Two genome sequences of the same bacterial strain, Gluconacetobacter diazotrophicus PAl 5, suggest a new standard in genome sequence submission

A deoxyviolacein-based transposon insertion vector for pigmented tracer studies

Pigments provide a simple means to rapidly visually ascertain the quantities or presence of specific microbes in a complex community. The selection of pigment-producing colonies that are simple to differentiate from common colony phenotypes provides a high degree of certainty for the identity of pigment-tagged strains. Successful employment of pigment production is dependent on various intrinsic factors related to proper levels of gene expression and pigment production that are not always easy to predict and vary within each microbe. We have constructed a simple transposon system that incorporates the genes for the production of deoxyviolacein, a pigment produced from intracellular reserves of the amino acid tryptophan, to randomly insert these genes throughout the genome. This tool allows the user to select from many thousands of potential sites throughout a bacterial genome for an ideal location to generate the desired amount of pigment. We have applied this system to a small selection of endophytes and other model bacteria to differentiate these strains from complex communities and confirm their presence after several weeks in natural environments. We provide two examples of applications using the pigments to trace strains following introduction into plant tissues or to produce a reporter strain for extracellular nitrogen compound sensing. We recognize that this tool could have far broader utility in other applications and microbes, and describe the methodology for use by the greater scientific community.

Enhanced extracellular ammonium release in the plant endophyte Gluconacetobacter diazotrophicus through genome editing

IMPORTANCE

Our results demonstrate increased extracellular ammonium release in the endophyte plant growth-promoting bacterium Gluconacetobacter diazotrophicus. Strains were constructed in a manner that leaves no antibiotic markers behind, such that these strains contain no transgenes. Levels of ammonium achieved by cultures of modified G. diazotrophicus strains reached concentrations of approximately 18 mM ammonium, while wild-type G. diazotrophicus remained much lower (below 50 µM). These findings demonstrate a strong potential for further improving the biofertilizer potential of this important microbe.

Gluconacetobacter diazotrophicus Gene Fitness during Diazotrophic Growth

Plant growth-promoting (PGP) bacteria are important to the development of sustainable agricultural systems. PGP microbes that fix atmospheric nitrogen (diazotrophs) could minimize the application of industrially derived fertilizers and function as a biofertilizer. The bacterium Gluconacetobacter diazotrophicus is a nitrogen-fixing PGP microbe originally discovered in association with sugarcane plants, where it functions as an endophyte. It also forms endophyte associations with a range of other agriculturally relevant crop plants. G. diazotrophicus requires microaerobic conditions for diazotrophic growth. We generated a transposon library for G. diazotrophicus and cultured the library under various growth conditions and culture medium compositions to measure fitness defects associated with individual transposon inserts (transposon insertion sequencing [Tn-seq]). Using this library, we probed more than 3,200 genes and ascertained the importance of various genes for diazotrophic growth of this microaerobic endophyte. We also identified a set of essential genes. IMPORTANCE Our results demonstrate a succinct set of genes involved in diazotrophic growth for G. diazotrophicus, with a lower degree of redundancy than what is found in other model diazotrophs. The results will serve as a valuable resource for those interested in biological nitrogen fixation and will establish a baseline data set for plant free growth, which could complement future studies related to the endophyte relationship.

Dissection and reconstitution provide insights into electron transport in the membrane-bound aldehyde dehydrogenase complex of Gluconacetobacter diazotrophicus

