Open-source software accelerates bioinformatics.
A report on the Wellcome Trust/Cold Spring Harbor Genome Informatics meeting, Cold Spring Harbor, USA, 7-11 May 2003.
Computational Biology↗
Biomedical subjects
John Quackenbush
Publications and source records attributed to John Quackenbush.
Regulation of caspase expression and apoptosis by adenomatous polyposis coli.
The adenomatous polyposis coli (APC) gene, a member of the WNT pathway, has been shown to assign intestinal epithelial cells to a program of proliferation or differentiation through regulation of the beta-catenin/TCF-4 complex. Wild-type APC, in certain cellular contexts, appears to induce differentiation and apoptosis, although mutant forms of APC, known to produce polyps and ultimately cancers, may suppress these events. Here, we show that mutant forms of APC can induce repression of select terminal caspases as a potential means of attenuating responses to apoptotic stimuli. Using gene expression profiling to interrogate the intact intestines of Apc(+/min) mice harboring numerous polyps, we identified a reduction in the mRNA expression of both caspases 3 and 7. We additionally identified a reduction in protein levels of caspase-3, caspase-7, and caspase-9 in human colon cancer specimens known to harbor APC mutations. A reduction in caspase protein levels resulted in resistance to apoptotic-inducing agents and restoration of caspase levels reinstated apoptotic capacities. Consistent with Wnt pathway involvement, dominant negative TCF/LEF induced caspase protein expression. These data provide support for the hypothesis that one of the functions of APC is the regulation of caspase activity and other apoptotic proteins by controlling their expression levels in the cell.
Adenomatous Polyposis Coli Protein↗
Shear-induced cyclooxygenase-2 via a JNK2/c-Jun-dependent pathway regulates prostaglandin receptor expression in chondrocytic cells.
Using cDNA microarrays coupled with bioinformatics tools, we elucidated a signaling cascade regulating cyclooxygenase-2 (COX-2), a pivotal pro-inflammatory enzyme expressed in rheumatic and osteoarthritic, but not normal, cartilage. Exposure of T/C-28a2 chondrocytic cells to fluid shear results in co-regulation of c-Jun N-terminal kinase2 (JNK2), c-Jun, and COX-2 as well as concomitant downstream expression of prostaglandin receptors EP2 and EP3a1. JNK2 transcript inhibition abrogated shear-induced COX-2, EP2, and EP3a1 mRNA up-regulation as well as c-Jun phosphorylation. Functional knock-out experiments using an antisense c-Jun oligonucleotide revealed the abolition of shear-induced COX-2, EP2, and EP3a1, but not JNK2, transcripts. Moreover, inhibition of COX-2 activity eliminated mRNA upregulation of EP2 and EP3a1 induced by shear. Hence, a biochemical pathway exists wherein fluid shear activates COX-2, via a JNK2/c-Jun-dependent pathway, which in turn elicits downstream EP2 and EP3a1 mRNA synthesis.
Cell Line↗
TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets.
TGICL is a pipeline for analysis of large Expressed Sequence Tags (EST) and mRNA databases in which the sequences are first clustered based on pairwise sequence similarity, and then assembled by individual clusters (optionally with quality values) to produce longer, more complete consensus sequences. The system can run on multi-CPU architectures including SMP and PVM.
Cluster Analysis↗
Integrating computationally assembled mouse transcript sequences with the Mouse Genome Informatics (MGI) database.
Databases of experimentally generated and computationally derived transcript sequences are valuable resources for genome analysis and annotation. The utility of such databases is enhanced when the sequences they contain are integrated with such biological information as genomic location, gene function, gene expression and phenotypic variation. We present the analysis and results of a semi-automated process of connecting transcript assemblies with highly curated biological information for mouse genes that is available through the Mouse Genome Informatics (MGI) database.
Animals↗
John Quackenbush talks about the clinical promise of genetic microarrays. Interviewed by Brian Vastag.
Explore the source record for details and available documents.
Forecasting↗
Assessing unmodified 70-mer oligonucleotide probe performance on glass-slide microarrays.
BACKGROUND: Long oligonucleotide microarrays are potentially more cost- and management-efficient than cDNA microarrays, but there is little information on the relative performance of these two probe types. The feasibility of using unmodified oligonucleotides to accurately measure changes in gene expression is also unclear. RESULTS: Unmodified sense and antisense 70-mer oligonucleotides representing 75 known rat genes and 10 Arabidopsis control genes were synthesized, printed and UV cross-linked onto glass slides. Printed alongside were PCR-amplified cDNA clones corresponding to the same genes, enabling us to compare the two probe types simultaneously. Our study was designed to evaluate the mRNA profiles of heart and brain, along with Arabidopsis cRNA spiked into the labeling reaction at different relative copy number. Hybridization signal intensity did not correlate with probe type but depended on the extent of UV irradiation. To determine the effect of oligonucleotide concentration on hybridization signal, 70-mers were serially diluted. No significant change in gene-expression ratio or loss in hybridization signal was detected, even at the lowest concentration tested (6.25 microm). In many instances, signal intensity actually increased with decreasing concentration. The correlation coefficient between oligonucleotide and cDNA probes for identifying differentially expressed genes was 0.80, with an average coefficient of variation of 13.4%. Approximately 8% of the genes showed discordant results with the two probe types, and in each case the cDNA results were more accurate, as determined by real-time PCR. CONCLUSIONS: Microarrays of UV cross-linked unmodified oligonucleotides provided sensitive and specific measurements for most of the genes studied.
Arabidopsis↗
The TIGR rice genome annotation resource: annotating the rice genome and creating resources for plant biologists.
Rice is not only a major food staple for the world's population but it also is a model species for a major group of flowering plants, the monocotyledonous plants. Draft genomic sequence of two subspecies of rice, Oryza sativa spp. japonica and indica ssp. are publicly available. To provide the community with a resource to data-mine the rice genome, we have constructed an annotation resource for rice (http://www.tigr.org/tdb/e2k1/osa1/). In this resource, we have annotated the rice genome for gene content, identified motifs/domains within the predicted genes, constructed a rice repeat database, identified related sequences in other plant species, and identified syntenic sequences between rice and maize. All of the data is available through web-based interfaces, FTP downloads, and a Distributed Annotation System.
Chromosomes, Artificial↗
Continued discovery of transcriptional units expressed in cells of the mouse mononuclear phagocyte lineage.
The current RIKEN transcript set represents a significant proportion of the mouse transcriptome but transcripts expressed in the innate and acquired immune systems are poorly represented. In the present study we have assessed the complexity of the transcriptome expressed in mouse macrophages before and after treatment with lipopolysaccharide, a global regulator of macrophage gene expression, using existing RIKEN 19K arrays. By comparison to array profiles of other cells and tissues, we identify a large set of macrophage-enriched genes, many of which have obvious functions in endocytosis and phagocytosis. In addition, a significant number of LPS-inducible genes were identified. The data suggest that macrophages are a complex source of mRNA for transcriptome studies. To assess complexity and identify additional macrophage expressed genes, cDNA libraries were created from purified populations of macrophage and dendritic cells, a functionally related cell type. Sequence analysis revealed a high incidence of novel mRNAs within these cDNA libraries. These studies provide insights into the depths of transcriptional complexity still untapped amongst products of inducible genes, and identify macrophage and dendritic cell populations as a starting point for sampling the inducible mammalian transcriptome.
Animals↗
Gene expression analyses of Arabidopsis chromosome 2 using a genomic DNA amplicon microarray.
The gene predictions and accompanying functional assignments resulting from the sequencing and annotation of a genome represent hypotheses that can be tested and used to develop a more complete understanding of the organism and its biology. In the model plant Arabidopsis thaliana, we developed a novel approach to constructing whole-genome microarrays based on PCR amplification of the 3' ends of each predicted gene from genomic DNA, and constructed an array representing more than 94% of the predicted genes and pseudogenes on chromosome 2. With this array, we examined various tissues and physiological conditions, providing expression-based validation for 84% of the gene predictions and providing clues as to the functions of many predicted genes. Further, by examining the distribution of expression along the physical chromosome, we were able to identify a region of repressed transcription that may represent a previously undescribed heterochromatic region.
Arabidopsis↗
Development and evaluation of an automated annotation pipeline and cDNA annotation system.
Manual curation has long been held to be the "gold standard" for functional annotation of DNA sequence. Our experience with the annotation of more than 20,000 full-length cDNA sequences revealed problems with this approach, including inaccurate and inconsistent assignment of gene names, as well as many good assignments that were difficult to reproduce using only computational methods. For the FANTOM2 annotation of more than 60,000 cDNA clones, we developed a number of methods and tools to circumvent some of these problems, including an automated annotation pipeline that provides high-quality preliminary annotation for each sequence by introducing an "uninformative filter" that eliminates uninformative annotations, controlled vocabularies to accurately reflect both the functional assignments and the evidence supporting them, and a highly refined, Web-based manual annotation tool that allows users to view a wide array of sequence analyses and to assign gene names and putative functions using a consistent nomenclature. The ultimate utility of our approach is reflected in the low rate of reassignment of automated assignments by manual curation. Based on these results, we propose a new standard for large-scale annotation, in which the initial automated annotations are manually investigated and then computational methods are iteratively modified and improved based on the results of manual curation.
Animals↗
Comparative analyses of potato expressed sequence tag libraries.
The cultivated potato (Solanum tuberosum) shares similar biology with other members of the Solanaceae, yet has features unique within the family, such as modified stems (stolons) that develop into edible tubers. To better understand potato biology, we have undertaken a survey of the potato transcriptome using expressed sequence tags (ESTs) from diverse tissues. A total of 61,940 ESTs were generated from aerial tissues, below-ground tissues, and tissues challenged with the late-blight pathogen (Phytophthora infestans). Clustering and assembly of these ESTs resulted in a total of 19,892 unique sequences with 8,741 tentative consensus sequences and 11,151 singleton ESTs. We were able to identify a putative function for 43.7% of these sequences. A number of sequences (48) were expressed throughout the libraries sampled, representing constitutively expressed sequences. Other sequences (13,068, 21%) were uniquely expressed and were detected only in a single library. Using hierarchal and k means clustering of the EST sequences, we were able to correlate changes in gene expression with major physiological events in potato biology. Using pair-wise comparisons of tuber-related tissues, we were able to associate genes with tuber initiation, dormancy, and sprouting. We also were able to identify a number of characterized as well as novel sequences that were unique to the incompatible interaction of late-blight pathogen, thereby providing a foundation for further understanding the mechanism of resistance.
Cluster Analysis↗
The genetics of innate immunity in the lung.
Explore the source record for details and available documents.
Animals↗
The Schistosoma mansoni gene index: gene discovery and biology by reconstruction and analysis of expressed gene sequences.
Expressed sequence tag (EST) sequencing and analysis is a primary research tool to identify and characterize the Schistosoma mansoni transcriptome. As part of our gene discovery effort, a total of 5,793 ESTs have been generated from clones selected randomly from complementary DNA (cDNA) libraries constructed from male and female adult worms. Assembly analysis of all the 16,813 public S. mansoni ESTs has identified 1,920 distinct tentative consensus sequences (TCs) and 5,571 nonoverlapping ESTs (singletons). Of these, 376 TCs (20%) and 1,449 singletons (26%) are unique to the SUNY/TIGR sequencing effort. Tentative consensus sequences and singletons were distributed into various categories of biological roles associated with cell structure, metabolism, protein fate, signal transduction, transcription, protein synthesis, transporters, and cell growth. The TCs and singletons represent transcripts that can be used as a resource for functional annotation of genomic sequence data, comparative sequence analysis, and cDNA clone selection for microarray projects. The utility of EST analysis is demonstrated by identifying new protease genes, which may be involved in hemoglobin degradation.
Amino Acid Sequence↗
Open source software for the analysis of microarray data.
DNA microarray assays represent the first widely used application that attempts to build upon the information provided by genome projects in the study of biological questions. One of the greatest challenges with working with microarrays is collecting, managing, and analyzing data. Although several commercial and noncommercial solutions exist, there is a growing body of freely available, open source software that allows users to analyze data using a host of existing techniques and to develop their own and integrate them within the system. Here we review three of the most widely used and comprehensive systems, the statistical analysis tools written in R through the Bioconductor project (http://www.bioconductor.org), the Java-based TM4 software system available from The Institute for Genomic Research (http://www.tigr.org/software), and BASE, the Web-based system developed at Lund University (http://base.thep.lu.se).
Database Management Systems↗
Within the fold: assessing differential expression measures and reproducibility in microarray assays.
BACKGROUND: 'Fold-change' cutoffs have been widely used in microarray assays to identify genes that are differentially expressed between query and reference samples. More accurate measures of differential expression and effective data-normalization strategies are required to identify high-confidence sets of genes with biologically meaningful changes in transcription. Further, the analysis of a large number of expression profiles is facilitated by a common reference sample, the construction of which must be carefully addressed. RESULTS: We carried out a series of 'self-self' hybridizations in which aliquots of the same RNA sample were labeled separately with Cy3 and Cy5 fluorescent dyes and co-hybridized to the same microarray. From this, we can analyze the intensity-dependent behavior of microarray data, define a statistically significant measure of differential expression that exploits the structure of the fluorescent signals, and measure the inherent reproducibility of the technique. We also devised a simple procedure for identifying and eliminating low-quality data for replicates within and between slides. We examine the properties required of a universal reference RNA sample and show how pooling a small number of samples with a diverse representation of expressed genes can outperform more complex mixtures as a reference sample. CONCLUSION: Analysis of cell-line samples can identify systematic structure in measured gene-expression levels. A general procedure for analyzing cDNA microarray data is proposed and validated. We show that pooled reference samples should be based not only on the expression of individual genes in each cell line but also on the expression levels of genes within cell lines.
Brain Neoplasms↗
Standards for microarray data.
Explore the source record for details and available documents.
Computational Biology↗
Identification of Src transformation fingerprint in human colon cancer.
We used a classical rodent model of transformation to understand the transcriptional processes, and hence the molecular and cellular events a given cell undergoes when progressing from a normal to a transformed phenotype. Src activation is evident in 80% of human colon cancer, yet the myriad of cellular processes effected at the level of gene expression has yet to be fully documented. We identified a Src 'transformation fingerprint' within the gene expression profiles of Src-transformed rat 3Y1 fibroblasts demonstrating a progression in transformation characteristics. To evaluate the role of this gene set in human cancer development and progression, we extracted the orthologous genes present on the Affymetrix Hu95A GeneChip (12k named genes) and compared expression profiles between the Src-induced rodent cell line model of transformation and staged colon tumors where Src is known to be activated. A similar gene expression pattern between the cell line model and staged colon tumors for components of the cell cycle, cytoskeletal associated proteins, transcription factors and lysosomal proteins suggests the need for co-regulation of several cellular processes in the progression of cancer. Genes not previously implicated in tumorigenesis were detected, as well as a set of 14 novel, highly conserved genes with here-to-fore unknown function. These studies define a set of transformation associated genes whose up-regulation has implications for understanding Src mediated transformation and strengthens the role of Src in the development and progression of human colon cancer. Supportive Supplemental Data can be viewed at http://pga.tigr.org/PGApubs.shtml.
Animals↗