Draft Genome Sequence of the Human Pathogen Halomonas stevensii S18214T

Halomonas stevensii is a Gram-negative, moderately halophilic bacterium causing environmental contamination and infections in a dialysis center. Here we present the 3.7-Mb draft genome sequence of the type strain (S18214(T)) of H. stevensii, which will give insight into the pathogenic potential of H. stevensii.     

Resequencing the Mycobacterium avium subsp. paratuberculosis K10 Genome: Improved Annotation and Revised Genome Sequence

We report the resequencing and revised annotation of the Mycobacterium avium subsp. paratuberculosis K10 genome. A total of 90 single-nucleotide errors and a 51-bp indel in the original K10 genome were corrected, and the whole genome annotation was revised. Correction of these sequencing errors resulted in 28 frameshift alterations. The amended genome sequence is accessible via the supplemental section of study SRR060191 in the NCBI Sequence Read Archive and will serve as a valuable reference genome for future studies.The American bovine isolate K10 remains the only Mycobacterium avium subsp. paratuberculosis genome to be fully sequenced and published to date (1). Although this 4.8-Mbp genome likely contains some assembly errors (3), it has provided, and will continue to provide, an invaluable resource for Mycobacterium research. The assembly errors were identified through optical mapping of related M. avium subsp. paratuberculosis strain ATCC 19698, which revealed a 648-kb inversion around the origin of replication and two additional copies of the insertion sequences IS1311 and IS_MAP03 (3). These findings were subsequently validated via PCR, Southern blotting, and (Sanger) sequence analysis in ATCC 19698 and were also confirmed to be present in K10 (3). We designate this interim corrected genome M. avium subsp. paratuberculosis K10′. To further improve this resource, we undertook a resequencing project of the original M. avium subsp. paratuberculosis K10 genome.Whole-genome sequencing was performed on the Illumina GAIIx platform using one flow cell lane with 36-cycle paired-end chemistry. Reads were variably trimmed at the 3′ end based on the Illumina Read Segment Quality Indicator (Illumina manual), and read pairs containing ambiguous bases were removed. Read mapping onto the K10′ genome sequence was performed using SHRiMP (ver. 1.3.2) (2), and single-nucleotide polymorphisms and indels (deletion and insertion polymorphisms [DIPs]) were called using Nesoni (ver. 0.29; Monash University Victorian Bioinformatics Consortium) with default parameters. Read mapping determined that the data set comprised an average sequence coverage of 72.6 across the K10′ genome. This high sequence coverage allowed differences between K10\K10′ and the resequenced version of the genome, designated K10", to be identified with high confidence.Ninety single-nucleotide differences and one 51-bp indel were identified in the K10" genome. As confirmation that these differences are likely to represent errors in the original genome sequence, we have also detected these polymorphisms in two additional bovine M. avium subsp. paratuberculosis genomes recently sequenced and assembled within our laboratory (data not shown). Seven of the 90 differences and the 51-bp indel were subjected to PCR and Sanger sequencing for verification. All of the polymorphisms were confirmed to be present in K10" compared to the original genome sequence.Thirty-six single-nucleotide deletions and four nucleotide insertions were identified in K10" compared to the reference. These DIPs resulted in 27 frameshift mutations of protein coding loci. As a consequence of these frameshifts, one complete coding sequence (CDS) feature was removed (MAPK_3751), one novel CDS was created (MAPK_2081b), and one pseudogene was repaired (MAPK_4158-4159). In almost all of the other cases, the frameshifts resulted in proteins which more closely resembled their orthologs in M. avium subsp. hominissuis and M. intracellulare. Other frameshifts of biological interest include the truncation of a PPE family protein (MAPK_1173) and the extension of an MCE (mammalian cell entry) family protein (MAPK_4086). Compared to the reference, K10" also had a 51-bp indel within a possible MCE family protein (MAPK_1575). This indel consisted of an 11-bp deletion (bases 2436510 to 2436520 in the original K10 sequence) and an insertion of 51 bp. The resulting protein sequence now more closely resembles orthologs of the MCE family in other Mycobacterium spp. In conclusion, the fact that so many of the amended bases have resulted in revised coding regions indicates the underlying importance of this exercise.     

