The limitations of draft assemblies for understanding prokaryotic adaptation and evolution

Dellaglioa spp. an underestimated genus isolated from high-oxygen modified-atmosphere packaged meat

The genus Dellaglioa (D.) actually comprises two species, i.e., D. algida and the recently described species D. carnosa. Both species are adapted to cold and have been typically recovered from meat products. However, their importance has thus far been underestimated, since routine culture-based analysis failed to support their growth. Furthermore, their occurrence on meat packed under high-oxygen MA conditions (HiOx-MAP) is controversial because they have been described as being oxygen-sensitive. In this study, we focused on the targeted isolation of Dellaglioa spp. from HiOx-MAP meat samples and the characterization of our isolates regarding their adaption to HiOx-MAP conditions, their spoilage potential, as well as food safety aspects. We used a medium recently developed specifically for strains of this genus and investigated ten meat batches from seven different suppliers. Our study confirms that the occurrence of Dellaglioa spp. on HiOx-MAP meat is non-sporadic, reaching cell counts ranging from log10 5.8-7.1 CFU/cm2 at a late stage of chilled storage. Autochthonous Dellaglioa spp. and Leuconostoc (L.) gasicomitatum dominated the microbiota of the beef steaks with similar growth behavior. Our results suggest that Dellaglioa spp. benefits from the heme-dependent respiration of oxygen by L. gasicomitatum. Furthermore, whole genome analysis revealed the presence of genes predictively involved in oxidative stress defense, survival, and adaptation in meat environments. Moreover, we predict a weak aminogenic potential of D. algida strains. Tyramine production from tyrosine seems to be a species-specific characteristic of D. carnosa. The extent to which D. algida and D. carnosa occurrence is influenced by or even dependent on the composition of the entire microbiota remains to be investigated.

Development of pBACpAK entrapment vector derivatives to detect intracellular transfer of mobile genetic elements within chloramphenicol resistant bacterial isolates

Antimicrobial resistance disseminates throughout bacterial populations via horizontal gene transfer, driven mainly by mobile genetic elements (MGEs). Entrapment vectors are key tools in determining MGE movement within a bacterial cell between different replicons or between sites within the same replicon. The pBACpAK entrapment vector has been previously used to study intracellular transfer in Gram-negative bacteria however since pBACpAK contains a chloramphenicol resistance gene, it cannot be used in bacterial isolates which are already resistant to chloramphenicol. Therefore, we developed new derivatives of the pBACpAK entrapment vector to determine intracellular transfer of MGEs in an Escherichia coli DH5α transconjugant containing the chloramphenicol resistance plasmid pD25466. The catA1 of pBACpAK was replaced by both mcr-1 in pBACpAK-COL and aph(3')-Ia in pBACpAK-KAN, allowing it to be used in chloramphenicol resistant strains. The plasmid constructs were verified and then used to transform the E. coli DH5α/pD25466 transconjugants in order to detect intracellular movement of the MGEs associated with the pD25466 plasmid. Here we report on the validation of the expanded suite of pBACpAK vectors which can be used to study the intracellular transfer of MGEs between, and within, replicons in bacteria with different antimicrobial resistance profiles.

Genomic analysis of 61 Chlamydia psittaci strains reveals extensive divergence associated with host preference

BACKGROUND

Chlamydia (C.) psittaci, the causative agent of avian chlamydiosis and human psittacosis, is a genetically heterogeneous species. Its broad host range includes parrots and many other birds, but occasionally also humans (via zoonotic transmission), ruminants, horses, swine and rodents. To assess whether there are genetic markers associated with host tropism we comparatively analyzed whole-genome sequences of 61 C. psittaci strains, 47 of which carrying a 7.6-kbp plasmid.

RESULTS

Following clean-up, reassembly and polishing of poorly assembled genomes from public databases, phylogenetic analyses using C. psittaci whole-genome sequence alignment revealed four major clades within this species. Clade 1 represents the most recent lineage comprising 40/61 strains and contains 9/10 of the psittacine strains, including type strain 6BC, and 10/13 of human isolates. Strains from different non-psittacine hosts clustered in Clades 2- 4. We found that clade membership correlates with typing schemes based on SNP types, ompA genotypes, multilocus sequence types as well as plasticity zone (PZ) structure and host preference. Genome analysis also revealed that i) sequence variation in the major outer membrane porin MOMP can result in 3D structural changes of immunogenic domains, ii) past host change of Clade 3 and 4 strains could be associated with loss of MAC/perforin in the PZ, rather than the large cytotoxin, iii) the distinct phylogeny of atypical strains (Clades 3 and 4) is also reflected in their repertoire of inclusion proteins (Inc family) and polymorphic membrane proteins (Pmps).

CONCLUSIONS

Our study identified a number of genomic features that can be correlated with the phylogeny and host preference of C. psittaci strains. Our data show that intra-species genomic divergence is associated with past host change and includes deletions in the plasticity zone, structural variations in immunogenic domains and distinct repertoires of virulence factors.

Diversity of the type VI secretion systems in the Neisseria spp

Complete Type VI Secretion Systems were identified in the genome sequence data of Neisseria subflava isolates sourced from throat swabs of human volunteers. The previous report was the first to describe two complete Type VI Secretion Systems in these isolates, both of which were distinct in terms of their gene organization and sequence homology. Since publication of the first report, Type VI Secretion System subtypes have been identified in Neisseria spp. The characteristics of each type in N. subflava are further investigated here and in the context of the other Neisseria spp., including identification of the lineages containing the different types and subtypes. Type VI Secretion Systems use VgrG for delivery of toxin effector proteins; several copies of vgrG and associated effector / immunity pairs are present in Neisseria spp. Based on sequence similarity between strains and species, these core Type VI Secretion System genes, vgrG, and effector / immunity genes may diversify via horizontal gene transfer, an instrument for gene acquisition and repair in Neisseria spp.

Exploring the mobilome and resistome of Enterococcus faecium in a One Health context across two continents

Enterococcus faecium is a ubiquitous opportunistic pathogen that is exhibiting increasing levels of antimicrobial resistance (AMR). Many of the genes that confer resistance and pathogenic functions are localized on mobile genetic elements (MGEs), which facilitate their transfer between lineages. Here, features including resistance determinants, virulence factors and MGEs were profiled in a set of 1273 E. faecium genomes from two disparate geographic locations (in the UK and Canada) from a range of agricultural, clinical and associated habitats. Neither lineages of E. faecium, type A and B, nor MGEs are constrained by geographic proximity, but our results show evidence of a strong association of many profiled genes and MGEs with habitat. Many features were associated with a group of clinical and municipal wastewater genomes that are likely forming a new human-associated ecotype within type A. The evolutionary dynamics of E. faecium make it a highly versatile emerging pathogen, and its ability to acquire, transmit and lose features presents a high risk for the emergence of new pathogenic variants and novel resistance combinations. This study provides a workflow for MGE-centric surveillance of AMR in Enterococcus that can be adapted to other pathogens.

Software evaluation for de novo detection of transposons

