Bacterial canker of stone-fruits: Infection experiments with Pseudomonas mors-prunorum and P. syringae

The effect was compared of inoculating Pseudomonas mors-prunorum and P. syringae at high inoculum concentration through wounds in plum stems and cherry branches, and through cherry-leaf scars. In the wound inoculations a North American cherry strain of P. syringae was considerably more virulent than any of three indigenous cherry strains of P. mors-prunorum, but two pear strains of P. syringae were less virulent. P. mors-prunorum showed a much greater capacity to invade through leaf scars, particularly towards the end of the leaf-fall period, when all three P. syringae strains were largely ineffective. Evidence is discussed which suggests that the greater infectivity of the P. mors-prunorum strains through leaf scars was related to their capacity to colonize the host tissues from small inocula. The P. syringae strains died out in cherry cankers earlier than P. mors-prunorum strains and they were less stable during host passage. The results of the experiments indicate that virulence in these organisms is a complex phenomenon and determined by several factors.     

Incongruence between multi‐locus sequence analysis (MLSA) and whole‐genome‐based phylogenies: Pseudomonas syringae pathovar pisi as a cautionary tale

Previous phylogenies, built using a subset of genomic loci, split Pseudomonas syringae pv. pisi into two well‐supported clades and implied convergence in host range for these lineages. The analysis of phenotypic and genotypic data within the context of this phylogenetic relationship implied further convergence at the level of virulence gene loss and acquisition. We generate draft genome assemblies for two additional P. syringae strains, isolated from diseased pea plants, and demonstrate incongruence between phylogenies created from a subset of the data compared with the whole genomes. Our whole‐genome analysis demonstrates that strains classified as pv. pisi actually form a coherent monophyletic clade, so that apparent convergence is actually the product of shared ancestry. We use this example to urge caution when making evolutionary inferences across closely related strains of P. syringae.     

Two native plasmids of Pseudomonas syringae pathovar tomato strain PT23 share a large amount of repeated DNA, including replication sequences

Strain PT23 of Pseudomonas syringae pv, tomato contains four native plasmids, designated A, B, C, and D. By DNA hybridization of genomic and plasmid DNA digests from the wild type and a plasmid-cured strain, we determined that c. 61 kb (c. 74%) of pPT23B is repeated in pPT23A and only c. 17 kb (c. 21%) is in single copy in strain PT23. pPT23B also contains DNA repeated in the chromosome that occurs in three DNA fragments of 0.6, 4.6, and 9.6 kb that might be transposable elements. Additionally, the 9.6 kb fragment also shares sequences with the three other plasmids of strain PT23. By DNA hybridization with the origin of replication from a native plasmid of P. syringae pv. syringae and in vivo replication tests, we identified the origins of replication of plasmids A, B, and D and showed that they cross-hybridize. The putative par region from pPT23 A has also been identified and is not conserved in the other three native plasmids from strain PT23. By using the defined minimal origin of replication from pPT23 A as a probe, we showed that it is highly conserved in 14 strains belonging to nine different pathovars of P. syringae and that as many as five different native plasmids with closely related origins of replication coexist in the same cell. The duplication and reorganization of plasmids might therefore occur at high frequency and could be responsible for the existence of large numbers of native plasmids in P. syringae strains.     

Global genomic analysis of Pseudomonas savastanoi pv. savastanoi plasmids

