Gene expression analysis of estrogenic compounds in the liver of common carp (Cyprinus carpio) using a custom cDNA microarray

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01787-9.

17α-ethinylestradiol and 4-tert-octylphenol concurrently disrupt the immune response of common carp

The aquatic environment is massively polluted with endocrine-disrupting compounds (EDCs) including synthetic estrogens (e.g. 17α-ethinylestradiol, EE2) and alkylphenols (e.g. 4-tert-octylphenol, 4t-OP). A major mechanism of action for estrogenic EDCs is their interaction with estrogen receptors and consequently their modulation of the action of enzymes involved in steroid conversion e.g. aromatase CYP19. We now studied the effects of EE2 and 4t-OP on the anti-bacterial immune response of common carp. We investigated effects on the number/composition of inflammatory leukocytes and on the gene expression of mediators that regulate inflammation and EDC binding. In vitro we found that high concentrations of both EE2 and 4t-OP down-regulated IFN-γ2 and IFN-γ-dependent immune responses in LPS-stimulated monocytes/macrophages. Similarly, during bacterial infection in fish, in vivo treated with EE2 and 4t-OP, decreased gene expression of il-12p35 and of ifn-γ2 was found in the focus of inflammation. Moreover, during A. salmonicida-induced infection in EE2-treated carp, but not in fish fed with 4t-OP-treated food, we found an enhanced inflammatory reaction manifested by high number of inflammatory peritoneal leukocytes, including phagocytes and higher expression of pro-inflammatory mediators (inos, il-1β, cxcl8_l2). Furthermore, in the liver, EE2 down-regulated the expression of acute phase proteins: CRPs and C3. Importantly, both in vitro and in vivo, EDCs altered the expression of estrogen receptors: nuclear (erα and erβ) and membrane (gpr30). EDCs also induced up-regulation of the cyp19b gene. Our findings reveal that contamination of the aquatic milieu with estrogenic EDCs, may considerably violate the subtle and particular allostatic interactions between the immune response and endogenous estrogens and this may have negative consequences for fish health.

Twenty years of transcriptomics, 17alpha-ethinylestradiol, and fish

In aquatic toxicology, perhaps no pharmaceutical has been investigated more intensely than 17alpha-ethinylestradiol (EE2), the active ingredient of the birth control pill. At the turn of the century, the fields of comparative endocrinology and endocrine disruption research witnessed the emergence of omics technologies, which were rapidly adapted to characterize potential hazards associated with exposures to environmental estrogens, such as EE2. Since then, significant advances have been made by the scientific community, and as a result, much has been learned about estrogen receptor signaling in fish from environmental xenoestrogens. Vitellogenin, the egg yolk precursor protein, was identified as a major estrogen-responsive gene, establishing itself as the premier biomarker for estrogenic exposures. Omics studies have identified a plethora of estrogen responsive genes, contributing to a wealth of knowledge on estrogen-mediated regulatory networks in teleosts. There have been ~40 studies that report on transcriptome responses to EE2 in a variety of fish species (e.g., zebrafish, fathead minnows, rainbow trout, pipefish, mummichog, stickleback, cod, and others). Data on the liver and testis transcriptomes dominate in the literature and have been the subject of many EE2 studies, yet there remain knowledge gaps for other tissues, such as the spleen, kidney, and pituitary. Inter-laboratory genomics studies have revealed transcriptional networks altered by EE2 treatment in the liver; networks related to amino acid activation and protein folding are increased by EE2 while those related to xenobiotic metabolism, immune system, circulation, and triglyceride storage are suppressed. EE2-responsive networks in other tissues are not as comprehensively defined which is a knowledge gap as regulated networks are expected to be tissue-specific. On the horizon, omics studies for estrogen-mediated effects in fish include: (1) Establishing conceptual frameworks for incorporating estrogen-responsive networks into environmental monitoring programs; (2) Leveraging in vitro and computational toxicology approaches to identify chemicals associated with estrogen receptor-mediated effects in fish (e.g., male vitellogenin production); (3) Discovering new tissue-specific estrogen receptor signaling pathways in fish; and (4) Developing quantitative adverse outcome pathway predictive models for estrogen signaling. As we look ahead, research into EE2 over the past several decades can serve as a template for the array of hormones and endocrine active substances yet to be fully characterized or discovered.

Metabolic signatures of bisphenol A and genistein in Atlantic salmon liver cells

Screening has revealed that aquafeeds with high inclusion of plant material may contain small amounts of endocrine disrupting agricultural pesticides. In this work, bisphenol A (BPA) and genistein (GEN) were selected as model endocrine disrupting toxicants with impact on DNA methylation in fish. Atlantic salmon hepatocytes were exposed in vitro to four concentrations of BPA and GEN (0.1, 1.0, 10 and 100 μM) for 48 h. Toxicity endpoints included cytotoxicity, global DNA methylation, targeted transcriptomics and metabolomic screening (100 μM). GEN was not cytotoxic in concentrations up to 100 μM, whereas one out of two cell viability assays indicated a cytotoxic response to 100 μM BPA. Compared to the control, significant global DNA hypomethylation was observed at 1.0 μM BPA. Both compounds upregulated cyp1a1 transcription at 100 μM, while estrogenic markers esr1 and vtg1 responded strongest at 10 μM. Dnmt3aa transcription was downregulated by both compounds at 100 μM. Metabolomic screening showed that BPA and GEN resulted in significant changes in numerous biochemical pathways consistent with alterations in carbohydrate metabolism, indicating perturbation in glucose homeostasis and energy generation, and glutamate metabolism. Pathway analysis showed that while the superpathway of methionine degradation was among the most strongly affected pathways by BPA, GEN induced changes to uridine and pyrimidine biosynthesis. In conclusion, this mechanistic study proposes metabolites associated with glucose and glutamate metabolism, glucuronidation detoxification, as well as cyp1a1, vtg1, esr1, ar, dnmt3aa, cdkn1b and insig1 as transcriptional markers for BPA and GEN exposure in fish liver cells.

Estrogen-responsive gene networks in the teleost liver: What are the key molecular indicators?

