Using transcriptomic profiles in the diatom Phaeodactylum tricornutum to identify and prioritize stressors

Potential cellular targets of platinum in the freshwater microalgae Chlamydomonas reinhardtii and Nitzschia palea revealed by transcriptomics

Platinum group element levels have increased in natural aquatic environments in the last few decades, in particular as a consequence of the use of automobile catalytic converters on a global scale. Concentrations of Pt over tens of μg L-1 have been observed in rivers and effluents. This raises questions regarding its possible impacts on aquatic ecosystems, as Pt natural background concentrations are extremely low to undetectable. Primary producers, such as microalgae, are of great ecological importance, as they are at the base of the food web. The purpose of this work was to better understand the impact of Pt on a cellular level for freshwater unicellular algae. Two species with different characteristics, a green alga C. reinhardtii and a diatom N. palea, were studied. The bioaccumulation of Pt as well as its effect on growth were quantified. Moreover, the induction or repression factors of 16 specific genes were determined and allowed for the determination of possible intracellular effects and pathways of Pt. Both species seemed to be experiencing copper deficiency as suggested by inductions of genes linked to copper transporters. This is an indication that Pt might be internalized through the Cu(I) metabolic pathway. Moreover, Pt could possibly be excreted using an efflux pump. Other highlights include a concentration-dependent negative impact of Pt on mitochondrial metabolism for C. reinhardtii which is not observed for N. palea. These findings allowed for a better understanding of some of the possible impacts of Pt on freshwater primary producers, and also lay the foundations for the investigation of pathways for Pt entry at the base of the aquatic food web.

Toxicity of herbicides to the marine microalgae Tisochrysis lutea and Tetraselmis sp

Pesticides are ubiquitous in the catchments of the Great Barrier Reef (GBR) and regularly discharge into the nearshore waters. Effective management of pesticides requires suitable water quality guideline values (WQGVs), and further ecotoxicological data for many pesticides are needed to improve the reliability of environmental risk assessments. To help address this issue, toxicity thresholds were determined to two species of tropical marine microalgae Tisochrysis lutea and Tetraselmis sp. for a suite of herbicides detected in the GBR. Photosystem II (PSII) herbicides significantly reduced growth with no effect concentration (NEC) and 10% effect concentration (EC10) values spanning two orders of magnitude from 0.60 µg L-1 for diuron to 60 µg L-1 for simazine across both species. However, growth was insensitive to the non-PSII herbicides. The NEC/EC10 thresholds for most herbicide-microalgae combinations were greater than recent WQGVs intended to protect 99% of species (PC99); however, metribuzin was toxic to T. lutea at concentrations lower than the current PC99 value, which may have to be revisited. The toxicity thresholds for alternative herbicides derived here further inform the development of national and GBR-specific WQGVs, but more toxicity data is needed to develop WQGVs for the > 50 additional pesticides detected in catchments of the GBR.

Copper toxicity leads to accumulation of free amino acids and polyphenols in Phaeodactylum tricornutum diatoms

This work is focused on the effect of lethal and sub-lethal copper (Cu) concentrations on the free amino acid and polyphenol production by the marine diatom Phaeodactylum tricornutum (P. tricornutum) after 12, 18, and 21 days of exposure. The concentrations of 10 amino acids (arginine, aspartic acid, glutamic acid, histidine, lysine, methionine, proline, valine, isoleucine, and phenylalanine) and 10 polyphenols (gallic acid, protocatechuic acid, p-coumaric acid, ferulic acid, catechin, vanillic acid, epicatechin syringic acid, rutin, and gentisic acid) were measured by RP-HPLC. Under lethal doses of Cu, free amino acids reached levels significantly higher than those in the control cells (up to 21.9 times higher), where histidine and methionine showed the highest increases (up to 37.4 and 65.8 times higher, respectively). The total phenolic content also increased up to 11.3 and 5.59 times higher compared to the reference cells, showing gallic acid the highest increase (45.8 times greater). The antioxidant activities of cells exposed to Cu were also enhanced with increasing doses of Cu(II). They were evaluated by 2,2-diphenyl-1-picrylhydrazyl (DPPH) free radical scavenging ability (RSA), cupric ion reducing antioxidant capacity (CUPRAC), and ferric reducing antioxidant power (FRAP) assays. Malonaldehyde (MDA) exhibited the same tendency: cells grown at the highest lethal Cu concentration yielded the highest MDA level. These findings reflect the involvement of amino acids and polyphenols in protective mechanisms to overcome the toxicity of copper in marine microalgae.

