Microarray data integration for genome-wide analysis of human tissue-selective gene expression

Background

Microarray gene expression data are accumulating in public databases. The expression profiles contain valuable information for understanding human gene expression patterns. However, the effective use of public microarray data requires integrating the expression profiles from heterogeneous sources.

Results

In this study, we have compiled a compendium of microarray expression profiles of various human tissue samples. The microarray raw data generated in different research laboratories have been obtained and combined into a single dataset after data normalization and transformation. To demonstrate the usefulness of the integrated microarray data for studying human gene expression patterns, we have analyzed the dataset to identify potential tissue-selective genes. A new method has been proposed for genome-wide identification of tissue-selective gene targets using both microarray intensity values and detection calls. The candidate genes for brain, liver and testis-selective expression have been examined, and the results suggest that our approach can select some interesting gene targets for further experimental studies.

Conclusion

A computational approach has been developed in this study for combining microarray expression profiles from heterogeneous sources. The integrated microarray data can be used to investigate tissue-selective expression patterns of human genes.

     

Orymold: ontology based gene expression data integration and analysis tool applied to rice

Background  

Integration and exploration of data obtained from genome wide monitoring technologies has become a major challenge for many bioinformaticists and biologists due to its heterogeneity and high dimensionality. A widely accepted approach to solve these issues has been the creation and use of controlled vocabularies (ontologies). Ontologies allow for the formalization of domain knowledge, which in turn enables generalization in the creation of querying interfaces as well as in the integration of heterogeneous data, providing both human and machine readable interfaces.     

Optimized between-group classification: a new jackknife-based gene selection procedure for genome-wide expression data

Background  

A recent publication described a supervised classification method for microarray data: Between Group Analysis (BGA). This method which is based on performing multivariate ordination of groups proved to be very efficient for both classification of samples into pre-defined groups and disease class prediction of new unknown samples. Classification and prediction with BGA are classically performed using the whole set of genes and no variable selection is required. We hypothesize that an optimized selection of highly discriminating genes might improve the prediction power of BGA.     

SplicerAV: a tool for mining microarray expression data for changes in RNA processing

Background  

Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets.     

ReplicationDomain: a visualization tool and comparative database for genome-wide replication timing data

Background  

Eukaryotic DNA replication is regulated at the level of large chromosomal domains (0.5–5 megabases in mammals) within which replicons are activated relatively synchronously. These domains replicate in a specific temporal order during S-phase and our genome-wide analyses of replication timing have demonstrated that this temporal order of domain replication is a stable property of specific cell types.     

Nonsense RNA: a tool for specifically inhibiting the expression of a gene in vivo

We describe a general technique to inhibit gene expression in eukaryotic cells. The gene we chose to inhibit was the E. coli LacZ gene (encoding beta-galactosidase), which has previously been cloned into a eukaryotic expression vector [1]. This plasmid is called pCH110. We constructed a variant of pCH110 in which we flipped a 2566 base pair 5' fragment of the LacZ gene into the antiparallel orientation. The plasmid containing this mutated LacZ gene is called pNSLacZ (NS signifies non-sense coding sequence). When equal amounts of pCH110 and pNSLacZ are co-transfected into 3T6 mouse fibroblasts, the beta-galactosidase activity is decreased by approximately a factor of ten. Increasing the ratio of pNSLacZ to pCH110 above 1:1 does not appreciably increase the level of inhibition. Next, we prove the specificity of the inhibition by adding a third gene to the transfection mixture. For this purpose, we used pSVneo beta, a plasmid which expresses a phosphotransferase. We found that even when the beta-galactosidase activity was diminished by a factor of 10, the phosphotransferase activity was unaffected. Therefore, we have demonstrated that: the presence of an antiparallel copy of the LacZ gene results in a significant and specific diminution of the LacZ gene's expression; only a fraction of the LacZ gene needs to be in the antiparallel orientation in order to observe this effect. These results suggest that this technique can serve as a tool to decrease the level of gene expression in order to study the function of specific genes, or as a therapeutic manoeuvre in the treatment of disorders of abnormal gene expression.     

