Limited sex-biased neural gene expression patterns across strains in Zebrafish (Danio rerio)

Zebrafish gender-specific anxiety-like behavioral and physiological reactions elicited by caffeine

Caffeine exerts a biphasic effect on zebrafish behavior. High doses of caffeine have been associated with increased stress and anxiety, whereas low doses have been found to enhance performance on tasks requiring focus and attention. However, the sex-specific nature of these biphasic effects on behavior and physiology remains unclear. This study assessed the behavioral responses and hormone levels in male and female zebrafish after acute exposure to caffeine ranging from 0.3 to 600 mg/L. The results showed no significant difference in caffeine intake between males and females after acute exposure at each concentration. Caffeine-induced behavioral and physiological responses indicated a threshold dosage existed between 30 and 300 mg/L. Female fish displayed increased anxiety-like behavioral phenotypes, i.e., latency to upper and freezing, whereas males exhibited more erratic movement following acute exposure to a high-dose treatment. In addition, females exhibited a significant increase in whole-body cortisol levels, while males experienced a testosterone elevation at 300 mg/L of caffeine acute exposure. There was a significant decrease in the duration of erratic movements in males treated with the androgen receptor antagonist flutamide compared to the control group. The transcriptome analysis uncovered 511 and 592 up-regulated and 761 and 922 down-regulated differential expression genes in males and females, respectively, compared to the control. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analysis revealed that caffeine has the potential to impact various pathways in zebrafish, including phototransduction and steroid hormone biosynthesis. Our findings demonstrate that testosterone and cortisol play a combined role in regulating stress responses in both behavior and physiology. Furthermore, our study highlights the significance of encompassing both male and female zebrafish as a model system.

Species variation in steroid hormone-related gene expression contributes to species diversity in sexually dimorphic communication in electric fishes

Sexually dimorphic behaviors are often regulated by gonadal steroid hormones. Species diversity in behavioral sex differences may arise as expression of genes mediating steroid action in brain regions controlling these behaviors evolves. The electric communication signals of apteronotid knifefishes are an excellent model for comparatively studying neuroendocrine regulation of sexually dimorphic behavior. These fish produce and detect weak electric organ discharges (EODs) for electrolocation and communication. EOD frequency (EODf), controlled by the medullary pacemaker nucleus (Pn), is sexually dimorphic and regulated by androgens and estrogens in some species, but is sexually monomorphic and unaffected by hormones in other species. We quantified expression of genes for steroid receptors, metabolizing enzymes, and cofactors in the Pn of two species with sexually dimorphic EODf (Apteronotus albifrons and Apteronotus leptorhynchus) and two species with sexually monomorphic EODf ("Apteronotus" bonapartii and Parapteronotus hasemani). The "A." bonapartii Pn expressed lower levels of androgen receptor (AR) genes than the Pn of species with sexually dimorphic EODf. In contrast, the P. hasemani Pn robustly expressed AR genes, but expressed lower levels of genes for 5α-reductases, which convert androgens to more potent metabolites, and higher levels of genes for 17β-hydroxysteroid dehydrogenases that oxidize androgens and estrogens to less potent forms. These findings suggest that sexual monomorphism of EODf arose convergently via two different mechanisms. In "A." bonapartii, reduced Pn expression of ARs likely results in insensitivity of EODf to androgens, whereas in P. hasemani, gonadal steroids may be metabolically inactivated in the Pn, reducing their potential to influence EODf.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

First neurotranscriptome of adults Tambaquis (Colossoma macropomum) with characterization and differential expression between males and females

The Tambaqui is one of the most representative Amazon fish species, being highly exploited in fisheries, aquaculture and as a research model. Nonetheless, data about functional genome are still required to evaluate reproductive and nutrition parameters as well as resistance to pathogens. The of next-generation sequencing has allows assessing the transcriptional processes in non-model species by providing comprehensive gene collections to be used as a database in further genomic applications and increased performance of captive populations. In this study, we relied on RNAseq approach to generate the first transcriptome of the telencephalon from adult males and females of Colossoma macropomum, resulting in a reference dataset for future functional studies. We retrieved 896,238 transcripts, including the identification of 267,785 contigs and 203,790 genes. From this total, 91 transcripts were differentially expressed, being 63 and 28 of them positively regulated for females and males, respectively. The functional annotation resulted in a library of 40 candidate genes for females and 20 for males. The functional enrichment classes comprised reproductive processes (GO:0,048,609; GO:0,003,006; GO:0,044,703; GO:0,032,504; GO:0,019,953) being related to sex differentiation (e.g., SAFB) and immune response (e.g., SLC2A6, AHNAK, NLRC3, NLRP3 and IgC MHC I alpha3), thus indicating that the genes in the neurotranscriptome of Tambaqui participate in sex differentiation and homeostasis of captive specimens. These data are useful to design the selection of genes related to sex determination and animal welfare in raising systems of Tambaqui.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting the hypothesis that the Chr4R transcriptome might be different from the rest of the genome. We conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes and identified four regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brain and liver validated few transcripts from Region-2 in somatic cells, but without sex-specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. In Region-2, protein-coding genes lack human orthologs; it has zinc finger genes expressed early in zygotic genome activation; it has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and an distinct set of repetitive elements. The colocalization of 1) genes silenced in ovaries but not in testes that are 2) expressed in embryos briefly at the onset of zygotic genome activation; 3) maternal-specific genes for translation machinery; 4) maternal-specific spliceosome components; and 4) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a Maternal-to-Zygotic-Transition Gene Regulatory Block.

Transcriptomic signatures associated with underlying rapid changes in the early phase brain of bi-directional sex change in Trimma okinawae

Teleost fish exhibit remarkable sexual plasticity and divergent developmental systems, including sequential hermaphroditism. One of the more fascinating models of sexual plasticity is socially controlled sex change, which is often observed in coral reef fish. The Okinawa rubble goby, Trimma okinawae, is a bi-directional sex-changing fish. It can rapidly change sex in either direction based on social circumstances. Although behavioural and neuroendocrine sex change occurs immediately and is believed to trigger gonadal changes, the underlying mechanisms remain poorly understood. In this study, we conducted a de novo transcriptome analysis of the T. okinawae brain and identified genes that are differentially expressed between the sexes and genes that were immediately controlled by social stimulation implicating sex change. Several genes showed concordant expression shifts regardless of the sex change direction and were associated with histone modification in nerve cells. These genes are known to function in the neuroendocrine control of reproduction in nerve cells. Overall, we identified genes associated with the initiation of sex change, which provides insight into the regulation of sex change and sexual plasticity.

