Increased comparability between RNA-Seq and microarray data by utilization of gene sets.

The field of transcriptomics uses and measures mRNA as a proxy of gene expression. There are currently two major platforms in use for quantifying mRNA, microarray and RNA-Seq. Many comparative studies have shown that their results are not always consistent. In this study we aim to find a robust method to increase comparability of both platforms enabling data analysis of merged data from both platforms. We transformed high dimensional transcriptomics data from two different platforms into a lower dimensional, and biologically relevant dataset by calculating enrichment scores based on gene set collections for all samples. We compared the similarity between data from both platforms based on the raw data and on the enrichment scores. We show that the performed data transforms the data in a biologically relevant way and filters out noise which leads to increased platform concordance. We validate the procedure using predictive models built with microarray based enrichment scores to predict subtypes of breast cancer using enrichment scores based on sequenced data. Although microarray and RNA-Seq expression levels might appear different, transforming them into biologically relevant gene set enrichment scores significantly increases their correlation, which is a step forward in data integration of the two platforms. The gene set collections were shown to contain biologically relevant gene sets. More in-depth investigation on the effect of the composition, size, and number of gene sets that are used for the transformation is suggested for future research.

A novel toxicogenomics-based approach to categorize (non-)genotoxic carcinogens.

Alternative methods to detect non-genotoxic carcinogens are urgently needed, as this class of carcinogens goes undetected in the current testing strategy for carcinogenicity under REACH. A complicating factor is that non-genotoxic carcinogens act through several distinctive modes of action, which makes prediction of their carcinogenic property difficult. We have recently demonstrated that gene expression profiling in primary mouse hepatocytes is a useful approach to categorize non-genotoxic carcinogens according to their modes of action. In the current study, we improved the methods used for analysis and added mouse embryonic stem cells as a second in vitro test system, because of their features complementary to hepatocytes. Our approach involved an unsupervised analysis based on the 30 most significantly up- and down-regulated genes per chemical. Mouse embryonic stem cells and primary mouse hepatocytes were exposed to a selected set of chemicals and subsequently subjected to gene expression profiling. We focused on non-genotoxic carcinogens, but also included genotoxic carcinogens and non-carcinogens to test the robustness of this approach. Application of the optimized comparison approach resulted in improved categorization of non-genotoxic carcinogens. Mouse embryonic stem cells were a useful addition, especially for genotoxic substances, but also for detection of non-genotoxic carcinogens that went undetected by primary hepatocytes. The approach presented here is an important step forward to categorize chemicals, especially those that are carcinogenic.

Staphylococcus aureus ST398 gene expression profiling during ex vivo colonization of porcine nasal epithelium.

BACKGROUND: Staphylococcus aureus is a common human and animal opportunistic pathogen. In humans nasal carriage of S. aureus is a risk factor for various infections. Methicillin-resistant S. aureus ST398 is highly prevalent in pigs in Europe and North America. The mechanism of successful pig colonization by MRSA ST398 is poorly understood. Previously, we developed a nasal colonization model of porcine nasal mucosa explants to identify molecular traits involved in nasal MRSA colonization of pigs. RESULTS: We report the analysis of changes in the transcription of MRSA ST398 strain S0462 during colonization on the explant epithelium. Major regulated genes were encoding metabolic processes and regulation of these genes may represent metabolic adaptation to nasal mucosa explants. Colonization was not accompanied by significant changes in transcripts of the main virulence associated genes or known human colonization factors. Here, we documented regulation of two genes which have potential influence on S. aureus colonization; cysteine extracellular proteinase (scpA) and von Willebrand factor-binding protein (vWbp, encoded on SaPIbov5). Colonization with isogenic-deletion strains (Δvwbp and ΔscpA) did not alter the ex vivo nasal S. aureus colonization compared to wild type. CONCLUSIONS: Our results suggest that nasal colonization with MRSA ST398 is a complex event that is accompanied with changes in bacterial gene expression regulation and metabolic adaptation.

In vivo murine hepatic microRNA and mRNA expression signatures predicting the (non-)genotoxic carcinogenic potential of chemicals.

There is a high need to improve the assessment of, especially non-genotoxic, carcinogenic features of chemicals. We therefore explored a toxicogenomics-based approach using genome-wide microRNA and mRNA expression profiles upon short-term exposure in mice. For this, wild-type mice were exposed for seven days to three different classes of chemicals, i.e., four genotoxic carcinogens (GTXC), seven non-genotoxic carcinogens (NGTXC), and five toxic non-carcinogens. Hepatic expression patterns of mRNA and microRNA transcripts were determined after exposure and used to assess the discriminative power of the in vivo transcriptome for GTXC and NGTXC. A final classifier set, discriminative for GTXC and NGTXC, was generated from the transcriptomic data using a tiered approach. This appeared to be a valid approach, since the predictive power of the final classifier set in three different classifier algorithms was very high for the original training set of chemicals. Subsequent validation in an additional set of chemicals revealed that the predictive power for GTXC remained high, in contrast to NGTXC, which appeared to be more troublesome. Our study demonstrated that the in vivo microRNA-ome has less discriminative power to correctly identify (non-)genotoxic carcinogen classes. The results generally indicate that single mRNA transcripts do have the potential to be applied in risk assessment, but that additional (genomic) strategies are necessary to correctly predict the non-genotoxic carcinogenic potential of a chemical.

