Comparative toxicological assessment of two bisphenols using a systems approach: Evaluation of the behavioral and transcriptomic responses of Danio rerio to Bisphenol A and Tetrabromobisphenol A.

The zebrafish (Danio rerio) is becoming a critical component of New Approach Methods (NAMs) in chemical risk assessment. As a whole organism in vitro NAM, the zebrafish model offers significant advantages over individual cell-line testing, including toxicokinetic and toxicodynamic competencies. A transcriptomic approach not only allows for insight into mechanism of action for both apical endpoints and unobservable adverse outcomes, but also changes in gene expression induced by lower, environmentally relevant concentrations. In this study, we used a larval zebrafish model to assess the behavioral and transcriptomic alterations caused by sub-phenotypic concentrations of two chemicals with the same structural backbone, the endocrine disrupting chemicals: Bisphenol A and Tetrabromobisphenol A. Following assessment of behavioral toxicity, we used a transcriptomic approach to identify molecular pathways associated with previously described phenotypes. We also determined the transcriptomic Point of Departure (POD) for each chemical by modelling gene expression changes as continuous systems which allows for the identification of a single concentration at which toxic effects can be predicted. This can then be investigated with confirmatory cell-based testing in an integrated approach to testing and assessment (IATA) to determine risk to human health and the environment with greater confidence. This paper demonstrates the impact of using a multi-faceted approach for evaluating the physiological and neurotoxic effects of exposure to structurally related chemicals. By comparing phenotypic effects with transcriptomic outcomes, we were able to differentiate, characterize and rank the toxicities of related bisphenols, which demonstrates methodological advantages unique to the larval zebrafish NAM.

The contribution of larval zebrafish transcriptomics to chemical risk assessment.

In Canada, the Canadian Environmental Protection Act (1999) requires human health and environmental risk assessments be conducted for new substances prior to their manufacture or import. While this toxicity data is historically obtained using rodents, in response to the international effort to eliminate animal testing, Health Canada is collaborating with the National Research Council (NRC) of Canada to develop a New Approach Method by refining existing NRC zebrafish models. The embryo/larval zebrafish model evaluates systemic (whole body) general toxicity which is currently unachievable with cell-based testing. The model is strengthened using behavioral, toxicokinetic and transcriptomic responses to assess non-visible indicators of toxicity following chemical exposure at sub-phenotypic concentrations. In this paper, the predictive power of zebrafish transcriptomics is demonstrated using two chemicals; Raloxifene and Resorcinol. Raloxifene exposure produced darkening of the liver and malformation of the nose/mandible, while Resorcinol exposure produced increased locomotor activity. Transcriptomic analysis correlated differentially expressed genes with the phenotypic effects and benchmark dose calculations determined that the transcriptomic Point of Departure (POD) occurred at subphenotypic concentrations. Correlating gene expression with apical (phenotypic) effects strengthens confidence in evaluation of chemical toxicity, thereby demonstrating the significant advancement that the larval zebrafish transcriptomics model represents in chemical risk assessment.

Fatty acid-binding protein genes of the ancient, air-breathing, ray-finned fish, spotted gar (Lepisosteus oculatus).

With the advent of high-throughput DNA sequencing technology, the genomic sequence of many disparate species has led to the relatively new discipline of genomics, the study of genome structure, function and evolution. Much work has been focused on the role of whole genome duplications (WGD) in the architecture of extant vertebrate genomes, particularly those of teleost fishes which underwent a WGD early in the teleost radiation >230 million years ago (mya). Our past work has focused on the fate of duplicated copies of a multigene family coding for the intracellular lipid-binding protein (iLBP) genes in the teleost fishes. To define the evolutionary processes that determined the fate of duplicated genes and generated the structure of extant fish genomes, however, requires comparative genomic analysis with a fish lineage that diverged before the teleost WGD, such as the spotted gar (Lepisosteus oculatus), an ancient, air-breathing, ray-finned fish. Here, we describe the genomic organization, chromosomal location and tissue-specific expression of a subfamily of the iLBP genes that code for fatty acid-binding proteins (Fabps) in spotted gar. Based on this work, we have defined the minimum suite of fabp genes prior to their duplication in the teleost lineages ~230-400 mya. Spotted gar, therefore, serves as an appropriate outgroup, or ancestral/ancient fish, that did not undergo the teleost-specific WGD. As such, analyses of the spatio-temporal regulation of spotted gar genes provides a foundation to determine whether the duplicated fabp genes have been retained in teleost genomes owing to either sub- or neofunctionalization.