Sex and gametogenesis stage are strong drivers of gene expression in Mytilus edulis exposed to environmentally relevant plasticiser levels and pH 7.7

Plastic pollution and changes in oceanic pH are both pressing environmental issues. Little emphasis, however, has been placed on the influence of sex and gametogenesis stage when investigating the effects of such stressors. Here, we examined histology and molecular biomarkers of blue mussels Mytilus edulis exposed for 7 days to a pH 7.7 scenario (- 0.4 units) in combination with environmentally relevant concentrations (0, 0.5 and 50 µg/L) of the endocrine disrupting plasticiser di-2-ethylhexyl phthalate (DEHP). Through a factorial design, we investigated the gametogenesis cycle and sex-related expression of genes involved in pH homeostasis, stress response and oestrogen receptor-like pathways after the exposure to the two environmental stressors. As expected, we found sex-related differences in the proportion of developing, mature and spawning gonads in histological sections. Male gonads also showed higher levels of the acid-base regulator CA2, but females had a higher expression of stress response-related genes (i.e. sod, cat, hsp70). We found a significant effect of DEHP on stress response-related gene expression that was dependent on the gametogenesis stage, but there was only a trend towards downregulation of CA2 in response to pH 7.7. In addition, differences in gene expression between males and females were most pronounced in experimental conditions containing DEHP and/or acidified pH but never the control, indicating that it is important to consider sex and gametogenesis stage when studying the response of mussels to diverse stressors.

The critical impact of sex on preclinical alcohol research - Insights from zebrafish

Sex is an important biological variable that is widely recognized in studies of alcohol-related effects. Complementing clinical and preclinical rodent research, the zebrafish (Danio rerio) is the second most used laboratory species, and a powerful model organism in biomedicine. Like clinical and rodent models, zebrafish demonstrate overt sex differences in alcohol-related responses. Collectively, this evidence shows that the zebrafish becomes a sensitive model species to further probe in-depth sex differences commonly reported in alcohol research.

Sex-specific differences in zebrafish brains

In this systematic review, we highlight the differences between the male and female zebrafish brains to understand their differentiation and their use in studying sex-specific neurological diseases. Male and female brains display subtle differences at the cellular level which may be important in driving sex-specific signaling. Sex differences in the brain have been observed in humans as well as in non-human species. However, the molecular mechanisms of brain sex differentiation remain unclear. The classical model of brain sex differentiation suggests that the steroid hormones derived from the gonads are the primary determinants in establishing male and female neural networks. Recent studies indicate that the developing brain shows sex-specific differences in gene expression prior to gonadal hormone action. Hence, genetic differences may also be responsible for differentiating the brain into male and female types. Understanding the signaling mechanisms involved in brain sex differentiation could help further elucidate the sex-specific incidences of certain neurological diseases. The zebrafish model could be appropriate for enhancing our understanding of brain sex differentiation and the signaling involved in neurological diseases. Zebrafish brains show sex-specific differences at the hormonal level, and recent advances in RNA sequencing have highlighted critical sex-specific differences at the transcript level. The differences are also evident at the cellular and metabolite levels, which could be important in organizing sex-specific neuronal signaling. Furthermore, in addition to having one ortholog for 70% of the human gene, zebrafish also shares brain structural similarities with other higher eukaryotes, including mammals. Hence, deciphering brain sex differentiation in zebrafish will help further enhance the diagnostic and pharmacological intervention of neurological diseases.

A Glitch in the Matrix: The Role of Extracellular Matrix Remodeling in Opioid Use Disorder

Opioid use disorder (OUD) and deaths from drug overdoses have reached unprecedented levels. Given the enormous impact of the opioid crisis on public health, a more thorough, in-depth understanding of the consequences of opioids on the brain is required to develop novel interventions and pharmacological therapeutics. In the brain, the effects of opioids are far reaching, from genes to cells, synapses, circuits, and ultimately behavior. Accumulating evidence implicates a primary role for the extracellular matrix (ECM) in opioid-induced plasticity of synapses and circuits, and the development of dependence and addiction to opioids. As a network of proteins and polysaccharides, including cell adhesion molecules, proteases, and perineuronal nets, the ECM is intimately involved in both the formation and structural support of synapses. In the human brain, recent findings support an association between altered ECM signaling and OUD, particularly within the cortical and striatal circuits involved in cognition, reward, and craving. Furthermore, the ECM signaling proteins, including matrix metalloproteinases and proteoglycans, are directly involved in opioid seeking, craving, and relapse behaviors in rodent opioid models. Both the impact of opioids on the ECM and the role of ECM signaling proteins in opioid use disorder, may, in part, depend on biological sex. Here, we highlight the current evidence supporting sex-specific roles for ECM signaling proteins in the brain and their associations with OUD. We emphasize knowledge gaps and future directions to further investigate the potential of the ECM as a therapeutic target for the treatment of OUD.

Behavioral genetics of alcohol's effects in three zebrafish (Danio rerio) populations

Alcohol abuse is one of the most dangerous and serious problems for patients and society. Interpopulation studies are important in understanding how genetic background contributes to the effects of alcohol. In this study, we applied a chronic alcohol exposure protocol in three zebrafish populations (Danio rerio; both sexes; AB, TU, and outbred fish - OB). We analyzed the behavioral responses and mRNA expression involved in neurotransmitter metabolism - th1, tph1, ache, ada1, gaba1, gad1b, and bdnf. Locomotion patterns were similar between populations (increased speed after acute alcohol and unaltered locomotion after chronic and withdrawal treatments). All populations exhibited increased expression of genes associated with locomotion (th1, gad1b, and gaba1) after acute alcohol exposure. Anxiety-like responses increased in AB and TU fish during withdrawal and decreased in AB fish after acute alcohol exposure. Genes related to anxiety-like behavior (tph1 and ada1) were overexpressed in AB and TU fish after acute and withdrawal treatments, while OB fish exhibited unaltered responses. Bdnf levels decreased during withdrawal in AB and OB fish, while TU showed upregulated levels in both chronic and withdrawal treatments. Our results suggest that zebrafish populations respond differently to alcohol exposure, which may contribute to understanding the mechanisms underlying alcohol use and dependence. Moreover, we found that a more diverse genetic background (OB) was related to higher variability in behavioral and mRNA expression, demonstrating that inbred populations (AB and TU) may be useful tools in identifying alcohol use and abuse mechanisms.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

Npas4a expression in the teleost forebrain is associated with stress coping style differences in fear learning

