Completing bacterial genome assemblies: strategy and performance comparisons

Full-length transcriptome characterization of Platycladus orientalis based on the PacBio platform

As a unique and native conifer in China, Platycladus orientalis is widely used in soil erosion control, garden landscapes, timber, and traditional Chinese medicine. However, due to the lack of reference genome and transcriptome, it is limited to the further molecular mechanism research and gene function mining. To develop a full-length reference transcriptome, tissues from five different parts of P. orientalis and four cone developmental stages were sequenced and analyzed by single-molecule real-time (SMRT) sequencing through the PacBio platform in this study. Overall, 37,111 isoforms were detected by PacBio with an N50 length of 2,317 nt, an average length of 1,999 bp, and the GC content of 41.81%. Meanwhile, 36,120 coding sequences, 5,645 simple sequence repeats (SSRs), 1,201 non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events with five types were identified using the results obtained from the PacBio transcript isoforms. Furthermore, 1,659 transcription factors (TFs) were detected and belonged to 51 TF families. A total of 35,689 transcripts (96.17%) were annotated through the NCBI nr, KOG, Swiss-Prot and KEGG databases, and 385 transcript isoforms related to 8 types of hormones were identified incorporated into plant hormone signal transduction pathways. The assembly and revelation of the full-length transcriptome of P. orientalis offer a pioneering insight for future investigations into gene function and genetic breeding within Platycladus species.

Comparative genomics of clinical Stenotrophomonas maltophilia isolates reveals genetic diversity which correlates with colonization and persistence in vivo

Stenotrophomonas maltophilia is a Gram-negative emerging opportunistic pathogen often present in people with respiratory diseases such as cystic fibrosis (CF). People with CF (pwCF) experience lifelong polymicrobial infections of the respiratory mucosa. Our prior work showed that Pseudomonas aeruginosa promotes persistence of S. maltophilia in mouse respiratory infections. As is typical for environmental opportunistic pathogens, S. maltophilia has a large genome and a high degree of genetic diversity. In this study, we evaluated the genomic content of S. maltophilia, combining short and long read sequencing to construct nearly complete genomes of 10 clinical isolates. The genomes of these isolates were then compared with all publicly available S. maltophilia genome assemblies, and each isolate was then evaluated for colonization/persistence in vivo, both alone and in coinfection with P. aeruginosa. We found that while the overall genome size and GC content were fairly consistent between strains, there was considerable variability in both genome structure and gene content. Similarly, there was significant variability in S. maltophilia colonization and persistence in experimental mouse respiratory infections in the presence or absence of P. aeruginosa. Ultimately, this study gives us a greater understanding of the genomic diversity of clinical S. maltophilia isolates, and how this genomic diversity relates to both interactions with other pulmonary pathogens and to host disease progression. Identifying the molecular determinants of infection with S. maltophilia can facilitate development of novel antimicrobial strategies for a highly drug-resistant pathogen.

Exploration of Bioinformatics on Microbial Fuel Cell Technology: Trends, Challenges, and Future Prospects

Microbial fuel cells (MFCs) are a cost-effective and environmentally friendly alternative energy method. MFC technology has gained much interest in recent decades owing to its effectiveness in remediating wastewater and generating bioelectricity. The microbial fuel cell generates energy mainlybecause of oxidation-reduction reactions. In this reaction, electrons were transferred between two reactants. Bioinformatics is expanding across a wide range of microbial fuel cell technology. Electroactive species in the microbial community were evaluated using bioinformatics methodologies in whole genome sequencing, RNA sequencing, transcriptomics, metagenomics, and phylogenetics. Technology advancements in microbial fuel cells primarily produce power from organic and inorganic waste from various sources. Reduced chemical oxygen demand and waste degradation are two added advantages for microbial fuel cells. From plants, bacteria, and algae, microbial fuel cells were developed. Due to the rapid advancement of sequencing techniques, bioinformatics approaches are currently widely used in the technology of microbial fuel cells. In addition, they play an important role in determining the composition of electroactive species in microorganisms. The metabolic pathway is also possible to determine with bioinformatics resources. A computational technique that reveals the nature of the mediators and the substrate was also used to predict the electrochemical properties. Computational strategies were used to tackle significant challenges in experimental procedures, such as optimization and understanding microbiological systems. The main focus of this review is on utilizing bioinformatics techniques to improve microbial fuel cell technology.

Genome Resource for Acinetobacter schindleri H4-3-C1: An Endophyte of Pseudostellaria heterophylla with Degradation Activity to Toxins Produced by Fungal Pathogens

Acinetobacter schindleri is an endophyte of Pseudostellaria heterophylla, a traditional Chinese herbal plant. It has high degradation activity to toxins produced by fungal pathogen Fusarium graminearum. Here, we deployed PacBio single-molecule real-time long-read sequencing technology to generate a complete genome assembly for the Acinetobacter schindleri H4-3-C1 strain and obtained 1.59 Gb of clean reads. These reads were assembled to a single circular DNA chromosome with a length of 3,265,024 bp, and no plasmid was found in the genome. Totals of 3,193 coding sequences, 91 transfer RNA, 21 ribosomal RNA, and 75 small RNAs were identified in the genome. This high-quality genome assembly and gene annotation resource will facilitate the excavation of the zearalenone degradation gene and provide valuable resources for preventing and controlling toxigenic fungal diseases of P. heterophylla. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Draft Genome Sequence of Clostridium butyricum Strain 16-3, Isolated from Neonatal Feces

Here, we report the genome sequence of Clostridium butyricum strain 16-3, which was isolated from infant feces. The genome contains circular contigs of 3,861,515 bp and 769,300 bp, with G+C contents of 28.8% and 28.3%, respectively.

Software Choice and Sequencing Coverage Can Impact Plastid Genome Assembly-A Case Study in the Narrow Endemic Calligonum bakuense

Most plastid genome sequences are assembled from short-read whole-genome sequencing data, yet the impact that sequencing coverage and the choice of assembly software can have on the accuracy of the resulting assemblies is poorly understood. In this study, we test the impact of both factors on plastid genome assembly in the threatened and rare endemic shrub Calligonum bakuense. We aim to characterize the differences across plastid genome assemblies generated by different assembly software tools and levels of sequencing coverage and to determine if these differences are large enough to affect the phylogenetic position inferred for C. bakuense compared to congeners. Four assembly software tools (FastPlast, GetOrganelle, IOGA, and NOVOPlasty) and seven levels of sequencing coverage across the plastid genome (original sequencing depth, 2,000x, 1,000x, 500x, 250x, 100x, and 50x) are compared in our analyses. The resulting assemblies are evaluated with regard to reproducibility, contig number, gene complement, inverted repeat length, and computation time; the impact of sequence differences on phylogenetic reconstruction is assessed. Our results show that software choice can have a considerable impact on the accuracy and reproducibility of plastid genome assembly and that GetOrganelle produces the most consistent assemblies for C. bakuense. Moreover, we demonstrate that a sequencing coverage between 500x and 100x can reduce both the sequence variability across assembly contigs and computation time. When comparing the most reliable plastid genome assemblies of C. bakuense, a sequence difference in only three nucleotide positions is detected, which is less than the difference potentially introduced through software choice.