Acetic acid bacteria catalyze the two-step oxidation of ethanol to acetic acid using the membrane-bound enzymes pyrroloquinoline quinone-dependent alcohol dehydrogenase and molybdopterin-dependent aldehyde dehydrogenase (ALDH). Although the reducing equivalents from the substrate are transferred to ubiquinone in the membrane, intramolecular electron transport in ALDH is not understood. Here, we purified the AldFGH complex, the membrane-bound ALDH that is physiologically relevant to acetic acid fermentation in Gluconacetobacter diazotrophicus strain PAL5. The purified AldFGH complex showed acetaldehyde:ubiquinone (Q 2 ) oxidoreductase activity. C -type cytochromes of the AldFGH complex (in the AldF subunit) were reduced by acetaldehyde. Then, we genetically dissected the AldFGH complex into AldGH and AldF units and reconstituted them. The AldGH subcomplex showed acetaldehyde:ferricyanide oxidoreductase activity, but not Q 2 reductase activity. The ALDH activity of AldGH was not found in membranes but in the soluble fraction of the recombinant strain, suggesting that the AldF subunit is responsible for membrane binding of the AldFGH complex. AldFGH complex reconstituted from the AldGH subcomplex and AldF showed Q 2 reductase activity. Absorption spectra of the purified AldGH subcomplex suggested the presence of an [Fe–S] cluster, which can be reduced by acetaldehyde. We propose a model in which electrons from the substrate are abstracted by a molybdopterin in the AldH subunit and transferred to [Fe–S] cluster(s) in the AldG subunit, followed by electron transport to c -type cytochrome centers in the AldF subunit, which is the site of ubiquinone reduction in the membrane.

Importance

Two membrane-bound enzymes of acetic acid bacteria—pyrroloquinoline quinone-dependent alcohol dehydrogenase and molybdopterin-dependent aldehyde dehydrogenase (ALDH)—are responsible for vinegar production. Upon oxidation of acetaldehyde, ALDH reduces ubiquinone in the cytoplasmic membrane. ALDH is an enzyme complex of three subunits. Here, we tried to understand how ALDH works by using a classical biochemical approach and genetic engineering to dissect the enzyme complex into soluble and membrane-bound parts. The soluble part had limited activity in vitro , and did not reduce ubiquinone. However, enzyme complex reconstituted from the soluble and membrane-bound parts showed ubiquinone reduction activity. The proposed working model of ALDH provides a better understanding of how the enzyme works in the vinegar fermentation process.

Whole genome analysis of Gluconacetobacter azotocaptans DS1 and its beneficial effects on plant growth

Plant-associated bacteria play an important role in the enhancement of plant growth and productivity. Gluconacetobacter azotocaptans is an exceptional bacterium considering that till today it has been isolated and reported only from Mexico and Canada. It is a plant growth-promoting bacterium and can be used as biofertilizer for different crops and vegetables. The objective of the current study was to evaluate the inoculation effect of Gluconacetobacter azotocaptans DS1, Pseudomonas putida CQ179, Azosprillium zeae N7, Azosprillium brasilense N8, and Azosprillium canadense DS2, on the growth of vegetables including cucumber, sweet pepper, radish, and tomato. All strains increased the vegetables' growth; however, G. azotocaptans DS1 showed better results as compared to other inoculated and control plants and significantly increased the plant biomass of all vegetables. Therefore, the whole genome sequence of G. azotocaptans DS1 was analyzed to predict genes involved in plant growth promotion, secondary metabolism, antibiotics resistance, and bioremediation of heavy metals. Results of genome analysis revealed that G. azotocaptans DS1 has a circular chromosome with a size of 4.3 Mbp and total 3898 protein-coding sequences. Based on functional analysis, genes for nitrogen fixation, phosphate solubilization, indole acetic acid, phenazine, siderophore production, antibiotic resistance, and bioremediation of heavy metals including copper, zinc, cobalt, and cadmium were identified. Collectively, our findings indicated that G. azotocaptans DS1 can be used as a biofertilizer and biocontrol agent for growth enhancement of different crops and vegetables.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02996-1.

Elimination of aromatic fusel alcohols as by-products of Saccharomyces cerevisiae strains engineered for phenylpropanoid production by 2-oxo-acid decarboxylase replacement