Virulence and Draft Genome Sequence Overview of Multiple Strains of the Swine Pathogen Haemophilus parasuis

Haemophilus parasuis is the cause of Glässer''s disease in swine, which is characterized by systemic infection resulting in polyserositis, meningitis, and arthritis. Investigation of this animal disease is complicated by the enormous differences in the severity of disease caused by H. parasuis strains, ranging from lethal systemic disease to subclinical carriage. To identify differences in genotype that could account for virulence phenotypes, we established the virulence of, and performed whole genome sequence analysis on, 11 H. parasuis strains. Virulence was assessed by evaluating morbidity and mortality following intranasal challenge of Caesarean-derived, colostrum-deprived (CDCD) pigs. Genomic DNA from strains Nagasaki (serotype 5), 12939 (serotype 1), SW140 (serotype 2), 29755 (serotype 5), MN-H (serotype 13), 84-15995 (serotype 15), SW114 (serotype 3), H465 (serotype 11), D74 (serotype 9), and 174 (serotype 7) was used to generate Illumina paired-end libraries for genomic sequencing and de novo assembly. H. parasuis strains Nagasaki, 12939, SH0165 (serotype 5), SW140, 29755, and MN-H exhibited a high level of virulence. Despite minor differences in expression of disease among these groups, all pigs challenged with these strains developed clinical signs consistent with Glässer''s disease between 1–7 days post-challenge. H. parasuis strains 84-15995 and SW114 were moderately virulent, in that approximately half of the pigs infected with each developed Glässer''s disease. H. parasuis strains H465, D74, and 174 were minimally virulent or avirulent in the CDCD pig model. Comparative genomic analysis among strains identified several noteworthy differences in coding regions. These coding regions include predicted outer membrane, metabolism, and pilin or adhesin related genes, some of which likely contributed to the differences in virulence and systemic disease observed following challenge. These data will be useful for identifying H. parasuis virulence factors and vaccine targets.     

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.     

Genome Sequence of the Emerging Pathogen Helicobacter canadensis

We determined the genome sequence of the type strain of Helicobacter canadensis, an emerging human pathogen with diverse animal reservoirs. Potential virulence determinants carried by the genome include systems for N-linked glycosylation and capsular export. A protein-based phylogenetic analysis places H. canadensis close to Wolinella succinogenes.Helicobacter canadensis is an emerging pathogen that has been isolated from four Canadian patients with diarrhea and an Australian patient with bacteremia (3, 10). Wild geese have been identified as a reservoir in Europe (12, 13), while in China, the organism has been isolated from the feces of wild rodents (5). H. canadensis has also been isolated from laboratory rabbits (11) and from Guinea fowl in France (8). Atypical isolates have been obtained from swine feces from The Netherlands and Denmark (6). To shed light on the virulence and colonization factors of H. canadensis and to reassess its phylogenetic status, we performed whole-genome sequencing of the type strain, H. canadensis strain MIT 98-5491/NCTC 13241.Single- and paired-end whole-genome shotgun sequencing was performed using 454 pyrosequencing technology, supplemented by Solexa sequencing. An initial assembly of the 454 single-end and paired-end data was created using a Newbler assembler (Roche), generating four scaffolds with an average size of 402 kb. The genome sequence was completed using a combination of BLASTX searches and analysis of the Solexa data, followed by confirmatory PCRs, PCR-assisted contig extension (1), and combinatorial PCR, with conventional and long-PCR protocols. Annotation was performed using GenDB (7).The genome of H. canadensis NCTC 13241 consists of a single circular chromosome 1,623,845 bp in length, with a G+C content of 34%. There are 1,535 protein-coding sequences (CDSs), 40 transfer RNAs, and three rRNA loci.Although phylogenetic analyses based on 16S rRNA gene sequences place H. canadensis in the genus Helicobacter (3), analysis of 23S rRNA sequences supports a clade containing H. canadensis and Wolinella succinogenes (2). Our own phylogenetic analysis using a concatenation of 482 conserved CDSs from the H. canadensis genome and related genomes provides strong support for a W. succinogenes/H. canadensis clade (data not shown), suggesting that the current taxonomy should be reevaluated.One hundred seventeen CDSs from H. canadensis have no detectable ortholog in eight other genome-sequenced epsilonproteobacteria. We found only one contiguous region of difference longer than 10 genes (ROD1; HCAN_0630-HCAN_0663). Most CDSs in ROD1 are of unknown function, although, curiously, the region carries three versions of asparagine synthetase and contains six homopolymeric tracts. A second region of difference (ROD2; HCAN_0479-HCAN_0496) shows homology to HHGI1, a pathogenicity island from Helicobacter hepaticus strain ATCC 51449 that contributes to virulence in H. hepaticus (4).The H. canadensis NCTC 13241 genome contains 29 potentially phase-variable genes with homopolymeric tracts, including several genes that, by homology, might be implicated in virulence, e.g., an immunoglobulin A protease (HCAN_0234) and two homologues of the vacuolating cytotoxin from Helicobacter pylori (HCAN_0457 and HCAN_714). The genome contains a capsular polysaccharide export locus similar to that in Campylobacter jejuni (9). Also, like C. jejuni, H. canadensis possesses genes encoding an N-linked glycosylation pathway, including two copies of PglB (HCAN_0729 and HCAN_0930).     

