AffyRNADegradation: control and correction of RNA quality effects in GeneChip expression data

The integration of large-scale public data and network analysis uncovers molecular characteristics of psoriasis

In recent years, a growing interest in the characterization of the molecular basis of psoriasis has been observed. However, despite the availability of a large amount of molecular data, many pathogenic mechanisms of psoriasis are still poorly understood. In this study, we performed an integrated analysis of 23 public transcriptomic datasets encompassing both lesional and uninvolved skin samples from psoriasis patients. We defined comprehensive gene co-expression network models of psoriatic lesions and uninvolved skin. Moreover, we curated and exploited a wide range of functional information from multiple public sources in order to systematically annotate the inferred networks. The integrated analysis of transcriptomics data and co-expression networks highlighted genes that are frequently dysregulated and show aberrant patterns of connectivity in the psoriatic lesion compared with the unaffected skin. Our approach allowed us to also identify plausible, previously unknown, actors in the expression of the psoriasis phenotype. Finally, we characterized communities of co-expressed genes associated with relevant molecular functions and expression signatures of specific immune cell types associated with the psoriasis lesion. Overall, integrating experimental driven results with curated functional information from public repositories represents an efficient approach to empower knowledge generation about psoriasis and may be applicable to other complex diseases.

Genome-scale meta-analysis of breast cancer datasets identifies promising targets for drug development

Background

Because of the highly heterogeneous nature of breast cancer, each subtype differs in response to several treatment regimens. This has limited the therapeutic options for metastatic breast cancer disease requiring exploration of diverse therapeutic models to target tumor specific biomarkers.

Methods

Differentially expressed breast cancer genes identified through extensive data mapping were studied for their interaction with other target proteins involved in breast cancer progression. The molecular mechanisms by which these signature genes are involved in breast cancer metastasis were also studied through pathway analysis. The potential drug targets for these genes were also identified.

Results

From 50 DEGs, 20 genes were identified based on fold change and p-value and the data curation of these genes helped in shortlisting 8 potential gene signatures that can be used as potential candidates for breast cancer. Their network and pathway analysis clarified the role of these genes in breast cancer and their interaction with other signaling pathways involved in the progression of disease metastasis. The miRNA targets identified through miRDB predictor provided potential miRNA targets for these genes that can be involved in breast cancer progression. Several FDA approved drug targets were identified for the signature genes easing the therapeutic options for breast cancer treatment.

Conclusion

The study provides a more clarified role of signature genes, their interaction with other genes as well as signaling pathways. The miRNA prediction and the potential drugs identified will aid in assessing the role of these targets in breast cancer.

Manually curated and harmonised transcriptomics datasets of psoriasis and atopic dermatitis patients

We present manually curated transcriptomics data of psoriasis and atopic dermatitis patients retrieved from the NCBI Gene Expression Omnibus and EBI ArrayExpress repositories. We collected 39 transcriptomics datasets, deriving from DNA microarrays and RNA-Sequencing technologies, for a total of 1677 samples. We provide quality-checked, homogenised and preprocessed gene expression matrices and their corresponding metadata tables along with the estimated surrogate variables. These data represent a ready-made valuable source of knowledge for translational researchers in the dermatology field.

Analysis of global gene expression at seven brain regions of patients with schizophrenia

Previous transcriptome analyses of brain samples provided several insights into the pathophysiology of schizophrenia. In this study, we aimed to re-investigate gene expression datasets from seven brain regions of patients with schizophrenia and healthy controls by adopting a unified approach. After adjustment for confounding factors, we detected gene expression changes in 2 out of 7 brain regions - the dorsolateral prefrontal cortex (DLPFC) and parietal cortex (PC). We found relatively small effect sizes, not exceeding absolute log fold changes of 1. Gene-set enrichment analysis revealed the following alterations: 1) down-regulation of GABAergic signaling (in DLPFC and PC); 2) up-regulation of interleukin-23 signaling together with up-regulation of transcription mediated by RUNX1 and RUNX3 as well as down-regulation of RUNX2 signaling (in DLPFC) and 3) up-regulation of genes associated with responses to metal ions and RUNX1 signaling (PC). The number of neurons was significantly lower and the number of astrocytes was significantly higher at both brain regions. In turn, the index of microglia was increased in DLPFC and decreased in PC. Finally, our unsupervised analysis demonstrated that cellular composition of the samples was a major confounding factor in the analysis of gene expression across all datasets. In conclusion, our analysis provides further evidence that small but significant changes in the expression of genes related to GABAergic signaling, brain development, neuroinflammation and responses to metal ions might be involved in the pathophysiology of schizophrenia. Cell sorting techniques need to be used by future studies to dissect the effect of cellular content.

Transcriptomics in Toxicogenomics, Part II: Preprocessing and Differential Expression Analysis for High Quality Data

