Global rank-invariant set normalization (GRSN) to reduce systematic distortions in microarray data

Predicting oncology drug-induced cardiotoxicity with donor-specific iPSC-CMs-a proof-of-concept study with doxorubicin

Many oncology drugs have been found to induce cardiotoxicity in a subset of patients, which significantly limits their clinical use and impedes the benefit of lifesaving anticancer treatments. Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) carry donor-specific genetic information and have been proposed for exploring the interindividual difference in oncology drug-induced cardiotoxicity. Herein, we evaluated the inter- and intraindividual variability of iPSC-CM-related assays and presented a proof of concept to prospectively predict doxorubicin (DOX)-induced cardiotoxicity (DIC) using donor-specific iPSC-CMs. Our findings demonstrated that donor-specific iPSC-CMs exhibited greater line-to-line variability than the intraindividual variability in impedance cytotoxicity and transcriptome assays. The variable and dose-dependent cytotoxic responses of iPSC-CMs resembled those observed in clinical practice and largely replicated the reported mechanisms. By categorizing iPSC-CMs into resistant and sensitive cell lines based on their time- and concentration-related phenotypic responses to DOX, we found that the sensitivity of donor-specific iPSC-CMs to DOX may predict in vivo DIC risk. Furthermore, we identified a differentially expressed gene, DND microRNA-mediated repression inhibitor 1 (DND1), between the DOX-resistant and DOX-sensitive iPSC-CMs. Our results support the utilization of donor-specific iPSC-CMs in assessing interindividual differences in DIC. Further studies will encompass a large panel of donor-specific iPSC-CMs to identify potential novel molecular and genetic biomarkers for predicting DOX and other oncology drug-induced cardiotoxicity.

Rapid Transient Transcriptional Adaptation to Hypergravity in Jurkat T Cells Revealed by Comparative Analysis of Microarray and RNA-Seq Data

Cellular responses to micro- and hypergravity are rapid and complex and appear within the first few seconds of exposure. Transcriptomic analyses are a valuable tool to analyze these genome-wide cellular alterations. For a better understanding of the cellular dynamics upon altered gravity exposure, it is important to compare different time points. However, since most of the experiments are designed as endpoint measurements, the combination of cross-experiment meta-studies is inevitable. Microarray and RNA-Seq analyses are two of the main methods to study transcriptomics. In the field of altered gravity research, both methods are frequently used. However, the generation of these data sets is difficult and time-consuming and therefore the number of available data sets in this research field is limited. In this study, we investigated the comparability of microarray and RNA-Seq data and applied the results to a comparison of the transcriptomics dynamics between the hypergravity conditions during two real flight platforms and a centrifuge experiment to identify temporal adaptation processes. We performed a comparative study on an Affymetrix HTA2.0 microarray and a paired-end RNA-Seq data set originating from the same Jurkat T cell RNA samples from a short-term hypergravity experiment. The overall agreeability was high, with better sensitivity of the RNA-Seq analysis. The microarray data set showed weaknesses on the level of single upregulated genes, likely due to its normalization approach. On an aggregated level of biotypes, chromosomal distribution, and gene sets, both technologies performed equally well. The microarray showed better performance on the detection of altered gravity-related splicing events. We found that all initially altered transcripts fully adapted after 15 min to hypergravity and concluded that the altered gene expression response to hypergravity is transient and fully reversible. Based on the combined multiple-platform meta-analysis, we could demonstrate rapid transcriptional adaptation to hypergravity, the differential expression of the ATPase subunits ATP6V1A and ATP6V1D, and the cluster of differentiation (CD) molecules CD1E, CD2AP, CD46, CD47, CD53, CD69, CD96, CD164, and CD226 in hypergravity. We could experimentally demonstrate that it is possible to develop methodological evidence for the meta-analysis of individual data.

Identification and validation of endogenous control miRNAs in plasma samples for normalization of qPCR data for Alzheimer's disease

Background

MicroRNAs (miRNAs) are noncoding RNAs that are highly relevant as disease biomarkers. Several studies that explored the role of miRNAs in Alzheimer’s disease (AD) demonstrated their usefulness in clinical identification. Nevertheless, miRNAs that may act as endogenous controls (ECs) have not yet been established. The identification of ECs would contribute to the standardization of these biomarkers in AD. The objective of the study was to identify miRNAs that can be used as ECs in AD.

Methods

We evaluated 145 patients divided into two different cohorts. One was a discovery cohort of 19 women diagnosed with mild to moderate AD (Mini-Mental State Examination (MMSE) score ≥ 20) and with confirmed pathologic levels of Aβ42 in CSF. The stability assessment cohort consisted of 126 individuals: 24 subjects without AD or any kind of dementia and negative for all core CSF biomarkers of AD, 25 subjects with MCI and negative for CSF biomarkers (MCI −), 22 subjects with MCI and positive for CSF biomarkers (MCI +), and 55 subjects with AD and positive for CSF biomarkers. In the discovery cohort, a profile of 384 miRNAs was determined in the plasma by TaqMan low-density array. The best EC candidates were identified by mean-centering and concordance correlation restricted normalization methods. The stability of the EC candidates was assessed using the GeNorm, BestKeeper, and NormFinder algorithms.

Results

Nine miRNAs (hsa-miR-324-5p, hsa-miR-22-5p, hsa-miR-103a-2-5p, hsa-miR-362-5p, hsa-miR-425-3p, hsa-miR-423-5p, hsa-let-7i-3p, hsa-miR-532-5p, and hsa-miR-1301-3p) were identified as EC candidates in the discovery cohort. The validation results indicated that hsa-miR-103a-2-5p was the best EC, followed by hsa-miR-22-5p, hsa-miR-1301-3p, and hsa-miR-425-3p, which had similar stability values in all three algorithms.

Conclusions

We identified a profile of four miRNAs as potential plasma ECs to be used for normalization of miRNA expression data in studies of subjects with cognitive impairment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13195-020-00735-x.

Enabling Routine MHC-II-Associated Peptide Proteomics for Risk Assessment of Drug-Induced Immunogenicity

Major histocompatibility complex-II (MHC-II)-Associated Peptide Proteomics (MAPPs) is a mass spectrometry-based approach to identify and relatively quantitate naturally processed and presented MHC-II-associated peptides that can potentially activate T cells and contribute to the immunogenicity of a drug. Acceptance of the MAPPs technology as an appropriate preclinical (and potentially clinical) immunogenicity risk assessment tool depends not only on its technical stability and robustness but also on the ability to compare results across experiments and donors. To this end, we developed a specialized MAPPs data processing pipeline, dataMAPPs, which presents complex mass spectrometric data sets in the form of heat maps (heatMAPPs), enabling rapid and convenient comparison between conditions and donors. A customized normalization procedure based on identified endogenous peptides standardizes signal intensities within and between donors and enables cross-experimental comparison. We evaluated the technical reproducibility of the MAPPs platform using tool compounds with respect to the most prominent experimental factors and found that the systematic biological differences across donors by far outweighed any technical source of variation. We illustrate the capability of the MAPPs platform to generate data that may be used for preclinical risk assessment of drug-induced immunogenicity and discuss its applicability in the clinics.