An overarching goal of environmental genomics is to leverage sensitive suites of markers that are robust and reliable to assess biological responses in a range of species inhabiting variable environments. The objective of this study was to identify core groups of transcripts and molecular signaling pathways that respond to 17alpha-ethylinestadiol (EE2), a ubiquitous estrogenic contaminant, using transcriptome datasets generated from six independent laboratories. We sought to determine which biomarkers and gene networks were those most robust and reliably detected in multiple laboratories. Six laboratories conducted microarray analysis in pieces of the same liver from male fathead minnows exposed to ∼15ng/L EE2 for 96h. There were common transcriptional networks identified in every dataset. These included down-regulation of gene networks associated with blood clotting, complement activation, triglyceride storage, and xenobiotic metabolism. Noteworthy was that more than ∼85% of the gene networks were suppressed by EE2. Leveraging both these data and those mined from the Comparative Toxicogenomics Database (CTD), we narrowed in on an EE2-responsive transcriptional network. All transcripts in this network responded ∼±5-fold or more to EE2, increasing reliability of detection. This network included estrogen receptor alpha, transferrin, myeloid cell leukemia 1, insulin like growth factor 1, insulin like growth factor binding protein 2, and methionine adenosyltransferase 2A. This estrogen-responsive interactome has the advantage over single markers (e.g. vitellogenin) in that these entities are directly connected to each other based upon evidence of expression regulation and protein binding. Thus, it represents an interacting functional suite of estrogenic markers. Vitellogenin, the gold standard for estrogenic exposures, can show high individual variability in its response to estrogens, and the use of a multi-gene approach for estrogenic chemicals is expected to improve sensitivity. In our case, the coefficient of variation was significantly lowered by the gene network (∼67%) compared to Vtg alone, supporting the use of this transcriptional network as a sensitive alternative for detecting estrogenic effluents and chemicals. We propose that screening chemicals for estrogenicity using interacting genes within a defined expression network will improve sensitivity, accuracy, and reduce the number of animals required for endocrine disruption assessments.

Zearalenone (ZEN) and Its Influence on Regulation of Gene Expression in Carp (Cyprinus carpio L.) Liver Tissue

Zearalenone (ZEN) is a frequently-occurring mycotoxin in both animal and fish feeds. In order to characterize its effects on carp, three groups of fish were fed for 28 days with feeds contaminated with three different levels of ZEN (low: 332 µg kg⁻¹, medium: 621 µg kg⁻¹, and high: 797 µg kg⁻¹ feed). The reversibility of the effects of ZEN was assessed by feeding all of the groups with uncontaminated feed for a further 14 days. Gene expression of immune genes in the liver tissue of the fish was analysed, revealing reduced expressions of immune, antioxidative, and estrogen-related genes after the fish had been exposed to ZEN. However, the expression of vacuole-type H⁺ ATPase increased substantially with ZEN exposure, thus supporting the previously-reported sensitivity of lysosomal functions to ZEN. Feeding the fish with a ZEN-free diet for a further two weeks changed the effects of ZEN on the expression of some genes, including the expressions of the cytokines IL-1β, IL-8, IL-10, and arginase 2, which were not influenced after four weeks of treatment, but showed lower values after the recovery phase in fish previously treated with ZEN compared with the control group. In summary, this study confirmed the broad effects of ZEN on different essential functions in carp and suggests that the current maximum allowable levels in compound feed are too high to prevent damage to fish.

Estrogen-dependent seasonal adaptations in the immune response of fish

Clinical and experimental evidence shows that estrogens affect immunity in mammals. Less is known about this interaction in the evolutionary older, non-mammalian, vertebrates. Fish form an excellent model to identify evolutionary conserved neuroendocrine-immune interactions: i) they are the earliest vertebrates with fully developed innate and adaptive immunity, ii) immune and endocrine parameters vary with season, and iii) physiology is constantly disrupted by increasing contamination of the aquatic environment. Neuro-immuno-endocrine interactions enable adaption to changing internal and external environment and are based on shared signaling molecules and receptors. The presence of specific estrogen receptors on/in fish leukocytes, implies direct estrogen-mediated immunoregulation. Fish leukocytes most probably are also capable to produce estrogens as they express the cyp19a and cyp19b - genes, encoding aromatase cytochrome P450, the enzyme critical for conversion of C19 steroids to estrogens. Immunoregulatory actions of estrogens, vary among animal species, and also with dose, target cell type, or physiological condition (e.g., infected/non-infected, reproductive status). They moreover are multifaceted. Interestingly, season-dependent changes in immune status correlate with changes in the levels of circulating sex hormones. Whereas E2 circulating in the bloodstream is perhaps the most likely candidate to be the physiological mediator of systemic immune-reproductive trade-offs, leukocyte-derived hormones are hypothesized to be mainly involved in local tuning of the immune response. Contamination of the aquatic environment with estrogenic EDCs may violate the delicate and precise allostatic interactions between the endogenous estrogen system and the immune system. This has negative effects on fish health, but will also affect the physiology of its consumers.

Genome Wide Identification, Phylogeny, and Expression of Aquaporin Genes in Common Carp (Cyprinus carpio)

Background

Aquaporins (Aqps) are integral membrane proteins that facilitate the transport of water and small solutes across cell membranes. Among vertebrate species, Aqps are highly conserved in both gene structure and amino acid sequence. These proteins are vital for maintaining water homeostasis in living organisms, especially for aquatic animals such as teleost fish. Studies on teleost Aqps are mainly limited to several model species with diploid genomes. Common carp, which has a tetraploidized genome, is one of the most common aquaculture species being adapted to a wide range of aquatic environments. The complete common carp genome has recently been released, providing us the possibility for gene evolution of aqp gene family after whole genome duplication.

Results

In this study, we identified a total of 37 aqp genes from common carp genome. Phylogenetic analysis revealed that most of aqps are highly conserved. Comparative analysis was performed across five typical vertebrate genomes. We found that almost all of the aqp genes in common carp were duplicated in the evolution of the gene family. We postulated that the expansion of the aqp gene family in common carp was the result of an additional whole genome duplication event and that the aqp gene family in other teleosts has been lost in their evolution history with the reason that the functions of genes are redundant and conservation. Expression patterns were assessed in various tissues, including brain, heart, spleen, liver, intestine, gill, muscle, and skin, which demonstrated the comprehensive expression profiles of aqp genes in the tetraploidized genome. Significant gene expression divergences have been observed, revealing substantial expression divergences or functional divergences in those duplicated aqp genes post the latest WGD event.

Conclusions

