Quality Control Test for Sequence-Phenotype Assignments

Relating a gene mutation to a phenotype is a common task in different disciplines such as protein biochemistry. In this endeavour, it is common to find false relationships arising from mutations introduced by cells that may be depurated using a phenotypic assay; yet, such phenotypic assays may introduce additional false relationships arising from experimental errors. Here we introduce the use of high-throughput DNA sequencers and statistical analysis aimed to identify incorrect DNA sequence-phenotype assignments and observed that 10–20% of these false assignments are expected in large screenings aimed to identify critical residues for protein function. We further show that this level of incorrect DNA sequence-phenotype assignments may significantly alter our understanding about the structure-function relationship of proteins. We have made available an implementation of our method at http://bis.ifc.unam.mx/en/software/chispas.     

Operon prediction in Pyrococcus furiosus

Identification of operons in the hyperthermophilic archaeon Pyrococcus furiosus represents an important step to understanding the regulatory mechanisms that enable the organism to adapt and thrive in extreme environments. We have predicted operons in P.furiosus by combining the results from three existing algorithms using a neural network (NN). These algorithms use intergenic distances, phylogenetic profiles, functional categories and gene-order conservation in their operon prediction. Our method takes as inputs the confidence scores of the three programs, and outputs a prediction of whether adjacent genes on the same strand belong to the same operon. In addition, we have applied Gene Ontology (GO) and KEGG pathway information to improve the accuracy of our algorithm. The parameters of this NN predictor are trained on a subset of all experimentally verified operon gene pairs of Bacillus subtilis. It subsequently achieved 86.5% prediction accuracy when applied to a subset of gene pairs for Escherichia coli, which is substantially better than any of the three prediction programs. Using this new algorithm, we predicted 470 operons in the P.furiosus genome. Of these, 349 were validated using DNA microarray data.     

Metabolic Engineering of Carotenoid Biosynthesis in Escherichia coli by Ordered Gene Assembly in Bacillus subtilis

We attempted to optimize the production of zeaxanthin in Escherichia coli by reordering five biosynthetic genes in the natural carotenoid cluster of Pantoea ananatis. Newly designed operons for zeaxanthin production were constructed by the ordered gene assembly in Bacillus subtilis (OGAB) method, which can assemble multiple genes in one step using an intrinsic B. subtilis plasmid transformation system. The highest level of production of zeaxanthin in E. coli (820 μg/g [dry weight]) was observed in the transformant with a plasmid in which the gene order corresponds to the order of the zeaxanthin metabolic pathway (crtE-crtB-crtI-crtY-crtZ), among a series of plasmids with circularly permuted gene orders. Although two of five operons using intrinsic zeaxanthin promoters failed to assemble in B. subtilis, the full set of operons was obtained by repressing operon expression during OGAB assembly with a p_R promoter-cI repressor system. This result suggests that repressing the expression of foreign genes in B. subtilis is important for their assembly by the OGAB method. For all tested operons, the abundance of mRNA decreased monotonically with the increasing distance of the gene from the promoter in E. coli, and this may influence the yield of zeaxanthin. Our results suggest that rearrangement of biosynthetic genes in the order of the metabolic pathway by the OGAB method could be a useful approach for metabolic engineering.     

The life-cycle of operons

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution may be driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysA and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although genes within operons are usually closely spaced because of a neutral bias toward deletion and because of selection against large overlaps, genes in highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution may be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.     

A fuzzy guided genetic algorithm for operon prediction

Motivation: The operon structure of the prokaryotic genome isa critical input for the reconstruction of regulatory networksat the whole genome level. As experimental methods for the detectionof operons are difficult and time-consuming, efforts are beingput into developing computational methods that can use availablebiological information to predict operons. Method: A genetic algorithm is developed to evolve a startingpopulation of putative operon maps of the genome into progressivelybetter predictions. Fuzzy scoring functions based on multiplecriteria are used for assessing the ‘fitness’ ofthe newly evolved operon maps and guiding their evolution. Results: The algorithm organizes the whole genome into operons.The fuzzy guided genetic algorithm-based approach makes it possibleto use diverse biological information like genome sequence data,functional annotations and conservation across multiple genomes,to guide the organization process. This approach does not requireany prior training with experimental operons. The predictionsfrom this algorithm for Escherchia coli K12 and Bacillus subtilisare evaluated against experimentally discovered operons forthese organisms. The accuracy of the method is evaluated usingan ROC (receiver operating characteristic) analysis. The areaunder the ROC curve is around 0.9, which indicates excellentaccuracy. Contact: roschen_csir{at}rediffmail.com     