Construction of a chromosome-level Japanese stickleback species genome using ultra-dense linkage analysis with single-cell sperm sequencing

It is still difficult to construct the genomes of higher organisms as their genome sequences must be extended to the length of the chromosome by linkage analysis. In this study, we attempted to provide an innovative alternative to conventional linkage analysis by devising a method to genotype sperm using 10× Genomics single-cell genome sequencing libraries to generate a linkage map without interbreeding individuals. A genome was assembled using sperm from the Japanese stickleback Gasterosteus nipponicus, with single-cell genotyping yielding 1 864 430 very dense hetero-SNPs and an average coverage per sperm cell of 0.13×. In total, 1665 sperm were used, which is an order of magnitude higher than the number of recombinations used for conventional linkage analysis. We then improved the linkage analysis tool scaffold extender with low depth linkage analysis (SELDLA) to analyze the data according to the characteristics of the single-cell genotyping data. Finally, we were able to determine the chromosomal location (97.1%) and orientation (64.4%) of the contigs in the 456 Mb genome of G. nipponicus, sequenced using nanopores. This method promises to be a useful tool for determining the genomes of non-model organisms for which breeding systems have not yet been established by linkage analysis.

Complete Genome Sequence of Lacticaseibacillus rhamnosus CAU 1365, Isolated from Kimchi

Here, we report the whole-genome sequence of Lacticaseibacillus rhamnosus CAU 1365, which was isolated from kimchi. The genome was composed of 1 contig with a total length of 2,991,039 bp and had 2,658 coding sequences, including 62 tRNA genes and 15 rRNA genes.

Mining the Microbiome and Microbiota-Derived Molecules in Inflammatory Bowel Disease

The intestinal microbiota is a complex community that consists of an ecosystem with a dynamic interplay between bacteria, fungi, archaea, and viruses. Recent advances in model systems have revealed that the gut microbiome is critical for maintaining homeostasis through metabolic digestive function, immune regulation, and intestinal barrier integrity. Taxonomic shifts in the intestinal microbiota are strongly correlated with a multitude of human diseases, including inflammatory bowel disease (IBD). However, many of these studies have been descriptive, and thus the understanding of the cause and effect relationship often remains unclear. Using non-human experimental model systems such as gnotobiotic mice, probiotic mono-colonization, or prebiotic supplementation, researchers have defined numerous species-level functions of the intestinal microbiota that have produced therapeutic candidates for IBD. Despite these advances, the molecular mechanisms responsible for the function of much of the microbiota and the interplay with host cellular processes remain areas of tremendous research potential. In particular, future research will need to unlock the functional molecular units of the microbiota in order to utilize this untapped resource of bioactive molecules for therapy. This review will highlight the advances and remaining challenges of microbiota-based functional studies and therapeutic discovery, specifically in IBD. One of the limiting factors for reviewing this topic is the nascent development of this area with information on some drug candidates still under early commercial development. We will also highlight the current and evolving strategies, including in the biotech industry, used for the discovery of microbiota-derived bioactive molecules in health and disease.

PacBio sequencing revealed variation in the microbiota diversity, species richness and composition between milk collected from healthy and mastitis cows

Mastitis is the economically most important disease of dairy cows. This study used PacBio single-molecule real-time sequencing technology to sequence the full-length 16S rRNAs from 27 milk samples (18 from mastitis and nine from healthy cows; the cows were at different stages of lactation). We observed that healthy or late stage milk microbiota had significantly higher microbial diversity and richness. The community composition of the microbiota of different groups also varied greatly. The healthy cow milk microbiota was predominantly comprised of Lactococcus lactis, Acinetobacter johnsonii, and Bacteroides dorei, while the milk from mastitis cows was predominantly comprised of Bacillus cereus. The prevalence of L. lactis and B. cereus in the milk samples was confirmed by digital droplets PCR. Differences in the milk microbiota diversity and composition could suggest an important role for some these microbes in protecting the host from mastitis while others associated with mastitis. The results of our research serve as useful references for designing strategies to prevent and treat mastitis.

De novo genome assembly of Bacillus altitudinis 19RS3 and Bacillus altitudinis T5S-T4, two plant growth-promoting bacteria isolated from Ilex paraguariensis St. Hil. (yerba mate)

Plant growth-promoting bacteria (PGPB) are a heterogeneous group of bacteria that can exert beneficial effects on plant growth directly or indirectly by different mechanisms. PGPB-based inoculant formulation has been used to replace chemical fertilizers and pesticides. In our previous studies, two endophytic endospore-forming bacteria identified as Bacillus altitudinis were isolated from roots of Ilex paraguariensis St. Hil. seedlings and selected for their plant growth-promoting (PGP) properties shown in vitro and in vivo. The purposes of this work were to assemble the genomes of B. altitudinis 19RS3 and T5S-T4, using different assemblers available for Windows and Linux and to select the best assembly for each strain. Both genomes were also automatically annotated to detect PGP genes and compare sequences with other genomes reported. Library construction and draft genome sequencing were performed by Macrogen services. Raw reads were filtered using the Trimmomatic tool. Genomes were assembled using SPAdes, ABySS, Velvet, and SOAPdenovo2 assemblers for Linux, and Geneious and CLC Genomics Workbench assemblers for Windows. Assembly evaluation was done by the QUAST tool. The parameters evaluated were the number of contigs ≥ 500 bp and ≥ 1000 bp, the length of the longest contig, and the N50 value. For genome annotation PROKKA, RAST, and KAAS tools were used. The best assembly for both genomes was obtained using Velvet. The B. altitudinis 19RS3 genome was assembled into 15 contigs with an N50 value of 1,943,801 bp. The B. altitudinis T5S-T4 genome was assembled into 24 contigs with an N50 of 344,151 bp. Both genomes comprise several genes related to PGP mechanisms, such as those for nitrogen fixation, iron metabolism, phosphate metabolism, and auxin biosynthesis. The results obtained offer the basis for a better understanding of B. altitudinis 19RS3 and T5S-T4 and make them promissory for bioinoculant development.

Description of Paenibacillus dokdonensis sp. nov., a new bacterium isolated from soil

