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

Genomic resources for the Yellowfin tuna Thunnus albacares

BACKGROUND

The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares.

METHODS AND RESULTS

Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology.

CONCLUSIONS

The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid 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.

Phenotypic and genomic characterization of Pseudomonas aeruginosa isolates recovered from catheter-associated urinary tract infections in an Egyptian hospital

Catheter-associated urinary tract infections (CAUTIs) represent one of the major healthcare-associated infections, and Pseudomonas aeruginosa is a common Gram-negative bacterium associated with catheter infections in Egyptian clinical settings. The present study describes the phenotypic and genotypic characteristics of 31 P. aeruginosa isolates recovered from CAUTIs in an Egyptian hospital over a 3 month period. Genomes of isolates were of good quality and were confirmed to be P. aeruginosa by comparison to the type strain (average nucleotide identity, phylogenetic analysis). Clonal diversity among the isolates was determined; eight different sequence types were found (STs 244, 357, 381, 621, 773, 1430, 1667 and 3765), of which ST357 and ST773 are considered to be high-risk clones. Antimicrobial resistance (AMR) testing according to European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines showed that the isolates were highly resistant to quinolones [ciprofloxacin (12/31, 38.7 %) and levofloxacin (9/31, 29 %) followed by tobramycin (10/31, 32.5 %)] and cephalosporins (7/31, 22.5 %). Genotypic analysis of resistance determinants predicted all isolates to encode a range of AMR genes, including those conferring resistance to aminoglycosides, β-lactamases, fluoroquinolones, fosfomycin, sulfonamides, tetracyclines and chloramphenicol. One isolate was found to carry a 422 938 bp pBT2436-like megaplasmid encoding OXA-520, the first report from Egypt of this emerging family of clinically important mobile genetic elements. All isolates were able to form biofilms and were predicted to encode virulence genes associated with adherence, antimicrobial activity, anti-phagocytosis, phospholipase enzymes, iron uptake, proteases, secretion systems and toxins. The present study shows how phenotypic analysis alongside genomic analysis may help us understand the AMR and virulence profiles of P. aeruginosa contributing to CAUTIs in Egypt.

De novo genome assembly resolving repetitive structures enables genomic analysis of 35 European Mycoplasmopsis bovis strains

Mycoplasmopsis (M.) bovis, the agent of mastitis, pneumonia, and arthritis in cattle, harbors a small genome of approximately 1 Mbp. Combining data from Illumina and Nanopore technologies, we sequenced and assembled the genomes of 35 European strains and isolate DL422_88 from Cuba. While the high proportion of repetitive structures in M. bovis genomes represent a particular challenge, implementation of our own pipeline Mycovista (available on GitHub www.github.com/sandraTriebel/mycovista ) in a hybrid approach enabled contiguous assembly of the genomes and, consequently, improved annotation rates considerably. To put our European strain panel in a global context, we analyzed the new genome sequences together with 175 genome assemblies from public databases. Construction of a phylogenetic tree based on core genes of these 219 strains revealed a clustering pattern according to geographical origin, with European isolates positioned on clades 4 and 5. Genomic data allowing assignment of strains to tissue specificity or certain disease manifestations could not be identified. Seven strains isolated from cattle with systemic circular condition (SCC), still a largely unknown manifestation of M. bovis disease, were located on both clades 4 and 5. Pairwise association analysis revealed 108 genomic elements associated with a particular clade of the phylogenetic tree. Further analyzing these hits, 25 genes are functionally annotated and could be linked to a M. bovis protein, e.g. various proteases and nucleases, as well as ten variable surface lipoproteins (Vsps) and other surface proteins. These clade-specific genes could serve as useful markers in epidemiological and clinical surveys.

A Novel Robust Screening Assay Identifies Pseudomonas Strains as Reliable Antagonists of the Root-Knot Nematode Meloidogyne incognita

Forty-four bacterial strains isolated from greenhouse soil and beetroots were tested for their antagonistic activity against the plant-parasitic root-knot nematode (RKN) Meloidogyne incognita, which causes significant yield losses in a number of important crops worldwide. Through a novel combination of in vitro and on planta screening assays, Pseudomonas spp. 105 and 108 were identified as the most promising bacterial isolates. Both strains were evaluated for their potential to control different RKN population densities and as root protectants against nematode infestation. Regardless of the application method, both strains significantly reduced root galling caused by M. incognita. These two strains were subjected to whole genome sequencing and de novo genome assembly as a basis for phylogenetic and future functional characterization. Phylogenetic analysis revealed that both Pseudomonas strains cluster within the Pseudomonas fluorescens clade among previously characterized RKN antagonists and Pseudomonas-based biocontrol agents of plant diseases.

Elucidation of the Bovine Intramammary Bacteriome and Resistome from healthy cows of Swiss dairy farms in the Canton Tessin

Healthy, untreated cows of nine dairy herds from the Swiss Canton Tessin were analyzed three times within one year to identify the most abundant species of the intramammary bacteriome. Aseptically collected milk samples were cultured and bacteria identified using MALDI-TOF. Of 256 cows analyzed, 96% were bacteriologically positive and 80% of the 1,024 quarters were positive for at least one bacterial species. 84.5% of the quarters were healthy with somatic cell counts (SCC) < 200,000 cells/mL, whereas 15.5% of the quarters showed a subclinical mastitis (SCC ≥ 200,000 cells/mL). We could assign 1,288 isolates to 104 different bacterial species including 23 predominant species. Non-aureus staphylococci and mammaliicocci (NASM) were most prevalent (14 different species; 73.5% quarters). Staphylococcus xylosus and Mammaliicoccus sciuri accounted for 74.7% of all NASM isolates. To describe the intramammary resistome, 350 isolates of the predominant species were selected and subjected to short-read whole genome sequencing (WGS) and phenotypic antibiotic resistance profiling. While complete genomes of eight type strains were available, the remaining 15 were de novo assembled with long reads as a resource for the community. The 23 complete genomes served for reference-based assembly of the Illumina WGS data. Both chromosomes and mobile genetic elements were examined for antibiotic resistance genes (ARGs) using in-house and online software tools. ARGs were then correlated with phenotypic antibiotic resistance data from minimum inhibitory concentration (MIC). Phenotypic and genomic antimicrobial resistance was isolate-specific. Resistance to clindamycin and oxacillin was most frequently observed (65 and 30%) in Staphylococcus xylosus but could not be linked to chromosomal or plasmid-borne ARGs. However, in several cases, the observed antimicrobial resistance could be explained by the presence of mobile genetic elements like tetK carried on small plasmids. This represents a possible mechanism of transfer between non-pathogenic bacteria and pathogens of the mammary gland within and between herds. The-to our knowledge-most extensive bacteriome reported and the first attempt to link it with the resistome promise to profoundly affect veterinary bacteriology in the future and are highly relevant in a One Health context, in particular for mastitis, the treatment of which still heavily relies on antibiotics.

Insights into the Diversity of Transcription Activator-Like Effectors in Indian Pathotype Strains of Xanthomonas oryzae pv. oryzae