Pseudomonas savastanoi pv. savastanoi strains harbor native plasmids belonging to the pPT23A plasmid family (PFPs) which are detected in all pathovars of the related species Pseudomonas syringae examined and contribute to the ecological and pathogenic fitness of their host. However, there is a general lack of information about the gene content of P. savastanoi pv. savastanoi plasmids and their role in the interaction of this pathogen with olive plants. We designed a DNA macroarray containing 135 plasmid-borne P. syringae genes to conduct a global genetic analysis of 32 plasmids obtained from 10 P. savastanoi pv. savastanoi strains. Hybridization results revealed that the number of PFPs per strain varied from one to four. Additionally, most strains contained at least one plasmid (designated non-PFP) that did not hybridize to the repA gene of pPT23A. Only three PFPs contained genes involved in the biosynthesis of the virulence factor indole-3-acetic acid (iaaM, iaaH, and iaaL). In contrast, ptz, a gene involved in the biosynthesis of cytokinins, was found in five PFPs and one non-PFP. Genes encoding a type IV secretion system (T4SS), type IVA, were found in both PFPs and non-PFPs; however, type IVB genes were found only on PFPs. Nine plasmids encoded both T4SSs, whereas seven other plasmids carried none of these genes. Most PFPs and non-PFPs hybridized to at least one putative type III secretion system effector gene and to a variety of additional genes encoding known P. syringae virulence factors and one or more insertion sequence transposase genes. These results indicate that non-PFPs may contribute to the virulence and fitness of the P. savastanoi pv. savastanoi host. The overall gene content of P. savastanoi pv. savastanoi plasmids, with their repeated information, mosaic arrangement, and insertion sequences, suggests a possible role in adaptation to a changing environment.     

Is the whole greater than the sum of its parts? De novo assembly strategies for bacterial genomes based on paired-end sequencing

Background

Whole genome sequence construction is becoming increasingly feasible because of advances in next generation sequencing (NGS), including increasing throughput and read length. By simply overlapping paired-end reads, we can obtain longer reads with higher accuracy, which can facilitate the assembly process. However, the influences of different library sizes and assembly methods on paired-end sequencing-based de novo assembly remain poorly understood.

Results

We used 250 bp Illumina Miseq paired-end reads of different library sizes generated from genomic DNA from Escherichia coli DH1 and Streptococcus parasanguinis FW213 to compare the assembly results of different library sizes and assembly approaches. Our data indicate that overlapping paired-end reads can increase read accuracy but sometimes cause insertion or deletions. Regarding genome assembly, merged reads only outcompete original paired-end reads when coverage depth is low, and larger libraries tend to yield better assembly results. These results imply that distance information is the most critical factor during assembly. Our results also indicate that when depth is sufficiently high, assembly from subsets can sometimes produce better results.

Conclusions

In summary, this study provides systematic evaluations of de novo assembly from paired end sequencing data. Among the assembly strategies, we find that overlapping paired-end reads is not always beneficial for bacteria genome assembly and should be avoided or used with caution especially for genomes containing high fraction of repetitive sequences. Because increasing numbers of projects aim at bacteria genome sequencing, our study provides valuable suggestions for the field of genomic sequence construction.

Electronic supplementary material

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

Identifying the inserted locus of randomly integrated expression plasmids by whole-genome sequencing of Aspergillus strains

ABSTRACT

Whole-genome sequencing was conducted on two Aspergillus oryzae strains used for the manufacturing of food enzymes, Acrylaway® and Shearzyme®, with the aim of identifying the inserted locus of randomly integrated expression plasmid and obtaining flanking sequences for safety assessment. Illumina paired-end sequencing was employed, and the obtained reads were mapped to two references: the public genome sequence of Aspergillus oryzae RIB40 and the in-house sequence of the used expression plasmid. Introducing the concept of linking-reads, one locus for each was successfully identified as the integrated site. In the case of Acrylaway®, the obtained sequences suggested that the expression plasmid had been integrated as multiple copies in tandem form. In the case of Shearzyme®, however, information on one edge of the insert was missing, which required extra polymerase chain reaction (PCR) cloning for safety assessment. A 4-kb deletion was detected at the integrated site. There was also evidence of rearrangement occurring in Shearzyme® strain.     

Complete assembly,annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.     

Inhibitory effect of Thymus vulgaris and Origanum vulgare essential oils on virulence factors of phytopathogenic Pseudomonas syringae strains