Engineered strains of the yeast Saccharomyces cerevisiae are intensively studied as production platforms for aromatic compounds such as hydroxycinnamic acids, stilbenoids and flavonoids. Heterologous pathways for production of these compounds use l-phenylalanine and/or l-tyrosine, generated by the yeast shikimate pathway, as aromatic precursors. The Ehrlich pathway converts these precursors to aromatic fusel alcohols and acids, which are undesirable by-products of yeast strains engineered for production of high-value aromatic compounds. Activity of the Ehrlich pathway requires any of four S. cerevisiae 2-oxo-acid decarboxylases (2-OADCs): Aro10 or the pyruvate-decarboxylase isoenzymes Pdc1, Pdc5, and Pdc6. Elimination of pyruvate-decarboxylase activity from S. cerevisiae is not straightforward as it plays a key role in cytosolic acetyl-CoA biosynthesis during growth on glucose. In a search for pyruvate decarboxylases that do not decarboxylate aromatic 2-oxo acids, eleven yeast and bacterial 2-OADC-encoding genes were investigated. Homologs from Kluyveromyces lactis (KlPDC1), Kluyveromyces marxianus (KmPDC1), Yarrowia lipolytica (YlPDC1), Zymomonas mobilis (Zmpdc1) and Gluconacetobacter diazotrophicus (Gdpdc1.2 and Gdpdc1.3) complemented a Pdc⁻ strain of S. cerevisiae for growth on glucose. Enzyme-activity assays in cell extracts showed that these genes encoded active pyruvate decarboxylases with different substrate specificities. In these in vitro assays, ZmPdc1, GdPdc1.2 or GdPdc1.3 had no substrate specificity towards phenylpyruvate. Replacing Aro10 and Pdc1,5,6 by these bacterial decarboxylases completely eliminated aromatic fusel-alcohol production in glucose-grown batch cultures of an engineered coumaric acid-producing S. cerevisiae strain. These results outline a strategy to prevent formation of an important class of by-products in ‘chassis’ yeast strains for production of non-native aromatic compounds.

•

Identification of pyruvate decarboxylases active with pyruvate but not with aromatic 2-oxo acids.

•

Zymomonas mobilis pyruvate decarboxylase can replace the native yeast enzymes.

•

Expression of Z. mobilis pyruvate decarboxylase removes formation of fusel alcohols.

•

Elimination of fusel alcohol by products improves formation of coumaric acid.

•

Decarboxylase swapping is a beneficial strategy for production of non-native aromatics.

Relocation of dehydroquinate dehydratase to the periplasmic space improves dehydroshikimate production with Gluconobacter oxydans strain NBRC3244

3-Dehydroshikimate (3-DHS) is a key intermediate for the synthesis of various compounds, including the antiviral drug oseltamivir. The Gluconobacter oxydans strain NBRC3244 intrinsically oxidizes quinate to produce 3-dehydroquinate (3-DHQ) in the periplasmic space. Even though a considerable activity is detected in the recombinant G. oxydans homologously overexpressing type II dehydroquinate dehydratase (DHQase) encoded in the aroQ gene at a pH where it grows, an alkaline shift of the culture medium is required for 3-DHS production in the middle of cultivation. Here, we attempted to adopt type I DHQase encoded in the aroD gene of Gluconacetobacter diazotrophicus strain PAL5 because the type I DHQase works optimally at weak acid, which is preferable for growth conditions of G. oxydans. In addition, we anticipated that subcellular localization of DHQase is the cytoplasm, and therefore, transports of 3-DHQ and 3-DHS across the cytoplasmic membrane are rate-limiting steps in the biotransformation. The Sec- and TAT-dependent signal sequences for secretion were attached to the N terminus of AroD to change the subcellular localization. G. oxydans that expresses the TAT-AroD derivative achieved 3-DHS production at a tenfold higher rate than the reference strain that expresses wild-type AroD even devoid of alkaline shift. Enzyme activity with the intact cell suspension and signal sequence cleavage supported the relocation of AroD to the periplasmic space. The present study suggests that the relocation of DHQase improves 3-DHS production in G. oxydans and represents a proof of concept for the potential of enzyme relocation in metabolic engineering. KEY POINTS: • Type-I dehydroquinate dehydratase (DHQase) was expressed in Gluconobacter oxydans. • Cytoplasmic DHQase was relocated to the periplasmic space in G. oxydans. • Relocation of DHQase in G. oxydans improved 3-dehydroshikimate production.

Phylogenomic Framework for Taxonomic Delineation of Paracoccus spp. and Exploration of Core-Pan Genome

