Quality determination and the repair of poor quality spots in array experiments

Background  

A common feature of microarray experiments is the occurence of missing gene expression data. These missing values occur for a variety of reasons, in particular, because of the filtering of poor quality spots and the removal of undefined values when a logarithmic transformation is applied to negative background-corrected intensities. The efficiency and power of an analysis performed can be substantially reduced by having an incomplete matrix of gene intensities. Additionally, most statistical methods require a complete intensity matrix. Furthermore, biases may be introduced into analyses through missing information on some genes. Thus methods for appropriately replacing (imputing) missing data and/or weighting poor quality spots are required.     

Variations among cell lines in the synthesis of sphingolipids in de novo and recycling pathways

There are several pathways for the incorporation of sugars into glycosphingolipids (GSL). Sugars can be added to ceramide that contains sphinganine (dihydrosphingosine) synthesized de novo (pathway 1), to ceramide synthesized from sphingoid bases produced by hydrolysis of sphingolipids (pathway 2), and into GSL recycling from the endosomal pathway through the Golgi (pathway 3). We reported previously the surprising observation that SW13 cells, a human adrenal carcinoma cell line, synthesize most of their GSL in pathway 2. We now present data on the synthesis of GSL in four additional cell lines. Approximately 90% of sugar incorporation took place in pathway 2, and 10% or less in pathway 1, in human foreskin fibroblasts and NB41A3 neuroblastoma cells. In contrast, approximately 50-90% of sugar incorporation took place in pathway 1 in C2C12 myoblasts. The C2C12 cells divide more rapidly and synthesize 10-14 times as much GSL as the other three cell lines. In C6 glioma cells, approximately 30% of sugar incorporation occurred in pathway 1 and 60% in pathway 2. There was no relation between the utilization of pathways for GSL and sphingomyelin synthesis in foreskin fibroblasts and C2C12 cells. In both cells pathways 1 and 2 each accounted for 50% of incorporation of choline into sphingomyelin. In five of the six cell lines that we have studied, most GSL synthesis takes place in pathway 2. We suggest that when the need for synthesis is relatively low, as in slowly dividing cells, GSL are synthesized predominantly from sphingoid bases salvaged from the hydrolytic pathway. When cells are dividing more rapidly, the need for increased synthesis is met by upregulating the de novo pathway.      

Single-strand conformation polymorphism analysis using capillary array electrophoresis for large-scale mutation detection

This protocol describes capillary array electrophoresis single-strand conformation polymorphism (CAE-SSCP), a screening method for detection of unknown and previously identified mutations. The method detects 98% of mutations in a sample material and can be applied to any organism where the goal is to determine genetic variation. This protocol describes how to screen for mutations in 192 singleplex or up to 768 multiplex samples over 3 days. The protocol is based on the principle of sequence-specific mobility of single-stranded DNA in a native polymer, and covers all stages in the procedure, from initial DNA purification to final CAE-SSCP data analysis, as follows: DNA is purified, followed by PCR amplification using fluorescent primers. After PCR amplification, double-stranded DNA is heat-denatured to separate the strands and subsequently cooled on ice to avoid reannealing. Finally, samples are analyzed by capillary electrophoresis and appropriate analysis software.     

Non-native interactions play an effective role in protein folding dynamics

Systematic Monte Carlo simulations of simple lattice models show that the final stage of protein folding is an ordered process where native contacts get locked (i.e., the residues come into contact and remain in contact for the duration of the folding process) in a well‐defined order. The detailed study of the folding dynamics of protein‐like sequences designed as to exhibit different contact energy distributions, as well as different degrees of sequence optimization (i.e., participation of non‐native interactions in the folding process), reveals significant differences in the corresponding locking scenarios—the collection of native contacts and their average locking times, which are largely ascribable to the dynamics of non‐native contacts. Furthermore, strong evidence for a positive role played by non‐native contacts at an early folding stage was also found. Interestingly, for topologically simple target structures, a positive interplay between native and non‐native contacts is observed also toward the end of the folding process, suggesting that non‐native contacts may indeed affect the overall folding process. For target models exhibiting clear two‐state kinetics, the relation between the nucleation mechanism of folding and the locking scenario is investigated. Our results suggest that the stabilization of the folding transition state can be achieved through the establishment of a very small network of native contacts that are the first to lock during the folding process.     

Statistical modeling of biomedical corpora: mining the Caenorhabditis Genetic Center Bibliography for genes related to life span

Background  

The statistical modeling of biomedical corpora could yield integrated, coarse-to-fine views of biological phenomena that complement discoveries made from analysis of molecular sequence and profiling data. Here, the potential of such modeling is demonstrated by examining the 5,225 free-text items in the Caenorhabditis Genetic Center (CGC) Bibliography using techniques from statistical information retrieval. Items in the CGC biomedical text corpus were modeled using the Latent Dirichlet Allocation (LDA) model. LDA is a hierarchical Bayesian model which represents a document as a random mixture over latent topics; each topic is characterized by a distribution over words.     

The shape of human gene family phylogenies

Background  

The shape of phylogenetic trees has been used to make inferences about the evolutionary process by comparing the shapes of actual phylogenies with those expected under simple models of the speciation process. Previous studies have focused on speciation events, but gene duplication is another lineage splitting event, analogous to speciation, and gene loss or deletion is analogous to extinction. Measures of the shape of gene family phylogenies can thus be used to investigate the processes of gene duplication and loss. We make the first systematic attempt to use tree shape to study gene duplication using human gene phylogenies.     

A Taxonomic Search Engine: Federating taxonomic databases using web services

Background  

The taxonomic name of an organism is a key link between different databases that store information on that organism. However, in the absence of a single, comprehensive database of organism names, individual databases lack an easy means of checking the correctness of a name. Furthermore, the same organism may have more than one name, and the same name may apply to more than one organism.     

tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities

Background

To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative genomic studies.

Results

Here, we present a method that allows for the identification of differential gene expression associated with DNA copy number changes over time. The temporal variation in gene expression is described by a generalized linear mixed model employing low-rank thin-plate splines. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, posteriors of hyperparameters and model parameters are estimated. The empirical Bayes procedure shrinks multiple dispersion-related parameters. Shrinkage leads to more stable estimates of the model parameters, better control of false positives and improvement of reproducibility. In addition, to make estimates of the DNA copy number more stable, model parameters are also estimated in a multivariate way using triplets of features, imposing a spatial prior for the copy number effect.

Conclusion

With the proposed method for analysis of time-course multilevel molecular data, more profound insight may be gained through the identification of temporal differential expression induced by DNA copy number abnormalities. In particular, in the analysis of an integrative oncogenomics study with a time-course set-up our method finds genes previously reported to be involved in cervical carcinogenesis. Furthermore, the proposed method yields improvements in sensitivity, specificity and reproducibility compared to existing methods. Finally, the proposed method is able to handle count (RNAseq) data from time course experiments as is shown on a real data set.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-327) contains supplementary material, which is available to authorized users.