Transposable elements (TEs) are major genomic components in most eukaryotic genomes and play an important role in genome evolution. However, despite their relevance the identification of TEs is not an easy task and a number of tools were developed to tackle this problem. To better understand how they perform, we tested several widely used tools for de novo TE detection and compared their performance on both simulated data and well curated genomic sequences. As expected, tools that build TE-models performed better than k-mer counting ones, with RepeatModeler beating competitors in most datasets. However, there is a tendency for most tools to identify TE-regions in a fragmented manner and it is also frequent that small TEs or fragmented TEs are not detected. Consequently, the identification of TEs is still a challenging endeavor and it requires a significant manual curation by an experienced expert. The results will be helpful for identifying common issues associated with TE-annotation and for evaluating how comparable are the results obtained with different tools.

Accelerating Complete Phytoplasma Genome Assembly by Immunoprecipitation-Based Enrichment and MinION-Based DNA Sequencing for Comparative Analyses

Phytoplasmas are uncultivated plant-pathogenic bacteria with agricultural importance. Those belonging to the 16SrII group, represented by 'Candidatus P. aurantifolia', have a wide range of plant hosts and cause significant yield losses in valuable crops, such as pear, sweet potato, peanut, and soybean. In this study, a method that combines immunoprecipitation-based enrichment and MinION long-read DNA sequencing was developed to solve the challenge of phytoplasma genome studies. This approach produced long reads with high mapping rates and high genomic coverage that can be combined with Illumina reads to produce complete genome assemblies with high accuracy. We applied this method to strain NCHU2014 and determined its complete genome sequence, which consists of one circular chromosome with 635,584 bp and one plasmid with 4,224 bp. Although 'Ca. P. aurantifolia' NCHU2014 has a small chromosome with only 471 protein-coding genes, it contains 33 transporter genes and 27 putative effector genes, which may contribute to obtaining nutrients from hosts and manipulating host developments for their survival and multiplication. Two effectors, the homologs of SAP11 and SAP54/PHYL1 identified in 'Ca. P. aurantifolia' NCHU2014, have the biochemical activities in destabilizing host transcription factors, which can explain the disease symptoms observed in infected plants. Taken together, this study provides the first complete genome available for the 16SrII phytoplasmas and contributes to the understanding of phytoplasma pathogenicity.

Genome-Inferred Correspondence between Phylogeny and Metabolic Traits in the Wild Drosophila Gut Microbiome

Annotated genome sequences provide valuable insight into the functional capabilities of members of microbial communities. Nevertheless, most studies on the microbiome in animal guts use metagenomic data, hampering the assignment of genes to specific microbial taxa. Here, we make use of the readily culturable bacterial communities in the gut of the fruit fly Drosophila melanogaster to obtain draft genome sequences for 96 isolates from wild flies. These include 81 new de novo assembled genomes, assigned to three orders (Enterobacterales, Lactobacillales, and Rhodospirillales) with 80% of strains identified to species level using average nucleotide identity and phylogenomic reconstruction. Based on annotations by the RAST pipeline, among-isolate variation in metabolic function partitioned strongly by bacterial order, particularly by amino acid metabolism (Rhodospirillales), fermentation, and nucleotide metabolism (Lactobacillales) and arginine, urea, and polyamine metabolism (Enterobacterales). Seven bacterial species, comprising 2-3 species in each order, were well-represented among the isolates and included ≥5 strains, permitting analysis of metabolic functions in the accessory genome (i.e., genes not present in every strain). Overall, the metabolic function in the accessory genome partitioned by bacterial order. Two species, Gluconobacter cerinus (Rhodospirillales) and Lactiplantibacillus plantarum (Lactobacillales) had large accessory genomes, and metabolic functions were dominated by amino acid metabolism (G. cerinus) and carbohydrate metabolism (La. plantarum). The patterns of variation in metabolic capabilities at multiple phylogenetic scales provide the basis for future studies of the ecological and evolutionary processes shaping the diversity of microorganisms associated with natural populations of Drosophila.

AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence

Antimicrobial resistance (AMR) is a significant public health threat. With the rise of affordable whole genome sequencing, in silico approaches to assessing AMR gene content can be used to detect known resistance mechanisms and potentially identify novel mechanisms. To enable accurate assessment of AMR gene content, as part of a multi-agency collaboration, NCBI developed a comprehensive AMR gene database, the Bacterial Antimicrobial Resistance Reference Gene Database and the AMR gene detection tool AMRFinder. Here, we describe the expansion of the Reference Gene Database, now called the Reference Gene Catalog, to include putative acid, biocide, metal, stress resistance genes, in addition to virulence genes and species-specific point mutations. Genes and point mutations are classified by broad functions, as well as more detailed functions. As we have expanded both the functional repertoire of identified genes and functionality, NCBI released a new version of AMRFinder, known as AMRFinderPlus. This new tool allows users the option to utilize only the core set of AMR elements, or include stress response and virulence genes, too. AMRFinderPlus can detect acquired genes and point mutations in both protein and nucleotide sequence. In addition, the evidence used to identify the gene has been expanded to include whether nucleotide or protein sequence was used, its location in the contig, and presence of an internal stop codon. These database improvements and functional expansions will enable increased precision in identifying AMR genes, linking AMR genotypes and phenotypes, and determining possible relationships between AMR, virulence, and stress response.

Comparative Genome Analysis of 33 Chlamydia Strains Reveals Characteristic Features of Chlamydia Psittaci and Closely Related Species

To identify genome-based features characteristic of the avian and human pathogen Chlamydia(C.) psittaci and related chlamydiae, we analyzed whole-genome sequences of 33 strains belonging to 12 species. Using a novel genome analysis tool termed Roary ILP Bacterial Annotation Pipeline (RIBAP), this panel of strains was shown to share a large core genome comprising 784 genes and representing approximately 80% of individual genomes. Analyzing the most variable genomic sites, we identified a set of features of C. psittaci that in its entirety is characteristic of this species: (i) a relatively short plasticity zone of less than 30,000 nt without a tryptophan operon (also in C. abortus, C. avium, C. gallinacea, C. pneumoniae), (ii) a characteristic set of of Inc proteins comprising IncA, B, C, V, X, Y (with homologs in C. abortus, C. caviae and C. felis as closest relatives), (iii) a 502-aa SinC protein, the largest among Chlamydia spp., and (iv) an elevated number of Pmp proteins of subtype G (14 in C. psittaci, 14 in Cand. C. ibidis). In combination with future functional studies, the common and distinctive criteria revealed in this study provide important clues for understanding the complexity of host-specific behavior of individual Chlamydia spp.

Lysogeny in Streptococcus pneumoniae

Bacterial viruses, or bacteriophages, are major contributors to the evolution, pathogenesis and overall biology of their host bacteria. During their life cycle, temperate bacteriophages form stable associations with their host by integrating into the chromosome, a process called lysogeny. Isolates of the human pathogen Streptococcus pneumoniae are frequently lysogenic, and genomic studies have allowed the classification of these phages into distinct phylogenetic groups. Here, we review the recent advances in the characterization of temperate pneumococcal phages, with a focus on their genetic features and chromosomal integration loci. We also discuss the contribution of phages, and specific phage-encoded features, to colonization and virulence. Finally, we discuss interesting research perspectives in this field.

Potential evolutionary impact of integrative and conjugative elements (ICEs) and genomic islands in the Ralstonia solanacearum species complex

