DNA microarray technology in toxicogenomics of aquatic models: methods and applications

OMICs Approaches in Diarrhetic Shellfish Toxins Research

Diarrhetic shellfish toxins (DSTs) are among the most prevalent marine toxins in Europe's and in other temperate coastal regions. These toxins are produced by several dinoflagellate species; however, the contamination of the marine trophic chain is often attributed to species of the genus Dinophysis. This group of toxins, constituted by okadaic acid (OA) and analogous molecules (dinophysistoxins, DTXs), are highly harmful to humans, causing severe poisoning symptoms caused by the ingestion of contaminated seafood. Knowledge on the mode of action and toxicology of OA and the chemical characterization and accumulation of DSTs in seafood species (bivalves, gastropods and crustaceans) has significantly contributed to understand the impacts of these toxins in humans. Considerable information is however missing, particularly at the molecular and metabolic levels involving toxin uptake, distribution, compartmentalization and biotransformation and the interaction of DSTs with aquatic organisms. Recent contributions to the knowledge of DSTs arise from transcriptomics and proteomics research. Indeed, OMICs constitute a research field dedicated to the systematic analysis on the organisms' metabolisms. The methodologies used in OMICs are also highly effective to identify critical metabolic pathways affecting the physiology of the organisms. In this review, we analyze the main contributions provided so far by OMICs to DSTs research and discuss the prospects of OMICs with regard to the DSTs toxicology and the significance of these toxins to public health, food safety and aquaculture.

Gene Expression Profiling in Fish Toxicology: A Review

In this review, we present an overview of transcriptomic responses to chemical exposures in a variety of fish species. We have discussed the use of several molecular approaches such as northern blotting, differential display reverse transcription-polymerase chain reaction (DDRT-PCR), suppression subtractive hybridization (SSH), real time quantitative PCR (RT-qPCR), microarrays, and next-generation sequencing (NGS) for measuring gene expression. These techniques have been mainly used to measure the toxic effects of single compounds or simple mixtures in laboratory conditions. In addition, only few studies have been conducted to examine the biological significance of differentially expressed gene sets following chemical exposure. Therefore, future studies should focus more under field conditions using a multidisciplinary approach (genomics, proteomics and metabolomics) to understand the synergetic effects of multiple environmental stressors and to determine the functional significance of differentially expressed genes. Nevertheless, recent developments in NGS technologies and decreasing costs of sequencing holds the promise to uncover the complexity of anthropogenic impacts and biological effects in wild fish populations.

Cloning and expression analysis of a HSP70 gene from Korean rockfish (Sebastes schlegeli)

The gene encoding HSP70 was isolated from Korean rockfish Sebastes schlegeli by homologous cloning and rapid amplification of cDNA ends (RACE). The full-length of HSP70 cDNA was composed of 2259 bp and encoded a polypeptide of 639 amino acids. BLAST analysis showed that HSP70 of S. schlegeli shared high identities with those of the Lates calcarifer, Oreochromis niloticus, Seriola quinqueradiata HSP70s (88-89%). Our current study also revealed that HSP70 of Korean rockfish was expressed in many tissues by RT-PCR under unstressed condition. Quantitative real-time PCR showed that the expression patterns of Korean rockfish HSP70 were developmental stage-dependency. The expression of HSP70 was measured by quantitative real-time PCR after different oxygen treatments. The results showed that expression of HSP70 increased significantly after exposure to hypoxia for 30 min in gill and ovary, and then decreased for 60 min, and the level in spleen and liver gradually increased and reached the highest at 60 min. In addition, in gill, spleen and liver, the HSP70 mRNA level reached the maximum in hypoxia group after one hour different oxygen concentration stress. Increased amounts of serum thyroxine (T4), and triiodothyronine (T3) were also found during 30 min hypoxia treatment and 60 min normoxia group in our study. All of the results provide information to further study the mechanism of physiology and immune function under stress conditions of ovoviviparous teleosts.

Applications for next-generation sequencing in fish ecotoxicogenomics

The new technologies for next-generation sequencing (NGS) and global gene expression analyses that are widely used in molecular medicine are increasingly applied to the field of fish biology. This has facilitated new directions to address research areas that could not be previously considered due to the lack of molecular information for ecologically relevant species. Over the past decade, the cost of NGS has decreased significantly, making it possible to use non-model fish species to investigate emerging environmental issues. NGS technologies have permitted researchers to obtain large amounts of raw data in short periods of time. There have also been significant improvements in bioinformatics to assemble the sequences and annotate the genes, thus facilitating the management of these large datasets.The combination of DNA sequencing and bioinformatics has improved our abilities to design custom microarrays and study the genome and transcriptome of a wide variety of organisms. Despite the promising results obtained using these techniques in fish studies, NGS technologies are currently underused in ecotoxicogenomics and few studies have employed these methods. These issues should be addressed in order to exploit the full potential of NGS in ecotoxicological studies and expand our understanding of the biology of non-model organisms.