Revised Sequence and Annotation of the Rhodobacter sphaeroides 2.4.1 Genome

The DNA sequences of chromosomes I and II of Rhodobacter sphaeroides strain 2.4.1 have been revised, and the annotation of the entire genomic sequence, including both chromosomes and the five plasmids, has been updated. Errors in the originally published sequence have been corrected, and ∼11% of the coding regions in the original sequence have been affected by the revised annotation.     

The Unique Structure of the Apicoplast Genome of the Rodent Malaria Parasite Plasmodium chabaudi chabaudi

The apicoplast, a non-photosynthetic plastid of apicomplexan species, has an extremely reduced but highly conserved genome. Here, the apicoplast genome of the rodent malaria parasite Plasmodium chabaudi chabaudi (Pcc) isolate CB was characterized. Although the set of genes in the genome is identical, the copy number of some tRNA genes differs between Pcc and other Plasmodium species because the Pcc DNA has only one rRNA/tRNA gene cluster, which is normally duplicated in other species. The location of the duplicated trnR(ACG) and trnM implies that one of the duplicated clusters in the ancestral molecule has been lost due to an intramolecular recombination event. The Pcc DNA occurs in two isoforms with an internal inversion between them. The presence of a unique variant in the duplicated trnT gene suggests that the two isoforms are interconvertible. This is the first report of the complete nucleotide sequence of a Plasmodium apicoplast DNA.     

Genome Sequence of the Pea Aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.     

De novo Genome Assembly of the Fungal Plant Pathogen Pyrenophora semeniperda

Pyrenophora semeniperda (anamorph Drechslera campulata) is a necrotrophic fungal seed pathogen that has a wide host range within the Poaceae. One of its hosts is cheatgrass (Bromus tectorum), a species exotic to the United States that has invaded natural ecosystems of the Intermountain West. As a natural pathogen of cheatgrass, P. semeniperda has potential as a biocontrol agent due to its effectiveness at killing seeds within the seed bank; however, few genetic resources exist for the fungus. Here, the genome of P. semeniperda isolate assembled from sequence reads of 454 pyrosequencing is presented. The total assembly is 32.5 Mb and includes 11,453 gene models encoding putative proteins larger than 24 amino acids. The models represent a variety of putative genes that are involved in pathogenic pathways typically found in necrotrophic fungi. In addition, extensive rearrangements, including inter- and intrachromosomal rearrangements, were found when the P. semeniperda genome was compared to P. tritici-repentis, a related fungal species.     

Genome Sequence of the Rice Pathogen Pseudomonas fuscovaginae CB98818

