Comparison of long- and short-read metagenomic assembly for low-abundance species and resistance genes.

Recent technological and computational advances have made metagenomic assembly a viable approach to achieving high-resolution views of complex microbial communities. In previous benchmarking, short-read (SR) metagenomic assemblers had the highest accuracy, long-read (LR) assemblers generated the most contiguous sequences and hybrid (HY) assemblers balanced length and accuracy. However, no assessments have specifically compared the performance of these assemblers on low-abundance species, which include clinically relevant organisms in the gut. We generated semi-synthetic LR and SR datasets by spiking small and increasing amounts of Escherichia coli isolate reads into fecal metagenomes and, using different assemblers, examined E. coli contigs and the presence of antibiotic resistance genes (ARGs). For ARG assembly, although SR assemblers recovered more ARGs with high accuracy, even at low coverages, LR assemblies allowed for the placement of ARGs within longer, E. coli-specific contigs, thus pinpointing their taxonomic origin. HY assemblies identified resistance genes with high accuracy and had lower contiguity than LR assemblies. Each assembler type's strengths were maintained even when our isolate was spiked in with a competing strain, which fragmented and reduced the accuracy of all assemblies. For strain characterization and determining gene context, LR assembly is optimal, while for base-accurate gene identification, SR assemblers outperform other options. HY assembly offers contiguity and base accuracy, but requires generating data on multiple platforms, and may suffer high misassembly rates when strain diversity exists. Our results highlight the trade-offs associated with each approach for recovering low-abundance taxa, and that the optimal approach is goal-dependent.

WeFaceNano: a user-friendly pipeline for complete ONT sequence assembly and detection of antibiotic resistance in multi-plasmid bacterial isolates.

BACKGROUND: Bacterial plasmids often carry antibiotic resistance genes and are a significant factor in the spread of antibiotic resistance. The ability to completely assemble plasmid sequences would facilitate the localization of antibiotic resistance genes, the identification of genes that promote plasmid transmission and the accurate tracking of plasmid mobility. However, the complete assembly of plasmid sequences using the currently most widely used sequencing platform (Illumina-based sequencing) is restricted due to the generation of short sequence lengths. The long-read Oxford Nanopore Technologies (ONT) sequencing platform overcomes this limitation. Still, the assembly of plasmid sequence data remains challenging due to software incompatibility with long-reads and the error rate generated using ONT sequencing. Bioinformatics pipelines have been developed for ONT-generated sequencing but require computational skills that frequently are beyond the abilities of scientific researchers. To overcome this challenge, the authors developed 'WeFaceNano', a user-friendly Web interFace for rapid assembly and analysis of plasmid DNA sequences generated using the ONT platform. WeFaceNano includes: a read statistics report; two assemblers (Miniasm and Flye); BLAST searching; the detection of antibiotic resistance- and replicon genes and several plasmid visualizations. A user-friendly interface displays the main features of WeFaceNano and gives access to the analysis tools. RESULTS: Publicly available ONT sequence data of 21 plasmids were used to validate WeFaceNano, with plasmid assemblages and anti-microbial resistance gene detection being concordant with the published results. Interestingly, the "Flye" assembler with "meta" settings generated the most complete plasmids. CONCLUSIONS: WeFaceNano is a user-friendly open-source software pipeline suitable for accurate plasmid assembly and the detection of anti-microbial resistance genes in (clinical) samples where multiple plasmids can be present.

Rapid modification of the pET-28 expression vector for ligation independent cloning using homologous recombination in Saccharomyces cerevisiae.

The ability to rapidly customize an expression vector of choice is a valuable tool for any researcher involved in high-throughput molecular cloning for protein overexpression. Unfortunately, it is common practice to amend or neglect protein targets if the gene that encodes the protein of interest is incompatible with the multiple-cloning region of a preferred expression vector. To address this issue, a method was developed to quickly exchange the multiple-cloning region of the popular expression plasmid pET-28 with a ligation-independent cloning cassette, generating pGAY-28. This cassette contains dual inverted restriction sites that reduce false positive clones by generating a linearized plasmid incapable of self-annealing after a single restriction-enzyme digest. We also establish that progressively cooling the vector and insert leads to a significant increase in ligation-independent transformation efficiency, demonstrated by the incorporation of a 10.3 kb insert into the vector. The method reported to accomplish plasmid reconstruction is uniquely versatile yet simple, relying on the strategic placement of primers combined with homologous recombination of PCR products in yeast.