Normalization Methods for the Analysis of Unbalanced Transcriptome Data: A Review

Dozens of normalization methods for correcting experimental variation and bias in high-throughput expression data have been developed during the last two decades. Up to 23 methods among them consider the skewness of expression data between sample states, which are even more than the conventional methods, such as loess and quantile. From the perspective of reference selection, we classified the normalization methods for skewed expression data into three categories, data-driven reference, foreign reference, and entire gene set. We separately introduced and summarized these normalization methods designed for gene expression data with global shift between compared conditions, including both microarray and RNA-seq, based on the reference selection strategies. To our best knowledge, this is the most comprehensive review of available preprocessing algorithms for the unbalanced transcriptome data. The anatomy and summarization of these methods shed light on the understanding and appropriate application of preprocessing methods.

CellExpress: a comprehensive microarray-based cancer cell line and clinical sample gene expression analysis online system

Database URL

http://cellexpress.cgm.ntu.edu.tw/.

The 10,000 Immunomes Project: Building a Resource for Human Immunology

There is increasing appreciation that the immune system plays critical roles not only in the traditional domains of infection and inflammation but also in many areas of biology, including tumorigenesis, metabolism, and even neurobiology. However, one of the major barriers for understanding human immunological mechanisms is that immune assays have not been reproducibly characterized for a sufficiently large and diverse healthy human cohort. Here, we present the 10,000 Immunomes Project (10KIP), a framework for growing a diverse human immunology reference, from ImmPort, a publicly available resource of subject-level immunology data. Although some measurement types are sparse in the presently deposited ImmPort database, the extant data allow for a diversity of robust comparisons. Using 10KIP, we describe variations in serum cytokines and leukocytes by age, race, and sex; define a baseline cell-cytokine network; and describe immunologic changes in pregnancy. All data in the resource are available for visualization and download at http://10kimmunomes.org/.

An Optimized Framework for Cancer Prediction Using Immunosignature

Background:

Cancer is a complex disease which can engages the immune system of the patient. In this regard, determination of distinct immunosignatures for various cancers has received increasing interest recently. However, prediction accuracy and reproducibility of the computational methods are limited. In this article, we introduce a robust method for predicting eight types of cancers including astrocytoma, breast cancer, multiple myeloma, lung cancer, oligodendroglia, ovarian cancer, advanced pancreatic cancer, and Ewing sarcoma.

Methods:

In the proposed scheme, at first, the database is normalized with a dictionary of normalization methods that are combined with particle swarm optimization (PSO) for selecting the best normalization method for each feature. Then, statistical feature selection methods are used to separate discriminative features and they were further improved by PSO with appropriate weights as the inputs of the classification system. Finally, the support vector machines, decision tree, and multilayer perceptron neural network were used as classifiers.

Results:

The performance of the hybrid predictor was assessed using the holdout method. According to this method, the minimum sensitivity, specificity, precision, and accuracy of the proposed algorithm were 92.4 ± 1.1, 99.1 ± 1.1, 90.6 ± 2.1, and 98.3 ± 1.0, respectively, among the three types of classification that are used in our algorithm.

Conclusion:

The proposed algorithm considers all the circumstances and works with each feature in its special way. Thus, the proposed algorithm can be used as a promising framework for cancer prediction with immunosignature.

FC1000: normalized gene expression changes of systematically perturbed human cells

The systematic study of transcriptional responses to genetic and chemical perturbations in human cells is still in its early stages. The largest available dataset to date is the newly released L1000 compendium. With its 1.3 million gene expression profiles of treated human cells it offers many opportunities for biomedical data mining, but also data normalization challenges of new dimensions. We developed a novel and practical approach to obtain accurate estimates of fold change response profiles from L1000, based on the RUV (Remove Unwanted Variation) statistical framework. Extending RUV to a big data setting, we propose an estimation procedure, in which an underlying RUV model is tuned by feedback through dataset specific statistical measures, reflecting p-value distributions and internal gene knockdown controls. Applying these metrics - termed evaluation endpoints - to disjoint data splits and integrating the results to select an optimal normalization, the procedure reduces bias and noise in the L1000 data, which in turn broadens the potential of this resource for pharmacological and functional genomic analyses. Our pipeline and normalization results are distributed as an R package (nelanderlab.org/FC1000.html).

An improved analysis methodology for translational profiling by microarray

Translational regulation plays a central role in the global gene expression of a cell, and detection of such regulation has allowed deciphering of critical biological mechanisms. Genome-wide studies of the regulation of translation (translatome) performed on microarrays represent a substantial proportion of studies, alongside with recent advances in deep-sequencing methods. However, there has been a lack of development in specific processing methodologies that deal with the distinct nature of translatome array data. In this study, we confirm that polysome profiling yields skewed data and thus violates the conventional transcriptome analysis assumptions. Using a comprehensive simulation of translatome array data varying the percentage and symmetry of deregulation, we show that conventional analysis methods (Quantile and LOESS normalizations) and statistical tests failed, respectively, to correctly normalize the data and to identify correctly deregulated genes (DEGs). We thus propose a novel analysis methodology available as a CRAN package; Internal Control Analysis of Translatome (INCATome) based on a normalization tied to a group of invariant controls. We confirm that INCATome outperforms the other normalization methods and allows a stringent identification of DEGs. More importantly, INCATome implementation on a biological translatome data set (cells silenced for splicing factor PSF) resulted in the best normalization performance and an improved validation concordance for identification of true positive DEGs. Finally, we provide evidence that INCATome is able to infer novel biological pathways with superior discovery potential, thus confirming the benefits for researchers of implementing INCATome for future translatome studies as well as for existing data sets to generate novel avenues for research.

Quantitative analysis of ChIP-seq data uncovers dynamic and sustained H3K4me3 and H3K27me3 modulation in cancer cells under hypoxia

Background

A comprehensive assessment of the epigenetic dynamics in cancer cells is the key to understanding the molecular mechanisms underlying cancer and to improving cancer diagnostics, prognostics and treatment. By combining genome-wide ChIP-seq epigenomics and microarray transcriptomics, we studied the effects of oxygen deprivation and subsequent reoxygenation on histone 3 trimethylation of lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in a breast cancer cell line, serving as a model for abnormal oxygenation in solid tumors. A priori, epigenetic markings and gene expression levels not only are expected to vary greatly between hypoxic and normoxic conditions, but also display a large degree of heterogeneity across the cell population. Where traditionally ChIP-seq data are often treated as dichotomous data, the model and experiment here necessitate a quantitative, data-driven analysis of both datasets.

