Genomes of Two New Ammonia-Oxidizing Archaea Enriched from Deep Marine Sediments

Ammonia-oxidizing archaea (AOA) are ubiquitous and abundant and contribute significantly to the carbon and nitrogen cycles in the ocean. In this study, we assembled AOA draft genomes from two deep marine sediments from Donghae, South Korea, and Svalbard, Arctic region, by sequencing the enriched metagenomes. Three major microorganism clusters belonging to Thaumarchaeota, Epsilonproteobacteria, and Gammaproteobacteria were deduced from their 16S rRNA genes, GC contents, and oligonucleotide frequencies. Three archaeal genomes were identified, two of which were distinct and were designated Ca. “Nitrosopumilus koreensis” AR1 and “Nitrosopumilus sediminis” AR2. AR1 and AR2 exhibited average nucleotide identities of 85.2% and 79.5% to N. maritimus, respectively. The AR1 and AR2 genomes contained genes pertaining to energy metabolism and carbon fixation as conserved in other AOA, but, conversely, had fewer heme-containing proteins and more copper-containing proteins than other AOA. Most of the distinctive AR1 and AR2 genes were located in genomic islands (GIs) that were not present in other AOA genomes or in a reference water-column metagenome from the Sargasso Sea. A putative gene cluster involved in urea utilization was found in the AR2 genome, but not the AR1 genome, suggesting niche specialization in marine AOA. Co-cultured bacterial genome analysis suggested that bacterial sulfur and nitrogen metabolism could be involved in interactions with AOA. Our results provide fundamental information concerning the metabolic potential of deep marine sedimentary AOA.     

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

Background

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

Methodology/Principal Findings

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

Conclusions/Significance

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.     

Comparison of Cytochrome P450 Genes from Six Plant Genomes

Plants depend on cytochrome P450 (CYP) enzymes for nearly every aspect of their biology. In several sequenced angiosperms, CYP genes constitute up to 1% of the protein coding genes. The angiosperm sequence diversity is encapsulated by 59 CYP families, of which 52 families form a widely distributed core set. In the 20 years since the first plant P450 was sequenced, 3,387 P450 sequences have been identified and annotated in plant databases. As no new angiosperm CYP families have been discovered since 2004, it is now apparent that the sampling of CYP diversity is beginning to plateau. This review presents a comparison of 1,415 cytochrome P450 sequences from the six sequenced genomes of Vitis vinifera (grape), Carica papaya (papaya), Populus trichocarpa (poplar), Oryza sativa (rice), Arabidopsis thaliana (Arabidopsis or mouse ear’s cress) and Physcomitrella patens (moss). An evolutionary analysis is presented that tracks land plant P450 innovation over time from the most ancient and conserved sequences to the newest dicot-specific families. The earliest or oldest P450 families are devoted to the essential biochemistries of sterol and carotenoid synthesis. The next evolutionary radiation of P450 families appears to mediate crucial adaptations to a land environment. And, the newest CYP families appear to have driven the diversity of angiosperms in mediating the synthesis of pigments, odorants, flavors and order-/genus-specific secondary metabolites. Family-by-family comparisons allow the visualization of plant genome plasticity by whole genome duplications and massive gene family expansions via tandem duplications. Molecular evidence of human domestication is quite apparent in the repeated P450 gene duplications occurring in the grape genome.     

Genomic Characterization of Campylobacter jejuni Strain M1

Campylobacter jejuni strain M1 (laboratory designation 99/308) is a rarely documented case of direct transmission of C. jejuni from chicken to a person, resulting in enteritis. We have sequenced the genome of C. jejuni strain M1, and compared this to 12 other C. jejuni sequenced genomes currently publicly available. Compared to these, M1 is closest to strain 81116. Based on the 13 genome sequences, we have identified the C. jejuni pan-genome, as well as the core genome, the auxiliary genes, and genes unique between strains M1 and 81116. The pan-genome contains 2,427 gene families, whilst the core genome comprised 1,295 gene families, or about two-thirds of the gene content of the average of the sequenced C. jejuni genomes. Various comparison and visualization tools were applied to the 13 C. jejuni genome sequences, including a species pan- and core genome plot, a BLAST Matrix and a BLAST Atlas. Trees based on 16S rRNA sequences and on the total gene families in each genome are presented. The findings are discussed in the background of the proven virulence potential of M1.     

Novel Virus Discovery and Genome Reconstruction from Field RNA Samples Reveals Highly Divergent Viruses in Dipteran Hosts

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected.     

