Search

High throughput RNA sequencing experiments are widely conducted and analyzed to identify differentially expressed genes (DEGs). The statistical models calculated for this task are often not clear to practitioners, and analyses may not be optimally tailored to the research hypothesis. Often, interaction effects (IEs) are the mathematical equivalent of the biological research question but are not considered for different reasons. We fill this gap by explaining and presenting the potential benefit of IEs in the search for DEGs using RNA-Seq data of mice that receive different diets for different time periods. Using an IE model leads to a smaller, but likely more biologically informative set of DEGs compared to a common approach that avoids the calculation of IEs., (© 2023. The Author(s).)

The in vitro test battery of the European research consortium ESNATS ('novel stem cell-based test systems') has been used to screen for potential human developmental toxicants. As part of this effort, the migration of neural crest (MINC) assay has been used to evaluate chemical effects on neural crest function. It identified some drug-like compounds in addition to known environmental toxicants. The hits included the HSP90 inhibitor geldanamycin, the chemotherapeutic arsenic trioxide, the flame-retardant PBDE-99, the pesticide triadimefon and the histone deacetylase inhibitors valproic acid and trichostatin A. Transcriptome changes triggered by these substances in human neural crest cells were recorded and analysed here to answer three questions: (1) can toxicants be individually identified based on their transcript profile; (2) how can the toxicity pattern reflected by transcript changes be compacted/dimensionality-reduced for practical regulatory use; (3) how can a reduced set of biomarkers be selected for large-scale follow-up? Transcript profiling allowed clear separation of different toxicants and the identification of toxicant types in a blinded test study. We also developed a diagrammatic system to visualize and compare toxicity patterns of a group of chemicals by giving a quantitative overview of altered superordinate biological processes (e.g. activation of KEGG pathways or overrepresentation of gene ontology terms). The transcript data were mined for potential markers of toxicity, and 39 transcripts were selected to either indicate general developmental toxicity or distinguish compounds with different modes-of-action in read-across. In summary, we found inclusion of transcriptome data to largely increase the information from the MINC phenotypic test.

The comprehensive transcriptomic analysis of clinically annotated human tissue has found widespread use in oncology, cell biology, immunology, and toxicology. In cancer research, microarray-based gene expression profiling has successfully been applied to subclassify disease entities, predict therapy response, and identify cellular mechanisms. Public accessibility of raw data, together with corresponding information on clinicopathological parameters, offers the opportunity to reuse previously analyzed data and to gain statistical power by combining multiple datasets. However, results and conclusions obviously depend on the reliability of the available information. Here, we propose gene expression-based methods for identifying sample misannotations in public transcriptomic datasets. Sample mix-up can be detected by a classifier that differentiates between samples from male and female patients. Correlation analysis identifies multiple measurements of material from the same sample. The analysis of 45 datasets (including 4913 patients) revealed that erroneous sample annotation, affecting 40 % of the analyzed datasets, may be a more widespread phenomenon than previously thought. Removal of erroneously labelled samples may influence the results of the statistical evaluation in some datasets. Our methods may help to identify individual datasets that contain numerous discrepancies and could be routinely included into the statistical analysis of clinical gene expression data.

Purpose: Global gene expression profiling has been widely used in lung cancer research to identify clinically relevant molecular subtypes as well as to predict prognosis and therapy response. So far, the value of these multigene signatures in clinical practice is unclear, and the biologic importance of individual genes is difficult to assess, as the published signatures virtually do not overlap., Experimental Design: Here, we describe a novel single institute cohort, including 196 non-small lung cancers (NSCLC) with clinical information and long-term follow-up. Gene expression array data were used as a training set to screen for single genes with prognostic impact. The top 450 probe sets identified using a univariate Cox regression model (significance level P < 0.01) were tested in a meta-analysis including five publicly available independent lung cancer cohorts (n = 860)., Results: The meta-analysis revealed 14 genes that were significantly associated with survival (P < 0.001) with a false discovery rate <1%. The prognostic impact of one of these genes, the cell adhesion molecule 1 (CADM1), was confirmed by use of immunohistochemistry on tissue microarrays from 2 independent NSCLC cohorts, altogether including 617 NSCLC samples. Low CADM1 protein expression was significantly associated with shorter survival, with particular influence in the adenocarcinoma patient subgroup., Conclusions: Using a novel NSCLC cohort together with a meta-analysis validation approach, we have identified a set of single genes with independent prognostic impact. One of these genes, CADM1, was further established as an immunohistochemical marker with a potential application in clinical diagnostics.

Developmental neurotoxicity (DNT) and many forms of reproductive toxicity (RT) often manifest themselves in functional deficits that are not necessarily based on cell death, but rather on minor changes relating to cell differentiation or communication. The fields of DNT/RT would greatly benefit from in vitro tests that allow the identification of toxicant-induced changes of the cellular proteostasis, or of its underlying transcriptome network. Therefore, the 'human embryonic stem cell (hESC)-derived novel alternative test systems (ESNATS)' European commission research project established RT tests based on defined differentiation protocols of hESC and their progeny. Valproic acid (VPA) and methylmercury (MeHg) were used as positive control compounds to address the following fundamental questions: (1) Does transcriptome analysis allow discrimination of the two compounds? (2) How does analysis of enriched transcription factor binding sites (TFBS) and of individual probe sets (PS) distinguish between test systems? (3) Can batch effects be controlled? (4) How many DNA microarrays are needed? (5) Is the highest non-cytotoxic concentration optimal and relevant for the study of transcriptome changes? VPA triggered vast transcriptional changes, whereas MeHg altered fewer transcripts. To attenuate batch effects, analysis has been focused on the 500 PS with highest variability. The test systems differed significantly in their responses (<20 % overlap). Moreover, within one test system, little overlap between the PS changed by the two compounds has been observed. However, using TFBS enrichment, a relatively large 'common response' to VPA and MeHg could be distinguished from 'compound-specific' responses. In conclusion, the ESNATS assay battery allows classification of human DNT/RT toxicants on the basis of their transcriptome profiles.

