GENCODE: producing a reference annotation for ENCODE.

BACKGROUND: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manual annotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results. RESULTS: The GENCODE gene features are divided into eight different categories of which only the first two (known and novel coding sequence) are confidently predicted to be protein-coding genes. 5' rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentally verify the initial annotation. Of the 420 coding loci tested, 229 RACE products have been sequenced. They supported 5' extensions of 30 loci and new splice variants in 50 loci. In addition, 46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15 putative transcripts. We assessed the comprehensiveness of the GENCODE annotation by attempting to validate all the predicted exon boundaries outside the GENCODE annotation. Out of 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only two of them in intergenic regions. CONCLUSION: In total, 487 loci, of which 434 are coding, have been annotated as part of the GENCODE reference set available from the UCSC browser. Comparison of GENCODE annotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained within the two sets, which is a reflection of the high number of alternative splice forms with unique exons annotated. Over 50% of coding loci have been experimentally verified by 5' RACE for EGASP and the GENCODE collaboration is continuing to refine its annotation of 1% human genome with the aid of experimental validation.

Análise do perfil de expressão dos genes da cana-de-açúcar envolvidos na interação com Leifsonia xyli subsp: xyli

Utilizando a técnica de macroarranjos de cDNA em membranas de náilon, analisou-se o perfil de expressão de 3.575 ESTs ("Expressed Sequence Tags") de cana-de-açúcar, oriundas do projeto SUCEST, em duas variedades, uma tolerante (SP80-0185) e outra suscetível (SP70-3370) ao Raquitismo da Soqueira. Foram analisadas amostras foliares de plantas inoculadas com Leifsonia xyli subsp. xyli., agente etiológico do Raquitismo, contrastadas com plantas não inoculadas (controle), para cada variedade, marcadas com sondas de cDNA e hibridizadas contra os macroarranjos. Após as hibridizações e análises estatísticas dos dados foi possível identificar 49 ESTs com expressão alterada, sendo 44 na variedade tolerante (41 ESTs induzidos e 3 reprimidos) e 5 na variedade suscetível (2 ESTs induzidos e 3 reprimidos). Os resultados obtidos sugerem que a tolerância da variedade SP80-0185 de cana-de-açúcar à bactéria fitopatogênica pode estar relacionada com a percepção de sinais extracelulares, visto que ESTs relacionados a vias de transdução de sinais apresentaram expressão gênica induzida na variedade tolerante, os quais codificam para uma EST com similaridade à H+-ATPase da membrana plasmática, fatores de transcrição G-box, OsNAC6, "DNA binding", família MYB e "Zinc Finger" e ainda uma EST com similaridade ao fator de ligação ao G-Box, o qual corresponde a uma seqüência de DNA cis presente em vários promotores de plantas e requerido para o reconhecimento de muitos estímulos ambientais. Na variedade suscetível foi reprimido uma EST com similaridade à lipase. Esta enzima, também de membrana, faz parte da síntese do jasmonato, o qual ativa as defesas vegetais contra patógenos de plantas. Possíveis funções para os genes induzidos ou reprimidos nas cultivares de cana tolerante ou resistente ao Raquitismo são discutidas neste trabalho.

Eucalyptus ESTs involved in mechanisms against plant pathogens and environmental stresses

Eucalypt plantation has high economical importance in Brazil; however, it has been attacked by various pathogens under different environmental stress conditions. Disease resistance and survival under unfavorable environmental conditions have revealed that the eucalypt has developed highly efficient defense systems. Here we show the results of the Eucalyptus ESTs Genome Project (FORESTs). Using the expressed sequence tags (ESTs) obtained by the Project, contigs of similar sequences from each cDNA library induced and not induced by stress agents were formed, and cDNA sequences similar to other already known molecules, such as plant-signaling molecules, phytoalexins, lignin biosynthesis pathways, PR-proteins and putative genes corresponding to enzymes involved in the detoxification of reactive oxygen species, were identified. We also present general considerations about the mechanisms of Eucalyptus defense against biotic and abiotic stresses. These data are of extreme importance for future eucalypt breeding programs aimed at developing plants with enhanced resistance against pathogens and environmental stresses.

