Common pitfalls in bioinformatics-based analyses: look before you leap

A Targeted Bioinformatics Assessment of Adrenocortical Carcinoma Reveals Prognostic Implications of GABA System Gene Expression

Adrenocortical carcinoma (ACC) is a rare but deadly cancer for which few treatments exist. Here, we have undertaken a targeted bioinformatics study of The Cancer Genome Atlas (TCGA) ACC dataset focusing on the 30 genes encoding the γ-aminobutyric acid (GABA) system—an under-studied, evolutionarily-conserved system that is an emerging potential player in cancer progression. Our analysis identified a subset of ACC patients whose tumors expressed a distinct GABA system transcriptome. Transcript levels of ABAT (encoding a key GABA shunt enzyme), were upregulated in over 40% of tumors, and this correlated with several favorable clinical outcomes including patient survival; while enrichment and ontology analysis implicated two cancer-related biological pathways involved in metastasis and immune response. The phenotype associated with ABAT upregulation revealed a potential metabolic heterogeneity among ACC tumors associated with enhanced mitochondrial metabolism. Furthermore, many GABA_A receptor subunit-encoding transcripts were expressed, including two (GABRB2 and GABRD) prognostic for patient survival. Transcripts encoding GABA_B receptor subunits and GABA transporters were also ubiquitously expressed. The GABA system transcriptome of ACC tumors is largely mirrored in the ACC NCI-H295R cell line, suggesting that this cell line may be appropriate for future functional studies investigating the role of the GABA system in ACC cell growth phenotypes and metabolism.

Patterns of selection and tissue-specific expression among maize domestication and crop improvement loci

The domestication of maize (Zea mays sp. mays) from its wild progenitors represents an opportunity to investigate the timing and genetic basis of morphological divergence resulting from artificial selection on target genes. We compared sequence diversity of 30 candidate selected and 15 reference loci between the three populations of wild teosintes, maize landraces, and maize inbred lines. We inferred an approximately equal ratio of genes selected during early domestication and genes selected during modern crop breeding. Using an expanded dataset of 48 candidate selected and 658 neutral reference loci, we tested the hypothesis that candidate selected genes in maize are more likely to have transcriptional functions than neutral reference genes, but there was no overrepresentation of regulatory genes in the selected gene dataset. Electronic northern analysis revealed that candidate genes are significantly overexpressed in the maize ear relative to vegetative tissues such as maize shoot, leaf, and root tissue. The maize ear underwent dramatic morphological alteration upon domestication and has been a continuing target of selection for maize yield. Therefore, we hypothesize that genes targeted by selection are more likely to be expressed in tissues that experienced high levels of morphological divergence during domestication and crop improvement.

Application of DNA array technology to mammalian embryos

Early embryogenesis depends on a tightly choreographed succession of gene expression patterns which define normal development. Fertilization and the first zygotic cleavage involve major changes to paternal and maternal chromatin and translation of maternal RNAs which have been sequestered in the oocyte during oogenesis. At a critical species-specific point known as the major onset of embryonic expression, there is a dramatic increase in expression from the new diploid genome. The advent of array technology has, for the first time, made possible to determine the transcriptional profile of all approximately 20,000 mammalian genes during embryogenesis, although the small amount of mRNA in a single embryo necessitates either pooling large numbers of embryos or a global amplification procedure to give sufficient labeled RNA for analysis. Following array hybridization, various bioinformatic tools must be employed to determine the expression level for each gene, often based on multiple oligonucleotide probes and complex background estimation protocols. The grouped analysis of clusters of genes which represent specific biological pathways provides the key to understanding embryonic development, embryonic stem cell proliferation and the reprogramming of gene expression after somatic cloning. Arrays are being developed to address specific biological questions related to embryonic development including DNA methylation and microRNA expression. Array technology in its various facets is an important diagnostic tool for the early detection of developmental aberrations; for improving the safety of assisted reproduction technologies for man; and for improving the efficiency of producing cloned and/or transgenic farm animals. This review discusses current approaches and limitations of DNA microarray technology with emphasis on bovine embryos.

