ParAlign: a parallel sequence alignment algorithm for rapid and sensitive database searches

There is a need for faster and more sensitive algorithms for sequence similarity searching in view of the rapidly increasing amounts of genomic sequence data available. Parallel processing capabilities in the form of the single instruction, multiple data (SIMD) technology are now available in common microprocessors and enable a single microprocessor to perform many operations in parallel. The ParAlign algorithm has been specifically designed to take advantage of this technology. The new algorithm initially exploits parallelism to perform a very rapid computation of the exact optimal ungapped alignment score for all diagonals in the alignment matrix. Then, a novel heuristic is employed to compute an approximate score of a gapped alignment by combining the scores of several diagonals. This approximate score is used to select the most interesting database sequences for a subsequent Smith–Waterman alignment, which is also parallelised. The resulting method represents a substantial improvement compared to existing heuristics. The sensitivity and specificity of ParAlign was found to be as good as Smith–Waterman implementations when the same method for computing the statistical significance of the matches was used. In terms of speed, only the significantly less sensitive NCBI BLAST 2 program was found to outperform the new approach. Online searches are available at http://dna.uio.no/search/

The prion protein gene: Identifying regulatory signals using marsupial sequence

The function of the prion protein gene (PRNP) and its normal product PrPC is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyTI, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrPC in signal transduction and synaptic plasticity. (c) 2004 Elsevier B.V. All rights reserved.

Identification of protein coding regions using the modified Gabor-wavelet transform

An important topic in genomic sequence analysis is the identification of protein coding regions. In this context, several coding DNA model-independent methods based on the occurrence of specific patterns of nucleotides at coding regions have been proposed. Nonetheless, these methods have not been completely suitable due to their dependence on an empirically predefined window length required for a local analysis of a DNA region. We introduce a method based on a modified Gabor-wavelet transform (MGWT) for the identification of protein coding regions. This novel transform is tuned to analyze periodic signal components and presents the advantage of being independent of the window length. We compared the performance of the MGWT with other methods by using eukaryote data sets. The results show that MGWT outperforms all assessed model-independent methods with respect to identification accuracy. These results indicate that the source of at least part of the identification errors produced by the previous methods is the fixed working scale. The new method not only avoids this source of errors but also makes a tool available for detailed exploration of the nucleotide occurrence.

Molecular cloning and characterization of a ripening-induced polygalacturonase related to papaya fruit softening

Pulp softening is one of the most remarkable changes during ripening of papaya (Carica papaya) fruit and it is a major cause for post-harvest losses. Although cell wall catabolism has a major influence on papaya fruit, quality information on the gene products involved in this process is limited. A full-length polygalacturonase cDNA (cpPG) was isolated from papaya pulp and used to study gene expression and enzyme activity during normal and ethylene-induced ripening and after exposure of the fruit to 1-MCP. Northern-blot analysis demonstrated that cpPG transcription was strongly induced during ripening and was highly ethylene-dependent. The accumulation of cpPG transcript was paralleled by enzyme activity, and inversely correlated to the pulp firmness. Preliminary in silica analysis of the cpPG genomic sequence revealed the occurrence of putative regulatory motifs in the promoter region that may help to explain the effects of plant hormones and non-abiotic stresses on papaya fruit firmness. This newly isolated cpPG is an important candidate for functional characterization and manipulation to control the process of pulp softening during papaya ripening. (C) 2009 Elsevier Masson SAS. All rights reserved.

Rapid evolution of cancer/testis genes on the X chromosome.

BACKGROUND: Cancer/testis (CT) genes are normally expressed only in germ cells, but can be activated in the cancer state. This unusual property, together with the finding that many CT proteins elicit an antigenic response in cancer patients, has established a role for this class of genes as targets in immunotherapy regimes. Many families of CT genes have been identified in the human genome, but their biological function for the most part remains unclear. While it has been shown that some CT genes are under diversifying selection, this question has not been addressed before for the class as a whole. RESULTS: To shed more light on this interesting group of genes, we exploited the generation of a draft chimpanzee (Pan troglodytes) genomic sequence to examine CT genes in an organism that is closely related to human, and generated a high-quality, manually curated set of human:chimpanzee CT gene alignments. We find that the chimpanzee genome contains homologues to most of the human CT families, and that the genes are located on the same chromosome and at a similar copy number to those in human. Comparison of putative human:chimpanzee orthologues indicates that CT genes located on chromosome X are diverging faster and are undergoing stronger diversifying selection than those on the autosomes or than a set of control genes on either chromosome X or autosomes. CONCLUSION: Given their high level of diversifying selection, we suggest that CT genes are primarily responsible for the observed rapid evolution of protein-coding genes on the X chromosome.