A tobacco cDNA reveals two different transcription patterns in vegetative and reproductive organs

In order to identify genes expressed in the pistil that may have a role in the reproduction process, we have established an expressed sequence tags project to randomly sequence clones from a Nicotiana tabacum stigma/style cDNA library. A cDNA clone (MTL-8) showing high sequence similarity to genes encoding glycine-rich RNA-binding proteins was chosen for further characterization. Based on the extensive identity of MTL-8 to the RGP-1a sequence of N. sylvestris, a primer was defined to extend the 5' sequence of MTL-8 by RT-PCR from stigma/style RNAs. The amplification product was sequenced and it was confirmed that MTL-8 corresponds to an mRNA encoding a glycine-rich RNA-binding protein. Two transcripts of different sizes and expression patterns were identified when the MTL-8 cDNA insert was used as a probe in RNA blots. The largest is 1,100 nucleotides (nt) long and markedly predominant in ovaries. The smaller transcript, with 600 nt, is ubiquitous to the vegetative and reproductive organs analyzed (roots, stems, leaves, sepals, petals, stamens, stigmas/styles and ovaries). Plants submitted to stress (wounding, virus infection and ethylene treatment) presented an increased level of the 600-nt transcript in leaves, especially after tobacco necrosis virus infection. In contrast, the level of the 1,100-nt transcript seems to be unaffected by the stress conditions tested. Results of Southern blot experiments have suggested that MTL-8 is present in one or two copies in the tobacco genome. Our results suggest that the shorter transcript is related to stress while the larger one is a flower predominant and nonstress-inducible messenger.

Pipeline for macro- and microarray analyses

The pipeline for macro- and microarray analyses (PMmA) is a set of scripts with a web interface developed to analyze DNA array data generated by array image quantification software. PMmA is designed for use with single- or double-color array data and to work as a pipeline in five classes (data format, normalization, data analysis, clustering, and array maps). It can also be used as a plugin in the BioArray Software Environment, an open-source database for array analysis, or used in a local version of the web service. All scripts in PMmA were developed in the PERL programming language and statistical analysis functions were implemented in the R statistical language. Consequently, our package is a platform-independent software. Our algorithms can correctly select almost 90% of the differentially expressed genes, showing a superior performance compared to other methods of analysis. The pipeline software has been applied to 1536 expressed sequence tags macroarray public data of sugarcane exposed to cold for 3 to 48 h. PMmA identified thirty cold-responsive genes previously unidentified in this public dataset. Fourteen genes were up-regulated, two had a variable expression and the other fourteen were down-regulated in the treatments. These new findings certainly were a consequence of using a superior statistical analysis approach, since the original study did not take into account the dependence of data variability on the average signal intensity of each gene. The web interface, supplementary information, and the package source code are available, free, to non-commercial users at http://ipe.cbmeg.unicamp.br/pub/PMmA.

In silico analysis of mitochondrial proteins

