Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

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

Discovery of Phytophthora infestans Genes Expressed in Planta through Mining of cDNA Libraries

Background

Phytophthora infestans (Mont.) de Bary causes late blight of potato and tomato, and has a broad host range within the Solanaceae family. Most studies of the Phytophthora – Solanum pathosystem have focused on gene expression in the host and have not analyzed pathogen gene expression in planta.

Methodology/Principal Findings

We describe in detail an in silico approach to mine ESTs from inoculated host plants deposited in a database in order to identify particular pathogen sequences associated with disease. We identified candidate effector genes through mining of 22,795 ESTs corresponding to P. infestans cDNA libraries in compatible and incompatible interactions with hosts from the Solanaceae family.

Conclusions/Significance

We annotated genes of P. infestans expressed in planta associated with late blight using different approaches and assigned putative functions to 373 out of the 501 sequences found in the P. infestans genome draft, including putative secreted proteins, domains associated with pathogenicity and poorly characterized proteins ideal for further experimental studies. Our study provides a methodology for analyzing cDNA libraries and provides an understanding of the plant – oomycete pathosystems that is independent of the host, condition, or type of sample by identifying genes of the pathogen expressed in planta.     

Genomic characterisation of the effector complement of the potato cyst nematode Globodera pallida

Background

The potato cyst nematode Globodera pallida has biotrophic interactions with its host. The nematode induces a feeding structure – the syncytium – which it keeps alive for the duration of the life cycle and on which it depends for all nutrients required to develop to the adult stage. Interactions of G. pallida with the host are mediated by effectors, which are produced in two sets of gland cells. These effectors suppress host defences, facilitate migration and induce the formation of the syncytium.

Results

The recent completion of the G. pallida genome sequence has allowed us to identify the effector complement from this species. We identify 128 orthologues of effectors from other nematodes as well as 117 novel effector candidates. We have used in situ hybridisation to confirm gland cell expression of a subset of these effectors, demonstrating the validity of our effector identification approach. We have examined the expression profiles of all effector candidates using RNAseq; this analysis shows that the majority of effectors fall into one of three clusters of sequences showing conserved expression characteristics (invasive stage nematode only, parasitic stage only or invasive stage and adult male only). We demonstrate that further diversity in the effector pool is generated by alternative splicing. In addition, we show that effectors target a diverse range of structures in plant cells, including the peroxisome. This is the first identification of effectors from any plant pathogen that target this structure.

Conclusion

This is the first genome scale search for effectors, combined to a life-cycle expression analysis, for any plant-parasitic nematode. We show that, like other phylogenetically unrelated plant pathogens, plant parasitic nematodes deploy hundreds of effectors in order to parasitise plants, with different effectors required for different phases of the infection process.

Electronic supplementary material

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

The genome of Eimeria falciformis - reduction and specialization in a single host apicomplexan parasite

Background

The phylum Apicomplexa comprises important unicellular human parasites such as Toxoplasma and Plasmodium. Eimeria is the largest and most diverse genus of apicomplexan parasites and some species of the genus are the causative agent of coccidiosis, a disease economically devastating in poultry. We report a complete genome sequence of the mouse parasite Eimeria falciformis. We assembled and annotated the genome sequence to study host-parasite interactions in this understudied genus in a model organism host.

Results

The genome of E. falciformis is 44 Mb in size and contains 5,879 predicted protein coding genes. Comparative analysis of E. falciformis with Toxoplasma gondii shows an emergence and diversification of gene families associated with motility and invasion mainly at the level of the Coccidia. Many rhoptry kinases, among them important virulence factors in T. gondii, are absent from the E. falciformis genome. Surface antigens are divergent between Eimeria species. Comparisons with T. gondii showed differences between genes involved in metabolism, N-glycan and GPI-anchor synthesis. E. falciformis possesses a reduced set of transmembrane transporters and we suggest an altered mode of iron uptake in the genus Eimeria.

Conclusions

Reduced diversity of genes required for host-parasite interaction and transmembrane transport allow hypotheses on host adaptation and specialization of a single host parasite. The E. falciformis genome sequence sheds light on the evolution of the Coccidia and helps to identify determinants of host-parasite interaction critical for drug and vaccine development.

Electronic supplementary material

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

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

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

Unlocking the mystery of the hard-to-sequence phage genome: PaP1 methylome and bacterial immunity

Background

Whole-genome sequencing is an important method to understand the genetic information, gene function, biological characteristics and survival mechanisms of organisms. Sequencing large genomes is very simple at present. However, we encountered a hard-to-sequence genome of Pseudomonas aeruginosa phage PaP1. Shotgun sequencing method failed to complete the sequence of this genome.

