An investigation of biomarkers derived from legacy microarray data for their utility in the RNA-seq era

Background

Gene expression microarray has been the primary biomarker platform ubiquitously applied in biomedical research, resulting in enormous data, predictive models, and biomarkers accrued. Recently, RNA-seq has looked likely to replace microarrays, but there will be a period where both technologies co-exist. This raises two important questions: Can microarray-based models and biomarkers be directly applied to RNA-seq data? Can future RNA-seq-based predictive models and biomarkers be applied to microarray data to leverage past investment?

Results

We systematically evaluated the transferability of predictive models and signature genes between microarray and RNA-seq using two large clinical data sets. The complexity of cross-platform sequence correspondence was considered in the analysis and examined using three human and two rat data sets, and three levels of mapping complexity were revealed. Three algorithms representing different modeling complexity were applied to the three levels of mappings for each of the eight binary endpoints and Cox regression was used to model survival times with expression data. In total, 240,096 predictive models were examined.

Conclusions

Signature genes of predictive models are reciprocally transferable between microarray and RNA-seq data for model development, and microarray-based models can accurately predict RNA-seq-profiled samples; while RNA-seq-based models are less accurate in predicting microarray-profiled samples and are affected both by the choice of modeling algorithm and the gene mapping complexity. The results suggest continued usefulness of legacy microarray data and established microarray biomarkers and predictive models in the forthcoming RNA-seq era.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0523-y) contains supplementary material, which is available to authorized users.     

Regeneration in the auditory system of nymphal and adult bush crickets Tettigonia viridissima

Regeneration and reestablishment of synaptic connections is an important topic in neurobiological research. In the present study, the regeneration of auditory afferents and the accompanying effects in the central nervous system are investigated in nymphs and adults of the bush cricket Tettigonia viridissima L. (Orthoptera: Tettigoniidae). In all animals in which the tympanal nerve is crushed, neuronal tracing shows a regrowth of the afferents into the prothoracic ganglion. This regeneration is seen in both adult and nymphal stages and starts 10–15 days after nerve crushing. Physiological recordings from the leg nerve indicate a recovery of tympanal fibres and a formation of functional connections to interneurones in the same time range. Electrophysiological recordings from the neck connective suggest additional contralateral sprouting of interneurones and the formation of aberrant connections. The regeneration processes of the tympanal nerve in nymphal stages and adults appear to be similar.     

Human cytochrome P-450 PB-1: A multigene family involved in mephenytoin and steroid oxidations that maps to chromosome 10

The cytochrome P-450 monooxygenase system possesses catalytic activity toward many exogenous compounds (e.g., drugs, insecticides, and polycyclic aromatic hydrocarbons) and endogenous compounds (e.g., steroids, fatty acids, and prostaglandins). Multiple forms of cytochrome P-450 with different substrate specificities have been isolated. In the present paper we report the isolation and sequence of a cDNA clone for the human hepatic cytochrome P-450 responsible for mephenytoin (an anticonvulsant) oxidation. The mephenytoin cytochrome P-450 is analogous to the rat cytochrome P-450 form termed PB-1 (family P450C2C). We also report that human PB-1 is encoded by one of a small family of related genes all of which map to human chromosome 10q24.1-10q24.3. The endogenous role of this enzyme appears to be in steroid oxidations. This cytochrome P-450 family does not correspond to any of the hepatic cytochrome P-450 gene families previously mapped in humans.     

Phylogenetic relationships of dasyuromorphian marsupials revisited

We reassessed the phylogenetic relationships of dasyuromorphians using a large molecular database comprising previously published and new sequences for both nuclear (nDNA) and mitochondrial (mtDNA) genes from the numbat (Myrmecobius fasciatus), most living species of Dasyuridae, and the recently extinct marsupial wolf, Thylacinus cynocephalus. Our molecular tree suggests that Thylacinidae is sister to Myrmecobiidae + Dasyuridae. We show robust support for the dasyurid intrafamilial classification proposed by Krajewski & Westerman as well as for placement of most dasyurid genera, which suggests substantial homoplasy amongst craniodental characters presently used to generate morphology‐based taxonomies. Molecular dating with relaxed molecular clocks suggests that dasyuromorphian cladogenesis began in the Eocene, and that all three dasyuromorphian families originated prior to the end of this epoch. Radiation within Thylacinidae and Dasyuridae had occurred by the middle to late Oligocene, consistent with recognition of primitive thylacinids (e.g. Badjcinus turnbulli) in the later Oligocene and of putative dasyurids (e.g. Barinya wangala) by the early Miocene. We propose that all four extant dasyurid tribes were in existence by the early Miocene and that most modern dasyurid genera/species were established before the later Miocene. This is in marked contrast to the popularly accepted advocation of their origins in the latest Miocene–early Pliocene. © 2015 The Linnean Society of London     

