Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq

Background

Massively parallel cDNA sequencing (RNA-seq) experiments are gradually superseding microarrays in quantitative gene expression profiling. However, many biologists are uncertain about the choice of differentially expressed gene (DEG) analysis methods and the validity of cost-saving sample pooling strategies for their RNA-seq experiments. Hence, we performed experimental validation of DEGs identified by Cuffdiff2, edgeR, DESeq2 and Two-stage Poisson Model (TSPM) in a RNA-seq experiment involving mice amygdalae micro-punches, using high-throughput qPCR on independent biological replicate samples. Moreover, we sequenced RNA-pools and compared their results with sequencing corresponding individual RNA samples.

Results

False-positivity rate of Cuffdiff2 and false-negativity rates of DESeq2 and TSPM were high. Among the four investigated DEG analysis methods, sensitivity and specificity of edgeR was relatively high. We documented the pooling bias and that the DEGs identified in pooled samples suffered low positive predictive values.

Conclusions

Our results highlighted the need for combined use of more sensitive DEG analysis methods and high-throughput validation of identified DEGs in future RNA-seq experiments. They indicated limited utility of sample pooling strategies for RNA-seq in similar setups and supported increasing the number of biological replicate samples.

Electronic supplementary material

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

Characterization of the flavin association in hexose oxidase from Chondrus crispus

Hexose oxidase (EC 1.1.3.5) from Hansenula polymorpha was found to exhibit a dual covalent association of FAD with His79 via an 8 alpha-histidyl linkage as well as a covalent association between Cys138 and C-6 of the isoalloxazine moiety of FAD. Spectral properties of the wild-type enzyme exhibited maxima at 364 nm and 437 nm as well as a distinct shoulder at 445 nm. An H79K mutant enzyme exhibited only one maximum at 437 nm. The difference absorption spectrum between an oxidized and a substrate-reduced enzyme preparation showed maxima at 360 nm and 445 nm corresponding to an apparent novel type of association. Hexose oxidase showed a low, pH-independent fluorescence at 525 nm when excited at 450 nm. Flavin was released from the holoenzyme by treatment with trypsin. Sequencing of the flavopeptide revealed two peptides comprising positions 74-91 and 132-157 associated with FAD in equimolar amounts. A homology model of hexose oxidase was constructed using the crystal structure of glucooligosaccharide oxidase from Acremonium strictum as template. The model placed both of the sequences found above in the close vicinity of the FAD cofactor, and suggests covalent bonds between both His79 and Cys138 and FAD, in accordance with the chemical evidence. Based on the results, hexose oxidase is identified as incorporating FAD with a double covalent association with His79 and Cys138 in the holoenzyme. A reaction mechanism involving the concerted action of Tyr488 and Asp409 in hexose oxidase is suggested as the initiator of the proton abstraction from the substrate molecule in the active site.     

Biophysical consequences of remodeling the neutral side chains of rhamnogalacturonan I in tubers of transgenic potatoes

Two lines of transgenic potato (Solanum tuberosum L.) plants modified in their cell wall structure were characterized and compared to wild type with regard to biomechanical properties in order to assign functional roles to the particular cell wall polysaccharides that were targeted by the genetic changes. The targeted polymer was rhamnogalacturonan I (RG-I), a complex pectic polysaccharide comprised of mainly neutral oligosaccharide side chains attached to a backbone of alternating rhamnosyl and galacturonosyl units. Tuber rhamnogalacturonan I molecules from the two transformed lines are reduced in linear galactans and branched arabinans, respectively. The transformed tuber tissues were found to be more brittle when subjected to uniaxial compression and the side-chain truncation was found to be correlated with the physical properties of the tissue. Interpretation of the force–deflection curves was aided by a mathematical model that describes the contribution of the cellulose microfibrils, and the results lead to the proposition that the pectic matrix plays a role in transmitting stresses to the load-bearing cellulose microfibrils and that even small changes to the rheological properties of the matrix have consequences for the biophysical properties of the wall.     

Virtual screening and bioassay study of novel inhibitors for dengue virus mRNA cap (nucleoside-2'O)-methyltransferase

We report high-throughput structure-based virtual screening of putative Flavivirus 2′-O-methyltransferase inhibitors together with results from subsequent bioassay tests of selected compounds. Potential inhibitors for the S-adenosylmethionine binding site were explored using 2D similarity searching, pharmacophore filtering and docking. The inhibitory activities of 15 top-ranking compounds from the docking calculations were tested on a recombinant methyltransferase with the RNA substrate ^7MeGpppAC₅. Local and global docking simulations were combined to estimate the ligand selectivity for the target site. The results of the combined computational and experimental screening identified a novel inhibitor, with a previously unknown scaffold, that has an IC₅₀ value of 60 μM.     

SORLA is required for insulin-induced expansion of the adipocyte precursor pool in visceral fat

Visceral adipose tissue shows remarkable plasticity, constantly replacing mature adipocytes from an inherent pool of adipocyte precursors. The number of precursors is set in the juvenile organism and remains constant in adult life. Which signals drive precursor pool expansion in juveniles and why they operate in visceral but not in subcutaneous white adipose tissue (WAT) are unclear. Using mouse models, we identified the insulin-sensitizing receptor SORLA as a molecular factor explaining the distinct proliferative capacity of visceral WAT. High levels of SORLA activity in precursors of juvenile visceral WAT prime these cells for nutritional stimuli provided through insulin, promoting mitotic expansion of the visceral precursor cell pool in overfed juvenile mice. SORLA activity is low in subcutaneous precursors, blunting their response to insulin and preventing diet-induced proliferation of this cell type. Our findings provide a molecular explanation for the unique proliferative properties of juvenile visceral WAT, and for the genetic association of SORLA with visceral obesity in humans.