Le rôle important joué par la mitochondrie dans la cellule eucaryote est admis depuis longtemps. Cependant, la composition exacte des mitochondries, ainsi que les processus biologiques qui sy déroulent restent encore largement inconnus. Deux facteurs principaux permettent dexpliquer pourquoi létude des mitochondries progresse si lentement : le manque defficacité des méthodes didentification des protéines mitochondriales et le manque de précision dans lannotation de ces protéines. En conséquence, nous avons développé un nouvel outil informatique, YimLoc, qui permet de prédire avec succès les protéines mitochondriales à partir des séquences génomiques. Cet outil intègre plusieurs indicateurs existants, et sa performance est supérieure à celle des indicateurs considérés individuellement. Nous avons analysé environ 60 génomes fongiques avec YimLoc afin de lever la controverse concernant la localisation de la bêta-oxydation dans ces organismes. Contrairement à ce qui était généralement admis, nos résultats montrent que la plupart des groupes de Fungi possèdent une bêta-oxydation mitochondriale. Ce travail met également en évidence la diversité des processus de bêta-oxydation chez les champignons, en corrélation avec leur utilisation des acides gras comme source dénergie et de carbone. De plus, nous avons étudié le composant clef de la voie de bêta-oxydation mitochondriale, lacyl-CoA déshydrogénase (ACAD), dans 250 espèces, couvrant les 3 domaines de la vie, en combinant la prédiction de la localisation subcellulaire avec la classification en sous-familles et linférence phylogénétique. Notre étude suggère que les gènes ACAD font partie dune ancienne famille qui a adopté des stratégies évolutionnaires innovatrices afin de générer un large ensemble denzymes susceptibles dutiliser la plupart des acides gras et des acides aminés. Finalement, afin de permettre la prédiction de protéines mitochondriales à partir de données autres que les séquences génomiques, nous avons développé le logiciel TESTLoc qui utilise comme données des Expressed Sequence Tags (ESTs). La performance de TESTLoc est significativement supérieure à celle de tout autre outil de prédiction connu. En plus de fournir deux nouveaux outils de prédiction de la localisation subcellulaire utilisant différents types de données, nos travaux démontrent comment lassociation de la prédiction de la localisation subcellulaire à dautres méthodes danalyse in silico permet daméliorer la connaissance des protéines mitochondriales. De plus, ces travaux proposent des hypothèses claires et faciles à vérifier par des expériences, ce qui présente un grand potentiel pour faire progresser nos connaissances des métabolismes mitochondriaux.

Implication des peptides RALFs dans les communications cellulaires lors du développement du gamétophyte femelle chez Solanum chacoense et Arabidopsis thaliana.

Chez les angiospermes, la reproduction passe par la double fécondation. Le tube pollinique délivre deux cellules spermatiques au sein du gamétophyte femelle. Une cellule féconde la cellule œuf pour produire un zygote; l’autre féconde la cellule centrale pour produire l’endosperme. Pour assurer un succès reproductif, le développement du gamétophyte femelle au sein de l’ovule doit établir un patron cellulaire qui favorise les interactions avec le tube pollinique et les cellules spermatiques. Pour ce faire, un dialogue doit s’établir entre les différentes cellules de l’ovule lors de son développement, de même que lors de la fécondation. D’ailleurs, plusieurs types de communications intercellulaires sont supposées suite à la caractérisation de plusieurs mutants développementaux. De même, ces communications semblent persister au sein du zygote et de l’endosperme pour permettre la formation d’un embryon viable au sein de la graine. Malgré les développements récents qui ont permis de trouver des molécules de signalisation supportant les modèles d’interactions cellulaires avancés par la communauté scientifique, les voies de signalisation sont de loin très incomplètes. Dans le but de caractériser des gènes encodant des protéines de signalisation potentiellement impliqués dans la reproduction chez Solanum chacoense, l’analyse d’expression des gènes de type RALF présents dans une banque d’ESTs (Expressed Sequence Tags) spécifiques à l’ovule après fécondation a été entreprise. RALF, Rapid Alcalinization Factor, est un peptide de 5 kDa qui fait partie de la superfamille des «protéines riches en cystéines (CRPs)», dont les rôles physiologiques au sein de la plante sont multiples. Cette analyse d’expression a conduit à une analyse approfondie de ScRALF3, dont l’expression au sein de la plante se limite essentiellement à l’ovule. L’analyse de plantes transgéniques d’interférence pour le gène ScRALF3 a révélé un rôle particulier lors de la mégagamétogénèse. Les plantes transgéniques présentent des divisions mitotiques anormales qui empêchent le développement complet du sac embryonnaire. Le positionnement des noyaux, de même que la synchronisation des divisions au sein du syncytium, semblent responsables de cette perte de progression lors de la mégagamétogénèse. L’isolement du promoteur de même que l’analyse plus précise d’expression au sein de l’ovule révèle une localisation sporophytique du transcrit. La voie de signalisation de l’auxine régule également la transcription de ScRALF3. De surcroît, ScRALF3 est un peptide empruntant la voie de sécrétion médiée par le réticulum endoplasmique et l’appareil de Golgi. En somme, ScRALF3 est un important facteur facilitant la communication entre le sporophyte et le gamétophyte pour amener à maturité le sac embryonnaire. L’identification d’un orthologue potentiel chez Arabidopsis thaliana a conduit à la caractérisation de AtRALF34. L’absence de phénotype lors du développement du sac embryonnaire suggère, cependant, de la redondance génétique au sein de la grande famille des gènes de type RALF. Néanmoins, les peptides RALFs apparaissent comme d’importants régulateurs lors de la reproduction chez Solanum chacoense et Arabidopsis thaliana.