Identification of robust hypoxia biomarker candidates from fin of medaka (Oryzias latipes)

Aquatic hypoxia caused by organic pollution and eutrophication is a pressing worldwide water pollution problem. Better methods for monitoring oxygen levels are needed to assist efforts to maintain and protect the health of natural aquatic environments. In this project, we used a Japanese ricefish (medaka, Oryzias latipes) 8K oligonucleotide array as a platform to identify potential hypoxic biomarkers in different organs (fin, gill, liver and brain) upon exposure to hypoxia. The microarray results were validated by qRT-PCR employing a subset of candidate biomarkers. Interestingly, the largest number and most significant of hypoxia responding array features were detected in hypoxia exposed fin tissues. We identified 173 array features that exhibited a significant response (over 2 fold change in expression) upon exposure to hypoxic conditions and validated a subset of these by quantitative RT-PCR. These gene targets were subjected to annotation and gene ontology mining. Positively identifiable gene targets that may be useful for development of a rapid and accurate biomarker test using fin clips are discussed in relation to previous reports on hypoxia responsive genes.

Transcriptomics for understanding marine fish larval development1This review is part of a virtual symposium on current topics in aquaculture of marine fish and shellfish

The larval phase is a crucial period in the life of marine fish. During this phase, the organism will acquire the phenotype of an adult fish through the development of tissues and organs and the maturation of some of the principal physiological functions. Many biological processes (differentiation, cellular proliferation, growth, etc.) are regulated during this period. These regulations take place at different biological levels and particularly concern the expression of genes involved in larval ontogenesis processes. The development of bioinformatic resources (DNA or cDNA sequences) and molecular tools enabling high throughput gene expression analysis (microarrays) have allowed the transcriptome of marine fish species to be studied. In the present review, we summarize the main findings from transcriptomic investigations of development of marine fish larvae. Special attention is paid to investigations of transcriptomic patterns during postembryonic development and to the impact of environmental or nutritional factors on the transcriptome of marine fish larvae. Transcriptomic approaches will be especially useful in the future for investigating the effect of temperature and water acidification (or pH) on the development of different fish species in the context of global climate change.

Constructing a fish metabolic network model

We report the construction of a genome-wide fish metabolic network model, MetaFishNet, and its application to analyzing high throughput gene expression data. This model is a stepping stone to broader applications of fish systems biology, for example by guiding study design through comparison with human metabolism and the integration of multiple data types. MetaFishNet resources, including a pathway enrichment analysis tool, are accessible at http://metafishnet.appspot.com.

DNA chip analysis in diverse organisms with unsequenced genomes

Whether for basic research or biotechnology, DNA microarrays have become indispensable tools for studying the transcriptome. Normally, analyses begin with a set of known cDNA sequences to prepare microarray chips specific for a target organism with an extensively sequenced and annotated genome. For many organisms, however, genome programs are not complete or have not been initiated. The present study demonstrates that, whether using homologous or heterologous arrays, the chances of seeing interesting differences are similar. When a specific DNA microarray is not available, the results indicate that a reverse approach based on a heterologous array can be used to probe for interesting differences in gene expression. This may be sufficient in many studies but, if necessary, the genes exhibiting the most significant changes subsequently could be identified by traditional molecular approaches. Such a reverse strategy can provide a convenient and inexpensive approach to probe for significant genetic changes in many diverse studies, to monitor or mine critical biological information for basic or applied research, long before complete sequence data are available.

Multi-gene expression-based statistical approaches to predicting patients' clinical outcomes and responses

Gene expression profiling technique now enables scientists to obtain a genome-wide picture of cellular functions on various human disease mechanisms which has also proven to be extremely valuable in forecasting patients' prognosis and therapeutic responses. A wide range of multivariate techniques have been employed in biomedical applications on such expression profiling data in order to identify expression biomarkers that are highly associated with patients' clinical outcome and to train multi-gene prediction models that can forecast various human disease outcome and drug toxicities. We provide here a brief overview on some of these approaches, succinctly summarizing relevant basic concepts, statistical algorithms, and several practical applications. We also introduce our recent in vitro molecular expression-based algorithm, the so-called COXEN technique, which uses specialized gene profile signatures as a Rosetta Stone for translating the information between two different biological systems or populations.