• Pseudomonas syringae is a phytopathogenic bacterium that causes lesions in leaves during the colonisation process. The damage is associated with production of many virulence factors, such as biofilm and phytotoxins. The essential oils of Thymus vulgaris (thyme) and Origanum vulgare (oregano) have been demonstrated to inhibit P. syringae. The aim of this study was to investigate the effects of T. vulgaris and O. vulgare essential oils on production of virulence factors of phytopathogenic P. syringae strains, including anti‐biofilm and anti‐toxins activities.
• The broth microdilution method was used for determination of MIC and biofilm inhibition assays. Coronatine, syringomycin and tabtoxin were pheno‐ and genotypically evaluated.
• Both oils showed good inhibitory activity against P. syringae, with MIC values from 1.43 to 11.5 mg·ml^?1 for thyme and 5.8 to 11.6 mg·ml^?1 for oregano. Biofilm formation, production of coronatine, syringomycin and tabtoxin were inhibited by thyme and oregano essential oil in most strains.
• The results presented here are promising, demonstrating the bactericidal activity and reduction of virulence factor production after treatment with thyme and oregano oil, providing insight into how they exert their antibacterial activity. These natural products could be considered in the future for the control of diseases caused by P. syringae.

     

Genome-Wide Microsatellite Identification in the Fungus Anisogramma anomala Using Illumina Sequencing and Genome Assembly

High-throughput sequencing has been dramatically accelerating the discovery of microsatellite markers (also known as Simple Sequence Repeats). Both 454 and Illumina reads have been used directly in microsatellite discovery and primer design (the “Seq-to-SSR” approach). However, constraints of this approach include: 1) many microsatellite-containing reads do not have sufficient flanking sequences to allow primer design, and 2) difficulties in removing microsatellite loci residing in longer, repetitive regions. In the current study, we applied the novel “Seq-Assembly-SSR” approach to overcome these constraints in Anisogramma anomala. In our approach, Illumina reads were first assembled into a draft genome, and the latter was then used in microsatellite discovery. A. anomala is an obligate biotrophic ascomycete that causes eastern filbert blight disease of commercial European hazelnut. Little is known about its population structure or diversity. Approximately 26 M 146 bp Illumina reads were generated from a paired-end library of a fungal strain from Oregon. The reads were assembled into a draft genome of 333 Mb (excluding gaps), with contig N₅₀ of 10,384 bp and scaffold N₅₀ of 32,987 bp. A bioinformatics pipeline identified 46,677 microsatellite motifs at 44,247 loci, including 2,430 compound loci. Primers were successfully designed for 42,923 loci (97%). After removing 2,886 loci close to assembly gaps and 676 loci in repetitive regions, a genome-wide microsatellite database of 39,361 loci was generated for the fungus. In experimental screening of 236 loci using four geographically representative strains, 228 (96.6%) were successfully amplified and 214 (90.7%) produced single PCR products. Twenty-three (9.7%) were found to be perfect polymorphic loci. A small-scale population study using 11 polymorphic loci revealed considerable gene diversity. Clustering analysis grouped isolates of this fungus into two clades in accordance with their geographic origins. Thus, the “Seq-Assembly-SSR” approach has proven to be a successful one for microsatellite discovery.     

Genome Sequence of Lactobacillus pentosus KCA1: Vaginal Isolate from a Healthy Premenopausal Woman

The vaginal microbiota, in particular Lactobacillus species, play an important role in female health through modulation of immunity, countering pathogens and maintaining a pH below 4.7. We report the isolation and genome sequence of Lactobacillus pentosus strain KCA1 (formally known as L. plantarum) from the vagina of a healthy Nigerian woman. The genome was sequenced using Illumina GA II technology. The resulting 16,920,226 paired-end reads were assembled with the Velvet tool. Contigs were annotated using the RAST server, and manually curated. A comparative analysis with the available genomes of L. pentosus IG1 and L. plantarum WCFS1 showed that over 15% of the predicted functional activities are found only in this strain. The strain has a chromosome sequence of 3,418,159 bp with a G+C content of 46.4%, and is devoid of plasmids. Novel gene clusters or variants of known genes relative to the reference genomes were found. In particular, the strain has loci encoding additional putative mannose phosphotransferase systems. Clusters of genes include those for utilization of hydantoin, isopropylmalate, malonate, rhamnosides, and genes for assimilation of polyglycans, suggesting the metabolic versatility of L. pentosus KCA1. Loci encoding putative phage defense systems were also found including clustered regularly interspaced short palindromic repeats (CRISPRs), abortive infection (Abi) systems and toxin-antitoxin systems (TA). A putative cluster of genes for biosynthesis of a cyclic bacteriocin precursor, here designated as pentocin KCA1 (penA) were identified. These findings add crucial information for understanding the genomic and geographic diversity of vaginal lactobacilli.     