To some extent, the gene families are also considered as a unique source for evolutionary studies. Moreover, the whole set of common carp aqp gene family provides an essential genomic resource for future biochemical, toxicological, physiological, and evolutionary studies in common carp.

Phylogeny and Evolution of Multiple Common Carp (Cyprinus carpio L.) Populations Clarified by Phylogenetic Analysis Based on Complete Mitochondrial Genomes

Common carp (Cyprinus carpio L.) is one of the oldest, most widely farmed commercially important freshwater fish in the world. However, many undetermined phylogenetic relationships and origins of common carp lineages remain, which are obstacles to conservation and genetic breeding of this species. Phylogenetic analyses based on molecular tools are helpful to distinguish the origin of species, understand and clarify their evolutionary history, and provide a genetic basis for selective breeding. In this study, we demonstrated a method to extract complete mitochondrial genome sequences from whole-genome resequencing data using the Illumina platform. The complete mitochondrial genome sequences of 26 individuals representing nine strains were obtained and subjected to a phylogenetic analysis. We reconstructed the phylogenetic topologies of the nine strains and analyzed the haplotypes. Results from both analyses suggested that the genome sequences belonged to two distinct subspecies from Europe and East Asia. We also estimated the time of divergence of the nine strains, which was up to 100 KYA. The phylogenetic results clarified the breeding history of Songpu mirror carp and suggest that this species may be hybrid of paternal European mirror carp and maternal Xingguo red carp. The results also support a previous hypothesis that koi may have originated from or have close ancestry with Oujiang color carp in China.

Dietary Exposure to Individual Polybrominated Diphenyl Ether Congeners BDE-47 and BDE-99 Alters Innate Immunity and Disease Susceptibility in Juvenile Chinook Salmon

Polybrominated diphenyl ethers (PBDEs), used as commercial flame-retardants, are bioaccumulating in threatened Pacific salmon. However, little is known of PBDE effects on critical physiological functions required for optimal health and survival. BDE-47 and BDE-99 are the predominant PBDE congeners found in Chinook salmon collected from the Pacific Northwest. In the present study, both innate immunity (phagocytosis and production of superoxide anion) and pathogen challenge were used to evaluate health and survival in groups of juvenile Chinook salmon exposed orally to either BDE-47 or BDE-99 at environmentally relevant concentrations. Head kidney macrophages from Chinook salmon exposed to BDE-99, but not those exposed to BDE-47, were found to have a reduced ability in vitro to engulf foreign particles. However, both congeners increased the in vitro production of superoxide anion in head kidney macrophages. Salmon exposed to either congener had reduced survival during challenge with the pathogenic marine bacteria Listonella anguillarum. The concentration response curves generated for these end points were nonmonotonic and demonstrated a requirement for using multiple environmentally relevant PBDE concentrations for effect studies. Consequently, predicting risk from toxicity reference values traditionally generated with monotonic concentration responses may underestimate PBDE effect on critical physiological functions required for optimal health and survival in salmon.

Influence of short-term exposure to low levels of 17α-ethynylestradiol on expression of genes involved in immunity and on immune parameters in rainbow trout, Oncorhynchus mykiss

Fish are exposed to endocrine-disrupting chemicals (EDC), which are well known to disturb not only the reproductive system but also the immune system in vertebrates. However, the mechanisms by which these compounds are able to modify fish immunity are not well understood. In order to test the EE2 effects on immunity in selected organs, we exposed rainbow trout male juveniles for 3 weeks to EE2 concentrations ranging from 0.01 to 1 μg/L. The results of this study suggest that EE2 affects the immunity of rainbow trout in a tissue dependent manner. This molecule affects both cellular and humoral immune systems. Indeed, blood leukocyte populations, as well as hepatic and plasma lysozyme, plasma MPO and renal complement activities, are modulated by EE2. Moreover, EE2 alters the gene expression of some mucus compounds, hepatic expression of complement sub-unit and lysozyme, or genes involved in the hepatic phagocytosis and transport of immunoglobulin across the liver.

Dietary exposure of 17-alpha ethinylestradiol modulates physiological endpoints and gene signaling pathways in female largemouth bass (Micropterus salmoides)

17Alpha-ethinylestradiol (EE2), used for birth control in humans, is a potent estrogen that is found in wastewater at low concentrations (ng/l). EE2 has the ability to interfere with the endocrine system of fish, affecting reproduction which can result in population level effects. The objective of this study was to determine if dietary exposure to EE2 would alter gene expression patterns and key pathways in the liver and ovary and whether these could be associated with reproductive endpoints in female largemouth bass during egg development. Female LMB received 70ng EE2/g feed (administered at 1% of body weight) for 60 days. EE2 dietary exposure significantly reduced plasma vitellogenin concentrations by 70%. Hepatosomatic and gonadosomatic indices were also decreased with EE2 feeding by 38.5% and 40%, respectively. Transcriptomic profiling revealed that there were more changes in steady state mRNA levels in the liver compared to the ovary. Genes associated with reproduction were differentially expressed, such as vitellogenin in the liver and aromatase in the gonad. In addition, a set of genes related with oxidative stress (e.g. glutathione reductase and glutathione peroxidase) were identified as altered in the liver and genes associated with the immune system (e.g. complement component 1, and macrophage-inducible C-type lectin) were altered in the gonad. In a follow-up study with 0.2ng EE2/g feed for 60 days, similar phenotypic and gene expression changes were observed that support these findings with the higher concentrations. This study provides new insights into how dietary exposure to EE2 interferes with endocrine signaling pathways in female LMB during a critical period of reproductive oogenesis.

Development and evaluation of the first high-throughput SNP array for common carp (Cyprinus carpio)

Background

A large number of single nucleotide polymorphisms (SNPs) have been identified in common carp (Cyprinus carpio) but, as yet, no high-throughput genotyping platform is available for this species. C. carpio is an important aquaculture species that accounts for nearly 14% of freshwater aquaculture production worldwide. We have developed an array for C. carpio with 250,000 SNPs and evaluated its performance using samples from various strains of C. carpio.

Results