Xanthomonas oryzae pv. oryzae (Xoo) is a major rice pathogen, and its genome harbors extensive inter-strain and inter-lineage variations. The emergence of highly virulent pathotypes of Xoo that can overcome major resistance (R) genes deployed in rice breeding programs is a grave threat to rice cultivation. The present study reports on a long-read Oxford nanopore-based complete genomic investigation of Xoo isolates from 11 pathotypes that are reported based on their reaction toward 10 R genes. The investigation revealed remarkable variation in the genome structure in the strains belonging to different pathotypes. Furthermore, transcription activator-like effector (TALE) proteins secreted by the type III secretion system display marked variation in content, genomic location, classes, and DNA-binding domain. We also found the association of tal genes in the vicinity of regions with genome structural variations. Furthermore, in silico analysis of the genome-wide rice targets of TALEs allowed us to understand the emergence of pathotypes compatible with major R genes. Long-read, cost-effective sequencing technologies such as nanopore can be a game changer in the surveillance of major and emerging pathotypes. The resource and findings will be invaluable in the management of Xoo and in appropriate deployment of R genes in rice breeding programs.

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.

MEGAnnotator2: a pipeline for the assembly and annotation of microbial genomes

The reconstruction of microbial genome sequences by bioinformatic pipelines and the consequent functional annotation of their genes' repertoire are fundamental activities aiming at unveiling their biological mechanisms, such as metabolism, virulence factors, and antimicrobial resistances. Here, we describe the development of the MEGAnnotator2 pipeline able to manage all next-generation sequencing methodologies producing short- and long-read DNA sequences. Starting from raw sequencing data, the updated pipeline can manage multiple analyses leading to the assembly of high-quality genome sequences and the functional classification of their genetic repertoire, providing the user with a useful report constituting features and statistics related to the microbial genome. The updated pipeline is fully automated from the installation to the delivery of the output, thus requiring minimal bioinformatics knowledge to be executed.

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.

Whole-Genome Sequence of a Clinical Enterococcus gallinarum Strain in a Positive Hemoculture from Siriraj Hospital, Thailand

Enterococcus gallinarum is a rare causative agent of hospital-acquired bacteremia. Here, we reported the complete genome of E. gallinarum strain WKB01, which was directly isolated from a positive hemoculture in Siriraj Hospital, Thailand. Comprehensive analysis demonstrated the clinically relevant tetM gene presenting in its genome.

Hybrid assembly using long reads resolves repeats and completes the genome sequence of a laboratory strain of Staphylococcus aureus subsp. Aureus RN4220

Staphylococcus aureus RN4220 has been extensively used by staphylococcal researchers as an intermediate strain for genetic manipulation due to its ability to accept foreign DNA. Despite its wide use in laboratories, its complete genome is not available. In this study, we used a hybrid genome assembly approach using minION long reads and Illumina short reads to sequence the complete genome of S. aureus RN4220. The comparative analysis of the annotated complete genome showed the presence of 39 genes fragmented in the previous assembly, many of which were located near the repeat regions. Using RNA-Seq reads, we showed that a higher number of reads could be mapped to the complete genome than the draft genome and the gene expression profile obtained using the complete genome also differs from that obtained from the draft genome. Furthermore, by comparative transcriptomic analysis, we showed the correlation between expression levels of staphyloxanthin biosynthetic genes and the production of yellow pigment. This study highlighted the importance of long reads in completing microbial genomes, especially those possessing repetitive elements.

•

S. aureus RN4220 is used as an intermediate strain for genetic manipulation.

•

We completed its genome and found 39 fragmented genes in previous genome assembly.

•

RNA-Seq analysis improved mapping of the reads with the use of complete genome.

•

Expression of staphyloxanthin biosynthetic genes was correlated with its production.

Oral microbiome diversity: The curious case of Corynebacterium sp. isolation

Oral microbiome sequencing efforts revealed the presence of hundreds of different microbes. Interindividual differences at strain and species resolution suggest that microbiome diversity could lead to mechanistically distinct gene regulation as well as species-related differences in phenotypes. Commonly, gene regulation and related phenotypes are studied in a few selected strains of a particular species with conclusions that are mostly generalized. The aim of this study was to isolate several species of Corynebacterium using an established protocol that led to the previous isolation of C. durum. Characterization of C. durum interspecies interactions revealed a specific mechanism for chain elongation in Streptococcus sanguinis that was the result of corynebacterial fatty acid production and secretion. While the protocol was successfully applied to isolate what we presumed to be additional Corynebacterium based on several phenotypic traits that seem to be identical to C. durum, genome sequencing of the newly isolated strains placed them closer to Actinomyces. Both Corynebacterium and Actinomyces are suborders of the Actinobacteridae and related species. Our study suggests to take several comprehensive strategies into consideration when taxonomically identifying closely related microorganisms. Furthermore, it seems to be important to test common core phenotypes in bacterial ecology to understand the behavior of specific groups of microbes, rather than simply relying upon genome sequence homology to establish relationships in the microbiome.

Molecular epidemiology and population genomics of tet(X4), blaNDM or mcr-1 positive Escherichia coli from migratory birds in southeast coast of China

The emergence of multidrug-resistant (MDR) bacteria harboring tet(X4), bla_NDM or mcr-1 posed a serious threat to public health. Wild birds, especially migratory birds, were considered as one of important transmission vectors for antibiotic resistance genes (ARGs) globally, however, few studies were performed on the genomic epidemiology of critical resistance genes among them. Isolates harboring tet(X4), mcr-1 or bla_NDM from migratory birds were identified and characterized by PCR, antimicrobial susceptibility testing, conjugation assays, whole genome sequencing and bioinformatics analysis. A total of 14 tet(X4)-bearing E. coli, 4 bla_NDM-bearing E. coli and 23 mcr-1-bearing E. coli isolates were recovered from 1060 fecal samples of migratory birds. All isolates were MDR bacteria and most plasmids carrying tet(X4), bla_NDM or mcr-1 were conjugative. We first identified an E. coli of migratory bird origin carrying bla_NDM-4, which was located on a conjugative IncHI2 plasmid and embedded on a novel MDR region flanked by IS26 that could generate the circular intermediate. The emergency of E. coli isolates co-harboring mcr-1 and bla_NDM-5 in migratory birds indicated the coexistence of ARGs in migratory birds was a novel threat. This study revealed the prevalence and molecular characteristics of three important ARGs in migratory birds, provided evidence that migratory birds were potential vectors of novel resistance genes and highlighted the monitoring of ARGs in migratory birds should be strengthened to prevent the spread of ARGs in a One Health strategy.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Benchmarking the empirical accuracy of short-read sequencing across the M. tuberculosis genome

MOTIVATION

Short-read whole-genome sequencing (WGS) is a vital tool for clinical applications and basic research. Genetic divergence from the reference genome, repetitive sequences and sequencing bias reduces the performance of variant calling using short-read alignment, but the loss in recall and specificity has not been adequately characterized. To benchmark short-read variant calling, we used 36 diverse clinical Mycobacterium tuberculosis (Mtb) isolates dually sequenced with Illumina short-reads and PacBio long-reads. We systematically studied the short-read variant calling accuracy and the influence of sequence uniqueness, reference bias and GC content.

RESULTS