Pseudomonas fuscovaginae is a phytopathogenic bacterium causing bacterial sheath brown rot of cereal crops. Here, we present the draft genome sequence of P. fuscovaginae CB98818, originally isolated from a diseased rice plant in China. The draft genome will aid in epidemiological studies, comparative genomics, and quarantine of this broad-host-range pathogen.     

The Complete Mitochondrial Genome of the Foodborne Parasitic Pathogen Cyclospora cayetanensis

Cyclospora cayetanensis is a human-specific coccidian parasite responsible for several food and water-related outbreaks around the world, including the most recent ones involving over 900 persons in 2013 and 2014 outbreaks in the USA. Multicopy organellar DNA such as mitochondrion genomes have been particularly informative for detection and genetic traceback analysis in other parasites. We sequenced the C. cayetanensis genomic DNA obtained from stool samples from patients infected with Cyclospora in Nepal using the Illumina MiSeq platform. By bioinformatically filtering out the metagenomic reads of non-coccidian origin sequences and concentrating the reads by targeted alignment, we were able to obtain contigs containing Eimeria-like mitochondrial, apicoplastic and some chromosomal genomic fragments. A mitochondrial genomic sequence was assembled and confirmed by cloning and sequencing targeted PCR products amplified from Cyclospora DNA using primers based on our draft assembly sequence. The results show that the C. cayetanensis mitochondrion genome is 6274 bp in length, with 33% GC content, and likely exists in concatemeric arrays as in Eimeria mitochondrial genomes. Phylogenetic analysis of the C. cayetanensis mitochondrial genome places this organism in a tight cluster with Eimeria species. The mitochondrial genome of C. cayetanensis contains three protein coding genes, cytochrome (cytb), cytochrome C oxidase subunit 1 (cox1), and cytochrome C oxidase subunit 3 (cox3), in addition to 14 large subunit (LSU) and nine small subunit (SSU) fragmented rRNA genes.     

Genome Sequence of the Rice Pathogen Dickeya zeae Strain ZJU1202

Dickeya zeae is a phytopathogenic bacterium causing soft rot diseases in a wide range of economically important crops. Here we present the draft genome sequence of strain ZJU1202, which is the causal agent of rice foot rot in China. The draft genome will contribute to epidemiological and comparative genomic studies and the quarantine of this devastating phytopathogen.     

Genome Sequence of Human Adenovirus Type 55, a Re-Emergent Acute Respiratory Disease Pathogen in China

Human adenovirus type 55 (HAdV-B55) is an acute respiratory disease (ARD) pathogen first completely characterized in China (2006). This is a unique Trojan horse microbe with the virus neutralization attribute of a renal pathogen and the cell tropism and clinical attributes of a respiratory pathogen, bypassing herd immunity. It appeared to be an uncommon pathogen, with earlier putative, sporadic occurrences in Spain (1969) and Turkey (2004); these isolates were incompletely characterized using only two epitopes. Reported here is the genome of a second recent isolate (China, 2011), indicating that it may occur more frequently. The availability of this HAdV-B55 genome provides a foundation for studying adenovirus molecular evolution, the dynamics of epidemics, and patterns of pathogen emergence and re-emergence. These data facilitate studies to predict genome recombination between adenoviruses, as well as sequence divergence rates and hotspots, all of which are important for vaccine development and because HAdVs are used for epitope and/or gene delivery vectors.     

Multi-Locus Sequence Typing of a Geographically and Temporally Diverse Sample of the Highly Clonal Human Pathogen Bartonella quintana

Bartonella quintana is a re-emerging pathogen and the causative agent of a variety of disease manifestations in humans including trench fever. Various typing methods have been developed for B. quintana, but these tend to be limited by poor resolution and, in the case of gel-based methods, a lack of portability. Multilocus sequence typing (MLST) has been used to study the molecular epidemiology of a large number of pathogens, including B. henselae, a close relative of B. quintana. We developed a MLST scheme for B. quintana based on the 7 MLST loci employed for B. henselae with two additional loci to cover underrepresented regions of the B. quintana chromosome. A total of 16 B. quintana isolates spanning over 60 years and three continents were characterized. Allelic variation was detected in five of the nine loci. Although only 8/4270 (0.002%) of the nucleotide sites examined were variable over all loci, these polymorphisms resolved the 16 isolates into seven sequence types (STs). We also demonstrate that MLST can be applied on uncultured isolates by direct PCR from cardiac valve tissue, and suggest this method presents a promising approach for epidemiological studies in this highly clonal organism. Phylogenetic and clustering analyses suggest that two of the seven STs form a distinct lineage within the population.     