Results

We first identified genomic regions with sustained epigenetic markings, which provided a sample-specific reference enabling quantitative ChIP-seq data analysis. Sustained H3K27me3 marking was located around centromeres and intergenic regions, while sustained H3K4me3 marking is associated with genes involved in RNA binding, translation and protein transport and localization. Dynamic marking with both H3K4me3 and H3K27me3 (hypoxia-induced bivalency) was found in CpG-rich regions at loci encoding factors that control developmental processes, congruent with observations in embryonic stem cells.

Conclusions

In silico-identified epigenetically sustained and dynamic genomic regions were confirmed through ChIP-PCR in vitro, and obtained results are corroborated by published data and current insights regarding epigenetic regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0090-4) contains supplementary material, which is available to authorized users.

Hypoxia increases genome-wide bivalent epigenetic marking by specific gain of H3K27me3

Background

Trimethylation at histone H3 lysine 4 (H3K4me3) and lysine 27 (H3K27me3) controls gene activity during development and differentiation. Whether H3K4me3 and H3K27me3 changes dynamically in response to altered microenvironmental conditions, including low-oxygen conditions commonly present in solid tumors, is relatively unknown. Demethylation of H3K4me3 and H3K27me3 is mediated by oxygen and 2-oxoglutarate dioxygenases enzymes, suggesting that oxygen deprivation (hypoxia) may influence histone trimethylation. Using the MCF7 breast epithelial adenocarcinoma cell model, we have determined the relationship between epigenomic and transcriptomic reprogramming as a function of fluctuating oxygen tension.

Results

We find that in MCF7, H3K4me3 and H3K27me3 marks rapidly increase at specific locations throughout the genome and are largely reversed upon reoxygenation. Whereas dynamic changes are relatively highest for H3K27me3 marking under hypoxic conditions, H3K4me3 occupation is identified as the defining epigenetic marker of transcriptional control. In agreement with the global increase of H3K27 trimethylation, we provide direct evidence that the histone H3K27me3 demethylase KDM6B/JMJD3 is inactivated by limited oxygen. In situ immunohistochemical analysis confirms a marked rise of histone trimethylation in hypoxic tumor areas. Acquisition of H3K27me3 at H3K4me3-marked loci results in a striking increase in “bivalent” epigenetic marking. Hypoxia-induced bivalency substantially overlaps with embryonal stem cell-associated genic bivalency and is retained at numerous loci upon reoxygenation. Transcriptional activity is selectively and progressively dampened at bivalently marked loci upon repeated exposure to hypoxia, indicating that this subset of genes uniquely maintains the potential for epigenetic regulation by KDM activity.

Conclusions

These data suggest that dynamic regulation of the epigenetic state within the tumor environment may have important consequences for tumor plasticity and biology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0086-0) contains supplementary material, which is available to authorized users.

A capillary electrophoresis coupled to mass spectrometry pipeline for long term comparable assessment of the urinary metabolome

Although capillary electrophoresis coupled to mass spectrometry (CE-MS) has potential application in the field of metabolite profiling, very few studies actually used CE-MS to identify clinically useful body fluid metabolites. Here we present an optimized CE-MS setup and analysis pipeline to reproducibly explore the metabolite content of urine. We show that the use of a beveled tip capillary improves the sensitivity of detection over a flat tip. We also present a novel normalization procedure based on the use of endogenous stable urinary metabolites identified in the combined metabolome of 75 different urine samples from healthy and diseased individuals. This method allows a highly reproducible comparison of the same sample analyzed nearly 130 times over a range of 4 years. To demonstrate the use of this pipeline in clinical research we compared the urinary metabolome of 34 newborns with ureteropelvic junction (UPJ) obstruction and 15 healthy newborns. We identified 32 features with differential urinary abundance. Combination of the 32 compounds in a SVM classifier predicted with 76% sensitivity and 86% specificity UPJ obstruction in a separate validation cohort of 24 individuals. Thus, this study demonstrates the feasibility to use CE-MS as a tool for the identification of clinically relevant urinary metabolites.

Filtrating colorectal cancer associated genes by integrated analyses of global DNA methylation and hydroxymethylation in cancer and normal tissue

Recently, 5-hydroxymethylcytosine patterning across the tumor genome was considered as a hallmark of cancer development and progression. However, locus-specific difference of hydroxymethylation between colorectal cancer and normal tissue is unknown. In this study, we performed a newly developed method, HMST-seq, to profile 726 aberrant methylated loci and 689 aberrant hydroxymethylated loci synchronously in genome wide of colorectal cancers, majority of which presented higher methylation or lower hydroxymethylationin than in normal group. Besides, abnormal hydroxymethylated modification was more frequently occur at proximal regions close to TSSs and TSSs regions than abnormal methylation. Subsequently, we screened four genes (ALOX15, GHRHR, TFPI2 and TKTL1) with aberrant methylation and aberrant hydroxymethylation at some genome position by functional enrichment analysis as candidate genes associated with colorectal cancer. Our results may allow us to select differentially epigenetically modified target genes implicated in colorectal cancer tumorigenesis.

Practical guidelines and consideration of using RRHP for 5hmC detection

5-hydroxymethylcytosine (5hmC) is an epigenetic modification, which has been associated with gene expression in many biological contexts. Reduced representation hydroxymethylation profiling was developed as an enzymatic-based method for genome-wide 5hmC detection. It exploits β-glucosyltransferase to inhibit enzymatic cleavage of adapters ligated to a genomic library, allowing only fragments with glucosylated 5hmC residues at adapter junctions to be amplified and sequenced. The simple workflow and high sensitivity make it an efficient assay for 5hmC mapping. In this review, we discuss some technical consideration in applying reduced representation hydroxymethylation profiling, such as the use of alternative restriction enzymes for increased genomic coverage in different species, application of control libraries and specifications for multiplexing, data processing and normalization.

Functional analysis and transcriptional output of the Göttingen minipig genome

BACKGROUND

In the past decade the Göttingen minipig has gained increasing recognition as animal model in pharmaceutical and safety research because it recapitulates many aspects of human physiology and metabolism. Genome-based comparison of drug targets together with quantitative tissue expression analysis allows rational prediction of pharmacology and cross-reactivity of human drugs in animal models thereby improving drug attrition which is an important challenge in the process of drug development.

RESULTS

