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

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

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system.

BACKGROUND: Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii. RESULTS: Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats. CONCLUSIONS: Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

Seasonal Dynamics of Glyphosate and AMPA in Lake Greifensee: Rapid Microbial Degradation in the Epilimnion During Summer.

Occurrence and fate of glyphosate, a widely used herbicide, and its main metabolite AMPA was investigated in Lake Greifensee, Switzerland. Monthly vertical concentration profiles in the lake showed an increase of glyphosate concentrations in the epilimnion from 15 ng/L in March to 145 ng/L in July, followed by a sharp decline to <5 ng/L in August. A similar pattern was observed for AMPA. Concentrations of glyphosate and AMPA in the two main tributaries generally were much higher than in the lake. Simulations using a numerical lake model indicated that a substantial amount of glyphosate and AMPA dissipated in the epilimnion, mainly in July and August, with half-lives of only ≈2-4 days which is ≫100 times faster than in the preceding months. Fast dissipation coincided with high water temperatures and phytoplankton densities, and low phosphate concentrations. This indicates that glyphosate might have been used as an alternative phosphorus source by bacterio- and phytoplankton. Metagenomic analysis of lake water revealed the presence of organisms known to be capable of degrading glyphosate and AMPA.

Erwinia gerundensis sp. nov., a cosmopolitan epiphyte originally isolated from pome fruit trees.

A survey to obtain potential antagonists of pome fruit tree diseases yielded two yellow epiphytic bacterial isolates morphologically similar to Pantoea agglomerans, but showing no biocontrol activity. Whole-cell MALDI-TOF mass spectrometry and analysis of 16S rRNA gene and gyrB sequences suggested the possibility of a novel species with a phylogenetic position in either the genus Pantoea or the genus Erwinia. Multi-locus sequence analysis (MLSA) placed the two strains in the genus Erwinia and supported their classification as a novel species. The strains showed general phenotypic characteristics typical of this genus and results of DNA-DNA hybridizations confirmed that they represent a single novel species. Both strains showed a DNA G+C content, as determined by HPLC, of 54.5 mol% and could be discriminated from phylogenetically related species of the genus Erwinia by their ability to utilize potassium gluconate, potassium 2-ketogluconate, maltose, melibiose and raffinose. Whole-genome sequencing of strain EM595T revealed the presence of a chromosomal carotenoid biosynthesis gene cluster similar to those found in species of the genera Cronobacter and Pantoea that explains the pigmentation of the strain, which is atypical for the genus Erwinia. Additional strains belonging to the same species were recovered from different plant hosts in three different continents, revealing the cosmopolitan nature of this epiphyte. The name Erwinia gerundensis sp. nov. is proposed, with EM595T ( = LMG 28990T = CCOS 903T) as the designated type strain.

Development and Validation of LNA-Based Quantitative Real-Time PCR Assays for Detection and Identification of the Root-Knot Nematode Meloidogyne enterolobii in Complex DNA Backgrounds.

Meloidogyne enterolobii is a quarantine root-knot nematode posing a major threat to agricultural production systems worldwide. It attacks many host plants, including important agricultural crops, ornamentals, and trees. M. enterolobii is a highly virulent and pathogenic root-knot nematode species, able to reproduce on plants resistant to other Meloidogyne spp. Significant crop damage has been reported in Asia, South America, Africa, the United States, France, and greenhouses in Switzerland. To identify potential introduction pathways and ensure appropriate phytosanitary measures and management strategies, accurate detection and identification tools are needed. Therefore, two real-time quantitative polymerase chain reaction (PCR) assays based on the second intergenic spacer region of the ribosomal DNA cistron and the cytochrome oxidase c subunit I (COI) gene using locked nucleic acid probes were developed and validated for fast and reliable detection and identification of M. enterolobii. Analytical specificity was confirmed with 16 M. enterolobii populations, 16 populations of eight closely related Meloidogyne spp., and four species from other nematode genera. Optimizing and testing the assays on two real-time PCR platforms revealed an analytical sensitivity of one juvenile in a background of 1,000 nematodes and the intended limit of detection of one juvenile per 100 ml of soil. Both assays performed equally well, with the COI-based assay showing a slightly better performance concerning detection of M. enterolobii target DNA in complex DNA backgrounds.

Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites.