Learning to anticipate potentially dangerous contexts is an adaptive behavioral response to coping with stressors. An animal's stress coping style (e.g. proactive-reactive axis) is known to influence how it encodes salient events. However, the neural and molecular mechanisms underlying these stress coping style differences in learning are unknown. Further, while a number of neuroplasticity-related genes have been associated with alternative stress coping styles, it is unclear if these genes may bias the development of conditioned behavioral responses to stressful stimuli, and if so, which brain regions are involved. Here, we trained adult zebrafish to associate a naturally aversive olfactory cue with a given context. Next, we investigated if expression of two neural plasticity and neurotransmission-related genes (npas4a and gabbr1a) were associated with the contextual fear conditioning differences between proactive and reactive stress coping styles. Reactive zebrafish developed a stronger conditioned fear response and showed significantly higher npas4a expression in the medial and lateral zones of the dorsal telencephalon (Dm, Dl), and the supracommissural nucleus of the ventral telencephalon (Vs). Our findings suggest that the expression of activity-dependent genes like npas4a may be differentially expressed across several interconnected forebrain regions in response to fearful stimuli and promote biases in fear learning among different stress coping styles.

Early social deprivation does not affect cortisol response to acute and chronic stress in zebrafish

Social isolation is a well-established technique for inducing early adversity but, in rodent models, the need of parental care makes it difficult to distinguish the effects of social deprivation from the consequences of nutritional deficiencies. Zebrafish do not require parental care, allowing separation of social deprivation from nutritional deprivation, and have emerged as a promising model to study ontogeny of normal and pathological behaviors relevant for human neuropsychological disorders. Previous reports of life-long isolation in zebrafish showed some consistency with mammalian literature, depicting later social deficits and locomotor hyperactivity. However, unlike reports of higher anxiety and stress behavior in isolated rodents and primates, behavioral responses were tapered in isolated fish. To examine whether life-long developmental isolation has a dampening effect on zebrafish endocrine stress response, we applied stressors to zebrafish siblings that were either isolated or socialized, and compared their whole-body cortisol levels with non-stressed control siblings kept in low-housing densities. Utilizing previously validated paradigms (exposure to novel tank and unpredictable chronic stress), we exposed separate groups (n = 9-14, mixed-sex) of social and isolated zebrafish to acute and chronic stressors and measured their cortisol levels. A univariate ANOVA and post-hoc Tukey's HSD tests confirmed that compared to socially raised control fish, developmental isolation did not increase baseline cortisol levels in zebrafish. Additionally, compared to the non-stressed condition, application of both acute and chronic stressors significantly increased cortisol levels in isolated fish and, to a similar degree, to socially raised fish. Our findings suggest that zebrafish isolation studies may help separate effects of social deprivation from nonsocial aspects of early adversity. These studies further substantiate the use of developmental isolation in zebrafish, particularly with acute and chronic stress paradigms, for modeling neuropsychological disorders.LAY SUMMARYA difficult childhood can make humans react more frequently or severely to later stress and modeling this effect in animals can help explain how and why early stress affects subsequent mental and physical health. Early social isolation does not affect later response to stressful situations in adult zebrafish, providing us with a model of psychiatric disorders that allows separation of effects of poor physical environments (lacking food, shelter, etc.) from poor social environments (lack of appropriate socialization).

Analysis across diverse fish species highlights no conserved transcriptome signature for proactive behaviour

Background

Consistent individual differences in behaviour, known as animal personalities, have been demonstrated within and across species. In fish, studies applying an animal personality approach have been used to resolve variation in physiological and molecular data suggesting a linkage, genotype-phenotype, between behaviour and transcriptome regulation. In this study, using three fish species (zebrafish; Danio rerio, Atlantic salmon; Salmo salar and European sea bass; Dicentrarchus labrax), we firstly address whether personality-specific mRNA transcript abundances are transferrable across distantly-related fish species and secondly whether a proactive transcriptome signature is conserved across all three species.

Results

Previous zebrafish transcriptome data was used as a foundation to produce a curated list of mRNA transcripts related to animal personality across all three species. mRNA transcript copy numbers for selected gene targets show that differential mRNA transcript abundance in the brain appears to be partially conserved across species relative to personality type. Secondly, we performed RNA-Seq using whole brains from S. salar and D. labrax scoring positively for both behavioural and molecular assays for proactive behaviour. We further enriched this dataset by incorporating a zebrafish brain transcriptome dataset specific to the proactive phenotype. Our results indicate that cross-species molecular signatures related to proactive behaviour are functionally conserved where shared functional pathways suggest that evolutionary convergence may be more important than individual mRNAs.

Conclusions

Our data supports the proposition that highly polygenic clusters of genes, with small additive effects, likely support the underpinning molecular variation related to the animal personalities in the fish used in this study. The polygenic nature of the proactive brain transcriptome across all three species questions the existence of specific molecular signatures for proactive behaviour, at least at the granularity of specific regulatory gene modules, level of genes, gene networks and molecular functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07317-z.

Comparative transcriptomics implicate mitochondrial and neurodevelopmental impairments in larval zebrafish (Danio rerio) exposed to two selective serotonin reuptake inhibitors (SSRIs)

Pharmaceuticals and personal care products are emerging contaminants that are increasingly detected in the environment worldwide. Certain classes of pharmaceuticals, such as selective serotonin reuptake inhibitors (SSRIs), are a major environmental concern due to their widespread use and the fact that these compounds are designed to have biological effects at low doses. A complication in predicting toxic effects of SSRIs in nontarget organisms is that their mechanism of action is not fully understood. To better understand the potential toxic effects of SSRIs, we employed an ultra-low input RNA-sequencing method to identify potential pathways that are affected by early exposure to two SSRIs (fluoxetine and paroxetine). We exposed wildtype zebrafish (Danio rerio) embryos to 100 μg/L of either fluoxetine or paroxetine for 6 days before extracting and sequencing mRNA from individual larval brains. Differential gene expression analysis identified 1550 genes that were significantly affected by SSRI exposure with a core set of 138 genes altered by both SSRIs. Weighted gene co-expression network analysis identified 7 modules of genes whose expression patterns were significantly correlated with SSRI exposure. Functional enrichment analysis of differentially expressed genes as well as network module genes repeatedly identified various terms associated with mitochondrial and neuronal structures, mitochondrial respiration, and neurodevelopmental processes. The enrichment of these terms indicates that toxic effects of SSRI exposure are likely caused by mitochondrial dysfunction and subsequent neurodevelopmental effects. To our knowledge, this is the first effort to study the tissue-specific transcriptomic effects of SSRIs in developing zebrafish, providing specific, high resolution molecular data regarding the sublethal effects of SSRI exposure.

Age Bias in Zebrafish Models of Epilepsy: What Can We Learn From Old Fish?