Reference-based Illumina variant calling demonstrated a maximum recall of 89.0% and minimum precision of 98.5% across parameters evaluated. The approach that maximized variant recall while still maintaining high precision (<99%) was tuning the mapping quality filtering threshold, i.e. confidence of the read mapping (recall = 85.8%, precision = 99.1%, MQ ≥ 40). Additional masking of repetitive sequence content is an alternative conservative approach to variant calling that increases precision at cost to recall (recall = 70.2%, precision = 99.6%, MQ ≥ 40). Of the genomic positions typically excluded for Mtb, 68% are accurately called using Illumina WGS including 52/168 PE/PPE genes (34.5%). From these results, we present a refined list of low confidence regions across the Mtb genome, which we found to frequently overlap with regions with structural variation, low sequence uniqueness and low sequencing coverage. Our benchmarking results have broad implications for the use of WGS in the study of Mtb biology, inference of transmission in public health surveillance systems and more generally for WGS applications in other organisms.

AVAILABILITY AND IMPLEMENTATION

All relevant code is available at https://github.com/farhat-lab/mtb-illumina-wgs-evaluation.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Disentangle genus microdiversity within a complex microbial community by using a multi-distance long-read binning method: example of Candidatus Accumulibacter

Complete genomes can be recovered from metagenomes by assembling and binning DNA sequences into metagenome assembled genomes (MAGs). Yet, the presence of microdiversity can hamper the assembly and binning processes, possibly yielding chimeric, highly fragmented and incomplete genomes. Here, the metagenomes of four samples of aerobic granular sludge bioreactors containing Candidatus (Ca.) Accumulibacter, a phosphate-accumulating organism of interest for wastewater treatment, were sequenced with both PacBio and Illumina. Different strategies of genome assembly and binning were investigated, including published protocols and a binning procedure adapted to the binning of long contigs (MuLoBiSC). Multiple criteria were considered to select the best strategy for Ca. Accumulibacter, whose multiple strains in every sample represent a challenging microdiversity. In this case, the best strategy relies on long-read only assembly and a custom binning procedure including MuLoBiSC in metaWRAP. Several high-quality Ca. Accumulibacter MAGs, including a novel species, were obtained independently from different samples. Comparative genomic analysis showed that MAGs retrieved in different samples harbour genomic rearrangements in addition to accumulation of point mutations. The microdiversity of Ca. Accumulibacter, likely driven by mobile genetic elements, causes major difficulties in recovering MAGs, but it is also a hallmark of the panmictic lifestyle of these bacteria.

A Comparative Genomic and Safety Assessment of Six Lactiplantibacillus plantarum subsp. argentoratensis Strains Isolated from Spontaneously Fermented Greek Wheat Sourdoughs for Potential Biotechnological Application

The comparative genome analysis of six Lactiplantibacillus plantarum subsp. argentoratensis strains previously isolated from spontaneously fermented Greek wheat sourdoughs is presented. Genomic attributes related to food safety have been studied according to the European Food Safety Authority (EFSA) suggestions for the use of lactic acid bacteria (LAB) in the production of foods. Bioinformatic analysis revealed a complete set of genes for maltose, sucrose, glucose, and fructose fermentation; conversion of fructose to mannitol; folate and riboflavin biosynthesis; acetoin production; conversion of citrate to oxaloacetate; and the ability to produce antimicrobial compounds (plantaricins). Pathogenic factors were absent but some antibiotic resistance genes were detected. CRISPR and cas genes were present as well as various mobile genetic elements (MGEs) such as plasmids, prophages, and insertion sequences. The production of biogenic amines by these strains was not possible due to the absence of key genes in their genome except lysine decarboxylase associated with cadaverine; however, potential degradation of these substances was identified due to the presence of a blue copper oxidase precursor and a multicopper oxidase protein family. Finally, comparative genomics and pan-genome analysis showed genetic differences between the strains (e.g., variable pln locus), and it facilitated the identification of various phenotypic and probiotic-related properties.

Application and Challenge of 3rd Generation Sequencing for Clinical Bacterial Studies

Over the past 25 years, the powerful combination of genome sequencing and bioinformatics analysis has played a crucial role in interpreting information encoded in bacterial genomes. High-throughput sequencing technologies have paved the way towards understanding an increasingly wide range of biological questions. This revolution has enabled advances in areas ranging from genome composition to how proteins interact with nucleic acids. This has created unprecedented opportunities through the integration of genomic data into clinics for the diagnosis of genetic traits associated with disease. Since then, these technologies have continued to evolve, and recently, long-read sequencing has overcome previous limitations in terms of accuracy, thus expanding its applications in genomics, transcriptomics and metagenomics. In this review, we describe a brief history of the bacterial genome sequencing revolution and its application in public health and molecular epidemiology. We present a chronology that encompasses the various technological developments: whole-genome shotgun sequencing, high-throughput sequencing, long-read sequencing. We mainly discuss the application of next-generation sequencing to decipher bacterial genomes. Secondly, we highlight how long-read sequencing technologies go beyond the limitations of traditional short-read sequencing. We intend to provide a description of the guiding principles of the 3rd generation sequencing applications and ongoing improvements in the field of microbial medical research.

Genome sequence data of the antagonistic soil-borne yeast Cyberlindnera sargentensis (SHA 17.2)

Cyberlindnera sargentensis strain SHA 17.2, isolated from a Swiss soil sample, exhibited strong antagonistic activity against several plant pathogenic fungi in vitro and was highly competitive against other yeasts in soil. As a basis for identifying the mechanisms underlying its strong antagonistic activity, we have sequenced the genome of C. sargentensis (SHA 17.2) by long- and short read sequencing, de novo assembled them into seven contigs/chromosomes and a mitogenome (total genome size 11.4 Mbp), and annotated 5455 genes. This high-quality genome is the reference for transcriptome and proteome analyses aiming at elucidating the mode of action of C. sargentensis against fungal plant pathogens. It will thus serve as a resource for identifying potential biocontrol genes and performing comparative genomics analyses of yeast genomes.

Emergence of the cfr Gene in Vibrio diabolicus of Seafood Origin

The emergence and transmission of multidrug resistance (MDR) gene cfr have incurred great public health concerns worldwide. Recently, Gram-negative pathogens were found to carry cfr by various mobile elements. Here, we investigated a cfr-positive Vibrio diabolicus isolate by phenotyping and genomic analysis and found cfr in a translocatable structure (IS26-hp-cfr-IS26) among the MDR region in pNV27-cfr-208K, an emerging MDR plasmid in Vibrio species. This study highlights the necessity of surveillance of cfr in bacteria of diverse origins.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Identification of bacterial antibiotic resistance genes in next-generation sequencing data (review of literature)

The spread of antibiotic-resistant human bacterial pathogens is a serious threat to modern medicine. Antibiotic susceptibility testing is essential for treatment regimens optimization and preventing dissemination of antibiotic resistance. Therefore, development of antibiotic susceptibility testing methods is a priority challenge of laboratory medicine. The aim of this review is to analyze the capabilities of the bioinformatics tools for bacterial whole genome sequence data processing. The PubMed database, Russian scientific electronic library eLIBRARY, information networks of World health organization and European Society of Clinical Microbiology and Infectious Diseases (ESCMID) were used during the analysis. In this review, the platforms for whole genome sequencing, which are suitable for detection of bacterial genetic resistance determinants, are described. The classic step of genetic resistance determinants searching is an alignment between the query nucleotide/protein sequence and the subject (database) nucleotide/protein sequence, which is performed using the nucleotide and protein sequence databases. The most commonly used databases are Resfinder, CARD, Bacterial Antimicrobial Resistance Reference Gene Database. The results of the resistance determinants searching in genome assemblies is more correct in comparison to results of the searching in contigs. The new resistance genes searching bioinformatics tools, such as neural networks and machine learning, are discussed in the review. After critical appraisal of the current antibiotic resistance databases we designed a protocol for predicting antibiotic resistance using whole genome sequence data. The designed protocol can be used as a basis of the algorithm for qualitative and quantitative antimicrobial susceptibility testing based on whole genome sequence data.