Transcriptome analysis for discovering candidate genes involve in embryogenesis in coconut (Cocos nucifera L.) through 454 pyrosequencing

Coconut, Cocos nucifera L. is a major plantation crop, which ensures income for millions of people in the tropical region. Detailed molecular studies on zygotic embryo development would provide valuable clues for the identification of molecular markers to improve somatic embryogenesis. Since there is no ongoing genome project for this species, coconut expressed sequence tags (EST) would be an interesting technique to identify important coconut embryo specific genes as well as other functional genes in different biochemical pathways. The goal of this study was to analyse the ESTs by examining the transcriptome data of the different embryo tissue types together with one somatic tissue. Here, four cDNA libraries from immature embryo, mature embryo, microspore derived embryo and mature leaves were constructed. cDNA was sequenced by the Roche-454 GS-FLX system and assembled into 32621 putative unigenes and 155017 singletons. Of these unigenes, 18651 had significant sequence similarities to non-redundant protein database, from which 16153 were assigned to one or more gene ontology categories. Homologue genes, which are responsible for embryo development such as chitinase, beta-1,3-glucanase, ATP synthase CF0 subunit, thaumatin-like protein and metallothionein-like protein were identified among the embryo EST collection. Of the unigenes, 6694 were mapped into 139 KEGG pathways including carbohydrate metabolism, energy metabolism, lipid metabolism, amino acid metabolism and nucleotide metabolism. This collection of 454-derived EST data generated from different tissue types provides a significant resource for genome wide studies and gene discovery of coconut, a non-model species.

Evidence of a Ca(2+)-(center dot)NO-cGMP signaling pathway controlling zoospore biogenesis in the aquatic fungus Blastocladiella emersonii