The taxonomic classification of metabolically versatile Paracoccus spp. has been so far performed using polyphasic approach. The topology of single gene phylogenies, however, has highlighted ambiguous species assignments. In the present study, genome based multi-gene phylogenies and overall genome related index were used for species threshold assessment. Comprehensive phylogenomic analysis of Paracoccus genomes (n = 103) showed concordant clustering of strains across multi-gene marker set phylogenies (nMC = 0.08-0.14); as compared to 16S rDNA phylogeny (nMC = 0.37-0.42) suggesting robustness of multi gene phylogenies in drawing phylogenetic inferences. Functional gene content distribution across the genus showed that only 1.7% gene content constitutes the core genome highlighting the significance of extensive genomic variability in the evolution of Paracoccus spp. Further, genome metrics were used to validate characterized strains, identifying classification anomalies (n = 13), and based on this, genome derived taxonomic amendments were notified in present study. Conclusively, validated metric tools can be employed on whole genome sequences, including draft assemblies, for the assessment and assignment of uncharacterized strains and species level ascription of newly isolated Paracoccus strains in future.

Major aldehyde dehydrogenase AldFGH of Gluconacetobacter diazotrophicus is independent of pyrroloquinoline quinone but dependent on molybdopterin for acetic acid fermentation

Acetic acid fermentation involves the oxidation of ethanol to acetic acid via acetaldehyde as the intermediate and is catalyzed by the membrane-bound alcohol dehydrogenase (ADH) and aldehyde dehydrogenase (ALDH) of acetic acid bacteria. Although ADH depends on pyrroloquinoline quinone (PQQ), the prosthetic group associated with ALDH remains a matter of debate. This study aimed to address the dependency of ALDH of Gluconacetobacter diazotrophicus strain PAL5 on PQQ and the physiological role of ALDH in acetic acid fermentation. We constructed deletion mutant strains for both the ALDH gene clusters of PAL5, aldFGH and aldSLC. In addition, the adhAB operon for ADH was eliminated, since it shows ALDH activity. The triple-deletion derivative ΔaldFGH ΔaldSLC ΔadhAB failed to show ALDH activity, which suggested that ALDH activity in PAL5 is derived from these three enzyme complexes. Since the single-gene cluster deletion derivative ΔaldFGH lost most ALDH activity, and accumulated much higher acetaldehyde than wild type under acetic acid fermentation conditions, we concluded that AldFGH functions as the major ALDH in PAL5. Furthermore, deletion of the PQQ biosynthesis gene cluster (pqqABCDE) abolished ADH activity completely, but did not affect ALDH activity. Instead, the molybdopterin biosynthesis gene deletion derivatives lost ALDH activity. Thus, we concluded that the AldFGH and AldSLC complexes of Ga. diazotrophicus PAL5 require a form of molybdopterin but not PQQ for ALDH activity. KEY POINTS: • AldFGH is the major aldehyde dehydrogenase in Gluconacetobacter diazotrophicus PAL5. • Acetaldehyde accumulated from ethanol in the absence of AldFGH. • Molybdopterin, rather than pyrroloquinoline quinone, is required for AldFGH.

Molecular Cloning and Exploration of the Biochemical and Functional Analysis of Recombinant Glucose-6-Phosphate Dehydrogenase from Gluconoacetobacter diazotrophicus PAL5

Gluconacetobacter diazotrophicus PAL5 (GDI) is an endophytic bacterium with potential biotechnological applications in industry and agronomy. The recent description of its complete genome and its principal metabolic enzymes suggests that glucose metabolism is accomplished through the pentose phosphate pathway (PPP); however, the enzymes participating in this pathway have not yet been characterized in detail. The objective of the present work was to clone, purify, and biochemically and physicochemically characterize glucose-6-phosphate dehydrogenase (G6PD) from GDI. The gene was cloned and expressed as a tagged protein in E. coli to be purified by affinity chromatography. The native state of the G6PD protein in the solution was found to be a tetramer with optimal activity at pH 8.8 and a temperature between 37 and 50 °C. The apparent Km values for G6P and nicotinamide adenine dinucleotide phosphate (NADP⁺) were 63 and 7.2 μM, respectively. Finally, from the amino acid sequence a three-dimensional (3D) model was obtained, which allowed the arrangement of the amino acids involved in the catalytic activity, which are conserved (RIDHYLGKE, GxGGDLT, and EKPxG) with those of other species, to be identified. This characterization of the enzyme could help to identify new environmental conditions for the knowledge of the plant–microorganism interactions and a better use of GDI in new technological applications.

