Apuã C.M. Paquola

Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni

Schistosoma mansoni is the primary causative agent of schistosomiasis, which affects 200 million individuals in 74 countries. We generated 163,000 expressed-sequence tags (ESTs) from normalized cDNA libraries from six selected developmental stages of the parasite, resulting in 31,000 assembled sequences and 92% sampling of an estimated 14,000 gene complement. By analyzing automated Gene Ontology assignments, we provide a detailed view of important S. mansoni biological systems, including characterization of metazoa-specific and eukarya-conserved genes. Phylogenetic analysis suggests an early divergence from other metazoa. The data set provides insights into the molecular mechanisms of tissue organization, development, signaling, sexual dimorphism, host interactions and immune evasion and identifies novel proteins to be investigated as vaccine candidates and potential drug targets.

Antisense intronic non-coding RNA levels correlate to the degree of tumor differentiation in prostate cancer

A large fraction of transcripts are expressed antisense to introns of known genes in the human genome. Here we show the construction and use of a cDNA microarray platform enriched in intronic transcripts to assess their biological relevance in pathological conditions. To validate the approach, prostate cancer was used as a model, and 27 patient tumor samples with Gleason scores ranging from 5 to 10 were analyzed. We find that a considerably higher fraction (6.6%, [23/346]) of intronic transcripts are significantly correlated (P⩽0.001) to the degree of prostate tumor differentiation (Gleason score) when compared to transcripts from unannotated genomic regions (1%, [6/539]) or from exons of known genes (2%, [27/1369]). Among the top twelve transcripts most correlated to tumor differentiation, six are antisense intronic messages as shown by orientation-specific RT-PCR or Northern blot analysis with strand-specific riboprobe. Orientation-specific real-time RT–PCR with six tumor samples, confirmed the correlation (P=0.024) between the low/high degrees of tumor differentiation and antisense intronic RASSF1 transcript levels. The need to use intron arrays to reveal the transcriptome profile of antisense intronic RNA in cancer has clearly emerged.

Large-scale Transcriptome Analyses Reveal New Genetic Marker Candidates of Head, Neck, and Thyroid Cancer

A detailed genome mapping analysis of 213,636 expressed sequence tags (EST) derived from nontumor and tumor tissues of the oral cavity, larynx, pharynx, and thyroid was done. Transcripts matching known human genes were identified; potential new splice variants were flagged and subjected to manual curation, pointing to 788 putatively new alternative splicing isoforms, the majority (75%) being insertion events. A subset of 34 new splicing isoforms (5% of 788 events) was selected and 23 (68%) were confirmed by reverse transcription-PCR and DNA sequencing. Putative new genes were revealed, including six transcripts mapped to well-studied chromosomes such as 22, as well as transcripts that mapped to 253 intergenic regions. In addition, 2,251 noncoding intronic RNAs, eventually involved in transcriptional regulation, were found. A set of 250 candidate markers for loss of heterozygosis or gene amplification was selected by identifying transcripts that mapped to genomic regions previously known to be frequently amplified or deleted in head, neck, and thyroid tumors. Three of these markers were evaluated by quantitative reverse transcription-PCR in an independent set of individual samples. Along with detailed clinical data about tumor origin, the information reported here is now publicly available on a dedicated Web site as a resource for further biological investigation. This first in silico reconstruction of the head, neck, and thyroid transcriptomes points to a wealth of new candidate markers that can be used for future studies on the molecular basis of these tumors. Similar analysis is warranted for a number of other tumors for which large EST data sets are available.

Characterization of the SOS Regulon of Caulobacter crescentus

The SOS regulon is a paradigm of bacterial responses to DNA damage. A wide variety of bacterial species possess homologs of lexA and recA, the central players in the regulation of the SOS circuit. Nevertheless, the genes actually regulated by the SOS have been determined only experimentally in a few bacterial species. In this work, we describe 37 genes regulated in a LexA-dependent manner in the alphaproteobacterium Caulobacter crescentus. In agreement with previous results, we have found that the direct repeat GTTCN7GTTC is the SOS operator of C. crescentus, which was confirmed by site-directed mutagenesis studies of the imuA promoter. Several potential promoter regions containing the SOS operator were identified in the genome, and the expression of the corresponding genes was analyzed for both the wild type and the lexA strain, demonstrating that the vast majority of these genes are indeed SOS regulated. Interestingly, many of these genes encode proteins with unknown functions, revealing the potential of this approach for the discovery of novel genes involved in cellular responses to DNA damage in prokaryotes, and illustrating the diversity of SOS-regulated genes among different bacterial species.

Laterally transferred genomic islands in Xanthomonadales related to pathogenicity and primary metabolism

Lateral gene transfer (LGT) is considered as one of the drivers in bacterial genome evolution, usually associated with increased fitness and/or changes in behavior, especially if one considers pathogenic vs. non-pathogenic bacterial groups. The genomes of two phytopathogens, Xanthomonas campestris pv. campestris and Xanthomonas axonopodis pv. citri, were previously inspected for genome islands originating from LGT events, and, in this work, potentially early and late LGT events were identified according to their altered nucleotide composition. The biological role of the islands was also assessed, and pathogenicity, virulence and secondary metabolism pathways were functions highly represented, especially in islands that were found to be recently transferred. However, old islands are composed of a high proportion of genes related to cell primary metabolic functions. These old islands, normally undetected by traditional atypical composition analysis, but confirmed as product of LGT by atypical phylogenetic reconstruction, reveal the role of LGT events by replacing core metabolic genes normally inherited by vertical processes.

