Systematic Identification and Classification of ß-Lactamases Based on Sequence Similarity Criteria: ß-Lactamase Annotation.

ß-lactamases, the enzymes responsible for resistance to ß-lactam antibiotics, are widespread among prokaryotic genera. However, current ß-lactamase classification schemes do not represent their present diversity. Here, we propose a workflow to identify and classify ß-lactamases. Initially, a set of curated sequences was used as a model for the construction of profiles Hidden Markov Models (HMM), specific for each ß-lactamase class. An extensive, nonredundant set of ß-lactamase sequences was constructed from 7 different resistance proteins databases to test the methodology. The profiles HMM were improved for their specificity and sensitivity and then applied to fully assembled genomes. Five hierarchical classification levels are described, and a new class of ß-lactamases with fused domains is proposed. Our profiles HMM provide a better annotation of ß-lactamases, with classes and subclasses defined by objective criteria such as sequence similarity. This classification offers a solid base to the elaboration of studies on the diversity, dispersion, prevalence, and evolution of the different classes and subclasses of this critical enzymatic activity.

ESTs from Seeds to Assist the Selective Breeding of Jatropha curcas L. for Oil and Active Compounds.

We report here on the characterization of a cDNA library from seeds of Jatropha curcas L. at three stages of fruit maturation before yellowing. We sequenced a total of 2200 clones and obtained a set of 931 non-redundant sequences (unigenes) after trimming and quality control, ie, 140 contigs and 791 singlets with PHRED quality ≥10. We found low levels of sequence redundancy and extensive metabolic coverage by homology comparison to GO. After comparison of 5841 non-redundant ESTs from a total of 13193 reads from GenBank with KEGG, we identified tags with nucleotide variations among J. curcas accessions for genes of fatty acid, terpene, alkaloid, quinone and hormone pathways of biosynthesis. More specifically, the expression level of four genes (palmitoyl-acyl carrier protein thioesterase, 3-ketoacyl-CoA thiolase B, lysophosphatidic acid acyltransferase and geranyl pyrophosphate synthase) measured by real-time PCR proved to be significantly different between leaves and fruits. Since the nucleotide polymorphism of these tags is associated to higher level of gene expression in fruits compared to leaves, we propose this approach to speed up the search for quantitative traits in selective breeding of J. curcas. We also discuss its potential utility for the selective breeding of economically important traits in J. curcas.

Sequence analysis of NAT2 gene in Brazilians: identification of undescribed single nucleotide polymorphisms and molecular modeling of the N-acetyltransferase 2 protein structure.

N-Acetyltransferase 2 (NAT2) metabolizes a variety of xenobiotics that includes many drugs, chemicals and carcinogens. This enzyme is genetically variable in human populations and polymorphisms in the NAT2 gene have been associated with drug toxicity and efficacy as well as cancer susceptibility. Here, we have focused on the identification of NAT2 variants in Brazilian individuals from two different regions, Rio de Janeiro and Goiás, by direct sequencing, and on the characterization of new haplotypes after cloning and re-sequencing. Upon analysis of DNA samples from 404 individuals, six new SNPs (c.29T>C, c.152G>T, c.203G>A, c.228C>T, c.458C>T and c.600A>G) and seven new NAT2 alleles were identified with different frequencies in Rio de Janeiro and Goiás. All new SNPs were found as singletons (observed only once in 808 genes) and were confirmed by three independent technical replicates. Molecular modeling and structural analysis suggested that p.Gly51Val variant may have an important effect on substrate recognition by NAT2. We also observed that amino acid change p.Cys68Tyr would affect acetylating activity due to the resulting geometric restrictions and incompatibility of the functional group in the Tyr side chain with the admitted chemical mechanism for catalysis by NATs. Moreover, other variants, such like p.Thr153Ile, p.Thr193Met, p.Pro228Leu and p.Val280Met, may lead to the presence of hydrophobic residues on NAT2 surface involved in protein aggregation and/or targeted degradation. Finally, the new alleles NAT2*6H and NAT2*5N, which showed the highest frequency in the Brazilian populations considered in this study, may code for a slow activity. Functional studies are needed to clarify the mechanisms by which new SNPs interfere with acetylation.

AnEnPi: identification and annotation of analogous enzymes.

BACKGROUND: Enzymes are responsible for the catalysis of the biochemical reactions in metabolic pathways. Analogous enzymes are able to catalyze the same reactions, but they present no significant sequence similarity at the primary level, and possibly different tertiary structures as well. They are thought to have arisen as the result of independent evolutionary events. A detailed study of analogous enzymes may reveal new catalytic mechanisms, add information about the origin and evolution of biochemical pathways and disclose potential targets for drug development. RESULTS: In this work, we have constructed and implemented a new approach, AnEnPi (the Analogous Enzyme Pipeline), using a combination of bioinformatics tools like BLAST, HMMer, and in-house scripts, to assist in the identification, annotation, comparison and study of analogous and homologous enzymes. The algorithm for the detection of analogy is based i) on the construction of groups of homologous enzymes and ii) on the identification of cases where a given enzymatic activity is performed by two or more proteins without significant similarity between their primary structures. We applied this approach to a dataset obtained from KEGG Comprising all annotated enzymes, which resulted in the identification of 986 EC classes where putative analogy was detected (40.5% of all EC classes). AnEnPi is of considerable value in the construction of initial datasets that can be further curated, particularly in gene and genome annotation, in studies involving molecular evolution and metabolism and in the identification of new potential drug targets. CONCLUSION: AnEnPi is an efficient tool for detection and annotation of analogous enzymes and other enzymes in whole genomes. It is available for academic use at: http://bioinfo.pdtis.fiocruz.br/AnEnPi/