Two strains isolated from soil samples were designated as YH-JAE5^T and YH-JAE2. The isolates were facultative anaerobic, Gram-stain-variable, motile, rod-shaped bacteria. Phylogenetic analysis indicated that the isolates belonged to the genus Paenibacillus, but the 16S rRNA gene sequence similarities were <98 % when compared with other species within the genus. Analysis of rpoB gene revealed the isolates formed a sub-cluster with P. chibensis. The only menaquinone identified was MK-7. The two isolates contained meso-diaminopimelic acid within their cell wall peptidoglycan. The major polar lipids were diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, phospholipid, aminophospholipids, and lipids. The major fatty acids were C_15 : 0 anteiso and C_15 : 0 iso. The average nucleotide identity, average amino acid identity, and digital DNA-DNA hybridization values between isolate YH-JAE5^T and the most closely related reference strain (Paenibacillus chibensis KCTC 3758^T) were 81.7, 84.8 and 23.4 %, respectively. The G+C content of the genomic DNA was 47.4 mol%. Thus, the polyphasic data revealed that YH-JAE2 (=KCTC 43239=JCM 34435) and YH-JAE5^T (=KCTC 43059=JCM 33533) represent a new species. The name Paenibacillus dokdonensis sp. nov. is proposed.

Streamlining data-intensive biology with workflow systems

As the scale of biological data generation has increased, the bottleneck of research has shifted from data generation to analysis. Researchers commonly need to build computational workflows that include multiple analytic tools and require incremental development as experimental insights demand tool and parameter modifications. These workflows can produce hundreds to thousands of intermediate files and results that must be integrated for biological insight. Data-centric workflow systems that internally manage computational resources, software, and conditional execution of analysis steps are reshaping the landscape of biological data analysis and empowering researchers to conduct reproducible analyses at scale. Adoption of these tools can facilitate and expedite robust data analysis, but knowledge of these techniques is still lacking. Here, we provide a series of strategies for leveraging workflow systems with structured project, data, and resource management to streamline large-scale biological analysis. We present these practices in the context of high-throughput sequencing data analysis, but the principles are broadly applicable to biologists working beyond this field.

Nonomuraea nitratireducens sp. nov., a new actinobacterium isolated from Suaeda australis Moq. rhizosphere

A novel actinomycete, designated WYY166T, was isolated from the rhizosphere of Suaeda australis Moq. collected in Dongfang, PR China. The taxonomic position of this strain was investigated using a polyphasic approach. Phylogenetic analysis based on its 16S rRNA gene referred strain WYY166T to the genus Nonomuraea , and it was most closely related to the type strains Nonomuraea candida HMC10T, Nonomuraea turkmeniaca DSM 43926T, Nonomuraea maritima NBRC 106687T and Nonomuraea polychroma DSM 43925T (98.35, 97.60, 97.36 and 97.30% sequence similarity, respectively). Genome sequencing revealed a genome size of 11.27 Mbp and a G+C content of 71.10 mol%. The genome average nucleotide identity (ANI) values and the digital DNA - DNA hybridization (dDDH) values between strain WYY166T and the other species of the genus were found to be low (ANI 81.63~85.23 %, dDDH 23.6~31.6 %), suggesting that it represented a new species. The physiological evaluation showed that it had remarkable nitrate reduction activity. The whole-cell hydrolysates contained meso-diaminopimelic acid and madurose. The N-acyl type of muramic acid was acetyl. The major menaquinones were MK-9 (H4) (86.9 %) and MK-9 (H2) (13.1 %). The predominant fatty acids were iso-C16 : 0 (53.2 %), 10-methyl C17 : 0 (10.7 %), C17 : 1  ω6c (8.3 %) and iso-C16 : 1 h (7.3 %). These physiological, biochemical and chemotaxonomic data suggested that strain WYY166T should be classified as representing a novel species of the genus Nonomuraea , for which the name Nonomuraea nitratireducens sp. nov. is proposed. The type strain is WYY166T (=MCCC 1K03779T=KCTC 49343T).

Complete Genome Sequence of Weissella cibaria Strain BM2, Isolated from Korean Kimchi

Weissella cibaria appears to have broad-spectrum health benefits. Here, we report the genome sequence of Weissella cibaria strain BM2, which was isolated from homemade kimchi; it consists of one circular chromosome of 2,462,443 bp and one plasmid of 11,067 bp. A total of 2,337 coding sequences were predicted, including 2,117 protein-coding sequences and a G+C content of 45.06%.

Weissella cibaria appears to have broad-spectrum health benefits. Here, we report the genome sequence of Weissella cibaria strain BM2, which was isolated from homemade kimchi; it consists of one circular chromosome of 2,462,443 bp and one plasmid of 11,067 bp. A total of 2,337 coding sequences were predicted, including 2,117 protein-coding sequences and a G+C content of 45.06%.

Complete Genome Sequence of Bacillus velezensis Strain DKU_NT_04, Isolated from a Traditional Korean Food Made from Soybeans (Cheonggukjang)

In the present work, we report the complete genome sequence of Bacillus velezensis DKU_NT_04, isolated from cheonggukjang, which is a traditional Korean fermented soybean paste. The final genome assembly consists of a 4.328-Mbp chromosome with 4,134 coding sequences and a G+C content of 45.21%.

Comprehensive evaluation of non-hybrid genome assembly tools for third-generation PacBio long-read sequence data

Long reads obtained from third-generation sequencing platforms can help overcome the long-standing challenge of the de novo assembly of sequences for the genomic analysis of non-model eukaryotic organisms. Numerous long-read-aided de novo assemblies have been published recently, which exhibited superior quality of the assembled genomes in comparison with those achieved using earlier second-generation sequencing technologies. Evaluating assemblies is important in guiding the appropriate choice for specific research needs. In this study, we evaluated 10 long-read assemblers using a variety of metrics on Pacific Biosciences (PacBio) data sets from different taxonomic categories with considerable differences in genome size. The results allowed us to narrow down the list to a few assemblers that can be effectively applied to eukaryotic assembly projects. Moreover, we highlight how best to use limited genomic resources for effectively evaluating the genome assemblies of non-model organisms.

Plasmid Reconstruction from Next-Gen Data: A Detailed Protocol for the Use of PLACNETw for the Reconstruction of Plasmids from WGS Datasets

Mobile Genetic Elements (MGE) play essential roles in adaptive bacterial evolution, facilitating genetic exchange for extrachromosomal DNA, especially antibiotic resistance genes and virulence factors. For this reason, high-throughput next-generation sequencing of bacteria is of great relevance, especially for clinical pathogenic bacteria. Accurate identification of MGE from whole-genome sequencing (WGS) datasets is one of the major challenges, still hindered by methodological limitations and high sequencing costs.This chapter encompasses the protocol used for plasmid reconstruction by applying the PLACNETw methodology, from raw reads to assembled plasmids and chromosome. PLACNETw is a graphical user-friendly interface to visualize and reconstruct MGE from short-read WGS datasets. No bioinformatic background or sophisticated computational resources are required and high precision and sensitivity are achieved.

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.