ESTWeb: bioinformatics services for EST sequencing projects.

ESTWeb is an internet based software package designed for uniform data processing and storage for large-scale EST sequencing projects. The package provides for: (a) reception of sequencing chromatograms; (b) sequence processing such as base-calling, vector screening, comparison with public databases; (c) storage of data and analysis in a relational database, (d) generation of a graphical report of individual sequence quality; and (e) issuing of reports with statistics of productivity and redundancy. The software facilitates real-time monitoring and evaluation of EST sequence acquisition progress along an EST sequencing project.

A quantitative view of the transcriptome of Schistosoma mansoni adult-worms using SAGE

BackgroundFive species of the genus Schistosoma, a parasitic trematode flatworm, are causative agents of Schistosomiasis, a disease that is endemic in a large number of developing countries, affecting millions of patients around the world. By using SAGE (Serial Analysis of Gene Expression) we describe here the first large-scale quantitative analysis of the Schistosoma mansoni transcriptome, one of the most epidemiologically relevant species of this genus.ResultsAfter extracting mRNA from pooled male and female adult-worms, a SAGE library was constructed and sequenced, generating 68,238 tags that covered more than 6,000 genes expressed in this developmental stage. An analysis of the ordered tag-list shows the genes of F10 eggshell protein, pol-polyprotein, HSP86, 14-3-3 and a transcript yet to be identified to be the five top most abundant genes in pooled adult worms. Whereas only 8% of the 100 most abundant tags found in adult worms of S. mansoni could not be assigned to transcripts of this parasite, 46.9% of the total ditags could not be mapped, demonstrating that the 3 sequence of most of the rarest transcripts are still to be identified. Mapping of our SAGE tags to S. mansoni genes suggested the occurrence of alternative-polyadenylation in at least 13 gene transcripts. Most of these events seem to shorten the 3 UTR of the mRNAs, which may have consequences over their stability and regulation.ConclusionSAGE revealed the frequency of expression of the majority of the S. mansoni genes. Transcriptome data suggests that alternative polyadenylation is likely to be used in the control of mRNA stability in this organism. When transcriptome was compared with the proteomic data available, we observed a correlation of about 50%, suggesting that both transcriptional and post-transcriptional regulation are important for determining protein abundance in S. mansoni. The generation of SAGE tags from other life-cycle stages should contribute to reveal the dynamics of gene expression in this important parasite.

Transcriptional profiles of unirradiated or UV-irradiated human cells expressing either the cancer-prone XPB/CS allele or the noncancer-prone XPB/TTD allele

Xeroderma pigmentosum (XP) and trichothiodystrophy (TTD) syndromes are characterized by deficiency in nucleotide excision repair pathway, but with distinguished clinical manifestations. While XP patients exhibit a high frequency of skin cancer, TTD patients are not cancer prone. The relation between lack of DNA repair and their clinical manifestations was investigated through analysis of the transcriptional profile of 12 600 transcripts in two isogenic cell lines with different capabilities of DNA repair. These cell lines result from a stable transfection of the XPB-TTD allele into XP complementation group B fibroblasts, from an XP patient who also have clinical abnormalities corresponding to Cockaynes syndrome (CS). The microarray assays performed under normal growth conditions showed the expression of distinct groups of genes in each cell line. The UVC-transcription modulation of these cells revealed the changes in 869 transcripts. Some of these transcripts had similar modulation pattern in both cells, although with eventually different time patterns for induction or repression. However, some different ‘UVC signature’ for each cell line was also found, that is, transcripts that were specifically UV regulated depending on the DNA repair status of the cell. These results provide a detailed portrait of expression profiles that may potentially unravel the causes of the different phenotypes of XP/CS and TTD patients.

Zerg: a very fast BLAST parser library

SUMMARY Zerg is a library of sub-routines that parses the output from all NCBI BLAST programs (Blastn, Blastp, Blastx, Tblastn and Tblastx) and returns the attributes of a BLAST report to the user. It is optimized for speed, being especially useful for large-scale genomic analysis. Benchmark tests show that Zerg is over two orders of magnitude faster than some widely used BLAST parsers. AVAILABILITY http://bioinfo.iq.usp.br/zerg

Horizontal Gene Transfer Building Prokaryote Genomes: Genes Related to Exchange Between Cell and Environment are Frequently Transferred

Horizontal gene transfer (HGT) has a major impact on the evolution of prokaryotic genomes, as it allows genes evolved in different contexts to be combined in a single genome, greatly enhancing the ways evolving organisms can explore the gene content space and adapt to the environment. A systematic analysis of HGT in a large number of genomes is of key importance in understanding the impact of HGT in the evolution of prokaryotes. We developed a method for the detection of genes that potentially originated by HGT based on the comparison of BLAST scores between homologous genes to 16S rRNA-based phylogenetic distances between the involved organisms. The approach was applied to 697 prokaryote genomes and estimated that in average approximately 15% of the genes in prokaryote genomes originated by HGT, with a clear correlation between the proportion of predicted HGT genes and the size of the genome. The methodology was strongly supported by evolutionary relationships, as tested by the direct phylogenetic reconstruction of many of the HGT candidates. Studies performed with Escherichia coli W3110 genome clearly show that HGT proteins have fewer interactions when compared to those predicted as vertical inherited, an indication that the number of protein partners imposes limitations to horizontal transfer. A detailed functional classification confirms that genes related to protein translation are vertically inherited, whereas interestingly, transport and binding proteins are strongly enriched among HGT genes. Because these genes are related to the cell exchange with their environment, their transfer most likely contributed to successful adaptation throughout evolution.