Here we present a new chromosome level based version of the Göttingen minipig genome together with a comparative transcriptional analysis of tissues with pharmaceutical relevance as basis for translational research. We relied on mapping and assembly of WGS (whole-genome-shotgun sequencing) derived reads to the reference genome of the Duroc pig and predict 19,228 human orthologous protein-coding genes. Genome-based prediction of the sequence of human drug targets enables the prediction of drug cross-reactivity based on conservation of binding sites. We further support the finding that the genome of Sus scrofa contains about ten-times less pseudogenized genes compared to other vertebrates. Among the functional human orthologs of these minipig pseudogenes we found HEPN1, a putative tumor suppressor gene. The genomes of Sus scrofa, the Tibetan boar, the African Bushpig, and the Warthog show sequence conservation of all inactivating HEPN1 mutations suggesting disruption before the evolutionary split of these pig species. We identify 133 Sus scrofa specific, conserved long non-coding RNAs (lncRNAs) in the minipig genome and show that these transcripts are highly conserved in the African pigs and the Tibetan boar suggesting functional significance. Using a new minipig specific microarray we show high conservation of gene expression signatures in 13 tissues with biomedical relevance between humans and adult minipigs. We underline this relationship for minipig and human liver where we could demonstrate similar expression levels for most phase I drug-metabolizing enzymes. Higher expression levels and metabolic activities were found for FMO1, AKR/CRs and for phase II drug metabolizing enzymes in minipig as compared to human. The variability of gene expression in equivalent human and minipig tissues is considerably higher in minipig organs, which is important for study design in case a human target belongs to this variable category in the minipig. The first analysis of gene expression in multiple tissues during development from young to adult shows that the majority of transcriptional programs are concluded four weeks after birth. This finding is in line with the advanced state of human postnatal organ development at comparative age categories and further supports the minipig as model for pediatric drug safety studies.

CONCLUSIONS

Genome based assessment of sequence conservation combined with gene expression data in several tissues improves the translational value of the minipig for human drug development. The genome and gene expression data presented here are important resources for researchers using the minipig as model for biomedical research or commercial breeding. Potential impact of our data for comparative genomics, translational research, and experimental medicine are discussed.

Analysis of Reverse Phase Protein Array Data: From Experimental Design towards Targeted Biomarker Discovery

Mastering the systematic analysis of tumor tissues on a large scale has long been a technical challenge for proteomics. In 2001, reverse phase protein arrays (RPPA) were added to the repertoire of existing immunoassays, which, for the first time, allowed a profiling of minute amounts of tumor lysates even after microdissection. A characteristic feature of RPPA is its outstanding sample capacity permitting the analysis of thousands of samples in parallel as a routine task. Until today, the RPPA approach has matured to a robust and highly sensitive high-throughput platform, which is ideally suited for biomarker discovery. Concomitant with technical advancements, new bioinformatic tools were developed for data normalization and data analysis as outlined in detail in this review. Furthermore, biomarker signatures obtained by different RPPA screens were compared with another or with that obtained by other proteomic formats, if possible. Options for overcoming the downside of RPPA, which is the need to steadily validate new antibody batches, will be discussed. Finally, a debate on using RPPA to advance personalized medicine will conclude this article.

GATA2 regulates Wnt signaling to promote primitive red blood cell fate

Primitive erythropoiesis is regulated in a non cell-autonomous fashion across evolution from frogs to mammals. In Xenopus laevis, signals from the overlying ectoderm are required to induce the mesoderm to adopt an erythroid fate. Previous studies in our lab identified the transcription factor GATA2 as a key regulator of this ectodermal signal. To identify GATA2 target genes in the ectoderm required for red blood cell formation in the mesoderm, we used microarray analysis to compare gene expression in ectoderm from GATA2 depleted and wild type embryos. Our analysis identified components of the non-canonical and canonical Wnt pathways as being reciprocally up- and down-regulated downstream of GATA2 in both mesoderm and ectoderm. We show that up-regulation of canonical Wnt signaling during gastrulation blocks commitment to a hematopoietic fate while down-regulation of non-canonical Wnt signaling impairs erythroid differentiation. Our results are consistent with a model in which GATA2 contributes to inhibition of canonical Wnt signaling, thereby permitting progenitors to exit the cell cycle and commit to a hematopoietic fate. Subsequently, activation of non-canonical Wnt signaling plays a later role in enabling these progenitors to differentiate as mature red blood cells.

Data Normalization of (1)H NMR Metabolite Fingerprinting Data Sets in the Presence of Unbalanced Metabolite Regulation

Data normalization is an essential step in NMR-based metabolomics. Conducted properly, it improves data quality and removes unwanted biases. The choice of the appropriate normalization method is critical and depends on the inherent properties of the data set in question. In particular, the presence of unbalanced metabolic regulation, where the different specimens and cohorts under investigation do not contain approximately equal shares of up- and down-regulated features, may strongly influence data normalization. Here, we demonstrate the suitability of the Shapiro-Wilk test to detect such unbalanced regulation. Next, employing a Latin-square design consisting of eight metabolites spiked into a urine specimen at eight different known concentrations, we show that commonly used normalization and scaling methods fail to retrieve true metabolite concentrations in the presence of increasing amounts of glucose added to simulate unbalanced regulation. However, by learning the normalization parameters on a subset of nonregulated features only, Linear Baseline Normalization, Probabilistic Quotient Normalization, and Variance Stabilization Normalization were found to account well for different dilutions of the samples without distorting the true spike-in levels even in the presence of marked unbalanced metabolic regulation. Finally, the methods described were applied successfully to a real world example of unbalanced regulation, namely, a set of plasma specimens collected from patients with and without acute kidney injury after cardiac surgery with cardiopulmonary bypass use.

Transcriptome dynamics of the microRNA inhibition response

We report a high-resolution time series study of transcriptome dynamics following antimiR-mediated inhibition of miR-9 in a Hodgkin lymphoma cell-line—the first such dynamic study of the microRNA inhibition response—revealing both general and specific aspects of the physiological response. We show miR-9 inhibition inducing a multiphasic transcriptome response, with a direct target perturbation before 4 h, earlier than previously reported, amplified by a downstream peak at ∼32 h consistent with an indirect response due to secondary coherent regulation. Predictive modelling indicates a major role for miR-9 in post-transcriptional control of RNA processing and RNA binding protein regulation. Cluster analysis identifies multiple co-regulated gene regulatory modules. Functionally, we observe a shift over time from mRNA processing at early time points to translation at later time points. We validate the key observations with independent time series qPCR and we experimentally validate key predicted miR-9 targets. Methodologically, we developed sensitive functional data analytic predictive methods to analyse the weak response inherent in microRNA inhibition experiments. The methods of this study will be applicable to similar high-resolution time series transcriptome analyses and provides the context for more accurate experimental design and interpretation of future microRNA inhibition studies.

Identification of epigenetic modifications that contribute to pathogenesis in therapy-related AML: Effective integration of genome-wide histone modification with transcriptional profiles

Background