A metatranscriptomic analysis of changing dynamics in the plankton communities adjacent to aquaculture leases in southern Tasmania, Australia

Aquaculture releases nitrogen to the marine environment, potentially changing dynamics of local plankton populations and causing adverse impacts. Metatranscriptomics have been used to study planktonic nutrient cycles and community dynamics. We hypothesised that the metatranscriptome could be used to monitor changing phytoplankton physiology near leases. To test this hypothesis, opportunistic samples were collected from one oceanic location in winter and one estuarine location in spring and analysed via RNASeq. Transcriptomes from different locations were found to have little overlap, due to different community compositions in the oceanic and estuarine locations. Transcript function was similar at each location. Proximity to the salmon pen had little influence over the transcriptome at the estuarine location. In the oceanic environment, diatom-based activity decreased near pens and dinoflagellate-based activity increased as demonstrated through the abundance of carbon fixation and nitrogen-acquisition-related transcripts. Our initial results suggest that the use of the metatranscriptome in monitoring is promising.

Development of High-Level Omega-3 Eicosapentaenoic Acid (EPA) Production from Phaeodactylum tricornutum

Phaeodactylum tricornutum is a lipid-rich marine diatom that contains a high level of omega-3 polyunsaturated fatty acids, especially eicosapentaenoic acid (EPA). In an effort to reduce costs for large-scale cultivation of this microalga, this study first established a New BBM medium (0.3 x strength BBM with only 3% of the initial phosphate level) to replace the traditional F/2 medium. Phaeodactylum tricornutum could grow in extremely low phosphate concentrations (25 µM), without compromising the EPA content. In the presence of sea salts, silicate addition was not necessary for high rate growth, high EPA content, or lipid accumulation in this species. Using urea as the sole nitrogen source tended to increase EPA contents per dry biomass (by 24.7%) while not affecting growth performance. The use of sea salts, rather than just sodium chloride, led to significantly improved biomass yields (20% increase) and EPA contents of total fatty acid (46-52% increase), most likely because it supplied sufficient essential elements such as magnesium. A salinity level of 35 led to significantly higher biomass yields compared with 20, but salinity had no significant influence on EPA content. EPA became the dominant fatty acid with average levels of 51.8% of total fatty acids during the exponential growth phase at 20 ppt in New BBM medium with sea salts.

Toxicity of ten herbicides to the tropical marine microalgae Rhodomonas salina

Herbicide contamination of nearshore tropical marine ecosystems is widespread and persistent; however, risks posed by most ‘alternative’ herbicides to tropical marine microalgae remain poorly understood. Experimental exposures of the important but understudied microalgae Rhodomonas salina to seven individual Photosystem II (PSII) inhibitor herbicides (diuron, metribuzin, hexazinone, tebuthiuron, bromacil, simazine, propazine) led to inhibition of effective quantum yield (ΔF/F_m′) and subsequent reductions in specific growth rates (SGR). The concentrations which reduced ΔF/F_m′ by 50% (EC₅₀) ranged from 1.71-59.2 µg L⁻¹, while the EC₅₀s for SGR were 4-times higher, ranging from 6.27-188 µg L⁻¹. Inhibition of ΔF/F_m′ indicated reduced photosynthetic capacity, and this correlated linearly with reduced SGR (R² = 0.89), supporting the application of ∆F/F_m’ inhibition as a robust and sensitive indicator of sub-lethal toxicity of PSII inhibitors for this microalga. The three non-PSII inhibitor herbicides (imazapic, haloxyfop and 2,4-Dichlorophenoxyacetic acid (2,4-D)) caused low or no toxic responses to the function of the PSII or growth at the highest concentrations tested suggesting these herbicides pose little risk to R. salina. This study highlights the suitability of including R. salina in future species sensitivity distributions (SSDs) to support water quality guideline development for the management of herbicide contamination in tropical marine ecosystems.