Erwinia amylovora causes a major disease of pome fruit trees worldwide, and is regulated as a quarantine organism in many countries. While some diversity of isolates has been observed, molecular epidemiology of this bacterium is hindered by a lack of simple molecular typing techniques with sufficiently high resolution. We report a molecular typing system of E. amylovora based on variable number of tandem repeats (VNTR) analysis. Repeats in the E. amylovora genome were identified with comparative genomic tools, and VNTR markers were developed and validated. A Multiple-Locus VNTR Analysis (MLVA) was applied to E. amylovora isolates from bacterial collections representing global and regional distribution of the pathogen. Based on six repeats, MLVA allowed the distinction of 227 haplotypes among a collection of 833 isolates of worldwide origin. Three geographically separated groups were recognized among global isolates using Bayesian clustering methods. Analysis of regional outbreaks confirmed presence of diverse haplotypes but also high representation of certain haplotypes during outbreaks. MLVA analysis is a practical method for epidemiological studies of E. amylovora, identifying previously unresolved population structure within outbreaks. Knowledge of such structure can increase our understanding on how plant diseases emerge and spread over a given geographical region.

Genome sequence data of the strongly antagonistic yeast Pichia kluyveri isolate APC 11.10 B as a foundation for analysing biocontrol mechanisms.

Pichia kluyveri strain APC 11.10 B was isolated from apple bark in Switzerland and exhibited strong antagonistic activity against plant pathogenic fungi in vitro (e.g., Botrytis, Fusarium or Monilinia isolates). In order to identify the mechanisms underlying this antagonism, we have sequenced the genome of this isolate by long- and short-read sequencing technologies. The sequence data were de novo assembled into nine scaffolds and a fully resolved circularized mitogenome. The total genome size was 10.9 Mbp and 7451 potential open reading frames (ORFs) and 202 tRNA genes were predicted. In comparison to two P. kluyveri genomes deposited at the NCBI (of strains X31-10 and CBA6002), the APC 11.10 B strain seemed to represent a hybrid because backmapping of sequencing reads resulted in a high rate of heterozygous and structural variants in the nuclear genome (this was not observed for the mitochondrial genome). The P. kluyveri (APC 11.10 B) draft genome represents a first step and resource for genome mining, comparative and functional genomics (e.g., identifying the biocontrol mode of action), and evolutionary studies. Since the genus Pichia comprises many biotechnologically relevant yeasts, the genome data may be used in a variety of fields and disciplines.

Complete genome sequence of clinical isolate Pantoea ananatis LMG 5342.

The enterobacterium Pantoea ananatis is an ecologically versatile species. It has been found in the environment, as plant epiphyte and endophyte, as an emerging phytopathogen, and as a presumptive, opportunistic human pathogen. Here, we report the complete genome sequence of P. ananatis LMG 5342, isolated from a human wound.

Next generation biosecurity: Towards genome based identification to prevent spread of agronomic pests and pathogens using nanopore sequencing.