Therapy-related, secondary acute myeloid leukemia (t-AML) is an increasingly frequent complication of intensive chemotherapy. This malignancy is often characterized by abnormalities of chromosome 7, including large deletions or chromosomal loss. A variety of studies suggest that decreased expression of the EZH2 gene located at 7q36.1 is critical in disease pathogenesis. This histone methyltransferase has been implicated in transcriptional repression through modifying histone H3 on lysine 27 (H3k27). However, the critical target genes of EZH2 and their regulatory roles remain unclear.

Method

To characterize the subset of EZH2 target genes that might contribute to t-AML pathogenesis, we developed a novel computational analysis to integrate tissue-specific histone modifications and genome-wide transcriptional regulation. Initial integrative analysis utilized a novel "seq2gene" strategy to explore largely the target genes of chromatin immuneprecipitation sequencing (ChIP-seq) enriched regions. By combining seq2gene with our Phenotype-Genotype-Network (PGNet) algorithm, we enriched genes with similar expression profiles and genomic or functional characteristics into "biomodules".

Results

Initial studies identified SEMA3A (semaphoring 3A) as a novel oncogenic candidate that is regulated by EZH2-silencing, using data derived from both normal and leukemic cell lines as well as murine cells deficient in EZH2. A microsatellite marker at the SEMA3A promoter has been associated with chemosensitivity and radiosensitivity. Notably, our subsequent studies in primary t-AML demonstrate an expected up-regulation of SEMA3A that is EZH2-modulated. Furthermore, we have identified three biomodules that are co-expressed with SEMA3A and up-regulated in t-AML, one of which consists of previously characterized EZH2-repressed gene targets. The other two biomodules include MAPK8 and TATA box targets. Together, our studies suggest an important role for EZH2 targets in t-AML pathogenesis that warrants further study.

Conclusion

These developed computational algorithms and systems biology strategies will enhance the knowledge discovery and hypothesis-driven analysis of multiple next generation sequencing data, for t-AML and other complex diseases.

Tissue elasticity regulated tumor gene expression: implication for diagnostic biomarkers of primitive neuroectodermal tumor

Background

The tumor microenvironment consists of both physical and chemical factors. Tissue elasticity is one physical factor contributing to the microenvironment of tumor cells. To test the importance of tissue elasticity in cell culture, primitive neuroectodermal tumor (PNET) stem cells were cultured on soft polyacrylamide (PAA) hydrogel plates that mimics the elasticity of brain tissue compared with PNET on standard polystyrene (PS) plates. We report the molecular profiles of PNET grown on either PAA or PS.

Methodology/Principal Findings

A whole-genome microarray profile of transcriptional expression between the two culture conditions was performed as a way to probe effects of substrate on cell behavior in culture. The results showed more genes downregulated on PAA compared to PS. This led us to propose microRNA (miRNA) silencing as a potential mechanism for downregulation. Bioinformatic analysis predicted a greater number of miRNA binding sites from the 3' UTR of downregulated genes and identified as specific miRNA binding sites that were enriched when cells were grown on PAA—this supports the hypothesis that tissue elasticity plays a role in influencing miRNA expression. Thus, Dicer was examined to determine if miRNA processing was affected by tissue elasticity. Dicer genes were downregulated on PAA and had multiple predicted miRNA binding sites in its 3' UTR that matched the miRNA binding sites found enriched on PAA. Many differentially regulated genes were found to be present on PS but downregulated on PAA were mapped onto intron sequences. This suggests expression of alternative polyadenylation sites within intron regions that provide alternative 3' UTRs and alternative miRNA binding sites. This results in tissue specific transcriptional downregulation of mRNA in humans by miRNA. We propose a mechanism, driven by the physical characteristics of the microenvironment by which downregulation of genes occur. We found that tissue elasticity-mediated cytokines (TGFβ2 and TNFα) signaling affect expression of ECM proteins.

Conclusions

Our results suggest that tissue elasticity plays important roles in miRNA expression, which, in turn, regulate tumor growth or tumorigenicity.

Integrated analyses of DNA methylation and hydroxymethylation reveal tumor suppressive roles of ECM1, ATF5, and EOMES in human hepatocellular carcinoma

Background

Differences in 5-hydroxymethylcytosine, 5hmC, distributions may complicate previous observations of abnormal cytosine methylation statuses that are used for the identification of new tumor suppressor gene candidates that are relevant to human hepatocarcinogenesis. The simultaneous detection of 5-methylcytosine and 5-hydroxymethylcytosine is likely to stimulate the discovery of aberrantly methylated genes with increased accuracy in human hepatocellular carcinoma.

Results

Here, we performed ultra-performance liquid chromatography/tandem mass spectrometry and single-base high-throughput sequencing, Hydroxymethylation and Methylation Sensitive Tag sequencing, HMST-seq, to synchronously measure these two modifications in human hepatocellular carcinoma samples. After identification of differentially methylated and hydroxymethylated genes in human hepatocellular carcinoma, we integrate DNA copy-number alterations, as determined using array-based comparative genomic hybridization data, with gene expression to identify genes that are potentially silenced by promoter hypermethylation.

Conclusions

We report a high enrichment of genes with epigenetic aberrations in cancer signaling pathways. Six genes were selected as tumor suppressor gene candidates, among which, ECM1, ATF5 and EOMES are confirmed via siRNA experiments to have potential anti-cancer functions.

Electronic supplementary material

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

A comprehensive comparison of normalization methods for loading control and variance stabilization of reverse-phase protein array data

Loading control (LC) and variance stabilization of reverse-phase protein array (RPPA) data have been challenging mainly due to the small number of proteins in an experiment and the lack of reliable inherent control markers. In this study, we compare eight different normalization methods for LC and variance stabilization. The invariant marker set concept was first applied to the normalization of high-throughput gene expression data. A set of “invariant” markers are selected to create a virtual reference sample. Then all the samples are normalized to the virtual reference. We propose a variant of this method in the context of RPPA data normalization and compare it with seven other normalization methods previously reported in the literature. The invariant marker set method performs well with respect to LC, variance stabilization and association with the immunohistochemistry/florescence in situ hybridization data for three key markers in breast tumor samples, while the other methods have inferior performance. The proposed method is a promising approach for improving the quality of RPPA data.

Species and gene divergence in Littorina snails detected by array comparative genomic hybridization

Background

Array comparative genomic hybridization (aCGH) is commonly used to screen different types of genetic variation in humans and model species. Here, we performed aCGH using an oligonucleotide gene-expression array for a non-model species, the intertidal snail Littorina saxatilis. First, we tested what types of genetic variation can be detected by this method using direct re-sequencing and comparison to the Littorina genome draft. Secondly, we performed a genome-wide comparison of four closely related Littorina species: L. fabalis, L. compressa, L. arcana and L. saxatilis and of populations of L. saxatilis found in Spain, Britain and Sweden. Finally, we tested whether we could identify genetic variation underlying “Crab” and “Wave” ecotypes of L. saxatilis.