Evolution of the duplicated intracellular lipid-binding protein genes of teleost fishes.

Increasing organismal complexity during the evolution of life has been attributed to the duplication of genes and entire genomes. More recently, theoretical models have been proposed that postulate the fate of duplicated genes, among them the duplication-degeneration-complementation (DDC) model. In the DDC model, the common fate of a duplicated gene is lost from the genome owing to nonfunctionalization. Duplicated genes are retained in the genome either by subfunctionalization, where the functions of the ancestral gene are sub-divided between the sister duplicate genes, or by neofunctionalization, where one of the duplicate genes acquires a new function. Both processes occur either by loss or gain of regulatory elements in the promoters of duplicated genes. Here, we review the genomic organization, evolution, and transcriptional regulation of the multigene family of intracellular lipid-binding protein (iLBP) genes from teleost fishes. Teleost fishes possess many copies of iLBP genes owing to a whole genome duplication (WGD) early in the teleost fish radiation. Moreover, the retention of duplicated iLBP genes is substantially higher than the retention of all other genes duplicated in the teleost genome. The fatty acid-binding protein genes, a subfamily of the iLBP multigene family in zebrafish, are differentially regulated by peroxisome proliferator-activated receptor (PPAR) isoforms, which may account for the retention of iLBP genes in the zebrafish genome by the process of subfunctionalization of cis-acting regulatory elements in iLBP gene promoters.

Tissue-specific transcriptional modulation of fatty acid-binding protein genes, fabp2, fabp3 and fabp6, by fatty acids and the peroxisome proliferator, clofibrate, in zebrafish (Danio rerio).

All fabp genes, except fabp2, fabp3 and fabp6, exist as duplicates in the zebrafish genome owing to a whole genome duplication event ~230-400 million years ago. Transcription of some duplicated fabp genes is modulated by fatty acids (FAs) and/or clofibrate, a peroxisome proliferator-activated receptor (PPAR) agonist. We had also shown previously that the steady-state level of acyl-CoA oxidase 1 (acox1) mRNA, a marker of PPARα activation, was elevated in liver, intestine, heart and muscle of fish fed clofibrate demonstrating that zebrafish, unlike some fishes, is responsive to this drug. acox1 transcripts were not induced in the brain of fish fed clofibrate, which suggests this drug may not cross the blood brain barrier. Here, we investigated the effect of dietary FAs and clofibrate on the transcription of single copy fabp genes, fabp2, fabp3 and fabp6, in five tissues of inbred zebrafish. The steady-state level of fabp2 transcripts increased in intestine, while fabp3 mRNA increased in liver of fish fed diets differing in FA content. In fish fed clofibrate, fabp3 mRNA in intestine, and fabp6 mRNA in intestine and heart, were elevated. Based on these findings, modulation of fabp2, fabp3 and fabp6 transcription by FAs and/or clofibrate in zebrafish implicates control of these genes by PPAR interaction with peroxisome proliferator response elements (PPRE) most likely in fabp promoters. Moreover, transcriptional induction of these fabp genes by dietary FAs and/or clofibrate is over-ridden by a tissue-specific mechanism(s), e.g., transcriptional activator or repressor proteins.

Comparative genomics and evolutionary diversification of the duplicated fabp6a and fabp6b genes in medaka and three-spined stickleback.

