Genome Sequence of Lactobacillus crispatus ST1

Lactobacillus crispatus is a common member of the beneficial microbiota present in the vertebrate gastrointestinal and human genitourinary tracts. Here, we report the genome sequence of L. crispatus ST1, a chicken isolate displaying strong adherence to vaginal epithelial cells.Lactobacillus crispatus can persist in the vertebrate gastrointestinal tract and is among the most prevalent species of the Lactobacillus-dominated human vaginal microbiota (2, 9, 13, 14). It belongs to the so-called acidophilus group (3), which has attracted interest because some of its species are important factors in the production of fermented foods (12) and some can, at least transiently, colonize the human host (2, 9, 13, 14). Moreover, some specific strains, mainly L. acidophilus NCFM and L. johnsonii NCC 533, have received prominence as intestinal-health-promoting microbes (4). Although the genomes of seven members of the acidophilus complex have been sequenced to date (12), the genome sequences of L. crispatus and other predominant lactobacillar species in the urogenital flora have mostly remained obscure. Vaginal lactobacilli can have an important role in controlling the health of the host (2, 14). They can, for example, positively influence and stabilize the host''s vaginal microbiota via the production of compounds that are acidic or exert a direct inhibiting action toward pathogenic bacteria (2, 14). In addition to the antimicrobial compounds, the competitive exclusion of pathogens is another mechanism by which the host''s microbiota can be balanced (2). L. crispatus ST1 was originally isolated from the crop of a chicken, and PCR profiling of L. crispatus isolates has verified it to be an abundant colonizer of the chicken crop (6, 8). It also displays a strong protein-dependent adhesion to the epithelial cells of the human vagina and has been shown to inhibit the adhesion of avian pathogenic Escherichia coli (6, 7).The genome was sequenced (18× coverage) using a 454 pyrosequencer with GS FLX chemistry (Roche). The contig order was confirmed and gaps were filled by sequencing PCR fragments from the genomic DNA template using ABI 3730 and Big Dye chemistry (Applied Biosystems). Genomic data were processed using the Staden Package (11) and gsAssembler (Roche). Coding sequences (CDSs) were predicted using Glimmer3 (5) followed by manual curation of the start sites. The remaining intergenic regions were reanalyzed for missed CDSs by using BlastX (1). Annotation transfer was performed based on a BlastP search, followed by Blannotator analysis using default settings (http://ekhidna.biocenter.helsinki.fi/poxo/blannotator) and manual verification. Orthologous groups between the different lactobacillar proteomes were identified using OrthoMCL (10).The genome of L. crispatus ST1 consists of a single circular chromosome 2.04 Mbp in size, with an overall G+C content of 37%, without any plasmids. There are 64 tRNA genes, 4 rRNA operons, and 2 CRISPR loci. Out of the 2,024 predicted CDSs, a putative function was assigned to 77%, whereas 10% of the CDSs were annotated as conserved and 13% as novel. Based on the orthologous grouping, 302 (15%) of the CDSs encoded by ST1 have no detectable homologs in any of the Lactobacillus proteomes published to date.     

Genome Sequence of Pantoea ananatis LMG20103, the Causative Agent of Eucalyptus Blight and Dieback

Pantoea ananatis is a Gram-negative plant pathogen that causes disease on a broad range of host plants, including pineapple, maize, rice, onion, melons, and Eucalyptus, and has been implicated in several cases of human disease. Here, we report the genome sequence of P. ananatis LMG20103 isolated from diseased Eucalyptus in South Africa.Pantoea ananatis belongs within the Enterobacteriaceae family and the genus Pantoea, which hosts a number of important plant pathogens and clinically relevant species. It has been isolated from a range of plant hosts worldwide, and its ability to cause disease in several important staple food crops, including maize, onion and rice, has led to its being recognized as an important emerging plant pathogen (3). In South Africa, P. ananatis causes severe blight and dieback of Eucalyptus hybrid and clone seedlings/cuttings, which affects the pulp and paper industry (4). For this reason, we chose to sequence the genome of the highly virulent Eucalyptus-pathogenic strain P. ananatis LMG20103. Previous analyses have indicated that this strain lacks genes encoding type II and type III secretion systems (3), which represent two major pathogenicity determinants in most plant- and animal-pathogenic bacteria (7). The genome sequence may resolve how P. ananatis LMG20103 causes disease in the absence of these secretion mechanisms. This also constitutes the first genome of a member of the genus Pantoea and the first bacterial phytopathogen genome sequenced in Africa.The genome was sequenced at Inqaba Biotec, South Africa, using a Roche 454 GS20 sequencer. A total of 991,246 reads, 80 to 120 nucleotides in length, yielding 97,152,414 nucleotides were sequenced to give a 20.9-times coverage. A draft assembly of 117 contigs was produced using Newbler Assembler 2.00.00. A scaffolding approach with 12 genomes of closely related members of the Enterobacteriaceae, combined with PCR gap closure, was utilized for final genome assembly. Protein-coding sequences (CDSs) were predicted using a combination of Glimmer 3.0 (5) and AMIGene (2). Annotation was performed using the automated annotation engines BASys (9) and MaGe (8) and validated manually against the NCBI (1) and COG databases (6).The genome of P. ananatis LMG20103 consists of a single chromosome ∼4.69 Mb in size, with a G+C content of 53.69%. There are 4,272 predicted CDSs with an average length of 969 nucleotides, giving a coding intensity of 88.27%. Seventy-two tRNAs and 6 rRNA operons were identified. Putative functions have been assigned to 73% of the encoded proteins. The presence of potential homologues to genes involved in pathogenesis in other plant- and animal-pathogenic bacteria suggests that the genome sequence of P. ananatis LMG20103 will provide considerable insights into the pathogenesis of this organism on Eucalyptus and into diseases caused by its plant- and animal-pathogenic relatives within the genus Pantoea.     

Complete Genome Sequence of the Representative γ-Hexachlorocyclohexane-Degrading Bacterium Sphingobium japonicum UT26

Sphingobium japonicum strain UT26 utilizes γ-hexachlorocyclohexane (γ-HCH), a man-made chlorinated pesticide that causes serious environmental problems due to its toxicity and long persistence, as a sole source of carbon and energy. Here, we report the complete genome sequence of UT26, which consists of two chromosomes and three plasmids. The 15 lin genes involved in γ-HCH degradation are dispersed on the two chromosomes and one of the three plasmids.γ-Hexachlorocylohexane (γ-HCH) is a man-made insecticide that has caused serious environmental problems worldwide (9). Although only several decades have passed since the initial release of γ-HCH into the environment, many γ-HCH-degrading bacterial strains have been isolated (9), suggesting that γ-HCH-degrading ability was acquired by such strains within a short period (14). Sphingobium japonicum UT26 was isolated from upland γ-HCH-contaminated soil in Japan and utilizes γ-HCH as its sole source of carbon and energy (8, 16). In UT26, 15 lin genes are involved in γ-HCH utilization (2, 11, 12). Our clarification of the complete genome sequence of this strain is expected to provide insights into the mechanisms by which bacteria adapt to xenobiotics.One S. japonicum UT26 colony was designated UT26S (NBRC 101211), and its complete genome sequence was determined by a whole-genome shotgun sequencing strategy using the Sanger method (10, 15, 17). The sequences of ca. 92,400 reads were assembled by the Phrap and CONSED assembly tools (4, 5, 7), and the gaps between contigs were closed by sequencing PCR products which were amplified from genomic DNA using the appropriate primers. The prediction and annotation of open reading frames (ORFs) were performed using Glimmer3 (1), BLASTP, the In Silico Molecular Cloning software suite (In Silico Biology, Inc.), and the GenomeMatcher software (13). Nontranslating genes were predicted using the Rfam, tRNAscan-SE, and ARAGORN programs.The genome of S. japonicum UT26S consists of two circular chromosomes (Chr), Chr 1 (3,514,822 bp, 64.8% G+C, 3,529 ORFs) and Chr 2 (681,892 bp, 65.9% G+C, 589 ORFs), and three circular plasmids, pCHQ1 (190,974 bp, 63.0% G+C, 224 ORFs), pUT1 (31,776 bp, 63.7% G+C, 44 ORFs), and pUT2 (5,398 bp, 61.0% G+C, 8 ORFs). Chr 1 and Chr 2 have one and two copies, respectively, of rRNA operons. Fifty-one and 4 tRNA genes were located on Chr 1 and Chr 2, respectively. One hundred ninety-six out of 206 bacterial essential genes proposed by Gil et al. (6) were all located on Chr 1, indicating that this is clearly a “main” chromosome. The 15 lin genes for γ-HCH degradation are dispersed on Chr 1, Chr 2, and pCHQ1. Comparison of the UT26S genome with those of five other sphingomonad (TM1) strains revealed that the lin genes (linA, linB, linC, linRED, and linF) specific for the conversion of γ-HCH to β-ketoadipate are located in the DNA regions unique to the UT26S genome. On the other hand, linGHIJ for the β-ketoadipate pathway (12) and linKLMN for the ABC transporter system (3) are located in conserved genomic regions of these sphingomonads. Based on these results, we propose a model in which UT26S was established by recruiting the specific lin genes into an ancestral strain having core functions of sphingomonads.     

Complete Genome Sequence of Halalkalicoccus jeotgali B3T,an Extremely Halophilic Archaeon

Halalkalicoccus jeotgali B3^T, isolated from salt-fermented seafood from South Korea, is an extremely halophilic archaeon belonging to the family Halobacteriaceae. Here, we present the complete genome sequence of the type strain H. jeotgali B3^T (3,698,650 bp, with a G+C content of 62.5%), which consists of one chromosome and six plasmids. This is the first complete genome sequence of the Halalkalicoccus species.Extremely halophilic archaea (haloarchaea) are adapted to hypersaline environments and grow optimally in NaCl solutions of 2.6 M or higher (12). These haloarchaea are classified within the family Halobacteriaceae in the order Halobacteriales; currently, this family comprises 28 genera (3), and only 11 complete genome sequences in Halobacteriaceae have been reported. In a study of archaeal diversity in salt-fermented small shrimp or shellfish from South Korea, our laboratory isolated and characterized 5 novel, extremely halophilic archaeal strains of Halobacteriaceae. These strains included Natronococcus jeotgali (9), Halalkalicoccus jeotgali (11), Halorubrum cibi (7), Haloterrigena jeotgali (10) and Haladaptatus cibarius (8). We have now sequenced the genome of Halalkalicoccus jeotgali B3^T; genome sequencing had not been completed or initiated for any strain in this genus when our sequencing project was begun. The genus Halalkalicoccus currently contains only two species, Halalkalicoccus tibetensis (13) and H. jeotgali, and these species exhibit 98.6% gene sequence similarity in their 16S rRNA. The genome of H. jeotgali B3^T is the first of this genus to be sequenced.The complete genome sequence of H. jeotgali B3^T was determined by a whole-genome shotgun strategy using Roche 454 GS (FLX Titanium) pyrosequencing (898,168 reads totaling ∼348 Mb; ∼94-fold coverage of the genome) and a fosmid library (514 reads totaling ∼680 kb) at the Genome Resource Center, KRIBB (Korea Research Institute of Bioscience and Biotechnology). Genome sequences from pyrosequencing were processed by Roche''s software according to the manufacturer''s instructions, and sequences from the fosmid library were processed by PESTAS (6). A total of 898,196 reads were assembled using Newbler Assembler 2.3 (454 Life Science), which generated 54 large contigs (>100 bp in size) with bases having quality scores of 40 and above. The gaps between contigs were closed by primer walking and sequencing of PCR products across the gaps. The annotation was done by merging results obtained from the RAST (Rapid Annotation using Subsystem Technology) pipeline (1), Glimmer 3.02 (2), tRNAscan-SE 1.21 (5), and RNAmmer 1.2 (4).The H. jeotgali B3^T genome is 3,698,650 bases long with a 62.5% G+C content. The chromosome consists of a single circular chromosome (2,809,118 bp, with a G+C content of 65.0%) and six plasmids (406,285 bp, 55.3%; 363,534 bp, 54.2%; 44,576 bp, 58.9%; 44,459 bp, 54.9%; 23,727 bp, 47.6%; 6,951 bp, 60.6%). The genome contains 3,860 predicted coding sequences and 52 RNA genes (determined using RAST). The chromosome is predicted to contain 3,101 coding sequences with a coding intensity of 90.0%, including 47 tRNA genes, 1 5S rRNA gene, 1 16S rRNA gene, and 1 23S rRNA gene. The largest plasmid contains 466 coding sequences with a coding intensity of 81.2% and 2 tRNA genes, while the other five plasmids contain 425, 44, 48, 29, and 5 coding sequences with coding intensities of 80.2%, 84.2%, 83.0%, 69.6%, and 22.8%, respectively (determined using Glimmer3). More detailed analysis of this genome and comparative analysis with other haloarchaea will provide further insight into the genomic differences and metabolism of the extremely halophilic archaea.     

Abundance of Novel and Diverse tfdA-Like Genes,Encoding Putative Phenoxyalkanoic Acid Herbicide-Degrading Dioxygenases,in Soil

Phenoxyalkanoic acid (PAA) herbicides are widely used in agriculture. Biotic degradation of such herbicides occurs in soils and is initiated by α-ketoglutarate- and Fe²⁺-dependent dioxygenases encoded by tfdA-like genes (i.e., tfdA and tfdAα). Novel primers and quantitative kinetic PCR (qPCR) assays were developed to analyze the diversity and abundance of tfdA-like genes in soil. Five primer sets targeting tfdA-like genes were designed and evaluated. Primer sets 3 to 5 specifically amplified tfdA-like genes from soil, and a total of 437 sequences were retrieved. Coverages of gene libraries were 62 to 100%, up to 122 genotypes were detected, and up to 389 genotypes were predicted to occur in the gene libraries as indicated by the richness estimator Chao1. Phylogenetic analysis of in silico-translated tfdA-like genes indicated that soil tfdA-like genes were related to those of group 2 and 3 Bradyrhizobium spp., Sphingomonas spp., and uncultured soil bacteria. Soil-derived tfdA-like genes were assigned to 11 clusters, 4 of which were composed of novel sequences from this study, indicating that soil harbors novel and diverse tfdA-like genes. Correlation analysis of 16S rRNA and tfdA-like gene similarity indicated that any two bacteria with D > 20% of group 2 tfdA-like gene-derived protein sequences belong to different species. Thus, data indicate that the soil analyzed harbors at least 48 novel bacterial species containing group 2 tfdA-like genes. Novel qPCR assays were established to quantify such new tfdA-like genes. Copy numbers of tfdA-like genes were 1.0 × 10⁶ to 65 × 10⁶ per gram (dry weight) soil in four different soils, indicating that hitherto-unknown, diverse tfdA-like genes are abundant in soils.Phenoxyalkanoic acid (PAA) herbicides such as MCPA (4-chloro-2-methyl-phenoxyacetic acid) and 2,4-D (2,4-dichlorophenoxyacetic acid) are widely used to control broad-leaf weeds in agricultural as well as nonagricultural areas (19, 77). Degradation occurs primarily under oxic conditions in soil, and microorganisms play a key role in the degradation of such herbicides in soil (62, 64). Although relatively rapidly degraded in soil (32, 45), both MCPA and 2,4-D are potential groundwater contaminants (10, 56, 70), accentuating the importance of bacterial PAA herbicide-degrading bacteria in soils (e.g., references 3, 5, 6, 20, 41, 59, and 78).Degradation can occur cometabolically or be associated with energy conservation (15, 54). The first step in the degradation of 2,4-D and MCPA is initiated by the product of cadAB or tfdA-like genes (29, 30, 35, 67), which constitutes an α-ketoglutarate (α-KG)- and Fe²⁺-dependent dioxygenase. TfdA removes the acetate side chain of 2,4-D and MCPA to produce 2,4-dichlorophenol and 4-chloro-2-methylphenol, respectively, and glyoxylate while oxidizing α-ketoglutarate to CO₂ and succinate (16, 17).Organisms capable of PAA herbicide degradation are phylogenetically diverse and belong to the Alpha-, Beta-, and Gammproteobacteria and the Bacteroidetes/Chlorobi group (e.g., references 2, 14, 29-34, 39, 60, 68, and 71). These bacteria harbor tfdA-like genes (i.e., tfdA or tfdAα) and are categorized into three groups on an evolutionary and physiological basis (34). The first group consists of beta- and gammaproteobacteria and can be further divided into three distinct classes based on their tfdA genes (30, 46). Class I tfdA genes are closely related to those of Cupriavidus necator JMP134 (formerly Ralstonia eutropha). Class II tfdA genes consist of those of Burkholderia sp. strain RASC and a few strains that are 76% identical to class I tfdA genes. Class III tfdA genes are 77% identical to class I and 80% identical to class II tfdA genes and linked to MCPA degradation in soil (3). The second group consists of alphaproteobacteria, which are closely related to Bradyrhizobium spp. with tfdAα genes having 60% identity to tfdA of group 1 (18, 29, 34). The third group also harbors the tfdAα genes and consists of Sphingomonas spp. within the alphaproteobacteria (30).Diverse PAA herbicide degraders of all three groups were identified in soil by cultivation-dependent studies (32, 34, 41, 78). Besides CadAB, TfdA and certain TfdAα proteins catalyze the conversion of PAA herbicides (29, 30, 35). All groups of tfdA-like genes are potentially linked to the degradation of PAA herbicides, although alternative primary functions of group 2 and 3 TfdAs have been proposed (30, 35). However, recent cultivation-independent studies focused on 16S rRNA genes or solely on group 1 tfdA sequences in soil (e.g., references 3-5, 13, and 41). Whether group 2 and 3 tfdA-like genes are also quantitatively linked to the degradation of PAA herbicides in soils is unknown. Thus, tools to target a broad range of tfdA-like genes are needed to resolve such an issue. Primers used to assess the diversity of tfdA-like sequences used in previous studies were based on the alignment of approximately 50% or less of available sequences to date (3, 20, 29, 32, 39, 47, 58, 73). Primers specifically targeting all major groups of tfdA-like genes to assess and quantify a broad diversity of potential PAA degraders in soil are unavailable. Thus, the objectives of this study were (i) to develop primers specific for all three groups of tfdA-like genes, (ii) to establish quantitative kinetic PCR (qPCR) assays based on such primers for different soil samples, and (iii) to assess the diversity and abundance of tfdA-like genes in soil.     

Occurrence and Molecular Characteristics of Methicillin-Resistant Staphylococcus aureus and Methicillin-Resistant Staphylococcus pseudintermedius in an Academic Veterinary Hospital

Recently, methicillin-resistant Staphylococcus aureus (MRSA) and methicillin-resistant Staphylococcus pseudintermedius (MRSP) have been increasingly isolated from veterinarians and companion animals. With a view to preventing the spread of MRSA and MRSP, we evaluated the occurrence and molecular characteristics of each in a veterinary college. MRSA and MRSP were isolated from nasal samples from veterinarians, staff members, and veterinary students affiliated with a veterinary hospital. Using stepwise logistic regression, we identified two factors associated with MRSA carriage: (i) contact with an identified animal MRSA case (odds ratio [OR], 6.9; 95% confidence interval [95% CI], 2.2 to 21.6) and (ii) being an employee (OR, 6.2; 95% CI, 2.0 to 19.4). The majority of MRSA isolates obtained from individuals affiliated with the veterinary hospital and dog patients harbored spa type t002 and a type II staphylococcal cassette chromosome mec (SCCmec), similar to the hospital-acquired MRSA isolates in Japan. MRSA isolates harboring spa type t008 and a type IV SCCmec were obtained from one veterinarian on three different sampling occasions and also from dog patients. MRSA carriers can also be a source of MRSA infection in animals. The majority of MRSP isolates (85.2%) carried hybrid SCCmec type II-III, and almost all the remaining MRSP isolates (11.1%) carried SCCmec type V. MRSA and MRSP were also isolated from environmental samples collected from the veterinary hospital (5.1% and 6.4%, respectively). The application of certain disinfection procedures is important for the prevention of nosocomial infection, and MRSA and MRSP infection control strategies should be adopted in veterinary medical practice.Methicillin-resistant Staphylococcus aureus (MRSA) is an important cause of nosocomial infections in human hospitals. The prevalence of hospital-acquired MRSA (HA-MRSA) infection among inpatients in intensive care units (ICUs) continues to increase steadily in Japan. Recently, cases of community-acquired MRSA (CA-MRSA) have been documented in persons without an established risk factor for HA-MRSA infection (14, 32, 36, 49).There has also been an increase in the number of reports of the isolation of MRSA from veterinarians and companion animals (5, 21, 23-26, 28, 31, 34, 38, 44, 50, 51, 53). Values reported for the prevalence of MRSA among veterinary staff include 17.9% in the United Kingdom (21), 10% in Japan (38), 3.9% in Scotland (13), and 3.0% in Denmark (28). Loeffler et al. reported that the prevalence of MRSA among dog patients and healthy dogs owned by veterinary staff members was 8.9% (21). In Japan, an MRSA isolate was detected in only one inpatient dog (3.8%) and could not be detected in any of 31 outpatient dogs (38). In the United States, MRSA isolates were detected in both dog (0.1%) and cat (0.1%) patients (31). The prevalence of MRSA among healthy dogs has been reported to be 0.7% (5). Hanselman et al. suggested that MRSA colonization may be an occupational risk for large-animal veterinarians (12). Recently, Burstiner et al. reported that the frequency of MRSA colonization among companion-animal veterinary personnel was equal to the frequency among large-animal veterinary personnel (6).In addition, other methicillin-resistant coagulase-positive staphylococci (MRCPS), such as methicillin-resistant Staphylococcus pseudintermedius (MRSP) and methicillin-resistant Staphylococcus schleiferi (MRSS), isolated from dogs, cats, and a veterinarian have been reported (11, 31, 38, 40, 52). MRSP isolates have also been detected among inpatient dogs (46.2%) and outpatient dogs (19.4%) in a Japanese veterinary teaching hospital (38). In Canada, however, MRSP and MRSS isolates were detected in only 2.1% and 0.5% of dog patients, respectively (11).Methicillin-resistant staphylococci produce penicillin-binding protein 2′, which reduces their affinity for β-lactam antibiotics. This protein is encoded by the mecA gene (48), which is carried on the staphylococcal cassette chromosome mec (SCCmec). SCCmec is a mobile genetic element characterized by the combination of the mec and ccr complexes (16), and it is classified into subtypes according to differences in the junkyard regions (43). SCCmec typing can be used as a molecular tool (22, 27, 30, 33, 36, 55) for examining the molecular epidemiology of methicillin-resistant staphylococci.In this study, we investigated the occurrence and characteristics of MRCPS isolates in a veterinary hospital in order to establish the transmission route of MRCPS in a veterinary hospital and with a view to preventing the spread of MRCPS infection. In addition, we evaluated the factors associated with MRCPS. Further, as Heller et al. have reported the distribution of MRSA within veterinary hospital environments and suggested the necessity to review cleaning protocols of hospital environments (13), we also attempted to isolate MRCPS from environmental samples collected in a veterinary hospital for an evaluation of MRSA transmission cycle though environmental surfaces in the veterinary hospital.     

Complete Genome Sequence of Bacillus thuringiensis Mutant Strain BMB171

Bacillus thuringiensis has been widely used as a biopesticide for a long time. Here we report the finished and annotated genome sequence of B. thuringiensis mutant strain BMB171, an acrystalliferous mutant strain with a high transformation frequency obtained and stocked in our laboratory.Bacillus thuringiensis is an insect pathogen which is widely used as a biopesticide due to its various endogenous crystal proteins and spores (12). To improve the virulence and practical effectiveness of B. thuringiensis, genetic transformation of different genes with beneficial traits is a fundamental procedure. Simultaneously, genetic transformation can facilitate functional genomic research. However, wild-type strains are not suitable to be used as recipient strains because of low transformation efficiency. This obstacle is mainly caused by the thick cell wall layer of B. thuringiensis together with multiple plasmids inside the cell, which harbor genes encoding insecticidal crystal proteins. We used the method of elevating the growth temperature and adding 0.05% sodium dodecyl sulfate to treat several parental strains and finally obtained mutant strain BMB171, with no resident plasmid, from wild-type crystalliferous strain YBT-1463 (9). The electrotransformation frequency of mutant BMB171 could reach up to 10⁷ transformants/μg DNA after optimization of the electrotransformation parameters (7), which was 4.8 × 10⁴-fold higher than that of the parental strain (8). Moreover, mutant strain BMB171 exhibited the same characteristics as YBT-1463, such as metabolic abilities and growth properties, as well as sensitivity to 10 antibiotics (8). Of course, BMB171 could produce parasporal crystals with characteristic geometric shapes through the expression of relevant cry genes carried by plasmids (7). Thus, B. thuringiensis mutant strain BMB171 has become a major recipient strain and is widely used for insecticidal crystal protein-encoding gene expression (14, 15), cell surface display (10, 13), gene function and regulation researches (2, 5), etc.The B. thuringiensis mutant strain BMB171 genome was sequenced by using a massive parallel pyrosequencing technology (454 GS-FLX). A total of 448,963 high-quality reads with an average read length of 391 bp were produced, providing about 32-fold coverage of the genome. Assembly was performed using the Newbler software of the 454 suite package (454 Life Sciences), which resulted in 193 large (defined as >500 bp) contigs. The relationship of contigs was determined by multiplex PCR, and gaps were filled through sequencing of PCR products by primer walking or shotgun sequencing with an ABI 3730 sequencer. The Phred/Phrap/Consed software package (3) was used for final sequence assembly and quality assessment. Protein-coding genes were predicted by combining the results of Glimmer 3.02 (1) and ZCURVE (4), followed by manual inspection. Both tRNA and rRNA genes were identified by tRNAscan-SE (11) and RNAmmer (6), respectively. Functional annotation was performed by searching against a protein database of the microbial genome developed in house.The 5.64-Mb genome of B. thuringiensis mutant strain BMB171 contains two replicons: a circular chromosome (5.33 Mb) encoding 5,088 open reading frames (ORFs) and a circular plasmid (0.31 Mb), which is named pBMB171, encoding 276 predicted ORFs. The G+C content of the chromosome is 35.3%, while that of the plasmid is 33.3%. The mutant strain BMB171 genome encodes 104 tRNAs and 14 rRNA operons. A previous study indicated that BMB171 is a plasmid-free mutant (9); however, our sequencing results demonstrated that a large plasmid still remains. The reason why the plasmid was not detected previously might be its large size and low copy number. We did not find any crystal protein genes in either chromosome or plasmid sequences, which was consistent with previous observations (9).In summary, the complete B. thuringiensis mutant strain BMB171 genome provides a better-defined genetic background for gene expression and regulation studies, especially crystal protein production and metabolic network construction.     

Poultry-Associated Salmonella enterica subsp. enterica Serovar 4,12:d:? Reveals High Clonality and a Distinct Pathogenicity Gene Repertoire

A European baseline survey during the years 2005 and 2006 has revealed that the monophasic Salmonella enterica subsp. enterica serovar 4,12:d:− was, with a prevalence of 23.6%, the most frequently isolated serovar in German broiler flocks. In Denmark and the United Kingdom, its serovar prevalences were 15.15% and 2.8%, respectively. Although poultry is a major source of human salmonellosis, serovar 4,12:d:− is rarely isolated in humans (approximately 0.09% per year). Molecular typing studies using pulsed-field gel electrophoresis and DNA microarray analysis show that the serovar is highly clonal and lacks genes with known contributions to pathogenicity. In contrast to other poultry-associated serovars, all strains were susceptible to 17 antimicrobial agents tested and did not encode any resistance determinant. Furthermore, serovar 4,12:d:− lacked the genes involved in galactonate metabolism and in the glycolysis and glyconeogenesis important for energy production in the cells. The conclusion of the study is that serovar 4,12:d:− seems to be primarily adapted to broilers and therefore causes only rare infections in humans.Salmonella spp. are major zoonotic food-borne pathogens which cause outbreaks and sporadic cases of gastroenteritis in humans worldwide (12). Depending on the serovar, cases of salmonellosis can differ substantially in severity (13). The primary sources of salmonellosis are food-producing animals, such as poultry, pigs, and cattle (30). The pathogen is spread by trade in animals and nonheated animal food products (10).A European baseline survey on the prevalence of Salmonella in commercial broiler flocks of Gallus gallus in 2005 and 2006 showed that in the European Union (EU), 23.7% of the broiler flocks were Salmonella positive (8). However, the Salmonella prevalences and serovar distribution varied widely among the EU member states. The five most frequently isolated Salmonella enterica serovars in Europe were those classically observed, like serovar Enteritidis (33.8%), serovar Infantis (22.0%), serovar Mbandaka (8.1%), serovar Hadar (3.7%), and serovar Typhimurium (3.0%). In Germany, the flock prevalence of Salmonella was 15.0% among the 377 broiler flocks investigated. In contrast to the well-known serovars described above, the predominating serovar was monophasic serovar 4,12:d:−, with a prevalence of 23.6%. This serovar was also isolated in Denmark and the United Kingdom, with prevalences of 15.2% and 2.8%, respectively.The German Salmonella National Reference Laboratory (NRL-Salmonella) has received 818 isolates of this serovar between 1998 and 2007, with peaks in 2001 (240 isolates) and 2004 (160 isolates), for diagnosis. Since 2005, the number has doubled annually, and the serovar obviously established itself well in poultry production lines. These isolates were found mostly in broilers (78%), occasionally in turkeys (11.6%) and feedstuff (8.4%), and rarely in pigs (1.3%) and cattle (0.6%). In contrast, infections in humans are only sporadic. During the last 10 years (1998 to 2007), the National Reference Centre for Salmonellae and Other Enterics located at the Robert-Koch Institute, Wernigerode branch, has received 55 strains of this serovar from sporadic human cases of salmonellosis and carriers in Germany (W. Rabsch, personal communication). Similarly, in Denmark, only two isolates in 1993 and in 2002 were isolated from humans (E. Møller Nielsen, Statens Serum Institut, Copenhagen, Denmark, personal communication).Subtyping food-borne pathogens is an approach often applied to facilitate the epidemiological investigation of outbreaks of gastrointestinal disease and to identify the source of entry into the food chain. Several molecular-based tools have been developed to type bacteria genotypically. Pulsed-field gel electrophoresis (PFGE) is currently the method of choice for the molecular subtyping of Salmonella serovars. It has been proven to be a useful discriminatory method which was standardized by the PulseNet Consortium (9). However, although this approach is certainly valuable, it does not reveal data on the gene repertoire and biological properties of a strain. To overcome this weakness, whole-genome DNA microarrays have successfully been applied in comparative genomic hybridizations for Salmonella (7, 24, 25). However, whole-genome arrays reflect only one genome of one strain. Because of many serovar or strain genome variations described for Salmonella, thematic arrays were developed, such as arrays specially targeting genes involved in resistance profiles (2, 17, 32), phage types (23), or serovars (33, 35). A condensed selection of 109 various Salmonella genetic markers comprising the detection of flagellar and somatic antigen-encoding genes, important virulence genes, phage-associated genes, and antibiotic resistance determinants have been used to show the usefulness of DNA microarrays for the discriminative characterization of Salmonella serovars (18).In this study, we elucidate the contradicting situation between the high prevalence in broilers and source attribution of broiler meat for humans and the low infection rates in humans of the serovar 4,12:d:− by genotypic characterization using PFGE and DNA microarray to determine the clonality, the pathogenic gene repertoire, and resistance determinants. These data give basic information to discuss the hazard potential of this serovar for humans. For that purpose, a new prototype of a Salmonella DNA microarray comprising 281 60-mer oligonucleotide probes was developed and validated in house.     

Completed Genome Sequence of the Anaerobic Iron-Oxidizing Bacterium Acidovorax ebreus Strain TPSY

Acidovorax ebreus strain TPSY is the first anaerobic nitrate-dependent Fe(II) oxidizer for which there is a completed genome sequence. Preliminary protein annotation revealed an organism optimized for survival in a complex environmental system. Here, we briefly report the completed and annotated genome sequence of strain TPSY.Microorganisms from diverse anoxic environments are capable of nitrate-dependent Fe(II) oxidation at circumneutral pH (4, 11, 17, 18, 20, 21). Despite their geochemical importance (22), little is known of the underlying biochemical and genetic mechanisms. Genome sequencing of several nitrate-dependent Fe(II) oxidizers will provide insight into this process. By comparing Fe(II) oxidation mechanisms in various organisms, we hope to identify both the conserved and disparate aspects of the metabolism. The genome of Acidovorax ebreus strain TPSY is the first of these to be sequenced.Strain TPSY is a motile, Gram-negative facultative anaerobe isolated from groundwater collected from the U.S. Department of Energy Integrated Field Research Challenge site at Oak Ridge, TN. Growth experiments performed as previously described (21) revealed TPSY''s incapacity for lithoautotrophic growth, which was supported by a lack of genes in the genome encoding any known CO₂ fixation pathways. TPSY did grow mixotrophically with Fe(II) as the electron donor and a 0.1 mM acetate carbon source. 16S rRNA gene sequence analysis placed TPSY in the class Betaproteobacteria with 99.8% similarity to Acidovorax sp. strain JS42 in the family Comamonadaceae.The completed genome consisted of a single circular chromosome of 3,796,573 bp with an average 66.8% G+C content. A total of 3,479 protein-encoding genes were predicted, and 34 (0.98%) had no similarity to public database sequences. Sequencing performed at the Department of Energy Joint Genome Institute (JGI) used Sanger sequencing and 454 pyrosequencing to a depth of 20× coverage. All JGI library construction and sequencing techniques can be found at http://www.jgi.doe.gov/. Sequence assembly, quality assessment, and annotation were performed using the software Phred/Phrap/Consed (www.phrap.com) (6-8), Dupfinisher (10), CRITICA (2), GLIMMER, and GENERATION (5) and the JGI Integrated Microbial Genomes site (12). The completed genome sequence contained 33,341 reads and had an average of ninefold coverage per base and an error rate of <1 in 100,000.TPSY was named in part for its meandering motility, and its genome confirmed the twitching phenotype with the presence of pilT, pilU, and a complete set of flagellar and chemotaxis genes. The ability of TPSY to oxidize simple alcohols and acids with oxygen or nitrate respiration was confirmed by the genome. In addition, biosynthetic pathways for all amino acids except tyrosine and phenylalanine were present. No homologues of chorismate mutase (EC 5.4.99.5), an enzyme required for tyrosine and phenylalanine anabolism, were identified. The genome contained both intact Embden-Meyerhof-Parnas and Entner-Doudoroff pathways, in addition to a pentose phosphate pathway and a trichloroacetic acid cycle.In support of its facultative anaerobicity, a complete set of genes for denitrification and three different terminal oxidases (cytochrome aa₃, cbb₃, and cytochrome d quinol oxidase) were present. The cbb₃ and cytochrome d oxidases, with their high oxygen affinity, putatively enable survival in microaerobic environments (14).TPSY had sequences encoding 30 transposases, 11 integrases, and 11 phage/prophage-related genes. A region of particular interest putatively conferred resistance to lead, arsenate, and mercury: pbrRATARTBC, arsRDAB, and merRPCADE. Evidence suggests horizontal transfer and insertion of this region, as it was flanked on the 5′ end by λ prophage-related genes and the 3′ end encoded a putative Tn21 transposase. Phenotypic studies by the method of Wang et al. (19) revealed MICs of 16 μM phenylmercuric acetate and 250 μM MgCl₂. TPSY was also capable of growth in the presence of arsenate (10 mM) but did not use it as an electron acceptor.Related to phage infection, one CRISPR (clustered, regularly interspaced, short palindromic repeats) region (3, 16) was predicted. The core proteins, the cas1 and cas2 genes, and a csn1 gene formed the CRISPR subtype Nmeni, which is associated with vertebrate pathogens and commensals (9). However, the lack of typical pathogenic type I or III secretion systems such as the hec cluster of Dickeya chrysanthemi (15) or the inv/spa system of Salmonella enterica serovar Typhimurium (13) indicated that TPSY would probably not exhibit a pathogenic lifestyle.     

Complete Genome Sequence of Mycoplasma hyorhinis Strain HUB-1

Mycoplasma hyorhinis is generally considered a swine pathogen yet is most commonly found infecting laboratory cell lines. An increasing body of evidence suggests that chronic infections with M. hyorhinis may cause oncogenic transformation. Here, we announce the complete genome sequence of M. hyorhinis strain HUB-1.Mycoplasma hyorhinis is generally considered to be a swine pathogen causing lung lesions, inflammation in the chest and abdominal lining, and arthritis (8). This agent also frequently contaminates laboratory cell cultures, impinging on many aspects of biological research (3). Recent studies have demonstrated that M. hyorhinis infections induce a malignant phenotype in human prostate (7) and gastric (4) cells, suggesting that M. hyorhinis infections are associated with oncogenic transformation. These properties of M. hyorhinis have increased its profile to researchers. The complete genome sequence of this microbe has yet to be determined.We sequenced the genome of M. hyorhinis strain HUB-1, a pathogenic strain isolated from the respiratory tract of swine. Whole-genome sequencing was performed by combining GS FLX (6) and Solexa paired-end sequencing technologies (1). Genomic libraries containing 3-kb inserts were constructed, and 308,604 reads (79.7% paired end) were produced using the GS FLX system, giving 65.9-fold coverage of the genome. About 93.4% of reads were assembled into one large scaffold using Newbler software (454 Life Sciences, Branford, CT). A total of 822,579 reads were generated using an Illumina Solexa Genome Analyzer IIx and were mapped to the scaffold using the Burrows-Wheeler alignment (BWA) tool (5). Gaps were filled by local assembly of the Solexa/Roche 454 reads or by sequencing PCR products by using an ABI 3730 capillary sequencer. Open reading frames containing more than 30 amino acid residues were predicted using Glimmer 3.0 (2) and verified by comparison with six other closely related genome sequences.The complete genome of M. hyorhinis HUB-1 consists of an 839,615-bp single circular chromosome with an average G+C content of 25.88%. A total of 654 protein-encoding genes are predicted. The average protein size is 364 amino acids, and the mean coding percentage is 85.2%. The genome includes 30 tRNA genes, and only a single copy of the 16S-23S rRNA operon can be found. The 5S rRNA operon is separate from the 16S-23S rRNA operon. Protein secretion occurs through a truncated membrane protein secretion system, consisting of SecA, SecD, SecY, PrsA, DnaK, Tig, and LepA. Additionally, 20 pseudogenes, which become truncated or inactivated, are identified in the genome.M. hyorhinis contains a special variable lipoprotein (Vlp) system that constitutes its major coat protein (9) and provides a mutational strategy for evasion of the host immune system. Different M. hyorhinis strains carry a variable number of vlp genes (9). M. hyorhinis HUB-1 is characterized to contain seven vlp genes displayed in the order 5′-vlpD-vlpE-vlpF-insertion sequence (IS)-vlpG-vlpA-IS-vlpB-vlpC-3′.This is the first complete genome sequence of M. hyorhinis, and its availability will provide a better-defined genetic background for future studies of gene expression and regulation.     

Optical Mapping Reveals a Large Genetic Inversion between Two Methicillin-Resistant Staphylococcus aureus Strains

Staphylococcus aureus is a highly versatile and evolving bacterium of great clinical importance. S. aureus can evolve by acquiring single nucleotide polymorphisms and mobile genetic elements and by recombination events. Identification and location of novel genomic elements in a bacterial genome are not straightforward, unless the whole genome is sequenced. Optical mapping is a new tool that creates a high-resolution, in situ ordered restriction map of a bacterial genome. These maps can be used to determine genomic organization and perform comparative genomics to identify genomic rearrangements, such as insertions, deletions, duplications, and inversions, compared to an in silico (virtual) restriction map of a known genome sequence. Using this technology, we report here the identification, approximate location, and characterization of a genetic inversion of ∼500 kb of a DNA element between the NRS387 (USA800) and FPR3757 (USA300) strains. The presence of the inversion and location of its junction sites were confirmed by site-specific PCR and sequencing. At both the left and right junction sites in NRS387, an IS1181 element and a 73-bp sequence were identified as inverted repeats, which could explain the possible mechanism of the inversion event.Staphylococcus aureus is a gram-positive bacterium of immense clinical importance. This opportunistic pathogen is capable of causing a wide range of diseases from skin and soft-tissue infections to life-threatening bacteremia, endocarditis, and osteomyelitis (14). Approximately 75% of the S. aureus genome is composed of a core genome that is common in all strains, and 25% of the genome is composed of variable regions which can differ between different strains (4, 16, 24-26). S. aureus evolves primarily by introducing single nucleotide polymorphisms in its core genome and by acquiring mobile genetic elements (MGEs) through horizontal gene transfer. These MGEs include pathogenicity/genomic islands, plasmids, transposons, and bacteriophages that become integrated in the chromosome (4, 11, 16, 31, 32). Despite being a heterogeneous organism, genetic recombination in S. aureus is proposed to be rather rare (20, 24, 29, 35). Its clones are more likely to evolve by point mutations than by recombination events (12). The MGEs contribute to the phenotypic and genotypic diversity seen with the S. aureus population. Acquisition of the staphylococcal cassette chromosome (SCCmec) elements through site-specific recombinases has led to the epidemic of methicillin-resistant S. aureus (MRSA) strains in hospitals and communities all over the world (6, 10, 15). In recent years, the integration of arginine catabolite mobile element in the USA300 lineage of MRSA has been proposed to give the pathogen its epidemiological advantage, including traits for surviving in low-pH conditions and oxygen tension environments (11). In addition, chromosomal replacements have been observed within lineages of sequence type 34 (ST34) and ST42 (34) and ST8 and ST30 (13).Genomic rearrangements, such as inversions, have been observed with genomes of enteric bacteria, such as Salmonella enterica, Shigella flexneri, Yersinia pestis KIM, Escherichia coli (K12 and O157:H7), and group A Streptococcus pyogenes (8, 9, 18, 27, 28, 30, 37). No genomic inversions in S. aureus have been reported to date. With the use of optical mapping, large genomic rearrangements, such as inversions, that would otherwise be missed with other comparative genotyping approaches, including microarray analysis, can be identified. Optical mapping uses high-resolution restriction maps (optical maps) of a bacterial genome to determine its genomic organization (5, 21, 23, 33, 36). These optical maps can be compared to an in silico (virtual) restriction map of a known genome sequence and can be used to identify gene rearrangements and their locations. Using optical mapping in conjunction with subsequent site-specific PCR and sequencing, we report the identification, approximate location, and partial characterization of an ∼500-kb DNA element in NRS387, a USA800 strain that was found to be inverted relative to USA300FPR3757. Identification of IS1181 elements and a novel 73-bp element at both ends of the ∼500-kb element in NRS387 could suggest the possibility of an inversion event in an ancestral strain of NRS387.     

Complete Genome Sequence of Staphylococcus aureus Strain JKD6159, a Unique Australian Clone of ST93-IV Community Methicillin-Resistant Staphylococcus aureus

Community methicillin-resistant Staphylococcus aureus (cMRSA) is an emerging issue that has resulted in multiple worldwide epidemics. We report the first complete genome sequence of an ST93-MRSA-IV clinical isolate that caused severe invasive infection and a familial outbreak of skin infection. This isolate is a representative of the most common Australian clone of cMRSA that is more distantly related to the previously sequenced genomes of S. aureus.Staphylococcus aureus is a major cause of both hospital- and community-acquired infections, with rapid emergence of antibiotic resistance, in particular methicillin resistance, adding complexity to the treatment of this organism (3). While previously a hospital problem, methicillin-resistant S. aureus (MRSA) is now being increasingly documented in healthy patients in the community, and these isolates are termed “community MRSA” (cMRSA). A number of cMRSA genomes have been sequenced; however, these are phylogenetically closely related to each other. In contrast, ST93-MRSA-IV, a unique Australian clone, is a singleton by multilocus sequence typing (MLST) eBURST analysis (4). It is now the dominant cMRSA clone in Australia and is associated with both skin infection and severe invasive infection, including necrotizing pneumonia, deep-seated abscesses, and septicemia (5, 10). JKD6159 is a representative ST93-MRSA-IV clinical isolate which caused septicemia and multifocal pulmonary and musculoskeletal abscesses in a previously well intravenous drug user and also resulted in a familial outbreak of skin infection.The genome sequence of S. aureus strain JKD6159 was determined by high-throughput whole-genome shotgun sequencing, using both Illumina GAII (Illumina, CA) and Roche GS FLX Titanium (Roche Diagnostics, Basel, Switzerland) sequencing technologies, producing approximately 164× and 32× coverage of the genome, respectively. The GS FLX Titanium reads were assembled using Newbler 2.0.01.12, resulting in 56 contigs totaling 2.8 Mbp (9). The paired GAII reads were aligned to the contigs using SHRiMP 1.3.2 to identify and correct 74 homopolymeric sequencing errors (11). Optical mapping was used to produce a high-resolution XbaI chromosome restriction map, and the contigs were ordered and oriented against this map using MapSolver 2.1.1 (Opgen). Gap closures were performed by PCR followed by Sanger sequencing of amplification products (3730S genetic analyzer sequencer; Applied Biosystems, CA). The finished sequence was validated by reference to the XbaI optical map, Roche GS FLX Titanium mate pair analysis, and Illumina paired-end-read analysis.Protein coding regions were predicted using GeneMarkS 4.6b, tRNA genes using tRNAscan-SE 1.23, and rRNA genes using RNAmmer 1.2 (2, 7, 8). Gene products were assigned using HMMER 3.0 against the Pfam database (release 23) and BLAST 2.2.23 against RefSeq Proteins (April 2010) and the Conserved Domain Database (v2.22) (1, 6). These automated analyses were followed by manual curation, including comparison with other completed S. aureus genomes.The genome of S. aureus strain JKD6159 consists of a circular 2,811,435-bp chromosome with 33% G+C content—similar to those of other staphylococci—and one circular plasmid of 20,730 bp. A total of 2,605 coding regions, 57 tRNA genes, and 5 rRNA loci were detected. Over 67% of genes were assigned to specific Clusters of Orthologous Groups (COG) Database functional groups, and 40% were assigned an enzyme classification number (12).Initial analysis of the whole-genome sequence of JKD6159 confirms that ST93-MRSA-IV is distantly related to other previously sequenced S. aureus genomes. ST93-MRSA-IV has a distinct accessory genome. There were a number of regions of difference in JKD6159 that contain coding sequences (CDS) not present in any other published S. aureus genomes. Additionally, the ssl gene cluster in JKD6159 appears distinct from other sequenced S. aureus isolates. Comparison with other S. aureus genomes also shows that although JKD6159 carries lukSF-PV (the genes encoding Panton-Valentine leukocidin), there is a relative paucity of virulence factors such as tst-1, genes encoding staphylococcal enterotoxins A to U, and the ACME locus. Further analysis of the genome is now under way to identify factors that might explain the emergence of this MRSA strain in the community.