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.      

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.     

Biomarkers for screening of lung cancer and pre-neoplastic lesions in a high risk Chilean population

Background

The mortality of lung cancer (LC), increases each year in the world, in spite of any advances, in development of new drugs to advance stages of LC. The high incidence of LC has been associated with smoking habit, genetic diversity and environmental pollution. Antofagasta region has been reported to have the highest LC mortality rate in Chile and its inhabitants were exposed to arsenic in their drinking water in concentrations as high as 870 μg/L. Non-invasive techniques such as biomarkers (Automatic Quantitative Cytometry: AQC and DR70) and Auto Fluorescence Bronchoscopy (AFB) might be potentially useful as a supplementary diagnostic approach and early detection. Early detection is one of the most important factors to intervene and prevent cancer progression in LC. This is a work of an ongoing prospective bimodality cancer surveillance study in high risk LC volunteers. Enrolment was done in subjects from Antofagasta and Metropolitan regions. In addition, we enrolled subjects who were suspected of having lung cancer. AQC, DR70 and AFB were used as tools in the detection of pre-neoplastic (PNL) and neoplastic lesions (NL).

Results

Half of the samples, classified as suspicious by AFB, were confirmed as metaplasia or dysplasia by histopathology. For LC, DR70 showed a higher sensitivity (95.8%) and specificity (91.9%) than AQC. However, for PNL AQC showed a higher sensitivity (91.9%) than DR70 (27.3%), although both with low PPV values. As a pre screener, both biomarkers might be employed as complementary tools to detect LC, especially as serially combined tests, with a sensitivity of 60% and a PPV of 65.2%. Additionally, the use of parallel combined tests might support the detection of PNL (sensitivity 91.2%; PPV 49.1%).

Conclusion

This work adds information on cellular and molecular biomarkers to complement imaging techniques for early detection of LC in Latin America that might contribute to formulate policies concerning screening of LC. Supported by INNOVA-CORFO, Chile.     

Connexin- and Pannexin-Based Channels in Normal Skeletal Muscles and Their Possible Role in Muscle Atrophy

Precursor cells of skeletal muscles express connexins 39, 43 and 45 and pannexin1. In these cells, most connexins form two types of membrane channels, gap junction channels and hemichannels, whereas pannexin1 forms only hemichannels. All these channels are low-resistance pathways permeable to ions and small molecules that coordinate developmental events. During late stages of skeletal muscle differentiation, myofibers become innervated and stop expressing connexins but still express pannexin1 hemichannels that are potential pathways for the ATP release required for potentiation of the contraction response. Adult injured muscles undergo regeneration, and connexins are reexpressed and form membrane channels. In vivo, connexin reexpression occurs in undifferentiated cells that form new myofibers, favoring the healing process of injured muscle. However, differentiated myofibers maintained in culture for 48 h or treated with proinflammatory cytokines for less than 3 h also reexpress connexins and only form functional hemichannels at the cell surface. We propose that opening of these hemichannels contributes to drastic changes in electrochemical gradients, including reduction of membrane potential, increases in intracellular free Ca²⁺ concentration and release of diverse metabolites (e.g., NAD⁺ and ATP) to the extracellular milieu, contributing to multiple metabolic and physiologic alterations that characterize muscles undergoing atrophy in several acquired and genetic human diseases. Consequently, inhibition of connexin hemichannels expressed by injured or denervated skeletal muscles might reduce or prevent deleterious changes triggered by conditions that promote muscle atrophy.     

Simple and reliable multiplex PCR assay for detection of blaTEM,bla(SHV) and blaOXA-1 genes in Enterobacteriaceae

Third-generation cephalosporin resistance is often mediated by TEM- and SHV-type beta-lactamases in Enterobacteriaceae. TEM-type and OXA-1 enzymes are the major plasmid-borne beta-lactamases implicated in amoxicillin-clavulanic acid resistance in Escherichia coli isolates. We have developed a rapid and simple multiplex polymerase chain reaction (PCR) which discriminates bla(TEM), bla(SHV) and bla(OXA-1) genes by generating fragments of 516, 392 and 619 bp respectively. Multiplex PCR analysis of 51 amoxicillin-clavulanate resistant E. coli isolates detected bla(TEM) and bla(SHV) genes in 45 and two strains, respectively, and only one strain harboured a bla(OXA-1) gene. Twenty-three of the 40 cefotaxime-resistant Enterobacteriaceae isolates produced amplicons with a size compatible with the presence of bla(TEM) (13 strains), bla(SHV) (six strains) genes or the association of both genes (four strains). These results were verified by colony hybridisation. Therefore, multiplex PCR is a suitable tool for initial rapid screening of bla genes in Enterobacteriaceae.     

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.