Phylogenomics of colistin-susceptible and resistant XDR Acinetobacter baumannii

Background

Acinetobacter baumannii is a healthcare-associated pathogen with high rates of carbapenem resistance. Colistin is now routinely used for treatment of infections by this pathogen. However, colistin use has been associated with development of resistance to this agent.

Objectives

To elucidate the phylogenomics of colistin-susceptible and -resistant A. baumannii strain pairs from a cohort of hospitalized patients at a tertiary medical centre in the USA.

Methods

WGS data from 21 pairs of colistin-susceptible and -resistant, XDR clinical strains were obtained and compared using phylogeny of aligned genome sequences, assessment of pairwise SNP differences and gene content.

Results

Fourteen patients had colistin-resistant strains that were highly genetically related to their own original susceptible strain with a median pairwise SNP distance of 5.5 (range 1-40 SNPs), while seven other strain pairs were divergent with ≥84 SNP differences. In addition, several strains from different patients formed distinct clusters on the phylogeny in keeping with closely linked transmission chains. The majority of colistin-resistant strains contained non-synonymous mutations within the pmrAB locus suggesting a central role for pmrAB mutations in colistin resistance. Excellent genotype-phenotype correlation was also observed for carbapenems, aminoglycosides and tetracyclines.

Conclusions

The findings suggest that colistin resistance in the clinical setting arises through both in vivo evolution from colistin-susceptible strains and reinfection by unrelated colistin-resistant strains, the latter of which may involve patient-to-patient transmission.

Completing Circular Bacterial Genomes With Assembly Complexity by Using a Sampling Strategy From a Single MinION Run With Barcoding

The Oxford Nanopore MinION is an affordable and portable DNA sequencer that can produce very long reads (tens of kilobase pairs), which enable de novo bacterial genome assembly. Although many algorithms and tools have been developed for base calling, read mapping, de novo assembly, and polishing, an automated pipeline is not available for one-stop analysis for circular bacterial genome reconstruction. In this paper, we present the pipeline CCBGpipe for completing circular bacterial genomes. Raw current signals are demultiplexed and base called to generate sequencing data. Sequencing reads are de novo assembled several times by using a sampling strategy to produce circular contigs that have a sequence in common between their start and end. The circular contigs are polished by using raw signals and sequencing reads; then, duplicated sequences are removed to form a linear representation of circular sequences. The circularized contigs are finally rearranged to start at the start position of dnaA/repA or a replication origin based on the GC skew. CCBGpipe implemented in Python is available at https://github.com/jade-nhri/CCBGpipe. Using sequencing data produced from a single MinION run, we obtained 48 circular sequences, comprising 12 chromosomes and 36 plasmids of 12 bacteria, including Acinetobacter nosocomialis, Acinetobacter pittii, and Staphylococcus aureus. With adequate quantities of sequencing reads (80×), CCBGpipe can provide a complete and automated assembly of circular bacterial genomes.

Genomic and Metagenomic Approaches for Predictive Surveillance of Emerging Pathogens and Antibiotic Resistance

Antibiotic-resistant organisms (AROs) are a major concern to public health worldwide. While antibiotics have been naturally produced by environmental bacteria for millions of years, modern widespread use of antibiotics has enriched resistance mechanisms in human-impacted bacterial environments. Antibiotic resistance genes (ARGs) continue to emerge and spread rapidly. To combat the global threat of antibiotic resistance, researchers must develop methods to rapidly characterize AROs and ARGs, monitor their spread across space and time, and identify novel ARGs and resistance pathways. We review how high-throughput sequencing-based methods can be combined with classic culture-based assays to characterize, monitor, and track AROs and ARGs. Then, we evaluate genomic and metagenomic methods for identifying ARGs and biosynthetic pathways for novel antibiotics from genomic data sets. Together, these genomic analyses can improve surveillance and prediction of emerging resistance threats and accelerate the development of new antibiotic therapies to combat resistance.

Sequencing-based methods and resources to study antimicrobial resistance

Antimicrobial resistance extracts high morbidity, mortality and economic costs yearly by rendering bacteria immune to antibiotics. Identifying and understanding antimicrobial resistance are imperative for clinical practice to treat resistant infections and for public health efforts to limit the spread of resistance. Technologies such as next-generation sequencing are expanding our abilities to detect and study antimicrobial resistance. This Review provides a detailed overview of antimicrobial resistance identification and characterization methods, from traditional antimicrobial susceptibility testing to recent deep-learning methods. We focus on sequencing-based resistance discovery and discuss tools and databases used in antimicrobial resistance studies.

TSD: A Computational Tool To Study the Complex Structural Variants Using PacBio Targeted Sequencing Data

PacBio sequencing is a powerful approach to study DNA or RNA sequences in a longer scope. It is especially useful in exploring the complex structural variants generated by random integration or multiple rearrangement of endogenous or exogenous sequences. Here, we present a tool, TSD, for complex structural variant discovery using PacBio targeted sequencing data. It allows researchers to identify and visualize the genomic structures of targeted sequences by unlimited splitting, alignment and assembly of long PacBio reads. Application to the sequencing data derived from an HBV integrated human cell line(PLC/PRF/5) indicated that TSD could recover the full profile of HBV integration events, especially for the regions with the complex human-HBV genome integrations and multiple HBV rearrangements. Compared to other long read analysis tools, TSD showed a better performance for detecting complex genomic structural variants. TSD is publicly available at: https://github.com/menggf/tsd.

The present and future of de novo whole-genome assembly

As the advent of next-generation sequencing (NGS) technology, various de novo assembly algorithms based on the de Bruijn graph have been developed to construct chromosome-level sequences. However, numerous technical or computational challenges in de novo assembly still remain, although many bright ideas and heuristics have been suggested to tackle the challenges in both experimental and computational settings. In this review, we categorize de novo assemblers on the basis of the type of de Bruijn graphs (Hamiltonian and Eulerian) and discuss the challenges of de novo assembly for short NGS reads regarding computational complexity and assembly ambiguity. Then, we discuss how the limitations of the short reads can be overcome by using a single-molecule sequencing platform that generates long reads of up to several kilobases. In fact, the long read assembly has caused a paradigm shift in whole-genome assembly in terms of algorithms and supporting steps. We also summarize (i) hybrid assemblies using both short and long reads and (ii) overlap-based assemblies for long reads and discuss their challenges and future prospects. This review provides guidelines to determine the optimal approach for a given input data type, computational budget or genome.