Genomic analysis of bacteria in the Acute Oak Decline pathobiome

The UK’s native oak is under serious threat from Acute Oak Decline (AOD). Stem tissue necrosis is a primary symptom of AOD and several bacteria are associated with necrotic lesions. Two members of the lesion pathobiome, Brenneria goodwinii and Gibbsiella quercinecans, have been identified as causative agents of tissue necrosis. However, additional bacteria including Lonsdalea britannica and Rahnella species have been detected in the lesion microbiome, but their role in tissue degradation is unclear. Consequently, information on potential genome-encoded mechanisms for tissue necrosis is critical to understand the role and mechanisms used by bacterial members of the lesion pathobiome in the aetiology of AOD. Here, the whole genomes of bacteria isolated from AOD-affected trees were sequenced, annotated and compared against canonical bacterial phytopathogens and non-pathogenic symbionts. Using orthologous gene inference methods, shared virulence genes that retain the same function were identified. Furthermore, functional annotation of phytopathogenic virulence genes demonstrated that all studied members of the AOD lesion microbiota possessed genes associated with phytopathogens. However, the genome of B. goodwinii was the most characteristic of a necrogenic phytopathogen, corroborating previous pathological and metatranscriptomic studies that implicate it as the key causal agent of AOD lesions. Furthermore, we investigated the genome sequences of other AOD lesion microbiota to understand the potential ability of microbes to cause disease or contribute to pathogenic potential of organisms isolated from this complex pathobiome. The role of these members remains uncertain but some such as G. quercinecans may contribute to tissue necrosis through the release of necrotizing enzymes and may help more dangerous pathogens activate and realize their pathogenic potential or they may contribute as secondary/opportunistic pathogens with the potential to act as accessory species for B. goodwinii. We demonstrate that in combination with ecological data, whole genome sequencing provides key insights into the pathogenic potential of bacterial species whether they be phytopathogens, part-contributors or stimulators of the pathobiome.

Conformationally Gated Electron Transfer in Nitrogenase. Isolation, Purification, and Characterization of Nitrogenase From Gluconacetobacter diazotrophicus

Nitrogenase is a complex, bacterial enzyme that catalyzes the ATP-dependent reduction of dinitrogen (N₂) to ammonia (NH₃). In its most prevalent form, it consists of two proteins, the catalytic molybdenum-iron protein (MoFeP) and its specific reductase, the iron protein (FeP). A defining feature of nitrogenase is that electron and proton transfer processes linked to substrate reduction are synchronized by conformational changes driven by ATP-dependent FeP-MoFeP interactions. Yet, despite extensive crystallographic, spectroscopic, and biochemical information on nitrogenase, the structural basis of the ATP-dependent synchronization mechanism is not understood in detail. In this chapter, we summarize some of our efforts toward obtaining such an understanding. Experimental investigations of the structure-function relationships in nitrogenase are challenged by the fact that it cannot be readily expressed heterologously in nondiazotrophic bacteria, and the purification protocols for nitrogenase are only known for a small number of diazotrophic organisms. Here, we present methods for purifying and characterizing nitrogenase from a new model organism, Gluconacetobacter diazotrophicus. We also describe procedures for observing redox-dependent conformational changes in G. diazotrophicus nitrogenase by X-ray crystallography and electron paramagnetic resonance spectroscopy, which have provided new insights into the redox-dependent conformational gating processes in nitrogenase.

Characterization of newly identified DnaA and DnaB proteins from Acetobacter