Evaluation of WGS performance for bacterial pathogen characterization with the Illumina technology optimized for time-critical situations

Whole genome sequencing (WGS) has become the reference standard for bacterial outbreak investigation and pathogen typing, providing a resolution unattainable with conventional molecular methods. Data generated with Illumina sequencers can however only be analysed after the sequencing run has finished, thereby losing valuable time during emergency situations. We evaluated both the effect of decreasing overall run time, and also a protocol to transfer and convert intermediary files generated by Illumina sequencers enabling real-time data analysis for multiple samples part of the same ongoing sequencing run, as soon as the forward reads have been sequenced. To facilitate implementation for laboratories operating under strict quality systems, extensive validation of several bioinformatics assays (16S rRNA species confirmation, gene detection against virulence factor and antimicrobial resistance databases, SNP-based antimicrobial resistance detection, serotype determination, and core genome multilocus sequence typing) for three bacterial pathogens (Mycobacterium tuberculosis , Neisseria meningitidis , and Shiga-toxin producing Escherichia coli ) was performed by evaluating performance in function of the two most critical sequencing parameters, i.e. read length and coverage. For the majority of evaluated bioinformatics assays, actionable results could be obtained between 14 and 22 h of sequencing, decreasing the overall sequencing-to-results time by more than half. This study aids in reducing the turn-around time of WGS analysis by facilitating a faster response in time-critical scenarios and provides recommendations for time-optimized WGS with respect to required read length and coverage to achieve a minimum level of performance for the considered bioinformatics assay(s), which can also be used to maximize the cost-effectiveness of routine surveillance sequencing when response time is not essential.

Towards population-scale long-read sequencing

Long-read sequencing technologies have now reached a level of accuracy and yield that allows their application to variant detection at a scale of tens to thousands of samples. Concomitant with the development of new computational tools, the first population-scale studies involving long-read sequencing have emerged over the past 2 years and, given the continuous advancement of the field, many more are likely to follow. In this Review, we survey recent developments in population-scale long-read sequencing, highlight potential challenges of a scaled-up approach and provide guidance regarding experimental design. We provide an overview of current long-read sequencing platforms, variant calling methodologies and approaches for de novo assemblies and reference-based mapping approaches. Furthermore, we summarize strategies for variant validation, genotyping and predicting functional impact and emphasize challenges remaining in achieving long-read sequencing at a population scale.

Identification of a Novel Plasmid-Mediated Carbapenemase-Encoding Gene, blaVMB-2, in Vibrio diabolicus

Here, we characterized a carbapenem-resistant Vibrio diabolicus strain of shrimp origin with various experiments and bioinformatics analysis. A novel metallo-β-lactamase (MBL) gene, blaVMB-2, that confers resistance to β-lactams, including meropenem and cephalosporins, was identified on a plasmid-borne composite transposon, ISShfr9-ISCR1-blaVMB-2-blaCARB-12-aadA1-ISShfr9, which is capable of generating a blaVMB-2-bearing circular intermediate. ISShfr9 was found to be disseminated in multidrug-resistant (MDR) pathogens, arousing the concern of further transmission of a blaVMB-2-bearing circular intermediate to clinical Enterobacterales strains via such insertion sequences, which warrants further investigations.

Plasmids of novel incompatibility group IncpRBL16 from Pseudomonas species

OBJECTIVES

To dissect genomic features of IncpRBL16 plasmids from Pseudomonas.

METHODS

An extensive genomic comparison was applied to all 17 available sequenced IncpRBL16 plasmids, including 8 sequenced in this study and another 2 sequenced in two of our previous studies.

RESULTS

Conserved IncpRBL16 backbone markers repAIncpRBL16 together with its iterons, parB2-parA, che, pil and ter were present in all 17 plasmids. At least 18 regions or sites across IncpRBL16 genomes exhibited major modular differences, including insertion of accessory modules, deletion of backbone regions surrounding insertion sites and substitution of multiple-gene backbone regions. Ten plasmids carried a sole IncpRBL16 replicon, while exogenous acquisition of an auxiliary replicon (located in an accessory module) besides the primary IncpRBL16 replicon was observed in each of the remaining seven plasmids. The 17 IncpRBL16 plasmids carried at least 71 different accessory modules, notably including Tn1403-related regions, Tn7-family transposons, Tn6571-family transposons, integrative and conjugative elements, and integrative and mobilizable elements. There were a total of 40 known resistance genes, which were involved in resistance to 15 categories of antibiotics and heavy metals, notably including blaIMP-9, blaIMP-45, blaVIM-2, blaDIM-2, blaOXA-246, blaPER-1, aphA and armA.

CONCLUSIONS

Different IncpRBL16 plasmids contain different profiles of accessory modules and thus diverse collections of resistance genes. To the best of our knowledge, this is the first report of fully sequenced blaOXA-246-carrying (p12939-PER) and blaPER-1-carrying (p12939-PER and pA681-IMP) IncpRBL16 plasmids and also that of 14 novel (first identified in this study) and additionally 31 newly named (first designated in this study, but with previously determined sequences) mobile elements.

Biosurfactants Produced by Phyllosphere-Colonizing Pseudomonads Impact Diesel Degradation but Not Colonization of Leaves of Gnotobiotic Arabidopsis thaliana

Biosurfactant production is a common trait in leaf surface-colonizing bacteria that has been associated with increased survival and movement on leaves. At the same time, the ability to degrade aliphatics is common in biosurfactant-producing leaf colonizers. Pseudomonads are common leaf colonizers and have been recognized for their ability to produce biosurfactants and degrade aliphatic compounds. In this study, we investigated the role of biosurfactants in four non-plant-pathogenic Pseudomonas strains by performing a series of experiments to characterize their surfactant properties and their role during leaf colonization and diesel degradation. The biosurfactants produced were identified using mass spectrometry. Two strains produced viscosin-like biosurfactants, and the other two produced massetolide A-like biosurfactants, which aligned with the phylogenetic relatedness between the strains. To further investigate the role of surfactant production, random Tn5 transposon mutagenesis was performed to generate knockout mutants. The knockout mutants were compared to their respective wild types with regard to their ability to colonize gnotobiotic Arabidopsis thaliana and to degrade diesel or dodecane. It was not possible to detect negative effects during plant colonization in direct competition or individual colonization experiments. When grown on diesel, knockout mutants grew significantly slower than their respective wild types. When grown on dodecane, knockout mutants were less impacted than during growth on diesel. By adding isolated wild-type biosurfactants, it was possible to complement the growth of the knockout mutants.IMPORTANCE Many leaf-colonizing bacteria produce surfactants and are able to degrade aliphatic compounds; however, whether surfactant production provides a competitive advantage during leaf colonization is unclear. Furthermore, it is unclear if leaf colonizers take advantage of the aliphatic compounds that constitute the leaf cuticle and cuticular waxes. Here, we tested the effect of surfactant production on leaf colonization, and we demonstrate that the lack of surfactant production decreases the ability to degrade aliphatic compounds. This indicates that leaf surface-dwelling, surfactant-producing bacteria contribute to degradation of environmental hydrocarbons and may be able to utilize leaf surface waxes. This has implications for plant-microbe interactions and future studies.

