Gene expression profiling in the neuroendocrine brain of male goldfish (Carassius auratus) exposed to 17alpha-ethinylestradiol

Gene expression responses in male fathead minnows exposed to binary mixtures of an estrogen and antiestrogen

Background

Aquatic organisms are continuously exposed to complex mixtures of chemicals, many of which can interfere with their endocrine system, resulting in impaired reproduction, development or survival, among others. In order to analyze the effects and mechanisms of action of estrogen/anti-estrogen mixtures, we exposed male fathead minnows (Pimephales promelas) for 48 hours via the water to 2, 5, 10, and 50 ng 17α-ethinylestradiol (EE₂)/L, 100 ng ZM 189,154/L (a potent antiestrogen known to block activity of estrogen receptors) or mixtures of 5 or 50 ng EE₂/L with 100 ng ZM 189,154/L. We analyzed gene expression changes in the gonad, as well as hormone and vitellogenin plasma levels.

Results

Steroidogenesis was down-regulated by EE₂ as reflected by the reduced plasma levels of testosterone in the exposed fish and down-regulation of genes in the steroidogenic pathway. Microarray analysis of testis of fathead minnows treated with 5 ng EE₂/L or with the mixture of 5 ng EE₂/L and 100 ng ZM 189,154/L indicated that some of the genes whose expression was changed by EE₂ were blocked by ZM 189,154, while others were either not blocked or enhanced by the mixture, generating two distinct expression patterns. Gene ontology and pathway analysis programs were used to determine categories of genes for each expression pattern.

Conclusion

Our results suggest that response to estrogens occurs via multiple mechanisms, including canonical binding to soluble estrogen receptors, membrane estrogen receptors, and other mechanisms that are not blocked by pure antiestrogens.

Quantitative proteomic profiles of androgen receptor signaling in the liver of fathead minnows (Pimephales promelas)

Androgenic chemicals are present in the environment at concentrations that impair reproductive processes in fish. The objective of this experiment was to identify proteins and cell processes mediated through androgen receptor signaling using an androgen receptor agonist (17beta-trenbolone) and antagonist (flutamide) in the liver. Female fathead minnows were exposed to nominal concentrations of either 17beta-trenbolone (0.05, 0.5, or 5 microg/L), flutamide (50, 150, or 500 microg/L), or a mixture (500 microg flutamide/L and 0.5 microg 17beta-trenbolone/L) for 48 h. The iTRAQ method was used to label peptides after protein extraction and trypsin-digestion from livers of untreated controls or from fish treated with 17beta-trenbolone (5 microg/L), flutamide (500 microg/L), or a mixture of both compounds. Forty-five proteins were differentially altered by one or more treatments (p<0.05). Many altered proteins were involved in cellular metabolism (e.g., glyceraldehyde 3-phosphate dehydrogenase, phosphoglycerate mutase), general and oxidative stress response (e.g., superoxide dismutase and heat shock proteins), and the regulation of translation (e.g., ribosomal proteins). Cellular pathway analysis identified additional signaling cascades activated or inhibited by flutamide that may not be androgen receptor mediated. We also compared changes in select proteins to changes in their mRNA levels and observed, in general, that proteins and mRNA changes did not correlate, suggesting complex regulation at the level of both the transcriptome and proteome. It is concluded that both transcriptomic and proteomic approaches offer unique and complementary insights into mechanisms of regulation. We demonstrate the utility of proteomic profiling for use on a model species with value to ecotoxicology but having limited genomic information.

Defining global neuroendocrine gene expression patterns associated with reproductive seasonality in fish

BACKGROUND

Many vertebrates, including the goldfish, exhibit seasonal reproductive rhythms, which are a result of interactions between external environmental stimuli and internal endocrine systems in the hypothalamo-pituitary-gonadal axis. While it is long believed that differential expression of neuroendocrine genes contributes to establishing seasonal reproductive rhythms, no systems-level investigation has yet been conducted.

METHODOLOGY/PRINCIPAL FINDINGS