Complete sequence of kenaf (Hibiscus cannabinus) mitochondrial genome and comparative analysis with the mitochondrial genomes of other plants

Plant mitochondrial (mt) genomes are species specific due to the vast of foreign DNA migration and frequent recombination of repeated sequences. Sequencing of the mt genome of kenaf (Hibiscus cannabinus) is essential for elucidating its evolutionary characteristics. In the present study, single-molecule real-time sequencing technology (SMRT) was used to sequence the complete mt genome of kenaf. Results showed that the complete kenaf mt genome was 569,915 bp long and consisted of 62 genes, including 36 protein-coding, 3 rRNA and 23 tRNA genes. Twenty-five introns were found among nine of the 36 protein-coding genes, and five introns were trans-spliced. A comparative analysis with other plant mt genomes showed that four syntenic gene clusters were conserved in all plant mtDNAs. Fifteen chloroplast-derived fragments were strongly associated with mt genes, including the intact sequences of the chloroplast genes psaA, ndhB and rps7. According to the plant mt genome evolution analysis, some ribosomal protein genes and succinate dehydrogenase genes were frequently lost during the evolution of angiosperms. Our data suggest that the kenaf mt genome retained evolutionarily conserved characteristics. Overall, the complete sequencing of the kenaf mt genome provides additional information and enhances our better understanding of mt genomic evolution across angiosperms.

Complete Genome Sequence of Bacillus subtilis Strain DKU_NT_02, Isolated from Traditional Korean Food Using Soybean (Chung-gook-jang) for High-Quality Poly-γ-Glutamic Acid Activity

The complete genome sequence of Bacillus subtilis strain DKU_NT_02, isolated from traditional Korean food using soybeans (chung-gook-jang), is presented here. This strain was chosen to help identify genetic factors with high-quality poly-γ-glutamic acid (γPGA) activity.

The complete genome sequence of Bacillus subtilis strain DKU_NT_02, isolated from traditional Korean food using soybeans (chung-gook-jang), is presented here. This strain was chosen to help identify genetic factors with high-quality poly-γ-glutamic acid (γPGA) activity.

Complete Genome Sequence of Bacillus subtilis Strain DKU_NT_03, Isolated from a Traditional Korean Food Using Soybean (Chung-gook-jang) for High-Quality Nattokinase Activity

We present here the complete genome sequence of Bacillus subtilis strain DKU_NT_03 isolated from the traditional Korean food chung-gook-jang, which is made from soybeans. This strain was chosen to identify genetic factors with high-quality nattokinase activity.

Recent Advances in Molecular Diagnosis of Pseudomonasaeruginosa Infection by State-of-the-Art Genotyping Techniques

Pseudomonas aeruginosa is a rod-shaped Gram-negative bacterium which is notably known as a pathogen in humans, animals, and plants. Infections caused by P. aeruginosa especially in hospitalized patients are often life-threatening and rapidly increasing worldwide throughout the years. Recently, multidrug-resistant P. aeruginosa has taken a toll on humans' health due to the inefficiency of antimicrobial agents. Therefore, the rapid and advanced diagnostic techniques to accurately detect this bacterium particularly in clinical samples are indeed necessary to ensure timely and effective treatments and to prevent outbreaks. This review aims to discuss most recent of state-of-the-art molecular diagnostic techniques enabling fast and accurate detection and identification of P. aeruginosa based on well-developed genotyping techniques, e.g., polymerase chain reaction, pulse-field gel electrophoresis, and next generation sequencing. The advantages and limitations of each of the methods are also reviewed.

Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes

Advancement of DNA sequencing technology allows the routine use of genome sequences in the various fields of microbiology. The information held in genome sequences proved to provide objective and reliable means in the taxonomy of prokaryotes. Here, we describe the minimal standards for the quality of genome sequences and how they can be applied for taxonomic purposes.

Genome characterization of two bile-isolated Vibrio fluvialis strains: an insight into pathogenicity and bile salt adaption

Vibrio fluvialis is recognized as an emerging pathogen. However, not much is known about the mechanism of its pathogenesis, and its adaptation to a special niche such as the gall bladder. Here we describe two V. fluvialis strains that cause acute cholecystitis. It is noteworthy that both strains were susceptible to all antibiotics tested, which is in contrast to previous studies, suggesting substantial genetic diversity among V. fluvialis isolates. In agreement with their survival and growth in the gall bladder, the genomes of strains 12605 and 3663 contain a considerable number of genes that confer resistance to bile, including toxR, ompU, tolC, cmeABC, rlpB, yrbK, rpoS, damX and gltK. Furthermore, integrative and conjugative elements (ICEs), virulence factors and prophage regions were also detected in strains 12605 and 3663, reflecting their flexibility in recombination during the evolution of pathogenicity. Comparative analysis of nine available genomes of V. fluvialis revealed a core genome consisting of 3,147 genes. Our results highlight the association of V. fluvialis with a rare disease profile and shed light on the evolution of pathogenesis and niche adaptation of V. fluvialis.

PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Although PacBio third-generation sequencers have improved the read lengths of genome sequencing which facilitates the assembly of complete genomes, no study has reported success in using PacBio data alone to completely sequence a two-chromosome bacterial genome from a single library in a single run. Previous studies using earlier versions of sequencing chemistries have at most been able to finish bacterial genomes containing only one chromosome with de novo assembly. In this study, we compared the robustness of PacBio RS II, using one SMRT cell and the latest P6-C4 chemistry, with Illumina HiSeq 1500 in sequencing the genome of Burkholderia pseudomallei, a bacterium which contains two large circular chromosomes, very high G+C content of 68–69%, highly repetitive regions and substantial genomic diversity, and represents one of the largest and most complex bacterial genomes sequenced, using a reference genome generated by hybrid assembly using PacBio and Illumina datasets with subsequent manual validation. Results showed that PacBio data with de novo assembly, but not Illumina, was able to completely sequence the B. pseudomallei genome without any gaps or mis-assemblies. The two large contigs of the PacBio assembly aligned unambiguously to the reference genome, sharing >99.9% nucleotide identities. Conversely, Illumina data assembled using three different assemblers resulted in fragmented assemblies (201–366 contigs), sharing only 92.2–100% and 92.0–100% nucleotide identities to chromosomes I and II reference sequences, respectively, with no indication that the B. pseudomallei genome consisted of two chromosomes with four copies of ribosomal operons. Among all assemblies, the PacBio assembly recovered the highest number of core and virulence proteins, and housekeeping genes based on whole-genome multilocus sequence typing (wgMLST). Most notably, assembly solely based on PacBio outperformed even hybrid assembly using both PacBio and Illumina datasets. Hybrid approach generated only 74 contigs, while the PacBio data alone with de novo assembly achieved complete closure of the two-chromosome B. pseudomallei genome without additional costly bench work and further sequencing. PacBio RS II using P6-C4 chemistry is highly robust and cost-effective and should be the platform of choice in sequencing bacterial genomes, particularly for those that are well-known to be difficult-to-sequence.