The Complete Sequence of the Mitochondrial Genome of Butomus umbellatus – A Member of an Early Branching Lineage of Monocotyledons

In order to study the evolution of mitochondrial genomes in the early branching lineages of the monocotyledons, i.e., the Acorales and Alismatales, we are sequencing complete genomes from a suite of key taxa. As a starting point the present paper describes the mitochondrial genome of Butomus umbellatus (Butomaceae) based on next-generation sequencing data. The genome was assembled into a circular molecule, 450,826 bp in length. Coding sequences cover only 8.2% of the genome and include 28 protein coding genes, four rRNA genes, and 12 tRNA genes. Some of the tRNA genes and a 16S rRNA gene are transferred from the plastid genome. However, the total amount of recognized plastid sequences in the mitochondrial genome is only 1.5% and the amount of DNA transferred from the nucleus is also low. RNA editing is abundant and a total of 557 edited sites are predicted in the protein coding genes. Compared to the 40 angiosperm mitochondrial genomes sequenced to date, the GC content of the Butomus genome is uniquely high (49.1%). The overall similarity between the mitochondrial genomes of Butomus and Spirodela (Araceae), the closest relative yet sequenced, is low (less than 20%), and the two genomes differ in size by a factor 2. Gene order is also largely unconserved. However, based on its phylogenetic position within the core alismatids Butomus will serve as a good reference point for subsequent studies in the early branching lineages of the monocotyledons.     

A Mesophilic,Autotrophic, Ammonia-Oxidizing Archaeon of Thaumarchaeal Group I.1a Cultivated from a Deep Oligotrophic Soil Horizon

Soil nitrification plays an important role in the reduction of soil fertility and in nitrate enrichment of groundwater. Various ammonia-oxidizing archaea (AOA) are considered to be members of the pool of ammonia-oxidizing microorganisms in soil. This study reports the discovery of a chemolithoautotrophic ammonia oxidizer that belongs to a distinct clade of nonmarine thaumarchaeal group I.1a, which is widespread in terrestrial environments. The archaeal strain MY2 was cultivated from a deep oligotrophic soil horizon. The similarity of the 16S rRNA gene sequence of strain MY2 to those of other cultivated group I.1a thaumarchaeota members, i.e., Nitrosopumilus maritimus and “Candidatus Nitrosoarchaeum koreensis,” is 92.9% for both species. Extensive growth assays showed that strain MY2 is chemolithoautotrophic, mesophilic (optimum temperature, 30°C), and neutrophilic (optimum pH, 7 to 7.5). The accumulation of nitrite above 1 mM inhibited ammonia oxidation, while ammonia oxidation itself was not inhibited in the presence of up to 5 mM ammonia. The genome size of strain MY2 was 1.76 Mb, similar to those of N. maritimus and “Ca. Nitrosoarchaeum koreensis,” and the repertoire of genes required for ammonia oxidation and carbon fixation in thaumarchaeal group I.1a was conserved. A high level of representation of conserved orthologous genes for signal transduction and motility in the noncore genome might be implicated in niche adaptation by strain MY2. On the basis of phenotypic, phylogenetic, and genomic characteristics, we propose the name “Candidatus Nitrosotenuis chungbukensis” for the ammonia-oxidizing archaeal strain MY2.     

An in vitro strategy for the selective isolation of anomalous DNA from prokaryotic genomes

In sequenced genomes of prokaryotes, anomalous DNA (aDNA) can be recognized, among others, by atypical clustering of dinucleotides. We hypothesized that atypical clustering of hexameric endonuclease recognition sites in aDNA allows the specific isolation of anomalous sequences in vitro. Clustering of endonuclease recognition sites in aDNA regions of eight published prokaryotic genome sequences was demonstrated. In silico digestion of the Neisseria meningitidis MC58 genome, using four selected endonucleases, revealed that out of 27 of the small fragments predicted (<5 kb), 21 were located in known genomic islands. Of the 24 calculated fragments (>300 bp and <5 kb), 22 met our criteria for aDNA, i.e. a high dinucleotide dissimilarity and/or aberrant GC content. The four enzymes also allowed the identification of aDNA fragments from the related Z2491 strain. Similarly, the sequenced genomes of three strains of Escherichia coli assessed by in silico digestion using XbaI yielded strain-specific sets of fragments of anomalous composition. In vitro applicability of the method was demonstrated by using adaptor-linked PCR, yielding the predicted fragments from the N.meningitidis MC58 genome. In conclusion, this strategy allows the selective isolation of aDNA from prokaryotic genomes by a simple restriction digest–amplification–cloning–sequencing scheme.     