Purpose: Although the central role of the immune system for tumor prognosis is generally accepted, a single robust marker is not yet available., Experimental Design: On the basis of receiver operating characteristic analyses, robust markers were identified from a 60-gene B cell-derived metagene and analyzed in gene expression profiles of 1,810 breast cancer; 1,056 non-small cell lung carcinoma (NSCLC); 513 colorectal; and 426 ovarian cancer patients. Protein and RNA levels were examined in paraffin-embedded tissue of 330 breast cancer patients. The cell types were identified with immunohistochemical costaining and confocal fluorescence microscopy., Results: We identified immunoglobulin κ C (IGKC) which as a single marker is similarly predictive and prognostic as the entire B-cell metagene. IGKC was consistently associated with metastasis-free survival across different molecular subtypes in node-negative breast cancer (n = 965) and predicted response to anthracycline-based neoadjuvant chemotherapy (n = 845; P < 0.001). In addition, IGKC gene expression was prognostic in NSCLC and colorectal cancer. No association was observed in ovarian cancer. IGKC protein expression was significantly associated with survival in paraffin-embedded tissues of 330 breast cancer patients. Tumor-infiltrating plasma cells were identified as the source of IGKC expression., Conclusion: Our findings provide IGKC as a novel diagnostic marker for risk stratification in human cancer and support concepts to exploit the humoral immune response for anticancer therapy. It could be validated in several independent cohorts and carried out similarly well in RNA from fresh frozen as well as from paraffin tissue and on protein level by immunostaining., (©2012 AACR.)

Objectives: We characterize typical problems encountered in microarray image analysis and present algorithmic approaches dealing with background estimation, spot identification and intensity extraction. Validation of the quality of resulting measurements is discussed., Methods: We describe sources for errors in microarray images and present algorithms that have been specifically developed to deal with such experimental imperfections., Results: For the image analysis of hybridization experiments, discriminating spot regions from a background is the most critical step. Spot shape detection algorithms, intensity histogram methods and hybrid approaches have been proposed. The correctness of final intensity estimates is difficult to verify. Nevertheless, the application of sophisticated algorithms provides a significant reduction of the possible information loss., Conclusions: The initial analysis step for array hybridization experiments is the estimation of expression intensities. The quality of this process is crucial for the validity of interpretations from subsequent analysis steps.

Objectives: We introduce methods for the exploratory analysis of microarray data, especially focusing on cluster algorithms. Benefits and problems are discussed., Methods: We describe application and suitability of unsupervised learning methods for the classification of gene expression data. Cluster algorithms are treated in more detail, including assessment of cluster quality., Results: When dealing with microarray data, most cluster algorithms must be applied with caution. As long as the structure of the true generating models of such data is not fully understood, the use of simple algorithms seems to be more appropriate than the application of complex black-box algorithms. New methods explicitly targeted to the analysis of microarray data are increasingly being developed in order to increase the amount of useful information extracted from the experiments., Conclusions: Unsupervised methods can be a helpful tool for the analysis of microarray data, but a critical choice of the algorithm and a careful interpretation of the results are required in order to avoid false conclusions.

Background: Image analysis is the first crucial step to obtain reliable results from microarray experiments. First, areas in the image belonging to single spots have to be identified. Then, those target areas have to be partitioned into foreground and background. Finally, two scalar values for the intensities have to be extracted. These goals have been tackled either by spot shape methods or intensity histogram methods, but it would be desirable to have hybrid algorithms which combine the advantages of both approaches., Results: A new robust and adaptive histogram type method is pixel clustering, which has been successfully applied for detecting and quantifying microarray spots. This paper demonstrates how the spot shape can be effectively integrated in this approach. Based on the clustering results, a bivalence mask is constructed. It estimates the expected spot shape and is used to filter the data, improving the results of the cluster algorithm. The quality measure 'stability' is defined and evaluated on a real data set. The improved clustering method is compared with the established Spot software on a data set with replicates., Conclusion: The new method presents a successful hybrid microarray image analysis solution. It incorporates both shape and histogram features and is specifically adapted to deal with typical microarray image characteristics. As a consequence of the filtering step pixels are divided into three groups, namely foreground, background and deletions. This allows a separate treatment of artifacts and their elimination from the further analysis.

Motivation: Microarray images challenge existing analytical methods in many ways given that gene spots are often comprised of characteristic imperfections. Irregular contours, donut shapes, artifacts, and low or heterogeneous expression impair corresponding values for red and green intensities as well as their ratio R/G. New approaches are needed to ensure accurate data extraction from these images., Results: Herein we introduce a novel method for intensity assessment of gene spots. The technique is based on clustering pixels of a target area into foreground and background. For this purpose we implemented two clustering algorithms derived from k-means and Partitioning Around Medoids (PAM), respectively. Results from the analysis of real gene spots indicate that our approach performs superior to other existing analytical methods. This is particularly true for spots generally considered as problematic due to imperfections or almost absent expression. Both PX(PAM) and PX(KMEANS) prove to be highly robust against various types of artifacts through adaptive partitioning, which more correctly assesses expression intensity values., Availability: The implementation of this method is a combination of two complementary tools Extractiff (Java) and Pixclust (free statistical language R), which are available upon request from the authors.