Ralstonia solanacearum, a soil-borne plant pathogen, encompasses a large number of strains known as R. solanacearum species complex (RSSC). Although it has been suggested that mobile genetic elements (MGEs) may play an important role in the RSSC genome, the evolutionary impact of these elements remains unknown. Here, we identified and analysed Integrative and Conjugative Elements (ICEs) and Genomic Islands (GIs) in the 121 genomes published for Ralstonia spp., including RSSC strains and three other non-plant pathogenic Ralstonia spp. Our results provided a dataset of 12 ICEs and 31 GIs distributed throughout Ralstonia spp. Four novel ICEs in RSSC were found. Some of these elements cover 5% of the host genome and carry accessory genes with a potential impact on the fitness and pathogenicity of RSSC. In addition, phylogenetic analysis revealed that these MGEs clustered to the same species, but there is evidence of strains from different countries that host the same element. Our results provide novel insight into the RSSC adaptation, opening new paths to a better understanding of how these elements affect this soil-borne plant pathogen.

Comparative genomics of multidrug-resistant Enterococcus spp. isolated from wastewater treatment plants

BACKGROUND

Wastewater treatment plants (WWTPs) are considered hotspots for the environmental dissemination of antimicrobial resistance (AMR) determinants. Vancomycin-Resistant Enterococcus (VRE) are candidates for gauging the degree of AMR bacteria in wastewater. Enterococcus faecalis and Enterococcus faecium are recognized indicators of fecal contamination in water. Comparative genomics of enterococci isolated from conventional activated sludge (CAS) and biological aerated filter (BAF) WWTPs was conducted.

RESULTS

VRE isolates, including E. faecalis (n = 24), E. faecium (n = 11), E. casseliflavus (n = 2) and E. gallinarum (n = 2) were selected for sequencing based on WWTP source, species and AMR phenotype. The pangenomes of E. faecium and E. faecalis were both open. The genomic fraction related to the mobilome was positively correlated with genome size in E. faecium (p < 0.001) and E. faecalis (p < 0.001) and with the number of AMR genes in E. faecium (p = 0.005). Genes conferring vancomycin resistance, including vanA and vanM (E. faecium), vanG (E. faecalis), and vanC (E. casseliflavus/E. gallinarum), were detected in 20 genomes. The most prominent functional AMR genes were efflux pumps and transporters. A minimum of 16, 6, 5 and 3 virulence genes were detected in E. faecium, E. faecalis, E. casseliflavus and E. gallinarum, respectively. Virulence genes were more common in E. faecalis and E. faecium, than E. casseliflavus and E. gallinarum. A number of mobile genetic elements were shared among species. Functional CRISPR/Cas arrays were detected in 13 E. faecalis genomes, with all but one also containing a prophage. The lack of a functional CRISPR/Cas arrays was associated with multi-drug resistance in E. faecium. Phylogenetic analysis demonstrated differential clustering of isolates based on original source but not WWTP. Genes related to phage and CRISPR/Cas arrays could potentially serve as environmental biomarkers.

CONCLUSIONS

There was no discernible difference between enterococcal genomes from the CAS and BAF WWTPs. E. faecalis and E. faecium have smaller genomes and harbor more virulence, AMR, and mobile genetic elements than other Enterococcus spp.

MinION sequencing of Streptococcus suis allows for functional characterization of bacteria by multilocus sequence typing and antimicrobial resistance profiling

In recent years, high-throughput sequencing has revolutionized disease diagnosis by its powerful ability to provide high resolution genomic information. The Oxford Nanopore MinION sequencer has unparalleled potential as a rapid disease diagnostic tool due to its high mobility, accessibility, and short turnaround time. However, there is a lack of rigorous quality assessment and control processes standardizing the testing on the MinION, which is necessary for incorporation into a diagnostic workflow. Thus, our study examined the use of the MinION sequencer for bacterial whole genome generation and characterization. Using Streptococcus suis as a model, we optimized DNA isolation and treatments to be used for MinION sequencing and standardized de novo assembly to quickly generate a full-length consensus sequence achieving a 99.4% average accuracy. The consensus genomes from MinION sequencing were able to accurately predict the multilocus sequence type in 8 out of 10 samples and identified antimicrobial resistance profiles for 100% of the samples, despite the concern of a high error rate. The inability to unequivocally predict sequence types was due to difficulty in differentiating high identity alleles, which was overcome by applying additional error correction methods to increase consensus accuracy. This manuscript provides methods for the use of MinION sequencing for identification of S. suis genome sequence, sequence type, and antibiotic resistance profile that can be used as a framework for identification and classification of other pathogens.

Multiplexed Non-barcoded Long-Read Sequencing and Assembling Genomes of Bacillus Strains in Error-Free Simulations

The generation of genomic data from microorganisms has revolutionized our abilities to understand their biology, but it is still challenging to obtain complete genome sequences of microbes in an automated high-throughput and cost-effective manner. While the advent of second-generation sequencing technologies provided significantly higher throughput, their shorter lengths and more pronounced sequence-context bias led to a shift towards resequencing applications. Recently, single molecule real-time (SMRT) DNA sequencing has been used to generate sequencing reads that are much longer than other sequencing platforms, facilitating de novo genome assembly and genome finishing. Here we introduced a novel multiplex strategy to make full use of the capacity and characteristics of SMRT sequencing in microbe genome assembly. We used error-free simulations to evaluate the practicability of assembling SMRT genomic sequencing data from multiple microbes into finished genomes once at a time. Then we compared the influence of two key factors, including sequencing coverage and read length, on multiplex assembling. Our results showed that long-read genomic sequencing inherently provided the ability to assemble genomic sequencing data from multiple microbes into finished genomes due to its long length. This approach might be helpful for the various groups of microbial genome projects or metagenomics research.

Complete hybrid genome assembly of clinical multidrug-resistant Bacteroides fragilis isolates enables comprehensive identification of antimicrobial-resistance genes and plasmids

Bacteroides fragilis constitutes a significant part of the normal human gut microbiota and can also act as an opportunistic pathogen. Antimicrobial resistance (AMR) and the prevalence of AMR genes are increasing, and prediction of antimicrobial susceptibility based on sequence information could support targeted antimicrobial therapy in a clinical setting. Complete identification of insertion sequence (IS) elements carrying promoter sequences upstream of resistance genes is necessary for prediction of AMR. However, de novo assemblies from short reads alone are often fractured due to repeat regions and the presence of multiple copies of identical IS elements. Identification of plasmids in clinical isolates can aid in the surveillance of the dissemination of AMR, and comprehensive sequence databases support microbiome and metagenomic studies. We tested several short-read, hybrid and long-lead assembly pipelines by assembling the type strain B. fragilis CCUG4856T (=ATCC25285=NCTC9343) with Illumina short reads and long reads generated by Oxford Nanopore Technologies (ONT) MinION sequencing. Hybrid assembly with Unicycler, using quality filtered Illumina reads and Filtlong filtered and Canu-corrected ONT reads, produced the assembly of highest quality. This approach was then applied to six clinical multidrug-resistant B. fragilis isolates and, with minimal manual finishing of chromosomal assemblies of three isolates, complete, circular assemblies of all isolates were produced. Eleven circular, putative plasmids were identified in the six assemblies, of which only three corresponded to a known cultured Bacteroides plasmid. Complete IS elements could be identified upstream of AMR genes; however, there was not complete correlation between the absence of IS elements and antimicrobial susceptibility. As our knowledge on factors that increase expression of resistance genes in the absence of IS elements is limited, further research is needed prior to implementing AMR prediction for B. fragilis from whole-genome sequencing.

Rapid Detection of Genetic Engineering, Structural Variation, and Antimicrobial Resistance Markers in Bacterial Biothreat Pathogens by Nanopore Sequencing