Results

We could reliably detect copy number variations, deletions and high sequence divergence (i.e. above 3%), but not single nucleotide polymorphisms. The overall hybridization pattern and number of significantly diverged genes were in close agreement with earlier phylogenetic reconstructions based on single genes. The trichotomy of L. arcana, L. compressa and L. saxatilis could not be resolved and we argue that these divergence events have occurred recently and very close in time. We found evidence for high levels of segmental duplication in the Littorina genome (10% of the transcripts represented on the array and up to 23% of the analyzed genomic fragments); duplicated genes and regions were mostly the same in all analyzed species. Finally, this method discriminated geographically distant populations of L. saxatilis, but we did not detect any significant genome divergence associated with ecotypes of L. saxatilis.

Conclusions

The present study provides new information on the sensitivity and the potential use of oligonucleotide arrays for genotyping of non-model organisms. Applying this method to Littorina species yields insights into genome evolution following the recent species radiation and supports earlier single-gene based phylogenies. Genetic differentiation of L. saxatilis ecotypes was not detected in this study, despite pronounced innate phenotypic differences. The reason may be that these differences are due to single-nucleotide polymorphisms.

Electronic supplementary material

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

Improved normalization of systematic biases affecting ion current measurements in label-free proteomics data

Normalization is an important step in the analysis of quantitative proteomics data. If this step is ignored, systematic biases can lead to incorrect assumptions about regulation. Most statistical procedures for normalizing proteomics data have been borrowed from genomics where their development has focused on the removal of so-called 'batch effects.' In general, a typical normalization step in proteomics works under the assumption that most peptides/proteins do not change; scaling is then used to give a median log-ratio of 0. The focus of this work was to identify other factors, derived from knowledge of the variables in proteomics, which might be used to improve normalization. Here we have examined the multi-laboratory data sets from Phase I of the NCI's CPTAC program. Surprisingly, the most important bias variables affecting peptide intensities within labs were retention time and charge state. The magnitude of these observations was exaggerated in samples of unequal concentrations or "spike-in" levels, presumably because the average precursor charge for peptides with higher charge state potentials is lower at higher relative sample concentrations. These effects are consistent with reduced protonation during electrospray and demonstrate that the physical properties of the peptides themselves can serve as good reporters of systematic biases. Between labs, retention time, precursor m/z, and peptide length were most commonly the top-ranked bias variables, over the standardly used average intensity (A). A larger set of variables was then used to develop a stepwise normalization procedure. This statistical model was found to perform as well or better on the CPTAC mock biomarker data than other commonly used methods. Furthermore, the method described here does not require a priori knowledge of the systematic biases in a given data set. These improvements can be attributed to the inclusion of variables other than average intensity during normalization.

Iron deficiency alters expression of dopamine-related genes in the ventral midbrain in mice

A clear link exists between iron deficiency (ID) and nigrostriatal dopamine malfunction. This link appears to play an important role in at least restless legs syndrome (RLS) if not several other neurological diseases. Yet, the underlying mechanisms remain unclear. The effects of ID on gene expression in the brain have not been studied extensively. Here, to better understand how exactly ID alters dopamine functioning, we investigated the effects of ID on gene expression in the brain, seeking to identify any potential transcription-based mechanisms. We used six strains of recombinant inbred mice (BXD type) known to differ in susceptibility to ID in the brain. Upon weaning, we subjected mice from each strain to either an iron-deficient or iron-adequate diet. After 100 days of dietary treatment, we measured the effects of ID on gene expression in the ventral midbrain, a region containing the substantia nigra. The substantia nigra is the base of the nigrostriatal dopamine pathway and a region particularly affected by iron loss in RLS. We screened for ID-induced changes in expression, including changes in that of both iron-regulating and dopamine-related genes. Results revealed a number of expression changes occurring in ID, with large strain-dependent differences in the genes involved and number of expression changes occurring. In terms of dopamine-related genes, results revealed ID-induced expression changes in three genes with direct ties to nigrostriatal dopamine functioning, two of which have never before been implicated in an iron-dopamine pathway. These were stromal cell-derived factor 1 (Cxcl12, or SDF-1), a ferritin regulator and potent dopamine neuromodulator, and hemoglobin, beta adult chain 1 (Hbb-b1), a gene recently shown to play a functional role in dopaminergic neurons. The extent of up-regulation of these genes varied by strain. This work not only demonstrates a wide genetic variation in the transcriptional response to ID in the brain, but also reveals two novel biochemical pathways by which iron may potentially alter dopamine function.

Circulating miRNA biomarkers for Alzheimer's disease

A minimally invasive diagnostic assay for early detection of Alzheimer's disease (AD) is required to select optimal patient groups in clinical trials, monitor disease progression and response to treatment, and to better plan patient clinical care. Blood is an attractive source for biomarkers due to minimal discomfort to the patient, encouraging greater compliance in clinical trials and frequent testing. MiRNAs belong to the class of non-coding regulatory RNA molecules of ∼22 nt length and are now recognized to regulate ∼60% of all known genes through post-transcriptional gene silencing (RNAi). They have potential as useful biomarkers for clinical use because of their stability and ease of detection in many tissues, especially blood. Circulating profiles of miRNAs have been shown to discriminate different tumor types, indicate staging and progression of the disease and to be useful as prognostic markers. Recently their role in neurodegenerative diseases, both as diagnostic biomarkers as well as explaining basic disease etiology has come into focus. Here we report the discovery and validation of a unique circulating 7-miRNA signature (hsa-let-7d-5p, hsa-let-7g-5p, hsa-miR-15b-5p, hsa-miR-142-3p, hsa-miR-191-5p, hsa-miR-301a-3p and hsa-miR-545-3p) in plasma, which could distinguish AD patients from normal controls (NC) with >95% accuracy (AUC of 0.953). There was a >2 fold difference for all signature miRNAs between the AD and NC samples, with p-values<0.05. Pathway analysis, taking into account enriched target mRNAs for these signature miRNAs was also carried out, suggesting that the disturbance of multiple enzymatic pathways including lipid metabolism could play a role in AD etiology.

Analysis of boutique arrays: a universal method for the selection of the optimal data normalization procedure