Dissecting modes of action of non-genotoxic carcinogens in primary mouse hepatocytes.

Under REACH, the European Community Regulation on chemicals, the testing strategy for carcinogenicity is based on in vitro and in vivo genotoxicity assays. Given that non-genotoxic carcinogens are negative for genotoxicity and chronic bioassays are no longer regularly performed, this class of carcinogens will go undetected. Therefore, test systems detecting non-genotoxic carcinogens, or even better their modes of action, are required. Here, we investigated whether gene expression profiling in primary hepatocytes can be used to distinguish different modes of action of non-genotoxic carcinogens. For this, primary mouse hepatocytes were exposed to 16 non-genotoxic carcinogens with diverse modes of action. Upon profiling, pathway analysis was performed to obtain insight into the biological relevance of the observed changes in gene expression. Subsequently, both a supervised and an unsupervised comparison approach were applied to recognize the modes of action at the transcriptomic level. These analyses resulted in the detection of three of eight compound classes, that is, peroxisome proliferators, metalloids and skin tumor promotors. In conclusion, gene expression profiles in primary hepatocytes, at least in rodent hepatocytes, appear to be useful to detect some, certainly not all, modes of action of non-genotoxic carcinogens.

Finding transcriptomics biomarkers for in vivo identification of (non-)genotoxic carcinogens using wild-type and Xpa/p53 mutant mouse models.

The carcinogenic potential of chemicals and pharmaceuticals is traditionally tested in the chronic, 2 year rodent bioassay. This assay is not only time consuming, expensive and often with a limited sensitivity and specificity but it also causes major distress to the experimental animals. A major improvement in carcinogenicity testing, especially regarding reduction and refinement of animal experimentation, could be the application of toxicogenomics. The ultimate aim of this study is to demonstrate a proof-of-principle for transcriptomics biomarkers in various tissues for identification of (subclasses of) carcinogenic compounds after short-term in vivo exposure studies. Both wild-type and DNA repair-deficient Xpa(-/-)/p53(+/-) (Xpa/p53) mice were exposed up to 14 days to compounds of three distinct classes: genotoxic carcinogens (GTXC), non-genotoxic carcinogens (NGTXC) and non-carcinogens. Subsequently, extensive transcriptomics analyses were performed on several tissues, and transcriptomics data were screened for potential biomarkers using advanced statistical learning techniques. For all tissues analyzed, we identified multigene gene-expression signatures that are, with a high confidence, predictive for GTXC and NGTXC exposures in both mouse genotypes. Xpa/p53 mice did not perform better in the short-term bioassay. We were able to achieve a proof-of-principle for the identification and use of transcriptomics biomarkers for GTXC or NGTXC. This supports the view that toxicogenomics with short-term in vivo exposure provides a viable tool for classifying (geno)toxic compounds.

The absence of Ser389 phosphorylation in p53 affects the basal gene expression level of many p53-dependent genes and alters the biphasic response to UV exposure in mouse embryonic fibroblasts.

Phosphorylation is important in p53-mediated DNA damage responses. After UV irradiation, p53 is phosphorylated specifically at murine residue Ser389. Phosphorylation mutant p53.S389A cells and mice show reduced apoptosis and compromised tumor suppression after UV irradiation. We investigated the underlying cellular processes by time-series analysis of UV-induced gene expression responses in wild-type, p53.S389A, and p53(-/-) mouse embryonic fibroblasts. The absence of p53.S389 phosphorylation already causes small endogenous gene expression changes for 2,253, mostly p53-dependent, genes. These genes showed basal gene expression levels intermediate to the wild type and p53(-/-), possibly to readjust the p53 network. Overall, the p53.S389A mutation lifts p53-dependent gene repression to a level similar to that of p53(-/-) but has lesser effect on p53-dependently induced genes. In the wild type, the response of 6,058 genes to UV irradiation was strictly biphasic. The early stress response, from 0 to 3 h, results in the activation of processes to prevent the accumulation of DNA damage in cells, whereas the late response, from 12 to 24 h, relates more to reentering the cell cycle. Although the p53.S389A UV gene response was only subtly changed, many cellular processes were significantly affected. The early response was affected the most, and many cellular processes were phase-specifically lost, gained, or altered, e.g., induction of apoptosis, cell division, and DNA repair, respectively. Altogether, p53.S389 phosphorylation seems essential for many p53 target genes and p53-dependent processes.

Overlapping gene expression profiles of model compounds provide opportunities for immunotoxicity screening.