In the present study, by analyzing multiple female goldfish brain microarray datasets, we have characterized global gene expression patterns for a seasonal cycle. A core set of genes (873 genes) in the hypothalamus were identified to be differentially expressed between May, August and December, which correspond to physiologically distinct stages that are sexually mature (prespawning), sexual regression, and early gonadal redevelopment, respectively. Expression changes of these genes are also shared by another brain region, the telencephalon, as revealed by multivariate analysis. More importantly, by examining one dataset obtained from fish in October who were kept under long-daylength photoperiod (16 h) typical of the springtime breeding season (May), we observed that the expression of identified genes appears regulated by photoperiod, a major factor controlling vertebrate reproductive cyclicity. Gene ontology analysis revealed that hormone genes and genes functionally involved in G-protein coupled receptor signaling pathway and transmission of nerve impulses are significantly enriched in an expression pattern, whose transition is located between prespawning and sexually regressed stages. The existence of seasonal expression patterns was verified for several genes including isotocin, ependymin II, GABA(A) gamma2 receptor, calmodulin, and aromatase b by independent samplings of goldfish brains from six seasonal time points and real-time PCR assays.

CONCLUSIONS/SIGNIFICANCE

Using both theoretical and experimental strategies, we report for the first time global gene expression patterns throughout a breeding season which may account for dynamic neuroendocrine regulation of seasonal reproductive development.

Effects of fluoxetine on the reproductive axis of female goldfish (Carassius auratus)

We investigated the effects of fluoxetine, a selective serotonin reuptake inhibitor, on neuroendocrine function and the reproductive axis in female goldfish. Fish were given intraperitoneal injections of fluoxetine twice a week for 14 days, resulting in five injections of 5 μg fluoxetine/g body wt. We measured the monoamine neurotransmitters serotonin, dopamine, and norepinephrine in addition to their metabolites with HPLC. Homovanillic acid, a metabolite in the dopaminergic pathway, increased significantly in the hypothalamus. Plasma estradiol levels were measured by radioimmunoassay and were significantly reduced approximately threefold after fluoxetine treatment. We found that fluoxetine also significantly reduced the expression of estrogen receptor (ER)β1 mRNA by 4-fold in both the hypothalamus and the telencephalon and ERα mRNA by 1.7-fold in the telencephalon. Fluoxetine had no effect on the expression of ERβ2 mRNA in the hypothalamus or telencephalon. Microarray analysis identified isotocin, a neuropeptide that stimulates reproductive behavior in fish, as a candidate gene affected by fluoxetine treatment. Real-time RT-PCR verified that isotocin mRNA was downregulated approximately sixfold in the hypothalamus and fivefold in the telencephalon. Intraperitoneal injection of isotocin (1 μg/g) increased plasma estradiol, providing a potential link between changes in isotocin gene expression and decreased circulating estrogen in fluoxetine-injected fish. Our results reveal targets of serotonergic modulation in the neuroendocrine brain and indicate that fluoxetine has the potential to affect sex hormones and modulate genes involved in reproductive function and behavior in the brain of female goldfish. We discuss these findings in the context of endocrine disruption because fluoxetine has been detected in the environment.

The goldfish (Carassius auratus) as a model for neuroendocrine signaling

Goldfish (Carassius auratus) are excellent model organisms for the neuroendocrine signaling and the regulation of reproduction in vertebrates. Goldfish also serve as useful model organisms in numerous other fields. In contrast to mammals, teleost fish do not have a median eminence; the anterior pituitary is innervated by numerous neuronal cell types and thus, pituitary hormone release is directly regulated. Here we briefly describe the neuroendocrine control of luteinizing hormone. Stimulation by gonadotropin-releasing hormone and a multitude of classical neurotransmitters and neuropeptides is opposed by the potent inhibitory actions of dopamine. The stimulatory actions of gamma-aminobutyric acid and serotonin are also discussed. We will focus on the development of a cDNA microarray composed of carp and goldfish sequences which has allowed us to examine neurotransmitter-regulated gene expression in the neuroendocrine brain and to investigate potential genomic interactions between these key neurotransmitter systems. We observed that isotocin (fish homologue of oxytocin) and activins are regulated by multiple neurotransmitters, which is discussed in light of their roles in reproduction in other species. We have also found that many novel and uncharacterized goldfish expressed sequence tags in the brain are also regulated by neurotransmitters. Their sites of production and whether they play a role in neuroendocrine signaling and control of reproduction remain to be determined. The transcriptomic tools developed to study reproduction could also be used to advance our understanding of neuroendocrine-immune interactions and the relationship between growth and food intake in fish.

Fish and chips: functional genomics of social plasticity in an African cichlid fish