The SNPs used on the array were selected from two resources: the transcribed sequences from RNA-seq data of four strains of C. carpio, and the genome re-sequencing data of five strains of C. carpio. The 250,000 SNPs on the resulting array are distributed evenly across the reference C.carpio genome with an average spacing of 6.6 kb. To evaluate the SNP array, 1,072 C. carpio samples were collected and tested. Of the 250,000 SNPs on the array, 185,150 (74.06%) were found to be polymorphic sites. Genotyping accuracy was checked using genotyping data from a group of full-siblings and their parents, and over 99.8% of the qualified SNPs were found to be reliable. Analysis of the linkage disequilibrium on all samples and on three domestic C.carpio strains revealed that the latter had the longer haplotype blocks. We also evaluated our SNP array on 80 samples from eight species related to C. carpio, with from 53,526 to 71,984 polymorphic SNPs. An identity by state analysis divided all the samples into three clusters; most of the C. carpio strains formed the largest cluster.

Conclusions

The Carp SNP array described here is the first high-throughput genotyping platform for C. carpio. Our evaluation of this array indicates that it will be valuable for farmed carp and for genetic and population biology studies in C. carpio and related species.

Estrogen-related receptor γ is an in vivo receptor of bisphenol A

Bisphenol A (BPA) is an endocrine disruptor that displays estrogenic activity. Several reports suggest that BPA may have estrogen receptor-independent effects. In zebrafish, 50 μM BPA exposure induces otic vesicle abnormalities, including otolith aggregation. The purpose of this study was to test if BPA action was mediated in vivo during zebrafish development by the orphan nuclear estrogen related receptor (ERR) γ. Combining pharmacological and functional approaches, we demonstrate that the zebrafish ERRγ mediates BPA-induced malformations in otoliths. Using different bisphenol derivatives, we show that different compounds can induce a similar otolith phenotype than BPA and that the binding affinity of these derivatives to the zebrafish ERRγ correlates with their ability to induce otolith malformations. Morpholino knockdown of ERRγ function suppresses the BPA effect on otoliths whereas overexpression of ERRγ led to a BPA-like otolith phenotype. Moreover, a subphenotypical dose of BPA (1 μM) combined with ERRγ overexpression led to a full-dose (50 μM) BPA otolith phenotype. We therefore conclude that ERRγ mediates the otic vesicle phenotype generated by BPA. Our results suggest that the range of pathways perturbed by this compound and its potential harmful effect are larger than expected.-Tohmé, M., Prud'homme, S. M., Boulahtouf, A., Samarut, E., Brunet, F., Bernard, L., Bourguet, W., Gibert, Y., Balaguer, P., Laudet, V. Estrogen-related receptor γ is an in vivo receptor of bisphenol A.

A dense genetic linkage map for common carp and its integration with a BAC-based physical map

Background

Common carp (Cyprinus carpio) is one of the most important aquaculture species with an annual global production of 3.4 million metric tons. It is also an important ornamental species as well as an important model species for aquaculture research. To improve the economically important traits of this fish, a number of genomic resources and genetic tools have been developed, including several genetic maps and a bacterial artificial chromosome (BAC)-based physical map. However, integrated genetic and physical maps are not available to study quantitative trait loci (QTL) and assist with fine mapping, positional cloning and whole genome sequencing and assembly. The objective of this study was to integrate the currently available BAC-based physical and genetic maps.

Results

The genetic map was updated with 592 novel markers, including 312 BAC-anchored microsatellites and 130 SNP markers, and contained 1,209 genetic markers on 50 linkage groups, spanning 3,565.9 cM in the common carp genome. An integrated genetic and physical map of the common carp genome was then constructed, which was composed of 463 physical map contigs and 88 single BACs. Combined lengths of the contigs and single BACs covered a physical length of 498.75 Mb, or around 30% of the common carp genome. Comparative analysis between common carp and zebrafish genomes was performed based on the integrated map, providing more insights into the common carp specific whole genome duplication and segmental rearrangements in the genome.

Conclusion

We integrated a BAC-based physical map to a genetic linkage map of common carp by anchoring BAC-associated genetic markers. The density of the genetic linkage map was significantly increased. The integrated map provides a tool for both genetic and genomic studies of common carp, which will help us to understand the genomic architecture of common carp and facilitate fine mapping and positional cloning of economically important traits for genetic improvement and modification.

Transcriptomic analyses of sexual dimorphism of the zebrafish liver and the effect of sex hormones

The liver is one of the most sex-dimorphic organs in both oviparous and viviparous animals. In order to understand the molecular basis of the difference between male and female livers, high-throughput RNA-SAGE (serial analysis of gene expression) sequencing was performed for zebrafish livers of both sexes and their transcriptomes were compared. Both sexes had abundantly expressed genes involved in translation, coagulation and lipid metabolism, consistent with the general function of the liver. For sex-biased transcripts, from in addition to the high enrichment of vitellogenin transcripts in spawning female livers, which constituted nearly 80% of total mRNA, it is apparent that the female-biased genes were mostly involved in ribosome/translation, estrogen pathway, lipid transport, etc, while the male-biased genes were enriched for oxidation reduction, carbohydrate metabolism, coagulation, protein transport and localization, etc. Sexual dimorphism on xenobiotic metabolism and anti-oxidation was also noted and it is likely that retinol x receptor (RXR) and liver x receptor (LXR) play central roles in regulating the sexual differences of lipid and cholesterol metabolisms. Consistent with high ribosomal/translational activities in the female liver, female-biased genes were significantly regulated by two important transcription factors, Myc and Mycn. In contrast, Male livers showed activation of transcription factors Ppargc1b, Hnf4a, and Stat4, which regulate lipid and glucose metabolisms and various cellular activities. The transcriptomic responses to sex hormones, 17β-estradiol (E2) or 11-keto testosterone (KT11), were also investigated in both male and female livers and we found that female livers were relatively insensitive to sex hormone disturbance, while the male livers were readily affected. E2 feminized male liver by up-regulating female-biased transcripts and down-regulating male-biased transcripts. The information obtained in this study provides comprehensive insights into the sexual dimorphism of zebrafish liver transcriptome and will facilitate further development of the zebrafish as a human liver disease model.

Genome-Wide SNP Discovery from Transcriptome of Four Common Carp Strains

Background

Single nucleotide polymorphisms (SNPs) have been used as genetic marker for genome-wide association studies in many species. Gene-associated SNPs could offer sufficient coverage in trait related research and further more could themselves be causative SNPs for traits. Common carp (Cyprinus carpio) is one of the most important aquaculture species in the world accounting for nearly 14% of freshwater aquaculture production. There are various strains of common carp with different economic traits, however, the genetic mechanism underlying the different traits have not been elucidated yet. In this project, we identified a large number of gene-associated SNPs from four strains of common carp using next-generation sequencing.