First report of two complete Clostridium chauvoei genome sequences and detailed in silico genome analysis

Clostridium (C.) chauvoei is a Gram-positive, spore forming, anaerobic bacterium. It causes black leg in ruminants, a typically fatal histotoxic myonecrosis. High quality circular genome sequences were generated for the C. chauvoei type strain DSM 7528^T (ATCC 10092^T) and a field strain 12S0467 isolated in Germany. The origin of replication (oriC) was comparable to that of Bacillus subtilis in structure with two regions containing DnaA boxes. Similar prophages were identified in the genomes of both C. chauvoei strains which also harbored hemolysin and bacterial spore formation genes. A CRISPR type I-B system with limited variations in the repeat number was identified. Sporulation and germination process related genes were homologous to that of the Clostridia cluster I group but novel variations for regulatory genes were identified indicative for strain specific control of regulatory events. Phylogenomics showed a higher relatedness to C. septicum than to other so far sequenced genomes of species belonging to the genus Clostridium. Comparative genome analysis of three C. chauvoei circular genome sequences revealed the presence of few inversions and translocations in locally collinear blocks (LCBs). The species genome also shows a large number of genes involved in proteolysis, genes for glycosyl hydrolases and metal iron transportation genes which are presumably involved in virulence and survival in the host. Three conserved flagellar genes (fliC) were identified in each of the circular genomes. In conclusion this is the first comparative analysis of circular genomes for the species C. chauvoei, enabling insights into genome composition and virulence factor variation.

Comparative genome analysis of fish pathogen Flavobacterium columnare reveals extensive sequence diversity within the species

Flavobacterium columnare is one of the deadliest fish pathogens causing devastating mortality in various freshwater fish species globally. To gain an insight into bacterial genomic contents and structures, comparative genome analyses were performed using the reference and newly sequenced genomes of F. columnare including genomovar I, II and I/II strains isolated from Thailand, Europe and the USA. Bacterial genomes varied in size from 3.09 to 3.39Mb (2714 to 3101 CDSs). The pan-genome analysis revealed open pan-genome nature of F. columnare strains, which possessed at least 4953 genes and tended to increase progressively with the addition of a new genome. Genomic islands (GIs) present in bacterial genomes were diverse, in which 65% (39 out of 60) of possible GIs were strain-specific. A CRISPR/cas investigation indicated at least two different CRISPR systems with varied spacer profiles. On the other hand, putative virulence genes, including those related to gliding motility, type IX secretion system (T9SS), outer membrane proteins (Omp), were equally distributed among F. columnare strains. The MLSA scheme categorized bacterial strains into nine different sequence types (ST 9-17). Phylogenetic analyses based on either 16S rRNA, MLSA and concatenated SNPs of core genome revealed the diversity of F. columnare strains. DNA homology analysis indicated that the estimated digital DNA-DNA hybridization (dDDH) between strains of genomovar I and II can be as low as 42.6%, while the three uniquely tilapia-originated strains from Thailand (1214, NK01 and 1215) were clearly dissimilar to other F. columnare strains as the dDDH values were only 27.7-30.4%. Collectively, this extensive diversity among bacterial strains suggested that species designation of F. columnare would potentially require re-emendation.

Comparative genomics of Burkholderia multivorans, a ubiquitous pathogen with a highly conserved genomic structure

The natural environment serves as a reservoir of opportunistic pathogens. A well-established method for studying the epidemiology of such opportunists is multilocus sequence typing, which in many cases has defined strains predisposed to causing infection. Burkholderia multivorans is an important pathogen in people with cystic fibrosis (CF) and its epidemiology suggests that strains are acquired from non-human sources such as the natural environment. This raises the central question of whether the isolation source (CF or environment) or the multilocus sequence type (ST) of B. multivorans better predicts their genomic content and functionality. We identified four pairs of B. multivorans isolates, representing distinct STs and consisting of one CF and one environmental isolate each. All genomes were sequenced using the PacBio SMRT sequencing technology, which resulted in eight high-quality B. multivorans genome assemblies. The present study demonstrated that the genomic structure of the examined B. multivorans STs is highly conserved and that the B. multivorans genomic lineages are defined by their ST. Orthologous protein families were not uniformly distributed among chromosomes, with core orthologs being enriched on the primary chromosome and ST-specific orthologs being enriched on the second and third chromosome. The ST-specific orthologs were enriched in genes involved in defense mechanisms and secondary metabolism, corroborating the strain-specificity of these virulence characteristics. Finally, the same B. multivorans genomic lineages occur in both CF and environmental samples and on different continents, demonstrating their ubiquity and evolutionary persistence.

A brief primer on genomic epidemiology: lessons learned from Mycobacterium tuberculosis

Genomics is now firmly established as a technique for the investigation and reconstruction of communicable disease outbreaks, with many genomic epidemiology studies focusing on revealing transmission routes of Mycobacterium tuberculosis. In this primer, we introduce the basic techniques underlying transmission inference from genomic data, using illustrative examples from M. tuberculosis and other pathogens routinely sequenced by public health agencies. We describe the laboratory and epidemiological scenarios under which genomics may or may not be used, provide an introduction to sequencing technologies and bioinformatics approaches to identifying transmission-informative variation and resistance-associated mutations, and discuss how variation must be considered in the light of available clinical and epidemiological information to infer transmission.

BAC-pool sequencing and analysis confirms growth-associated QTLs in the Asian seabass genome

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

Navigating Microbiological Food Safety in the Era of Whole-Genome Sequencing

The epidemiological investigation of a foodborne outbreak, including identification of related cases, source attribution, and development of intervention strategies, relies heavily on the ability to subtype the etiological agent at a high enough resolution to differentiate related from nonrelated cases. Historically, several different molecular subtyping methods have been used for this purpose; however, emerging techniques, such as single nucleotide polymorphism (SNP)-based techniques, that use whole-genome sequencing (WGS) offer a resolution that was previously not possible. With WGS, unlike traditional subtyping methods that lack complete information, data can be used to elucidate phylogenetic relationships and disease-causing lineages can be tracked and monitored over time. The subtyping resolution and evolutionary context provided by WGS data allow investigators to connect related illnesses that would be missed by traditional techniques. The added advantage of data generated by WGS is that these data can also be used for secondary analyses, such as virulence gene detection, antibiotic resistance gene profiling, synteny comparisons, mobile genetic element identification, and geographic attribution. In addition, several software packages are now available to generate in silico results for traditional molecular subtyping methods from the whole-genome sequence, allowing for efficient comparison with historical databases. Metagenomic approaches using next-generation sequencing have also been successful in the detection of nonculturable foodborne pathogens. This review addresses state-of-the-art techniques in microbial WGS and analysis and then discusses how this technology can be used to help support food safety investigations. Retrospective outbreak investigations using WGS are presented to provide organism-specific examples of the benefits, and challenges, associated with WGS in comparison to traditional molecular subtyping techniques.