DNA microarrays, which are among the most popular genomic tools, are widely applied in biology and medicine. Boutique arrays, which are small, spotted, dedicated microarrays, constitute an inexpensive alternative to whole-genome screening methods. The data extracted from each microarray-based experiment must be transformed and processed prior to further analysis to eliminate any technical bias. The normalization of the data is the most crucial step of microarray data pre-processing and this process must be carefully considered as it has a profound effect on the results of the analysis. Several normalization algorithms have been developed and implemented in data analysis software packages. However, most of these methods were designed for whole-genome analysis. In this study, we tested 13 normalization strategies (ten for double-channel data and three for single-channel data) available on R Bioconductor and compared their effectiveness in the normalization of four boutique array datasets. The results revealed that boutique arrays can be successfully normalized using standard methods, but not every method is suitable for each dataset. We also suggest a universal seven-step workflow that can be applied for the selection of the optimal normalization procedure for any boutique array dataset. The described workflow enables the evaluation of the investigated normalization methods based on the bias and variance values for the control probes, a differential expression analysis and a receiver operating characteristic curve analysis. The analysis of each component results in a separate ranking of the normalization methods. A combination of the ranks obtained from all the normalization procedures facilitates the selection of the most appropriate normalization method for the studied dataset and determines which methods can be used interchangeably.

Integrated detection of both 5-mC and 5-hmC by high-throughput tag sequencing technology highlights methylation reprogramming of bivalent genes during cellular differentiation

5-methylcytosine (5-mC) can be oxidized to 5-hydroxymethylcytosine (5-hmC). Genome-wide profiling of 5-hmC thus far indicates 5-hmC may not only be an intermediate form of DNA demethylation but could also constitute an epigenetic mark per se. Here we describe a cost-effective and selective method to detect both the hydroxymethylation and methylation status of cytosines in a subset of cytosines in the human genome. This method involves the selective glucosylation of 5-hmC residues, short-Sequence tag generation and high-throughput sequencing. We tested this method by screening H9 human embryonic stem cells and their differentiated embroid body cells, and found that differential hydroxymethylation preferentially occurs in bivalent genes during cellular differentiation. Especially, our results support hydroxymethylation can regulate key transcription regulators with bivalent marks through demethylation and affect cellular decision on choosing active or inactive state of these genes upon cellular differentiation. Future application of this technology would enable us to uncover the status of methylation and hydroxymethylation in dynamic biological processes and disease development in multiple biological samples.

Coex-Rank: An approach incorporating co-expression information for combined analysis of microarray data

Microarrays have been widely used to study differential gene expression at the genomic level. They can also provide genome-wide co-expression information. Biologically related datasets from independent studies are publicly available, which requires robust combined approaches for integration and validation. Previously, meta-analysis has been adopted to solve this problem. As an alternative to meta-analysis, for microarray data with high similarity in biological experimental design, a more direct combined approach is possible. Gene-level normalization across datasets is motivated by the different scale and distribution of data due to separate origins. However, there has been limited discussion about this point in the past. Here we describe a combined approach for microarray analysis, including gene-level normalization and Coex-Rank approach. After normalization, a linear modeling process is used to identify lists of differentially expressed genes. The Coex-Rank approach incorporates co-expression information into a rank-aggregation procedure. We applied this computational approach to our data, which illustrated an improvement in statistical power and a complementary advantage of the Coex-Rank approach from a biological perspective. Our combined approach for microarray data analysis (Coex-rank) is based on normalization, which is naturally driven. The Coex-rank process not only takes advantage of merging the power of multiple methods regarding normalization but also assists in the discovery of functional clusters of genes.

Systems genetic analysis of the effects of iron deficiency in mouse brain

Iron regulation in the brain is both necessary and highly complex. Too little or too much iron can compromise neurological function, yet we still do not know all of the regulatory processes. In our research, we seek to identify genes and gene networks underlying individual differences in brain iron regulation. To this end, we fed mice from 20+ inbred strains a diet low in iron from weaning to 4 months of age. At sacrifice, we measured iron content in the ventral midbrain (VMB). The VMB contains the substantia nigra, a region particularly vulnerable to iron imbalance. The results showed high, inter-strain variability in dietary iron reduction, from almost no loss to more than 40 % vs. control. When we performed quantitative trait loci (QTL) analysis, we observed a significant area on chromosome 2. Within this QTL, we selected glial high-affinity glutamate transporter 1 (Glt1) as the leading candidate. Expression of this gene is both correlated with VMB iron and is also cis-modulated by local sequence variants that segregate in the BXD family. VMB expression differences of Glt1 in six strains covary with differential susceptibility to VMB iron loss.

Preservation of ranking order in the expression of human Housekeeping genes

Housekeeping (HK) genes fulfill the basic needs for a cell to survive and function properly. Their ubiquitous expression, originally thought to be constant, can vary from tissue to tissue, but this variation remains largely uncharacterized and it could not be explained by previously identified properties of HK genes such as short gene length and high GC content. By analyzing microarray expression data for human genes, we uncovered a previously unnoted characteristic of HK gene expression, namely that the ranking order of their expression levels tends to be preserved from one tissue to another. Further analysis by tensor product decomposition and pathway stratification identified three main factors of the observed ranking preservation, namely that, compared to those of non-HK (NHK) genes, the expression levels of HK genes show a greater degree of dispersion (less overlap), stableness (a smaller variation in expression between tissues), and correlation of expression. Our results shed light on regulatory mechanisms of HK gene expression that are probably different for different HK genes or pathways, but are consistent and coordinated in different tissues.

Normalization of high dimensional genomics data where the distribution of the altered variables is skewed

Genome-wide analysis of gene expression or protein binding patterns using different array or sequencing based technologies is now routinely performed to compare different populations, such as treatment and reference groups. It is often necessary to normalize the data obtained to remove technical variation introduced in the course of conducting experimental work, but standard normalization techniques are not capable of eliminating technical bias in cases where the distribution of the truly altered variables is skewed, i.e. when a large fraction of the variables are either positively or negatively affected by the treatment. However, several experiments are likely to generate such skewed distributions, including ChIP-chip experiments for the study of chromatin, gene expression experiments for the study of apoptosis, and SNP-studies of copy number variation in normal and tumour tissues. A preliminary study using spike-in array data established that the capacity of an experiment to identify altered variables and generate unbiased estimates of the fold change decreases as the fraction of altered variables and the skewness increases. We propose the following work-flow for analyzing high-dimensional experiments with regions of altered variables: (1) Pre-process raw data using one of the standard normalization techniques. (2) Investigate if the distribution of the altered variables is skewed. (3) If the distribution is not believed to be skewed, no additional normalization is needed. Otherwise, re-normalize the data using a novel HMM-assisted normalization procedure. (4) Perform downstream analysis. Here, ChIP-chip data and simulated data were used to evaluate the performance of the work-flow. It was found that skewed distributions can be detected by using the novel DSE-test (Detection of Skewed Experiments). Furthermore, applying the HMM-assisted normalization to experiments where the distribution of the truly altered variables is skewed results in considerably higher sensitivity and lower bias than can be attained using standard and invariant normalization methods.

Kernel density weighted loess normalization improves the performance of detection within asymmetrical data

Background