Results

Transcriptome sequencing of four strains of common carp (mirror carp, purse red carp, Xingguo red carp, Yellow River carp) was performed with Solexa HiSeq2000 platform. De novo assembled transcriptome was used as reference for alignments, and SNP calling was done through BWA and SAMtools. A total of 712,042 Intra-strain SNPs were discovered in four strains, of which 483,276 SNPs for mirror carp, 486,629 SNPs for purse red carp, 478,028 SNPs for Xingguo red carp and 488,281 SNPs for Yellow River carp were discovered, respectively. Besides, 53,893 inter-SNPs were identified. Strain-specific SNPs of four strains were 53,938, 53,866, 48,701, 40,131 in mirror carp, purse red carp, Xingguo red carp and Yellow River carp, respectively. GO and KEGG pathway analysis were done to reveal strain-specific genes affected by strain-specific non-synonymous SNPs. Validation of selected SNPs revealed that 48% percent of SNPs (12 of 25) were tested to be true SNPs.

Conclusions

Transcriptome analysis of common carp using RNA-Seq is a cost-effective way of generating numerous reads for SNP discovery. After validation of identified SNPs, these data will provide a solid base for SNP array designing and genome-wide association studies.

Exploring androgen-regulated pathways in teleost fish using transcriptomics and proteomics

In the environment, there are aquatic pollutants that disrupt androgen signaling in fish. Laboratory and field-based experiments have utilized omics technologies to characterize the molecular mechanisms underlying androgen-receptor agonism/antagonism. Transcriptomics and proteomics studies with 17β-trenbolone, a growth-promoting pharmaceutical found in water systems surrounding cattle feed lots, and androgens such as 17α-methyltestosterone and 17α-methyldihydrotestosterone, have been conducted in ovary and liver of fish that include the fathead minnow (FHM) (Pimephales promelas), common carp (Cyprinus carpio), Qurt medaka (Oryzias latipes), and zebrafish (Danio rerio). In this mini-review, we survey recent omics studies in fish and reveal that, despite the diversity of species and tissues examined, there are common cellular responses that are observed with waterborne androgenic treatments. Recurring themes in gene ontology include apoptosis, transport and oxidation of lipids, synthesis and transport of hormones, immune response, protein metabolism, and cell proliferation. However, we also discuss other mechanisms other than androgen receptor (AR) activation, such as responses to toxicant stress, estrogen receptor agonism, aromatization of androgens into estrogens, and inhibitory feedback mechanisms by high levels of androgens that may also explain molecular responses in fish. To further explore androgen-responsive protein networks, a sub-network enrichment analysis was performed on protein data collected from the livers of female FHMs exposed to 17β-trenbolone. We construct a putative AR-regulated protein/cell process network in the liver that includes B-lymphocyte differentiation, xenobiotic clearance, low-density lipoprotein oxidation, proliferation of smooth muscle cells, and permeability of blood vessels. We demonstrate that construction of protein networks can offer insight into cell processes that are potentially regulated by androgens.

Characterization of common carp transcriptome: sequencing, de novo assembly, annotation and comparative genomics

BACKGROUND

Common carp (Cyprinus carpio) is one of the most important aquaculture species of Cyprinidae with an annual global production of 3.4 million tons, accounting for nearly 14% of the freshwater aquaculture production in the world. Due to the economical and ecological importance of common carp, genomic data are eagerly needed for genetic improvement purpose. However, there is still no sufficient transcriptome data available. The objective of the project is to sequence transcriptome deeply and provide well-assembled transcriptome sequences to common carp research community.

RESULT

Transcriptome sequencing of common carp was performed using Roche 454 platform. A total of 1,418,591 clean ESTs were collected and assembled into 36,811 cDNA contigs, with average length of 888 bp and N50 length of 1,002 bp. Annotation was performed and a total of 19,165 unique proteins were identified from assembled contigs. Gene ontology and KEGG analysis were performed and classified all contigs into functional categories for understanding gene functions and regulation pathways. Open Reading Frames (ORFs) were detected from 29,869 (81.1%) contigs with an average ORF length of 763 bp. From these contigs, 9,625 full-length cDNAs were identified with sequence length from 201 bp to 9,956 bp. Comparative analysis revealed that 27,693(75.2%) contigs have significant similarity to zebrafish Refseq proteins, and 24,371(66.2%), 24,501(66.5%) and 25,025(70.0%) to teraodon, medaka and three-spined stickleback refseq proteins. A total of 2,064 microsatellites were initially identified from 1,730 contigs, and 1,639 unique sequences had sufficient flanking sequences on both sides for primer design.

CONCLUSION

The transcriptome of common carp had been deep sequenced, de novo assembled and characterized, providing the valuable resource for better understanding of common carp genome. The transcriptome data will facilitate future functional studies on common carp genome, and gradually apply in breeding programs of common carp, as well as closely related other Cyprinids.

Lumiestrone is Photochemically Derived from Estrone and may be Released to the Environment without Detection

Endocrine disrupting chemicals are adversely affecting the reproductive health and metabolic status of aquatic vertebrates. Estrone is often the dominant natural estrogen in urban sewage, yet little is known about its environmental fate and biological effects. Increased use of UV-B radiation for effluent treatments, and exposure of effluents to sunlight in holding ponds led us to examine the effects of environmentally relevant levels of UV-B radiation on the photodegradation potential of estrone. Surprisingly, UV-B-mediated degradation leads to the photoproduction of lumiestrone, a little known 13α-epimer form of estrone. We show for the first time that lumiestrone possesses novel biological activity. In vivo treatment with estrone stimulated estrogen receptor (ER) α mRNA production in the male goldfish liver, whereas lumiestrone was without effect, suggesting a total loss of estrogenicity. In contrast, results from in vitro ER-dependent reporter gene assays indicate that lumiestrone showed relatively higher estrogenic potency with the zebrafish ERβ2 than zfERα, suggesting that it may act through an ERβ-selectivity. Lumiestrone also activated human ERs. Microarray analysis of male goldfish liver following in vivo treatments showed that lumiestrone respectively up- and down-regulated 20 and 69 mRNAs, which was indicative of metabolic upsets and endocrine activities. As a photodegradation product from a common estrogen of both human and farm animal origin, lumiestrone is present in sewage effluent, is produced from estrone upon exposure to natural sunlight and should be considered as a new environmental contaminant.