Widespread release of Bacillus anthracis (anthrax) or Yersinia pestis (plague) would prompt a public health emergency. During an exposure event, high-quality whole genome sequencing (WGS) can identify genetic engineering, including the introduction of antimicrobial resistance (AMR) genes. Here, we developed rapid WGS laboratory and bioinformatics workflows using a long-read nanopore sequencer (MinION) for Y. pestis (6.5 h) and B. anthracis (8.5 h) and sequenced strains with different AMR profiles. Both salt-precipitation and silica-membrane extracted DNA were suitable for MinION WGS using both rapid and field library preparation methods. In replicate experiments, nanopore quality metrics were defined for genome assembly and mutation analysis. AMR markers were correctly detected and >99% coverage of chromosomes and plasmids was achieved using 100,000 raw sequencing reads. While chromosomes and large and small plasmids were accurately assembled, including novel multimeric forms of the Y. pestis virulence plasmid, pPCP1, MinION reads were error-prone, particularly in homopolymer regions. MinION sequencing holds promise as a practical, front-line strategy for on-site pathogen characterization to speed the public health response during a biothreat emergency.

Pushing the limits of de novo genome assembly for complex prokaryotic genomes harboring very long, near identical repeats

Generating a complete, de novo genome assembly for prokaryotes is often considered a solved problem. However, we here show that Pseudomonas koreensis P19E3 harbors multiple, near identical repeat pairs up to 70 kilobase pairs in length, which contained several genes that may confer fitness advantages to the strain. Its complex genome, which also included a variable shufflon region, could not be de novo assembled with long reads produced by Pacific Biosciences’ technology, but required very long reads from Oxford Nanopore Technologies. Importantly, a repeat analysis, whose results we release for over 9600 prokaryotes, indicated that very complex bacterial genomes represent a general phenomenon beyond Pseudomonas. Roughly 10% of 9331 complete bacterial and a handful of 293 complete archaeal genomes represented this ‘dark matter’ for de novo genome assembly of prokaryotes. Several of these ‘dark matter’ genome assemblies contained repeats far beyond the resolution of the sequencing technology employed and likely contain errors, other genomes were closed employing labor-intense steps like cosmid libraries, primer walking or optical mapping. Using very long sequencing reads in combination with assembly algorithms capable of resolving long, near identical repeats will bring most prokaryotic genomes within reach of fast and complete de novo genome assembly.

Genomic Characterization Provides New Insights for Detailed Phage- Resistant Mechanism for Brucella abortus

As the causative agent of cattle brucellosis, Brucella abortus commonly exhibits smooth phenotype (by virtue of colony morphology) that is characteristically sensitive to specific Brucella phages, playing until recently a major role in taxonomical classification of the Brucella species by the phage typing approach. We previously reported the discrepancy between traditional phenotypic typing and MLVA results of a smooth phage-resistant (SPR) strain Bab8416 isolated from a 45-year-old custodial worker with brucellosis in a cattle farm. Here, we performed whole genome sequencing and further obtained a complete genome sequence of strain Bab8416 by a combination of multiple NGS technologies and routine PCR sequencing. The detailed genetic differences between B. abortus SPR Bab8416 and large smooth phage-sensitive (SPS) strains were investigated in a comprehensively comparative genomic study. The large indels between B. abortus SPS strains and Bab8416 showed possible divergence between two evolutionary branches at a far phylogenetic node. Compared to B. abortus SPS strain 9-941 (Bab9-941), the specific re-arrangement event in Bab8416 displaying a closer linear relationship with B. melitensis 16M than other B. abortus strains resulted in the truncation of c-di-GMP synthesis, and 3 c-di-GMP-metabolizing genes, were present in Bab8416 and B. melitensis 16M, but absent in Bab9-941 and other B. abortus strains, indicating potential SPR-associated key determinants and novel molecular mechanisms. Moreover, despite almost completely intact smooth LPS related genes, only one mutated OmpA family protein of Bab8416, functionally related to flagellar and efflux pump, was newly identified. Several point mutations were identified to be Bab8416 specific while a majority of them were verified to be B. abortus ST2 characteristic. In conclusion, our study therefore identifies new SPR-associated factors that could play a role in refining and updating Brucella taxonomic schemes and provides resources for further detailed analysis of mechanism for Brucella phage resistance.

Accessory genome of the multi-drug resistant ocular isolate of Pseudomonas aeruginosa PA34

Bacteria can acquire an accessory genome through the horizontal transfer of genetic elements from non-parental lineages. This leads to rapid genetic evolution allowing traits such as antibiotic resistance and virulence to spread through bacterial communities. The study of complete genomes of bacterial strains helps to understand the genomic traits associated with virulence and antibiotic resistance. We aimed to investigate the complete accessory genome of an ocular isolate of Pseudomonas aeruginosa strain PA34. We obtained the complete genome of PA34 utilising genome sequence reads from Illumina and Oxford Nanopore Technology followed by PCR to close any identified gaps. In-depth genomic analysis was performed using various bioinformatics tools. The susceptibility to heavy metals and cytotoxicity was determined to confirm expression of certain traits. The complete genome of PA34 includes a chromosome of 6.8 Mbp and two plasmids of 95.4 Kbp (pMKPA34-1) and 26.8 Kbp (pMKPA34-2). PA34 had a large accessory genome of 1,213 genes and had 543 unique genes not present in other strains. These exclusive genes encoded features related to metal and antibiotic resistance, phage integrase and transposons. At least 24 genomic islands (GIs) were predicated in the complete chromosome, of which two were integrated into novel sites. Eleven GIs carried virulence factors or replaced pathogenic genes. A bacteriophage carried the aminoglycoside resistance gene (AAC(3)-IId). The two plasmids carried other six antibiotic resistance genes. The large accessory genome of this ocular isolate plays a large role in shaping its virulence and antibiotic resistance.

Genomic and Transcriptomic Changes that Mediate Increased Platinum Resistance in Cupriavidus metallidurans

The extensive anthropogenic use of platinum, a rare element found in low natural abundance in the Earth's continental crust and one of the critical raw materials in the EU innovation partnership framework, has resulted in increased concentrations in surface environments. To minimize its spread and increase its recovery from the environment, biological recovery via different microbial systems is explored. In contrast, studies focusing on the effects of prolonged exposure to Pt are limited. In this study, we used the metal-resistant Cupriavidus metallidurans NA4 strain to explore the adaptation of environmental bacteria to platinum exposure. We used a combined Nanopore⁻Illumina sequencing approach to fully resolve all six replicons of the C. metallidurans NA4 genome, and compared them with the C. metallidurans CH34 genome, revealing an important role in metal resistance for its chromid rather than its megaplasmids. In addition, we identified the genomic and transcriptomic changes in a laboratory-evolved strain, displaying resistance to 160 µM Pt4+. The latter carried 20 mutations, including a large 69.9 kb deletion in its plasmid pNA4_D (89.6 kb in size), and 226 differentially-expressed genes compared to its parental strain. Many membrane-related processes were affected, including up-regulation of cytochrome c and a lytic transglycosylase, down-regulation of flagellar and pili-related genes, and loss of the pNA4_D conjugative machinery, pointing towards a significant role in the adaptation to platinum.

Comparative genomics of clinical strains of Pseudomonas aeruginosa strains isolated from different geographic sites