Development and annotation of shotgun sequence libraries from New World monkeys

The draft genome sequences of several primates are available, providing insights into evolutionary and anthropological research. However, genomic resources from New World monkeys are conspicuously lacking. To date, the genomes of only two platyrrhine species, the common marmoset and the Bolivian squirrel monkey, have been fully sequenced. This is especially limiting for comparative genomics research, considering that New World monkeys are the most speciose primate group, and platyrrhine genetic diversity is comparable to that of the catarrhines (i.e. apes and Old World monkeys). Here, we present the generation and annotation of numerous sequence reads from the genomes of Spider monkey (Ateles belzebuth), Owl monkey (Aotus lemurinus) and Uakari (Cacajao calvus), representing the three platyrrhine families, Atelidae, Cebidae and Pitheciidae, respectively. These sequencing reads were developed from gDNA shotgun libraries containing over 3000 individual sequences with an average length of 726 bps. Of these sequences, 1220 contain <20% repeats, and thus are potentially highly useful phylogenetic markers for other platyrrhine species. Among them, a large number of sequencing reads were found to match unique regions within the human (2462 sequences) and the marmoset (2829 sequences) genomes. In particular, the majority of these sequencing reads are from putatively neutrally evolving intergenic regions. Thus, they are likely to be highly informative for inferring neutral evolutionary patterns and genomic evolution for other New World monkeys.     

Isotopic signatures of N2O produced by ammonia-oxidizing archaea from soils

N₂O gas is involved in global warming and ozone depletion. The major sources of N₂O are soil microbial processes. Anthropogenic inputs into the nitrogen cycle have exacerbated these microbial processes, including nitrification. Ammonia-oxidizing archaea (AOA) are major members of the pool of soil ammonia-oxidizing microorganisms. This study investigated the isotopic signatures of N₂O produced by soil AOA and associated N₂O production processes. All five AOA strains (I.1a, I.1a-associated and I.1b clades of Thaumarchaeota) from soil produced N₂O and their yields were comparable to those of ammonia-oxidizing bacteria (AOB). The levels of site preference (SP), δ¹⁵N^bulk and δ¹⁸O -N₂O of soil AOA strains were 13–30%, −13 to −35% and 22–36%, respectively, and strains MY1–3 and other soil AOA strains had distinct isotopic signatures. A ¹⁵N-NH₄⁺-labeling experiment indicated that N₂O originated from two different production pathways (that is, ammonia oxidation and nitrifier denitrification), which suggests that the isotopic signatures of N₂O from AOA may be attributable to the relative contributions of these two processes. The highest N₂O production yield and lowest site preference of acidophilic strain CS may be related to enhanced nitrifier denitrification for detoxifying nitrite. Previously, it was not possible to detect N₂O from soil AOA because of similarities between its isotopic signatures and those from AOB. Given the predominance of AOA over AOB in most soils, a significant proportion of the total N₂O emissions from soil nitrification may be attributable to AOA.     

Enrichment and Genome Sequence of the Group I.1a Ammonia-Oxidizing Archaeon “Ca. Nitrosotenuis uzonensis” Representing a Clade Globally Distributed in Thermal Habitats

The discovery of ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota and the high abundance of archaeal ammonia monooxygenase subunit A encoding gene sequences in many environments have extended our perception of nitrifying microbial communities. Moreover, AOA are the only aerobic ammonia oxidizers known to be active in geothermal environments. Molecular data indicate that in many globally distributed terrestrial high-temperature habits a thaumarchaeotal lineage within the Nitrosopumilus cluster (also called “marine” group I.1a) thrives, but these microbes have neither been isolated from these systems nor functionally characterized in situ yet. In this study, we report on the enrichment and genomic characterization of a representative of this lineage from a thermal spring in Kamchatka. This thaumarchaeote, provisionally classified as “Candidatus Nitrosotenuis uzonensis”, is a moderately thermophilic, non-halophilic, chemolithoautotrophic ammonia oxidizer. The nearly complete genome sequence (assembled into a single scaffold) of this AOA confirmed the presence of the typical thaumarchaeotal pathways for ammonia oxidation and carbon fixation, and indicated its ability to produce coenzyme F₄₂₀ and to chemotactically react to its environment. Interestingly, like members of the genus Nitrosoarchaeum, “Candidatus N. uzonensis” also possesses a putative artubulin-encoding gene. Genome comparisons to related AOA with available genome sequences confirmed that the newly cultured AOA has an average nucleotide identity far below the species threshold and revealed a substantial degree of genomic plasticity with unique genomic regions in “Ca. N. uzonensis”, which potentially include genetic determinants of ecological niche differentiation.     

