Development of a lateral flow test for the rapid detection of <Emphasis Type="Italic">Avibacterium paragallinarum</Emphasis> in chickens suspected of having infectious coryza

Background

Infectious coryza (IC) is an acute respiratory disease of growing chickens and layers caused by Avibacterium paragallinarum. The development of tools that allow rapid pathogen detection is necessary in order to avoid disease dissemination and economic losses in poultry. An Av. paragallinarum-specific Ma-4 epitope of the TonB-dependent transporter (TBDT) was selected using bioinformatic tools in order to immunize a BalbC mouse and to produce monoclonal antibodies to be used in a lateral flow test (LFT) developed for Av. paragallinarum detection in chicken nasal mucus samples.

Results

The 1G7G8 monoclonal antibody was able to detect TBDT in Av. paragallinarum cultures (serogroups: A, B and C) by Western blot and indirect ELISA assay. Consequently, we developed a self-pairing prototype LFT. The limit of detection of the prototype LFT using Av. paragallinarum cultures was 1?×?10⁴ colony-forming units (CFU)/mL. Thirty-five nasal mucus samples from chickens suspected of having infectious coryza were evaluated for the LFT detection capacity and compared with bacterial isolation (B.I) and polymerase chain reaction (PCR). Comparative indicators such as sensitivity (Se), specificity (Sp), positive predictive value (PPV), negative predictive values (NPV) and the kappa index (K) were obtained. The values were 100.0% Se, 50% Sp, 65.4% PPV, 100% NPV, and 0.49?K and 83.9% Se, 100% Sp, 100% PPV, 44.4% NPV, and 0.54?K for the comparison of the LFT with B.I and PCR, respectively. Additionally, the LFT allowed the detection of Av. paragallinarum from coinfection cases of Av. paragallinarum with Gallibacterium anatis.

Conclusions

The results indicate that the self-pairing prototype LFT is suitable for the detection of TBDT in Av. paragallinarum cultures as well as in field samples such as nasal mucus from Av. paragallinarum-infected chickens. Therefore, this prototype LFT could be considered a rapid and promising tool to be used in farm conditions for Av. paragallinarum diagnosis.

     

Resident bacteria contribute to opportunistic infections of the respiratory tract

Opportunistic pathogens frequently cause volatile infections in hosts with compromised immune systems or a disrupted normal microbiota. The commensalism of diverse microorganisms contributes to colonization resistance, which prevents the expansion of opportunistic pathogens. Following microbiota disruption, pathogens promptly adapt to altered niches and obtain growth advantages. Nevertheless, whether and how resident bacteria modulate the growth dynamics of invasive pathogens and the eventual outcome of such infections are still unclear. Here, we utilized birds as a model animal and observed a resident bacterium exacerbating the invasion of Avibacterium paragallinarum (previously Haemophilus paragallinarum) in the respiratory tract. We first found that negligibly abundant Staphylococcus chromogenes, rather than Staphylococcus aureus, played a dominant role in Av. paragallinarum-associated infectious coryza in poultry based on epidemic investigations and in vitro analyses. Furthermore, we determined that S. chromogenes not only directly provides the necessary nutrition factor nicotinamide adenine dinucleotide (NAD⁺) but also accelerates its biosynthesis and release from host cells to promote the survival and growth of Av. paragallinarum. Last, we successfully intervened in Av. paragallinarum-associated infections in animal models using antibiotics that specifically target S. chromogenes. Our findings show that opportunistic pathogens can hijack commensal bacteria to initiate infection and expansion and suggest a new paradigm to ameliorate opportunistic infections by modulating the dynamics of resident bacteria.     

Enhanced De Novo Assembly of High Throughput Pyrosequencing Data Using Whole Genome Mapping

Despite major advances in next-generation sequencing, assembly of sequencing data, especially data from novel microorganisms or re-emerging pathogens, remains constrained by the lack of suitable reference sequences. De novo assembly is the best approach to achieve an accurate finished sequence, but multiple sequencing platforms or paired-end libraries are often required to achieve full genome coverage. In this study, we demonstrated a method to assemble complete bacterial genome sequences by integrating shotgun Roche 454 pyrosequencing with optical whole genome mapping (WGM). The whole genome restriction map (WGRM) was used as the reference to scaffold de novo assembled sequence contigs through a stepwise process. Large de novo contigs were placed in the correct order and orientation through alignment to the WGRM. De novo contigs that were not aligned to WGRM were merged into scaffolds using contig branching structure information. These extended scaffolds were then aligned to the WGRM to identify the overlaps to be eliminated and the gaps and mismatches to be resolved with unused contigs. The process was repeated until a sequence with full coverage and alignment with the whole genome map was achieved. Using this method we were able to achieved 100% WGRM coverage without a paired-end library. We assembled complete sequences for three distinct genetic components of a clinical isolate of Providencia stuartii: a bacterial chromosome, a novel bla _NDM-1 plasmid, and a novel bacteriophage, without separately purifying them to homogeneity.     