Epidemic Territorial Spread of IncP-2-Type VIM-2 Carbapenemase-Encoding Megaplasmids in Nosocomial Pseudomonas aeruginosa Populations

In 2003 to 2004, the first five VIM-2 metallo-β-lactamase (MBL)-producing Pseudomonas aeruginosa (MPPA) isolates with an In4-like integron, In461 (aadB-blaVIM-2-aadA6), on conjugative plasmids were identified in three hospitals in Poland. In 2005 to 2015, MPPA expanded much in the country, and as many as 80 isolates in a collection of 454 MPPA (∼18%) had In461, one of the two most common MBL-encoding integrons. The organisms occurred in 49 hospitals in 33 cities of 11/16 main administrative regions. Pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST) classified them into 55 pulsotypes and 35 sequence types (STs), respectively, revealing their remarkable genetic diversity overall, with only a few small clonal clusters. S1 nuclease/hybridization assays and mating of 63 representative isolates showed that ∼85% of these had large In461-carrying plasmids, ∼350 to 550 kb, usually self-transmitting with high efficiency (∼10-1 to 10-2 per donor cell). The plasmids from 19 isolates were sequenced and subjected to structural and single-nucleotide-polymorphism (SNP)-based phylogenetic analysis. These formed a subgroup within a family of IncP-2-type megaplasmids, observed worldwide in pseudomonads from various environments and conferring resistance/tolerance to multiple stress factors, including antibiotics. Their microdiversity in Poland arose mainly from acquisition of different accessory fragments, as well as new resistance genes and multiplication of these. Short-read sequence and/or PCR mapping confirmed the In461-carrying plasmids in the remaining isolates to be the IncP-2 types. The study demonstrated a large-scale epidemic spread of multidrug resistance plasmids in P. aeruginosa populations, creating an epidemiological threat. It contributes to the knowledge on IncP-2 types, which are interesting research objects in resistance epidemiology, environmental microbiology, and biotechnology.

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.

A comprehensive review of scaffolding methods in genome assembly

In the field of genome assembly, scaffolding methods make it possible to obtain a more complete and contiguous reference genome, which is the cornerstone of genomic research. Scaffolding methods typically utilize the alignments between contigs and sequencing data (reads) to determine the orientation and order among contigs and to produce longer scaffolds, which are helpful for genomic downstream analysis. With the rapid development of high-throughput sequencing technologies, diverse types of reads have emerged over the past decade, especially in long-range sequencing, which have greatly enhanced the assembly quality of scaffolding methods. As the number of scaffolding methods increases, biology and bioinformatics researchers need to perform in-depth analyses of state-of-the-art scaffolding methods. In this article, we focus on the difficulties in scaffolding, the differences in characteristics among various kinds of reads, the methods by which current scaffolding methods address these difficulties, and future research opportunities. We hope this work will benefit the design of new scaffolding methods and the selection of appropriate scaffolding methods for specific biological studies.

Antarctic desert soil bacteria exhibit high novel natural product potential, evaluated through long-read genome sequencing and comparative genomics

Actinobacteria and Proteobacteria are important producers of bioactive natural products (NP), and these phyla dominate in the arid soils of Antarctica, where metabolic adaptations influence survival under harsh conditions. Biosynthetic gene clusters (BGCs) which encode NPs, are typically long and repetitious high G + C regions difficult to sequence with short-read technologies. We sequenced 17 Antarctic soil bacteria from multi-genome libraries, employing the long-read PacBio platform, to optimize capture of BGCs and to facilitate a comprehensive analysis of their NP capacity. We report 13 complete bacterial genomes of high quality and contiguity, representing 10 different cold-adapted genera including novel species. Antarctic BGCs exhibited low similarity to known compound BGCs (av. 31%), with an abundance of terpene, non-ribosomal peptide and polyketide-encoding clusters. Comparative genome analysis was used to map BGC variation between closely related strains from geographically distant environments. Results showed the greatest biosynthetic differences to be in a psychrotolerant Streptomyces strain, as well as a rare Actinobacteria genus, Kribbella, while two other Streptomyces spp. were surprisingly similar to known genomes. Streptomyces and Kribbella BGCs were predicted to encode antitumour, antifungal, antibacterial and biosurfactant-like compounds, and the synthesis of NPs with antibacterial, antifungal and surfactant properties was confirmed through bioactivity assays.

An integrated model system to gain mechanistic insights into biofilm-associated antimicrobial resistance in Pseudomonas aeruginosa MPAO1

Pseudomonas aeruginosa MPAO1 is the parental strain of the widely utilized transposon mutant collection for this important clinical pathogen. Here, we validate a model system to identify genes involved in biofilm growth and biofilm-associated antibiotic resistance. Our model employs a genomics-driven workflow to assemble the complete MPAO1 genome, identify unique and conserved genes by comparative genomics with the PAO1 reference strain and genes missed within existing assemblies by proteogenomics. Among over 200 unique MPAO1 genes, we identified six general essential genes that were overlooked when mapping public Tn-seq data sets against PAO1, including an antitoxin. Genomic data were integrated with phenotypic data from an experimental workflow using a user-friendly, soft lithography-based microfluidic flow chamber for biofilm growth and a screen with the Tn-mutant library in microtiter plates. The screen identified hitherto unknown genes involved in biofilm growth and antibiotic resistance. Experiments conducted with the flow chamber across three laboratories delivered reproducible data on P. aeruginosa biofilms and validated the function of both known genes and genes identified in the Tn-mutant screens. Differential protein abundance data from planktonic cells versus biofilm confirmed the upregulation of candidates known to affect biofilm formation, of structural and secreted proteins of type VI secretion systems, and provided proteogenomic evidence for some missed MPAO1 genes. This integrated, broadly applicable model promises to improve the mechanistic understanding of biofilm formation, antimicrobial tolerance, and resistance evolution in biofilms.

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.

NanoGalaxy: Nanopore long-read sequencing data analysis in Galaxy

BACKGROUND

Long-read sequencing can be applied to generate very long contigs and even completely assembled genomes at relatively low cost and with minimal sample preparation. As a result, long-read sequencing platforms are becoming more popular. In this respect, the Oxford Nanopore Technologies-based long-read sequencing "nanopore" platform is becoming a widely used tool with a broad range of applications and end-users. However, the need to explore and manipulate the complex data generated by long-read sequencing platforms necessitates accompanying specialized bioinformatics platforms and tools to process the long-read data correctly. Importantly, such tools should additionally help democratize bioinformatics analysis by enabling easy access and ease-of-use solutions for researchers.

RESULTS

The Galaxy platform provides a user-friendly interface to computational command line-based tools, handles the software dependencies, and provides refined workflows. The users do not have to possess programming experience or extended computer skills. The interface enables researchers to perform powerful bioinformatics analysis, including the assembly and analysis of short- or long-read sequence data. The newly developed "NanoGalaxy" is a Galaxy-based toolkit for analysing long-read sequencing data, which is suitable for diverse applications, including de novo genome assembly from genomic, metagenomic, and plasmid sequence reads.