Although chromosomal replication is an essential feature of the bacterial life cycle, the replication mechanism and involved molecular players have never been properly characterized in the Acetobacter genera. Thanks to whole-genome sequencing, the unknown replication proteins from Acetobacter pasteurianus and Acetobacter orleanensis, DnaA-like and DnaB-like, could be identified. Despite the low nucleotide or amino acid similarity to the respective orthologs from Escherichia coli, their involvement during replication regulation was corroborated by artificial microRNA. In the Acetobacter genome, a novel replication origin, oriAo, was detected with three 9-nucleotide-long DnaA boxes to which DnaA-like proteins bind actively. Bacterial two-hybrid systems and co-immunoprecipitation confirmed the homologous and heterologous interactions between DnaA-like and DnaB-like proteins with their E. coli orthologs. This communication is due to the conserved tryptophan at position 6 for E. coli or 25 for Acetobacter that unables DnaA-like proteins to form oligomeric protein structures after its substitution. Altogether, these results provide novel insights into the genome replication mechanism in Acetobacter.

Isolation and characterization of active promoters from Gluconacetobacter diazotrophicus strain PAL5 using a promoter-trapping plasmid

Gluconacetobacter diazotrophicus is a nitrogen-fixing, endophytic bacterium that has the potential to promote plant growth and increase yield. Genetically modified strains might get more benefits to host plants, including through expression of useful proteins, such as Cry toxins from B. thuringiensis, or enzymes involved in phytohormone production, proteins with antagonistic activity for phytopathogens, or that improve nutrient utilization by the plant. For that, expression systems for G. diazotrophicus are needed, which requires active promoters fused to foreign (or innate) genes. This article describes the construction of a G. diazotrophicus PAL5 promoter library using a promoter-less lacZ-bearing vector, and the identification of six active promoters through β-galactosidase activity assays, sequencing and localization in the bacterial genome. The characterized promoters, which are located on distinct regions of the bacterial genome and encoding either sense or antisense transcripts, present variable expression strengths and might be used in the future for expressing useful proteins.

De novo assembly of Dekkera bruxellensis: a multi technology approach using short and long-read sequencing and optical mapping

BACKGROUND

It remains a challenge to perform de novo assembly using next-generation sequencing (NGS). Despite the availability of multiple sequencing technologies and tools (e.g., assemblers) it is still difficult to assemble new genomes at chromosome resolution (i.e., one sequence per chromosome). Obtaining high quality draft assemblies is extremely important in the case of yeast genomes to better characterise major events in their evolutionary history. The aim of this work is two-fold: on the one hand we want to show how combining different and somewhat complementary technologies is key to improving assembly quality and correctness, and on the other hand we present a de novo assembly pipeline we believe to be beneficial to core facility bioinformaticians. To demonstrate both the effectiveness of combining technologies and the simplicity of the pipeline, here we present the results obtained using the Dekkera bruxellensis genome.

METHODS

In this work we used short-read Illumina data and long-read PacBio data combined with the extreme long-range information from OpGen optical maps in the task of de novo genome assembly and finishing. Moreover, we developed NouGAT, a semi-automated pipeline for read-preprocessing, de novo assembly and assembly evaluation, which was instrumental for this work.

RESULTS

We obtained a high quality draft assembly of a yeast genome, resolved on a chromosomal level. Furthermore, this assembly was corrected for mis-assembly errors as demonstrated by resolving a large collapsed repeat and by receiving higher scores by assembly evaluation tools. With the inclusion of PacBio data we were able to fill about 5 % of the optical mapped genome not covered by the Illumina data.

CVTree3 Web Server for Whole-genome-based and Alignment-free Prokaryotic Phylogeny and Taxonomy