Earth Mover’s Distance (EMD): A True Metric for Comparing Biomarker Expression Levels in Cell Populations

Changes in the frequencies of cell subsets that (co)express characteristic biomarkers, or levels of the biomarkers on the subsets, are widely used as indices of drug response, disease prognosis, stem cell reconstitution, etc. However, although the currently available computational “gating” tools accurately reveal subset frequencies and marker expression levels, they fail to enable statistically reliable judgements as to whether these frequencies and expression levels differ significantly between/among subject groups. Here we introduce flow cytometry data analysis pipeline which includes the Earth Mover’s Distance (EMD) metric as solution to this problem. Well known as an informative quantitative measure of differences between distributions, we present three exemplary studies showing that EMD 1) reveals clinically-relevant shifts in two markers on blood basophils responding to an offending allergen; 2) shows that ablative tumor radiation induces significant changes in the murine colon cancer tumor microenvironment; and, 3) ranks immunological differences in mouse peritoneal cavity cells harvested from three genetically distinct mouse strains.     

Treatment of Premature Labour with Orciprenaline

The use of orciprenaline, a β-adrenergic agent, in 30 patients with probable premature labour proved to be effective and safe to both mother and fetus under the conditions employed.     

Cultivator for NMR studies of suspended cell cultures

When nuclear magnetic resonance (NMR) spectroscopy is employed for physiological experiments with suspended cells, providing for adequate nutrient and oxygen delivery is particularly important, because the inherent insensitivity of NMR requires that concentrated cell suspensions be used. In addition, it is desirable to be able to manipulate the growth rate of cells during a NMR experiment. To address these concerns, a continuous cell cultivator that provides convective oxygen and nutrient transport has been constructed for NMR experiments. The NMR detector coil is located within the cultivator volume. The location is advantageous because the rapid exchange of cells in and out of the coil leads to a small apparent spin lattice relaxation time, thus allowing for rapid pulsing and fast signal averaging. In this article we present the physical principles on which the cultivator's design is based. (31)P spectra showing the response of continuously cultivated Saccharomyces cerevisiae cultures to a phosphate bolus and growth rate shift are then given. (c) 1992 John Wiley & Sons, Inc.     

Oocysts of Isospora ernsti n. sp. and Isospora blagburni n. sp. (Apicomplexa: Eimeriidae) from the black-capped bulbul,Pycnonotus xanthopygos

Oocysts of Isospora ernsti n. sp. and Isospora blagburni n. sp. are described from the black-capped bulbul Pycnonotus xanthopygos from Lincoln Park Zoo, Chicago, Illinois. The bird came from southwestern Africa seven years earlier. I. ernsti oocysts are ellipsoidal to bluntly ovoid, 28–38 × 23–31m (mean 34 × 28 m) and have a single-layered oocyst wall. Micropyle, oocyst residuum and polar granules are absent. Sporocysts are elongate ovoid, 24–30 × 11–16 m (mean 27×13 m). Stieda and substiedal bodies and sporocyst residuum are present. I. blagburni oocysts are spherical to subspherical. 21–28 × 19–26 m (mean 25 × 23 m) and have a single oocyst wall. Sporocysts are ovoid and 17–23 × 10–13 m (mean 20 × 12 m). Stieda and substiedal bodies and sporocyst residuum are present.     

Polycyclic aromatic hydrocarbons physically intercalate into duplex regions of denatured DNA

We have investigated the physical binding of pyrene and benzo[a]pyrene derivatives to denatured DNA. These compounds exhibit a red shift in their absorbance spectra of 9 nm when bound to denatured calf thymus DNA, compared to a shift of 10 nm when binding occurs to native DNA. Fluorescence from the hydrocarbons is severely quenched when bound to both native and denatured DNA. Increasing sodium ion concentration decreases binding of neutral polycyclic aromatic hydrocarbons to native DNA and increases binding to denatured DNA. The direct relationship between binding to denatured DNA and salt concentration appears to be a general property of neutral polycyclic aromatic hydrocarbons. Absorption measurements at 260 nm were used to determine the duplex content of denatured DNA. When calculated on the basis of duplex binding sites, equilibrium constants for binding of 7,8,9,10-tetrahydroxy-7,8,9,10-tetrahydro-benzo[a]pyrene to denatured DNA are an order of magnitude larger than for binding to native DNA. The effect of salt on the binding constant was used to calculate the sodium ion release per bound ligand, which was 0.36 for both native and denatured DNA. Increasing salt concentration increases the duplex content of denatured DNA, and it appears that physical binding of polycyclic aromatic hydrocarbons consists of intercalation into these sites.