High accuracy operon prediction method based on STRING database scores

We present a simple and highly accurate computational method for operon prediction, based on intergenic distances and functional relationships between the protein products of contiguous genes, as defined by STRING database (Jensen,L.J., Kuhn,M., Stark,M., Chaffron,S., Creevey,C., Muller,J., Doerks,T., Julien,P., Roth,A., Simonovic,M. et al. (2009) STRING 8–a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res., 37, D412–D416). These two parameters were used to train a neural network on a subset of experimentally characterized Escherichia coli and Bacillus subtilis operons. Our predictive model was successfully tested on the set of experimentally defined operons in E. coli and B. subtilis, with accuracies of 94.6 and 93.3%, respectively. As far as we know, these are the highest accuracies ever obtained for predicting bacterial operons. Furthermore, in order to evaluate the predictable accuracy of our model when using an organism''s data set for the training procedure, and a different organism''s data set for testing, we repeated the E. coli operon prediction analysis using a neural network trained with B. subtilis data, and a B. subtilis analysis using a neural network trained with E. coli data. Even for these cases, the accuracies reached with our method were outstandingly high, 91.5 and 93%, respectively. These results show the potential use of our method for accurately predicting the operons of any other organism. Our operon predictions for fully-sequenced genomes are available at http://operons.ibt.unam.mx/OperonPredictor/.     

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.     

Engineering EGFP reporter constructs into a 200 kb human β-globin BAC clone using GET Recombination

GET Recombination, a simple inducible homologous recombination system for Escherichia coli, was used to target insertion of an EGFP cassette between the start and termination codons of the β-globin gene in a 200 kb BAC clone. The high degree of homology between the promoter regions of the β- and δ-globin genes also allowed the simultaneous generation of a δ-globin reporter construct with the deletion of 8.8 kb of intervening sequences. Both constructs expressed EGFP after transient transfection of MEL cells. Similarly, targeting of the EGFP cassette between the promoter regions of the γ-globin genes and the termination codon of the β-globin gene enabled the generation of reporter constructs for both ^Aγ- and ^Gγ-globin genes, involving specific deletions of 24 and 29 kb of genomic sequence, respectively. Finally the EGFP cassette was also inserted between the - and β-globin genes, with the simultaneous deletion of 44 kb of intervening sequence. The modified constructs were generated at high efficiency, illustrating the usefulness of GET Recombination to generate large deletions of specific sequences in BACs for functional studies. The establishment of stable erythropoietic cell lines with these globin constructs will facilitate the search for therapeutic agents that modify the expression of the individual globin genes in a physiologically relevant manner.