Methods of assessment of transfer of a gentamycin-resistance gene (aacC1) from genetically engineered microorganisms into agriculturally important soil bacteria

In order to assess the risk associated with the deliberate release of genetically engineered microorganisms (GEMs) into the agricultural environment, the transfer of plasmids between bacterial strains was investigated under laboratory conditions. Genetically modified Rhizobium leguminosarum and Agrobacterium tumefaciens strains carrying the gentamycin acetyltransferase resistance gene (aacC1) on various plasmids were investigated for their ability to transfer the aacC1 gene to their wild-type (w.t.) counterparts, as well as to Pseudomonas syringae. Conjugation experiments between the various strains, were carried out after the relevant characteristics and conditions for selective growth of each bacterial strain had been ascertained. After conjugations on filters had been completed, the putative transconjugants were grown in media containing antibiotics and assessed for the presence of aacC1 gene by: (a) DNA plasmid profile; (b) expression of AAC(3)-I enzyme activity; (c) colony hybridization using a ³²P-labelled DNA probe complementary to the aacC1 gene. The results obtained indicate that transfer of the aacC1 gene from genetically modified strains of R. leguminosarum into a plasmid-free strain of A. tumefaciens occurred via self-transmissible plasmids. Alternatively, genetically modified A. tumefaciens bearing the aacC1 gene on plasmids acquired from R. leguminosarum strains, transferred it ineffectively to a hardly detectable frequency. No transfer of the aacC1 gene from genetically modified R. leguminosarum or A. tumefaciens strains into P. syringae has been observed. These data indicate that in the absence of the RP4 element, genetically modified A. tumefaciens is not able to efficiently transfer aacC1 into w.t. R. leguminosarum and P. syringae. Correspondence to: A. S. Tsiftsoglou     

An integrative plasmid and multiple-sized plasmids of Pseudomonas syringae pv. phaseolicola have extensive homology

Summary Plasmids isolated from five strains of the bean pathogen Pseudomonas syringae pv. phaseolicola were characterized by restriction endonuclease and filter hybridization analyses. BamHI and EcoRI restriction patterns revealed that total plasmid DNA from each strain had a high level of sequence homology with pMC7105, a 148 kbp integrative plasmid found in a sixth strain. Only six BamHI fragments from the eight plasmids in these strains failed to hybridize with pMC7105 probe. Four of these fragments, three from pPP6520 and one from pPP6525 of strain PP652, hybridized strongly to plasmid DNA from a closely-related pathovar, P. syringae pv. glycinea. BamHI fragment 8, which is involved in the integration of pMC7105 into the host chromosome, contains a repeat sequence that was present on all the plasmids except pPP6120 (6.8 kbp), pPP6310 (40 kbp) and pPP6520 (45 kbp). Every plasmid but pPP6520 had fragments that showed weak hybridization to the small plasmid, pPP6120. This homology suggests that a second repetitive sequence is common to these plasmids. The large plasmids (148 to 151 kbp) were essentially identical to pMC7105. The intermediate plasmids (122 to 128 kbp) appeared to be derived mainly from pMC7105 or a related plasmid, whereas the smaller plasmids (6.8 to 45 kbp) appear to have been derived in part from sequences not present in pMC7105.     

Polishing the Oxford Nanopore long-read assemblies of bacterial pathogens with Illumina short reads to improve genomic analyses