Oxford Nanopore MinION Sequencing and Genome Assembly

The revolution of genome sequencing is continuing after the successful second-generation sequencing (SGS) technology. The third-generation sequencing (TGS) technology, led by Pacific Biosciences (PacBio), is progressing rapidly, moving from a technology once only capable of providing data for small genome analysis, or for performing targeted screening, to one that promises high quality de novo assembly and structural variation detection for human-sized genomes. In 2014, the MinION, the first commercial sequencer using nanopore technology, was released by Oxford Nanopore Technologies (ONT). MinION identifies DNA bases by measuring the changes in electrical conductivity generated as DNA strands pass through a biological pore. Its portability, affordability, and speed in data production makes it suitable for real-time applications, the release of the long read sequencer MinION has thus generated much excitement and interest in the genomics community. While de novo genome assemblies can be cheaply produced from SGS data, assembly continuity is often relatively poor, due to the limited ability of short reads to handle long repeats. Assembly quality can be greatly improved by using TGS long reads, since repetitive regions can be easily expanded into using longer sequencing lengths, despite having higher error rates at the base level. The potential of nanopore sequencing has been demonstrated by various studies in genome surveillance at locations where rapid and reliable sequencing is needed, but where resources are limited.

Third generation sequencing technologies applied to diagnostic microbiology: benefits and challenges in applications and data analysis

INTRODUCTION

The diagnosis of infectious diseases is among the most successful areas of application of new generation sequencing technologies. The field has seen the development of numerous experimental and analytical approaches for the detection and the fine description of pathogenic and non-pathogenic microorganisms.

AREAS COVERED

Without claiming to be exhaustive with respect to all applications and methods developed over the years, this review focuses on the advantages and the issues brought by the new technologies, with an eye in particular to third generation sequencing methods. Both experimental procedures and algorithmic strategies are presented, following the most relevant publications which have led to progress in our ability of detecting infectious agents. Expert commentary: The technical advance brought by third generation sequencing platforms has the potential to significantly expand the range of diagnostic tools that will be available to clinicians. Nonetheless, the implementation of these technologies in clinical practice is still far from being actionable and will temporally follow the path undertaken by second generation methods, which still require the setup of standardized pipelines in both wet and dry laboratory procedures.

Comparative genomic and functional analyses: unearthing the diversity and specificity of nematicidal factors in Pseudomonas putida strain 1A00316

We isolated Pseudomonas putida (P. putida) strain 1A00316 from Antarctica. This bacterium has a high efficiency against Meloidogyne incognita (M. incognita) in vitro and under greenhouse conditions. The complete genome of P. putida 1A00316 was sequenced using PacBio single molecule real-time (SMRT) technology. A comparative genomic analysis of 16 Pseudomonas strains revealed that although P. putida 1A00316 belonged to P. putida, it was phenotypically more similar to nematicidal Pseudomonas fluorescens (P. fluorescens) strains. We characterized the diversity and specificity of nematicidal factors in P. putida 1A00316 with comparative genomics and functional analysis, and found that P. putida 1A00316 has diverse nematicidal factors including protein alkaline metalloproteinase AprA and two secondary metabolites, hydrogen cyanide and cyclo-(l-isoleucyl-l-proline). We show for the first time that cyclo-(l-isoleucyl-l-proline) exhibit nematicidal activity in P. putida. Interestingly, our study had not detected common nematicidal factors such as 2,4-diacetylphloroglucinol (2,4-DAPG) and pyrrolnitrin in P. putida 1A00316. The results of the present study reveal the diversity and specificity of nematicidal factors in P. putida strain 1A00316.

Comparison of Sample Preparation Methods Used for the Next-Generation Sequencing of Mycobacterium tuberculosis

The advent and widespread application of next-generation sequencing (NGS) technologies to the study of microbial genomes has led to a substantial increase in the number of studies in which whole genome sequencing (WGS) is applied to the analysis of microbial genomic epidemiology. However, microorganisms such as Mycobacterium tuberculosis (MTB) present unique problems for sequencing and downstream analysis based on their unique physiology and the composition of their genomes. In this study, we compare the quality of sequence data generated using the Nextera and TruSeq isolate preparation kits for library construction prior to Illumina sequencing-by-synthesis. Our results confirm that MTB NGS data quality is highly dependent on the purity of the DNA sample submitted for sequencing and its guanine-cytosine content (or GC-content). Our data additionally demonstrate that the choice of library preparation method plays an important role in mitigating downstream sequencing quality issues. Importantly for MTB, the Illumina TruSeq library preparation kit produces more uniform data quality than the Nextera XT method, regardless of the quality of the input DNA. Furthermore, specific genomic sequence motifs are commonly missed by the Nextera XT method, as are regions of especially high GC-content relative to the rest of the MTB genome. As coverage bias is highly undesirable, this study illustrates the importance of appropriate protocol selection when performing NGS studies in order to ensure that sound inferences can be made regarding mycobacterial genomes.

Evaluation and Validation of Assembling Corrected PacBio Long Reads for Microbial Genome Completion via Hybrid Approaches

Despite the ever-increasing output of next-generation sequencing data along with developing assemblers, dozens to hundreds of gaps still exist in de novo microbial assemblies due to uneven coverage and large genomic repeats. Third-generation single-molecule, real-time (SMRT) sequencing technology avoids amplification artifacts and generates kilobase-long reads with the potential to complete microbial genome assembly. However, due to the low accuracy (~85%) of third-generation sequences, a considerable amount of long reads (>50X) are required for self-correction and for subsequent de novo assembly. Recently-developed hybrid approaches, using next-generation sequencing data and as few as 5X long reads, have been proposed to improve the completeness of microbial assembly. In this study we have evaluated the contemporary hybrid approaches and demonstrated that assembling corrected long reads (by runCA) produced the best assembly compared to long-read scaffolding (e.g., AHA, Cerulean and SSPACE-LongRead) and gap-filling (SPAdes). For generating corrected long reads, we further examined long-read correction tools, such as ECTools, LSC, LoRDEC, PBcR pipeline and proovread. We have demonstrated that three microbial genomes including Escherichia coli K12 MG1655, Meiothermus ruber DSM1279 and Pdeobacter heparinus DSM2366 were successfully hybrid assembled by runCA into near-perfect assemblies using ECTools-corrected long reads. In addition, we developed a tool, Patch, which implements corrected long reads and pre-assembled contigs as inputs, to enhance microbial genome assemblies. With the additional 20X long reads, short reads of S. cerevisiae W303 were hybrid assembled into 115 contigs using the verified strategy, ECTools + runCA. Patch was subsequently applied to upgrade the assembly to a 35-contig draft genome. Our evaluation of the hybrid approaches shows that assembling the ECTools-corrected long reads via runCA generates near complete microbial genomes, suggesting that genome assembly could benefit from re-analyzing the available hybrid datasets that were not assembled in an optimal fashion.