The large and complex genome of Pseudomonas aeruginosa, which consists of significant portions (up to 20%) of transferable genetic elements contributes to the rapid development of antibiotic resistance. The whole genome sequences of 22 strains isolated from eye and cystic fibrosis patients in Australia and India between 1992 and 2007 were used to compare genomic divergence and phylogenetic relationships as well as genes for antibiotic resistance and virulence factors. Analysis of the pangenome indicated a large variation in the size of accessory genome amongst 22 stains and the size of the accessory genome correlated with number of genomic islands, insertion sequences and prophages. The strains were diverse in terms of sequence type and dissimilar to that of global epidemic P. aeruginosa clones. Of the eye isolates, 62% clustered together within a single lineage. Indian eye isolates possessed genes associated with resistance to aminoglycoside, beta-lactams, sulphonamide, quaternary ammonium compounds, tetracycline, trimethoprims and chloramphenicols. These genes were, however, absent in Australian isolates regardless of source. Overall, our results provide valuable information for understanding the genomic diversity of P. aeruginosa isolated from two different infection types and countries.

Cupriavidus metallidurans Strains with Different Mobilomes and from Distinct Environments Have Comparable Phenomes

Cupriavidus metallidurans has been mostly studied because of its resistance to numerous heavy metals and is increasingly being recovered from other environments not typified by metal contamination. They host a large and diverse mobile gene pool, next to their native megaplasmids. Here, we used comparative genomics and global metabolic comparison to assess the impact of the mobilome on growth capabilities, nutrient utilization, and sensitivity to chemicals of type strain CH34 and three isolates (NA1, NA4 and H1130). The latter were isolated from water sources aboard the International Space Station (NA1 and NA4) and from an invasive human infection (H1130). The mobilome was expanded as prophages were predicted in NA4 and H1130, and a genomic island putatively involved in abietane diterpenoids metabolism was identified in H1130. An active CRISPR-Cas system was identified in strain NA4, providing immunity to a plasmid that integrated in CH34 and NA1. No correlation between the mobilome and isolation environment was found. In addition, our comparison indicated that the metal resistance determinants and properties are conserved among these strains and thus maintained in these environments. Furthermore, all strains were highly resistant to a wide variety of chemicals, much broader than metals. Only minor differences were observed in the phenomes (measured by phenotype microarrays), despite the large difference in mobilomes and the variable (shared by two or three strains) and strain-specific genomes.

Identification and Classification for the Lactobacillus casei Group

Lactobacillus casei, Lactobacillus paracasei, and Lactobacillus rhamnosus are phenotypically and genotypically closely related, and together comprise the L. casei group. Although the strains of this group are commercially valuable as probiotics, the taxonomic status and nomenclature of the L. casei group have long been contentious because of the difficulties in identifying these three species by using the most frequently used genotypic methodology of 16S rRNA gene sequencing. Long used as the gold standard for species classification, DNA–DNA hybridization is laborious, requires expert skills, and is difficult to use routinely in laboratories. Currently, genome-based comparisons, including average nucleotide identity (ANI) and digital DNA–DNA hybridization (dDDH), are commonly applied to bacterial taxonomy as alternatives to the gold standard method for the demarcating phylogenetic relationships. To establish quick and accurate methods for identifying strains in the L. casei group at the species and subspecies levels, we developed species- and subspecies-specific identification methods based on housekeeping gene sequences and whole-cell matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) spectral pattern analysis. By phylogenetic analysis based on concatenated housekeeping gene sequences (dnaJ, dnaK, mutL, pheS, and yycH), 53 strains were separated into four clusters corresponding to the four species: L. casei, L. paracasei and L. rhamnosus, and Lactobacillus chiayiensis sp. nov. A multiplex minisequencing assay using single nucleotide polymorphism (SNP)-specific primers based on the dnaK gene sequences and species-specific primers based on the mutL gene sequences provided high resolution that enabled the strains at the species level to be identified as L. casei, L. paracasei, and L. rhamnosus. By MALDI-TOF MS analysis coupled with an internal database and ClinProTools software, species- and subspecies-level L. casei group strains were identified based on reliable scores and species- and subspecies-specific MS peaks. The L. paracasei strains were distinguished clearly at the subspecies level based on subspecies-specific MS peaks. This article describes the rapid and accurate methods used for identification and classification of strains in the L. casei group based on housekeeping gene sequences and MALDI-TOF MS analysis as well as the novel speciation of this group including L. chiayiensis sp. nov. and ‘Lactobacillus zeae’ by genome-based methods.

Unexpected diversity in the mobilome of a Pseudomonas aeruginosa strain isolated from a dental unit waterline revealed by SMRT Sequencing

The Gram-negative bacterium Pseudomonas aeruginosa is found in several habitats, both natural and human-made, and is particularly known for its recurrent presence as a pathogen in the lungs of patients suffering from cystic fibrosis, a genetic disease. Given its clinical importance, several major studies have investigated the genomic adaptation of P. aeruginosa in lungs and its transition as acute infections become chronic. However, our knowledge about the diversity and adaptation of the P. aeruginosa genome to non-clinical environments is still fragmentary, in part due to the lack of accurate reference genomes of strains from the numerous environments colonized by the bacterium. Here, we used PacBio long-read technology to sequence the genome of PPF-1, a strain of P. aeruginosa isolated from a dental unit waterline. Generating this closed genome was an opportunity to investigate genomic features that are difficult to accurately study in a draft genome (contigs state). It was possible to shed light on putative genomic islands, some shared with other reference genomes, new prophages, and the complete content of insertion sequences. In addition, four different group II introns were also found, including two characterized here and not listed in the specialized group II intron database.

Comparative Genomics of Completely Sequenced Lactobacillus helveticus Genomes Provides Insights into Strain-Specific Genes and Resolves Metagenomics Data Down to the Strain Level

Although complete genome sequences hold particular value for an accurate description of core genomes, the identification of strain-specific genes, and as the optimal basis for functional genomics studies, they are still largely underrepresented in public repositories. Based on an assessment of the genome assembly complexity for all lactobacilli, we used Pacific Biosciences' long read technology to sequence and de novo assemble the genomes of three Lactobacillus helveticus starter strains, raising the number of completely sequenced strains to 12. The first comparative genomics study for L. helveticus—to our knowledge—identified a core genome of 988 genes and sets of unique, strain-specific genes ranging from about 30 to more than 200 genes. Importantly, the comparison of MiSeq- and PacBio-based assemblies uncovered that not only accessory but also core genes can be missed in incomplete genome assemblies based on short reads. Analysis of the three genomes revealed that a large number of pseudogenes were enriched for functional Gene Ontology categories such as amino acid transmembrane transport and carbohydrate metabolism, which is in line with a reductive genome evolution in the rich natural habitat of L. helveticus. Notably, the functional Clusters of Orthologous Groups of proteins categories “cell wall/membrane biogenesis” and “defense mechanisms” were found to be enriched among the strain-specific genes. A genome mining effort uncovered examples where an experimentally observed phenotype could be linked to the underlying genotype, such as for cell envelope proteinase PrtH3 of strain FAM8627. Another possible link identified for peptidoglycan hydrolases will require further experiments. Of note, strain FAM22155 did not harbor a CRISPR/Cas system; its loss was also observed in other L. helveticus strains and lactobacillus species, thus questioning the value of the CRISPR/Cas system for diagnostic purposes. Importantly, the complete genome sequences proved to be very useful for the analysis of natural whey starter cultures with metagenomics, as a larger percentage of the sequenced reads of these complex mixtures could be unambiguously assigned down to the strain level.