Oxford Nanopore sequencing has been widely used to achieve complete genomes of bacterial pathogens. However, the error rates of Oxford Nanopore long reads are high. Various polishing algorithms using Illumina short reads to correct the errors in Oxford Nanopore long-read assemblies have been developed. The impact of polishing the Oxford Nanopore long-read assemblies of bacterial pathogens with Illumina short reads on improving genomic analyses was evaluated using both simulated and real reads. Ten species (10 strains) were selected for simulated reads, while real reads were tested on 11 species (11 strains). Oxford Nanopore long reads were assembled with Unicycler to produce a draft assembly, followed by three rounds of polishing with Illumina short reads using two polishing tools, Pilon and NextPolish. One round of NextPolish polishing generated genome completeness and accuracy parameters similar to the reference genomes, whereas two or three rounds of Pilon polishing were needed, though contiguity remained unchanged after polishing. The polished assemblies of Escherichia coli O157:H7, Salmonella Typhimurium, and Cronobacter sakazakii with simulated reads did not provide accurate plasmid identifications. One round of NextPolish polishing was needed for accurately identifying plasmids in Staphylococcus aureus and E. coli O26:H11 with real reads, whereas one and two rounds of Pilon polishing were necessary for these two strains, respectively. Polishing failed to provide an accurate antimicrobial resistance (AMR) genotype for S. aureus with real reads. One round of polishing recovered an accurate AMR genotype for Klebsiella pneumoniae with real reads. The reference genome and draft assembly of Citrobacter braakii with real reads differed, which carried blaCMY-83 and fosA6, respectively, while both genes were present after one round of polishing. However, polishing did not improve the assembly of E. coli O26:H11 with real reads to achieve numbers of virulence genes similar to the reference genome. The draft and polished assemblies showed a phylogenetic tree topology comparable with the reference genomes. For multilocus sequence typing and pan-genome analyses, one round of NextPolish polishing was sufficient to obtain accurate results, while two or three rounds of Pilon polishing were needed. Overall, NextPolish outperformed Pilon for polishing the Oxford Nanopore long-read assemblies of bacterial pathogens, though both polishing strategies improved genomic analyses compared to the draft assemblies.     

Sequencing and Genetic Variation of Multidrug Resistance Plasmids in Klebsiella pneumoniae

Background

The development of multidrug resistance is a major problem in the treatment of pathogenic microorganisms by distinct antimicrobial agents. Characterizing the genetic variation among plasmids from different bacterial species or strains is a key step towards understanding the mechanism of virulence and their evolution.

Results

We applied a deep sequencing approach to 206 clinical strains of Klebsiella pneumoniae collected from 2002 to 2008 to understand the genetic variation of multidrug resistance plasmids, and to reveal the dynamic change of drug resistance over time. First, we sequenced three plasmids (70 Kb, 94 Kb, and 147 Kb) from a clonal strain of K. pneumoniae using Sanger sequencing. Using the Illumina sequencing technology, we obtained more than 17 million of short reads from two pooled plasmid samples. We mapped these short reads to the three reference plasmid sequences, and identified a large number of single nucleotide polymorphisms (SNPs) in these pooled plasmids. Many of these SNPs are present in drug-resistance genes. We also found that a significant fraction of short reads could not be mapped to the reference sequences, indicating a high degree of genetic variation among the collection of K. pneumoniae isolates. Moreover, we identified that plasmid conjugative transfer genes and antibiotic resistance genes are more likely to suffer from positive selection, as indicated by the elevated rates of nonsynonymous substitution.

Conclusion

These data represent the first large-scale study of genetic variation in multidrug resistance plasmids and provide insight into the mechanisms of plasmid diversification and the genetic basis of antibiotic resistance.     

Comparative genomics reveals genes significantly associated with woody hosts in the plant pathogen Pseudomonas syringae