In order to investigate immunotoxic effects of a set of model compounds in mice, a toxicogenomics approach was combined with information on macroscopical and histopathological effects on spleens and on modulation of immune function. Bis(tri-n-butyltin)oxide (TBTO), cyclosporin A (CsA), and benzo[a]pyrene (B[a]P) were administered to C57BL/6 mice at immunosuppressive dose levels. Acetaminophen (APAP) was included in the study since indications of immunomodulating properties of this compound have appeared in the literature. TBTO exposure caused the most pronounced effect on gene expression and also resulted in the most severe reduction of body weight gain and induction of splenic irregularities. All compounds caused inhibition of cell division in the spleen as shown by microarray analysis as well as by suppression of lymphocyte proliferation after application of a contact sensitizer as demonstrated in an immune function assay that was adapted from the local lymph node assay. The immunotoxicogenomics approach applied in this study thus pointed to immunosuppression through cell cycle arrest as a common mechanism of action of immunotoxicants, including APAP. Genes related to cell division such as Ccna2, Brca1, Birc5, Incenp, and Cdkn1a (p21) were identified as candidate genes to indicate anti-proliferative effects of xenobiotics in immune cells for future screening assays. The results of our experiments also show the value of group wise pathway analysis for detection of more subtle transcriptional effects and the potency of evaluation of effects in the spleen to demonstrate immunotoxicity.

Primary nasal epithelium exposed to house dust mite extract shows activated expression in allergic individuals.

Nasal epithelial cells form the outermost protective layer against environmental factors. However, this defense is not just physical; it has been shown that epithelial cells respond by the production of inflammatory mediators that may affect local immune responses. In this research we set out to characterize potential differences between the responses of nasal epithelium from healthy and allergic individuals to house dust mite (HDM) allergen. These differences will help us to define local mechanisms that could contribute to allergic disease expression. Epithelial cells were cultured from nasal biopsies taken from five healthy and five allergic individuals. These cultures were exposed for 24 hours to culture medium containing HDM allergen, or to culture medium alone. Isolated RNA was used for microarray analysis. Gene-ontology of the response in healthy epithelium revealed mainly up-regulation of chemokines, growth factors, and structural proteins. Moreover, we saw increased expression of two transcription factors (NF-kappaB and AP-1) and their regulatory members. The expression pattern of epithelium from allergic individuals in the absence of the HDM stimulus suggests that it is already in an activated state. Most striking is that, while the already activated NF-kappaB regulatory pathway remained unchanged in allergic epithelium, the AP-1 pathway is down-regulated upon exposure to HDM allergen; this is contrary to what we see in healthy epithelium. Clear differences in the expression pattern exist between epithelial cells isolated from healthy and allergic individuals at baseline and between their responses to allergen exposure; these differences may contribute to the inflammatory response.

Delayed expression of apoptotic and cell-cycle control genes in carcinogen-exposed bladders of mice lacking p53.S389 phosphorylation.

Mice with non-phosphorylated serine 389 in p53 are susceptible for bladder tumors induced by 2-acetylaminofluorene (2-AAF). Since p53 is a transcription factor, this might well be preceded by differences in the regulation of gene expression. Microarray analysis was used to determine early transcriptional changes that might underlie this cancer-prone phenotype. Interestingly, lack of Ser389 phosphorylation led to endogenously different gene expression levels. The number of genes affected was, however, rather small. Conversely, after short-term exposure to 2-AAF, wild-type and p53.S389A bladders demonstrated a significant number of differentially expressed genes. Differences between wild-type and p53.S389A could mainly be attributed to a delayed, rather than complete absence of, transcriptional response of a group of genes, including well-known p53 target genes involved in apoptosis and cell-cycle control like Bax, Perp and P21. An analysis of differentially expressed genes in non-tumorigenic tissue and bladder tumors of p53.S389A after long-term exposure to 2-AAF revealed 319 genes. Comparison of these with those found after short-term exposure resulted in 23 transcripts. These possible marker genes might be useful for the early prediction of bladder tumor development. In conclusion, our data indicate that lack of Ser389 phosphorylation results in aberrant expression of genes needed to execute vital responses to DNA damage. Post-translational modifications, like Ser389 phosphorylation, seem crucial for fine-tuning the transcription of a specific set of genes and do not appear to give rise to major changes in transcription patterns. As such, Ser389 phosphorylation is needed for some, but certainly not all, p53 functions.

Toxicity of binary mixtures of cadmium-copper and carbendazim-copper to the nematode Caenorhabditis elegans.

For ecological risk assessment, the additive model may be used to empirically predict toxic mixture effects. Detailed toxicity tests were performed to determine whether effects of mixtures of copper-cadmium and copper-carbendazim on Caenorhabditis elegans were similar to the effects of the individual compounds. Effects on the course of reproduction, the length of the juvenile period, the length of the reproductive period, and body length were analyzed. Dose-response data were compared to the additive model and tested for four deviation patterns from additivity: No deviation, synergistic/antagonistic deviation, dose ratio-dependent deviation, dose level-dependent deviation. During the exposure, the cadmium-copper effect on reproduction changed from a synergistic, to a dose ratio-dependent deviation from additivity. More cadmium in the mixture decreased the toxicity and more copper increased the toxicity. The effect of copper-carbendazim on reproduction was synergistic at low dose levels and antagonistic at high dose levels and independent of time. Mixture effects on the juvenile and reproductive period were similar to single component effects. It was concluded that the observed time-dependence of toxic interactions was small and that interactions on the timing of reproduction were not found. The additive model underestimated mixture effects on reproduction and body length.