STEM: a tool for the analysis of short time series gene expression data

Background  

Time series microarray experiments are widely used to study dynamical biological processes. Due to the cost of microarray experiments, and also in some cases the limited availability of biological material, about 80% of microarray time series experiments are short (3–8 time points). Previously short time series gene expression data has been mainly analyzed using more general gene expression analysis tools not designed for the unique challenges and opportunities inherent in short time series gene expression data.     

The ribosomal RNA gene cluster in aneuploid chickens: evidence for increased gene dosage and regulation of gene expression

In the chicken, the nucleolus organizer regions, or sites of the genes encoding 18S, 5.8S, and 28S ribosomal RNA (rRNA), map to one pair of microchromosomes that can be identified by silver nitrate cytochemistry. This nucleolar organizer chromosome also contains the major histocompatibility complex. Chickens aneuploid for this chromosome have been identified and reproduced for over seven generations. Crossing two trisomic parents results in the production of viable disomic, trisomic, and tetrasomic progeny, showing two, three, and four nucleoli and nucleolar organizers per cell, respectively. A molecular analysis of rRNA genes was undertaken to establish the gene copy numbers in the aneuploid genotypes, and to determine if elevated numbers of rRNA genes are stably maintained and inherited over multiple generations. Gene copy numbers were determined using hybridization analysis of erythrocyte DNA obtained from individuals comprising a family which segregated disomic, trisomic, and tetrasomic genotypes. The values obtained were 290, 420, and 570 rDNA repeats per cell for disomic, trisomic, and tetrasomic animals, respectively. These results provide molecular confirmation of the two aneuploid states and show that elevated gene copy numbers have been maintained over multiple generations. Fibroblasts derived from disomic and tetrasomic embryos were found to grow at similar rates in culture, and mature rRNA levels in chicken embryo fibroblasts from disomic, trisomic and tetrasomic embryos were also found to have similar levels of mature rRNA. Therefore, despite the increase in rDNA content, the level of rRNA is regulated to diploid amounts in aneuploid fibroblasts.     

Kleisli: a new tool for data integration in biology.

One of the central problems in bioinformatics is data retrieval and integration. The existing biological databases are geographically distributed across the Internet, complex and heterogeneous in data types and data structures, and constantly changing. With the current rapid growth of biomedical data, the challenge is how large volumes of data retrieved from multiple databases can be transformed and integrated automatically and flexibly. This article describes a powerful new tool, the Kleisli system, for complex queries across multiple databases and data integration.     

基因鉴定集成法:全基因组基因表达研究的新策略

人类基因组包含的核苷酸数目庞大,基因鉴定（识别）的技术策略是基因克隆研究至为重要的基础。在全基因组基因表达分析策略方面,已相继建立了ｍＲＮＡ差异显示、代表性差异分析、抑制性消减杂交、基因表达系列分析和ｃＤＮＡ微阵列等技术。基因鉴定集成法是新近在综合上述技术的优缺点的基础上建立的全基因组分析新策略,具有充分利用生物基因信息数据库进行基因鉴定（识别）,并能提高稀有拷贝基因鉴定效率的优点。本文简要介绍其     

SPI: a tool for incorporating gene expression data into a four-dimensional database of Caenorhabditis elegans embryogenesis