The diversification of lineages within Pseudomonas syringae has involved a number of adaptive shifts from herbaceous hosts onto various species of tree, resulting in the emergence of highly destructive diseases such as bacterial canker of kiwi and bleeding canker of horse chestnut. This diversification has involved a high level of gene gain and loss, and these processes are likely to play major roles in the adaptation of individual lineages onto their host plants. In order to better understand the evolution of P. syringae onto woody plants, we have generated de novo genome sequences for 26 strains from the P. syringae species complex that are pathogenic on a range of woody species, and have looked for statistically significant associations between gene presence and host type (i.e. woody or herbaceous) across a phylogeny of 64 strains. We have found evidence for a common set of genes associated with strains that are able to colonize woody plants, suggesting that divergent lineages have acquired similarities in genome composition that may form the genetic basis of their adaptation to woody hosts. We also describe in detail the gain, loss and rearrangement of specific loci that may be functionally important in facilitating this adaptive shift. Overall, our analyses allow for a greater understanding of how gene gain and loss may contribute to adaptation in P. syringae.     

Different impacts of pMP19 on the virulence of Melissococcus plutonius strains with different genetic backgrounds

Virulence factors responsible for bacterial pathogenicity are often encoded by plasmids. In Melissococcus plutonius, the causative agent of European foulbrood of honey bees, a putative virulence plasmid (pMP19) possessing mtxA, which encodes a putative insecticidal toxin, was found by comparative genome analyses. However, as the role of pMP19 in the pathogenesis of European foulbrood remains to be elucidated, we generated pMP19 cured-M. plutonius from representative strains of the three genetically distinct groups (CC3, CC12 and CC13) and compared their virulence against Apis mellifera larvae using our in vitro infection model. Under the conditions tested, the loss of pMP19 abrogated the pathogenicity in CC3 strains, and > 94% of pMP19-cured CC3 strain-infected larvae became adult bees, suggesting that pMP19 is a virulence determinant of CC3 strains. However, introduction of mtxA on its own did not increase the virulence of pMP19-cured strains. In contrast to CC3 strains, the representative CC12 strain remained virulent even in the absence of pMP19, whereas the representative CC13 strain was avirulent even in the presence of the plasmid. Thus, pMP19 plays a role in the virulence of M. plutonius; however, its impact on the virulence varies among strains with different genetic backgrounds.     

Predictive modeling of Pseudomonas syringae virulence on bean using gradient boosted decision trees

Pseudomonas syringae is a genetically diverse bacterial species complex responsible for numerous agronomically important crop diseases. Individual P. syringae isolates are assigned pathovar designations based on their host of isolation and the associated disease symptoms, and these pathovar designations are often assumed to reflect host specificity although this assumption has rarely been rigorously tested. Here we developed a rapid seed infection assay to measure the virulence of 121 diverse P. syringae isolates on common bean (Phaseolus vulgaris). This collection includes P. syringae phylogroup 2 (PG2) bean isolates (pathovar syringae) that cause bacterial spot disease and P. syringae phylogroup 3 (PG3) bean isolates (pathovar phaseolicola) that cause the more serious halo blight disease. We found that bean isolates in general were significantly more virulent on bean than non-bean isolates and observed no significant virulence difference between the PG2 and PG3 bean isolates. However, when we compared virulence within PGs we found that PG3 bean isolates were significantly more virulent than PG3 non-bean isolates, while there was no significant difference in virulence between PG2 bean and non-bean isolates. These results indicate that PG3 strains have a higher level of host specificity than PG2 strains. We then used gradient boosting machine learning to predict each strain’s virulence on bean based on whole genome k-mers, type III secreted effector k-mers, and the presence/absence of type III effectors and phytotoxins. Our model performed best using whole genome data and was able to predict virulence with high accuracy (mean absolute error = 0.05). Finally, we functionally validated the model by predicting virulence for 16 strains and found that 15 (94%) had virulence levels within the bounds of estimated predictions. This study strengthens the hypothesis that P. syringae PG2 strains have evolved a different lifestyle than other P. syringae strains as reflected in their lower level of host specificity. It also acts as a proof-of-principle to demonstrate the power of machine learning for predicting host specific adaptation.     

Pseudo-Sanger sequencing: massively parallel production of long and near error-free reads using NGS technology

Background