Dissecting quantitative trait nucleotides by saturation genome editing.

Genome editing technologies have the potential to transform our understanding of how genetic variation gives rise to complex traits through the systematic engineering and phenotypic characterization of genetic variants. However, there has yet to be a system with sufficient efficiency, fidelity, and throughput to comprehensively identify causal variants at the genome scale. Here we explored the ability of templated CRISPR editing systems to install natural variants genome-wide in budding yeast. We optimized several approaches to enhance homology-directed repair (HDR) with donor DNA templates, including donor recruitment to target sites, single-stranded donor production by bacterial retrons, and in vivo plasmid assembly. We uncovered unique advantages of each system that we integrated into a single superior system named MAGESTIC 3.0. We used MAGESTIC 3.0 to dissect causal variants residing in 112 quantitative trait loci across 32 environmental conditions, revealing an enrichment for missense variants and loci with multiple causal variants. MAGESTIC 3.0 will facilitate the functional analysis of the genome at single-nucleotide resolution and provides a roadmap for improving template-based genome editing systems in other organisms.

Multidrug resistance and virulence genes carried by mobile genomic elements in Salmonella enterica isolated from live food animals, processed, and retail meat in North Carolina, 2018-2019.

An estimated 1000,000 domestic salmonellosis cases are attributed to food as a vehicle of exposure. Among Food Safety and Inspection Service (FSIS)-regulated products, approximately 360,000 salmonellosis cases are associated with consumption of meat, poultry, and egg products. Salmonella vaccination programs instituted in U.S. poultry, cattle, and swine production have effectively reduced the prevalence of common Salmonella enterica serotypes Typhimurium, Enteritidis, Choleraesuis (swine), and Dublin (cattle) in the past several years, with some evidence of cross-immunity to other serovars. This study investigated S. enterica (n = 741) from live food animals, meat carcasses at production, and retail meat in North Carolina collected January 2018 to December 2019. Whole-genome sequencing (WGS) and bioinformatics were used to molecularly characterize and compare AMR profiles, virulence, and phylogeny of Salmonella at three stages of food processing. Multidrug-resistant (MDR) plasmids identified also contained the integrase recombinase virulence factor int associated with mobile integrons, qacE conferred quaternary ammonia resistance, and diverse AMR profiles. MDR Plasmid IncFIB(K)_1_Kpn3_JN233704, with virulence factor int had 51 different AMR profiles within poultry S. enterica Infantis isolates. Plasmid-mediated virulence factors also appear to provide a fitness advantage, as the dominant S. enterica serotype Kentucky in chicken retail meat held the greatest diversity of plasmid-mediated colicin virulence genes which are often upregulated by environmental stressors and confer a competitive advantage. Mobile genetic element recombination is increasing pathogen fitness in the food chain through the dissemination of virulence factors and resistance genes to clinically important antibiotics, posing a significant threat to human health.

pCAT vectors overcome inefficient electroporation of Cupriavidus necator H16.

Cupriavidus necator H16 is a chemolithoautotroph with a range of industrial biotechnological applications. Advanced metabolic engineering in the bacterium, however, is impeded by low transformation efficiency, making it difficult to introduce and screen new genetic functions rapidly. This study systematically characterized the broad host range plasmids pBHR1, pBBR1MCS-2 and pKT230 used frequently for C. necator engineering. Kanamycin resistance cassette (KanR) and a truncated sequence of the replication origin (Rep) are contributing factors to C. necator low electroporation transformation efficiency. Consequently, a series of modular minimal plasmids, named pCAT, were constructed. pCAT vectors transform C. necator H16 with a > 3000-fold higher efficiency (up to 107 CFU/µg DNA) compared to control plasmids. Further, pCAT vectors are highly stable, expressing reporter proteins over several days of serial cultivation in the absence of selection pressure. Finally, they can be assembled rapidly from PCR or synthesized DNA fragments, and restriction-ligation reactions can be efficiently electroporated directly into C. necator, circumventing the requirement to use Escherichia coli for plasmid maintenance or propagation. This study demonstrates that an understanding of the behaviour of the constituent parts of plasmids in a host is key to efficient propagation of genetic information, while offering new methods for engineering a bacterium with desirable industrial biotechnological features.

