RESUMO
The analysis and comparison of genomes rely on diï¬erent tools for tasks such as annotation, orthology prediction, and phylogenetic inference. Most tools are specialized for a single task, and additional eï¬orts are necessary to integrate and visualize the results. To ï¬ll this gap, we developed zDB, an application integrating a Nextï¬ow analysis pipeline and a Python visualization platform built on the Django framework. The application is available on GitHub (https://github.com/metagenlab/zDB) and from the bioconda channel. Starting from annotated Genbank ï¬les, zDB identiï¬es orthologs and infers a phylogeny for each orthogroup. A species phylogeny is also constructed from shared single-copy orthologs. The results can be enriched with Pfam protein domain prediction, Cluster of Orthologs Genes and Kyoto Encyclopedia of Genes and Genomes annotations, and Swissprot homologs. The web application allows searching for speciï¬c genes or annotations, running Blast queries, and comparing genomic regions and whole genomes. The metabolic capacities of organisms can be compared at either the module or pathway levels. Finally, users can run queries to examine the conservation of speciï¬c genes or annotations across a chosen subset of genomes and display the results as a list of genes, Venn diagram, or heatmaps. Those features make zDB useful for both bioinformaticians and researchers more accustomed to laboratory research.IMPORTANCEGenome comparison and analysis rely on many independent tools, leaving to scientists the burden to integrate and visualize their results for interpretation. To alleviate this burden, we have built zDB, a comparative genomics tool that includes both an analysis pipeline and a visualization platform. The analysis pipeline automates gene annotation, orthology prediction, and phylogenetic inference, while the visualization platform allows scientists to easily explore the results in a web browser. Among other features, the interface allows users to visually compare whole genomes and targeted regions, assess the conservation of genes or metabolic pathways, perform Blast searches, or look for speciï¬c annotations. Altogether, this tool will be useful for a broad range of applications in comparative studies between two and hundred genomes. Furthermore, it is designed to allow sharing of data sets easily at a local or international scale, thereby supporting exploratory analyses for non-bioinformaticians on the genome of their favorite organisms.
Assuntos
Genoma Bacteriano , Genômica , Filogenia , Software , Genômica/métodos , Genoma Bacteriano/genética , Bactérias/genética , Bactérias/classificação , Anotação de Sequência Molecular/métodos , Biologia Computacional/métodos , Bases de Dados GenéticasRESUMO
The International Committee on Systematics of Prokaryotes (ICSP) discussed and rejected in 2020 a proposal to modify the International Code of Nomenclature of Prokaryotes to allow the use of gene sequences as type for naming prokaryotes. An alternative nomenclatural code, the Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode), which considers genome sequences as type material for naming species, was published in 2022. Members of the ICSP subcommittee for the taxonomy of the phylum Chlamydiae (Chlamydiota) consider that the use of gene sequences as type would benefit the taxonomy of microorganisms that are difficult to culture such as the chlamydiae and other strictly intracellular bacteria. We recommend the registration of new names of uncultured prokaryotes in the SeqCode registry.
RESUMO
Genome sequences from evolving infectious pathogens allow quantification of case introductions and local transmission dynamics. We sequenced 11,357 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from Switzerland in 2020-the sixth largest effort globally. Using a representative subset of these data, we estimated viral introductions to Switzerland and their persistence over the course of 2020. We contrasted these estimates with simple null models representing the absence of certain public health measures. We show that Switzerland's border closures decoupled case introductions from incidence in neighboring countries. Under a simple model, we estimate an 86 to 98% reduction in introductions during Switzerland's strictest border closures. Furthermore, the Swiss 2020 partial lockdown roughly halved the time for sampled introductions to die out. Last, we quantified local transmission dynamics once introductions into Switzerland occurred using a phylodynamic model. We found that transmission slowed 35 to 63% upon outbreak detection in summer 2020 but not in fall. This finding may indicate successful contact tracing over summer before overburdening in fall. The study highlights the added value of genome sequencing data for understanding transmission dynamics.
Assuntos
COVID-19 , SARS-CoV-2 , Humanos , SARS-CoV-2/genética , COVID-19/genética , Saúde Pública , Suíça/epidemiologia , Controle de Doenças Transmissíveis , Genoma Viral/genética , FilogeniaRESUMO
Background: The need for effective public health surveillance systems to track virus spread for targeted interventions was highlighted during the COVID-19 pandemic. It spurred an interest in the use of spatiotemporal clustering and genomic analyses to identify high-risk areas and track the spread of the SARS-CoV-2 virus. However, these two approaches are rarely combined in surveillance systems to complement each one's limitations; spatiotemporal clustering approaches usually consider only one source of virus transmission (i.e., the residential setting) to detect case clusters, while genomic studies require significant resources and processing time that can delay decision-making. Here, we clarify the differences and possible synergies of these two approaches in the context of infectious disease surveillance systems by investigating to what extent geographically-defined clusters are confirmed as transmission clusters based on genome sequences, and how genomic-based analyses can improve the epidemiological investigations associated with spatiotemporal cluster detection. Methods: For this purpose, we sequenced the SARS-CoV-2 genomes of 172 cases that were part of a collection of spatiotemporal clusters found in a Swiss state (Vaud) during the first epidemic wave. We subsequently examined intra-cluster genetic similarities and spatiotemporal distributions across virus genotypes. Results: Our results suggest that the congruence between the two approaches might depend on geographic features of the area (rural/urban) and epidemic context (e.g., lockdown). We also identified two potential superspreading events that started from cases in the main urban area of the state, leading to smaller spreading events in neighboring regions, as well as a large spreading in a geographically-isolated area. These superspreading events were characterized by specific mutations assumed to originate from Mulhouse and Milan, respectively. Our analyses propose synergistic benefits of using two complementary approaches in public health surveillance, saving resources and improving surveillance efficiency.
Assuntos
COVID-19 , Humanos , SARS-CoV-2/genética , Pandemias , Controle de Doenças Transmissíveis , Genômica , Análise por ConglomeradosRESUMO
SARS-CoV-2 is currently causing an unprecedented pandemic. While vaccines are massively deployed, we still lack effective large-scale antiviral therapies. In the quest for antivirals targeting conserved structures, we focused on molecules able to bind viral RNA secondary structures. Aminoglycosides are a class of antibiotics known to interact with the ribosomal RNA of both prokaryotes and eukaryotes and have previously been shown to exert antiviral activities by interacting with viral RNA. Here we show that the aminoglycoside geneticin is endowed with antiviral activity against all tested variants of SARS-CoV-2, in different cell lines and in a respiratory tissue model at non-toxic concentrations. The mechanism of action is an early inhibition of RNA replication and protein expression related to a decrease in the efficiency of the -1 programmed ribosomal frameshift (PRF) signal of SARS-CoV-2. Using in silico modeling, we have identified a potential binding site of geneticin in the pseudoknot of frameshift RNA motif. Moreover, we have selected, through virtual screening, additional RNA binding compounds, interacting with the same site with increased potency.
Assuntos
Tratamento Farmacológico da COVID-19 , Mudança da Fase de Leitura do Gene Ribossômico , Humanos , SARS-CoV-2/genética , Antivirais/farmacologia , Antivirais/química , RNA Viral/metabolismoRESUMO
Administration of plasma therapy may contribute to viral control and survival of COVID-19 patients receiving B-cell-depleting agents that impair humoral immunity. However, little is known on the impact of anti-CD20 pre-exposition on the kinetics of SARS-CoV-2-specific antibodies. Here, we evaluated the relationship between anti-spike immunoglobulin G (IgG) kinetics and the clinical status or intra-host viral evolution after plasma therapy in 36 eligible hospitalized COVID-19 patients, pre-exposed or not to B-cell-depleting treatments. The majority of anti-CD20 pre-exposed patients (14/17) showed progressive declines of anti-spike IgG titres following plasma therapy, contrasting with the 4/19 patients who had not received B-cell-depleting agents (p = 0.0006). Patients with antibody decay also depicted prolonged clinical symptoms according to the World Health Organization (WHO) severity classification (p = 0.0267) and SARS-CoV-2 viral loads (p = 0.0032) before complete virus clearance. Moreover, they had higher mutation rates than patients able to mount an endogenous humoral response (p = 0.015), including three patients with one to four spike mutations, potentially associated with immune escape. No relevant differences were observed between patients treated with plasma from convalescent and/or mRNA-vaccinated donors. Our study emphasizes the need for an individualized clinical care and follow-up in the management of COVID-19 patients with B-cell lymphopenia.
Assuntos
COVID-19 , Humanos , COVID-19/terapia , SARS-CoV-2 , Formação de Anticorpos , Imunização Passiva , Anticorpos Antivirais , Imunoglobulina GRESUMO
The rapid emergence of antibiotic-resistant bacteria poses a serious threat to public health worldwide. Antimicrobial peptides (AMPs) are promising antibiotic alternatives; however, little is known about bacterial mechanisms of AMP resistance and the interplay between AMP resistance and the bacterial response to other antimicrobials. In this study, we identified Escherichia coli mutants resistant to the TAT-RasGAP317-326 antimicrobial peptide and found that resistant bacteria show collateral sensitivity to other AMPs and antibacterial agents. We determined that resistance to TAT-RasGAP317-326 peptide arises through mutations in the histidine kinase EnvZ, a member of the EnvZ/OmpR two-component system responsible for osmoregulation in E. coli. In particular, we found that TAT-RasGAP317-326 binding and entry is compromised in E. coli peptide-resistant mutants. We showed that peptide resistance is associated with transcriptional regulation of a number of pathways and EnvZ-mediated resistance is dependent on the OmpR response regulator but is independent of the OmpC and OmpF outer membrane porins. Our findings provide insight into the bacterial mechanisms of TAT-RasGAP317-326 resistance and demonstrate that resistance to this AMP is associated with collateral sensitivity to other antibacterial agents. IMPORTANCE Antimicrobial peptides (AMP) are promising alternatives to classical antibiotics in the fight against antibiotic resistance. Resistance toward antimicrobial peptides can occur, but little is known about the mechanisms driving this phenomenon. Moreover, there is limited knowledge on how AMP resistance relates to the bacterial response to other antimicrobial agents. Here, we address these questions in the context of the antimicrobial peptide TAT-RasGAP317-326. We show that resistant Escherichia coli strains can be selected and do not show resistance to other antimicrobial agents. Resistance is caused by a mutation in a regulatory pathway, which lowers binding and entry of the peptide in E. coli. Our results highlight a mechanism of resistance that is specific to TAT-RasGAP317-326. Further research is required to characterize these mechanisms and to evaluate the potential of antimicrobial combinations to curb the development of antimicrobial resistance.
Assuntos
Proteínas de Bactérias , Proteínas de Escherichia coli , Escherichia coli , Transativadores , Antibacterianos/farmacologia , Proteínas da Membrana Bacteriana Externa/metabolismo , Proteínas de Bactérias/metabolismo , Sensibilidade Colateral a Medicamentos , Escherichia coli/genética , Escherichia coli/metabolismo , Proteínas de Escherichia coli/genética , Proteínas de Escherichia coli/metabolismo , Proteínas Ativadoras de GTPase , Complexos Multienzimáticos/metabolismo , Fragmentos de Peptídeos , Porinas/genética , Porinas/metabolismo , Proteínas Ativadoras de ras GTPase/metabolismo , Transativadores/metabolismoRESUMO
SARS-CoV-2 is currently causing an unprecedented pandemic. While vaccines are massively deployed, we still lack effective large-scale antiviral therapies. In the quest for antivirals targeting conserved structures, we focused on molecules able to bind viral RNA secondary structures. Aminoglycosides are a class of antibiotics known to interact with the ribosomal RNA of both prokaryotes and eukaryotes and have previously been shown to exert antiviral activities by interacting with viral RNA. Here we show that the aminoglycoside geneticin is endowed with antiviral activity against all tested variants of SARS-CoV-2, in different cell lines and in a respiratory tissue model at non-toxic concentrations. The mechanism of action is an early inhibition of RNA replication and protein expression related to a decrease in the efficiency of the -1 programmed ribosomal frameshift (PRF) signal of SARS-CoV-2. Using in silico modelling, we have identified a potential binding site of geneticin in the pseudoknot of frameshift RNA motif. Moreover, we have selected, through virtual screening, additional RNA binding compounds, interacting with the same site with increased potency.
RESUMO
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is currently causing an unprecedented pandemic. Although vaccines and antivirals are limiting the spread, SARS-CoV-2 is still under selective pressure in human and animal populations, as demonstrated by the emergence of variants of concern. To better understand the driving forces leading to new subtypes of SARS-CoV-2, we infected an ex vivo cell model of the human upper respiratory tract with Alpha and Omicron BA.1 variants for one month. Although viral RNA was detected during the entire course of the infection, infectious virus production decreased over time. Sequencing analysis did not show any adaptation in the spike protein, suggesting a key role for the adaptive immune response or adaptation to other anatomical sites for the evolution of SARS-CoV-2.
Assuntos
COVID-19 , SARS-CoV-2 , Animais , Humanos , Antivirais , Nariz , SARS-CoV-2/genética , Glicoproteína da Espícula de Coronavírus/genética , Traqueia , Evolução MolecularRESUMO
Antibiotic resistance is an increasing threat for public health, underscoring the need for new antibacterial agents. Antimicrobial peptides (AMPs) represent an alternative to classical antibiotics. TAT-RasGAP317-326 is a recently described AMP effective against a broad range of bacteria, but little is known about the conditions that may influence its activity. Using RNA-sequencing and screening of mutant libraries, we show that Escherichia coli and Pseudomonas aeruginosa respond to TAT-RasGAP317-326 by regulating metabolic and stress response pathways, possibly implicating two-component systems. Our results also indicate that bacterial surface properties, in particular integrity of the lipopolysaccharide layer, influence peptide binding and entry. Finally, we found differences between bacterial species with respect to their rate of resistance emergence against this peptide. Our findings provide the basis for future investigation on the mode of action of TAT-RasGAP317-326, which may help developing antimicrobial treatments based on this peptide.
RESUMO
Although many laboratories worldwide have developed their sequencing capacities in response to the need for SARS-CoV-2 genome-based surveillance of variants, only a few reported some quality criteria to ensure sequence quality before lineage assignment and submission to public databases. Hence, we aimed here to provide simple quality control criteria for SARS-CoV-2 sequencing to prevent erroneous interpretation of low-quality or contaminated data. We retrospectively investigated 647 SARS-CoV-2 genomes obtained over 10 tiled amplicons sequencing runs. We extracted 26 potentially relevant metrics covering the entire workflow from sample selection to bioinformatics analysis. Based on data distribution, critical values were established for 11 selected metrics to prompt further quality investigations for problematic samples, in particular those with a low viral RNA quantity. Low-frequency variants (<70% of supporting reads) can result from PCR amplification errors, sample cross contaminations, or presence of distinct SARS-CoV2 genomes in the sample sequenced. The number and the prevalence of low-frequency variants can be used as a robust quality criterion to identify possible sequencing errors or contaminations. Overall, we propose 11 metrics with fixed cutoff values as a simple tool to evaluate the quality of SARS-CoV-2 genomes, among which are cycle thresholds, mean depth, proportion of genome covered at least 10×, and the number of low-frequency variants combined with mutation prevalence data.
Assuntos
COVID-19 , SARS-CoV-2 , Genoma Viral , Humanos , RNA Viral , Estudos RetrospectivosRESUMO
The rapid spread of the SARS-CoV-2 lineages B.1.1.7 (N501Y.V1) throughout the UK, B.1.351 (N501Y.V2) in South Africa, and P.1 (B.1.1.28.1; N501Y.V3) in Brazil has led to the definition of variants of concern (VoCs) and recommendations for lineage specific surveillance. In Switzerland, during the last weeks of December 2020, we established a nationwide screening protocol across multiple laboratories, focusing first on epidemiological and microbiological definitions. In January 2021, we validated and implemented an N501Y-specific PCR to rapidly screen for VoCs, which are then confirmed using amplicon sequencing or whole genome sequencing (WGS). A total of 13,387 VoCs have been identified since the detection of the first Swiss case in October 2020, with 4194 being B.1.1.7, 172 B.1.351, and 7 P.1. The remaining 9014 cases of VoCs have been described without further lineage specification. Overall, all diagnostic centers reported a rapid increase of the percentage of detected VOCs, with a range of 6 to 46% between 25 to 31 of January 2021 increasing towards 41 to 82% between 22 to 28 of February. A total of 739 N501Y positive genomes were analysed and show a broad range of introduction events to Switzerland. In this paper, we describe the nationwide coordination and implementation process across laboratories, public health institutions, and researchers, the first results of our N501Y-specific variant screening, and the phylogenetic analysis of all available WGS data in Switzerland, that together identified the early introduction events and subsequent community spreading of the VoCs.
RESUMO
The order Chlamydiales includes obligate intracellular pathogens capable of infecting mammals, fishes and amoeba. Unlike other intracellular bacteria for which intracellular adaptation led to the loss of glycogen metabolism pathway, all chlamydial families maintained the nucleotide-sugar dependent glycogen metabolism pathway i.e. the GlgC-pathway with the notable exception of both Criblamydiaceae and Waddliaceae families. Through detailed genome analysis and biochemical investigations, we have shown that genome rearrangement events have resulted in a defective GlgC-pathway and more importantly we have evidenced a distinct trehalose-dependent GlgE-pathway in both Criblamydiaceae and Waddliaceae families. Altogether, this study strongly indicates that the glycogen metabolism is retained in all Chlamydiales without exception, highlighting the pivotal function of storage polysaccharides, which has been underestimated to date. We propose that glycogen degradation is a mandatory process for fueling essential metabolic pathways that ensure the survival and virulence of extracellular forms i.e. elementary bodies of Chlamydiales.
Assuntos
Chlamydiales/metabolismo , Glicogênio/metabolismo , Glicogenólise , Polissacarídeos Bacterianos/metabolismo , Chlamydiales/genética , Chlamydiales/patogenicidade , Evolução Molecular , Regulação Bacteriana da Expressão Gênica , Genoma Bacteriano , Cinética , Filogenia , VirulênciaRESUMO
The Mycobacterium kansasii species comprises six subtypes that were recently classified into six closely related species; Mycobacterium kansasii (formerly M. kansasii subtype 1), Mycobacterium persicum (subtype 2), Mycobacterium pseudokansasii (subtype 3), Mycobacterium ostraviense (subtype 4), Mycobacterium innocens (subtype 5) and Mycobacterium attenuatum (subtype 6). Together with Mycobacterium gastri, they form the M. kansasii complex. M. kansasii is the most frequent and most pathogenic species of the complex. M. persicum is classically associated with diseases in immunosuppressed patients, and the other species are mostly colonizers, and are only very rarely reported in ill patients. Comparative genomics was used to assess the genetic determinants leading to the pathogenicity of members of the M. kansasii complex. The genomes of 51 isolates collected from patients with and without disease were sequenced and compared with 24 publicly available genomes. The pathogenicity of each isolate was determined based on the clinical records or public metadata. A comparative genomic analysis showed that all M. persicum, M. ostraviense, M innocens and M. gastri isolates lacked the ESX-1-associated EspACD locus that is thought to play a crucial role in the pathogenicity of M. tuberculosis and other non-tuberculous mycobacteria. Furthermore, M. kansasii was the only species exhibiting a 25-Kb-large genomic island encoding for 17 type-VII secretion system-associated proteins. Finally, a genome-wide association analysis revealed that two consecutive genes encoding a hemerythrin-like protein and a nitroreductase-like protein were significantly associated with pathogenicity. These two genes may be involved in the resistance to reactive oxygen and nitrogen species, a required mechanism for the intracellular survival of bacteria. Three non-pathogenic M. kansasii lacked these genes likely due to two distinct distributive conjugal transfers (DCTs) between M. attenuatum and M. kansasii, and one DCT between M. persicum and M. kansasii. To our knowledge, this is the first study linking DCT to reduced pathogenicity.
RESUMO
Chronic infections caused by obligate intracellular bacteria belonging to the Chlamydiales order are related to the formation of persistent developmental forms called aberrant bodies (ABs), which undergo DNA replication without cell division. These enlarged bacteria develop and persist upon exposure to different stressful conditions such as ß-lactam antibiotics, iron deprivation and interferon-γ. However, the mechanisms behind ABs biogenesis remain uncharted. Using an RNA-sequencing approach, we compared the transcriptional profile of ABs induced by iron starvation to untreated bacteria in the Chlamydia-related species Waddliachondrophila, a potential agent of abortion in ruminants and miscarriage in humans. Consistent with the growth arrest observed following iron depletion, our results indicate a significant reduction in the expression of genes related to energy production, carbohydrate and amino acid metabolism and cell wall/envelope biogenesis, compared to untreated, actively replicating bacteria. Conversely, three putative toxin-antitoxin modules were among the most up-regulated genes upon iron starvation, suggesting that their activation might be involved in growth arrest in adverse conditions, an uncommon feature in obligate intracellular bacteria. Our work represents the first complete transcriptomic profile of a Chlamydia-related species in stressful conditions and sets the grounds for further investigations on the mechanisms underlying chlamydial persistence.
Assuntos
Genômica/educação , Técnicas Microbiológicas , Biologia Computacional , Resistência Microbiana a Medicamentos/genética , Educação Médica Continuada , Sequenciamento de Nucleotídeos em Larga Escala , Metagenômica/educação , Tipagem Molecular , Filogenia , Alinhamento de Sequência , Análise de Sequência , Virulência/genéticaRESUMO
Whole genome sequencing (WGS) enables high resolution typing of bacteria up to the single nucleotide polymorphism (SNP) level. WGS is used in clinical microbiology laboratories for infection control, molecular surveillance and outbreak analyses. Given the large palette of WGS reagents and bioinformatics tools, the Swiss clinical bacteriology community decided to conduct a ring trial (RT) to foster harmonization of NGS-based bacterial typing. The RT aimed at assessing methicillin-susceptible Staphylococcus aureus strain relatedness from WGS and epidemiological data. The RT was designed to disentangle the variability arising from differences in sample preparation, SNP calling and phylogenetic methods. Nine laboratories participated. The resulting phylogenetic tree and cluster identification were highly reproducible across the laboratories. Cluster interpretation was, however, more laboratory dependent, suggesting that an increased sharing of expertise across laboratories would contribute to further harmonization of practices. More detailed bioinformatic analyses unveiled that while similar clusters were found across laboratories, these were actually based on different sets of SNPs, differentially retained after sample preparation and SNP calling procedures. Despite this, the observed number of SNP differences between pairs of strains, an important criterion to determine strain relatedness given epidemiological information, was similar across pipelines for closely related strains when restricting SNP calls to a common core genome defined by S. aureus cgMLST schema. The lessons learned from this pilot study will serve the implementation of larger-scale RT, as a mean to have regular external quality assessments for laboratories performing WGS analyses in a clinical setting.