Usually, next generation sequencing (NGS) technology has the property of ultra-high throughput but the read length is remarkably short compared to conventional Sanger sequencing. Paired-end NGS could computationally extend the read length but with a lot of practical inconvenience because of the inherent gaps. Now that Illumina paired-end sequencing has the ability of read both ends from 600 bp or even 800 bp DNA fragments, how to fill in the gaps between paired ends to produce accurate long reads is intriguing but challenging.

Results

We have developed a new technology, referred to as pseudo-Sanger (PS) sequencing. It tries to fill in the gaps between paired ends and could generate near error-free sequences equivalent to the conventional Sanger reads in length but with the high throughput of the Next Generation Sequencing. The major novelty of PS method lies on that the gap filling is based on local assembly of paired-end reads which have overlaps with at either end. Thus, we are able to fill in the gaps in repetitive genomic region correctly. The PS sequencing starts with short reads from NGS platforms, using a series of paired-end libraries of stepwise decreasing insert sizes. A computational method is introduced to transform these special paired-end reads into long and near error-free PS sequences, which correspond in length to those with the largest insert sizes. The PS construction has 3 advantages over untransformed reads: gap filling, error correction and heterozygote tolerance. Among the many applications of the PS construction is de novo genome assembly, which we tested in this study. Assembly of PS reads from a non-isogenic strain of Drosophila melanogaster yields an N50 contig of 190 kb, a 5 fold improvement over the existing de novo assembly methods and a 3 fold advantage over the assembly of long reads from 454 sequencing.

Conclusions

Our method generated near error-free long reads from NGS paired-end sequencing. We demonstrated that de novo assembly could benefit a lot from these Sanger-like reads. Besides, the characteristic of the long reads could be applied to such applications as structural variations detection and metagenomics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-711) contains supplementary material, which is available to authorized users.     

Pseudomonas caspiana sp. nov., a citrus pathogen in the Pseudomonas syringae phylogenetic group

In a screening by multilocus sequence analysis of Pseudomonas strains isolated from diverse origins, 4 phylogenetically closely related strains (FBF58, FBF102^T, FBF103, and FBF122) formed a well-defined cluster in the Pseudomonas syringae phylogenetic group. The strains were isolated from citrus orchards in northern Iran with disease symptoms in the leaves and stems and its pathogenicity against citrus plants was demonstrated. The whole genome of the type strain of the proposed new species (FBF102^T = CECT 9164^T = CCUG 69273^T) was sequenced and characterized. Comparative genomics with the 14 known Pseudomonas species type strains of the P. syringae phylogenetic group demonstrated that this strain belonged to a new genomic species, different from the species described thus far. Genome analysis detected genes predicted to be involved in pathogenesis, such as an atypical type 3 secretion system and two type 6 secretion systems, together with effectors and virulence factors. A polyphasic taxonomic characterization demonstrated that the 4 plant pathogenic strains represented a new species, for which the name Pseudomonas caspiana sp. nov. is proposed.     

Characterization of pPT23B, the Plasmid Involved in Syringolide Production by Pseudomonas syringae pv. tomato PT23

Avirulence gene D (avrD) in strain PT23 of Pseudomonas syringae pv. tomato (Pst) specifies the production of syringolides, which are elicitors of plant defense reactions. An 83-kb indigenous plasmid (pPT23B) that carries avrD has been mapped and characterized and a putative par region was identified. pPT23B contains a large amount of DNA that is repeated in other native plasmids in PT23. A putative mobile insertion element that occurs on plasmid pPT23A as well as on the chromosome was also identified in strain PT23. New broad-host-range expression vectors that functioned in Pst were constructed for overexpression of the cloned avrD gene and high-level production of the syringolides. Introduction of an avrD overexpression plasmid into PT23 or plasmid-cured strains led to identical syringolide peaks on HPLC with no new peaks observed. These results suggested that neither pPT23B nor other indigenous plasmids in Pst carry additional genes required for syringolide production or metabolism. Pst strains lacking pPT23B were not impaired in virulence on tomato plants.