ReRep: computational detection of repetitive sequences in genome survey sequences (GSS).

BACKGROUND: Genome survey sequences (GSS) offer a preliminary global view of a genome since, unlike ESTs, they cover coding as well as non-coding DNA and include repetitive regions of the genome. A more precise estimation of the nature, quantity and variability of repetitive sequences very early in a genome sequencing project is of considerable importance, as such data strongly influence the estimation of genome coverage, library quality and progress in scaffold construction. Also, the elimination of repetitive sequences from the initial assembly process is important to avoid errors and unnecessary complexity. Repetitive sequences are also of interest in a variety of other studies, for instance as molecular markers. RESULTS: We designed and implemented a straightforward pipeline called ReRep, which combines bioinformatics tools for identifying repetitive structures in a GSS dataset. In a case study, we first applied the pipeline to a set of 970 GSSs, sequenced in our laboratory from the human pathogen Leishmania braziliensis, the causative agent of leishmaniosis, an important public health problem in Brazil. We also verified the applicability of ReRep to new sequencing technologies using a set of 454-reads of an Escheria coli. The behaviour of several parameters in the algorithm is evaluated and suggestions are made for tuning of the analysis. CONCLUSION: The ReRep approach for identification of repetitive elements in GSS datasets proved to be straightforward and efficient. Several potential repetitive sequences were found in a L. braziliensis GSS dataset generated in our laboratory, and further validated by the analysis of a more complete genomic dataset from the EMBL and Sanger Centre databases. ReRep also identified most of the E. coli K12 repeats prior to assembly in an example dataset obtained by automated sequencing using 454 technology. The parameters controlling the algorithm behaved consistently and may be tuned to the properties of the dataset, in particular to the length of sequencing reads and the genome coverage. ReRep is freely available for academic use at http://bioinfo.pdtis.fiocruz.br/ReRep/.

Squid - a simple bioinformatics grid.

BACKGROUND: BLAST is a widely used genetic research tool for analysis of similarity between nucleotide and protein sequences. This paper presents a software application entitled "Squid" that makes use of grid technology. The current version, as an example, is configured for BLAST applications, but adaptation for other computing intensive repetitive tasks can be easily accomplished in the open source version. This enables the allocation of remote resources to perform distributed computing, making large BLAST queries viable without the need of high-end computers. RESULTS: Most distributed computing / grid solutions have complex installation procedures requiring a computer specialist, or have limitations regarding operating systems. Squid is a multi-platform, open-source program designed to "keep things simple" while offering high-end computing power for large scale applications. Squid also has an efficient fault tolerance and crash recovery system against data loss, being able to re-route jobs upon node failure and recover even if the master machine fails. Our results show that a Squid application, working with N nodes and proper network resources, can process BLAST queries almost N times faster than if working with only one computer. CONCLUSION: Squid offers high-end computing, even for the non-specialist, and is freely available at the project web site. Its open-source and binary Windows distributions contain detailed instructions and a "plug-n-play" instalation containing a pre-configured example.

Genotypic variation and stability of four variable-number tandem repeats and their suitability for discriminating strains of Mycobacterium leprae.

It has not been possible to distinguish different strains of Mycobacterium leprae according to their genetic sequence. However, the genome contains several variable-number tandem repeats (VNTR), which have been used effectively in strain typing of other bacteria. To determine their suitability for differentiating M. leprae, we developed PCR systems to amplify 5 different VNTR loci and examined a battery of 12 M. leprae strains derived from patients in different regions of the United States, Brazil, Mexico, and the Philippines, as well as from wild armadillos and a sooty mangabey monkey. We found diversity at four VNTR (D = 0.74), but one system (C(16)G(8)) failed to yield reproducible results. Alleles for the GAA VNTR varied in length from 10 to 16 copies, those for AT(17) varied in length from 10 to 15 copies, those for GTA varied in length from 9 to 12 copies, and those for TA(18) varied in length from 13 to 20 copies. Relatively little variation was seen with interspecies transfer of bacilli or during short-term passage of strains in nude mice or armadillos. The TA(18) locus was more polymorphic than other VNTR, and genotypic variation was more common after long-term expansion in armadillos. Most strain genotypes remained fairly stable in passage, but strain Thai-53 showed remarkable variability. Statistical cluster analysis segregated strains and passage samples appropriately but did not reveal any particular genotype associable with different regions or hosts of origin. VNTR polymorphisms can be used effectively to discriminate M. leprae strains. Inclusion of additional loci and other elements will likely lead to a robust typing system that can be used in community-based epidemiological studies and select clinical applications.