Next-generation phylogenomics using a Target Restricted Assembly Method

Next-generation sequencing technologies are revolutionizing the field of phylogenetics by making available genome scale data for a fraction of the cost of traditional targeted sequencing. One challenge will be to make use of these genomic level data without necessarily resorting to full-scale genome assembly and annotation, which is often time and labor intensive. Here we describe a technique, the Target Restricted Assembly Method (TRAM), in which the typical process of genome assembly and annotation is in essence reversed. Protein sequences of phylogenetically useful genes from a species within the group of interest are used as targets in tblastn searches of a data set from a lane of Illumina reads for a related species. Resulting blast hits are then assembled locally into contigs and these contigs are then aligned against the reference “cDNA” sequence to remove portions of the sequences that include introns. We illustrate the Target Restricted Assembly Method using genomic scale datasets for 20 species of lice (Insecta: Psocodea) to produce a test phylogenetic data set of 10 nuclear protein coding gene sequences. Given the advantages of using DNA instead of RNA, this technique is very cost effective and feasible given current technologies.     

Performance comparison of second- and third-generation sequencers using a bacterial genome with two chromosomes

Background

The availability of diverse second- and third-generation sequencing technologies enables the rapid determination of the sequences of bacterial genomes. However, identifying the sequencing technology most suitable for producing a finished genome with multiple chromosomes remains a challenge. We evaluated the abilities of the following three second-generation sequencers: Roche 454 GS Junior (GS Jr), Life Technologies Ion PGM (Ion PGM), and Illumina MiSeq (MiSeq) and a third-generation sequencer, the Pacific Biosciences RS sequencer (PacBio), by sequencing and assembling the genome of Vibrio parahaemolyticus, which consists of a 5-Mb genome comprising two circular chromosomes.

Results

We sequenced the genome of V. parahaemolyticus with GS Jr, Ion PGM, MiSeq, and PacBio and performed de novo assembly with several genome assemblers. Although GS Jr generated the longest mean read length of 418 bp among the second-generation sequencers, the maximum contig length of the best assembly from GS Jr was 165 kbp, and the number of contigs was 309. Single runs of Ion PGM and MiSeq produced data of considerably greater sequencing coverage, 279× and 1,927×, respectively. The optimized result for Ion PGM contained 61 contigs assembled from reads of 77× coverage, and the longest contig was 895 kbp in size. Those for MiSeq were 34 contigs, 58× coverage, and 733 kbp, respectively. These results suggest that higher coverage depth is unnecessary for a better assembly result. We observed that multiple rRNA coding regions were fragmented in the assemblies from the second-generation sequencers, whereas PacBio generated two exceptionally long contigs of 3,288,561 and 1,875,537 bps, each of which was from a single chromosome, with 73× coverage and mean read length 3,119 bp, allowing us to determine the absolute positions of all rRNA operons.

Conclusions

PacBio outperformed the other sequencers in terms of the length of contigs and reconstructed the greatest portion of the genome, achieving a genome assembly of “finished grade” because of its long reads. It showed the potential to assemble more complex genomes with multiple chromosomes containing more repetitive sequences.

Electronic supplementary material

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

Genome-wide analysis of the role of GlnR in Streptomyces venezuelae provides new insights into global nitrogen regulation in actinomycetes

Background

Knowledge of the origins, distribution, and inheritance of variation in the malaria parasite (Plasmodium falciparum) genome is crucial for understanding its evolution; however the 81% (A+T) genome poses challenges to high-throughput sequencing technologies. We explore the viability of the Roche 454 Genome Sequencer FLX (GS FLX) high throughput sequencing technology for both whole genome sequencing and fine-resolution characterization of genetic exchange in malaria parasites.

Results

We present a scheme to survey recombination in the haploid stage genomes of two sibling parasite clones, using whole genome pyrosequencing that includes a sliding window approach to predict recombination breakpoints. Whole genome shotgun (WGS) sequencing generated approximately 2 million reads, with an average read length of approximately 300 bp. De novo assembly using a combination of WGS and 3 kb paired end libraries resulted in contigs ≤ 34 kb. More than 8,000 of the 24,599 SNP markers identified between parents were genotyped in the progeny, resulting in a marker density of approximately 1 marker/3.3 kb and allowing for the detection of previously unrecognized crossovers (COs) and many non crossover (NCO) gene conversions throughout the genome.