CONCLUSIONS

A range of best-practice tools and workflows for long-read sequence genome assembly has been integrated into a NanoGalaxy platform to facilitate easy access and use of bioinformatics tools for researchers. NanoGalaxy is freely available at the European Galaxy server https://nanopore.usegalaxy.eu with supporting self-learning training material available at https://training.galaxyproject.org.

Long-read only assembly of Drechmeria coniospora genomes reveals widespread chromosome plasticity and illustrates the limitations of current nanopore methods

BACKGROUND

Long-read sequencing is increasingly being used to determine eukaryotic genomes. We used nanopore technology to generate chromosome-level assemblies for 3 different strains of Drechmeria coniospora, a nematophagous fungus used extensively in the study of innate immunity in Caenorhabditis elegans.

RESULTS

One natural geographical isolate demonstrated high stability over decades, whereas a second isolate not only had a profoundly altered genome structure but exhibited extensive instability. We conducted an in-depth analysis of sequence errors within the 3 genomes and established that even with state-of-the-art tools, nanopore methods alone are insufficient to generate eukaryotic genome sequences of sufficient accuracy to merit inclusion in public databases.

CONCLUSIONS

Although nanopore long-read sequencing is not accurate enough to produce publishable eukaryotic genomes, in our case, it has revealed new information about genome plasticity in D. coniospora and provided a backbone that will permit future detailed study to characterize gene evolution in this important model fungal pathogen.

Harnessing the Microbiomes of Suppressive Composts for Plant Protection: From Metagenomes to Beneficial Microorganisms and Reliable Diagnostics

Soil-borne diseases cause significant yield losses worldwide, are difficult to treat and often only limited options for disease management are available. It has long been known that compost amendments, which are routinely applied in organic and integrated farming as a part of good agricultural practice to close nutrient cycles, can convey a protective effect. Yet, the targeted use of composts against soil-borne diseases is hampered by the unpredictability of the efficacy. Several studies have identified and/or isolated beneficial microorganisms (i.e., bacteria, oomycetes, and fungi) from disease suppressive composts capable of suppressing pathogens (e.g., Pythium and Fusarium) in various crops (e.g., tomato, lettuce, and cucumber), and some of them have been developed into commercial products. Yet, there is growing evidence that synthetic or complex microbial consortia can be more effective in controlling diseases than single strains, but the underlying molecular mechanisms are poorly understood. Currently, a major bottleneck concerns the lack of functional assays to identify the most potent beneficial microorganisms and/or key microbial consortia from complex soil and compost microbiomes, which can harbor tens of thousands of species. This focused review describes microorganisms, which have been isolated from, amended to or found to be abundant in disease-suppressive composts and for which a beneficial effect has been documented. We point out opportunities to increasingly harness compost microbiomes for plant protection through an integrated systems approach that combines the power of functional assays to isolate biocontrol and plant growth promoting strains and further prioritize them, with functional genomics approaches that have been successfully applied in other fields of microbiome research. These include detailed metagenomics studies (i.e., amplicon and shotgun sequencing) to achieve a better understanding of the complex system compost and to identify members of taxa enriched in suppressive composts. Whole-genome sequencing and complete assembly of key isolates and their subsequent functional profiling can elucidate the mechanisms of action of biocontrol strains. Integrating the benefits of these approaches will bring the long-term goals of employing microorganisms for a sustainable control of plant pathogens and developing reliable diagnostic assays to assess the suppressiveness of composts within reach.

Complete genome sequences of Streptococcus pyogenes type strain reveal 100%-match between PacBio-solo and Illumina-Oxford Nanopore hybrid assemblies

We present the first complete, closed genome sequences of Streptococcus pyogenes strains NCTC 8198T and CCUG 4207T, the type strain of the type species of the genus Streptococcus and an important human pathogen that causes a wide range of infectious diseases. S. pyogenes NCTC 8198T and CCUG 4207T are derived from deposit of the same strain at two different culture collections. NCTC 8198T was sequenced, using a PacBio platform; the genome sequence was assembled de novo, using HGAP. CCUG 4207T was sequenced and a de novo hybrid assembly was generated, using SPAdes, combining Illumina and Oxford Nanopore sequence reads. Both strategies yielded closed genome sequences of 1,914,862 bp, identical in length and sequence identity. Combining short-read Illumina and long-read Oxford Nanopore sequence data circumvented the expected error rate of the nanopore sequencing technology, producing a genome sequence indistinguishable to the one determined with PacBio. Sequence analyses revealed five prophage regions, a CRISPR-Cas system, numerous virulence factors and no relevant antibiotic resistance genes. These two complete genome sequences of the type strain of S. pyogenes will effectively serve as valuable taxonomic and genomic references for infectious disease diagnostics, as well as references for future studies and applications within the genus Streptococcus.

Complete genome sequence and epigenetic profile of Bacillus velezensis UCMB5140 used for plant and crop protection in comparison with other plant-associated Bacillus strains

The application of biocontrol biopesticides based on plant growth-promoting rhizobacteria (PGPR), particularly members of the genus Bacillus, is considered a promising perspective to make agricultural practices sustainable and ecologically safe. Recent advances in genome sequencing by third-generation sequencing technologies, e.g., Pacific Biosciences' Single Molecule Real-Time (PacBio SMRT) platform, have allowed researchers to gain deeper insights into the molecular and genetic mechanisms of PGPR activities, and to compare whole genome sequences and global patterns of epigenetic modifications. In the current work, this approach was used to sequence and compare four Bacillus strains that exhibited various PGPR activities including the strain UCMB5140, which is used in the commercial biopesticide Phytosubtil. Whole genome comparison and phylogenomic inference assigned the strain UCMB5140 to the species Bacillus velezensis. Strong biocontrol activities of this strain were confirmed in several bioassays. Several factors that affect the evolution of active PGPR B. velezensis strains were identified: (1) horizontal acquisition of novel non-ribosomal peptide synthetases (NRPS) and adhesion genes; (2) rearrangements of functional modules of NRPS genes leading to strain specific combinations of their encoded products; (3) gain and loss of methyltransferases that can cause global alterations in DNA methylation patterns, which eventually may affect gene expression and regulate transcription. Notably, we identified a horizontally transferred NRPS operon encoding an uncharacterized polypeptide antibiotic in B. velezensis UCMB5140. Other horizontally acquired genes comprised a possible adhesin and a methyltransferase, which may explain the strain-specific methylation pattern of the chromosomal DNA of UCMB5140. KEY POINTS: • Whole genome sequence of the active PGPR Bacillus velezensis UCMB5140. • Identification of genetic determinants responsible for PGPR activities. • Role of methyltransferases and epigenetic mechanisms in evolution of bacteria.

Chromosomal toxin-antitoxin systems in Pseudomonas putida are rather selfish than beneficial