Carbon fixation gene expression in Skeletonema marinoi in nitrogen-, phosphate-, silicate-starvation, and low-temperature stress exposure

Diatoms are unicellular algae with a set of extraordinary genes, metabolic pathways, and physiological functions acquired by secondary endosymbiosis, especially for their efficient photosynthetic carbon fixation mechanisms, which can be a reason for their successful environmental adaptation and great contribution to primary production. Based on the available genomic information, the expression patterns of carbon fixation genes were analyzed using transcriptomic sequencing and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) in Skeletonema marinoi. Meanwhile, suitable reference genes applying to specific experimental treatments were selected. In our results, carbon fixation genes were standardized by actin and TATA box-binding protein-coding genes in growth phase samples and stress conditions, respectively. It was found that a series of carbon fixation genes, such as the pyruvate orthophosphate dikinase (PPDK)-coding gene, had significantly up-regulated expression in nitrogen-starvation, phosphate-starvation, and low-temperature conditions, but consistently down-regulated in silicate-starvation treatment. These carbon fixation genes exhibited variable expression levels in different conditions and will be useful for investigating gene expression mechanisms in S. marinoi and improve our understanding of diatom carbon fixation pathways.

PhaeoNet: A Holistic RNAseq-Based Portrait of Transcriptional Coordination in the Model Diatom Phaeodactylum tricornutum

Transcriptional coordination is a fundamental component of prokaryotic and eukaryotic cell biology, underpinning the cell cycle, physiological transitions, and facilitating holistic responses to environmental stress, but its overall dynamics in eukaryotic algae remain poorly understood. Better understanding of transcriptional partitioning may provide key insights into the primary metabolism pathways of eukaryotic algae, which frequently depend on intricate metabolic associations between the chloroplasts and mitochondria that are not found in plants. Here, we exploit 187 publically available RNAseq datasets generated under varying nitrogen, iron and phosphate growth conditions to understand the co-regulatory principles underpinning transcription in the model diatom Phaeodactylum tricornutum. Using WGCNA (Weighted Gene Correlation Network Analysis), we identify 28 merged modules of co-expressed genes in the P. tricornutum genome, which show high connectivity and correlate well with previous microarray-based surveys of gene co-regulation in this species. We use combined functional, subcellular localization and evolutionary annotations to reveal the fundamental principles underpinning the transcriptional co-regulation of genes implicated in P. tricornutum chloroplast and mitochondrial metabolism, as well as the functions of diverse transcription factors underpinning this co-regulation. The resource is publically available as PhaeoNet, an advanced tool to understand diatom gene co-regulation.

Salinity-dependent nanostructures and composition of cell surface and its relation to Cd toxicity in an estuarine diatom

The interactions between metal and phytoplankton are affected by salinity in estuarine environments. While water chemistry is an important factor regulating the metal bioavailability in phytoplankton, the physiological adaptation of the algae cells may also change their intrinsic response to metals. In this work, we tried to interpret the salinity-dependent Cd toxicity in a pennate diatom Nitzchia closterium from a biological side. As with many studies, we observed Cd toxicity to the diatom increased with decreasing salinity. However, changing free Cd ion concentrations may be partly responsible for the enhanced Cd toxicity. Multiple evidences showed that diatom cells acclimated at low salinity had stronger intrinsic Cd adsorption capacity. Salinity significantly affected not only the nanostructures but also the biochemical composition in the cell surface of the diatom. Diatom cells grown at lower salinity had a lower surface potential, higher specific surface area, and more sulfur-containing groups in the cell wall, leading to stronger Cd binding capacity in the cells. Meanwhile, more Si was present as poly-silicic acid when the salinity decreased. The change of Si content and speciation in the cell wall are also considered a major reason for the variations of Cd surface binding. Our study provided new clues for the salinity-dependent metal toxicity in marine diatoms.

Lipid metabolism and pro-oxidant/antioxidant balance of Halamphora oceanica from the Gulf of Mexico exposed to water accommodated fraction of Maya crude oil