Zebrafish are a powerful tool for investigating epilepsy. Mammalian seizures can be recapitulated molecularly, behaviorally, and electrophysiologically, using a fraction of the resources required for experiments in mammals. Larval zebrafish offer exceptionally economical and high-throughput approaches and are amenable to state-of-the-art genetic engineering techniques, providing valuable transgenic models of human diseases. For these reasons, larvae tend to be chosen for studying epilepsy, but the value of adult zebrafish may be underappreciated. Zebrafish exhibit transient larval - adult duality. The incompletely developed neural system of larval zebrafish may limit the translation of complex neurological disorders. Larval zebrafish go through dynamic changes during ontogenesis, whereas adult zebrafish are physiologically more stable. Adult zebrafish have a full range of complex brain structures and functions, such as an endothelial blood-brain barrier and adult neurogenesis, both are significant factors in epilepsy research. This review highlights the differences between larval and adult zebrafish that should be considered in pathophysiological and pharmacological studies of epilepsy.

Differential effects of ethanol on behavior and GABAA receptor expression in adult zebrafish (Danio rerio) with alternative stress coping styles

Variation in stress responses between individuals are linked to factors ranging from stress coping styles to sensitivity of neurotransmitter systems. Many anxiolytic compounds (e.g. ethanol) can increase stressor engagement through modulation of neurotransmitter systems and are used to investigate stress response mechanisms. There are two alternative suites of correlated behavioral and physiological responses to stressors (stress coping styles) that differ in exploration tendencies: proactive and reactive stress coping styles. By chronically treating individuals differing in stress coping style with ethanol, a GABA-acting drug, we assessed the role of the GABAergic system on the behavioral stress response. Specifically, we investigated resulting changes in stress-related behavior (i.e. exploratory behavior) and whole-brain GABAA receptor subunits (gabra1, gabra2, gabrd, & gabrg2) in response to a novelty stressor. We found that ethanol-treated proactive individuals showed lower stress-related behaviors than their reactive counterparts. Proactive individuals showed significantly higher expression of gabra1, gabra2, and gabrg2 compared to reactive individuals and ethanol treatment resulted in upregulation of gabra1 and gabrg2 in both stress coping styles. These results suggest that impacts of ethanol on stress-related behaviors vary by stress coping style and that expression of select GABAA receptor subunits may be one of the underlying mechanisms.

Environmental Enrichment Improved Learning and Memory, Increased Telencephalic Cell Proliferation, and Induced Differential Gene Expression in Colossoma macropomum

Fish use spatial cognition based on allocentric cues to navigate, but little is known about how environmental enrichment (EE) affects learning and memory in correlation with hematological changes or gene expression in the fish brain. Here we investigated these questions in Colossoma macropomum (Teleostei). Fish were housed for 192 days in either EE or in an impoverished environment (IE) aquarium. EE contained toys, natural plants, and a 12-h/day water stream for voluntary exercise, whereas IE had no toys, plants, or water stream. A third plus maze aquarium was used for spatial and object recognition tests. Compared with IE, the EE fish showed greater learning rates, body length, and body weight. After behavioral tests, whole brain tissue was taken, stored in RNA-later, and then homogenized for DNA sequencing after conversion of isolated RNA. To compare read mapping and gene expression profiles across libraries for neurotranscriptome differential expression, we mapped back RNA-seq reads to the C. macropomum de novo assembled transcriptome. The results showed significant differential behavior, cell counts and gene expression in EE and IE individuals. As compared with IE, we found a greater number of cells in the telencephalon of individuals maintained in EE but no significant difference in the tectum opticum, suggesting differential plasticity in these areas. A total of 107,669 transcripts were found that ultimately yielded 64 differentially expressed transcripts between IE and EE brains. Another group of adult fish growing in aquaculture conditions were either subjected to exercise using running water flow or maintained sedentary. Flow cytometry analysis of peripheral blood showed a significantly higher density of lymphocytes, and platelets but no significant differences in erythrocytes and granulocytes. Thus, under the influence of contrasting environments, our findings showed differential changes at the behavioral, cellular, and molecular levels. We propose that the differential expression of selected transcripts, number of telencephalic cell counts, learning and memory performance, and selective hematological cell changes may be part of Teleostei adaptive physiological responses triggered by EE visuospatial and somatomotor stimulation. Our findings suggest abundant differential gene expression changes depending on environment and provide a basis for exploring gene regulation mechanisms under EE in C. macropomum.

Large-scale analysis of zebrafish (Danio rerio) transcriptomes identifies functional modules associated with phenotypes