Chromosomal toxin-antitoxin (TA) systems are widespread genetic elements among bacteria, yet, despite extensive studies in the last decade, their biological importance remains ambivalent. The ability of TA-encoded toxins to affect stress tolerance when overexpressed supports the hypothesis of TA systems being associated with stress adaptation. However, the deletion of TA genes has usually no effects on stress tolerance, supporting the selfish elements hypothesis. Here, we aimed to evaluate the cost and benefits of chromosomal TA systems to Pseudomonas putida. We show that multiple TA systems do not confer fitness benefits to this bacterium as deletion of 13 TA loci does not influence stress tolerance, persistence or biofilm formation. Our results instead show that TA loci are costly and decrease the competitive fitness of P. putida. Still, the cost of multiple TA systems is low and detectable in certain conditions only. Construction of antitoxin deletion strains showed that only five TA systems code for toxic proteins, while other TA loci have evolved towards reduced toxicity and encode non-toxic or moderately potent proteins. Analysis of P. putida TA systems’ homologs among fully sequenced Pseudomonads suggests that the TA loci have been subjected to purifying selection and that TA systems spread among bacteria by horizontal gene transfer.

Linking Comparative Genomics of Nine Potato-Associated Pseudomonas Isolates With Their Differing Biocontrol Potential Against Late Blight

For plants, the advantages of associating with beneficial bacteria include plant growth promotion, reduction of abiotic and biotic stresses and enhanced protection against various pests and diseases. Beneficial bacteria rightly equipped for successful plant colonization and showing antagonistic activity toward plant pathogens seem to be actively recruited by plants. To gain more insights into the genetic determinants responsible for plant colonization and antagonistic activities, we first sequenced and de novo assembled the complete genomes of nine Pseudomonas strains that had exhibited varying antagonistic potential against the notorious oomycete Phytophthora infestans, placed them into the phylogenomic context of known Pseudomonas biocontrol strains and carried out a comparative genomic analysis to define core, accessory (i.e., genes found in two or more, but not all strains) and unique genes. Next, we assessed the colonizing abilities of these strains and used bioassays to characterize their inhibitory effects against different stages of P. infestans' lifecycle. The phenotype data were then correlated with genotype information, assessing over three hundred genes encoding known factors for plant colonization and antimicrobial activity as well as secondary metabolite biosynthesis clusters predicted by antiSMASH. All strains harbored genes required for successful plant colonization but also distinct arsenals of antimicrobial compounds. We identified genes coding for phenazine, hydrogen cyanide, 2-hexyl, 5-propyl resorcinol and pyrrolnitrin synthesis, as well as various siderophores, pyocins and type VI secretion systems. Additionally, the comparative genomic analysis revealed about a hundred accessory genes putatively involved in anti-Phytophthora activity, including a type II secretion system (T2SS), several peptidases and a toxin. Transcriptomic studies and mutagenesis are needed to further investigate the putative involvement of the novel candidate genes and to identify the various mechanisms involved in the inhibition of P. infestans by different Pseudomonas strains.

The Hypercomplex Genome of an Insect Reproductive Parasite Highlights the Importance of Lateral Gene Transfer in Symbiont Biology

The biology of many bacteria is critically dependent on genes carried on plasmid and phage mobile elements. These elements shuttle between microbial species, thus providing an important source of biological innovation across taxa. It has recently been recognized that mobile elements are also important in symbiotic bacteria, which form long-lasting interactions with their host. In this study, we report a bacterial symbiont genome that carries a highly complex array of these elements. Arsenophonus nasoniae is the son-killer microbe of the parasitic wasp Nasonia vitripennis and exists with the wasp throughout its life cycle. We completed its genome with the aid of recently developed long-read technology. This assembly contained over 50 chromosomal regions of phage origin and 17 extrachromosomal elements within the genome, encoding many important traits at the host-microbe interface. Thus, the biology of this symbiont is enabled by a complex array of mobile elements.

Mobile elements—plasmids and phages—are important components of microbial function and evolution via traits that they encode and their capacity to shuttle genetic material between species. We here report the unusually rich array of mobile elements within the genome of Arsenophonus nasoniae, the son-killer symbiont of the parasitic wasp Nasonia vitripennis. This microbe’s genome has the highest prophage complement reported to date, with over 50 genomic regions that represent either intact or degraded phage material. Moreover, the genome is predicted to include 17 extrachromosomal genetic elements, which carry many genes predicted to be important at the microbe-host interface, derived from a diverse assemblage of insect-associated gammaproteobacteria. In our system, this diversity was previously masked by repetitive mobile elements that broke the assembly derived from short reads. These findings suggest that other complex bacterial genomes will be revealed in the era of long-read sequencing.

Ultraplexing: increasing the efficiency of long-read sequencing for hybrid assembly with k-mer-based multiplexing

Hybrid genome assembly has emerged as an important technique in bacterial genomics, but cost and labor requirements limit large-scale application. We present Ultraplexing, a method to improve per-sample sequencing cost and hands-on time of Nanopore sequencing for hybrid assembly by at least 50% compared to molecular barcoding while maintaining high assembly quality. Ultraplexing requires the availability of Illumina data and uses inter-sample genetic variability to assign reads to isolates, which obviates the need for molecular barcoding. Thus, Ultraplexing can enable significant sequencing and labor cost reductions in large-scale bacterial genome projects.

A megaplasmid family driving dissemination of multidrug resistance in Pseudomonas

Multidrug resistance (MDR) represents a global threat to health. Here, we used whole genome sequencing to characterise Pseudomonas aeruginosa MDR clinical isolates from a hospital in Thailand. Using long-read sequence data we obtained complete sequences of two closely related megaplasmids (>420 kb) carrying large arrays of antibiotic resistance genes located in discrete, complex and dynamic resistance regions, and revealing evidence of extensive duplication and recombination events. A comprehensive pangenomic and phylogenomic analysis indicates that: 1) these large plasmids comprise an emerging family present in different members of the Pseudomonas genus, and associated with multiple sources (geographical, clinical or environmental); 2) the megaplasmids encode diverse niche-adaptive accessory traits, including multidrug resistance; 3) the accessory genome of the megaplasmid family is highly flexible and diverse. The history of the megaplasmid family, inferred from our analysis of the available database, suggests that members carrying multiple resistance genes date back to at least the 1970s.

A Proteogenomic Resource Enabling Integrated Analysis of Listeria Genotype-Proteotype-Phenotype Relationships

Listeria monocytogenes is an opportunistic foodborne pathogen responsible for listeriosis, a potentially fatal foodborne disease. Many different Listeria strains and serotypes exist, but a proteogenomic resource that bridges the gap in our molecular understanding of the relationships between the Listeria genotypes and phenotypes via proteotypes is still missing. Here, we devised a next-generation proteogenomics strategy that enables the community to rapidly proteotype Listeria strains and relate this information back to the genotype. Based on sequencing and de novo assembly of the two most commonly used Listeria model strains, EGD-e and ScottA, we established two comprehensive Listeria proteogenomic databases. A genome comparison established core- and strain-specific genes potentially responsible for virulence differences. Next, we established a DIA/SWATH-based proteotyping strategy, including a new and robust sample preparation workflow, that enables the reproducible, sensitive, and relative quantitative measurement of Listeria proteotypes. This reusable and publicly available DIA/SWATH library covers 70% of open reading frames of Listeria and represents the most extensive spectral library for Listeria proteotype analysis to date. We used these two new resources to investigate the Listeria proteotype in states mimicking the upper gastrointestinal passage. Exposure of Listeria to bile salts at 37 °C, which simulates conditions encountered in the duodenum, showed significant proteotype perturbations including an increase of FlaA, the structural protein of flagella. Given that Listeria is known to lose its flagella above 30 °C, this was an unexpected finding. The formation of flagella, which might have implications on infectivity, was validated by parallel reaction monitoring and light and scanning electron microscopy. flaA transcript levels did not change significantly upon exposure to bile salts at 37 °C, suggesting regulation at the post-transcriptional level. Together, these analyses provide a comprehensive proteogenomic resource and toolbox for the Listeria community enabling the analysis of Listeria genotype-proteotype-phenotype relationships.