Normalization of gene expression data has been studied for many years and various strategies have been formulated to deal with various types of data. Most normalization algorithms rely on the assumption that the number of up-regulated genes and the number of down-regulated genes are roughly the same. However, the well-known Golden Spike experiment presents a unique situation in which differentially regulated genes are biased toward one direction, thereby challenging the conclusions of previous bench mark studies.

Results

This study proposes two novel approaches, KDL and KDQ, based on kernel density estimation to improve upon the basic idea of invariant set selection. The key concept is to provide various importance scores to data points on the MA plot according to their proximity to the cluster of the null genes under the assumption that null genes are more densely distributed than those that are differentially regulated. The comparison is demonstrated in the Golden Spike experiment as well as with simulation data using the ROC curves and compression rates. KDL and KDQ in combination with GCRMA provided the best performance among all approaches.

Conclusions

This study determined that methods based on invariant sets are better able to resolve the problem of asymmetry. Normalization, either before or after expression summary for probesets, improves performance to a similar degree.

Extensive increase of microarray signals in cancers calls for novel normalization assumptions

When using microarray data for studying a complex disease such as cancer, it is a common practice to normalize data to force all arrays to have the same distribution of probe intensities regardless of the biological groups of samples. The assumption underlying such normalization is that in a disease the majority of genes are not differentially expressed genes (DE genes) and the numbers of up- and down-regulated genes are roughly equal. However, accumulated evidences suggest gene expressions could be widely altered in cancer, so we need to evaluate the sensitivities of biological discoveries to violation of the normalization assumption. Here, we analyzed 7 large Affymetrix datasets of pair-matched normal and cancer samples for cancers collected in the NCBI GEO database. We showed that in 6 of these 7 datasets, the medians of perfect match (PM) probe intensities increased in cancer state and the increases were significant in three datasets, suggesting the assumption that all arrays have the same median probe intensities regardless of the biological groups of samples might be misleading. Then, we evaluated the effects of three currently most widely used normalization algorithms (RMA, MAS5.0 and dChip) on the selection of DE genes by comparing them with LVS which relies less on the above-mentioned assumption. The results showed using RMA, MAS5.0 and dChip may produce lots of false results of down-regulated DE genes while missing many up-regulated DE genes. At least for cancer study, normalizing all arrays to have the same distribution of probe intensities regardless of the biological groups of samples might be misleading. Thus, most current normalizations based on unreliable assumptions may distort biological differences between normal and cancer samples. The LVS algorithm might perform relatively well due to that it relies less on the above-mentioned assumption. Also, our results indicate that genes may be widely up-regulated in most human cancer.

Experimental design, preprocessing, normalization and differential expression analysis of small RNA sequencing experiments

Prior to the advent of new, deep sequencing methods, small RNA (sRNA) discovery was dependent on Sanger sequencing, which was time-consuming and limited knowledge to only the most abundant sRNA. The innovation of large-scale, next-generation sequencing has exponentially increased knowledge of the biology, diversity and abundance of sRNA populations. In this review, we discuss issues involved in the design of sRNA sequencing experiments, including choosing a sequencing platform, inherent biases that affect sRNA measurements and replication. We outline the steps involved in preprocessing sRNA sequencing data and review both the principles behind and the current options for normalization. Finally, we discuss differential expression analysis in the absence and presence of biological replicates. While our focus is on sRNA sequencing experiments, many of the principles discussed are applicable to the sequencing of other RNA populations.

GRISSOM platform: enabling distributed processing and management of biological data through fusion of grid and web technologies

Transcriptomic technologies have a critical impact in the revolutionary changes that reshape biological research. Through the recruitment of novel high-throughput instrumentation and advanced computational methodologies, an unprecedented wealth of quantitative data is produced. Microarray experiments are considered high-throughput, both in terms of data volumes (data intensive) and processing complexity (computationally intensive). In this paper, we present grids for in silico systems biology and medicine (GRISSOM), a web-based application that exploits GRID infrastructures for distributed data processing and management, of DNA microarrays (cDNA, Affymetrix, Illumina) through a generic, consistent, computational analysis framework. GRISSOM performs versatile annotation and integrative analysis tasks, through the use of third-party application programming interfaces, delivered as web services. In parallel, by conforming to service-oriented architectures, it can be encapsulated in other biomedical processing workflows, with the help of workflow enacting software, like Taverna Workbench, thus rendering access to its algorithms, transparent and generic. GRISSOM aims to set a generic paradigm of efficient metamining that promotes translational research in biomedicine, through the fusion of grid and semantic web computing technologies.

Cold-induced silencing by long antisense transcripts of an Arabidopsis Polycomb target

Transcription in eukaryotic genomes generates an extensive array of non-protein-coding RNA, the functional significance of which is mostly unknown. We are investigating the link between non-coding RNA and chromatin regulation through analysis of FLC - a regulator of flowering time in Arabidopsis and a target of several chromatin pathways. Here we use an unbiased strategy to characterize non-coding transcripts of FLC and show that sense/antisense transcript levels correlate in a range of mutants and treatments, but change independently in cold-treated plants. Prolonged cold epigenetically silences FLC in a Polycomb-mediated process called vernalization. Our data indicate that upregulation of long non-coding antisense transcripts covering the entire FLC locus may be part of the cold-sensing mechanism. Induction of these antisense transcripts occurs earlier than, and is independent of, other vernalization markers and coincides with a reduction in sense transcription. We show that addition of the FLC antisense promoter sequences to a reporter gene is sufficient to confer cold-induced silencing of the reporter. Our data indicate that cold-induced FLC antisense transcripts have an early role in the epigenetic silencing of FLC, acting to silence FLC transcription transiently. Recruitment of the Polycomb machinery then confers the epigenetic memory. Antisense transcription events originating from 3' ends of genes might be a general mechanism to regulate the corresponding sense transcription in a condition/stage-dependent manner.

Gene expression profiling in the rhesus macaque: methodology, annotation and data interpretation

Gene microarray analyses represent potentially effective means for high-throughput gene expression profiling in non-human primates. In the companion article, we emphasize effective experimental design based on the in vivo physiology of the rhesus macaque, whereas this article emphasizes considerations for gene annotation and data interpretation using gene microarray platforms from Affymetrix. Initial annotation of the rhesus genome array was based on Affymetrix human GeneChips. However, annotation revisions improve the precision with which rhesus transcripts are identified. Annotation of the rhesus GeneChip is under continuous revision with large percentages of probesets under multiple annotation systems having undergone multiple reassignments between March 2007 and November 2008. It is also important to consider that quantitation and comparison of gene expression levels across multiple chips requires appropriate normalization. External corroboration of microarray results using PCR-based methodology also requires validation of appropriate internal reference genes for normalization of expression values. Many tools are now freely available to aid investigators with microarray normalization and selection of internal reference genes to be used for independent corroboration of microarray results.