A faithful phylogeny and an objective taxonomy for prokaryotes should agree with each other and ultimately follow the genome data. With the number of sequenced genomes reaching tens of thousands, both tree inference and detailed comparison with taxonomy are great challenges. We now provide one solution in the latest Release 3.0 of the alignment-free and whole-genome-based web server CVTree3. The server resides in a cluster of 64 cores and is equipped with an interactive, collapsible, and expandable tree display. It is capable of comparing the tree branching order with prokaryotic classification at all taxonomic ranks from domains down to species and strains. CVTree3 allows for inquiry by taxon names and trial on lineage modifications. In addition, it reports a summary of monophyletic and non-monophyletic taxa at all ranks as well as produces print-quality subtree figures. After giving an overview of retrospective verification of the CVTree approach, the power of the new server is described for the mega-classification of prokaryotes and determination of taxonomic placement of some newly-sequenced genomes. A few discrepancies between CVTree and 16S rRNA analyses are also summarized with regard to possible taxonomic revisions. CVTree3 is freely accessible to all users at http://tlife.fudan.edu.cn/cvtree3/ without login requirements.

Complete genome and gene expression analyses of Asaia bogorensis reveal unique responses to culture with mammalian cells as a potential opportunistic human pathogen

Asaia bogorensis, a member of acetic acid bacteria (AAB), is an aerobic bacterium isolated from flowers and fruits, as well as an opportunistic pathogen that causes human peritonitis and bacteraemia. Here, we determined the complete genomic sequence of the As. bogorensis type strain NBRC 16594, and conducted comparative analyses of gene expression under different conditions of co-culture with mammalian cells and standard AAB culture. The genome of As. bogorensis contained 2,758 protein-coding genes within a circular chromosome of 3,198,265 bp. There were two complete operons encoding cytochrome bo3-type ubiquinol terminal oxidases: cyoABCD-1 and cyoABCD-2. The cyoABCD-1 operon was phylogenetically common to AAB genomes, whereas the cyoABCD-2 operon belonged to a lineage distinctive from the cyoABCD-1 operon. Interestingly, cyoABCD-1 was less expressed under co-culture conditions than under the AAB culture conditions, whereas the converse was true for cyoABCD-2. Asaia bogorensis shared pathogenesis-related genes with another pathogenic AAB, Granulibacter bethesdensis, including a gene coding pathogen-specific large bacterial adhesin and additional genes for the inhibition of oxidation and antibiotic resistance. Expression alteration of the respiratory chain and unique hypothetical genes may be key traits that enable the bacterium to survive under the co-culture conditions.

Comparative Genomics Reveal That Host-Innate Immune Responses Influence the Clinical Prevalence of Legionella pneumophila Serogroups

Legionella pneumophila is the primary etiologic agent of legionellosis, a potentially fatal respiratory illness. Amongst the sixteen described L. pneumophila serogroups, a majority of the clinical infections diagnosed using standard methods are serogroup 1 (Sg1). This high clinical prevalence of Sg1 is hypothesized to be linked to environmental specific advantages and/or to increased virulence of strains belonging to Sg1. The genetic determinants for this prevalence remain unknown primarily due to the limited genomic information available for non-Sg1 clinical strains. Through a systematic attempt to culture Legionella from patient respiratory samples, we have previously reported that 34% of all culture confirmed legionellosis cases in Ontario (n = 351) are caused by non-Sg1 Legionella. Phylogenetic analysis combining multiple-locus variable number tandem repeat analysis and sequence based typing profiles of all non-Sg1 identified that L. pneumophila clinical strains (n = 73) belonging to the two most prevalent molecular types were Sg6. We conducted whole genome sequencing of two strains representative of these sequence types and one distant neighbour. Comparative genomics of the three L. pneumophila Sg6 genomes reported here with published L. pneumophila serogroup 1 genomes identified genetic differences in the O-antigen biosynthetic cluster. Comparative optical mapping analysis between Sg6 and Sg1 further corroborated this finding. We confirmed an altered O-antigen profile of Sg6, and tested its possible effects on growth and replication in in vitro biological models and experimental murine infections. Our data indicates that while clinical Sg1 might not be better suited than Sg6 in colonizing environmental niches, increased bloodstream dissemination through resistance to the alternative pathway of complement mediated killing in the human host may explain its higher prevalence.

Enhanced de novo assembly of high throughput pyrosequencing data using whole genome mapping

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

Rapid genomic-scale analysis of Escherichia coli O104:H4 by using high-resolution alternative methods to next-generation sequencing