Validation of in silico prediction by in vitro immunoserological results of fine epitope mapping on citrate synthase specific autoantibodies

In silico antibody-antigen binding predictions are generally employed in research to rationalize epitope development. These techniques are widely spread despite their technical limitations. To validate the results of these bioinformatic calculations evidence based comparative in vitro studies are necessary. We have used a well-conserved mitochondrial inner membrane antigen-citrate synthase to develop a model for comparative analysis of the predicted and the immunoserologically verified epitopes of circulating autoantibodies. Epitopes were predicted using accepted tools: the GCG Wisconsin package and TEPITOPE 2000. An overlapping multipin ELISA assay--covering 49% of the citrate synthase molecule--was developed to map autoantibody epitopes of individuals (healthy, systemic autoimmune, and heart transplanted) in different immunopathological conditions. From the 40 synthesized decapeptides 34 were predicted in silico and 27 were validated in vitro. Thirty-two percent of epitopes were recognized by majority of sera 47% by at least one sera. False positive predictions were 21%. There was major difference in the recognized epitope pattern under different immunopathological conditions. Our results suggest that special databases are needed for training and weighing prediction methods by clinically well-characterized samples, due to the differences in the immune response under different health status. The development of these special algorithms needs a new approach. A high number of samples under these special immunological conditions are to be mapped and then used for the "fine tuning" of different prediction algorithms.

Cloning of a novel signaling molecule, AMSH-2, that potentiates transforming growth factor beta signaling

BACKGROUND

Transforming growth factor-betas (TGF-betas), bone morphogenetic proteins (BMPs) and activins are important regulators of developmental cell growth and differentiation. Signaling by these factors is mediated chiefly by the Smad family of latent transcription factors.

RESULTS

There are a large number of uncharacterized cDNA clones that code for novel proteins with homology to known signaling molecules. We have identified a novel molecule from the HUGE database that is related to a previously known molecule, AMSH (associated molecule with the SH3 domain of STAM), an adapter shown to be involved in BMP signaling. Both of these molecules contain a coiled-coil domain located within the amino-terminus region and a JAB (Domain in Jun kinase activation domain binding protein and proteasomal subunits) domain at the carboxy-terminus. We show that this novel molecule, which we have designated AMSH-2, is widely expressed and its overexpression potentiates activation of TGF-beta-dependent promoters. Coimmunoprecipitation studies indicated that Smad7 and Smad2, but not Smad3 or 4, interact with AMSH-2. We show that overexpression of AMSH-2 decreases the inhibitory effect of Smad7 on TGF-beta signaling. Finally, we demonstrate that knocking down AMSH-2 expression by RNA interference decreases the activation of 3TP-lux reporter in response to TGF-beta.

CONCLUSIONS

This report implicates AMSH and AMSH-2 as a novel family of molecules that positively regulate the TGF-beta signaling pathway. Our results suggest that this effect could be partially explained by AMSH-2 mediated decrease of the action of Smad7 on TGF-beta signaling pathway.

From biological databases to platforms for biomedical discovery

The use of high-throughput DNA sequencing and proteomic methods has led to an unprecedented increase in the amount of genomic and proteomic data. Application of computing technologies and development of computational tools to analyze and present these data has not kept pace with the accumulation of information. Here, we discuss the use of different database systems to store biological information and mention some of the key emerging computing technologies that are likely to have a key role in the future of bioinformatics.

A reassessment of the translation initiation codon in vertebrates

More than two decades ago Marilyn Kozak proposed the scanning model of translation initiation, whereby translation is initiated at the first AUG codon that is in a particular context. In this article, we re-examine the context of initiator codons using a large dataset of curated human transcripts. We find that more than 40% of transcripts contain AUG codons upstream of the actual start codon and that most authentic AUGs contain three or more mismatches from the consensus sequence, CCACCaugG. Also, in a large fraction of transcripts, the sequences surrounding the initiator codon deviate more from the consensus than those surrounding upstream AUGs, indicating that translation initiation from downstream AUGs is more common than generally believed.