Transcriptional responses indicate attenuated oxidative stress in the springtail Folsomia candida exposed to mixtures of cadmium and phenanthrene.

Since the 'omics revolution', the assessment of toxic chemical mixtures has incorporated approaches where phenotypic endpoints are connected to a mechanistic understanding of toxicity. In this study we determined the effect of binary mixtures of cadmium and phenanthrene on the reproduction of Folsomia candida and investigated the cellular mechanisms underlying this response. Mixture toxicity modeling showed an antagonistic deviation from concentration addition for reproduction effects of the mixtures. Subsequent transcriptional response analysis was done using five mixtures at the modeled 50 % effect level for reproduction. The transcription profiles of 86 high throughput RT-qPCR assays were studied by means of partial least squares regression analysis. The first and second principal components (PCs) were correlated with global responses to cadmium and phenanthrene, while correlations with the mixture treatments were found in the higher PCs. Specifically associated with the mixture treatments were a biotransformation phase II gene, four mitochondrial related genes and a gene involved in the biosynthesis of antioxidant selenoproteins. Membrane integrity related gene inductions were correlated with the single phenanthrene treatment but not with the mixtures. Immune and inflammatory response assays did not correlate with any of the mixtures. These results suggest moderated oxidative stress, a higher mitochondrial maintenance and less compromised membrane function in the mixture exposed samples compared to the separate cadmium or phenanthrene exposures. The antagonism found for inhibition of reproduction may partially originate from these differences. Mechanistic studies on mixture toxicity can ultimately aid risk assessment by defining relevant toxicity pathways in organisms exposed to real-world mixture exposures present in the field.

Functional environmental genomics of a municipal landfill soil.

We investigated the toxicity of soil samples derived from a former municipal landfill site in the South of the Netherlands, where a bioremediation project is running aiming at reusing the site for recreation. Both an organic soil extract and the original soil sample was investigated using the ISO standardized Folsomia soil ecotoxicological testing and gene expression analysis. The 28 day survival/reproduction test revealed that the ecologically more relevant original soil sample was more toxic than the organic soil extract. Microarray analysis showed that the more toxic soil samples induced gene regulatory changes in twice as less genes compared to the soil extract. Consequently gene regulatory changes were highly dependent on sample type, and were to a lesser extent caused by exposure level. An important biological process shared among the two sample types was the detoxification pathway for xenobiotics (biotransformation I, II, and III) suggesting a link between compound type and observed adverse effects. Finally, we were able to retrieve a selected group of genes that show highly significant dose-dependent gene expression and thus were tightly linked with adverse effects on reproduction. Expression of four cytochrome P450 genes showed highest correlation values with reproduction, and maybe promising genetic markers for soil quality. However, a more elaborate set of environmental soil samples is needed to validate the correlation between gene expression induction and adverse phenotypic effects.

High throughput nano-liter RT-qPCR to classify soil contamination using a soil arthropod.

BACKGROUND: To incorporate genomics data into environmental assessments a mechanistic perspective of interactions between chemicals and induced biological processes needs to be developed. Since chemical compounds with structural similarity often induce comparable biological responses in exposed animals, gene expression signatures can serve as a starting point for the assessment of chemicals and their toxicity, but only when relevant and stable gene panels are available. To design such a panel, we isolated differentially expressed gene fragments from the soil arthropod Folsomia candida, a species often used for ecotoxicological testing. Animals were exposed to two chemically distinct compounds, being a metal (cadmium) and a polycyclic aromatic hydrocarbon (phenanthrene). We investigated the affected molecular responses resulting from either treatment and developed and validated 44 qPCR assays for their responses using a high throughput nano-liter RT-qPCR platform for the analysis of the samples. RESULTS: Suppressive subtractive hybridization (SSH) was used to retrieve stress-related gene fragments. SSH libraries revealed pathways involved in mitochondrial dysfunction and protein degradation for cadmium and biotransformation for phenanthrene to be overrepresented. Amongst a small cluster of SSH-derived cadmium responsive markers were an inflammatory response protein and an endo-glucanase. Conversely, cytochrome P450 family 6 or 9 was specifically induced by phenanthrene. Differential expressions of these candidate biomarkers were also highly significant in the independently generated test sample set. Toxicity levels in different training samples were not reflected by any of the markers' intensity of expressions. Though, a model based on partial least squares differential analysis (PLS-DA) (with RMSEPs between 9 and 22% and R(2)s between 0.82 and 0.97) using gene expressions of 25 important qPCR assays correctly predicted the nature of exposures of test samples. CONCLUSIONS: For the application of molecular bio-indication in environmental assessments, multivariate analyses obviously have an added value over univariate methods. Our results suggest that compound discrimination can be achieved by PLS-DA, based on a hard classification of the within-class rankings of samples from a test set. This study clearly shows that the use of high throughput RT-qPCR could be a valuable tool in ecotoxicology combining high throughput with analytical sensitivity.