Conclusions

By sequencing the 23 Mb genomes of two haploid progeny clones derived from a genetic cross at more than 30× coverage, we captured high resolution information on COs, NCOs and genetic variation within the progeny genomes. This study is the first to resequence progeny clones to examine fine structure of COs and NCOs in malaria parasites.     

Rapid Genome Mapping in Nanochannel Arrays for Highly Complete and Accurate De Novo Sequence Assembly of the Complex Aegilops tauschii Genome

Next-generation sequencing (NGS) technologies have enabled high-throughput and low-cost generation of sequence data; however, de novo genome assembly remains a great challenge, particularly for large genomes. NGS short reads are often insufficient to create large contigs that span repeat sequences and to facilitate unambiguous assembly. Plant genomes are notorious for containing high quantities of repetitive elements, which combined with huge genome sizes, makes accurate assembly of these large and complex genomes intractable thus far. Using two-color genome mapping of tiling bacterial artificial chromosomes (BAC) clones on nanochannel arrays, we completed high-confidence assembly of a 2.1-Mb, highly repetitive region in the large and complex genome of Aegilops tauschii, the D-genome donor of hexaploid wheat (Triticum aestivum). Genome mapping is based on direct visualization of sequence motifs on single DNA molecules hundreds of kilobases in length. With the genome map as a scaffold, we anchored unplaced sequence contigs, validated the initial draft assembly, and resolved instances of misassembly, some involving contigs <2 kb long, to dramatically improve the assembly from 75% to 95% complete.     

Sequencing and Analysis of the Pseudomonas fluorescens GcM5-1A Genome: A Pathogen Living in the Surface Coat of Bursaphelenchus xylophilus

It is known that several bacteria are adherent to the surface coat of pine wood nematode (Bursaphelenchus xylophilus), but their function and role in the pathogenesis of pine wilt disease remains debatable. The Pseudomonas fluorescens GcM5-1A is a bacterium isolated from the surface coat of pine wood nematodes. In previous studies, GcM5-1A was evident in connection with the pathogenicity of pine wilt disease. In this study, we report the de novo sequencing of the GcM5-1A genome. A 600-Mb collection of high-quality reads was obtained and assembled into sequence contigs spanning a 6.01-Mb length. Sequence annotation predicted 5,413 open reading frames, of which 2,988 were homologous to genes in the other four sequenced P. fluorescens isolates (SBW25, WH6, Pf0-1 and Pf-5) and 1,137 were unique to GcM5-1A. Phylogenetic studies and genome comparison revealed that GcM5-1A is more closely related to SBW25 and WH6 isolates than to Pf0-1 and Pf-5 isolates. Towards study of pathogenesis, we identified 79 candidate virulence factors in the genome of GcM5-1A, including the Alg, Fl, Waa gene families, and genes coding the major pathogenic protein fliC. In addition, genes for a complete T3SS system were identified in the genome of GcM5-1A. Such systems have proved to play a critical role in subverting and colonizing the host organisms of many gram-negative pathogenic bacteria. Although the functions of the candidate virulence factors need yet to be deciphered experimentally, the availability of this genome provides a basic platform to obtain informative clues to be addressed in future studies by the pine wilt disease research community.     

Relationship between homoeologous regulatory and structural genes in allopolyploid genome – A case study in bread wheat

Background

Whole-genome physical maps facilitate genome sequencing, sequence assembly, mapping of candidate genes, and the design of targeted genetic markers. An automated protocol was used to construct a Vitis vinifera 'Cabernet Sauvignon' physical map. The quality of the result was addressed with regard to the effect of high heterozygosity on the accuracy of contig assembly. Its usefulness for the genome-wide mapping of genes for disease resistance, which is an important trait for grapevine, was then assessed.

Results

The physical map included 29,727 BAC clones assembled into 1,770 contigs, spanning 715,684 kbp, and corresponding to 1.5-fold the genome size. Map inflation was due to high heterozygosity, which caused either the separation of allelic BACs in two different contigs, or local mis-assembly in contigs containing BACs from the two haplotypes. Genetic markers anchored 395 contigs or 255,476 kbp to chromosomes. The fully automated assembly and anchorage procedures were validated by BAC-by-BAC blast of the end sequences against the grape genome sequence, unveiling 7.3% of chimerical contigs. The distribution across the physical map of candidate genes for non-host and host resistance, and for defence signalling pathways was then studied. NBS-LRR and RLK genes for host resistance were found in 424 contigs, 133 of them (32%) were assigned to chromosomes, on which they are mostly organised in clusters. Non-host and defence signalling genes were found in 99 contigs dispersed without a discernable pattern across the genome.