Diatoms play key roles in primary production and carbon fixation at a global scale and in some cases these species live on marine ecosystems impacted by crude oil (CO) spills. Halamphora oceanica, a new diatom species from the Southwest of the Gulf of Mexico was isolated and cultured in the laboratory and was exposed to water accommodated fraction (WAF) of different Maya CO loads at 0.01, 0.1, 1 and 10g/L by 96h. A battery of biomarkers involved in oxidative stress (O₂•, H₂O₂, TBARS, ROOH, RC=O, SOD, CAT, GPx), biotransformation and conjugation (total CYP450 activity and GST) moreover fatty acid (FA) metabolism (FA levels, fatty-acid synthase and acyl-CoA oxidase) were measured. Obtained results suggest that increases of PAHs in the medium (below to EC₅₀) acts as external forces able to turn-on regulatory mechanisms on H. oceanica involved in both, on the PAHs uptake and changing its aerobic metabolism to anaerobic metabolism. However, the growth of this microalgae species evaluated as chlorophyll "a" and pheophytin levels increased as the WAF concentration indicating that PAHs and other hydrosoluble hydrocarbons were used as carbon and energy sources by unidentified enzymes not evaluated in the current study. Our hypothesis was also corroborated by IBRv2. In the current study, we suppose the change from aerobic to anaerobic metabolism as a strategy for Halamphora oceanica survival exposed to petroleum hydrocarbons.

The Combination of RNA and Protein Profiling Reveals the Response to Nitrogen Depletion in Thalassiosira pseudonana

Nitrogen (N) is essential for the growth of algae, and its concentration varies greatly in the ocean, which has been regarded as a limitation for phytoplankton growth. Despite its great importance, most of the existing studies on the mechanisms underlying the effects of N on diatoms have focused on physiology, biochemistry and a few target genes and have rarely involved whole genomic analyses. Therefore, in this study, we integrated physiological data with RNA and protein profiling data to reveal the response strategy of Thalassiosira pseudonana under N-depleted conditions. Physiological measurements indicated that the cell growth capacity and chlorophyll content of the cells decreased, as did the expression of photosynthesis- and chlorophyll biosynthesis-related genes or proteins. The RNA-Seq profile results showed that T. pseudonana responded to N deprivation through increases in glycolysis, the TCA cycle and N metabolism as well as down-regulation in the Calvin cycle, gluconeogenesis, pentose phosphate, oxidative phosphorylation and lipid synthesis. These results provide a basic understanding for further research addressing how N affects phytoplankton in terms of genomics.

A Community Multi-Omics Approach towards the Assessment of Surface Water Quality in an Urban River System

A multi-omics approach was applied to an urban river system (the Brisbane River (BR), Queensland, Australia) in order to investigate surface water quality and characterize the bacterial population with respect to water contaminants. To do this, bacterial metagenomic amplicon-sequencing using Illumina next-generation sequencing (NGS) of the V5-V6 hypervariable regions of the 16S rRNA gene and untargeted community metabolomics using gas chromatography coupled with mass spectrometry (GC-MS) were utilized. The multi-omics data, in combination with fecal indicator bacteria (FIB) counts, trace metal concentrations (by inductively coupled plasma mass spectrometry (ICP-MS)) and in-situ water quality measurements collected from various locations along the BR were then used to assess the health of the river ecosystem. Sites sampled represented the transition from less affected (upstream) to polluted (downstream) environments along the BR. Chemometric analysis of the combined datasets indicated a clear separation between the sampled environments. Burkholderiales and Cyanobacteria were common key factors for differentiation of pristine waters. Increased sugar alcohol and short-chain fatty acid production was observed by Actinomycetales and Rhodospirillaceae that are known to form biofilms in urban polluted and brackish waters. Results from this study indicate that a multi-omics approach enables a deep understanding of the health of an aquatic ecosystem, providing insight into the bacterial diversity present and the metabolic output of the population when exposed to environmental contaminants.

Proteome response of Phaeodactylum tricornutum, during lipid accumulation induced by nitrogen depletion