Zebrafish (Danio rerio) is an excellent model for biomedicine research due to its genetic accessibility and optical transparency. A large number of microarray based transcriptomes of zebrafish have been profiled in various cell types, tissues, development stages, toxicological exposures and other conditions. However, there is still no easy-to-use web tool to explore those precious data. We downloaded 1434 microarray data from National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO), constructed weighted gene co-expression network, and identified 50 modules of co-expressed genes that correspond to different cell types, tissues, development stages, and other experimental conditions. These modules were associated with experiments/traits, and may serve signature modules for phenotypes. Hub genes were screened by intra-modular connectivity. Higher-order module networks analysis suggested that nucleus and cell cycle modules are densely connected. Module-based gene function identification may help to discover novel gene function. Our web tool provides a new resource for gene function study in zebrafish (http://bioinformatics.fafu.edu.cn/zebrafish/).

Transcriptome of tambaqui Colossoma macropomum during gonad differentiation: Different molecular signals leading to sex identity

Tambaqui (Colossoma macropomum) is the major native species in Brazilian aquaculture, and we have shown that females exhibit a higher growth compared to males, opening up the possibility for the production of all-female population. To date, there is no information on the sex determination and differentiation molecular mechanisms of tambaqui. In the present study, transcriptome sequencing of juvenile trunks was performed to understand the molecular network involved in the gonadal sex differentiation. The results showed that before differentiation, components of the Wnt/β-catenin pathway, fox and fst genes imprint female sex development, whereas antagonistic pathways (gsk3b, wt1 and fgfr2), sox9 and genes for androgen synthesis indicate male differentiation. Hence, in undifferentiated tambaqui, the Wnt/β-catenin exerts a role on sex differentiation, either upregulated in female-like individuals, or antagonized in male-like individuals.

Transcriptome profiling analysis of sex-based differentially expressed mRNAs and lncRNAs in the brains of mature zebrafish (Danio rerio)

Background

Similar to humans, the zebrafish brain plays a central role in regulating sexual reproduction, maturation and sexual behavior. However, systematic studies of the dimorphic patterns of gene expression in the brain of male and female zebrafish are lacking.

Results

In this study, the mRNA and lncRNA expression profiles were obtained from the brain tissue samples of the three male and three female zebrafish by high-throughput transcriptome sequencing. We identified a total of 108 mRNAs and 50 lncRNAs with sex-based differential expression. We randomly selected four differentially expressed genes for RT-qPCR verification and the results certified that the expression pattern showed a similar trend between RNA-seq and RT-qPCR results. Protein-protein interaction network analysis, Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to obtain the biological significance of differentially expressed mRNA in the brain dimorphism of zebrafish. Finally, a Pearson correlation analysis was performed to construct the co-expression network of the mRNAs and lncRNAs.

Conclusions

We found that 12 new lncRNAs not only have significant gender specificity in the brain of zebrafish, and this finding may provide a clue to further study of the functional difference between male and female zebrafish brain.

Identification of sex differences in zebrafish (Danio rerio) brains during early sexual differentiation and masculinization using 17α-methyltestoterone

Sexual behavior in teleost fish is highly plastic. It can be attributed to the relatively few sex differences found in adult brain transcriptomes. Environmental and hormonal factors can influence sex-specific behavior. Androgen treatment stimulates behavioral masculinization. Sex dimorphic gene expression in developing teleost brains and the molecular basis for androgen-induced behavioral masculinization are poorly understood. In this study, juvenile zebrafish (Danio rerio) were treated with 100 ng/L of 17 alpha-methyltestosterone (MT) during sexual development from 20 days post fertilization to 40 days and 60 days post fertilization. We compared brain gene expression patterns in MT-treated zebrafish with control males and females using RNA-Seq to shed light on the dynamic changes in brain gene expression during sexual development and how androgens affect brain gene expression leading to behavior masculinization. We found modest differences in gene expression between juvenile male and female zebrafish brains. Brain aromatase (cyp19a1b), prostaglandin 3a synthase (ptges3a), and prostaglandin reductase 1 (ptgr1) were among the genes with sexually dimorphic expression patterns. MT treatment significantly altered gene expression relative to both male and female brains. Fewer differences were found among MT-treated brains and male brains compared to female brains, particularly at 60 dpf. MT treatment upregulated the expression of hydroxysteroid 11-beta dehydrogenase 2 (hsd11b2), deiodinase, iodothyronine, type II (dio2), and gonadotrophin releasing hormones (GnRH) 2 and 3 (gnrh2 and gnrh3) suggesting local synthesis of 11-ketotestosterone, triiodothyronine, and GnRHs in zebrafish brains which are influenced by androgens. Androgen, estrogen, prostaglandin, thyroid hormone, and GnRH signaling pathways likely interact to modulate teleost sexual behavior.

Effects of fasting on the central expression of appetite-regulating and reproductive hormones in wild-type and Casper zebrafish (Danio rerio)

Appetite and reproduction are closely related functions that are both regulated by brain hormones. Appetite stimulators include orexin and neuropeptide Y (NPY), and reproductive hormones include gonadotropin-releasing hormone (GnRH), gonadotropin-inhibitory hormone (GnIH), kisspeptin, and neurokinin B (NKB). GnRH stimulates the secretion of pituitary gonadotropes, and kisspeptin and GnIH modulate this action. Kisspeptin secretion is further controlled by neurokinin B (NKB) and dynorphin A (Dyn). To better understand the mechanisms regulating appetite and reproduction in fish, we examined the effects of fasting, reproductive stage, gender, and strain on the brain mRNA expression of appetite (orexin and NPY) and reproductive (GnRH, kisspeptin, GnIH, and NKB) hormones in zebrafish. In order to compare strains, we used both wild-type and transparent Casper zebrafish. In female wild-type zebrafish, fasting increased the expression of all hormones investigated, with the exception of Kiss2. Only NPY and Kiss2 were increased in male wild-type zebrafish during fasting. In Casper zebrafish, only GnIH and NKB in males were affected by fasting, suggesting that Casper fish may be more resistant to fasting than wild fish. Fasting increased expressions of orexin, GnRH2, Kiss1, GnIH and NKB in wild-type females with more eggs or larger eggs relative to body weight, compared to those with fewer or smaller eggs, suggesting that more mature females are more affected by fasting. No significant interactions of fasting and reproductive stage were noted in female Casper fish. To investigate whether differences between Casper and wild-type fish were due to genes involved in pigmentation, we compared the brain mRNA expressions of enzymes involved in melanin synthesis (tyrosinase and tyrosine hydroxylase - TH), melanocortin receptors (MC3R and MC4R), and the melanocortin precursor (proopiomelanocortin - POMC) between the two strains. Casper zebrafish had lower levels of MC3R, tyrosinase, TH1, TH2, and POMC than wild-type fish. Overall, our results suggest the existence of gender- and reproductive stage-specific, as well as strain-specific variations in the mechanisms regulating feeding and reproduction in zebrafish, and that the melanocortin system and melanin pathways may be in part responsible for these differences between strains.

The Genetics and Epigenetics of Sex Change in Fish

Fish show extraordinary sexual plasticity, changing sex naturally as part of their life cycle or reversing sex because of environmental stressors. This plasticity shows that sexual fate is not an irreversible process but the result of an ongoing tug-of-war for supremacy between male and female signaling networks. The behavioral, gonadal, and morphological changes involved in this process are well described, yet the molecular events that underpin those changes remain poorly understood. Epigenetic modifications emerge as a critical link between environmental stimuli, the onset of sex change, and subsequent maintenance of sexual phenotype. Here we synthesize current knowledge of sex change, focusing on the genetic and epigenetic processes that are likely involved in the initiation and regulation of sex change. We anticipate that better understanding of sex change in fish will shed new light on sex determination and development in vertebrates and on how environmental perturbations affect sexual fate.

Contextual fear learning and memory differ between stress coping styles in zebrafish

Animals frequently overcome stressors and the ability to learn and recall these salient experiences is essential to an individual’s survival. As part of an animal’s stress coping style, behavioral and physiological responses to stressors are often consistent across contexts and time. However, we are only beginning to understand how cognitive traits can be biased by different coping styles. Here we investigate learning and memory differences in zebrafish (Danio rerio) displaying proactive and reactive stress coping styles. We assessed learning rate and memory duration using an associative fear conditioning paradigm that trained zebrafish to associate a context with exposure to a natural olfactory alarm cue. Our results show that both proactive and reactive zebrafish learn and remember this fearful association. However, we note significant interaction effects between stress coping style and cognition. Zebrafish with the reactive stress coping style acquired the fear memory at a significantly faster rate than proactive fish. While both stress coping styles showed equal memory recall one day post-conditioning, reactive zebrafish showed significantly stronger recall of the conditioned context relative to proactive fish four days post-conditioning. Through understanding how stress coping strategies promote biases in processing salient information, we gain insight into mechanisms that can constrain adaptive behavioral responses.

Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change

Bluehead wrasses undergo dramatic, socially cued female-to-male sex change. We apply transcriptomic and methylome approaches in this wild coral reef fish to identify the primary trigger and subsequent molecular cascade of gonadal metamorphosis. Our data suggest that the environmental stimulus is exerted via the stress axis and that repression of the aromatase gene (encoding the enzyme converting androgens to estrogens) triggers a cascaded collapse of feminizing gene expression and identifies notable sex-specific gene neofunctionalization. Furthermore, sex change involves distinct epigenetic reprogramming and an intermediate state with altered epigenetic machinery expression akin to the early developmental cells of mammals. These findings reveal at a molecular level how a normally committed developmental process remains plastic and is reversed to completely alter organ structures.

Sex differences in behavior and neuropharmacology of zebrafish

Sex is an important variable in biomedical research. The zebrafish (Danio rerio) is increasingly utilized as a powerful new model organism in translational neuroscience and pharmacology. Mounting evidence indicates important sex differences in zebrafish behavioral and neuropharmacological responses. Here, we discuss the role of sex in zebrafish central nervous system (CNS) models, their molecular mechanisms, recent findings and the existing challenges in this field. We also emphasize the growing utility of zebrafish models in translational neuropharmacological research of sex differences, fostering future CNS drug discovery and the search for novel sex-specific therapies. Finally, we highlight the interplay between sex and environment in zebrafish models of sex-environment correlations as an important strategy of CNS disease modeling using this aquatic organism.

Differences in stress reactivity between zebrafish with alternative stress coping styles

Animals experience stress in a variety of contexts and the behavioural and neuroendocrine responses to stress can vary among conspecifics. The responses across stressors often covary within an individual and are consistently different between individuals, which represent distinct stress coping styles (e.g. proactive and reactive). While studies have identified differences in peak glucocorticoid levels, less is known about how cortisol levels differ between stress coping styles at other time points of the glucocorticoid stress response. Here we quantified whole-body cortisol levels and stress-related behaviours (e.g. depth preference, movement) at time points representing the rise and recovery periods of the stress response in zebrafish lines selectively bred to display the proactive and reactive coping style. We found that cortisol levels and stress behaviours are significantly different between the lines, sexes and time points. Further, individuals from the reactive line showed significantly higher cortisol levels during the rising phase of the stress response compared with those from the proactive line. We also observed a significant correlation between individual variation of cortisol levels and depth preference but only in the reactive line. Our results show that differences in cortisol levels between the alternative stress coping styles extend to the rising phase of the endocrine stress response and that cortisol levels may explain variation in depth preferences in the reactive line. Differences in the timing and duration of cortisol levels may influence immediate behavioural displays and longer lasting neuromolecular mechanisms that modulate future responses.

DNA Methylation Patterns in the Round Goby Hypothalamus Support an On-The-Spot Decision Scenario for Territorial Behavior

The question as to how early life experiences are stored on a molecular level and affect traits later in life is highly topical in ecology, medicine, and epigenetics. In this study, we use a fish model to investigate whether DNA methylation mediates early life experiences and predetermines a territorial male reproductive phenotype. In fish, adult reproductive phenotypes frequently depend on previous life experiences and are often associated with distinct morphological traits. DNA methylation is an epigenetic mechanism which is both sensitive to environmental conditions and stably inherited across cell divisions. We therefore investigate early life predisposition in the round goby Neogobius melanostomus by growth back-calculations and then study DNA methylation by MBD-Seq in the brain region controlling vertebrate reproductive behavior, the hypothalamus. We find a link between the territorial reproductive phenotype and high growth rates in the first year of life. However, hypothalamic DNA methylation patterns reflect the current behavioral status independently of early life experiences. Together, our data suggest a non-predetermination scenario in the round goby, in which indeterminate males progress to a non-territorial status in the spawning season, and in which some males then assume a specialized territorial phenotype if current conditions are favorable.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes (e.g., translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr.

Telencephalon transcriptome analysis of chronically stressed adult zebrafish

Chronic stress leads to disruptions in learning and memory processes. The effects of chronic stress experience on the adult zebrafish brain, particularly the memory associated telencephalon brain region, is unclear. The goal of this study was to identify gene expression changes in the adult zebrafish brain triggered by chronic unpredictable stress. Transcriptome analysis of the telencephalon revealed 155 differentially expressed genes. Of these genes, some are critical genes involved in learning and memory, such as cdk5 and chrna7, indicating effects of chronic unpredictable stress on zebrafish memory. Interestingly, several genes were annotated in the Orange domain, which is an amino acid sequence present in eukaryotic DNA-binding transcription repressors. Furthermore, we identified hsd11b2, a cortisol inactivating gene, as chronic stress-responsive in the whole zebrafish brain. Collectively, these findings suggest that memory associated gene expression changes in adult zebrafish telencephalon are affected by chronic stress experience.

Genomic analysis of Sparus aurata reveals the evolutionary dynamics of sex-biased genes in a sequential hermaphrodite fish

Sexual dimorphism is a fascinating subject in evolutionary biology and mostly results from sex-biased expression of genes, which have been shown to evolve faster in gonochoristic species. We report here genome and sex-specific transcriptome sequencing of Sparus aurata, a sequential hermaphrodite fish. Evolutionary comparative analysis reveals that sex-biased genes in S. aurata are similar in number and function, but evolved following strikingly divergent patterns compared with gonochoristic species, showing overall slower rates because of stronger functional constraints. Fast evolution is observed only for highly ovary-biased genes due to female-specific patterns of selection that are related to the peculiar reproduction mode of S. aurata, first maturing as male, then as female. To our knowledge, these findings represent the first genome-wide analysis on sex-biased loci in a hermaphrodite vertebrate species, demonstrating how having two sexes in the same individual profoundly affects the fate of a large set of evolutionarily relevant genes.

Repeatability and reliability of exploratory behavior in proactive and reactive zebrafish, Danio rerio

Behavioral responses to novel situations often vary and can belong to a suite of correlated behaviors. Characteristic behaviors of different personality types (e.g. stress coping styles) are generally consistent across contexts and time. Here, we compare the repeatability and reliability of exploratory behaviors between zebrafish strains selectively bred to display contrasting behavioral responses to stressors that represent the proactive-reactive axis. Specifically, we measure exploratory behavior of individual fish in an open field test over five weeks. We quantified the stationary time, average swimming speed and time spent by a fish in the center area. We found a number of strain differences for each behavioral measure. Stationary time was the most repeatable and reliable measure for assessing proactive-reactive behavioral differences. Reactive zebrafish generally showed the highest reliability and repeatability of exploratory behavior compared to proactive zebrafish and a separate wild caught strain. Given the increased interest in the evolutionary consequences and proximate mechanisms of consistent individual differences, it will be important to continue to investigate how different selective pressures may influence expression of stress coping styles and their effects on the consistency of an animal’s behavior.

Gender-Specific Effects of Selection for Drinking in the Dark on the Network Roles of Coding and Noncoding RNAs

BACKGROUND

Transcriptional differences between heterogeneous stock mice and high drinking-in-the-dark selected mouse lines have previously been described based on microarray technology coupled with network-based analysis. The network changes were reproducible in 2 independent selections and largely confined to 2 distinct network modules; in contrast, differential expression appeared more specific to each selected line. This study extends these results by utilizing RNA-Seq technology, allowing evaluation of the relationship between genetic risk and transcription of noncoding RNA (ncRNA); we additionally evaluate sex-specific transcriptional effects of selection.

METHODS

Naïve mice (N = 24/group and sex) were utilized for gene expression analysis in the ventral striatum; the transcriptome was sequenced with the Illumina HiSeq platform. Differential gene expression and the weighted gene co-expression network analysis were implemented largely as described elsewhere, resulting in the identification of genes that change expression level or (co)variance structure.

RESULTS

Across both sexes, we detect selection effects on the extracellular matrix and synaptic signaling, although the identity of individual genes varies. A majority of nc RNAs cluster in a single module of relatively low density in both the male and female network. The most strongly differentially expressed transcript in both sexes was Gm22513, a small nuclear RNA with unknown function. Associated with selection, we also found a number of network hubs that change edge strength and connectivity. At the individual gene level, there are many sex-specific effects; however, at the annotation level, results are more concordant.

CONCLUSIONS

In addition to demonstrating sex-specific effects of selection on the transcriptome, the data point to the involvement of extracellular matrix genes as being associated with the binge drinking phenotype.

Baseline mRNA expression differs widely between common laboratory strains of zebrafish

Common strains of wildtype zebrafish (Danio rerio) have unique genomic features including SNPs and CNV, but strain information often goes unreported in the literature. As a result, the confounding effects of interstrain variation makes repetition of studies in zebrafish challenging. Here we analyze hepatic mRNA expression patterns between three common zebrafish strains (AB, Tuebingen (TU), and WIK) using Agilent 4 × 44 K gene expression microarrays to establish baseline mRNA expression across strains and between sexes. We observed wide variation in sex-specific gene expression within AB and WIK strains (141 genes in AB and 67 genes in WIK), but no significant variation between sexes within TU. After partitioning the dataset into male and female subsets, we detected 421 unique mRNA transcripts with statistically significant differential expression; 269 mRNA transcripts varied between males, 212 mRNA transcripts varied between females, and 59 mRNA transcripts varied across the three strains, regardless of sex. It is not surprising that mRNA expression profiles differ between sexes and strains, but it is imperative to characterize the differences. These results highlight the complexity of variation within zebrafish and underscore the value of this model system as a valid representation of normal variation present in other species, including humans.

Concentration and timing of application reveal strong fungistatic effect of tebuconazole in a Daphnia-microparasitic yeast model

Given the importance of pollutant effects on host-parasite relationships and disease spread, the main goal of this study was to assess the influence of different exposure scenarios for the fungicide tebuconazole (concentration×timing of application) on a Daphnia-microparasitic yeast experimental system. Previous results had demonstrated that tebuconazole is able to suppress Metschnikowia bicuspidata infection at ecologically-relevant concentrations; here, we aimed to obtain an understanding of the mechanism underlying the anti-parasitic (fungicidal or fungistatic) action of tebuconazole. We exposed the Daphnia-yeast system to four nominal tebuconazole concentrations at four timings of application (according to the predicted stage of parasite development), replicated on two Daphnia genotypes, in a fully crossed experiment. An "all-or-nothing" effect was observed, with tebuconazole completely suppressing infection from 13.5μgl^-1 upwards, independent of the timing of tebuconazole application. A follow-up experiment confirmed that the suppression of infection occurred within a narrow range of tebuconazole concentrations (3.65-13.5μgl^-1), although a later application of the fungicide had to be compensated for by a slight increase in concentration to elicit the same anti-parasitic effect. The mechanism behind this anti-parasitic effect seems to be the inhibition of M. bicuspidata sporulation, since tebuconazole was effective in preventing ascospore production even when applied at a later time. However, this fungicide also seemed to affect the vegetative growth of the yeast, as demonstrated by the enhanced negative effect of the parasite (increasing mortality in one of the host genotypes) at a later time of application of tebuconazole, when no signs of infection were observed. Fungicide contamination can thus affect the severity and spread of disease in natural populations, as well as the inherent co-evolutionary dynamics in host-parasite systems.

Genes linked to species diversity in a sexually dimorphic communication signal in electric fish

Sexually dimorphic behaviors are often regulated by androgens and estrogens. Steroid receptors and metabolism are control points for evolutionary changes in sexual dimorphism. Electric communication signals of South American knifefishes are a model for understanding the evolution and physiology of sexually dimorphic behavior. These signals are regulated by gonadal steroids and controlled by a simple neural circuit. Sexual dimorphism of the signals varies across species. We used transcriptomics to examine mechanisms for sex differences in electric organ discharges (EODs) of two closely related species, Apteronotus leptorhynchus and Apteronotus albifrons, with reversed sexual dimorphism in their EODs. The pacemaker nucleus (Pn), which controls EOD frequency (EODf), expressed transcripts for steroid receptors and metabolizing enzymes, including androgen receptors, estrogen receptors, aromatase, and 5α-reductase. The Pn expressed mRNA for ion channels likely to regulate the high-frequency activity of Pn neurons and for neuromodulator and neurotransmitter receptors that may regulate EOD modulations used in aggression and courtship. Expression of several ion channel genes, including those for Kir3.1 inward-rectifying potassium channels and sodium channel β1 subunits, was sex-biased or correlated with EODf in ways consistent with EODf sex differences. Our findings provide a basis for future studies to characterize neurogenomic mechanisms by which sex differences evolve.

Both Male-Biased and Female-Biased Genes Evolve Faster in Fish Genomes

Males and females often display extensive phenotypic differences, and many of these sexual dimorphisms are thought to result from differences between males and females in expression of genes present in both sexes. Sex-biased genes have been shown to exhibit accelerated rates of evolution in a wide array of species, however the cause of this remains enigmatic. In this study, we investigate the extent and evolutionary dynamics of sex-biased gene expression in zebrafish. Our results indicate that both male-biased genes and female-biased genes exhibit accelerated evolution at the protein level. In order to differentiate between adaptive and nonadaptive causes, we tested for codon usage bias and signatures of different selective regimes in our sequence data. Our results show that both male- and female-biased genes show signatures consistent with adaptive evolution. In order to test the generality of our findings across fish, we also analyzed publicly available data on sticklebacks, and found results consistent with our findings in zebrafish.

Sex differences in DNA methylation and expression in zebrafish brain: a test of an extended 'male sex drive' hypothesis

The sex drive hypothesis predicts that stronger selection on male traits has resulted in masculinization of the genome. Here we test whether such masculinizing effects can be detected at the level of the transcriptome and methylome in the adult zebrafish brain. Although methylation is globally similar, we identified 914 specific differentially methylated CpGs (DMCs) between males and females (435 were hypermethylated and 479 were hypomethylated in males compared to females). These DMCs were prevalent in gene body, intergenic regions and CpG island shores. We also discovered 15 distinct CpG clusters with striking sex-specific DNA methylation differences. In contrast, at transcriptome level, more female-biased genes than male-biased genes were expressed, giving little support for the male sex drive hypothesis. Our study provides genome-wide methylome and transcriptome assessment and sheds light on sex-specific epigenetic patterns and in zebrafish for the first time.

Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish

Background

Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored.

Methods

RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females.

Results

Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed.

Conclusions

Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.

Genes located in a chromosomal inversion are correlated with territorial song in white-throated sparrows

The genome of the white-throated sparrow (Zonotrichia albicollis) contains an inversion polymorphism on chromosome 2 that is linked to predictable variation in a suite of phenotypic traits including plumage color, aggression and parental behavior. Differences in gene expression between the two color morphs, which represent the two common inversion genotypes (ZAL2/ZAL2 and ZAL2/ZAL2(m) ), may therefore advance our understanding of the molecular underpinnings of these phenotypes. To identify genes that are differentially expressed between the two morphs and correlated with behavior, we quantified gene expression and terrirorial aggression, including song, in a population of free-living white-throated sparrows. We analyzed gene expression in two brain regions, the medial amygdala (MeA) and hypothalamus. Both regions are part of a 'social behavior network', which is rich in steroid hormone receptors and previously linked with territorial behavior. Using weighted gene co-expression network analyses, we identified modules of genes that were correlated with both morph and singing behavior. The majority of these genes were located within the inversion, showing the profound effect of the inversion on the expression of genes captured by the rearrangement. These modules were enriched with genes related to retinoic acid signaling and basic cellular functioning. In the MeA, the most prominent pathways were those related to steroid hormone receptor activity. Within these pathways, the only gene encoding such a receptor was ESR1 (estrogen receptor 1), a gene previously shown to predict song rate in this species. The set of candidate genes we identified may mediate the effects of a chromosomal inversion on territorial behavior.

Neurotranscriptome profiles of multiple zebrafish strains

Behavioral displays or physiological responses are often influenced by intrinsic and extrinsic mechanisms in the context of the organism's evolutionary history. Understanding differences in transcriptome profiles can give insight into adaptive or pathological responses. We utilize high throughput sequencing (RNA-sequencing) to characterize the neurotranscriptome profiles in both males and females across four strains of zebrafish (Danio rerio). Strains varied by previously documented differences in stress and anxiety-like behavioral responses, and generations removed from wild-caught individuals. Here we describe detailed methodologies and quality controls in generating the raw RNA-sequencing reads that are publically available in NCBI's Gene Expression Omnibus database (GSE61108).

Zebrafish sexual behavior: role of sex steroid hormones and prostaglandins

Background

Mating behavior differ between sexes and involves gonadal hormones and possibly sexually dimorphic gene expression in the brain. Sex steroids and prostaglandin E₂ (PGE₂) have been shown to regulate mammalian sexual behavior. The present study was aimed at determining whether exposure to sex steroids and prostaglandins could alter zebrafish sexual mating behavior.

Methods

Mating behavior and successful spawning was recorded following exposure to 17β-estradiol (E2), 11-ketotestosterone (11-KT), prostaglandin D₂ (PGD₂) and PGE₂ via the water. qRT-PCR was used to analyze transcript levels in the forebrain, midbrain, and hindbrain of male and female zebrafish and compared to animals exposed to E2 via the water.

Results

Exposure of zebrafish to sex hormones resulted in alterations in behavior and spawning when male fish were exposed to E2 and female fish were exposed to 11-KT. Exposure to PGD_2, and PGE₂ did not alter mating behavior or spawning success. Determination of gene expression patterns of selected genes from three brain regions using qRT-PCR analysis demonstrated that the three brain regions differed in gene expression pattern and that there were differences between the sexes. In addition, E2 exposure also resulted in altered gene transcription profiles of several genes.

Conclusions

Exposure to sex hormones, but not prostaglandins altered mating behavior in zebrafish. The expression patterns of the studied genes indicate that there are large regional and gender-based differences in gene expression and that E2 treatment alter the gene expression pattern in all regions of the brain.

Electronic supplementary material

The online version of this article (doi:10.1186/s12993-015-0068-6) contains supplementary material, which is available to authorized users.

Characterizing the neurotranscriptomic states in alternative stress coping styles

BACKGROUND

Animals experience stress in many contexts and often successfully cope. Individuals exhibiting the proactive versus reactive stress coping styles display qualitatively different behavioral and neuroendocrine responses to stressors. The predisposition to exhibiting a particular coping style is due to genetic and environmental factors. In this study we explore the neurotranscriptomic and gene network biases that are associated with differences between zebrafish (Danio rerio) lines selected for proactive and reactive coping styles and reared in a common garden environment.

RESULTS

Using RNA-sequencing we quantified the basal transcriptomes from the brains of wild-derived zebrafish lines selectively bred to exhibit the proactive or reactive stress coping style. We identified 1953 genes that differed in baseline gene expression levels. Weighted gene coexpression network analyses identified one gene module associated with line differences. Together with our previous pharmacological experiment, we identified a core set of 62 genes associated with line differences. Gene ontology analyses reveal that many of these core genes are implicated in neurometabolism (e.g. organic acid biosynthetic and fatty acid metabolic processes).

CONCLUSIONS

Our results show that proactive and reactive stress coping individuals display distinct basal neurotranscriptomic states. Differences in baseline expression of select genes or regulation of specific gene modules are linked to the magnitude of the behavioral response and the display of a coping style, respectively. Our results expand the molecular mechanisms of stress coping from one focused on the neurotransmitter systems to a more complex system that involves an organism's capability to handle neurometabolic loads and allows for comparisons with other animal taxa to uncover potential conserved mechanisms.