We describe the evolutionary diversification of the duplicated ileal fatty acid-binding protein genes (fabp6a and fabp6b) from Japanese ricefish (Oryzias latipes; medaka) and three-spined stickleback (Gasterosteus aculeatus). The fabp6a and fabp6b genes from medaka and three-spined stickleback encode polypeptides of 125-127 amino acids, which share highest sequence identity with their orthologs in teleost fishes and tetrapods. All Fabp6a and Fabp6b from different species cluster together in a distinct clade in phylogenetic analysis and the topology of the tree suggests that fabp6a and fabp6b from medaka and three-spined stickleback are most likely duplicated genes of an ancestral FABP6 owing to teleost-specific whole-genome duplication. However, the topology of an alternate phylogenetic tree revealed that the duplication of the ancestral FABP6 that gave rise to the extant fabp6a and fabp6b possibly occurred before the divergence of tetrapods and fishes. Conserved gene synteny was evident between the teleost fabp6a and fabp6b genes and the human FABP6 gene. The tissue-specific distribution of fabp6a transcripts suggests the retention of ancestral function(s) of the fabp6a gene in medaka and three-spined stickleback with acquisition of new function(s) in different tissues. However, the tissue-specific regulation of the fabp6b gene has diverged markedly in medaka and three-spined stickleback since the duplication of the fabp6 gene.

Tissue-specific differential induction of duplicated fatty acid-binding protein genes by the peroxisome proliferator, clofibrate, in zebrafish (Danio rerio).

BACKGROUND: Force, Lynch and Conery proposed the duplication-degeneration-complementation (DDC) model in which partitioning of ancestral functions (subfunctionalization) and acquisition of novel functions (neofunctionalization) were the two primary mechanisms for the retention of duplicated genes. The DDC model was tested by analyzing the transcriptional induction of the duplicated fatty acid-binding protein (fabp) genes by clofibrate in zebrafish. Clofibrate is a specific ligand of the peroxisome proliferator-activated receptor (PPAR); it activates PPAR which then binds to a peroxisome proliferator response element (PPRE) to induce the transcriptional initiation of genes primarily involved in lipid homeostasis. Zebrafish was chosen as our model organism as it has many duplicated genes owing to a whole genome duplication (WGD) event that occurred ~230-400 million years ago in the teleost fish lineage. We assayed the steady-state levels of fabp mRNA and heterogeneous nuclear RNA (hnRNA) transcripts in liver, intestine, muscle, brain and heart for four sets of duplicated fabp genes, fabp1a/fabp1b.1/fabp1b.2, fabp7a/fabp7b, fabp10a/fabp10b and fabp11a/fabp11b in zebrafish fed different concentrations of clofibrate. RESULT: Electron microscopy showed an increase in the number of peroxisomes and mitochondria in liver and heart, respectively, in zebrafish fed clofibrate. Clofibrate also increased the steady-state level of acox1 mRNA and hnRNA transcripts in different tissues, a gene with a functional PPRE. These results demonstrate that zebrafish is responsive to clofibrate, unlike some other fishes. The levels of fabp mRNA and hnRNA transcripts for the four sets of duplicated fabp genes was determined by reverse transcription, quantitative polymerase chain reaction (RT-qPCR). The level of hnRNA coded by a gene is an indirect estimate of the rate of transcriptional initiation of that gene. Clofibrate increased the steady-state level of fabp mRNAs and hnRNAs for both the duplicated copies of fabp1a/fabp1b.1, and fabp7a/fabp7b, but in different tissues. Clofibrate also increased the steady-state level of fabp10a and fabp11a mRNAs and hnRNAs in liver, but not for fabp10b and fabp11b. CONCLUSION: Some duplicated fabp genes have, most likely, retained PPREs, but induction by clofibrate is over-ridden by an, as yet, unknown tissue-specific mechanism(s). Regardless of the tissue-specific mechanism(s), transcriptional control of duplicated zebrafish fabp genes by clofibrate has markedly diverged since the WGD event.

The evolutionary relationship of the transcriptionally active fabp11a (intronless) and fabp11b genes of medaka with fabp11 genes of other teleost fishes.