Stability of the hybrid plasmid pIM138 and its curing by some eliminating agents

Hybrid plasmid pIM138 was constructed by insertion of a chromosomal fragment with the threonme operon fromEscherichia coli into the pBR322 vector. Molar mass of pIM138 was 2.8 Mg/mol. Heteroduplexes between pBR322 vector and pIM138 hybrid DNA molecules were prepared. The hybrid plasmid shows a high stability against the curing effect of rifampicin and clorobiocm inE. coli SK1590thr host.     

The dnd operon for DNA phosphorothioation modification system in Escherichia coli is located in diverse genomic islands

Background

Strains of Escherichia coli that are non-typeable by pulsed-field gel electrophoresis (PFGE) due to in-gel degradation can influence their molecular epidemiological data. The DNA degradation phenotype (Dnd⁺) is mediated by the dnd operon that encode enzymes catalyzing the phosphorothioation of DNA, rendering the modified DNA susceptible to oxidative cleavage during a PFGE run. In this study, a PCR assay was developed to detect the presence of the dnd operon in Dnd⁺E. coli strains and to improve their typeability. Investigations into the genetic environments of the dnd operon in various E. coli strains led to the discovery that the dnd operon is harboured in various diverse genomic islands.

Results

The dndBCDE genes (dnd operon) were detected in all Dnd⁺E. coli strains by PCR. The addition of thiourea improved the typeability of Dnd⁺E. coli strains to 100% using PFGE and the Dnd⁺ phenotype can be observed in both clonal and genetically diverse E. coli strains.Genomic analysis of 101 dnd operons from genome sequences of Enterobacteriaceae revealed that the dnd operons of the same bacterial species were generally clustered together in the phylogenetic tree. Further analysis of dnd operons of 52 E. coli genomes together with their respective immediate genetic environments revealed a total of 7 types of genetic organizations, all of which were found to be associated with genomic islands designated dnd-encoding GIs. The dnd-encoding GIs displayed mosaic structure and the genomic context of the 7 islands (with 1 representative genome from each type of genetic organization) were also highly variable, suggesting multiple recombination events. This is also the first report where two dnd operons were found within a strain although the biological implication is unknown. Surprisingly, dnd operons were frequently found in pathogenic E. coli although their link with virulence has not been explored.

Conclusion

Genomic islands likely play an important role in facilitating the horizontal gene transfer of the dnd operons in E. coli with 7 different types of islands discovered so far.

Electronic supplementary material

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

Secret-AAR: a web server to assess the antigenic density of proteins and homology search against bacterial and parasite secretome proteins

The secretome refers to all the Excreted/Secreted (ES) proteins of a cell, and these are involved in critical biological processes, such as cell-cell communication, and host immune responses. Recently, we introduced the Abundance of Antigenic Aegions (AAR) value to assess the protein antigenic density and to evaluate the antigenic potential of secretomes. Here, to facilitate the AAR calculation, we implemented it as a user-friendly webserver. We extended the webserver capabilities implementing a sequence-based tool for searching homologous proteins across secretomes, including experimental and predicted secretomes of Mycobacterium tuberculosis and Taenia solium. Additionally, twelve secretomes of helminths, five of Mycobacterium and two of Gram-negative bacteria are also available. Our webserver is a useful tool for researchers working on immunoinformatics and reverse vaccinology, aiming at discovering candidate proteins for new vaccines or diagnostic tests, and it can be used to prioritize the experimental analysis of proteins for druggability assays. The Secret-AAR web server is available at http://microbiomics.ibt.unam.mx/tools/aar/.     

Escherichia coli Enoyl-Acyl Carrier Protein Reductase (FabI) Supports Efficient Operation of a Functional Reversal of the β-Oxidation Cycle

