-
Volume 7,
Issue 2,
2021
Volume 7, Issue 2, 2021
- Pathogens and Epidemiology
-
-
-
A comprehensive and high-quality collection of Escherichia coli genomes and their genes
Escherichia coli is a highly diverse organism that includes a range of commensal and pathogenic variants found across a range of niches and worldwide. In addition to causing severe intestinal and extraintestinal disease, E. coli is considered a priority pathogen due to high levels of observed drug resistance. The diversity in the E. coli population is driven by high genome plasticity and a very large gene pool. All these have made E. coli one of the most well-studied organisms, as well as a commonly used laboratory strain. Today, there are thousands of sequenced E. coli genomes stored in public databases. While data is widely available, accessing the information in order to perform analyses can still be a challenge. Collecting relevant available data requires accessing different sources, where data may be stored in a range of formats, and often requires further manipulation and processing to apply various analyses and extract useful information. In this study, we collated and intensely curated a collection of over 10 000 E. coli and Shigella genomes to provide a single, uniform, high-quality dataset. Shigella were included as they are considered specialized pathovars of E. coli . We provide these data in a number of easily accessible formats that can be used as the foundation for future studies addressing the biological differences between E. coli lineages and the distribution and flow of genes in the E. coli population at a high resolution. The analysis we present emphasizes our lack of understanding of the true diversity of the E. coli species, and the biased nature of our current understanding of the genetic diversity of such a key pathogen.
-
-
-
-
Neisseria gonorrhoeae clustering to reveal major European whole-genome-sequencing-based genogroups in association with antimicrobial resistance
Neisseria gonorrhoeae , the bacterium responsible for the sexually transmitted disease gonorrhoea, has shown an extraordinary ability to develop antimicrobial resistance (AMR) to multiple classes of antimicrobials. With no available vaccine, managing N. gonorrhoeae infections demands effective preventive measures, antibiotic treatment and epidemiological surveillance. The latter two are progressively being supported by the generation of whole-genome sequencing (WGS) data on behalf of national and international surveillance programmes. In this context, this study aims to perform N. gonorrhoeae clustering into genogroups based on WGS data, for enhanced prospective laboratory surveillance. Particularly, it aims to identify the major circulating WGS-genogroups in Europe and to establish a relationship between these and AMR. Ultimately, it enriches public databases by contributing with WGS data from Portuguese isolates spanning 15 years of surveillance. A total of 3791 carefully inspected N. gonorrhoeae genomes from isolates collected across Europe were analysed using a gene-by-gene approach (i.e. using cgMLST). Analysis of cluster composition and stability allowed the classification of isolates into a two-step hierarchical genogroup level determined by two allelic distance thresholds revealing cluster stability. Genogroup clustering in general agreed with available N. gonorrhoeae typing methods [i.e. MLST (multilocus sequence typing), NG-MAST ( N. gonorrhoeae multi-antigen sequence typing) and PubMLST core-genome groups], highlighting the predominant genogroups circulating in Europe, and revealed that the vast majority of the genogroups present a dominant AMR profile. Additionally, a non-static gene-by-gene approach combined with a more discriminatory threshold for potential epidemiological linkage enabled us to match data with previous reports on outbreaks or transmission chains. In conclusion, this genogroup assignment allows a comprehensive analysis of N. gonorrhoeae genetic diversity and the identification of the WGS-based genogroups circulating in Europe, while facilitating the assessment (and continuous monitoring) of their frequency, geographical dispersion and potential association with specific AMR signatures. This strategy may benefit public-health actions through the prioritization of genogroups to be controlled, the identification of emerging resistance carriage, and the potential facilitation of data sharing and communication.
-
-
-
Comparative genomics of Staphylococcus epidermidis from prosthetic-joint infections and nares highlights genetic traits associated with antimicrobial resistance, not virulence
More LessThere is increased awareness of the worldwide spread of specific epidemic multidrug-resistant (MDR) lineages of the human commensal Staphylococcus epidermidis . Here, using bioinformatic analyses accounting for population structure, we determined genomic traits (genes, SNPs and k-mers) that distinguish S. epidermidis causing prosthetic-joint infections (PJIs) from commensal isolates from nares, by analysing whole-genome sequencing data from S. epidermidis from PJIs prospectively collected over 10 years in Sweden, and contemporary S. epidermidis from the nares of patients scheduled for arthroplasty surgery. Previously suggested virulence determinants and the presence of genes and mutations linked to antimicrobial resistance (AMR) were also investigated. Publicly available S. epidermidis sequences were used for international extrapolation and validation of findings. Our data show that S. epidermidis causing PJIs differed from nasal isolates not by virulence but by traits associated with resistance to compounds used in prevention of PJIs: β-lactams, aminoglycosides and chlorhexidine. Almost a quarter of the PJI isolates did not belong to any of the previously described major nosocomial lineages, but the AMR-related traits were also over-represented in these isolates, as well as in international S. epidermidis isolates originating from PJIs. Genes previously associated with virulence in S. epidermidis were over-represented in individual lineages, but failed to reach statistical significance when adjusted for population structure. Our findings suggest that the current strategies for prevention of PJIs select for nosocomial MDR S. epidermidis lineages that have arisen from horizontal gene transfer of AMR-related traits into multiple genetic backgrounds.
-
-
-
Colonization with Staphylococcus aureus and Klebsiella pneumoniae causes infections in a Vietnamese intensive care unit
Pre-existing colonization with Staphylococcus aureus or Klebsiella pneumoniae has been found to increase the risk of infection in intensive care patients. We previously conducted a longitudinal study to characterize colonization of these two organisms in patients admitted to intensive care in a hospital in southern Vietnam. Here, using genomic and phylogenetic analyses, we aimed to assess the contribution these colonizing organisms made to infections. We found that in the majority of patients infected with S. aureus or K. pneumoniae , the sequence type of the disease-causing (infecting) isolate was identical to that of corresponding colonizing organisms in the respective patient. Further in-depth analysis revealed that in patients infected by S. aureus ST188 and by K. pneumoniae ST17, ST23, ST25 and ST86, the infecting isolate was closely related to and exhibited limited genetic variation relative to pre-infection colonizing isolates. Multidrug-resistant S. aureus ST188 was identified as the predominant agent of colonization and infection. Colonization and infection by K. pneumoniae were characterized by organisms with limited antimicrobial resistance profiles but extensive repertoires of virulence genes. Our findings augment the understanding of the link between bacterial colonization and infection in a low-resource setting, and could facilitate the development of novel evidence-based approaches to prevent and treat infections in high-risk patients in intensive care.
-
-
-
Phylogeny of Salmonella enterica subspecies arizonae by whole-genome sequencing reveals high incidence of polyphyly and low phase 1 H antigen variability
More LessSalmonella enterica subspecies arizonae is frequently associated with animal reservoirs, particularly reptiles, and can cause illness in some mammals, including humans. Using whole-genome sequencing data, core genome phylogenetic analyses were performed using 112 S . enterica subsp. arizonae isolates, representing 46 of 102 described serovars. Nearly one-third of these are polyphyletic, including two serovars that appear in four and five distinct evolutionary lineages. Subspecies arizonae has a monophasic H antigen. Among the 46 serovars investigated, only 8 phase 1 H antigens were identified, demonstrating high conservation for this antigen. Prophages and plasmids were found throughout this subspecies, including five novel prophages. Polyphyly was also reflected in prophage content, although some clade-specific enrichment for some phages was observed. IncFII(S) was the most frequent plasmid replicon identified and was found in a quarter of S. enterica subsp. arizonae genomes. Salmonella pathogenicity islands (SPIs) 1 and 2 are present across all Salmonella , including this subspecies, although effectors sipA, sptP and arvA in SPI-1 and sseG and ssaI in SPI-2 appear to be lost in this lineage. SPI-20, encoding a type VI secretion system, is exclusive to this subspecies and is well maintained in all genomes sampled. A number of fimbral operons were identified, including the sas operon that appears to be a synapomorphy for this subspecies, while others exhibited more clade-specific patterns. This work reveals evolutionary patterns in S. enterica subsp. arizonae that make this subspecies a unique lineage within this very diverse species.
-
-
-
Epidemiological analysis of pneumococcal strains isolated at Yangon Children’s Hospital in Myanmar via whole-genome sequencing-based methods
Streptococcus pneumoniae causes over one million deaths from lower respiratory infections per annum worldwide. Although mortality is very high in Southeast Asian countries, molecular epidemiological information remains unavailable for some countries. In this study, we report, for the first time, the whole-genome sequences and genetic profiles of pneumococcal strains isolated in Myanmar. We isolated 60 streptococcal strains from 300 children with acute respiratory infection at Yangon Children’s Hospital in Myanmar. We obtained whole-genome sequences and identified the species, serotypes, sequence types, antimicrobial resistance (AMR) profiles, virulence factor profiles and pangenome structure using sequencing-based analysis. Average nucleotide identity analysis indicated that 58 strains were S. pneumoniae and the other 2 strains were Streptococcus mitis . The major serotype was 19F (11 strains), followed by 6E (6B genetic variant; 7 strains) and 15 other serotypes; 5 untypable strains were also detected. Multilocus sequence typing analysis revealed 39 different sequence types, including 11 novel ones. In addition, genetic profiling indicated that AMR genes and mutations spread among pneumococcal strains in Myanmar. A minimum inhibitory concentration assay indicated that several pneumococcal strains had acquired azithromycin and tetracycline resistance, whereas no strains were found to be resistant against levofloxacin and high-dose penicillin G. Phylogenetic and pangenome analysis showed various pneumococcal lineages and that the pneumococcal strains contain a rich and mobile gene pool, providing them with the ability to adapt to selective pressures. This molecular epidemiological information can help in tracking global infection and supporting AMR control in addition to public health interventions in Myanmar.
-
- Functional Genomics and Microbe–Niche Interactions
-
-
-
Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history
Mireia Coscolla, Sebastien Gagneux, Fabrizio Menardo, Chloé Loiseau, Paula Ruiz-Rodriguez, Sonia Borrell, Isaac Darko Otchere, Adwoa Asante-Poku, Prince Asare, Leonor Sánchez-Busó, Florian Gehre, C. N’Dira Sanoussi, Martin Antonio, Dissou Affolabi, Janet Fyfe, Patrick Beckert, Stefan Niemann, Abraham S. Alabi, Martin P. Grobusch, Robin Kobbe, Julian Parkhill, Christian Beisel, Lukas Fenner, Erik C. Böttger, Conor J. Meehan, Simon R. Harris, Bouke C. de Jong, Dorothy Yeboah-Manu and Daniela BritesHuman tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum . Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum , and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.
-
-
-
-
Pathogenomic analyses of Mycobacterium microti, an ESX-1-deleted member of the Mycobacterium tuberculosis complex causing disease in various hosts
Mycobacterium microti is an animal-adapted member of the Mycobacterium tuberculosis complex (MTBC), which was originally isolated from voles, but has more recently also been isolated from other selected mammalian hosts, including occasionally from humans. Here, we have generated and analysed the complete genome sequences of five representative vole and clinical M. microti isolates using PacBio- and Illumina-based technologies, and have tested their virulence and vaccine potential in SCID (severe combined immune deficient) mouse and/or guinea pig infection models. We show that the clinical isolates studied here cluster separately in the phylogenetic tree from vole isolates and other clades from publicly available M. microti genome sequences. These data also confirm that the vole and clinical M. microti isolates were all lacking the specific RD1mic region, which in other tubercle bacilli encodes the ESX-1 type VII secretion system. Biochemical analysis further revealed marked phenotypic differences between isolates in type VII-mediated secretion of selected PE and PPE proteins, which in part were attributed to specific genetic polymorphisms. Infection experiments in the highly susceptible SCID mouse model showed that the clinical isolates were significantly more virulent than the tested vole isolates, but still much less virulent than the M. tuberculosis H37Rv control strain. The strong attenuation of the ATCC 35872 vole isolate in immunocompromised mice, even compared to the attenuated BCG (bacillus Calmette–Guérin) vaccine, and its historic use in human vaccine trials encouraged us to test this strain’s vaccine potential in a guinea pig model, where it demonstrated similar protective efficacy as a BCG control, making it a strong candidate for vaccination of immunocompromised individuals in whom BCG vaccination is contra-indicated. Overall, we provide new insights into the genomic and phenotypic variabilities and particularities of members of an understudied clade of the MTBC, which all share a recent common ancestor that is characterized by the deletion of the RD1mic region.
-
-
-
Acquisition of pcnB [poly(A) polymerase I] genes via horizontal transfer from the β, γ-Proteobacteria
More LessPoly(A) polymerases (PAPs) and tRNA nucleotidyltransferases belong to a superfamily of nucleotidyltransferases and modify RNA 3′-ends. The product of the pcnB gene, PAP I, has been characterized in a few β-, γ- and δ- Proteobacteria . Using the PAP I signature sequence, putative PAPs were identified in bacterial species from the α- and ε- Proteobacteria and from four other bacterial phyla ( Firmicutes , Actinobacteria , Bacteroidetes and Aquificae ). Phylogenetic analysis, alien index and G+C content calculations strongly suggest that the PAPs in the species identified in this study arose by horizontal gene transfer from the β- and γ- Proteobacteria .
-
-
-
Whole-genome comparative analysis of Malaysian Burkholderia pseudomallei clinical isolates
More LessBurkholderia pseudomallei , a soil-dwelling Gram-negative bacterium, is the causative agent of the endemic tropical disease melioidosis. Clinical manifestations of B. pseudomallei infection range from acute or chronic localized infection in a single organ to fulminant septicaemia in multiple organs. The diverse clinical manifestations are attributed to various factors, including the genome plasticity across B. pseudomallei strains. We previously characterized B. pseudomallei strains isolated in Malaysia and noted different levels of virulence in model hosts. We hypothesized that the difference in virulence might be a result of variance at the genome level. In this study, we sequenced and assembled four Malaysian clinical B. pseudomallei isolates, UKMR15, UKMPMC2000, UKMD286 and UKMH10. Phylogenomic analysis showed that Malaysian subclades emerged from the Asian subclade, suggesting that the Malaysian strains originated from the Asian region. Interestingly, the low-virulence strain, UKMH10, was the most distantly related compared to the other Malaysian isolates. Genomic island (GI) prediction analysis identified a new island of 23 kb, GI9c, which is present in B. pseudomallei and Burkholderia mallei , but not Burkholderia thailandensis . Genes encoding known B. pseudomallei virulence factors were present across all four genomes, but comparative analysis of the total gene content across the Malaysian strains identified 104 genes that are absent in UKMH10. We propose that these genes may encode novel virulence factors, which may explain the reduced virulence of this strain. Further investigation on the identity and role of these 104 proteins may aid in understanding B. pseudomallei pathogenicity to guide the design of new therapeutics for treating melioidosis.
-
- Genomic Methodologies
-
-
-
BSGatlas: a unified Bacillus subtilis genome and transcriptome annotation atlas with enhanced information access
A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis , as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5′ and 3′ UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR–Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.
-
-
- Evolution and Responses to Interventions
-
-
-
Rapid molecular evolution of Spiroplasma symbionts of Drosophila
Spiroplasma is a genus of Mollicutes whose members include plant pathogens, insect pathogens and endosymbionts of animals. Spiroplasma phenotypes have been repeatedly observed to be spontaneously lost in Drosophila cultures, and several studies have documented a high genomic turnover in Spiroplasma symbionts and plant pathogens. These observations suggest that Spiroplasma evolves quickly in comparison to other insect symbionts. Here, we systematically assess evolutionary rates and patterns of Spiroplasma poulsonii , a natural symbiont of Drosophila. We analysed genomic evolution of sHy within flies, and sMel within in vitro culture over several years. We observed that S. poulsonii substitution rates are among the highest reported for any bacteria, and around two orders of magnitude higher compared with other inherited arthropod endosymbionts. The absence of mismatch repair loci mutS and mutL is conserved across Spiroplasma , and likely contributes to elevated substitution rates. Further, the closely related strains sMel and sHy (>99.5 % sequence identity in shared loci) show extensive structural genomic differences, which potentially indicates a higher degree of host adaptation in sHy, a protective symbiont of Drosophila hydei. Finally, comparison across diverse Spiroplasma lineages confirms previous reports of dynamic evolution of toxins, and identifies loci similar to the male-killing toxin Spaid in several Spiroplasma lineages and other endosymbionts. Overall, our results highlight the peculiar nature of Spiroplasma genome evolution, which may explain unusual features of its evolutionary ecology.
-
-
Most Read This Month