Generation of the first BAC-based physical map of the common carp genome

Background

Common carp (Cyprinus carpio), a member of Cyprinidae, is the third most important aquaculture species in the world with an annual global production of 3.4 million metric tons, accounting for nearly 14% of the all freshwater aquaculture production in the world. Apparently genomic resources are needed for this species in order to study its performance and production traits. In spite of much progress, no physical maps have been available for common carp. The objective of this project was to generate a BAC-based physical map using fluorescent restriction fingerprinting.

Result

The first generation of common carp physical map was constructed using four- color High Information Content Fingerprinting (HICF). A total of 72,158 BAC clones were analyzed that generated 67,493 valid fingerprints (5.5 × genome coverage). These BAC clones were assembled into 3,696 contigs with the average length of 476 kb and a N50 length of 688 kb, representing approximately 1.76 Gb of the common carp genome. The largest contig contained 171 BAC clones with the physical length of 3.12 Mb. There are 761 contigs longer than the N50, and these contigs should be the most useful resource for future integrations with linkage map and whole genome sequence assembly. The common carp physical map is available at http://genomics.cafs.ac.cn/fpc/WebAGCoL/Carp/WebFPC/.

Conclusion

The reported common carp physical map is the first physical map of the common carp genome. It should be a valuable genome resource facilitating whole genome sequence assembly and characterization of position-based genes important for aquaculture traits.

A genomic and ecotoxicological perspective of DNA array studies in aquatic environmental risk assessment

Ecotoxicogenomics is developing into a key tool for the assessment of environmental impacts and environmental risk assessment for aquatic ecosystems. This review aims to report achievements and drawbacks of this technique and to explore potential conceptual and experimental procedures to improve future investigations. Ecotoxicogenomic literature evidences the ability of genomic technologies to characterize toxicant specific gene transcriptome patterns that can be used to identify major toxicants affecting aquatic species. They also contribute decisively to the development of new molecular biomarkers and, in many cases, to the determination of new possible gene targets. Primary transcriptomic responses obtained after short exposures provided more information of putative gene targets than secondary responses obtained after long, chronic exposures, which in turn are usually more accurate to describe actual environmental impacts in natural populations. Several problems need to be addressed in future investigations: the lack of studies (and genomic information) on key ecological species and taxa, the need to better understand the different transcriptomic responses to high and low doses and, especially, short and long exposures, and the need to improve experimental designs to minimize false transcriptome interpretations of target genes.

Genes are differentially expressed at transcriptional level of Neocaridina denticulata following short-term exposure to nonylphenol

To assess the toxicity of nonylphenol towards aquatic crustaceans, Neocaridina denticulata were exposed short-term to sublethal concentration (0.001-0.5 mg/L). Following treatment, differentially expressed genes were identified using suppression subtractive hybridization on samples prepared from whole specimens. There were 20 differentially expressed sequence tags that corresponded to known genes and could be divided into six functional classes: defence, translation, metabolism, ribosomal gene expression, respiration, and genes involved in the stress response. Using semi-quantitative RT-PCR, we found that 14 of the differentially expressed sequence tags significantly responded to nonylphenol, including six at a nominal concentration of 0.01 mg/L; among them, 12 genes were down-regulated. These results suggest that under non-lethal concentrations of nonylphenol, the polluted aquatic environment may still present a potential risk to N. denticulata.

Construction and characterization of the BAC library for common carp Cyprinus carpio L. and establishment of microsynteny with zebrafish Danio rerio

A bacterial artificial chromosome (BAC) library of common carp Cyprinus carpio L. was constructed as a part of ongoing common carp genome project, which is aiming assembly of common carp genome. The library, containing a total of 92,160 BAC clones with an average insert size of 141 kb, was constructed into the restriction site of Hind III on BAC vector CopyControl pCC1BAC, covering 7.7 X haploid genome equivalents. Three dimension pools and superpools of the BAC library were established and 23 positive clones of 14 targets were identified from one-fifth of the BAC library. Pilot project of BAC end sequencing was conducted on 2,688 BAC ends from 1,344 clones and harvested 2,522 high-quality Q20 sequences with average length of 677 bp. The sequencing success rate was 93.8% and pair-end success rate was 92.3%. A total of 212 microsyntenies had been established between common carp and zebrafish genomes as a trial for genome-wide comparative genomics in these two closely related species.

Endocrine disrupting compounds: can they target the immune system of fish?

Endocrine disruption, in particular disruption by estrogen-active compounds, has been identified as an important ecotoxicological hazard in the aquatic environment. Research on the impact of endocrine disrupting compounds (EDCs) on wildlife has focused on disturbances of the reproductive system. However, there is increasing evidence that EDCs affect a variety of physiological systems other than the reproductive system. Here, we discuss if EDCs may be able to affect the immune system of fish, as this would have direct implications for individual fitness and population growth. Evidence suggesting an immunomodulatory role of estrogens in fish comes from the following findings: (a) estrogen receptors are expressed in piscine immune organs, (b) immune gene expression is modulated by estrogen exposure, and (c) pathogen susceptibility of fish increases under estrogen exposure.

Methoxychlor affects multiple hormone signaling pathways in the largemouth bass (Micropterus salmoides) liver

Methoxychlor (MXC) is an organochlorine pesticide that has been shown to have estrogenic activity by activating estrogen receptors and inducing vitellogenin production in male fish. Previous studies report that exposure to MXC induces changes in mRNA abundance of reproductive genes in the liver and testes of largemouth bass (Micropterus salmoides). The objective of the present study was to better characterize the mode of action of MXC by measuring the global transcriptomic response in the male largemouth liver using an oligonucleotide microarray. Microarray analysis identified highly significant changes in the expression of 37 transcripts (p<0.001) (20 induced and 17 decreased) in the liver after MXC injection and a total of 900 expression changes (p<0.05) in transcripts with high homology to known genes. Largemouth bass estrogen receptor alpha (esr1) and androgen receptor (ar) were among the transcripts that were increased in the liver after MXC treatment. Functional enrichment analysis identified the molecular functions of steroid binding and androgen receptor activity as well as steroid hormone receptor activity as being significantly over-represented gene ontology terms. Pathway analysis identified c-fos signaling as being putatively affected through both estrogen and androgen signaling. This study provides evidence that MXC elicits transcriptional effects through the estrogen receptor as well as androgen receptor-mediated pathways in the liver.