Current millennium biotechniques for biomedical research on parasites and host-parasite interactions

The development of biotechnology in the last three decades has generated the feeling that the newest scientific achievements will deliver high standard quality of life through abundance of food and means for successfully combating diseases. Where the new biotechnologies give access to genetic information, there is a common belief that physiological and pathological processes result from subtle modifications of gene expression. Trustfully, modern genetics has produced genetic maps, physical maps and complete nucleotide sequences from 141 viruses, 51 organelles, two eubacteria, one archeon and one eukaryote (Saccharomices cerevisiae). In addition, during the Centennial Commemoration of the Oswaldo Cruz Institute the nearly complete human genome map was proudly announced, whereas the latest Brazilian key stone contribution to science was the publication of the Shillela fastidiosa genomic sequence highlythed on a Nature cover issue. There exists a belief among the populace that further scientific accomplishments will rapidly lead to new drugs and methodological approaches to cure genetic diseases and other incurable ailments. Yet, much evidence has been accumulated, showing that a large information gap exists between the knowledge of genome sequence and our knowledge of genome function. Now that many genome maps are available, people wish to know what are we going to do with them. Certainly, all these scientific accomplishments will shed light on many more secrets of life. Nevertheless, parsimony in the weekly announcements of promising scientific achievements is necessary. We also need many more creative experimental biologists to discover new, as yet un-envisaged biotechnological approaches, and the basic resource needed for carrying out mile stone research necessary for leading us to that "promised land"often proclaimed by the mass media.

Optimització d'una aplicació bioinformàtica d'aliniament de seqüències executada en processadors many-core (GPUs)

Las herramientas de análisis de secuencias genómicas permiten a los biólogos identificar y entender regiones fundamentales que tienen implicación en enfermedades genéticas. Actualmente existe una necesidad de dotar al ámbito científico de herramientas de análisis eficientes. Este proyecto lleva a cabo una caracterización y análisis del rendimiento de algoritmos utilizados en la comparación de secuencias genómicas completas, y ejecutadas en arquitecturas MultiCore y ManyCore. A partir del análisis se evalúa la idoneidad de este tipo de arquitecturas para resolver el problema de comparar secuencias genómicas. Finalmente se propone una serie de modificaciones en las implementaciones de estos algoritmos con el objetivo de mejorar el rendimiento.

Single-stranded oligonucleotide-mediated in vivo gene repair in the rd1 retina.