Nitrogen stress is a common strategy employed to stimulate lipid accumulation in microalgae, a biofuel feedstock of topical interest. Although widely investigated, the underlying mechanism of this strategy is still poorly understood. We examined the proteome response of lipid accumulation in the model diatom, Phaeodactylum tricornutum (CCAP 1055/1), at an earlier stage of exposure to selective nitrogen exclusion than previously investigated, and at a time point when changes would reflect lipid accumulation more than carbohydrate accumulation. In total 1043 proteins were confidently identified (≥ 2 unique peptides) with 645 significant (p < 0.05) changes observed, in the LC-MS/MS based iTRAQ investigation. Analysis of significant changes in KEGG pathways and individual proteins showed that under nitrogen starvation P. tricornutum reorganizes its proteome in favour of nitrogen scavenging and reduced lipid degradation whilst rearranging the central energy metabolism that deprioritizes photosynthetic pathways. By doing this, this species appears to increase nitrogen availability inside the cell and limit its use to the pathways where it is needed most. Compared to previously published proteomic analysis of nitrogen starvation in Chlamydomonas reinhardtii, central energy metabolism and photosynthesis appear to be affected more in the diatom, whilst the green algae appears to invest its energy in reorganizing respiration and the cellular organization pathways.

Challenges for using quantitative PCR test batteries as a TIE-type approach to identify metal exposure in benthic invertebrates

The epibenthic amphipod Melita plumulosa shows unique gene expression profiles when exposed to different contaminants. We hypothesized that specific changes in transcript abundance could be used in a battery of quantitative polymerase chain reaction (qPCR) assays as a toxicity identification evaluation (TIE)-like approach to identify the most relevant stressor in field-contaminated sediments. To test this hypothesis, seven candidate transcriptomic markers were selected, and their specificity following metal exposure was confirmed. The performance of these markers across different levels of added metals was verified. The ability of these transcripts to act as markers was tested by exposing amphipods to metal-contaminated field-collected sediments and measuring changes in transcript abundance via qPCR. For two of the three sediments tested, at least some of the transcriptomic patterns matched our predictions, suggesting that they would be effective in helping to identify metal exposure in field sediments. However, following exposure to the third sediment, transcriptomic patterns were unlike our predictions. These results suggest that the seven transcripts may be insufficient to discern individual contaminants from complex mixtures and that microarray or RNA-Seq global gene expression profiles may be more effective for TIE. Changes in transcriptomics based on laboratory exposures to single compounds should be carefully validated before the results are used to analyze mixtures.

Proteomics studies on stress responses in diatoms

Diatoms are a highly diverse group of eukaryotic phytoplankton that are distributed throughout marine and freshwater environments and are believed to be responsible for approximately 40% of the total marine primary productivity. The ecological success of diatoms suggests that they have developed a range of strategies to cope with various biotic and abiotic stress factors. It is of great interest to understand the adaptive responses of diatoms to different stresses in the marine environment. Proteomic technologies have been applied to the adaptive responses of marine diatoms under different growth conditions in recent years such as nitrogen starvation, iron limitation and phosphorus deficiency. These studies have provided clues to elucidate the sophisticated sensing mechanisms that control their adaptive responses. Although only a very limited number of proteomic studies were conducted in diatoms, the obtained data have led to a better understanding of the biochemical processes that contribute to their ecological success. This review presents the current status of proteomic studies of diatom stress responses and discusses the novel developments and applications for the analysis of protein post-translational modification in diatoms. The potential future application of proteomics could contribute to a better understanding of the physiological mechanisms underlying diatom acclimation to a given stress and the acquisition of an enhanced diatom stress tolerance. Future challenges and research opportunities in the proteomics studies of diatoms are also discussed.

454 pyrosequencing-based analysis of gene expression profiles in the amphipod Melita plumulosa: transcriptome assembly and toxicant induced changes