Complete, closed and curated genome sequences of Photobacterium damselae subsp. piscicida isolates from Australia indicate mobilome-driven localized evolution and novel pathogenicity determinants.

Despite the recent advances in sequencing technologies, the complete assembly of multi-chromosome genomes of the Vibrionaceae, often containing several plasmids, remains challenging. Using a combination of Oxford Nanopore MinION long reads and short Illumina reads, we fully sequenced, closed and curated the genomes of two strains of a primary aquatic pathogen Photobacterium damselae subsp. piscicida isolated in Australia. These are also the first genome sequences of P. damselae subsp. piscicida isolated in Oceania and, to our knowledge, in the Southern hemisphere. We also investigated the phylogenetic relationships between Australian and overseas isolates, revealing that Australian P. damselae subsp. piscicida are more closely related to the Asian and American strains rather than to the European ones. We investigated the mobilome and present new evidence showing that a host specialization process and progressive adaptive evolution to fish are ongoing in P. damselae subsp. piscicida, and are largely mediated by transposable elements, predominantly in chromosome 2, and by plasmids. Finally, we identified two novel potential virulence determinants in P. damselae subsp. piscicida - a chorismate mutase gene, which is ubiquitously retained and co-localized with the AIP56 apoptogenic toxin-encoding gene on the pPHDP10 plasmid, and transfer-messenger RNA gene ssrA located on the main chromosome, homologous to a critical-to-virulence determinant in Yersinia pseudotuberculosis. Our study describes, to our knowledge, the only fully closed and manually curated genomes of P. damselae subsp. piscicida available to date, offering new insights into this important fish pathogen and its evolution.

An improved, versatile and efficient modular plasmid assembly system for expression analyses of genes in Xanthomonas oryzae.

Xanthomonas oryzae pathovars oryzae (Xoo) and oryzicola (Xoc) infect rice, causing bacterial blight and bacterial leaf streak, respectively, which are two economically important bacterial diseases in paddy fields. The interactions of Xoo and Xoc with rice can be used as models for studying fundamental aspects of bacterial pathogenesis and host tissue specificity. However, an improved vector system for gene expression analysis is desired for Xoo and Xoc because some broad host range vectors that can replicate stably in X. oryzae pathovars are low-copy number plasmids. To overcome this limitation, we developed a modular plasmid assembly system to transfer the functional DNA modules from the entry vectors into the pHM1-derived backbone vectors on a high-copy number basis. We demonstrated the feasibility of our vector system for protein detection, and quantification of virulence gene expression under laboratory conditions and in association with host rice and nonhost tobacco cells. This system also allows execution of a mutant complementation equivalent to the single-copy chromosomal integration system and tracing of pathogens in rice leaf. Based on this assembly system, we constructed a series of protein expression and promoter-probe vectors suitable for classical double restriction enzyme cloning. These vector systems enable cloning of all genes or promoters of interest from Xoo and Xoc strains. Our modular assembly system represents a versatile and highly efficient toolkit for gene expression analysis that will accelerate studies on interactions of X. oryzae with rice.

Multiple site-directed mutagenesis via simple cloning by prolonged overlap extension.

We describe the application of simple cloning by prolonged overlap extension for multiple site-directed mutagenesis in the same plasmid. We show that it is possible to use this technique with very short PCR templates. The technique is ideally suited for the generation of longer donor DNA sequences for CRISPR/Cas9-mediated homologous repair.

Comparison of long-read sequencing technologies in the hybrid assembly of complex bacterial genomes.