The unintentional movement of agronomic pests and pathogens is steadily increasing due to the intensification of global trade. Being able to identify accurately and rapidly early stages of an invasion is critical for developing successful eradication or management strategies. For most invasive organisms, molecular diagnostics is today the method of choice for species identification. However, the currently implemented tools are often developed for certain taxa and need to be adapted for new species, making them ill-suited to cope with the current constant increase in new invasive species. To alleviate this impediment, we developed a fast and accurate sequencing tool allowing to modularly obtain genetic information at different taxonomical levels. Using whole genome amplification (WGA) followed by Oxford nanopore MinION sequencing, our workflow does not require any a priori knowledge on the investigated species and its classification. While mainly focusing on harmful plant pathogenic insects, we also demonstrate the suitability of our workflow for the molecular identification of bacteria (Erwinia amylovora and Escherichia coli), fungi (Cladosporium herbarum, Colletotrichum salicis, Neofabraea alba) and nematodes (Globodera rostochiensis). On average, the pairwise identity between the generated consensus sequences and best GenBank BLAST matches was 99.6 ± 0.6%. Additionally, assessing the generated insect genomic dataset, the potential power of the workflow to detect pesticide resistance genes, as well as arthropod-infecting viruses and endosymbiotic bacteria is demonstrated.

Different gut microbial communities correlate with efficacy of albendazole-ivermectin against soil-transmitted helminthiases.

Soil-transmitted helminth infections represent a large burden with over a quarter of the world's population at risk. Low cure rates are observed with standard of care (albendazole); therefore, a more effective combination therapy (albendazole and ivermectin) is being investigated but showed variable treatment efficacies without evidence of intrinsic parasite resistance. Here, we analyzed the microbiome of Trichuris trichiura and hookworm-infected patients and found an association of different enterotypes with treatment efficacy. 80 T. trichiura-infected patients with hookworm co-infections from Pak-Khan, Laos, received either albendazole (n = 41) or albendazole and ivermectin combination therapy (n = 39). Pre-/post-treatment stool samples were collected to monitor treatment efficacy and microbial communities were profiled using 16S rRNA gene sequencing, qPCR, and shotgun sequencing. We identified three bacterial enterotypes and show that pre-treatment enterotype is associated with efficacy of the combination treatment for both T. trichiura (CRET1 = 5.8%; CRET2 = 16.6%; CRET3 = 68.8%) and hookworm (CRET1 = 31.3%; CRET2 = 16.6%; CRET3 = 78.6%). This study shows that pre-treatment enterotype enables predicting treatment outcome of combination therapy for T. trichiura and hookworm infections.Trial registration: ClinicalTrials.gov, NCT03527732. Registered 17 May 2018, https://clinicaltrials.gov/ct2/show/NCT03527732 .

Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families.

Long-read DNA sequencing technologies require high molecular weight (HMW) DNA of adequate purity and integrity, which can be difficult to isolate from plant material. Plant leaves usually contain high levels of carbohydrates and secondary metabolites that can impact DNA purity, affecting downstream applications. Several protocols and kits are available for HMW DNA extraction, but they usually require a high amount of input material and often lead to substantial DNA fragmentation, making sequencing suboptimal in terms of read length and data yield. We here describe a protocol for plant HMW DNA extraction from low input material (0.1 g) which is easy to follow and quick (2.5 h). This method successfully enabled us to extract HMW from four species from different families (Orchidaceae, Poaceae, Brassicaceae, Asteraceae). In the case of recalcitrant species, we show that an additional purification step is sufficient to deliver a clean DNA sample. We demonstrate the suitability of our protocol for long-read sequencing on the Oxford Nanopore Technologies PromethION® platform, with and without the use of a short fragment depletion kit.

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.

Genome sequence of an Erwinia amylovora strain with pathogenicity restricted to Rubus plants.

Here, we present the genome of a strain of Erwinia amylovora, the fire blight pathogen, with pathogenicity restricted to Rubus spp. Comparative genomics of ATCC BAA-2158 with E. amylovora strains from non-Rubus hosts identified significant genetic differences but support the inclusion of this strain within the species E. amylovora.

Dynamics of the Apple Fruit Microbiome after Harvest and Implications for Fruit Quality.