The sporulation stage of the aquatic fungus Blastocladiella emersonii culminates with the formation and release to the medium of a number of zoospores, which are motile cells responsible for the dispersal of the fungus. The presence in the sporulation solution of 1H-[1,2,4]Oxadiazolo[4,3-a]quinoxalin-1-one (ODQ), a potent and selective inhibitor of nitric oxide-sensitive guanylyl cyclases, completely prevented biogenesis of the zoospores. In addition, this compound was able to significantly reduce cGMP levels, which increase drastically during late sporulation, suggesting the existence of a nitric oxide-dependent mechanism for cGMP synthesis. Furthermore, increased levels of nitric oxide-derived products were detected during sporulation by fluorescence assays using DAF-2 DA, whose signal was drastically reduced in the presence of the nitric oxide synthase inhibitor N omega-Nitro-L-arginine methyl ester (L-NAME). These results were confirmed by quantitative chemiluminescent determination of the intracellular levels of nitric oxide-derived products. A putative nitric oxide synthase (NOS) activity was detected throughout sporulation, and this enzyme activity decreased significantly when L-NAME and 1-[2-(Trifluoromethyl)phenyl]imidazole (TRIM) were added to the assays. NOS assays carried out in the presence of EGTA showed decreased enzyme activity, suggesting the involvement of calcium ions in enzyme activation. Additionally, expressed sequence tags (ESTs) encoding putative guanylyl cyclases and a cGMP-phosphodiesterase were found in B. emersonii EST database (http://blasto.iq.usp.br), and the mRNA levels of the corresponding genes were observed to increase during sporulation. Altogether, data presented here revealed the presence and expression of guanylyl cyclase and cGMP phosphodiesterase genes in B. emersonii and provided evidence of a Ca(2+)-(center dot)NO-cGMP signaling pathway playing a role in zoospore biogenesis. (C) 2009 Elsevier Inc. All rights reserved.

The Aedes aegypti larval transcriptome: a comparative perspective with emphasis on trypsins and the domain structure of peritrophins

The genome sequence of Aedes aegypti was recently reported. A significant amount of Expressed Sequence Tags (ESTs) were sequenced to aid in the gene prediction process. In the present work we describe an integrated analysis of the genomic and EST data, focusing on genes with preferential expression in larvae (LG), adults (AG) and in both stages (SG). A total of 913 genes (5.4% of the transcript complement) are LG, including ion transporters and cuticle proteins that are important for ion homeostasis and defense. From a starting set of 245 genes encoding the trypsin domain, we identified 66 putative LG, AG, and SG trypsins by manual curation. Phylogenetic analyses showed that AG trypsins are divergent from their larval counterparts (LG), grouping with blood-induced trypsins from Anopheles gambiae and Simulium vittatum. These results support the hypothesis that blood-feeding arose only once, in the ancestral Culicomorpha. Peritrophins are proteins that interlock chitin fibrils to form the peritrophic membrane (PM) that compartmentalizes the food in the midgut. These proteins are recognized by having chitin-binding domains with 6 conserved Cys and may also present mucin-like domains (regions expected to be highly O-glycosylated). PM may be formed by a ring of cells (type 2, seen in Ae. aegypti larvae and Drosophila melanogaster) or by most midgut cells (type 1, found in Ae. aegypti adult and Tribolium castaneum). LG and D. melanogaster peritrophins have more complex domain structures than AG and T. castaneum peritrophins. Furthermore, mucin-like domains of peritrophins from T. castaneum (feeding on rough food) are lengthier than those of adult Ae. aegypti (blood-feeding). This suggests, for the first time, that type 1 and type 2 PM may have variable molecular architectures determined by different peritrophins and/or ancillary proteins, which may be partly modulated by diet.

Transcriptome analysis of Taenia solium cysticerci using Open Reading Frame ESTs (ORESTES)

Background: Human infection by the pork tapeworm Taenia solium affects more than 50 million people worldwide, particularly in underdeveloped and developing countries. Cysticercosis which arises from larval encystation can be life threatening and difficult to treat. Here, we investigate for the first time the transcriptome of the clinically relevant cysticerci larval form. Results: Using Expressed Sequence Tags (ESTs) produced by the ORESTES method, a total of 1,520 high quality ESTs were generated from 20 ORESTES cDNA mini-libraries and its analysis revealed fragments of genes with promising applications including 51 ESTs matching antigens previously described in other species, as well as 113 sequences representing proteins with potential extracellular localization, with obvious applications for immune-diagnosis or vaccine development. Conclusion: The set of sequences described here will contribute to deciphering the expression profile of this important parasite and will be informative for the genome assembly and annotation, as well as for studies of intra- and inter-specific sequence variability. Genes of interest for developing new diagnostic and therapeutic tools are described and discussed.

Análise do perfil de expressão dos genes da cana-de-açúcar envolvidos na interação com Leifsonia xyli subsp: xyli

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genes diferencialmente expressos em cana-de-açúcar inoculada com Xanthomonas albilineans, o agente causal da escaldadura da folha

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Eucalyptus ESTs involved in mechanisms against plant pathogens and environmental stresses

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Comparative transcriptome analysis of Paracoccidioides brasiliensis during in vitro adhesion to type I collagen and fibronectin: identification of potential adhesins

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)