Behavior and physiology are regulated by both environment and social context. A central goal in the study of the social control of behavior is to determine the underlying physiological, cellular and molecular mechanisms in the brain. The African cichlid fish Astatotilapia burtoni has long been used as a model system to study how social interactions regulate neural and behavioral plasticity. In this species, males are either socially dominant and reproductively active or subordinate and reproductively suppressed. This phenotypic difference is reversible. Using an integrative approach that combines quantitative behavioral measurements, functional genomics and bioinformatic analyses, we examine neural gene expression in dominant and subordinate males as well as in brooding females. We confirm the role of numerous candidate genes that are part of neuroendocrine pathways and show that specific co-regulated gene sets (modules), as well as specific functional gene ontology categories, are significantly associated with either dominance or reproductive state. Finally, even though the dominant and subordinate phenotypes are robustly defined, we find a surprisingly high degree of individual variation in the transcript levels of the very genes that are differentially regulated between these phenotypes. The results of the present study demonstrate the molecular complexity in the brain underlying social behavior, identify novel targets for future studies, validate many candidate genes and exploit individual variation in order to gain biological insights.

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.

Using generalized procrustes analysis (GPA) for normalization of cDNA microarray data

BACKGROUND

Normalization is essential in dual-labelled microarray data analysis to remove non-biological variations and systematic biases. Many normalization methods have been used to remove such biases within slides (Global, Lowess) and across slides (Scale, Quantile and VSN). However, all these popular approaches have critical assumptions about data distribution, which is often not valid in practice.

RESULTS

In this study, we propose a novel assumption-free normalization method based on the Generalized Procrustes Analysis (GPA) algorithm. Using experimental and simulated normal microarray data and boutique array data, we systemically evaluate the ability of the GPA method in normalization compared with six other popular normalization methods including Global, Lowess, Scale, Quantile, VSN, and one boutique array-specific housekeeping gene method. The assessment of these methods is based on three different empirical criteria: across-slide variability, the Kolmogorov-Smirnov (K-S) statistic and the mean square error (MSE). Compared with other methods, the GPA method performs effectively and consistently better in reducing across-slide variability and removing systematic bias.

CONCLUSION

The GPA method is an effective normalization approach for microarray data analysis. In particular, it is free from the statistical and biological assumptions inherent in other normalization methods that are often difficult to validate. Therefore, the GPA method has a major advantage in that it can be applied to diverse types of array sets, especially to the boutique array where the majority of genes may be differentially expressed.

Microarray analysis in the zebrafish (Danio rerio) liver and telencephalon after exposure to low concentration of 17alpha-ethinylestradiol

17alpha-ethinylestradiol (EE2) is detected in sewage effluent at concentrations that can disrupt normal reproductive function in fish. The objectives of this study were to identify novel genomic responses to EE2 exposure using microarray and real-time RT-PCR analysis in the liver and telencephalon of male zebrafish. Zebrafish were exposed to an environmentally relevant nominal concentration of 10ng/L EE2 for a 21-day period. In the liver, common biomarkers for estrogenic exposure such as vitellogenin 1 and 3 (vtg1; vtg3), estrogen receptor alpha (esr1), and apolipoprotein A1 (apoA1) mRNA were identified by microarray analysis as being differentially regulated. Real-time RT-PCR confirmed that vtg1 was induced approximately 700-fold, vtg3 was induced approximately 100-fold and esr1 was induced approximately 20-fold. As determined by microarray analysis, ATPase Na+/K+ alpha 1a.4 (atp1a1a.4) and ATPase Na+/K+ beta 1a (atp1b1a) mRNA were down-regulated in the liver. Gene ontology (GO) analysis revealed that there were common biological processes and molecular functions regulated by EE2 in both tissues (e.g. electron transport and cell communication) but there were tissue specific changes in gene categories. For example, genes involved in protein metabolism, carbohydrate metabolism were down-regulated in the liver but were induced in the telencephalon. This study demonstrates that (1) tissues exhibit different gene responses to low EE2 exposure; (2) there are pronounced genomic effects in the liver and (3) multi-tissue gene profiling is needed to improve understanding of the effects of human pharmaceuticals on aquatic organisms.

Effective use of microarrays in neuroendocrine research

The development of microarray technology makes it possible to simultaneously assay the expression level of hundreds to tens of thousands of mRNA transcripts in one experiment. Genome-wide transcriptional analysis has increasing importance for many areas of neuroendocrinology research. The expense and technical complexity of microarray experiments can make it difficult to navigate the terrain of rival platforms and technologies. In this review, we provide a practical view and comparison of various microarray technologies. Affymetrix arrays, high-density cDNA arrays, membrane arrays and experimental design and data analysis are all discussed by researchers currently using these techniques to study gene regulation in neuroendocrine tissues.