Comparative genomic analysis of two Burkholderia glumae strains from different geographic origins reveals a high degree of plasticity in genome structure associated with genomic islands

Burkholderia glumae is the major causal agent of bacterial panicle blight of rice, a growing disease problem in global rice production. To better understand its genome-scale characteristics, the genome of the highly virulent B. glumae strain 336gr-1 isolated from Louisiana, USA was sequenced using the Illumina Genome Analyser II system. De novo assembled 336gr-1 contigs were aligned and compared with the previously sequenced genome of B. glumae strain BGR1, which was isolated from an infected rice plant in South Korea. Comparative analysis of the whole genomes of B. glumae 336gr-1 and B. glumae BGR1 revealed numerous unique genomic regions present only in one of the two strains. These unique regions contained accessory genes including mobile elements and phage-related genes, and some of the unique regions in B. glumae BGR1 corresponded to predicted genomic islands. In contrast, little variation was observed in known and potential virulence genes between the two genomes. The considerable amount of plasticity largely based on accessory genes and genome islands observed from the comparison of the genomes of these two strains of B. glumae may explain the versatility of this bacterial species in various environmental conditions and geographic locations.     

Generation of whole genome sequences of new Cryptosporidium hominis and Cryptosporidium parvum isolates directly from stool samples

Background

Whole genome sequencing (WGS) of Cryptosporidium spp. has previously relied on propagation of the parasite in animals to generate enough oocysts from which to extract DNA of sufficient quantity and purity for analysis. We have developed and validated a method for preparation of genomic Cryptosporidium DNA suitable for WGS directly from human stool samples and used it to generate 10 high-quality whole Cryptosporidium genome assemblies. Our method uses a combination of salt flotation, immunomagnetic separation (IMS), and surface sterilisation of oocysts prior to DNA extraction, with subsequent use of the transposome-based Nextera XT kit to generate libraries for sequencing on Illumina platforms. IMS was found to be superior to caesium chloride density centrifugation for purification of oocysts from small volume stool samples and for reducing levels of contaminant DNA.

Results

The IMS-based method was used initially to sequence whole genomes of Cryptosporidium hominis gp60 subtype IbA10G2 and Cryptosporidium parvum gp60 subtype IIaA19G1R2 from small amounts of stool left over from diagnostic testing of clinical cases of cryptosporidiosis. The C. parvum isolate was sequenced to a mean depth of 51.8X with reads covering 100 % of the bases of the C. parvum Iowa II reference genome (Bioproject PRJNA 15586), while the C. hominis isolate was sequenced to a mean depth of 34.7X with reads covering 98 % of the bases of the C. hominis TU502 v1 reference genome (Bioproject PRJNA 15585).The method was then applied to a further 17 stools, successfully generating another eight new whole genome sequences, of which two were C. hominis (gp60 subtypes IbA10G2 and IaA14R3) and six C. parvum (gp60 subtypes IIaA15G2R1 from three samples, and one each of IIaA17G1R1, IIaA18G2R1, and IIdA22G1), demonstrating the utility of this method to sequence Cryptosporidium genomes directly from clinical samples. This development is especially important as it reduces the requirement to propagate Cryptosporidium oocysts in animal models prior to genome sequencing.

Conclusion

This represents the first report of high-quality whole genome sequencing of Cryptosporidium isolates prepared directly from human stool samples.     

Cultivation of a highly enriched ammonia-oxidizing archaeon of thaumarchaeotal group I.1b from an agricultural soil