Systems biology: leading the revolution in ecotoxicology

The rapid development of new technologies such as transcriptomics, proteomics, and metabolomics (Omics) are changing the way ecotoxicology is practiced. The data deluge has begun with genomes of over 65 different aquatic species that are currently being sequenced, and many times that number with at least some level of transcriptome sequencing. Integrating these top-down methodologies is an essential task in the field of systems biology. Systems biology is a biology-based interdisciplinary field that focuses on complex interactions in biological systems, with the intent to model and discover emergent properties of the system. Recent studies demonstrate that Omics technologies provide valuable insight into ecotoxicity, both in laboratory exposures with model organisms and with animals exposed in the field. However, these approaches require a context of the whole animal and population to be relevant. Powerful approaches using reverse engineering to determine interacting networks of genes, proteins, or biochemical reactions are uncovering unique responses to toxicants. Modeling efforts in aquatic animals are evolving to interrelate the interacting networks of a system and the flow of information linking these elements. Just as is happening in medicine, systems biology approaches that allow the integration of many different scales of interaction and information are already driving a revolution in understanding the impacts of pollutants on aquatic systems.

An integrated transcriptomic and proteomic approach characterizing estrogenic and metabolic effects of 17 alpha-ethinylestradiol in zebrafish (Danio rerio)

Nowadays there is much concern about the presence of endocrine disrupting compounds (EDCs) in the environment due to their ability to interfere with the endocrine system. In the presented study, adult zebrafish (Danio rerio) were exposed to 30 ng L(-1) 17alpha-ethinylestradiol (EE2) for 4 and 28 days. The underlying molecular mechanisms of EE2 were studied in the zebrafish liver by applying a combined transcriptomics and proteomics approach. In addition, we assessed the added value of such an integrated-omics approach. Oligo microarrays, spotted with 3479 zebrafish-specific oligos, were employed to generate differential gene expression levels. The proteomic responses were evaluated by means of differential in-gel electrophoresis (DiGE), combined with MALDI-tandem mass spectrometry. Assessment of the major biological functions of the differentially expressed transcripts and proteins illustrated that both individual platforms could profile a clear estrogenic interference, next to numerous metabolism-related effects and stress responses. Cross-comparison of both transcriptomics and proteomics datasets displayed limited concordance, though, thorough revision of the results illustrated that transcriptional effects were projected on protein level as downstream effects of affected signalling pathways. Overall, this study demonstrated that a proteomics approach can lift the biological interpretation of microarrays to a higher level, and moreover, opens a window for identification of possible new biomarkers.

Towards functional genomics in fish using quantitative proteomics

Microarray and gene expression analysis have been key in our understanding of molecular pathways underlying physiological responses. Arguably, a large number of microarray based studies in fish have examined steroid nuclear receptor signaling (e.g., estrogens, androgens) in the context of both physiology and toxicology. Following close behind the advances in gene expression analysis, novel proteomic tools are available that have been under utilized in fish endocrinology studies. Quantitative proteomic approaches include both gel based (e.g., 2D gel electrophoresis, 2-D Fluorescence Difference Gel Electrophoresis; DIGE) and non-gel based methods that can be separated further into labeling approaches such as stable isotope labeling (SILAC), isotope coded affinity tags (ICAT), and isobaric tagging (iTRAQ) and label-free approaches (e.g., spectral counting and absolute quantitation). This review summarizes quantitative proteomic approaches and describes a successful application of iTRAQ) to study changes in the liver proteome in fathead minnows in response to the androgen, 17beta-trenbolone. The challenge remains to integrate molecular datasets in such a manner as to be able to consider temporal effects and complex regulation at the level of the genome and proteome.

Stimulation of growth and changes in the hepatic transcriptome by 17β-estradiol in the yellow perch (Perca flavescens)

The effects of dietary 17β-estradiol (E2) on growth and liver transcriptomics were investigated in the yellow perch ( Perca flavescens). After a 3-mo treatment, E2 significantly stimulated an increase in length and weight of juvenile male and female perch relative to control animals. The increase was significantly greater in females compared with males. Separate, unnormalized cDNA libraries were constructed from equal quantities of RNA from 6 male and 6 female livers of E2-treated and control perch, and 3,546 and 3,719 expressed sequence tags (ESTs) were obtained, respectively. To characterize E2-regulated transcripts, EST frequencies between libraries were calculated within contiguous sequences that were assembled from the combined ESTs of both libraries. Frequencies were also determined in EST transcript groupings produced by aligning all of the ESTs from both libraries at the nucleotide level. From these analyses, there were 28 annotated transcripts that were regulated by 75% between libraries and for which there were at least 5 ESTs of the same transcript between libraries. Regulation of a subset ( 14 ) of these transcripts was confirmed by quantitative reverse transcription-polymerase chain reaction (QPCR). Transcripts that were upregulated by E2 included reproduction-related proteins, binding proteins, and proteases and protease inhibitors. While not part of the transcript frequency analysis, QPCR showed significant upregulation of estrogen receptor esr1 and of insulin-like growth factor I (IGF-I) in E2 livers. E2-downregulated transcripts represented a variety of functional categories including components of the respiratory chain, lipid transport and metabolism, glycolysis, amino acid and nitrogen metabolism, binding proteins, a hydrolytic enzyme, and a transcriptional regulator. In perch it appears that exogenous estrogen drastically shifts liver metabolism toward the production of lipoproteins and carbohydrate binding proteins, and that the growth-promoting action may involve an increase in hepatic IGF-I production.

Effects of environmental relevant doses of pollutants from offshore oil production on Atlantic cod (Gadus morhua)

The release of produced water (PW), a by-product of offshore oil production, has increased in Norwegian waters in recent years. Alkylphenols (AP), a major component of PW, have been shown to have endocrine disrupting effects on several fish species. In the present study, four groups of female Atlantic cod (Gadus morhua) were orally exposed for 20 weeks to two different concentrations of a mixture of C4-C7 APs, PW or 17beta-estradiol. The transcriptional responses in the liver of Atlantic female cod were studied using a custom-made cDNA microarray. The largest transcriptional effects were seen in cod exposed to the lowest dose of APs. Several biological processes such as glycolysis, apoptosis and the general stress response were affected by exposure to APs. In addition, genes coding for the detoxification enzymes CYP1A and sulfotransferase 2 were up-regulated in the low exposure group. Significant reduction in gonadosomatic index (GSI) and the concentration of plasma vitellogenin were seen in both AP and 17beta-estradiol exposed cod. Exposure to PW had little effect on GSI and the regulation of stress responsive genes. The findings indicate that chronic exposure to low levels of APs may cause a stress response and delayed maturation in female cod.