We recently used a synthetic/bottom-up approach to establish the identity of the four enzymes composing an engineered functional reversal of the β-oxidation cycle for fuel and chemical production in Escherichia coli (J. M. Clomburg, J. E. Vick, M. D. Blankschien, M. Rodriguez-Moya, and R. Gonzalez, ACS Synth Biol 1:541–554, 2012, http://dx.doi.org/10.1021/sb3000782). While native enzymes that catalyze the first three steps of the pathway were identified, the identity of the native enzyme(s) acting as the trans-enoyl coenzyme A (CoA) reductase(s) remained unknown, limiting the amount of product that could be synthesized (e.g., 0.34 g/liter butyrate) and requiring the overexpression of a foreign enzyme (the Euglena gracilis trans-enoyl-CoA reductase [EgTER]) to achieve high titers (e.g., 3.4 g/liter butyrate). Here, we examine several native E. coli enzymes hypothesized to catalyze the reduction of enoyl-CoAs to acyl-CoAs. Our results indicate that FabI, the native enoyl-acyl carrier protein (enoyl-ACP) reductase (ENR) from type II fatty acid biosynthesis, possesses sufficient NADH-dependent TER activity to support the efficient operation of a β-oxidation reversal. Overexpression of FabI proved as effective as EgTER for the production of butyrate and longer-chain carboxylic acids. Given the essential nature of fabI, we investigated whether bacterial ENRs from other families were able to complement a fabI deletion without promiscuous reduction of crotonyl-CoA. These characteristics from Bacillus subtilis FabL enabled ΔfabI complementation experiments that conclusively established that FabI encodes a native enoyl-CoA reductase activity that supports the β-oxidation reversal in E. coli.     

Non-contiguous finished genome sequence and description of Bartonella florenciae sp. nov.

Bartonella florenciae sp. nov. strain R4^T is the type strain of B. florenciae sp. nov., a new species within the genus Bartonella. This strain, whose genome is described here, was isolated in France from the spleen of the shrew Crocidura russula. B. florenciae is an aerobic, rod-shaped, Gram-negative bacterium. Here we describe the features of this organism, together with the complete genome sequence and its annotation. The 2,010,844 bp-long genome contains 1,909 protein-coding and 46 RNA genes, including two rRNA operons.     

Engineered ribosomal RNA operon copy-number variants of E. coli reveal the evolutionary trade-offs shaping rRNA operon number

Ribosomal RNA (rrn) operons, characteristically present in several copies in bacterial genomes (7 in E. coli), play a central role in cellular physiology. We investigated the factors determining the optimal number of rrn operons in E. coli by constructing isogenic variants with 5–10 operons. We found that the total RNA and protein content, as well as the size of the cells reflected the number of rrn operons. While growth parameters showed only minor differences, competition experiments revealed a clear pattern: 7–8 copies were optimal under conditions of fluctuating, occasionally rich nutrient influx and lower numbers were favored in stable, nutrient-limited environments. We found that the advantages of quick adjustment to nutrient availability, rapid growth and economic regulation of ribosome number all contribute to the selection of the optimal rrn operon number. Our results suggest that the wt rrn operon number of E. coli reflects the natural, ‘feast and famine’ life-style of the bacterium, however, different copy numbers might be beneficial under different environmental conditions. Understanding the impact of the copy number of rrn operons on the fitness of the cell is an important step towards the creation of functional and robust genomes, the ultimate goal of synthetic biology.     

Genome sequence and description of Nesterenkonia massiliensis sp. nov. strain NP1T

Nesterenkonia massiliensis sp. nov., strain NP1^T, is the type strain of Nesterenkonia massiliensis sp. nov., a new species within the genus Nesterenkonia. This strain, whose genome is described here, was isolated from the feces of a 32-year-old French woman suffering from AIDS and living in Marseille. Nesterenkonia massiliensis is a Gram-positive aerobic coccus. Here, we describe the features of this bacterium, together with the complete genome sequencing and annotation. The 2,726,371 bp long genome (one chromosome but no plasmid) contains 2,663 protein-coding and 51 RNA genes, including 1 rRNA operon.     

Non-contiguous finished genome sequence and description of Oceanobacillus massiliensis sp. nov.

Oceanobacillus massiliensis strain N’Diop^T sp. nov. is the type strain of O. massiliensis sp. nov., a new species within the genus Oceanobacillus. This strain, whose genome is described here, was isolated from the fecal flora of a healthy patient. O. massiliensis is an aerobic rod. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 3,532,675 bp long genome contains 3,519 protein-coding genes and 72 RNA genes, including between 6 and 8 rRNA operons.