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.     

Endoplasmic reticulum stress reduces the export from the ER and alters the architecture of post-ER compartments

In eukaryotic cells several physiologic and pathologic conditions generate the accumulation of unfolded proteins in the endoplasmic reticulum (ER), leading to ER stress. To restore normal function, some ER transmembrane proteins sense the ER stress and activate coordinated signalling pathways collectively called the Unfolded Protein Response (UPR). Little is known on how the UPR relates to post-ER compartments and to the export from the ER of newly synthesized proteins. Here, we report that the ER stress response induced by either thapsigargin or nitric oxide modifies the dynamics of the intracellular distribution of ERGIC-53 and GM130, two markers of the ER Golgi Intermediate Compartment and of the cis-Golgi, respectively. In addition, induction of ER stress alters the morphology of the ERGIC and the Golgi complex and interferes with the reformation of both compartments. Moreover, ER stress rapidly reduces the transport to the Golgi complex of the temperature sensitive mutant of the Vesicular Stomatitis Virus G Glycoprotein (VSV-G) fused with the Green Fluorescent Protein (ts045G), without apparently decreasing the amount of the protein competent for export. Interestingly, a parallel rapid reduction of the number of Sec31 labelled fluorescent puncta on the ER membranes does occur, thus suggesting that the ER stress alters the ER export and the dynamic of post-ER compartments by rapidly targeting the formation of COPII-coated transport intermediates.     

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.     

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.     

Radiographic joint damage in rheumatoid arthritis is associated with differences in cartilage turnover and can be predicted by serum biomarkers: an evaluation from 1 to 4 years after diagnosis

Introduction  

The objective of this study was to determine whether serum biomarkers for degradation and synthesis of the extracellular matrix of cartilage are associated with, and can predict, radiographic damage in patients with rheumatoid arthritis (RA).     

The Urokinase Receptor Takes Control of Cell Migration by Recruiting Integrins and FPR1 on the Cell Surface

The receptor (uPAR) of the urokinase-type plasminogen activator (uPA) is crucial in cell migration since it concentrates uPA proteolytic activity at the cell surface, binds vitronectin and associates to integrins. uPAR cross-talk with receptors for the formylated peptide fMLF (fMLF-Rs) has been reported; however, cell-surface uPAR association to fMLF-Rs on the cell membrane has never been explored in detail.We now show that uPAR co-localizes at the cell-surface and co-immunoprecipitates with the high-affinity fMLF-R, FPR1, in uPAR-transfected HEK-293 (uPAR-293) cells. uPAR/β1 integrin and FPR1/β1 integrin co-localization was also observed. Serum or the WKYMVm peptide (W Pep), a FPR1 ligand, strongly increased all observed co-localizations in uPAR-293 cells, including FPR1/β1 integrin co-localization. By contrast, a low FPR1/β1 integrin co-localization was observed in uPAR-negative vector-transfected HEK-293 (V-293) cells, that was not increased by serum or W Pep stimulations.The role of uPAR interactions in cell migration was then explored. Both uPAR-293 and V-293 control cells efficiently migrated toward serum or purified EGF. However, cell treatments impairing uPAR interactions with fMLF-Rs or integrins, or inhibiting specific cell-signaling mediators abrogated uPAR-293 cell migration, without exerting any effect on V-293 control cells.Accordingly, uPAR depletion by a uPAR-targeting siRNA or uPAR blocking with an anti-uPAR polyclonal antibody in cells constitutively expressing high uPAR levels totally impaired their migration toward serum.Altogether, these results suggest that both uPAR-positive and uPAR-negative cells are able to migrate toward serum; however, uPAR expression renders cell migration totally and irreversibly uPAR-dependent, since it is completely inhibited by uPAR blocking.We propose that uPAR takes control of cell migration by recruiting fMLF-Rs and β1 integrins, thus promoting their co-localization at the cell-surface and driving pro-migratory signaling pathways.     

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.