PURPOSE: The aim of this study was to test whether oligonucleotide-targeted gene repair can correct the point mutation in genomic DNA of PDE6b(rd1) (rd1) mouse retinas in vivo. METHODS: Oligonucleotides (ODNs) of 25 nucleotide length and complementary to genomic sequence subsuming the rd1 point mutation in the gene encoding the beta-subunit of rod photoreceptor cGMP-phosphodiesterase (beta-PDE), were synthesized with a wild type nucleotide base at the rd1 point mutation position. Control ODNs contained the same nucleotide bases as the wild type ODNs but with varying degrees of sequence mismatch. We previously developed a repeatable and relatively non-invasive technique to enhance ODN delivery to photoreceptor nuclei using transpalpebral iontophoresis prior to intravitreal ODN injection. Three such treatments were performed on C3H/henJ (rd1) mouse pups before postnatal day (PN) 9. Treatment outcomes were evaluated at PN28 or PN33, when retinal degeneration was nearly complete in the untreated rd1 mice. The effect of treatment on photoreceptor survival was evaluated by counting the number of nuclei of photoreceptor cells and by assessing rhodopsin immunohistochemistry on flat-mount retinas and sections. Gene repair in the retina was quantified by allele-specific real time PCR and by detection of beta-PDE-immunoreactive photoreceptors. Confirmatory experiments were conducted using independent rd1 colonies in separate laboratories. These experiments had an additional negative control ODN that contained the rd1 mutant nucleotide base at the rd1 point mutation site such that the sole difference between treatment with wild type and control ODN was the single base at the rd1 point mutation site. RESULTS: Iontophoresis enhanced the penetration of intravitreally injected ODNs in all retinal layers. Using this delivery technique, significant survival of photoreceptors was observed in retinas from eyes treated with wild type ODNs but not control ODNs as demonstrated by cell counting and rhodopsin immunoreactivity at PN28. Beta-PDE immunoreactivity was present in retinas from eyes treated with wild type ODN but not from those treated with control ODNs. Gene correction demonstrated by allele-specific real time PCR and by counts of beta-PDE-immunoreactive cells was estimated at 0.2%. Independent confirmatory experiments showed that retinas from eyes treated with wild type ODN contained many more rhodopsin immunoreactive cells compared to retinas treated with control (rd1 sequence) ODN, even when harvested at PN33. CONCLUSIONS: Short ODNs can be delivered with repeatable efficiency to mouse photoreceptor cells in vivo using a combination of intravitreal injection and iontophoresis. Delivery of therapeutic ODNs to rd1 mouse eyes resulted in genomic DNA conversion from mutant to wild type sequence, low but observable beta-PDE immunoreactivity, and preservation of rhodopsin immunopositive cells in the outer nuclear layer, suggesting that ODN-directed gene repair occurred and preserved rod photoreceptor cells. Effects were not seen in eyes treated with buffer or with ODNs having the rd1 mutant sequence, a definitive control for this therapeutic approach. Importantly, critical experiments were confirmed in two laboratories by several different researchers using independent mouse colonies and ODN preparations from separate sources. These findings suggest that targeted gene repair can be achieved in the retina following enhanced ODN delivery.

A new approach for potential drug target discovery through in silico metabolic pathway analysis using Trypanosoma cruzi genome information

The current drug options for the treatment of chronic Chagas disease have not been sufficient and high hopes have been placed on the use of genomic data from the human parasite Trypanosoma cruzi to identify new drug targets and develop appropriate treatments for both acute and chronic Chagas disease. However, the lack of a complete assembly of the genomic sequence and the presence of many predicted proteins with unknown or unsure functions has hampered our complete view of the parasite's metabolic pathways. Moreover, pinpointing new drug targets has proven to be more complex than anticipated and has revealed large holes in our understanding of metabolic pathways and their integrated regulation, not only for this parasite, but for many other similar pathogens. Using an in silicocomparative study on pathway annotation and searching for analogous and specific enzymes, we have been able to predict a considerable number of additional enzymatic functions in T. cruzi. Here we focus on the energetic pathways, such as glycolysis, the pentose phosphate shunt, the Krebs cycle and lipid metabolism. We point out many enzymes that are analogous to those of the human host, which could be potential new therapeutic targets.

Towards an understanding of the epigenetics of schistosomes: a comparative epigenomic study

As in perhaps all eukaryotes, schistosomes use a supplementary information transmitting system, the epigenetic inheritance system, to shape genetic information and to produce different phenotypes. In contrast to other important parasites, the study of epigenetic phenomena in schistosomes is still in its infancy. Nevertheless, we are beginning to grasp what goes on behind the epigenetic scene in this parasite. We have developed techniques of native chromatin immunoprecipitation (N-ChIP) and associated the necessary bioinformatics tools that allow us to run genome-wide comparative chromatin studies on Schistosoma mansoni at different stages of its life cycle, on different strains and on different sexes. We present here an application of such an approach to study the genetic and epigenetic basis for a phenotypic trait, the compatibility of S. mansoni with its invertebrate host Biomphalaria glabrata. We have applied the ChIP procedure to two strains that are either compatible or incompatible with their intermediate host. The precipitated DNA was sequenced and aligned to a reference genome and this information was used to determine regions in which both strands differ in their genomic sequence and/or chromatin structure. This procedure allowed us to identify candidate genes that display either genetic or epigenetic difference between the two strains.

Genetic diversity of NS3 protease from Brazilian HCV isolates and possible implications for therapy with direct-acting antiviral drugs