Genomic Diversity in the Endosymbiotic Bacterium Rhizobium leguminosarum

Rhizobium leguminosarum bv. viciae is a soil α-proteobacterium that establishes a diazotrophic symbiosis with different legumes of the Fabeae tribe. The number of genome sequences from rhizobial strains available in public databases is constantly increasing, although complete, fully annotated genome structures from rhizobial genomes are scarce. In this work, we report and analyse the complete genome of R. leguminosarum bv. viciae UPM791. Whole genome sequencing can provide new insights into the genetic features contributing to symbiotically relevant processes such as bacterial adaptation to the rhizosphere, mechanisms for efficient competition with other bacteria, and the ability to establish a complex signalling dialogue with legumes, to enter the root without triggering plant defenses, and, ultimately, to fix nitrogen within the host. Comparison of the complete genome sequences of two strains of R. leguminosarum bv. viciae, 3841 and UPM791, highlights the existence of different symbiotic plasmids and a common core chromosome. Specific genomic traits, such as plasmid content or a distinctive regulation, define differential physiological capabilities of these endosymbionts. Among them, strain UPM791 presents unique adaptations for recycling the hydrogen generated in the nitrogen fixation process.

Approaches for in silico finishing of microbial genome sequences

The introduction of next-generation sequencing (NGS) had a significant effect on the availability of genomic information, leading to an increase in the number of sequenced genomes from a large spectrum of organisms. Unfortunately, due to the limitations implied by the short-read sequencing platforms, most of these newly sequenced genomes remained as "drafts", incomplete representations of the whole genetic content. The previous genome sequencing studies indicated that finishing a genome sequenced by NGS, even bacteria, may require additional sequencing to fill the gaps, making the entire process very expensive. As such, several in silico approaches have been developed to optimize the genome assemblies and facilitate the finishing process. The present review aims to explore some free (open source, in many cases) tools that are available to facilitate genome finishing.

From isolate to answer: how whole genome sequencing is helping us rapidly characterise nosocomial bacterial outbreaks

The occurrence of highly resistant bacterial pathogens has risen in recent years, causing immense strain on the healthcare industry. Hospital-acquired infections are arguably of most concern, as bacterial outbreaks in clinical settings provide an ideal environment for proliferation among vulnerable populations. Understanding these outbreaks beyond what can be determined with traditional clinical diagnostics and implementing these new techniques routinely in the hospital environment has now become a major focus. This brief review will discuss the three main whole genome sequence techniques available today, and how they are being used to further discriminate bacterial outbreaks in nosocomial settings.

Assessing Illumina technology for the high-throughput sequencing of bacteriophage genomes

Bacteriophages are the most abundant biological entities on the planet, playing crucial roles in the shaping of bacterial populations. Phages have smaller genomes than their bacterial hosts, yet there are currently fewer fully sequenced phage than bacterial genomes. We assessed the suitability of Illumina technology for high-throughput sequencing and subsequent assembly of phage genomes. In silico datasets reveal that 30× coverage is sufficient to correctly assemble the complete genome of ~98.5% of known phages, with experimental data confirming that the majority of phage genomes can be assembled at 30× coverage. Furthermore, in silico data demonstrate it is possible to co-sequence multiple phages from different hosts, without introducing assembly errors.

Exploration of the Drosophila buzzatii transposable element content suggests underestimation of repeats in Drosophila genomes

Background

Many new Drosophila genomes have been sequenced in recent years using new-generation sequencing platforms and assembly methods. Transposable elements (TEs), being repetitive sequences, are often misassembled, especially in the genomes sequenced with short reads. Consequently, the mobile fraction of many of the new genomes has not been analyzed in detail or compared with that of other genomes sequenced with different methods, which could shed light into the understanding of genome and TE evolution. Here we compare the TE content of three genomes: D. buzzatii st-1, j-19, and D. mojavensis.

Results

We have sequenced a new D. buzzatii genome (j-19) that complements the D. buzzatii reference genome (st-1) already published, and compared their TE contents with that of D. mojavensis. We found an underestimation of TE sequences in Drosophila genus NGS-genomes when compared to Sanger-genomes. To be able to compare genomes sequenced with different technologies, we developed a coverage-based method and applied it to the D. buzzatii st-1 and j-19 genome. Between 10.85 and 11.16 % of the D. buzzatii st-1 genome is made up of TEs, between 7 and 7,5 % of D. buzzatii j-19 genome, while TEs represent 15.35 % of the D. mojavensis genome. Helitrons are the most abundant order in the three genomes.

Conclusions

TEs in D. buzzatii are less abundant than in D. mojavensis, as expected according to the genome size and TE content positive correlation. However, TEs alone do not explain the genome size difference. TEs accumulate in the dot chromosomes and proximal regions of D. buzzatii and D. mojavensis chromosomes. We also report a significantly higher TE density in D. buzzatii and D. mojavensis X chromosomes, which is not expected under the current models. Our easy-to-use correction method allowed us to identify recently active families in D. buzzatii st-1 belonging to the LTR-retrotransposon superfamily Gypsy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2648-8) contains supplementary material, which is available to authorized users.

PacBio SMRT assembly of a complex multi-replicon genome reveals chlorocatechol degradative operon in a region of genome plasticity

We have sequenced a Burkholderia genome that contains multiple replicons and large repetitive elements that would make it inherently difficult to assemble by short read sequencing technologies. We illustrate how the integrated long read correction algorithms implemented through the PacBio Single Molecule Real-Time (SMRT) sequencing technology successfully provided a de novo assembly that is a reasonable estimate of both the gene content and genome organization without making any further modifications. This assembly is comparable to related organisms assembled by more labour intensive methods. Our assembled genome revealed regions of genome plasticity for further investigation, one of which harbours a chlorocatechol degradative operon highly homologous to those previously identified on globally ubiquitous plasmids. In an ideal world, this assembly would still require experimental validation to confirm gene order and copy number of repeated elements. However, we submit that particularly in instances where a polished genome is not the primary goal of the sequencing project, PacBio SMRT sequencing provides a financially viable option for generating a biologically relevant genome estimate that can be utilized by other researchers for comparative studies.

A Site-Specific Integrative Plasmid Found in Pseudomonas aeruginosa Clinical Isolate HS87 along with A Plasmid Carrying an Aminoglycoside-Resistant Gene

Plasmids play critical roles in bacterial fitness and evolution of Pseudomonas aeruginosa. Here two plasmids found in a drug-resistant P. aeruginosa clinical isolate HS87 were completely sequenced. The pHS87b plasmid (11.2 kb) carries phage-related genes and function-unknown genes. Notably, pHS87b encodes an integrase and has an adjacent tRNA^Thr-associated attachment site. A corresponding integrated form of pHS87b at the tRNA^Thr locus was identified on the chromosome of P. aeruginosa, showing that pHS87b is able to site-specifically integrate into the 3’-end of the tRNA^Thr gene. The pHS87a plasmid (26.8 kb) displays a plastic structure containing a putative replication module, stability factors and a variable region. The RepA of pHS87a shows significant similarity to the replication proteins of pPT23A-family plasmids. pHS87a carries a transposon Tn6049, a truncated insertion sequence ΔIS1071 and a Tn402-like class 1 integron which contains an aacA4 cassette that may confer aminoglycoside resistance. Thus, pHS87b is a site-specific integrative plasmid whereas pHS87a is a plastic antibiotic resistance plasmid. The two native plasmids may promote the fitness and evolution of P. aeruginosa.