Intragenomic polymorphisms in the ITS region of high-quality genomes of the Hypoxylaceae (Xylariales, Ascomycota)

The internal transcribed spacer (ITS) region of the ribosomal DNA (rDNA) has been established (and is generally accepted) as a primary “universal” genetic barcode for fungi for many years, but the actual value for taxonomy has been heavily disputed among mycologists. Recently, twelve draft genome sequences, mainly derived from type species of the family Hypoxylaceae (Xylariales, Ascomycota) and the ex-epitype strain of Xylaria hypoxylon have become available during the course of a large phylogenomic study that was primarily aimed at establishing a correlation between the existing multi-gene-based genealogy with a genome-based phylogeny and the discovery of novel biosynthetic gene clusters encoding for secondary metabolites. The genome sequences were obtained using combinations of Illumina and Oxford nanopore technologies or PacBio sequencing, respectively, and resulted in high-quality sequences with an average N50 of 3.2 Mbp. While the main results will be published concurrently in a separate paper, the current case study was dedicated to the detection of ITS nrDNA copies in the genomes, in an attempt to explain certain incongruities and apparent mismatches between phenotypes and genotypes that had been observed during previous polyphasic studies. The results revealed that all of the studied strains had at least three copies of rDNA in their genomes, with Hypoxylon fragiforme having at least 19 copies of the ITS region, followed by Xylaria hypoxylon with at least 13 copies. Several of the genomes contained 2–3 copies that were nearly identical, but in some cases drastic differences, below 97% identity were observed. In one case, ascribable to the presence of a pseudogene, the deviations of the ITS sequences from the same genome resulted in only ca. 90% of overall homology. These results are discussed in the scope of the current trends to use ITS data for species recognition and segregation of fungi. We propose that additional genomes should be checked for such ITS polymorphisms to reassess the validity of this non-coding part of the fungal DNA for molecular identification.

Extensive Genomic Rearrangements along with Distinct Mobilome and TALome are Associated with Extreme Pathotypes of a Rice Pathogen

Xanthomonas oryzae pv. oryzae (Xoo) is a serious pathogen of rice which displays tremendous interstrain variation. The emergence of highly-virulent strains of Xoo is a major threat to rice cultivation. Evolutionary insights into genome dynamics of highly virulent strains as compared with the less-virulent ones are crucial for understanding the molecular basis of exceptional success of Xoo as a highly evolved plant pathogen. In the present study, we report complete genome sequence of Xoo strains with extreme-virulent pathotypes (XVPs) characterized based on their reaction toward ten resistance (Xa) genes. One strain, IXO1088, can overcome resistance mediated by all the ten resistance genes while the other strain IXO704 cannot overcome any of them. Interestingly, our investigation revealed that XVPs display dramatic variation in the genome structure with numerous rearrangements/inversions. Moreover, XVPs also possess distinct transposon content and prophage elements that may provide genomic flux required for the acquisition of novel gene cassettes and structural changes in the genome. Interestingly, analysis of transcription activator-like effector proteins, which are major virulence determinants of Xanthomonas pathogen show marked variation in the transcription activator-like effector content and DNA binding domain of tal genes. Overall, the present study indicates the possible role of mobilomes and repetitive elements in major structural and sequence alterations, which may be leading to the emergence of novel and extreme pathotypes. The knowledge and resource of XVPs will be invaluable in the further systematic understanding of evolution and management of variant pathotypes of Xoo.

High quality 3C de novo assembly and annotation of a multidrug resistant ST-111 Pseudomonas aeruginosa genome: Benchmark of hybrid and non-hybrid assemblers

Genotyping methods and genome sequencing are indispensable to reveal genomic structure of bacterial species displaying high level of genome plasticity. However, reconstruction of genome or assembly is not straightforward due to data complexity, including repeats, mobile and accessory genetic elements of bacterial genomes. Moreover, since the solution to this problem is strongly influenced by sequencing technology, bioinformatics pipelines, and selection criteria to assess assemblers, there is no systematic way to select a priori the optimal assembler and parameter settings. To assembly the genome of Pseudomonas aeruginosa strain AG1 (PaeAG1), short reads (Illumina) and long reads (Oxford Nanopore) sequencing data were used in 13 different non-hybrid and hybrid approaches. PaeAG1 is a multiresistant high-risk sequence type 111 (ST-111) clone that was isolated from a Costa Rican hospital and it was the first report of an isolate of P. aeruginosa carrying both blaVIM-2 and blaIMP-18 genes encoding for metallo-β-lactamases (MBL) enzymes. To assess the assemblies, multiple metrics regard to contiguity, correctness and completeness (3C criterion, as we define here) were used for benchmarking the 13 approaches and select a definitive assembly. In addition, annotation was done to identify genes (coding and RNA regions) and to describe the genomic content of PaeAG1. Whereas long reads and hybrid approaches showed better performances in terms of contiguity, higher correctness and completeness metrics were obtained for short read only and hybrid approaches. A manually curated and polished hybrid assembly gave rise to a single circular sequence with 100% of core genes and known regions identified, >98% of reads mapped back, no gaps, and uniform coverage. The strategy followed to obtain this high-quality 3C assembly is detailed in the manuscript and we provide readers with an all-in-one script to replicate our results or to apply it to other troublesome cases. The final 3C assembly revealed that the PaeAG1 genome has 7,190,208 bp, a 65.7% GC content and 6,709 genes (6,620 coding sequences), many of which are included in multiple mobile genomic elements, such as 57 genomic islands, six prophages, and two complete integrons with blaVIM-2 and blaIMP-18 MBL genes. Up to 250 and 60 of the predicted genes are anticipated to play a role in virulence (adherence, quorum sensing and secretion) or antibiotic resistance (β-lactamases, efflux pumps, etc). Altogether, the assembly and annotation of the PaeAG1 genome provide new perspectives to continue studying the genomic diversity and gene content of this important human pathogen.

Tandem repeats lead to sequence assembly errors and impose multi-level challenges for genome and protein databases

The widespread occurrence of repetitive stretches of DNA in genomes of organisms across the tree of life imposes fundamental challenges for sequencing, genome assembly, and automated annotation of genes and proteins. This multi-level problem can lead to errors in genome and protein databases that are often not recognized or acknowledged. As a consequence, end users working with sequences with repetitive regions are faced with 'ready-to-use' deposited data whose trustworthiness is difficult to determine, let alone to quantify. Here, we provide a review of the problems associated with tandem repeat sequences that originate from different stages during the sequencing-assembly-annotation-deposition workflow, and that may proliferate in public database repositories affecting all downstream analyses. As a case study, we provide examples of the Atlantic cod genome, whose sequencing and assembly were hindered by a particularly high prevalence of tandem repeats. We complement this case study with examples from other species, where mis-annotations and sequencing errors have propagated into protein databases. With this review, we aim to raise the awareness level within the community of database users, and alert scientists working in the underlying workflow of database creation that the data they omit or improperly assemble may well contain important biological information valuable to others.