Next generation sequencing using Roche's 454 pyrosequencing platform can be used to generate genomic information for non-model organisms, although there are bioinformatic challenges associated with these studies. These challenges are compounded by a lack of a standardized protocol to either assemble data or to evaluate the quality of a de novo transcriptome. This study presents an assembly of the control and toxicant responsive transcriptome of Melita plumulosa, an Australian amphipod commonly used in ecotoxicological studies. RNA was harvested from control amphipods, juvenile amphipods, and from amphipods exposed to either metal or diesel contaminated sediments. This RNA was used as the basis for a 454 based transcriptome sequencing effort. Sequencing generated 1.3 million reads from control, juvenile, metal-exposed and diesel-exposed amphipods. Different read filtering and assembly protocols were evaluated to generate an assembly that (i) had an optimal number of contigs; (ii) had long contigs; (iii) contained a suitable representation of conserved genes; and (iv) had long ortholog alignment lengths relative to the length of each contig. A final assembly, generated using fixed-length trimming based on the sequence quality scores, followed by assembly using the MIRA algorithm, produced the best results. The 26,625 contigs generated via this approach were annotated using Blast2GO, and the differential expression between treatments and control was determined by mapping with BWA followed by DESeq. Although the mapping generated low coverage, many differentially expressed contigs, including some with known developmental or toxicological function, were identified. This study demonstrated that 454 pyrosequencing is an effective means of generating reference transcriptome information for organisms, such as the amphipod M. plumulosa, that have no genomic information available in databases or in closely related sequenced species. It also demonstrated how optimization of read filtering protocols and assembly approaches changes the utility of results obtained from next generation sequencing studies, and establishes criteria to determine the quality of a de novo assembly in species lacking a reference genome. This new transcriptomic knowledge provides the genomic foundation for the creation of microarray and qPCR assays, serving as a reference transcriptome in future RNAseq studies, and allowing both the biology and ecotoxicology of this organism to be better understood. This approach will allow genomics-based methodology to be applied to a wider range of environmentally relevant species.

The role of biomarkers in the assessment of aquatic ecosystem health

Ensuring the health of aquatic ecosystems and identifying species at risk from the detrimental effects of environmental contaminants can be facilitated by integrating analytical chemical analysis with carefully selected biological endpoints measured in tissues of species of concern. These biological endpoints include molecular, biochemical, and physiological markers (i.e., biomarkers) that when integrated, can clarify issues of contaminant bioavailability, bioaccumulation, and ecological effects while enabling a better understanding of the effects of nonchemical stressors. In the case of contaminant stressors, an understanding of chemical modes of toxicity can be incorporated with diagnostic markers of aquatic animal physiology to help understand the health status of aquatic organisms in the field. Furthermore, new approaches in functional genomics and bioinformatics can help discriminate individual chemicals, or groups of chemicals among complex mixtures that may contribute to adverse biological effects. Although the use of biomarkers is not a new paradigm, such approaches have been underused in the context of ecological risk assessment and natural resource damage assessment. From a regulatory standpoint, these approaches can help better assess the complex effects from coastal development activities to assessing ecosystem integrity pre- and post development or site remediation.

Effects of silver nitrate and silver nanoparticles on a planktonic community: general trends after short-term exposure

Among metal pollutants silver ions are one of the most toxic forms, and have thus been assigned to the highest toxicity class. Its toxicity to a wide range of microorganisms combined with its low toxicity to humans lead to the development of a wealth of silver-based products in many bactericidal applications accounting to more than 1000 nano-technology-based consumer products. Accordingly, silver is a widely distributed metal in the environment originating from its different forms of application as metal, salt and nanoparticle. A realistic assessment of silver nanoparticle toxicity in natural waters is, however, problematic and needs to be linked to experimental approaches. Here we apply metatranscriptome sequencing allowing for elucidating reactions of whole communities present in a water sample to stressors. We compared the toxicity of ionic silver and ligand-free silver nanoparticles by short term exposure on a natural community of aquatic microorganisms. We analyzed the effects of the treatments on metabolic pathways and species composition on the eukaryote metatranscriptome level in order to describe immediate molecular responses of organisms using a community approach. We found significant differences between the samples treated with 5 µg/L AgNO3 compared to the controls, but no significant differences in the samples treated with AgNP compared to the control samples. Statistical analysis yielded 126 genes (KO-IDs) with significant differential expression with a false discovery rate (FDR) <0.05 between the control (KO) and AgNO3 (NO3) groups. A KEGG pathway enrichment analysis showed significant results with a FDR below 0.05 for pathways related to photosynthesis. Our study therefore supports the view that ionic silver rather than silver nanoparticles are responsible for silver toxicity. Nevertheless, our results highlight the strength of metatranscriptome approaches for assessing metal toxicity on aquatic communities.