Multiple recent horizontal transfers of a large genomic region in cheese making fungi

While the extent and impact of horizontal transfers in prokaryotes are widely acknowledged, their importance to the eukaryotic kingdom is unclear and thought by many to be anecdotal. Here we report multiple recent transfers of a huge genomic island between Penicillium spp. found in the food environment. Sequencing of the two leading filamentous fungi used in cheese making, P. roqueforti and P. camemberti, and comparison with the penicillin producer P. rubens reveals a 575 kb long genomic island in P. roqueforti—called Wallaby—present as identical fragments at non-homologous loci in P. camemberti and P. rubens. Wallaby is detected in Penicillium collections exclusively in strains from food environments. Wallaby encompasses about 250 predicted genes, some of which are probably involved in competition with microorganisms. The occurrence of multiple recent eukaryotic transfers in the food environment provides strong evidence for the importance of this understudied and probably underestimated phenomenon in eukaryotes.

Horizontal gene transfers are known to play an important role in prokaryote evolution but their impact and prevalence in eukaryotes is less clear. Here, the authors sequence the genomes of cheese making fungi P. roqueforti and P. camemberti, and provide evidence for recent horizontal transfers of a large genomic region.

The rebirth of culture in microbiology through the example of culturomics to study human gut microbiota

Bacterial culture was the first method used to describe the human microbiota, but this method is considered outdated by many researchers. Metagenomics studies have since been applied to clinical microbiology; however, a "dark matter" of prokaryotes, which corresponds to a hole in our knowledge and includes minority bacterial populations, is not elucidated by these studies. By replicating the natural environment, environmental microbiologists were the first to reduce the "great plate count anomaly," which corresponds to the difference between microscopic and culture counts. The revolution in bacterial identification also allowed rapid progress. 16S rRNA bacterial identification allowed the accurate identification of new species. Mass spectrometry allowed the high-throughput identification of rare species and the detection of new species. By using these methods and by increasing the number of culture conditions, culturomics allowed the extension of the known human gut repertoire to levels equivalent to those of pyrosequencing. Finally, taxonogenomics strategies became an emerging method for describing new species, associating the genome sequence of the bacteria systematically. We provide a comprehensive review on these topics, demonstrating that both empirical and hypothesis-driven approaches will enable a rapid increase in the identification of the human prokaryote repertoire.

Development and validation of an rDNA operon based primer walking strategy applicable to de novo bacterial genome finishing

Advances in sequencing technology have drastically increased the depth and feasibility of bacterial genome sequencing. However, little information is available that details the specific techniques and procedures employed during genome sequencing despite the large numbers of published genomes. Shotgun approaches employed by second-generation sequencing platforms has necessitated the development of robust bioinformatics tools for in silico assembly, and complete assembly is limited by the presence of repetitive DNA sequences and multi-copy operons. Typically, re-sequencing with multiple platforms and laborious, targeted Sanger sequencing are employed to finish a draft bacterial genome. Here we describe a novel strategy based on the identification and targeted sequencing of repetitive rDNA operons to expedite bacterial genome assembly and finishing. Our strategy was validated by finishing the genome of Paenibacillus polymyxa strain CR1, a bacterium with potential in sustainable agriculture and bio-based processes. An analysis of the 38 contigs contained in the P. polymyxa strain CR1 draft genome revealed 12 repetitive rDNA operons with varied intragenic and flanking regions of variable length, unanimously located at contig boundaries and within contig gaps. These highly similar but not identical rDNA operons were experimentally verified and sequenced simultaneously with multiple, specially designed primer sets. This approach also identified and corrected significant sequence rearrangement generated during the initial in silico assembly of sequencing reads. Our approach reduces the required effort associated with blind primer walking for contig assembly, increasing both the speed and feasibility of genome finishing. Our study further reinforces the notion that repetitive DNA elements are major limiting factors for genome finishing. Moreover, we provided a step-by-step workflow for genome finishing, which may guide future bacterial genome finishing projects.

A precise chloroplast genome of Nelumbo nucifera (Nelumbonaceae) evaluated with Sanger, Illumina MiSeq, and PacBio RS II sequencing platforms: insight into the plastid evolution of basal eudicots

BACKGROUND

The chloroplast genome is important for plant development and plant evolution. Nelumbo nucifera is one member of relict plants surviving from the late Cretaceous. Recently, a new sequencing platform PacBio RS II, known as 'SMRT (Single Molecule, Real-Time) sequencing', has been developed. Using the SMRT sequencing to investigate the chloroplast genome of N. nucifera will help to elucidate the plastid evolution of basal eudicots.

RESULTS

The sizes of the de novo assembled complete chloroplast genome of N. nucifera were 163,307 bp, 163,747 bp and 163,600 bp with average depths of coverage of 7×, 712× and 105× sequenced by Sanger, Illumina MiSeq and PacBio RS II, respectively. The precise chloroplast genome of N. nucifera was obtained from PacBio RS II data proofread by Illumina MiSeq reads, with a quadripartite structure containing a large single copy region (91,846 bp) and a small single copy region (19,626 bp) separated by two inverted repeat regions (26,064 bp). The genome contains 113 different genes, including four distinct rRNAs, 30 distinct tRNAs and 79 distinct peptide-coding genes. A phylogenetic analysis of 133 taxa from 56 orders indicated that Nelumbo with an age of 177 million years is a sister clade to Platanus, which belongs to the basal eudicots. Basal eudicots began to emerge during the early Jurassic with estimated divergence times at 197 million years using MCMCTree. IR expansions/contractions within the basal eudicots seem to have occurred independently.

CONCLUSIONS

Because of long reads and lack of bias in coverage of AT-rich regions, PacBio RS II showed a great promise for highly accurate 'finished' genomes, especially for a de novo assembly of genomes. N. nucifera is one member of basal eudicots, however, evolutionary analyses of IR structural variations of N. nucifera and other basal eudicots suggested that IR expansions/contractions occurred independently in these basal eudicots or were caused by independent insertions and deletions. The precise chloroplast genome of N. nucifera will present new information for structural variation of chloroplast genomes and provide new insight into the evolution of basal eudicots at the primary sequence and structural level.

The complete genome sequence of Escherichia coli EC958: a high quality reference sequence for the globally disseminated multidrug resistant E. coli O25b:H4-ST131 clone

Escherichia coli ST131 is now recognised as a leading contributor to urinary tract and bloodstream infections in both community and clinical settings. Here we present the complete, annotated genome of E. coli EC958, which was isolated from the urine of a patient presenting with a urinary tract infection in the Northwest region of England and represents the most well characterised ST131 strain. Sequencing was carried out using the Pacific Biosciences platform, which provided sufficient depth and read-length to produce a complete genome without the need for other technologies. The discovery of spurious contigs within the assembly that correspond to site-specific inversions in the tail fibre regions of prophages demonstrates the potential for this technology to reveal dynamic evolutionary mechanisms. E. coli EC958 belongs to the major subgroup of ST131 strains that produce the CTX-M-15 extended spectrum β-lactamase, are fluoroquinolone resistant and encode the fimH30 type 1 fimbrial adhesin. This subgroup includes the Indian strain NA114 and the North American strain JJ1886. A comparison of the genomes of EC958, JJ1886 and NA114 revealed that differences in the arrangement of genomic islands, prophages and other repetitive elements in the NA114 genome are not biologically relevant and are due to misassembly. The availability of a high quality uropathogenic E. coli ST131 genome provides a reference for understanding this multidrug resistant pathogen and will facilitate novel functional, comparative and clinical studies of the E. coli ST131 clonal lineage.