Results

After persevering for 10 years and going over three generations of sequencing techniques, we successfully completed the sequence of the PaP1 genome with a length of 91,715 bp. Single-molecule real-time sequencing results revealed that this genome contains 51 N-6-methyladenines and 152 N-4-methylcytosines. Three significant modified sequence motifs were predicted, but not all of the sites found in the genome were methylated in these motifs. Further investigations revealed a novel immune mechanism of bacteria, in which host bacteria can recognise and repel modified bases containing inserts in a large scale. This mechanism could be accounted for the failure of the shotgun method in PaP1 genome sequencing. This problem was resolved using the nfi^- mutant of Escherichia coli DH5α as a host bacterium to construct a shotgun library.

Conclusions

This work provided insights into the hard-to-sequence phage PaP1 genome and discovered a new mechanism of bacterial immunity. The methylome of phage PaP1 is responsible for the failure of shotgun sequencing and for bacterial immunity mediated by enzyme Endo V activity; this methylome also provides a valuable resource for future studies on PaP1 genome replication and modification, as well as on gene regulation and host interaction.

Electronic supplementary material

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

Different dynamics of genome content shuffling among host-specificity groups of the symbiotic actinobacterium Frankia

Background

Frankia is a genus of soil actinobacteria forming nitrogen-fixing root-nodule symbiotic relationships with non-leguminous woody plant species, collectively called actinorhizals, from eight dicotyledonous families. Frankia strains are classified into four host-specificity groups (HSGs), each of which exhibits a distinct host range. Genome sizes of representative strains of Alnus, Casuarina, and Elaeagnus HSGs are highly diverged and are positively correlated with the size of their host ranges.

Results

The content and size of 12 Frankia genomes were investigated by in silico comparative genome hybridization and pulsed-field gel electrophoresis, respectively. Data were collected from four query strains of each HSG and compared with those of reference strains possessing completely sequenced genomes. The degree of difference in genome content between query and reference strains varied depending on HSG. Elaeagnus query strains were missing the greatest number (22–32%) of genes compared with the corresponding reference genome; Casuarina query strains lacked the fewest (0–4%), with Alnus query strains intermediate (14–18%). In spite of the remarkable gene loss, genome sizes of Alnus and Elaeagnus query strains were larger than would be expected based on total length of the absent genes. In contrast, Casuarina query strains had smaller genomes than expected.

Conclusions

The positive correlation between genome size and host range held true across all investigated strains, supporting the hypothesis that size and genome content differences are responsible for observed diversity in host plants and host plant biogeography among Frankia strains. In addition, our results suggest that different dynamics of shuffling of genome content have contributed to these symbiotic and biogeographic adaptations. Elaeagnus strains, and to a lesser extent Alnus strains, have gained and lost many genes to adapt to a wide range of environments and host plants. Conversely, rather than acquiring new genes, Casuarina strains have discarded genes to reduce genome size, suggesting an evolutionary orientation towards existence as specialist symbionts.

Electronic supplementary material

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

Whole genome capture of vector-borne pathogens from mixed DNA samples: a case study of Borrelia burgdorferi

Background

Rapid and accurate retrieval of whole genome sequences of human pathogens from disease vectors or animal reservoirs will enable fine-resolution studies of pathogen epidemiological and evolutionary dynamics. However, next generation sequencing technologies have not yet been fully harnessed for the study of vector-borne and zoonotic pathogens, due to the difficulty of obtaining high-quality pathogen sequence data directly from field specimens with a high ratio of host to pathogen DNA.

Results

We addressed this challenge by using custom probes for multiplexed hybrid capture to enrich for and sequence 30 Borrelia burgdorferi genomes from field samples of its arthropod vector. Hybrid capture enabled sequencing of nearly the complete genome (~99.5 %) of the Borrelia burgdorferi pathogen with 132-fold coverage, and identification of up to 12,291 single nucleotide polymorphisms per genome.

Conclusions

The proprosed culture-independent method enables efficient whole genome capture and sequencing of pathogens directly from arthropod vectors, thus making population genomic study of vector-borne and zoonotic infectious diseases economically feasible and scalable. Furthermore, given the similarities of invertebrate field specimens to other mixed DNA templates characterized by a high ratio of host to pathogen DNA, we discuss the potential applicabilty of hybrid capture for genomic study across diverse study systems.

Electronic supplementary material

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

Biogeography and Virulence of Staphylococcus aureus

Background

Staphylococcus aureus is commonly carried asymptomatically in the human anterior nares and occasionally enters the bloodstream to cause invasive disease. Much of the global diversity of S. aureus remains uncharacterised, and is not clear how disease propensity varies between strains, and between host populations.