Comparison of gene expression responses to hypoxia in viviparous (Xiphophorus) and oviparous (Oryzias) fishes using a medaka microarray

Gene expression profiling using DNA microarray technology is a useful tool for assessing gene transcript level responses after an organism is exposed to environmental stress. Herein, we detail results from studies using an 8 k medaka (Oryzias latipes) microarray to assess modulated gene expression patterns upon hypoxia exposure of the live-bearing aquaria fish, Xiphophorus maculatus. To assess the reproducibility and reliability of using the medaka array in cross-genus hybridization, a two-factor ANOVA analysis of gene expression was employed. The data show the tissue source of the RNA used for array hybridization contributed more to the observed response of modulated gene targets than did the species source of the RNA. In addition, hierarchical clustering via heat map analyses of groupings of tissues and species (Xiphophorus and medaka) suggests that hypoxia induced similar responses in the same tissues from these two diverse aquatic model organisms. Our Xiphophorus results indicate 206 brain, 37 liver, and 925 gill gene targets exhibit hypoxia induced expression changes. Analysis of the Xiphophorus data to determine those features exhibiting a significant (p<0.05)+/-3 fold change produced only two gene targets within brain tissue and 80 features within gill tissue. Of these 82 characterized features, 39 were identified via homology searching (cut-off E-value of 1 x 10(-5)) and placed into one or more biological process gene ontology groups. Among these 39 genes, metabolic energy changes and manipulation was the most affected biological pathway (13 genes).

Chronic Mycobacterium marinum infection acts as a tumor promoter in Japanese Medaka (Oryzias latipes)

An accumulating body of research indicates there is an increased cancer risk associated with chronic infections. The genus Mycobacterium contains a number of species, including M. tuberculosis, which mount chronic infections and have been implicated in higher cancer risk. Several non-tuberculosis mycobacterial species, including M. marinum, are known to cause chronic infections in fish and like human tuberculosis, often go undetected. The elevated carcinogenic potential for fish colonies infected with Mycobacterium spp. could have far reaching implications because fish models are widely used to study human diseases. Japanese medaka (Oryzias latipes) is an established laboratory fish model for toxicology, mutagenesis, and carcinogenesis; and produces a chronic tuberculosis-like disease when infected by M. marinum. We examined the role that chronic mycobacterial infections play in cancer risk for medaka. Experimental M. marinum infections of medaka alone did not increase the mutational loads or proliferative lesion incidence in all tissues examined. However, we showed that chronic M. marinum infections increased hepatocellular proliferative lesions in fish also exposed to low doses of the mutagen benzo[a]pyrene. These results indicate that chronic mycobacterial infections of medaka are acting as tumor promoters and thereby suggest increased human risks for cancer promotion in human populations burdened with chronic tuberculosis infections.

Sensitive and robust gene expression changes in fish exposed to estrogen--a microarray approach

Background

Vitellogenin is a well established biomarker for estrogenic exposure in fish. However, effects on gonadal differentiation at concentrations of estrogen not sufficient to give rise to a measurable vitellogenin response suggest that more sensitive biomarkers would be useful. Induction of zona pellucida genes may be more sensitive but their specificities are not as clear. The objective of this study was to find additional sensitive and robust candidate biomarkers of estrogenic exposure.

Results

Hepatic mRNA expression profiles were characterized in juvenile rainbow trout exposed to a measured concentration of 0.87 and 10 ng ethinylestradiol/L using a salmonid cDNA microarray. The higher concentration was used to guide the subsequent identification of generally more subtle responses at the low concentration not sufficient to induce vitellogenin. A meta-analysis was performed with data from the present study and three similar microarray studies using different fish species and platforms. Within the generated list of presumably robust responses, several well-known estrogen-regulated genes were identified. Two genes, confirmed by quantitative RT-PCR (qPCR), fulfilled both the criteria of high sensitivity and robustness; the induction of the genes encoding zona pellucida protein 3 and a nucleoside diphosphate kinase (nm23).

Conclusion

The cross-species, cross-platform meta-analysis correctly identified several robust responses. This adds confidence to our approach used for identifying candidate biomarkers. Specifically, we propose that analyses of an nm23 gene together with zona pellucida genes may increase the possibilities to detect an exposure to low levels of estrogenic compounds in fish.