Dissolved and particulate copper exposure induces differing gene expression profiles and mechanisms of toxicity in the deposit feeding amphipod Melita plumulosa

Uptake of metals via ingestion is an important route of exposure for many invertebrates, and it has been suggested that the toxic response to metals accumulated via food differs from that of metals accumulated via the dissolved phase. To test this hypothesis, the deposit-feeding epibenthic amphipod Melita plumulosa was exposed to nontoxic or reproductively toxic concentrations of copper via the overlying water, via ingestion of sediment, or via a combination of the two. Rates of copper uptake from the two exposure routes were predicted using a biokinetic model. Gene expression profiles were measured via microarray analysis and confirmed via quantitative polymerase chain reaction. Differences in expression profiles were related to the exposure route more than to individual or combined rates of copper uptake. Chitinase and digestive protease transcript expression levels correlated to the copper uptake rate from sediment, rather than from the dissolved phase or combined total uptake rate. Overall, this study supports the hypothesis that metals accumulated via ingestion have a different mode of toxic action than metals taken up from water. Consequently, guidelines that only consider dissolved metal exposure, including equilibrium-partitioning-based guidelines, may underestimate the potential effects from deposited or resuspended metal-contaminated sediments.

Assessing mechanisms of toxicant response in the amphipod Melita plumulosa through transcriptomic profiling

This study describes the function of transcripts with altered abundance in the epibenthic amphipod, Melita plumulosa, following whole-sediment exposure to a series of common environmental contaminants. M. plumulosa were exposed for 48 h to sediments spiked and equilibrated with the following contaminants at concentrations predicted to cause sublethal effects to reproduction: porewater ammonia 30 mg L(-1); bifenthrin at 100 μg kg(-1); fipronil at 50 μg kg(-1); 0.6% diesel; 0.3% crude oil; 250 mg Cu kg(-1); 400 mg Ni kg(-1); and 400 mg Zn kg(-1). RNA was extracted and hybridized against a custom Agilent microarray developed for this species. Although the microarray represented a partial transcriptome and not all features on the array could be annotated, unique transcriptomic profiles were generated for each of the contaminant exposures. Hierarchical clustering grouped the expression profiles together by contaminant class, with copper and zinc, the petroleum products and nickel, and the pesticides each forming a distinct cluster. Many of the transcriptional changes observed were consistent with patterns previously described in other crustaceans. The changes in the transcriptome demonstrated that contaminant exposure caused changes in digestive function, growth and moulting, and the cytoskeleton following metal exposure, whereas exposure to petroleum products caused changes in carbohydrate metabolism, xenobiotic metabolism and hormone cycling. Functional analysis of these gene expression profiles can provide a better understanding of modes of toxic action and permits the prediction of mixture effects within contaminated ecosystems.

Selection of reliable biomarkers from PCR array analyses using relative distance computational model: methodology and proof-of-concept study

It is increasingly evident about the difficulty to monitor chemical exposure through biomarkers as almost all the biomarkers so far proposed are not specific for any individual chemical. In this proof-of-concept study, adult male zebrafish (Danio rerio) were exposed to 5 or 25 µg/L 17β-estradiol (E2), 100 µg/L lindane, 5 nM 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) or 15 mg/L arsenic for 96 h, and the expression profiles of 59 genes involved in 7 pathways plus 2 well characterized biomarker genes, vtg1 (vitellogenin1) and cyp1a1 (cytochrome P450 1A1), were examined. Relative distance (RD) computational model was developed to screen favorable genes and generate appropriate gene sets for the differentiation of chemicals/concentrations selected. Our results demonstrated that the known biomarker genes were not always good candidates for the differentiation of pair of chemicals/concentrations, and other genes had higher potentials in some cases. Furthermore, the differentiation of 5 chemicals/concentrations examined were attainable using expression data of various gene sets, and the best combination was the set consisting of 50 genes; however, as few as two genes (e.g. vtg1 and hspa5 [heat shock protein 5]) were sufficient to differentiate the five chemical/concentration groups in the present test. These observations suggest that multi-parameter arrays should be more reliable for biomonitoring of chemical exposure than traditional biomarkers, and the RD computational model provides an effective tool for the selection of parameters and generation of parameter sets.