MOTIVATION: A comprehensive gene expression database is essential for computer modeling and simulation of biological phenomena, including development. Development is a four-dimensional (4D; 3D structure and time course) phenomenon. We are constructing a 4D database of gene expression for the early embryogenesis of the nematode Caenorhabditis elegans. As a framework of the 4D database, we have constructed computer graphics (CG), into which we will incorporate the expression data of a number of genes at the subcellular level. However, the assignment of 3D distribution of gene products (protein, mRNA), of embryos at various developmental stages, is both difficult and tedious. We need to automate this process. For this purpose, we developed a new system, named SPI after superimposing fluorescent confocal microscopic data onto a CG framework. RESULTS: The scheme of this system comprises the following: (1) acquirement of serial sections (40 slices) of fluorescent confocal images of three colors (4',6'-diamino-2-phenylindole (DAPI) for nuclei, indodicarbocyanine (Cy-3) for the internal marker, which is a germline-specific protein POS-1 and indocarbocyanine (Cy-5) for the gene product to be examined); (2) identification of several features of the stained embryos, such as contour, developmental stage and position of the internal marker; (3) selection of CG images of the corresponding stage for template matching; (4) superimposition of serial sections onto the CG; (5) assignment of the position of superimposed gene products. The Snakes algorithm identified the embryo contour. The detection accuracy of embryo contours was 92.1% when applied to 2- to 28-cell-stage embryos. The accuracy of the developmental stage prediction method was 81.2% for 2- to 8-cell-stage embryos. We manually judged only the later stage embryos because the accuracy for embryos at the later stages was unsatisfactory due to experimental noise effects. Finally, our system chose the optimal CG and performed the superposition and assignment of gene product distribution. We established an initial 4D gene expression database with 56 maternal gene products. AVAILABILITY: This system is available at http://anti.lab.nig.ac.jp/spi/ and http://anti.lab.nig.ac.jp/4ddb/     

Genephony: a knowledge management tool for genome-wide research

Background  

One of the consequences of the rapid and widespread adoption of high-throughput experimental technologies is an exponential increase of the amount of data produced by genome-wide experiments. Researchers increasingly need to handle very large volumes of heterogeneous data, including both the data generated by their own experiments and the data retrieved from publicly available repositories of genomic knowledge. Integration, exploration, manipulation and interpretation of data and information therefore need to become as automated as possible, since their scale and breadth are, in general, beyond the limits of what individual researchers and the basic data management tools in normal use can handle. This paper describes Genephony, a tool we are developing to address these challenges.     

GCTA: a tool for genome-wide complex trait analysis

For most human complex diseases and traits, SNPs identified by genome-wide association studies (GWAS) explain only a small fraction of the heritability. Here we report a user-friendly software tool called genome-wide complex trait analysis (GCTA), which was developed based on a method we recently developed to address the "missing heritability" problem. GCTA estimates the variance explained by all the SNPs on a chromosome or on the whole genome for a complex trait rather than testing the association of any particular SNP to the trait. We introduce GCTA's five main functions: data management, estimation of the genetic relationships from SNPs, mixed linear model analysis of variance explained by the SNPs, estimation of the linkage disequilibrium structure, and GWAS simulation. We focus on the function of estimating the variance explained by all the SNPs on the X chromosome and testing the hypotheses of dosage compensation. The GCTA software is a versatile tool to estimate and partition complex trait variation with large GWAS data sets.     

Molecular cloning into Tn5 and integration in the Pseudomonas aeruginosa chromosome: a tool for heterologous gene expression

The DNA primase gene of the promiscuous IncP-1 conjugative plasmid RP1, encoding two polypeptides of 118 and 80 kDa, was inserted into the transposon Tn5 in Escherichia coli. The derivative transposon, Tn2523, was then transposed to a temperature-sensitive replication mutant of the promiscuous IncP-1 conjugative plasmid R68 at permissive temperature and the plasmid transferred to Pseudomonas aeruginosa strain PAO. The latter strain was then grown at non-permissive temperature to identify transposition of Tn2523 into the P. aeruginosa chromosome. Immunological and enzymic analysis showed the expression of functional primase polypeptides in the constructed P. aeruginosa strain. This strain also restored wild-type conjugational transfer proficiency, by complementation, to mutants of the IncP-1 plasmid R18 affected in transfer from P. aeruginosa to P. stutzeri or to Acinetobacter calcoaceticus due to transposon Tn7 insertion mutations in the primase gene. This strategy of cloning into a transposon and integration into the bacterial chromosome should facilitate genetic manipulation and studies of gene expression in a range of Gram-negative bacteria.