The contribution of the apple microbiome to the production chain of apple was so far largely unknown. Here, we describe the apple fruit microbiome and influences on its composition by parameters such as storage season, storage duration, storage technology, apple variety, and plant protection schemes. A combined culturing and metabarcoding approach revealed significant differences in the abundance, composition, and diversity of the apple fruit microbiome. We showed that relatively few genera contribute a large portion of the microbiome on fruit and that the fruit microbiome changes during the storage season depending on the storage conditions. In addition, we show that the plant protection regime has an influence on the diversity of the fruit microbiome and on the dynamics of pathogenic fungal genera during the storage season. For the genus Neofabraea, the quantitative results from the metabarcoding approach were validated with real-time PCR. In conclusion, we identified key parameters determining the composition and temporal changes of the apple fruit microbiome, and the main abiotic driving factors of microbiome diversity on apple fruit were characterized.

Ultralong Oxford Nanopore Reads Enable the Development of a Reference-Grade Perennial Ryegrass Genome Assembly.

Despite the progress made in DNA sequencing over the last decade, reconstructing telomere-to-telomere genome assemblies of large and repeat-rich eukaryotic genomes is still difficult. More accurate basecalls or longer reads could address this issue, but no current sequencing platform can provide both simultaneously. Perennial ryegrass (Lolium perenne L.) is an example of an important species for which the lack of a reference genome assembly hindered a swift adoption of genomics-based methods into breeding programs. To fill this gap, we optimized the Oxford Nanopore Technologies' sequencing protocol, obtaining sequencing reads with an N50 of 62 kb-a very high value for a plant sample. The assembly of such reads produced a highly complete (2.3 of 2.7 Gb), correct (QV 45), and contiguous (contig N50 and N90 11.74 and 3.34 Mb, respectively) genome assembly. We show how read length was key in determining the assembly contiguity. Sequence annotation revealed the dominance of transposable elements and repeated sequences (81.6% of the assembly) and identified 38,868 protein coding genes. Almost 90% of the bases could be anchored to seven pseudomolecules, providing the first high-quality haploid reference assembly for perennial ryegrass. This protocol will enable producing longer Oxford Nanopore Technology reads for more plant samples and ushering forage grasses into modern genomics-assisted breeding programs.

Genome sequence of the biocontrol agent Pantoea vagans strain C9-1.

Pantoea vagans is a Gram-negative enterobacterial plant epiphyte of a broad range of plants. Here we report the 4.89-Mb genome sequence of P. vagans strain C9-1 (formerly Pantoea agglomerans), which is commercially registered for biological control of fire blight, a disease of pear and apple trees caused by Erwinia amylovora.

Complete genome sequence of the fire blight pathogen Erwinia amylovora CFBP 1430 and comparison to other Erwinia spp.

Fire blight, caused by the enterobacterium Erwinia amylovora, is a devastating disease of rosaceous plants that has global economic importance for apple and pear production and trade. The complete genome of E. amylovora CFBP 1430 was sequenced, annotated, and compared with the genomes of other Erwinia spp. Several singleton and shared features of the E. amylovora CFBP 1430 genome were identified that offer a first view into evolutionary aspects within the genus Erwinia. Comparative genomics identified or clarified virulence and fitness determinants and secretion systems. Novel insights revealed in the genome of E. amylovora CFBP 1430 hold potential for exploitation to improve the design of more effective fire blight control strategies.

Complete genome sequence of the fire blight pathogen Erwinia pyrifoliae DSM 12163T and comparative genomic insights into plant pathogenicity.