Illumina sequencing allows rapid, cheap and accurate whole genome bacterial analyses, but short reads (<300 bp) do not usually enable complete genome assembly. Long-read sequencing greatly assists with resolving complex bacterial genomes, particularly when combined with short-read Illumina data (hybrid assembly). However, it is not clear how different long-read sequencing methods affect hybrid assembly accuracy. Relative automation of the assembly process is also crucial to facilitating high-throughput complete bacterial genome reconstruction, avoiding multiple bespoke filtering and data manipulation steps. In this study, we compared hybrid assemblies for 20 bacterial isolates, including two reference strains, using Illumina sequencing and long reads from either Oxford Nanopore Technologies (ONT) or SMRT Pacific Biosciences (PacBio) sequencing platforms. We chose isolates from the family Enterobacteriaceae, as these frequently have highly plastic, repetitive genetic structures, and complete genome reconstruction for these species is relevant for a precise understanding of the epidemiology of antimicrobial resistance. We de novo assembled genomes using the hybrid assembler Unicycler and compared different read processing strategies, as well as comparing to long-read-only assembly with Flye followed by short-read polishing with Pilon. Hybrid assembly with either PacBio or ONT reads facilitated high-quality genome reconstruction, and was superior to the long-read assembly and polishing approach evaluated with respect to accuracy and completeness. Combining ONT and Illumina reads fully resolved most genomes without additional manual steps, and at a lower consumables cost per isolate in our setting. Automated hybrid assembly is a powerful tool for complete and accurate bacterial genome assembly.

Large Plasmid Complement Resolved: Complete Genome Sequencing of Lactobacillus plantarum MF1298, a Candidate Probiotic Strain Associated with Unfavorable Effect.

Considerable attention has been given to the species Lactobacillus plantarum regarding its probiotic potential. L. plantarum strains have shown health benefits in several studies, and even nonstrain-specific claims are allowed in certain markets. L. plantarum strain MF1298 was considered a candidate probiotic, demonstrating in vitro probiotic properties and the ability to survive passage through the human intestinal tract. However, the strain showed an unfavorable effect on symptoms in subjects with irritable bowel syndrome in a clinical trial. The properties and the genome of this strain are thus of general interest. Obtaining the complete genome of strain MF1298 proved difficult due to its large plasmid complement. Here, we exploit a combination of sequencing approaches to obtain the complete chromosome and plasmid assemblies of MF1298. The Oxford Nanopore Technologies MinION long-read sequencer was particularly useful in resolving the unusually large number of plasmids in the strain, 14 in total. The complete genome sequence of 3,576,440 basepairs contains 3272 protein-encoding genes, of which 315 are located on plasmids. Few unique regions were found in comparison with other L. plantarum genomes. Notably, however, one of the plasmids contains genes related to vitamin B12 (cobalamin) turnover and genes encoding bacterial reverse transcriptases, features not previously reported for L. plantarum. The extensive plasmid information will be important for future studies with this strain.

Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches.

This study aimed to assess the feasibility of using the Oxford Nanopore Technologies (ONT) MinION long-read sequencer in reconstructing fully closed plasmid sequences from eight Enterobacteriaceae isolates of six different species with plasmid populations of varying complexity. Species represented were Escherichia coli, Klebsiella pneumoniae, Citrobacter freundii, Enterobacter cloacae, Serratia marcescens and Klebsiella oxytoca, with plasmid populations ranging from 1-11 plasmids with sizes of 2-330 kb. Isolates were sequenced using Illumina (short-read) and ONT's MinION (long-read) platforms, and compared with fully resolved PacBio (long-read) sequence assemblies for the same isolates. We compared the performance of different assembly approaches including SPAdes, plasmidSPAdes, hybridSPAdes, Canu, Canu+Pilon (canuPilon) and npScarf in recovering the plasmid structures of these isolates by comparing with the gold-standard PacBio reference sequences. Overall, canuPilon provided consistently good quality assemblies both in terms of assembly statistics (N50, number of contigs) and assembly accuracy [presence of single nucleotide polymorphisms (SNPs)/indels with respect to the reference sequence]. For plasmid reconstruction, Canu recovered 70 % of the plasmids in complete contigs, and combining three assembly approaches (Canu or canuPilon, hybridSPAdes and plasmidSPAdes) resulted in a total 78 % recovery rate for all the plasmids. The analysis demonstrated the potential of using MinION sequencing technology to resolve important plasmid structures in Enterobacteriaceae species independent of and in conjunction with Illumina sequencing data. A consensus assembly derived from several assembly approaches could present significant benefit in accurately resolving the greatest number of plasmid structures.