Revised Genome Sequence of Burkholderia thailandensis MSMB43 with Improved Annotation

There is growing interest in discovery of novel bioactive natural products from Burkholderia thailandensis. Here we report a significantly improved genome sequence and reannotation of Burkholderia thailandensis MSMB43, which will facilitate the discovery of new natural products through genome mining and studies of the metabolic versatility of this bacterium.     

Complete Plastid Genome Sequence of the Brown Alga Undaria pinnatifida

In this study, we fully sequenced the circular plastid genome of a brown alga, Undaria pinnatifida. The genome is 130,383 base pairs (bp) in size; it contains a large single-copy (LSC, 76,598 bp) and a small single-copy region (SSC, 42,977 bp), separated by two inverted repeats (IRa and IRb: 5,404 bp). The genome contains 139 protein-coding, 28 tRNA, and 6 rRNA genes; none of these genes contains introns. Organization and gene contents of the U. pinnatifida plastid genome were similar to those of Saccharina japonica. There is a co-linear relationship between the plastid genome of U. pinnatifida and that of three previously sequenced large brown algal species. Phylogenetic analyses of 43 taxa based on 23 plastid protein-coding genes grouped all plastids into a red or green lineage. In the large brown algae branch, U. pinnatifida and S. japonica formed a sister clade with much closer relationship to Ectocarpus siliculosus than to Fucus vesiculosus. For the first time, the start codon ATT was identified in the plastid genome of large brown algae, in the atpA gene of U. pinnatifida. In addition, we found a gene-length change induced by a 3-bp repetitive DNA in ycf35 and ilvB genes of the U. pinnatifida plastid genome.     

Complete Genome Sequence of the Aerobic CO-Oxidizing Thermophile Thermomicrobium roseum

In order to enrich the phylogenetic diversity represented in the available sequenced bacterial genomes and as part of an “Assembling the Tree of Life” project, we determined the genome sequence of Thermomicrobium roseum DSM 5159. T. roseum DSM 5159 is a red-pigmented, rod-shaped, Gram-negative extreme thermophile isolated from a hot spring that possesses both an atypical cell wall composition and an unusual cell membrane that is composed entirely of long-chain 1,2-diols. Its genome is composed of two circular DNA elements, one of 2,006,217 bp (referred to as the chromosome) and one of 919,596 bp (referred to as the megaplasmid). Strikingly, though few standard housekeeping genes are found on the megaplasmid, it does encode a complete system for chemotaxis including both chemosensory components and an entire flagellar apparatus. This is the first known example of a complete flagellar system being encoded on a plasmid and suggests a straightforward means for lateral transfer of flagellum-based motility. Phylogenomic analyses support the recent rRNA-based analyses that led to T. roseum being removed from the phylum Thermomicrobia and assigned to the phylum Chloroflexi. Because T. roseum is a deep-branching member of this phylum, analysis of its genome provides insights into the evolution of the Chloroflexi. In addition, even though this species is not photosynthetic, analysis of the genome provides some insight into the origins of photosynthesis in the Chloroflexi. Metabolic pathway reconstructions and experimental studies revealed new aspects of the biology of this species. For example, we present evidence that T. roseum oxidizes CO aerobically, making it the first thermophile known to do so. In addition, we propose that glycosylation of its carotenoids plays a crucial role in the adaptation of the cell membrane to this bacterium''s thermophilic lifestyle. Analyses of published metagenomic sequences from two hot springs similar to the one from which this strain was isolated, show that close relatives of T. roseum DSM 5159 are present but have some key differences from the strain sequenced.