hybridSPAdes: an algorithm for hybrid assembly of short and long reads

MOTIVATION

Recent advances in single molecule real-time (SMRT) and nanopore sequencing technologies have enabled high-quality assemblies from long and inaccurate reads. However, these approaches require high coverage by long reads and remain expensive. On the other hand, the inexpensive short reads technologies produce accurate but fragmented assemblies. Thus, a hybrid approach that assembles long reads (with low coverage) and short reads has a potential to generate high-quality assemblies at reduced cost.

RESULTS

We describe hybridSPAdes algorithm for assembling short and long reads and benchmark it on a variety of bacterial assembly projects. Our results demonstrate that hybridSPAdes generates accurate assemblies (even in projects with relatively low coverage by long reads) thus reducing the overall cost of genome sequencing. We further present the first complete assembly of a genome from single cells using SMRT reads.

AVAILABILITY AND IMPLEMENTATION

hybridSPAdes is implemented in C++ as a part of SPAdes genome assembler and is publicly available at http://bioinf.spbau.ru/en/spades

CONTACT

d.antipov@spbu.ru

SUPPLEMENTARY INFORMATION

supplementary data are available at Bioinformatics online.

PacBio Sequencing and Its Applications

Single-molecule, real-time sequencing developed by Pacific BioSciences offers longer read lengths than the second-generation sequencing (SGS) technologies, making it well-suited for unsolved problems in genome, transcriptome, and epigenetics research. The highly-contiguous de novo assemblies using PacBio sequencing can close gaps in current reference assemblies and characterize structural variation (SV) in personal genomes. With longer reads, we can sequence through extended repetitive regions and detect mutations, many of which are associated with diseases. Moreover, PacBio transcriptome sequencing is advantageous for the identification of gene isoforms and facilitates reliable discoveries of novel genes and novel isoforms of annotated genes, due to its ability to sequence full-length transcripts or fragments with significant lengths. Additionally, PacBio’s sequencing technique provides information that is useful for the direct detection of base modifications, such as methylation. In addition to using PacBio sequencing alone, many hybrid sequencing strategies have been developed to make use of more accurate short reads in conjunction with PacBio long reads. In general, hybrid sequencing strategies are more affordable and scalable especially for small-size laboratories than using PacBio Sequencing alone. The advent of PacBio sequencing has made available much information that could not be obtained via SGS alone.

Mind the gap; seven reasons to close fragmented genome assemblies

Like other domains of life, research into the biology of filamentous microbes has greatly benefited from the advent of whole-genome sequencing. Next-generation sequencing (NGS) technologies have revolutionized sequencing, making genomic sciences accessible to many academic laboratories including those that study non-model organisms. Thus, hundreds of fungal genomes have been sequenced and are publically available today, although these initiatives have typically yielded considerably fragmented genome assemblies that often lack large contiguous genomic regions. Many important genomic features are contained in intergenic DNA that is often missing in current genome assemblies, and recent studies underscore the significance of non-coding regions and repetitive elements for the life style, adaptability and evolution of many organisms. The study of particular types of genetic elements, such as telomeres, centromeres, repetitive elements, effectors, and clusters of co-regulated genes, but also of phenomena such as structural rearrangements, genome compartmentalization and epigenetics, greatly benefits from having a contiguous and high-quality, preferably even complete and gapless, genome assembly. Here we discuss a number of important reasons to produce gapless, finished, genome assemblies to help answer important biological questions.

Bordetella pertussis evolution in the (functional) genomics era

The incidence of whooping cough caused by Bordetella pertussis in many developed countries has risen dramatically in recent years. This has been linked to the use of an acellular pertussis vaccine. In addition, it is thought that B. pertussis is adapting under acellular vaccine mediated immune selection pressure, towards vaccine escape. Genomics-based approaches have revolutionized the ability to resolve the fine structure of the global B. pertussis population and its evolution during the era of vaccination. Here, we discuss the current picture of B. pertussis evolution and diversity in the light of the current resurgence, highlight import questions raised by recent studies in this area and discuss the role that functional genomics can play in addressing current knowledge gaps.

Single-Molecule Real-Time Sequencing Combined with Optical Mapping Yields Completely Finished Fungal Genome

Next-generation sequencing (NGS) technologies have increased the scalability, speed, and resolution of genomic sequencing and, thus, have revolutionized genomic studies. However, eukaryotic genome sequencing initiatives typically yield considerably fragmented genome assemblies. Here, we assessed various state-of-the-art sequencing and assembly strategies in order to produce a contiguous and complete eukaryotic genome assembly, focusing on the filamentous fungus Verticillium dahliae. Compared with Illumina-based assemblies of the V. dahliae genome, hybrid assemblies that also include PacBio-generated long reads establish superior contiguity. Intriguingly, provided that sufficient sequence depth is reached, assemblies solely based on PacBio reads outperform hybrid assemblies and even result in fully assembled chromosomes. Furthermore, the addition of optical map data allowed us to produce a gapless and complete V. dahliae genome assembly of the expected eight chromosomes from telomere to telomere. Consequently, we can now study genomic regions that were previously not assembled or poorly assembled, including regions that are populated by repetitive sequences, such as transposons, allowing us to fully appreciate an organism’s biological complexity. Our data show that a combination of PacBio-generated long reads and optical mapping can be used to generate complete and gapless assemblies of fungal genomes.

Studying whole-genome sequences has become an important aspect of biological research. The advent of next-generation sequencing (NGS) technologies has nowadays brought genomic science within reach of most research laboratories, including those that study nonmodel organisms. However, most genome sequencing initiatives typically yield (highly) fragmented genome assemblies. Nevertheless, considerable relevant information related to genome structure and evolution is likely hidden in those nonassembled regions. Here, we investigated a diverse set of strategies to obtain gapless genome assemblies, using the genome of a typical ascomycete fungus as the template. Eventually, we were able to show that a combination of PacBio-generated long reads and optical mapping yields a gapless telomere-to-telomere genome assembly, allowing in-depth genome analyses to facilitate functional studies into an organism’s biology.