Methodology

We compared 147 isolates recovered from five kindergartens in Chengdu, China, with 51 isolates contemporaneously recovered from cases of pediatric infection from the main hospital serving this community. The samples were characterised by MLST, the presence/absence of PVL, and antibiotic resistance profiling.

Principal Findings

Genotype frequencies within individual kindergartens differ, but the sample recovered from cases of disease shows a general enrichment of certain MLST genotypes and PVL positive isolates. Genotypes under-represented in the disease sample tend to correspond to a single sequence cluster, and this cluster is more common in China than in other parts of the world.

Conclusions/Significance

Virulence propensity likely reflects a synergy between variation in the core genome (MLST) and accessory genome (PVL). By combining evidence form biogeography and virulence we demonstrate the existence of a “native” clade in West China which has lowered virulence, possibility due to acquired host immunity.     

Analysis of the Pantoea ananatis pan-genome reveals factors underlying its ability to colonize and interact with plant,insect and vertebrate hosts

Background

Pantoea ananatis is found in a wide range of natural environments, including water, soil, as part of the epi- and endophytic flora of various plant hosts, and in the insect gut. Some strains have proven effective as biological control agents and plant-growth promoters, while other strains have been implicated in diseases of a broad range of plant hosts and humans. By analysing the pan-genome of eight sequenced P. ananatis strains isolated from different sources we identified factors potentially underlying its ability to colonize and interact with hosts in both the plant and animal Kingdoms.

Results

The pan-genome of the eight compared P. ananatis strains consisted of a core genome comprised of 3,876 protein coding sequences (CDSs) and a sizeable accessory genome consisting of 1,690 CDSs. We estimate that ~106 unique CDSs would be added to the pan-genome with each additional P. ananatis genome sequenced in the future. The accessory fraction is derived mainly from integrated prophages and codes mostly for proteins of unknown function. Comparison of the translated CDSs on the P. ananatis pan-genome with the proteins encoded on all sequenced bacterial genomes currently available revealed that P. ananatis carries a number of CDSs with orthologs restricted to bacteria associated with distinct hosts, namely plant-, animal- and insect-associated bacteria. These CDSs encode proteins with putative roles in transport and metabolism of carbohydrate and amino acid substrates, adherence to host tissues, protection against plant and animal defense mechanisms and the biosynthesis of potential pathogenicity determinants including insecticidal peptides, phytotoxins and type VI secretion system effectors.

Conclusions

P. ananatis has an ‘open’ pan-genome typical of bacterial species that colonize several different environments. The pan-genome incorporates a large number of genes encoding proteins that may enable P. ananatis to colonize, persist in and potentially cause disease symptoms in a wide range of plant and animal hosts.

Electronic supplementary material

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

Nutrient requirements differ in two Pedicularis species in the absence of a host plant: implication for driving forces in the evolution of host preference of root hemiparasitic plants

Background and Aims

Facultative root hemiparasitic plants generally have a wide host range, but in most cases show an obvious host preference. The reasons for the marked difference in growth performance of hemiparasites when attached to different hosts are not fully understood. In this study, the hypothesis was tested that hemiparasites showing a preference for different hosts have different nutrient requirements.

Methods

Two facultative root hemiparasitic Pedicularis species (P. rex and P. tricolor) with a different host dependency and preference were used to test their responses to inorganic solutes. The effects of nitrogen, phosphorus and potassium on growth of the hemiparasitic plants not attached to a host were determined, using an orthogonal design in pot cultivation under greenhouse conditions. Variables including biomass, shoot nutrient concentration, root:shoot (R:S) ratios and the number of haustoria were measured.

Key Results

As in autotrophic plants, nutrient deficiency reduced dry weight (DW) and nutrient concentrations in the root hemiparasites. Nitrogen and phosphorus significantly influenced growth of both Pedicularis species, while potassium availability influenced only shoot DW of P. rex. Nitrogen had far more effect on growth of P. rex than on P. tricolor, while phosphorus deficiency caused more marked growth depression in P. tricolor than in P. rex. Pedicularis rex grew faster than P. tricolor in a range of nutrient supplies. Different patterns of biomass allocation between the two Pedicularis species were observed. While P. rex invested more into roots (particularly fine rootlets) than P. tricolor, the number of haustoria produced by P. rex was relatively much lower than that produced by P. tricolor, which had a much smaller root system.

Conclusions

The two Pedicularis species differ in nutrient requirements and biomass allocation. Distinct interspecific traits in growth and nutrient requirements can be driving forces for the differential interactions between hemiparasites and their hosts.     