The hepatitis C virus (HCV) NS3 protease has been one of the molecular targets of new therapeutic approaches. Its genomic sequence variability in Brazilian HCV isolates is poorly documented. To obtain more information on the magnitude of its genetic diversity, 114 Brazilian HCV samples were sequenced and analysed together with global reference sequences. Genetic distance (d) analyses revealed that subtype 1b had a higher degree of heterogeneity (d = 0.098) than subtypes 1a (d = 0.060) and 3a (d = 0.062). Brazilian isolates of subtype 1b were distributed in the phylogenetic tree among sequences from other countries, whereas most subtype 1a and 3a sequences clustered into a single branch. Additional characterisation of subtype 1a in clades 1 and 2 revealed that all but two Brazilian subtype 1a sequences formed a distinct and strongly supported (approximate likelihood-ratio test = 93) group of sequences inside clade 1. Moreover, this subcluster inside clade 1 presented an unusual phenotypic characteristic in relation to the presence of resistance mutations for macrocyclic inhibitors. In particular, the mutation Q80K was found in the majority of clade 1 sequences, but not in the Brazilian isolates. These data demonstrate that Brazilian HCV subtypes display a distinct pattern of genetic diversity and reinforce the importance of sequence information in future therapeutic approaches.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Genomics of secondary metabolite production by Pseudomonas fluorescens Pf-5

The complete sequence of the 7.07 Mb genome of the biological control agent Pseudomonas fluorescens Pf-5 is now available, providing a new opportunity to advance knowledge of biological control through genomics. P. fluorescens Pf-5 is a rhizosphere bacterium that suppresses seedling emergence diseases and produces a spectrum of antibiotics toxic to plant-pathogenic fungi and oomycetes. In addition to six known secondary metabolites produced by Pf-5, three novel secondary metabolite biosynthesis gene clusters identified in the genome could also contribute to biological control. The genomic sequence provides numerous clues as to mechanisms used by the bacterium to survive in the spermosphere and rhizosphere. These features include broad catabolic and transport capabilities for utilizing seed and root exudates, an expanded collection of efflux systems for defense against environmental stress and microbial competition, and the presence of 45 outer membrane receptors that should allow for the uptake of iron from a wide array of siderophores produced by soil microorganisms. As expected for a bacterium with a large genome that lives in a rapidly changing environment, Pf-5 has an extensive collection of regulatory genes, only some of which have been characterized for their roles in regulation of secondary metabolite production or biological control. Consistent with its commensal lifestyle, Pf-5 appears to lack a number of virulence and pathogenicity factors found in plant pathogen.

Methods for integration site distribution analyses in animal cell genomes.

The question of where retroviral DNA becomes integrated in chromosomes is important for understanding (i) the mechanisms of viral growth, (ii) devising new anti-retroviral therapy, (iii) understanding how genomes evolve, and (iv) developing safer methods for gene therapy. With the completion of genome sequences for many organisms, it has become possible to study integration targeting by cloning and sequencing large numbers of host-virus DNA junctions, then mapping the host DNA segments back onto the genomic sequence. This allows statistical analysis of the distribution of integration sites relative to the myriad types of genomic features that are also being mapped onto the sequence scaffold. Here we present methods for recovering and analyzing integration site sequences.

Method for retinal gene repair in neonatal mouse.

Gene correction at the site of the mutation in the chromosome is the absolute way to really cure a genetic disease. The oligonucleotide (ODN)-mediated gene repair technology uses an ODN perfectly complementary to the genomic sequence except for a mismatch at the base that is mutated. The endogenous repair machinery of the targeted cell then mediates substitution of the desired base in the gene, resulting in a completely normal sequence. Theoretically, it avoids potential gene silencing or random integration associated with common viral gene augmentation approaches and allows an intact regulation of expression of the therapeutic protein. The eye is a particularly attractive target for gene repair because of its unique features (small organ, easily accessible, low diffusion into systemic circulation). Moreover therapeutic effects on visual impairment could be obtained with modest levels of repair. This chapter describes in details the optimized method to target active ODNs to the nuclei of photoreceptors in neonatal mouse using (1) an electric current application at the eye surface (saline transpalpebral iontophoresis), (2) combined with an intravitreous injection of ODNs, as well as the experimental methods for (3) the dissection of adult neural retinas, (4) their immuno-labelling, and (5) flat-mounting for direct observation of photoreceptor survival, a relevant criteria of treatment outcomes for retinal degeneration.