Preprocessing of transcriptomics data plays a pivotal role in the development of toxicogenomics-driven tools for chemical toxicity assessment. The generation and exploitation of large volumes of molecular profiles, following an appropriate experimental design, allows the employment of toxicogenomics (TGx) approaches for a thorough characterisation of the mechanism of action (MOA) of different compounds. To date, a plethora of data preprocessing methodologies have been suggested. However, in most cases, building the optimal analytical workflow is not straightforward. A careful selection of the right tools must be carried out, since it will affect the downstream analyses and modelling approaches. Transcriptomics data preprocessing spans across multiple steps such as quality check, filtering, normalization, batch effect detection and correction. Currently, there is a lack of standard guidelines for data preprocessing in the TGx field. Defining the optimal tools and procedures to be employed in the transcriptomics data preprocessing will lead to the generation of homogeneous and unbiased data, allowing the development of more reliable, robust and accurate predictive models. In this review, we outline methods for the preprocessing of three main transcriptomic technologies including microarray, bulk RNA-Sequencing (RNA-Seq), and single cell RNA-Sequencing (scRNA-Seq). Moreover, we discuss the most common methods for the identification of differentially expressed genes and to perform a functional enrichment analysis. This review is the second part of a three-article series on Transcriptomics in Toxicogenomics.

Inference of cell type content from human brain transcriptomic datasets illuminates the effects of age, manner of death, dissection, and psychiatric diagnosis

Psychiatric illness is unlikely to arise from pathology occurring uniformly across all cell types in affected brain regions. Despite this, transcriptomic analyses of the human brain have typically been conducted using macro-dissected tissue due to the difficulty of performing single-cell type analyses with donated post-mortem brains. To address this issue statistically, we compiled a database of several thousand transcripts that were specifically-enriched in one of 10 primary cortical cell types in previous publications. Using this database, we predicted the relative cell type content for 833 human cortical samples using microarray or RNA-Seq data from the Pritzker Consortium (GSE92538) or publicly-available databases (GSE53987, GSE21935, GSE21138, CommonMind Consortium). These predictions were generated by averaging normalized expression levels across transcripts specific to each cell type using our R-package BrainInABlender (validated and publicly-released on github). Using this method, we found that the principal components of variation in the datasets strongly correlated with the predicted neuronal/glial content of the samples. This variability was not simply due to dissection–the relative balance of brain cell types appeared to be influenced by a variety of demographic, pre- and post-mortem variables. Prolonged hypoxia around the time of death predicted increased astrocytic and endothelial gene expression, illustrating vascular upregulation. Aging was associated with decreased neuronal gene expression. Red blood cell gene expression was reduced in individuals who died following systemic blood loss. Subjects with Major Depressive Disorder had decreased astrocytic gene expression, mirroring previous morphometric observations. Subjects with Schizophrenia had reduced red blood cell gene expression, resembling the hypofrontality detected in fMRI experiments. Finally, in datasets containing samples with especially variable cell content, we found that controlling for predicted sample cell content while evaluating differential expression improved the detection of previously-identified psychiatric effects. We conclude that accounting for cell type can greatly improve the interpretability of transcriptomic data.

Transcriptomic Studies of Malaria: a Paradigm for Investigation of Systemic Host-Pathogen Interactions

Transcriptomics, the analysis of genome-wide RNA expression, is a common approach to investigate host and pathogen processes in infectious diseases. Technical and bioinformatic advances have permitted increasingly thorough analyses of the association of RNA expression with fundamental biology, immunity, pathogenesis, diagnosis, and prognosis. Transcriptomic approaches can now be used to realize a previously unattainable goal, the simultaneous study of RNA expression in host and pathogen, in order to better understand their interactions. This exciting prospect is not without challenges, especially as focus moves from interactions in vitro under tightly controlled conditions to tissue- and systems-level interactions in animal models and natural and experimental infections in humans. Here we review the contribution of transcriptomic studies to the understanding of malaria, a parasitic disease which has exerted a major influence on human evolution and continues to cause a huge global burden of disease. We consider malaria a paradigm for the transcriptomic assessment of systemic host-pathogen interactions in humans, because much of the direct host-pathogen interaction occurs within the blood, a readily sampled compartment of the body. We illustrate lessons learned from transcriptomic studies of malaria and how these lessons may guide studies of host-pathogen interactions in other infectious diseases. We propose that the potential of transcriptomic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in study design rather than as a consequence of technological constraints. Further advances will require the integration of transcriptomic data with analytical approaches from other scientific disciplines, including epidemiology and mathematical modeling.

Immunological landscape of consensus clusters in colorectal cancer

Recent, large-scale expression–based subtyping has advanced our understanding of the genomic landscape of colorectal cancer (CRC) and resulted in a consensus molecular classification that enables the categorization of most CRC tumors into one of four consensus molecular subtypes (CMS). Currently, major progress in characterization of immune landscape of tumor-associated microenvironment has been made especially with respect to microsatellite status of CRCs. While these studies profoundly improved the understanding of molecular and immunological profile of CRCs heterogeneity less is known about repertoire of the tumor infiltrating immune cells of each CMS.