Here we describe the structure of the fatty acid-binding protein 11a and 11b genes (fabp11a and fabp11b) in medaka, and their evolutionary relationship to fabp11 genes from other teleost fishes. Initial studies indicated that the medaka fabp11a gene is intronless, but the fabp11b gene consists of four exons separated by three introns, a genomic organization that is characteristic of most members of the intracellular lipid-binding protein family. Based on genomic sequence, we conclude that the intronless fabp11a gene most likely arose as a result of reverse transcription of its mRNA transcript into cDNA followed by integration into chromosomal DNA. The ancestral intron-containing fabp11a gene was presumably lost from the medaka genome. The duplicated fabp11 genes extant in medaka encode polypeptides of 134 amino acids, which share highest sequence identity and similarity, and cluster in a distinct phylogenetic clade, with their orthologs in other teleost fishes. The fabp11a and fabp11b genes in medaka are therefore orthologs of the fabp11a and fabp11b genes, respectively, of other teleost fishes. No conserved gene synteny was found between medaka fabp11a and fabp11a genes from other teleost fishes, supporting our suggestion as to how this intronless gene arose. However, conserved gene synteny was evident between medaka fabp11b and fabp11b genes from other teleost fishes. The tissue-specific distribution of transcripts for medaka and zebrafish fabp11a and fabp11b genes revealed acquisition of a new function(s) in various tissues by the medaka fabp11b gene, which may explain the retention of sister duplicates of fabp11 in the medaka genome.

The duplicated retinol-binding protein 7 (rbp7) genes are differentially transcribed in embryos and adult zebrafish (Danio rerio).

Genomic and cDNA sequences coding for two cellular retinol-binding proteins (rbp) in zebrafish were retrieved from DNA sequence databases. Phylogenetic analysis revealed that these proteins were most similar to mammalian RBP7/Rbp7 proteins. Hence, the genes coding for these proteins were named rbp7a and rbp7b. Using a radiation hybrid panel, rbp7a and rbp7b were mapped to the zebrafish chromosomes 23 and 6, respectively. Conserved gene synteny indicated that these genes most likely arose as a result of a fish-specific whole-genome duplication event that had occurred 230-400 million years ago. Whole-mount in situ hybridization to zebrafish embryos detected rbp7a transcripts from the sphere stage (4h post-fertilization (hpf)) in the forerunner cells and the yolk syncytial layer, as well as in Kuppfer's vesicle and the periderm at 12 hpf. The transcripts of rbp7b were seen primarily in the somite stages (10-24 hpf) of zebrafish embryos, but also in the floor plate and hypochord, and did not overlap with the distribution of rbp7a transcripts in embryos. The hybridization signal for rbp7a and rbp7b transcripts was not detected in embryos after 12 hpf and 24 hpf, respectively. While transcripts for both rbp7a and rbp7b were found in all adult tissues assayed by RT-qPCR, the steady-state level of rbp7a transcripts were significantly higher than that of rbp7b transcripts in gill and ovary, whereas rbp7b transcripts were significantly higher than rbp7a transcripts in muscle and brain. The distribution of rbp7a and rbp7b transcripts in embryos and adult zebrafish indicate that the cis-elements that control the transcriptional regulation of the rbp7a and rbp7b genes have diverged considerably since their duplication.

Differential tissue-specific distribution of transcripts for the duplicated fatty acid-binding protein 10 (fabp10) genes in embryos, larvae and adult zebrafish (Danio rerio).

Genomic and cDNA sequences coding for a fatty acid-binding protein (FABP) in zebrafish were retrieved from DNA sequence databases. The cDNA codes for a protein of 14.7 kDa (pI = 5.94), and the gene consists of four exons, properties characteristic of most vertebrate FABP genes. Phylogenetic analyses using vertebrate FABPs indicated that this protein is most similar to zebrafish Fabp10. Currently, only one fabp10 gene is annotated in the zebrafish genome. In this article, the notations 'fabp10a' and 'fabp10b' are used to refer to the duplicate copies of fabp10. The zebrafish fabp10a and fabp10b genes were assigned by radiation hybrid mapping to chromosomes 16 and 19, respectively. On the basis of conserved gene synteny with chicken FABP10 on chromosome 23, zebrafish fabp10a and fabp10b are duplicates resulting from a whole-genome duplication event early in the ray-finned fish lineage some 230-400 million years ago. Whole-mount in situ hybridization detected fabp10b transcripts only in the olfactory vesicles of embryos and larvae, whereas fabp10a transcripts have been shown previously to be present only in the liver of embryos and larvae. In adults, RT-PCR detected fabp10b transcripts in all tissues assayed. By contrast, fabp10a transcripts were detected only in adult liver, intestine and testis. This differential tissue distribution of transcripts for the duplicated fabp10 genes suggests considerable divergence of their cis-acting regulatory elements since their duplication.