A critical analysis of the biological impacts of plasticizers on wildlife

This review provides a critical analysis of the biological effects of the most widely used plasticizers, including dibutyl phthalate, diethylhexyl phthalate, dimethyl phthalate, butyl benzyl phthalate and bisphenol A (BPA), on wildlife, with a focus on annelids (both aquatic and terrestrial), molluscs, crustaceans, insects, fish and amphibians. Moreover, the paper provides novel data on the biological effects of some of these plasticizers in invertebrates, fish and amphibians. Phthalates and BPA have been shown to affect reproduction in all studied animal groups, to impair development in crustaceans and amphibians and to induce genetic aberrations. Molluscs, crustaceans and amphibians appear to be especially sensitive to these compounds, and biological effects are observed at environmentally relevant exposures in the low ng l(-1) to microg l(-1) range. In contrast, most effects in fish (except for disturbance in spermatogenesis) occur at higher concentrations. Most plasticizers appear to act by interfering with the functioning of various hormone systems, but some phthalates have wider pathways of disruption. Effect concentrations of plasticizers in laboratory experiments coincide with measured environmental concentrations, and thus there is a very real potential for effects of these chemicals on some wildlife populations. The most striking gaps in our current knowledge on the impacts of plasticizers on wildlife are the lack of data for long-term exposures to environmentally relevant concentrations and their ecotoxicity when part of complex mixtures. Furthermore, the hazard of plasticizers has been investigated in annelids, molluscs and arthropods only, and given the sensitivity of some invertebrates, effects assessments are warranted in other invertebrate phyla.

Revealing genes associated with vitellogenesis in the liver of the zebrafish (Danio rerio) by transcriptome profiling

Background

In oviparous vertebrates, including fish, vitellogenesis consists of highly regulated pathways involving 17β-estradiol (E2). Previous studies focused on a relatively small number of hepatic expressed genes during vitellogenesis. This study aims to identify hepatic genes involved in vitellogenesis and regulated by E2, by using zebrafish microarray gene expression profiling, and to provide information on functional distinctive genes expressed in the liver of a vitellogenic female, using zebrafish as a model fish.

Results

Genes associated with vitellogenesis were revealed by the following paired t-tests (SAM) comparisons: a) two-month old vitellogenic (Vit2) females were compared with non-vitellogenic (NV) females, showing 825 differentially expressed transcripts during early stages of vitellogenesis, b) four-month old vitellogenic (Vit4) females were compared with NV females, showing 1,046 differentially expressed transcripts during vitellogenesis and c) E2-treated males were compared with control males, showing 1,828 differentially expressed transcripts regulated by E2. A Venn diagram revealed 822 common transcripts in the three groups, indicating that these transcripts were involved in vitellogenesis and putatively regulated by E2. In addition, 431 transcripts were differentially expressed in Vit2 and Vit4 females but not in E2-treated males, indicating that they were putatively not up-regulated by E2. Correspondence analysis showed high similarity in expression profiles of Vit2 with Vit4 and of NV females with control males. The E2-treated males differed from the other groups. The repertoire of genes putatively regulated by E2 in vitellogenic females included genes associated with protein synthesis and reproduction. Genes associated with the immune system processes and biological adhesion, were among the genes that were putatively not regulated by E2. E2-treated males expressed a large array of transcripts that were not associated with vitellogenesis.

The study revealed several genes that were not reported before as being regulated by E2. Also, the hepatic expression of several genes was reported here for the first time.

Conclusion

Gene expression profiling of liver samples revealed 1,046 differentially expressed transcripts during vitellogenesis of which at least ~64% were regulated by E2. The results raise the question on the regulation pattern and temporal pleiotropic expression of hepatic genes in vitellogenic females.

Construction of a robust microarray from a non-model species (largemouth bass) using pyrosequencing technology

A novel custom microarray for largemouth bass (Micropterus salmoides) was designed with sequences obtained from a normalized cDNA library using the 454 Life Sciences GS-20 pyrosequencer. This approach yielded in excess of 58 million bases of high-quality sequence. The sequence information was combined with 2,616 reads obtained by traditional suppressive subtractive hybridizations to derive a total of 31,391 unique sequences. Annotation and coding sequences were predicted for these transcripts where possible. 16,350 annotated transcripts were selected as target sequences for the design of the custom largemouth bass oligonucleotide microarray. The microarray was validated by examining the transcriptomic response in male largemouth bass exposed to 17beta-oestradiol. Transcriptomic responses were assessed in liver and gonad, and indicated gene expression profiles typical of exposure to oestradiol. The results demonstrate the potential to rapidly create the tools necessary to assess large scale transcriptional responses in non-model species, paving the way for expanded impact of toxicogenomics in ecotoxicology.

Identification of WSB1 gene as an important regulator in the development of zebrafish embryo during midblastula transition

To uncover novel genes potentially involved in embryo development, especially at the midblastula transition (MBT) phase in the developing embryo of zebrafish, Affymetrix zebrafish GeneChip microarray analysis was carried out on the expression of 14,900 gene transcripts. The results of the analysis showed that 360 genes were clearly up-regulated and 119 genes were markedly down-regulated. Many of these genes were involved in transcription factor activity, nucleic acid binding, and cell growth. The present study showed that significant changes in transcript abundance occurred during the MBT phase. The expression of eight of these 479 genes was identified by reverse transcription-polymerase chain reaction analysis, confirming the microarray results. The WSB1 gene, found to be down-regulated by the microarray and reverse transcription-polymerase chain reaction analyses, was selected for further study. Sequence analysis of the WSB1 gene showed that it encodes a protein with 75% identity to the corresponding active human orthologs. In addition, WSB1 gene expression was detected at a higher level at 2 h post fertilization and at a lower level at 4 h post fertilization, consistent with the chip results. Overexpression of the WSB1 gene can result in the formation of abnormalities in embryos, as determined by fluorescence-activated cell sorting. The present study showed unequivocally that the occurrence of WSB1 expression is an important event during the MBT phase in the development of zebrafish embryos.