In order to comprehensively characterize the immune landscape of CRC we re-analyzed a total of 15 CRC genome-wide expression data sets encompassing 1597 tumors and 125 normal adjacent colon tissues. After quality filtering, CRC clusters were discovered using a combination of multiple clustering algorithms and multiple validity metrics. CIBERSORT algorithm was used to compute relative proportions of 22 human leukocyte subpopulations across CRC clusters and normal colon tissue. Subsequently, differential expression specific to tumor epithelial cells was calculated to characterize mechanisms of tumor escape from immune surveillance occurring in particular CRC clusters.

Our results not only characterize the common and cluster-specific influx of immune cells into CRCs but also identify several deregulated gene targets that may contribute to improvement of immunotherapeutic strategies in CRC.

Reduced number of peripheral natural killer cells in schizophrenia but not in bipolar disorder

Overwhelming evidence indicates that subthreshold inflammatory state might be implicated in the pathophysiology of schizophrenia (SCZ) and bipolar disorder (BPD). It has been reported that both groups of patients might be characterized by abnormal lymphocyte counts. However, little is known about alterations in lymphocyte proportions that may differentiate SCZ and BPD patients. Therefore, in this study we investigated blood cell proportions quantified by means of microarray expression deconvolution using publicly available data from SCZ and BPD patients. We found significantly lower counts of natural killer (NK) cells in drug-naïve and medicated SCZ patients compared to healthy controls across all datasets. In one dataset from SCZ patients, there were no significant differences in the number of NK cells between acutely relapsed and remitted SCZ patients. No significant difference in the number of NK cells between BPD patients and healthy controls was observed in all datasets. Our results indicate that SCZ patients, but not BPD patients, might be characterized by reduced counts of NK cells. Future studies looking at lymphocyte counts in SCZ should combine the analysis of data obtained using computational deconvolution and flow cytometry techniques.

Microarray experiments and factors which affect their reliability

Oligonucleotide microarrays belong to the basic tools of molecular biology and allow for simultaneous assessment of the expression level of thousands of genes. Analysis of microarray data is however very complex, requiring sophisticated methods to control for various factors that are inherent to the procedures used. In this article we describe the individual steps of a microarray experiment, highlighting important elements and factors that may affect the processes involved and that influence the interpretation of the results. Additionally, we describe methods that can be used to estimate the influence of these factors, and to control the way in which they affect the expression estimates. A comprehensive understanding of the experimental protocol used in a microarray experiment aids the interpretation of the obtained results. By describing known factors which affect expression estimates this article provides guidelines for appropriate quality control and pre-processing of the data, additionally applicable to other transcriptome analysis methods that utilize similar sample handling protocols.

Gene expression profiling of sequential metastatic biopsies for biomarker discovery in breast cancer

The feasibility of longitudinal metastatic biopsies for gene expression profiling in breast cancer is unexplored. Dynamic changes in gene expression can potentially predict efficacy of targeted cancer drugs. Patients enrolled in a phase III trial of metastatic breast cancer with docetaxel monotherapy versus combination of docetaxel + sunitinib were offered to participate in a translational substudy comprising longitudinal fine needle aspiration biopsies and Positron Emission Tomography imaging before (T1) and two weeks after start of treatment (T2). Aspirated tumor material was used for microarray analysis, and treatment-induced changes (T2 versus T1) in gene expression and standardized uptake values (SUV) were investigated and correlated to clinical outcome measures. Gene expression profiling yielded high-quality data at both time points in 14/18 patients. Unsupervised clustering revealed specific patterns of changes caused by monotherapy vs. combination therapy (p = 0.021, Fisher's exact test). A therapy-induced reduction of known proliferation and hypoxia metagene scores was prominent in the combination arm. Changes in a previously reported hypoxia metagene score were strongly correlated to the objective responses seen by conventional radiology assessments after 6 weeks in the combination arm, Spearman's ρ = 1 (p = 0.017) but not in monotherapy, ρ = -0.029 (p = 1). Similarly, the Predictor Analysis of Microarrays 50 (PAM50) proliferation metagene correlated to tumor changes merely in the combination arm at 6 and 12 weeks (ρ = 0.900, p = 0.083 and ρ = 1, p = 0.017 respectively). Reductions in mean SUV were a reliable early predictor of objective response in monotherapy, ρ = 0.833 (p = 0.008), but not in the combination arm ρ = -0.029 (p = 1). Gene expression profiling of longitudinal metastatic aspiration biopsies was feasible, demonstrated biological validity and provided predictive information.

Variation of RNA Quality and Quantity Are Major Sources of Batch Effects in Microarray Expression Data

The great utility of microarrays for genome-scale expression analysis is challenged by the widespread presence of batch effects, which bias expression measurements in particular within large data sets. These unwanted technical artifacts can obscure biological variation and thus significantly reduce the reliability of the analysis results. It is largely unknown which are the predominant technical sources leading to batch effects. We here quantitatively assess the prevalence and impact of several known technical effects on microarray expression results. Particularly, we focus on important factors such as RNA degradation, RNA quantity, and sequence biases including multiple guanine effects. We find that the common variation of RNA quality and RNA quantity can not only yield low-quality expression results, but that both factors also correlate with batch effects and biological characteristics of the samples.