Two technologies, involving DNA microarray and optical mapping, were used to quickly assess gene content and genomic architecture of recent emergent Escherichia coli O104:H4 and related strains. In real-time outbreak investigations, these technologies can provide congruent perspectives on strain, serotype, and pathotype relationships. Our data demonstrated clear discrimination between clinically, temporally, and geographically distinct O104:H4 isolates and rapid characterization of strain differences.

Rapid whole genome optical mapping of Plasmodium falciparum

Background

Immune evasion and drug resistance in malaria have been linked to chromosomal recombination and gene copy number variation (CNV). These events are ideally studied using comparative genomic analyses; however in malaria these analyses are not as common or thorough as in other infectious diseases, partly due to the difficulty in sequencing and assembling complete genome drafts. Recently, whole genome optical mapping has gained wide use in support of genomic sequence assembly and comparison. Here, a rapid technique for producing whole genome optical maps of Plasmodium falciparum is described and the results of mapping four genomes are presented.

Methods

Four laboratory strains of P. falciparum were analysed using the Argus™ optical mapping system to produce ordered restriction fragment maps of all 14 chromosomes in each genome. Plasmodium falciparum DNA was isolated directly from blood culture, visualized using the Argus™ system and assembled in a manner analogous to next generation sequence assembly into maps (AssemblyViewer™, OpGen Inc.^®). Full coverage maps were generated for P. falciparum strains 3D7, FVO, D6 and C235. A reference P. falciparum in silico map was created by the digestion of the genomic sequence of P. falciparum with the restriction enzyme AflII, for comparisons to genomic optical maps. Maps were then compared using the MapSolver™ software.

Results

Genomic variation was observed among the mapped strains, as well as between the map of the reference strain and the map derived from the putative sequence of that same strain. Duplications, deletions, insertions, inversions and misassemblies of sizes ranging from 3,500 base pairs up to 78,000 base pairs were observed. Many genomic events occurred in areas of known repetitive sequence or high copy number genes, including var gene clusters and rifin complexes.

Conclusions

This technique for optical mapping of multiple malaria genomes allows for whole genome comparison of multiple strains and can assist in identifying genetic variation and sequence contig assembly. New protocols and technology allowed us to produce high quality contigs spanning four P. falciparum genomes in six weeks for less than $1,000.00 per genome. This relatively low cost and quick turnaround makes the technique valuable compared to other genomic sequencing technologies for studying genetic variation in malaria.

Whole-genome sequencing and phenotypic analysis of Bacillus subtilis mutants following evolution under conditions of relaxed selection for sporulation

Little is known about how genetic variation at the nucleotide level contributes to competitive fitness within species. During a 6,000-generation study of Bacillus subtilis evolved under relaxed selection for sporulation, a new strain, designated WN716, emerged with significantly different colony and cell morphologies; loss of sporulation, competence, acetoin production, and motility; multiple auxotrophies; and increased competitive fitness (H. Maughan and W. L. Nicholson, Appl. Environ. Microbiol. 77:4105-4118, 2011). The genome of WN716 was analyzed by OpGen optical mapping, whole-genome 454 pyrosequencing, and the CLC Genomics Workbench. No large chromosomal rearrangements were found; however, 34 single-nucleotide polymorphisms (SNPs) and +1 frameshifts were identified in WN716 that resulted in amino acid changes in coding sequences of annotated genes, and 11 SNPs were located in intergenic regions. Several classes of genes were affected, including biosynthetic pathways, sporulation, competence, and DNA repair. In several cases, attempts were made to link observed phenotypes of WN716 with the discovered mutations, with various degrees of success. For example, a +1 frameshift was identified at codon 13 of sigW, the product of which (SigW) controls a regulon of genes involved in resistance to bacteriocins and membrane-damaging antibiotics. Consistent with this finding, WN716 exhibited sensitivity to fosfomycin and to a bacteriocin produced by B. subtilis subsp. spizizenii and exhibited downregulation of SigW-dependent genes on a transcriptional microarray, consistent with WN716 carrying a knockout of sigW. The results suggest that propagation of B. subtilis for less than 2,000 generations in a nutrient-rich environment where sporulation is suppressed led to rapid initiation of genomic erosion.