Conclusion

Despite some limitations that interfere with the correct assembly of heterozygous clones into contigs, the 'Cabernet Sauvignon' physical map is a useful and reliable intermediary step between a genetic map and the genome sequence. This tool was successfully exploited for a quick mapping of complex families of genes, and it strengthened previous clues of co-localisation of major NBS-LRR clusters and disease resistance loci in grapevine.     

Projector: automatic contig mapping for gap closure purposes

Projector was designed for automatic positioning of contigs from an unfinished prokaryotic genome onto a template genome of a closely related strain or species. Projector mapped 84 contigs of Lactococcus lactis MG1363 (corresponding to 81% of the assembly nucleotides) against the genome of L.lactis IL1403. Ninety three percent of subsequent gap closure PCRs were successful. Moreover, a significant improvement in the N50 and N80 values (describing the assembly quality) was observed after the use of Projector. Because increasing numbers of bacterial genomes are being sequenced, Projector provides an efficient method to close a significant number of remaining gaps in the late stages of a genome sequencing project.     

Pyrosequencing of the northern red oak (Quercus rubra L.) chloroplast genome reveals high quality polymorphisms for population management

Given the low intraspecific chloroplast diversity detected in northern red oak (Quercus rubra L.), more powerful genetic tools are necessary to accurately characterize Q. rubra chloroplast diversity and structure. We report the sequencing, assembly, and annotation of the chloroplast genome of northern red oak via pyrosequencing and a combination of de novo and reference-guided assembly (RGA). Chloroplast DNA from 16 individuals was separated into four MID-tagged pools for a Genome Sequencer 20 quarter-run (Roche Life Sciences, Indianapolis, IN, USA). A four-step assembly method was used to generate the Q. rubra chloroplast consensus sequence: (1) reads were assembled de novo into contigs, (2) de novo contigs were aligned to a reference genome and merged to produce a consensus sequence, (3) the consensus sequence was aligned to the reference sequence and gaps between contigs were filled with reference sequence to generate a "pseudoreference", and (4) reads were mapped to the pseudoreference using RGA to generate the draft chloroplast genome. One hundred percent of the pseudoreference sequence was covered with a minimum coverage of 2× and an average coverage of 43.75×. The 161,304-bp Q. rubra chloroplast genome draft sequence contained 137 genes and one rps19 pseudogene. The sequence was compared to that of Quercus robur and Q. nigra with 951 and 186 insertion/deletion or SNP polymorphisms detected, respectively. A total of 51 intraspecific polymorphisms were detected among four northern red oak individuals. The fully sequenced and annotated Q. rubra chloroplast genome containing locations of interspecific and intraspecific polymorphisms will be essential for studying population differentiation, phylogeography, and evolutionary history of this species as well as meeting management goals such as monitoring reintroduced populations, tracking wood products, and certifying seed lots and forests.     

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.     

Decoding the oak genome: public release of sequence data,assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole‐genome shotgun (WGS) approach, without the use of costly and time‐consuming methods, such as fosmid or BAC clone‐based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS‐FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired‐end and mate‐pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired‐end reads and contaminants detected, resulting in a total of 17 910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.     

Simplifier: a web tool to eliminate redundant NGS contigs

Modern genomic sequencing technologies produce a large amount of data with reduced cost per base; however, this data consists of short reads. This reduction in the size of the reads, compared to those obtained with previous methodologies, presents new challenges, including a need for efficient algorithms for the assembly of genomes from short reads and for resolving repetitions. Additionally after abinitio assembly, curation of the hundreds or thousands of contigs generated by assemblers demands considerable time and computational resources. We developed Simplifier, a stand-alone software that selectively eliminates redundant sequences from the collection of contigs generated by ab initio assembly of genomes. Application of Simplifier to data generated by assembly of the genome of Corynebacterium pseudotuberculosis strain 258 reduced the number of contigs generated by ab initio methods from 8,004 to 5,272, a reduction of 34.14%; in addition, N50 increased from 1 kb to 1.5 kb. Processing the contigs of Escherichia coli DH10B with Simplifier reduced the mate-paired library 17.47% and the fragment library 23.91%. Simplifier removed redundant sequences from datasets produced by assemblers, thereby reducing the effort required for finalization of genome assembly in tests with data from Prokaryotic organisms.

Availability

Simplifier is available at http://www.genoma.ufpa.br/rramos/softwares/simplifier.xhtmlIt requires Sun jdk 6 or higher.