BACKGROUND: Erwinia pyrifoliae is a newly described necrotrophic pathogen, which causes fire blight on Asian (Nashi) pear and is geographically restricted to Eastern Asia. Relatively little is known about its genetics compared to the closely related main fire blight pathogen E. amylovora. RESULTS: The genome of the type strain of E. pyrifoliae strain DSM 12163T, was sequenced using both 454 and Solexa pyrosequencing and annotated. The genome contains a circular chromosome of 4.026 Mb and four small plasmids. Based on their respective role in virulence in E. amylovora or related organisms, we identified several putative virulence factors, including type III and type VI secretion systems and their effectors, flagellar genes, sorbitol metabolism, iron uptake determinants, and quorum-sensing components. A deletion in the rpoS gene covering the most conserved region of the protein was identified which may contribute to the difference in virulence/host-range compared to E. amylovora. Comparative genomics with the pome fruit epiphyte Erwinia tasmaniensis Et1/99 showed that both species are overall highly similar, although specific differences were identified, for example the presence of some phage gene-containing regions and a high number of putative genomic islands containing transposases in the E. pyrifoliae DSM 12163T genome. CONCLUSIONS: The E. pyrifoliae genome is an important addition to the published genome of E. tasmaniensis and the unfinished genome of E. amylovora providing a foundation for re-sequencing additional strains that may shed light on the evolution of the host-range and virulence/pathogenicity of this important group of plant-associated bacteria.

Broad spectrum microarray for fingerprint-based bacterial species identification.

BACKGROUND: Microarrays are powerful tools for DNA-based molecular diagnostics and identification of pathogens. Most target a limited range of organisms and are based on only one or a very few genes for specific identification. Such microarrays are limited to organisms for which specific probes are available, and often have difficulty discriminating closely related taxa. We have developed an alternative broad-spectrum microarray that employs hybridisation fingerprints generated by high-density anonymous markers distributed over the entire genome for identification based on comparison to a reference database. RESULTS: A high-density microarray carrying 95,000 unique 13-mer probes was designed. Optimized methods were developed to deliver reproducible hybridisation patterns that enabled confident discrimination of bacteria at the species, subspecies, and strain levels. High correlation coefficients were achieved between replicates. A sub-selection of 12,071 probes, determined by ANOVA and class prediction analysis, enabled the discrimination of all samples in our panel. Mismatch probe hybridisation was observed but was found to have no effect on the discriminatory capacity of our system. CONCLUSIONS: These results indicate the potential of our genome chip for reliable identification of a wide range of bacterial taxa at the subspecies level without laborious prior sequencing and probe design. With its high resolution capacity, our proof-of-principle chip demonstrates great potential as a tool for molecular diagnostics of broad taxonomic groups.

Genotypic comparison of Pantoea agglomerans plant and clinical strains.

BACKGROUND: Pantoea agglomerans strains are among the most promising biocontrol agents for a variety of bacterial and fungal plant diseases, particularly fire blight of apple and pear. However, commercial registration of P. agglomerans biocontrol products is hampered because this species is currently listed as a biosafety level 2 (BL2) organism due to clinical reports as an opportunistic human pathogen. This study compares plant-origin and clinical strains in a search for discriminating genotypic/phenotypic markers using multi-locus phylogenetic analysis and fluorescent amplified fragment length polymorphisms (fAFLP) fingerprinting. RESULTS: Majority of the clinical isolates from culture collections were found to be improperly designated as P. agglomerans after sequence analysis. The frequent taxonomic rearrangements underwent by the Enterobacter agglomerans/Erwinia herbicola complex may be a major problem in assessing clinical associations within P. agglomerans. In the P. agglomerans sensu stricto (in the stricter sense) group, there was no discrete clustering of clinical/biocontrol strains and no marker was identified that was uniquely associated to clinical strains. A putative biocontrol-specific fAFLP marker was identified only in biocontrol strains. The partial ORF located in this band corresponded to an ABC transporter that was found in all P. agglomerans strains. CONCLUSION: Taxonomic mischaracterization was identified as a major problem with P. agglomerans, and current techniques removed a majority of clinical strains from this species. Although clear discrimination between P. agglomerans plant and clinical strains was not obtained with phylogenetic analysis, a single marker characteristic of biocontrol strains was identified which may be of use in strain biosafety determinations. In addition, the lack of Koch's postulate fulfilment, rare retention of clinical strains for subsequent confirmation, and the polymicrobial nature of P. agglomerans clinical reports should be considered in biosafety assessment of beneficial strains in this species.