Nitrification of excess ammonia in soil causes eutrophication of water resources and emission of atmospheric N₂O gas. The first step of nitrification, ammonia oxidation, is mediated by Archaea as well as Bacteria. The physiological reactions mediated by ammonia‐oxidizing archaea (AOA) and their contribution to soil nitrification are still unclear. Results of non‐culture‐based studies have shown the thaumarchaeotal group I.1b lineage of AOA to be dominant over both AOA of group I.1a and ammonia‐oxidizing bacteria in various soils. We obtained from an agricultural soil a highly enriched ammonia‐oxidizing culture dominated by a single archaeal population [c. 90% of total cells, as determined microscopically (by fluorescence in situ hybridization) and by quantitative PCR of its 16S rRNA gene]. The archaeon (termed ‘strain JG1’) fell within thaumarchaeotal group I.1b and was related to the moderately thermophilic archaeon, Candidatus Nitrososphaera gargensis, and the mesophilic archaeon, Ca. Nitrososphaera viennensis with 97.0% and 99.1% 16S rRNA gene sequence similarity respectively. Strain JG1 was neutrophilic (growth range pH 6.0–8.0) and mesophilic (growth range temperature 25–40°C). The optimum temperature of strain JG1 (35–40°C) is > 10°C higher than that of ammonia‐oxidizing bacteria (AOB). Membrane analysis showed that strain JG1 contained a glycerol dialkyl glycerol tetraether, GDGT‐4, and its regioisomer as major core lipids; this crenarchaeol regioisomer was previously detected in similar abundance in the thermophile, Ca. N. gargensis and has been frequently observed in tropical soils. Substrate uptake assays showed that the affinity of strain JG1 for ammonia and oxygen was much higher than those of AOB. These traits may give a competitive advantage to AOA related to strain JG1 in oligotrophic environments. ¹³C‐bicarbonate incorporation into archaeal lipids of strain JG1 established its ability to grow autotrophically. Strain JG1 produced a significant amount of N₂O gas – implicating AOA as a possible source of N₂O emission from soils. Sequences of archaeal amoA and 16S rRNA genes closely related to those of strain JG1 have been retrieved from various terrestrial environments in which lineage of strain JG1 is likely engaged in autotrophic nitrification.     

High-throughput sequencing of three Lemnoideae (duckweeds) chloroplast genomes from total DNA

Background

Chloroplast genomes provide a wealth of information for evolutionary and population genetic studies. Chloroplasts play a particularly important role in the adaption for aquatic plants because they float on water and their major surface is exposed continuously to sunlight. The subfamily of Lemnoideae represents such a collection of aquatic species that because of photosynthesis represents one of the fastest growing plant species on earth.

Methods

We sequenced the chloroplast genomes from three different genera of Lemnoideae, Spirodela polyrhiza, Wolffiella lingulata and Wolffia australiana by high-throughput DNA sequencing of genomic DNA using the SOLiD platform. Unfractionated total DNA contains high copies of plastid DNA so that sequences from the nucleus and mitochondria can easily be filtered computationally. Remaining sequence reads were assembled into contiguous sequences (contigs) using SOLiD software tools. Contigs were mapped to a reference genome of Lemna minor and gaps, selected by PCR, were sequenced on the ABI3730xl platform.

Conclusions

This combinatorial approach yielded whole genomic contiguous sequences in a cost-effective manner. Over 1,000-time coverage of chloroplast from total DNA were reached by the SOLiD platform in a single spot on a quadrant slide without purification. Comparative analysis indicated that the chloroplast genome was conserved in gene number and organization with respect to the reference genome of L. minor. However, higher nucleotide substitution, abundant deletions and insertions occurred in non-coding regions of these genomes, indicating a greater genomic dynamics than expected from the comparison of other related species in the Pooideae. Noticeably, there was no transition bias over transversion in Lemnoideae. The data should have immediate applications in evolutionary biology and plant taxonomy with increased resolution and statistical power.     

A Novel Construction of Genome Space with Biological Geometry

A genome space is a moduli space of genomes. In this space, each point corresponds to a genome. The natural distance between two points in the genome space reflects the biological distance between these two genomes. Currently, there is no method to represent genomes by a point in a space without losing biological information. Here, we propose a new graphical representation for DNA sequences. The breakthrough of the subject is that we can construct the moment vectors from DNA sequences using this new graphical method and prove that the correspondence between moment vectors and DNA sequences is one-to-one. Using these moment vectors, we have constructed a novel genome space as a subspace in R^N. It allows us to show that the SARS-CoV is most closely related to a coronavirus from the palm civet not from a bird as initially suspected, and the newly discovered human coronavirus HCoV-HKU1 is more closely related to SARS than to any other known member of group 2 coronavirus. Furthermore, we reconstructed the phylogenetic tree for 34 lentiviruses (including human immunodeficiency virus) based on their whole genome sequences. Our genome space will provide a new powerful tool for analyzing the classification of genomes and their phylogenetic relationships.     

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.     

2× genomes - depth does matter

Background  

Given the availability of full genome sequences, mapping gene gains, duplications, and losses during evolution should theoretically be straightforward. However, this endeavor suffers from overemphasis on detecting conserved genome features, which in turn has led to sequencing multiple eutherian genomes with low coverage rather than fewer genomes with high-coverage and more even distribution in the phylogeny. Although limitations associated with analysis of low coverage genomes are recognized, they have not been quantified.