Reference genes for QRT-PCR tested under various stress conditions in Folsomia candida and Orchesella cincta (Insecta, Collembola).

BACKGROUND: Genomic studies measuring transcriptional responses to changing environments and stress currently make their way into the field of evolutionary ecology and ecotoxicology. To investigate a small to medium number of genes or to confirm large scale microarray studies, Quantitative Reverse Transcriptase PCR (QRT-PCR) can achieve high accuracy of quantification when key standards, such as normalization, are carefully set. In this study, we validated potential reference genes for their use as endogenous controls under different chemical and physical stresses in two species of soil-living Collembola, Folsomia candida and Orchesella cincta. Treatments for F. candida were cadmium exposure, phenanthrene exposure, desiccation, heat shock and pH stress, and for O. cincta cadmium, desiccation, heat shock and starvation. RESULTS: Eight potential reference genes for F. candida and seven for O. cincta were ranked by their stability per stress factor using the programs geNorm and Normfinder. For F. candida the succinate dehydrogenase (SDHA) and eukaryotic transcription initiation factor 1A (ETIF) genes were found the most stable over the different treatments, while for O. cincta, the beta actin (ACTb) and tyrosine 3-monooxygenase (YWHAZ) genes were the most stable. CONCLUSION: We present a panel of reference genes for two emerging ecological genomic model species tested under a variety of treatments. Within each species, different treatments resulted in differences in the top stable reference genes. Moreover, the two species differed in suitable reference genes even when exposed to similar stresses. This might be attributed to dissimilarity of physiology. It is vital to rigorously test a panel of reference genes for each species and treatment, in advance of relative quantification of QRT-PCR gene expression measurements.

Gene expression analysis of collembola in cadmium containing soil.

Increasing concern about pollution of our environment calls for advanced and rapid methods to estimate ecological toxicity. The use of gene expression microarrays in environmental studies can potentially meet this challenge. We present a novel method to examine soil toxicity. We exposed the collembolan Folsomia candida to soil containing an ecologically relevant cadmium concentration, and found a cumulative total of 1586 differentially expressed transcripts across three exposure durations, including transcripts involved in stress response, detoxification, and hypoxia. Additional enrichment analysis of gene ontology (GO) terms revealed that antibiotic biosynthesis is important at all time points examined. Interestingly, genes involved in the "penicillin and cephalosporin biosynthesis pathway" have never been identified in animals before, but are expressed in F. candida's tissue. The synthesis of antibiotics can possibly be a response to increased cadmium-induced susceptibility to invading pathogens, which might be caused by repression of genes involved in the immune-system (C-type lectins and Toll receptor). This study presents a first global view on the environmental stress response of an arthropod species exposed to contaminated soil, and provides a mechanistic basis for the development of a gene expression soil quality test.

Collembase: a repository for springtail genomics and soil quality assessment.

BACKGROUND: Environmental quality assessment is traditionally based on responses of reproduction and survival of indicator organisms. For soil assessment the springtail Folsomia candida (Collembola) is an accepted standard test organism. We argue that environmental quality assessment using gene expression profiles of indicator organisms exposed to test substrates is more sensitive, more toxicant specific and significantly faster than current risk assessment methods. To apply this species as a genomic model for soil quality testing we conducted an EST sequencing project and developed an online database. DESCRIPTION: Collembase is a web-accessible database comprising springtail (F. candida) genomic data. Presently, the database contains information on 8686 ESTs that are assembled into 5952 unique gene objects. Of those gene objects approximately 40% showed homology to other protein sequences available in GenBank (blastx analysis; non-redundant (nr) database; expect-value < 10-5). Software was applied to infer protein sequences. The putative peptides, which had an average length of 115 amino-acids (ranging between 23 and 440) were annotated with Gene Ontology (GO) terms. In total 1025 peptides (approximately 17% of the gene objects) were assigned at least one GO term (expect-value < 10-25). Within Collembase searches can be conducted based on BLAST and GO annotation, cluster name or using a BLAST server. The system furthermore enables easy sequence retrieval for functional genomic and Quantitative-PCR experiments. Sequences are submitted to GenBank (Accession numbers: EV473060 - EV481745). CONCLUSION: Collembase http://www.collembase.org is a resource of sequence data on the springtail F. candida. The information within the database will be linked to a custom made microarray, based on the Agilent platform, which can be applied for soil quality testing. In addition, Collembase supplies information that is valuable for related scientific disciplines such as molecular ecology, ecogenomics, molecular evolution and phylogenetics.