Elucidation of insertion elements carried on plasmids and in vitro construction of shuttle vectors from the toxic cyanobacterium Planktothrix

Several gene clusters that are responsible for toxin synthesis in bloom-forming cyanobacteria have been found to be associated with transposable elements (TEs). In particular, insertion sequence (IS) elements were shown to play a role in the inactivation or recombination of the genes responsible for cyanotoxin synthesis. Plasmids have been considered important vectors of IS element distribution to the host. In this study, we aimed to elucidate the IS elements propagated on the plasmids and the chromosome of the toxic cyanobacterium Planktothrix agardhii NIVA-CYA126/8 by means of high-throughput sequencing. In total, five plasmids (pPA5.5, pPA14, pPA50, pPA79, and pPA115, of 5, 6, 50, 79, and 120 kbp, respectively) were elucidated, and two plasmids (pPA5.5, pPA115) were found to propagate full IS element copies. Large stretches of shared DNA information between plasmids were constituted of TEs. Two plasmids (pPA5.5, pPA14) were used as candidates to engineer shuttle vectors (named pPA5.5SV and pPA14SV, respectively) in vitro by PCR amplification and the subsequent transposition of the Tn5 cat transposon containing the R6Kγ origin of replication of Escherichia coli. While pPA5.5SV was found to be fully segregated, pPA14SV consistently co-occurred with its wild-type plasmid even under the highest selective pressure. Interestingly, the Tn5 cat transposon became transferred by homologous recombination into another plasmid, pPA50. The availability of shuttle vectors is considered to be of relevance in investigating genome plasticity as a consequence of homologous recombination events. Combining the potential of high-throughput sequencing and in vitro production of shuttle vectors makes it simple to produce species-specific shuttle vectors for many cultivable prokaryotes.

Value of a newly sequenced bacterial genome

Next-generation sequencing (NGS) technologies have made high-throughput sequencing available to medium- and small-size laboratories, culminating in a tidal wave of genomic information. The quantity of sequenced bacterial genomes has not only brought excitement to the field of genomics but also heightened expectations that NGS would boost antibacterial discovery and vaccine development. Although many possible drug and vaccine targets have been discovered, the success rate of genome-based analysis has remained below expectations. Furthermore, NGS has had consequences for genome quality, resulting in an exponential increase in draft (partial data) genome deposits in public databases. If no further interests are expressed for a particular bacterial genome, it is more likely that the sequencing of its genome will be limited to a draft stage, and the painstaking tasks of completing the sequencing of its genome and annotation will not be undertaken. It is important to know what is lost when we settle for a draft genome and to determine the “scientific value” of a newly sequenced genome. This review addresses the expected impact of newly sequenced genomes on antibacterial discovery and vaccinology. Also, it discusses the factors that could be leading to the increase in the number of draft deposits and the consequent loss of relevant biological information.

A polyphasic strategy incorporating genomic data for the taxonomic description of novel bacterial species

Currently, bacterial taxonomy relies on a polyphasic approach based on the combination of phenotypic and genotypic characteristics. However, the current situation is paradoxical in that the genetic criteria that are used, including DNA–DNA hybridization, 16S rRNA gene sequence nucleotide similarity and phylogeny, and DNA G+C content, have significant limitations, but genome sequences that contain the whole genetic information of bacterial strains are not used for taxonomic purposes, despite the decreasing costs of sequencing and the increasing number of available genomes. Recently, we diversified bacterial culture conditions with the aim of isolating uncultivated bacteria. To classify the putative novel species that we cultivated, we used a polyphasic strategy that included phenotypic as well as genomic criteria (genome characteristics as well as genomic sequence similarity). Herein, we review the pros and cons of genome sequencing for taxonomy and propose that the incorporation of genome sequences in taxonomic studies has the advantage of using reliable and reproducible data. This strategy, which we name taxono-genomics, may contribute to the taxonomic classification of bacteria.

Toward a statistically explicit understanding of de novo sequence assembly

MOTIVATION

Draft de novo genome assemblies are now available for many organisms. These assemblies are point estimates of the true genome sequences. Each is a specific hypothesis, drawn from among many alternative hypotheses, of the sequence of a genome. Assembly uncertainty, the inability to distinguish between multiple alternative assembly hypotheses, can be due to real variation between copies of the genome in the sample, errors and ambiguities in the sequenced data and assumptions and heuristics of the assemblers. Most assemblers select a single assembly according to ad hoc criteria, and do not yet report and quantify the uncertainty of their outputs. Those assemblers that do report uncertainty take different approaches to describing multiple assembly hypotheses and the support for each.

RESULTS

Here we review and examine the problem of representing and measuring uncertainty in assemblies. A promising recent development is the implementation of assemblers that are built according to explicit statistical models. Some new assembly methods, for example, estimate and maximize assembly likelihood. These advances, combined with technical advances in the representation of alternative assembly hypotheses, will lead to a more complete and biologically relevant understanding of assembly uncertainty. This will in turn facilitate the interpretation of downstream analyses and tests of specific biological hypotheses.

Rapid quantification of sequence repeats to resolve the size, structure and contents of bacterial genomes

BACKGROUND

The numerous classes of repeats often impede the assembly of genome sequences from the short reads provided by new sequencing technologies. We demonstrate a simple and rapid means to ascertain the repeat structure and total size of a bacterial or archaeal genome without the need for assembly by directly analyzing the abundances of distinct k-mers among reads.

RESULTS

The sensitivity of this procedure to resolve variation within a bacterial species is demonstrated: genome sizes and repeat structure of five environmental strains of E. coli from short Illumina reads were estimated by this method, and total genome sizes corresponded well with those obtained for the same strains by pulsed-field gel electrophoresis. In addition, this approach was applied to read-sets for completed genomes and shown to be accurate over a wide range of microbial genome sizes.

CONCLUSIONS

Application of these procedures, based solely on k-mer abundances in short read data sets, allows aspects of genome structure to be resolved that are not apparent from conventional short read assemblies. This knowledge of the repetitive content of genomes provides insights into genome evolution and diversity.

Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data

We present a hierarchical genome-assembly process (HGAP) for high-quality de novo microbial genome assemblies using only a single, long-insert shotgun DNA library in conjunction with Single Molecule, Real-Time (SMRT) DNA sequencing. Our method uses the longest reads as seeds to recruit all other reads for construction of highly accurate preassembled reads through a directed acyclic graph-based consensus procedure, which we follow with assembly using off-the-shelf long-read assemblers. In contrast to hybrid approaches, HGAP does not require highly accurate raw reads for error correction. We demonstrate efficient genome assembly for several microorganisms using as few as three SMRT Cell zero-mode waveguide arrays of sequencing and for BACs using just one SMRT Cell. Long repeat regions can be successfully resolved with this workflow. We also describe a consensus algorithm that incorporates SMRT sequencing primary quality values to produce de novo genome sequence exceeding 99.999% accuracy.