Plasmodium falciparum antigens on the surface of the gametocyte-infected erythrocyte

Background

The asexual blood stages of the human malaria parasite Plasmodium falciparum produce highly immunogenic polymorphic antigens that are expressed on the surface of the host cell. In contrast, few studies have examined the surface of the gametocyte-infected erythrocyte.

Methodology/Principal Findings

We used flow cytometry to detect antibodies recognising the surface of live cultured erythrocytes infected with gametocytes of P. falciparum strain 3D7 in the plasma of 200 Gambian children. The majority of children had been identified as carrying gametocytes after treatment for malaria, and each donated blood for mosquito-feeding experiments. None of the plasma recognised the surface of erythrocytes infected with developmental stages of gametocytes (I–IV), but 66 of 194 (34.0%) plasma contained IgG that recognised the surface of erythrocytes infected with mature (stage V) gametocytes. Thirty-four (17.0%) of 200 plasma tested recognised erythrocytes infected with trophozoites and schizonts, but there was no association with recognition of the surface of gametocyte-infected erythrocytes (odds ratio 1.08, 95% C.I. 0.434–2.57; P = 0.851). Plasma antibodies with the ability to recognise gametocyte surface antigens (GSA) were associated with the presence of antibodies that recognise the gamete antigen Pfs 230, but not Pfs48/45. Antibodies recognising GSA were associated with donors having lower gametocyte densities 4 weeks after antimalarial treatment.

Conclusions/Significance

We provide evidence that GSA are distinct from antigens detected on the surface of asexual 3D7 parasites. Our findings suggest a novel strategy for the development of transmission-blocking vaccines.     

Trans-cellular introduction of HIV-1 protein Nef induces pathogenic response in Caenorhabditis elegans

Background

Caenorhabditis elegans has emerged as a very powerful model for studying the host pathogen interactions. Despite the absence of a naturally occurring viral infection for C. elegans, the model is now being exploited experimentally to study the basic aspects of virus-host interplay. The data generated from recent studies suggests that the virus that infects mammalian cells does infect, replicate and accumulate in C. elegans.

Methodology/Principal Findings

We took advantage of the easy-to-achieve protein introduction in C. elegans and employing the methodology, we administered HIV-1 protein Nef into live worms. Nef is known to be an important protein for exacerbating HIV-1 pathogenesis in host by enhancing viral replication. The deletion of nef from the viral genome has been reported to inhibit its replication in the host, thereby leading to delayed pathogenesis. Our studies, employing Nef introduction into C. elegans, led to creation of an in-vivo model that allowed us to study, whether or not, the protein induces effect in the whole organism. We observed a marked lipodystrophy, effect on neuromuscular function, impaired fertility and reduced longevity in the worms exposed to Nef. The observed effects resemble to those observed in Nef transgenic mice and most interestingly the effects also relate to some of the pathogenic aspects exhibited by human AIDS patients.

Conclusions/Significance

Our studies underline the importance of this in vivo model for studying the interactions of Nef with host proteins, which could further be used for identifying possible inhibitors of such interactions.     

Inoculation with arbuscular mycorrhizal fungi suppresses initiation of haustoria in the root hemiparasite Pedicularis tricolor

Background and Aims

Plant parasitism and arbuscular mycorrhizal (AM) associations have many parallels and share a number of regulatory pathways. Despite a rapid increase in investigations addressing the roles of AM fungi in regulating interactions between parasitic plants and their hosts, few studies have tested the effect of AM fungi on the initiation and differentiation of haustoria, the parasite-specific structures exclusively responsible for host attachment and nutrient transfer. In this study, we tested the influence of AM fungi on haustorium formation in a root hemiparasitic plant.

Methods

Using a facultative root hemiparasitic species (Pedicularis tricolor) with the potential to form AM associations, the effects of inoculation were tested with two AM fungal species, Glomus mosseae and Glomus intraradices, on haustorium initiation in P. tricolor grown alone or with Hordeum vulgare ‘Fleet’ (barley) as the host plant. This study consisted of two greenhouse pot experiments.

Key Results

Both AM fungal species dramatically suppressed intraspecific haustorium initiation in P. tricolor at a very low colonization level. The suppression over-rode inductive effects of the parasite''s host plant on haustoria production and caused significant growth depression of P. tricolor.

Conclusions

AM fungi had strong and direct suppressive effects on haustorium formation in the root hemiparasite. The significant role of AM fungi in haustorium initiation of parasitic plants was demonstrated for the first time. This study provides new clues for the regulation of haustorium formation and a route to development of new biocontrol strategies in management of parasitic weeds.