Research | Open | Published:
Expression profiles for six zebrafish genes during gonadal sex differentiation
Reproductive Biology and Endocrinologyvolume 6, Article number: 25 (2008)
The mechanism of sex determination in zebrafish is largely unknown and neither sex chromosomes nor a sex-determining gene have been identified. This indicates that sex determination in zebrafish is mediated by genetic signals from autosomal genes. The aim of this study was to determine the precise timing of expression of six genes previously suggested to be associated with sex differentiation in zebrafish. The current study investigates the expression of all six genes in the same individual fish with extensive sampling dates during sex determination and -differentiation.
In the present study, we have used quantitative real-time PCR to investigate the expression of ar, sox9a, dmrt1, fig alpha, cyp19a1a and cyp19a1b during the expected sex determination and gonadal sex differentiation period. The expression of the genes expected to be high in males (ar, sox9a and dmrt1a) and high in females (fig alpha and cyp19a1a) was segregated in two groups with more than 10 times difference in expression levels. All of the investigated genes showed peaks in expression levels during the time of sex determination and gonadal sex differentiation. Expression of all genes was investigated on cDNA from the same fish allowing comparison of the high and low expressers of genes that are expected to be highest expressed in either males or females. There were 78% high or low expressers of all three "male" genes (ar, sox9a and dmrt1) in the investigated period and 81% were high or low expressers of both "female" genes (fig alpha and cyp19a1a). When comparing all five genes with expected sex related expression 56% show expression expected for either male or female. Furthermore, the expression of all genes was investigated in different tissue of adult male and female zebrafish.
In zebrafish, the first significant peak in gene expression during the investigated period (2–40 dph) was dmrt1 at 10 dph which indicates involvement of this gene in the early gonadal sex differentiation of males.
The mechanism underlying sex determination in zebrafish is largely unknown. Currently, no morphological differences in the chromosomes of the two sexes have been identified neither by classical karyotyping nor from the zebrafish genome sequencing project. This indicates that sex determination in zebrafish is mediated by genetic signals from autosomal genes and not by XY/XX or ZZ/ZW sex determination [1–3]. Furthermore, the mammalian sex determining factor Sry (located on the Y chromosome) and the sex determining gene dmy (dmrt1y) located on the Y chromosome in Japanese medaka Oryzias latipes [4, 5] does not appear to be present in the zebrafish genome. Current knowledge on sex determination and differentiation in zebrafish includes that zebrafish are sexually mature after three months, and that separate sexes can be detected by gonad histology at approximately 40 days post hatch (dph). Furthermore, all zebrafish develop ovary-like gonads regardless of genetic background prior to sex differentiation . Ovarian development is initiated at approximately 10 dph and progresses until 20 dph. At 20 dph until approximately 30 dph testis development is initiated in males simultaneously with ovarian cell apoptosis [2, 3].
A number of genes (dmrt1, sox9a, amh, wt1, ftz-f1, gata) have previously been associated with the process of sex determination or differentiation in zebrafish . None of these genes are suggested to be the single factor responsible for specifying sex in zebrafish . However, the expression pattern and function of these genes suggest that they are part of a signalling network responsible for the development of sex specific gonads . Current research in medaka Oryzias latipes suggests a potential signalling pathway in XY individuals of: dmy, sox9a and dmrt1 with dmy blocking meiosis, sox9a regulating testicular tubule development and dmrt1 being important in spermatogenesis . Also, the research in medaka Oryzias latipes suggests an important role of fig α as well as cyp19a1a and cyp19a1b in sexual development of gonads in XX individuals .
The Dmrt1 gene belongs to a group of multiple genes containing a zinc-finger-like DNA-binding motif (DM domain) that has been identified in both invertebrates and vertebrates . Interestingly, the medaka sex determining gene DMY (dmrt1y) originates from a dmrt1 duplication, and sex determining homologues have been identified in the fruit fly Drosophila melanogaster (doublesex gene) and the worm Caenorhabditis elegans (mab-3 gene) . In zebrafish, in situ hybridization to gonads showed that dmrt1 was expressed in developing germ cells of both testis and ovary, suggesting that the dmrt1 gene is not only associated with testis development, but may also be important in ovary differentiation of zebrafish . In diverse species including frog, turtle, alligator, bird, and mouse, the dmrt1 gene is expressed at higher levels in males compared to females, suggesting that high expression is necessary for testicular differentiation, whereas lower expression is compatible with ovarian differentiation [9–12]. The Sox (SRY-related genes containing a HMG box) gene family encodes an important group of developmental regulators involved in sex determination. The HMG (high mobility group) box that characterises Sox proteins is a DNA-binding domain and proteins encoded by Sox genes act as transcription factors. Sry, the founder member of the Sox gene family, is the Y chromosomal male determinant in most mammals [13–15]Sry is a poorly conserved gene that appears to be exclusive to mammals. In contrast, Sox 9 is a conserved gene present in all vertebrate types. Like Sry, Sox 9 is required for testis development in mammals, and Sox9 deficiency can result in sex reversal in human males. The expression of Sox9 during gonadal differentiation is up-regulated in testis and down-regulated in ovaries in mammals, birds and turtles . However, the organisation and function of the Sox gene family is less understood in other types of vertebrates and despite the wide distribution of Sox9 genes in fish, only few have been investigated [16, 17]. In zebrafish, two sox9 genes are present sox9a and sox9b) and their expression patterns indicate that they have unique functions during development . In adult zebrafish the sox9a transcript was observed in testis but not in ovary. Conversely, sox9b transcripts were detected in ovary, but not in testis .
Androgens play a key role in male sex differentiation and sex maturation in vertebrates, including teleosts . The most important androgens in teleost fish are 11-ketotestosterone and testosterone and their action is mediated through specific nuclear androgen receptors (ar) [19–21]. In sea bass (Dicentrarchus labrax) expression of the ar was measured in the gonads during development of a male-dominated population and a female-dominated population . The expression patterns of the two populations were different, with a peak in the ar expression in the male-dominated population coinciding with the time of sex differentiation in sea bass . We recently identified and conducted the initial characterisation of a novel androgen receptor from zebrafish  and in another recent study, expression of this receptor was higher in the transforming testis compared to ovary suggesting a role during male gonadal differentiation . A key enzyme balancing the ratio of steroid hormones is cyp19 (aromatase) which is the terminal enzyme in the steroidogenic pathway. It converts androgens (e.g., testosterone) into estrogens (e.g., estradiol). Regulation of this gene dictates the ratio of androgens to estrogens; therefore, appropriate expression of this enzyme is critical for sex differentiation and reproduction in vertebrates. Most vertebrates have a single CYP19 gene that is regulated by multiple tissue-specific promoter regions. However, the zebrafish has two genes (cyp19a1a and cyp19a1b) encoding different proteins and each possessing its own regulatory mechanism . In general, the gonadal form cyp19a1a is more abundant than the brain form (cyp19a1b). The expression of the two cyp19 genes has previously been investigated from 0–41 days post fertilisation (dpf) which is the expected time of sex determination and differentiation in zebrafish. cyp19a1a expression was highest shortly after hatch from 4–8 dpf. The pattern of cyp19a1b expression was segregated into two populations, suggesting an association with sex differentiation . fig α is a germ cell-specific transcription factor required for ovarian follicle formation. fig α is involved in the coordinate expression of the zona pellucida (zpc) genes. The expression of fig α and zpc coincides with the onset of gonadal differentiation in zebrafish at 22 dpf . Previous in situ hybridisation studies of fig α expression in zebrafish have shown that fig α is expressed abundantly in ovaries of adult fish whereas no fig α signal could be detected in adult zebrafish testes . Furthermore, at 30 dpf the expression of fig α and zpc was investigated in eight fish. In five fish both genes were expressed whereas no expression was detected in the last three fish. This indicates female restricted expression of fig α in zebrafish .
Since previous studies indicate involvement of ar, sox9a, dmrt1, cyp19a1b, cyp19a1a and fig α in sex differentiation in vertebrates, including fish, we found it interesting to investigate the expression of these six genes during sex determination and differentiation in zebrafish. The idea is to measure the expression of all genes in the same individual, thereby allowing comparison of expression. This might lead to identification of putative male and female zebrafish and indicate whether one or more of these genes could function as an early genetic sex marker in zebrafish. Therefore, the aim of the current study is to determine the expression of the ar, sox9a, dmrt1, cyp19a1b, cyp19a1a and fig α during the developmental period in which sex determination and gonadal differentiation takes place in zebrafish.
Juvenile zebrafish originated from a brood population of fish. In the evening breeding boxes were placed in an aquarium with parent fish and eggs were collected the following morning. The eggs were sorted into fertilised and non-fertilised and the fertilised eggs were placed in 900 ml glass beakers and kept at 26 +/- 1°C and a light-dark period of 14:10 h. In the interval 3–22 dph the larvae were fed two times daily with powdered dry food (Sera Micron) and one time daily with newly hatched artemia sp. nauplii (Intér Ryba GmbH, Germany). For a period of three days (23–25 dph) the powdered dry food was given in combination with TetraMin Baby Powder Food. From 26 dph and onwards the dry food consisted solely of TetraMin Baby Powder Food. Artemia was still given once daily. At 2, 4, 6, 8, 10, 12, 14, 16, 18, 19, 20, 21, 22, 25, 30 and 40 dph zebrafish were frozen individually in liquid nitrogen and stored at -80°C until further analysis.
Total RNA was purified from whole fish for the AR expression during sex determination and differentiation experiment. Zebrafish or zebrafish tissue was homogenised for 20 sec. using an Ultra-Turrax homogenizer (IKA-Werke). For purification of RNA, Total RNA Isolation Kit (Macherey-Nagel) was used according to manufacturer's instructions. RNA concentration was measured spectrophotometrically (Gene-Quant, Pharmacia Biotech), checked by gel electrophoresis (1.2% agarose gel) and stored at -80°C until further use. cDNA was obtained from 0.3 μg RNA using SuperScript II reverse transcriptase kit with Oligo dT primer (Invitrogen) according to manufacturers instructions.
Expression of Danio rerio ar, sox9a, dmrt1, cyp19a1b, fig α and cyp19a1a was analysed by qRT-PCR using a real-time light-cycler (Roche). As zebrafish is an established model-organism, the genome has been sequenced and the seventh assembly is available. All genes investigated have previously been identified and sequences were available in GenBank. Primers for qRT-PCR analysis were designed using the Primer3 program  (Table 1). The final PCR reactions contained: 0.4 mM of each primer; 0.25 × SYBR Green (Invitrogen); 4 mM MgCl2 and as template 5 μl of cDNA reverse transcribed from a standardized amount of total RNA (0.3 μg). qRT-PCR was performed using Hotstart Taq polymerase (Qiagen) in a final volume of 20 μl. All quantitative reactions were subjected to: 95°C for 15 min followed by 45 cycles at 94°C for 15 s, 59°C 15 s and 72°C 15 s. Melting curve analysis was applied to all reactions to ensure homogeneity of the reaction product. In addition, the amplicon size was checked electrophoretically for each primer set and subsequent sequencing revealed that it corresponded to the zebrafish ar, sox9a, dmrt1, cyp19a1b, fig α and cyp19a1a, thus verifying the identity of the genes. Potential contamination was assessed by including non-reverse transcribed total RNA (genomic DNA contamination) and no-template controls. No products were observed in these reactions. Dilution curves generated by serial dilutions (1:10) of cDNA were used to calculate amplification efficiencies according to the Roche protocol. All assays were quantitative with standard curve (mean threshold cycle [Ct] vs. log cDNA dilution) slopes of -3.96 (β-actin), -2.10 (sox9a), -3.39 (dmrt1), -2.88 (ar), -3.00 (cyp19a1a), -4.22 (cyp19a1b) and -3.68 (fig α), translating to relatively high PCR efficiencies (E) of 1.79 (β-actin), 3.01 (sox9a), 1.97 (dmrt1), 2.23 (ar), 2.15 (cyp19a1a), 1.73 (cyp19a1b) and 1.87 (fig α). over the detection range, the linear correlation (R2) between the mean Ct and the log cDNA dilution was > 0.99 in each case. Transcript levels of the target genes were normalized to D. rerio β-actin after correcting for differences in amplification efficiencies.
Adult zebrafish, 5 males and 5 females were anaesthetized in a buffered solution of MS-222 (0.1 g/l) and quickly dissected into brain, gonads, liver, eyes, spleen, heart, gut, gall bladder, muscle and gills. The dissected tissues were immediately frozen in liquid nitrogen and stored at -80°C. The RNA purification and quantitative qRT-PCR was conducted as described for the gene expression during sex determination and differentiation experiment except that approximately 5 mg of zebrafish tissue were used for total RNA purification. Also, 1 μg total RNA was reverse transcribed using SuperScript II reverse transcriptase kit with Oligo dT primer (Invitrogen) according to manufacturers instructions. Following the qRT-PCR, the reactions were spinned down and loaded on a 0.8% agarose gel.
Data handling and statistical analysis
To analyse the expression patterns of ar, sox9a, dmrt1, fig α, cyp19a1a and cyp19a1b during the investigated period, it was necessary to discriminate between individual fish with high and low gene expression. This was done using the GraphPad Prism program makes a cut-off value between the high and low expressers as indicated on the graph, which is generated automatically in the program when making a two-segment graph. The investigated fish were a mixture of males and females and since large differences in expression were observed during the investigated period, it was difficult to separate the different effects from each other without this distinction between high and low expressers.
As expression of all the six investigated genes has been measured on individual fish, it is possible to compare the expression patterns. First, ar, sox9a and dmrt1 genes with expected high expression in males were divided into high and low expressers, then for each individual fish the expression levels (high or low) of the three genes were compared. Next, the percentage of individuals with high or low expression of the three genes was calculated. Secondly, the same was calculated for the two genes (fig α and cyp19a1a) with expected high expression in females. Third, the percentage of individuals with expression patterns as expected for male or female for all five genes was calculated.
Gene expressions in zebrafish during the sex determination and differentiation period were statistically analysed for variance. Data were normally distributed with similar variance after log transformation and a two-way ANOVA (variables: high/low expressers and dph) analysis was followed by Tukey test with a significance level of (p < 0.05).
To analyse the expression patterns of ar, sox9a, dmrt1a, fig α, cyp19a1a and cyp19a1b during the expected sex determination and differentiation period, individuals were divided into groups of high and low expressers based on the method described previously where the distribution between investigated individuals and the cut-off value between high and low expressers are shown in Figure 1 for each of the investigated genes. The expression pattern of ar, sox9a, dmrt1, fig α and cyp19a1a were compared in each individual zebrafish, three genes with expected high expression in males (ar, sox9a, dmrt1) and two genes with expected high expression in females (fig α, cyp19a1a). The data show that for ar, sox9a and dmrt1 expression 78% were either high or low expressers of all three genes. For fig α and cyp19a1a expression, 81% were either high or low expressers of both genes. Individuals with expression pattern of all five genes that were as expected for either male or female corresponded to 56% and square symbols were used for the tentative identification of males and triangles were used for the tentative identification of females. Individuals in which the expression of all five genes was not as expected for either female or male corresponded to 44% and were called intermediate and a round symbol was used. cyp19a1b was not included in the tentative identification of sex and therefore the round symbol was used. There was always more than 10 times difference in expression between the two groups and statistical analysis showed that for all investigated genes there was a significant difference (p < 0.05) between high and low expresser groups. The ratio between genes expected to be primarily "female" and genes expected to be primarily "male" were calculated (Table 2) together with ratios for fig α/dmrt1.
The expression levels of ar were low for both the high and low expresser group during the first period (2–8 dph) (Figure 2), however, with two distinct groups of high and low expressers. The ar expression has a peak at 16 dph, which is significantly different (p < 0.05) compared to the expression on the day before. The ar expression reaches the highest level at 22 dph for the high-expresser group, also significantly different (p < 0.05) from the expression on the day before. The 22 dph peak corresponds to the time period right after the bipotential gonad is starting to differentiate into ovary or testes. In the low expresser group no increase in ar levels was seen. In the period 19–20 dph there is a marked decrease in the expression level for the high-expresser group, which coincides with the time of juvenile ovary-to-testes transformation in zebrafish.
The expression of sox9a is very low for the low expresser group during the entire investigated period. In the high expresser group there are three peaks, at 4, 18 and 22 dph (Figure 3). The peak in the high expresser group at 4 dph is very early in development and the expression levels at 2 dph and 4 dph are not significantly different (p > 0.05). Interestingly, the peak at 18 dph which is significantly different compared to the expression of sox9a at 16 dph is just before the time of oocyte apoptosis (20–23 dph) in individuals expected to develop into males and might therefore represent a signal in the male developmental pathway. The peak at 18 dph is followed by a rapid decrease and very low sox9a expression at 19 dph. At 22 dph the expression level in the high-expresser group peaks again with expression level significantly different (p < 0.05) from that at 21 dph. The 22 dph peak in sox9a expression coincide with an expression peak in ar in the high expresser group.
The high-expresser group of dmrt1 has two distinct expression peaks (at 10 and 14 dph) during the investigated period (Figure 4). Both peaks are significantly different (p < 0.05) compared to the expression during previous days in the high expresser group. The peak in dmrt1 expression at 10 dph is early in development and might represent an early marker of male development. The peak at 14 dph is also quite early in development and these results could indicate that dmrt1 expression is important in juvenile sex determination and/or differentiation. The 14 dph peak is followed by several days (21–22 dph) with low dmrt1 expression which is coinciding with low expression of ar and sox9a in the high expresser group. The low-expresser group shows a relatively low expression of dmrt1 during the entire period, but there might be an indication that at 10 and 16 dph the dmrt1 expression decreases in the low-expresser group which is opposite to the expression in the high-expresser group.
In the group of high-expressers of cyp19a1b, one distinct expression peak is seen at 14 dph, the level being significantly different (p < 0.05) compared to expression at 12 dph (Figure 5). During the period 18–30 dph small oscillations in expression are seen in the high-expresser group and from 30–40 dph the cyp19a1b expression increases. The levels of cyp19a1b expression in the low-expresser group are decreasing in the period 6–21 dph. Thereafter the expression level remains very low.
fig α expression
The expression of fig α is low with small oscillations in the high expresser groups until 22 dph and the expression in the low expresser groups is almost undetectable during the entire experimental period. A very distinct peak in fig α expression in the high expresser group is seen at 25 dph, significantly different (p < 0.05) compared to the expression at 22 dph (Figure 6). The peak in fig α expression coincides with the expected onset of gonadal differentiation to ovary in female zebrafish.
The cyp19a1a expression in the low expressers is low during the entire investigated period whereas two peaks in expression are seen in the high-expresser groups, at 18 dph and 30 dph, respectively (Figure 7). For both peaks the expression is significantly different (p < 0.05) compared to the expression on the previous day. The peak at 18 dph is just prior to the initiation of gonadal differentiation in zebrafish and the peak at 30 dph is just after the expected time of gonadal differentiation and might reflect an increased estrogen demand in the gonads.
Tissue expression in adult zebrafish
The tissue specific expression of ar, sox9a dmrt1, fig α, cyp19a1a and cyp19a1b in adult male and female zebrafish was investigated by qRT-PCR to determine whether the genes were uniformly expressed or restricted to specific tissues or sex. Investigated tissues included liver, heart, brain, spleen, muscle, gall bladder, gut, eyes, gill and gonad (Figure 8). ar was expressed in all of the ten different tissues investigated in males and females. The sox9a expression was similar in the two sexes except in gonads and gill where it was low in female fish compared to the expression in male fish indicating that this "male" gene still has a role in the adult phase. dmrt1 was expressed in spleen and gonads in both sexes, but at highest level in males. The only tissue with expression of the "female" gene fig α was female gonads. cyp19a1a and cyp19a1b were differentially expressed in different organs of male and female zebrafish, liver expression of cyp19a1a being higher in male than in female, and higher in female gonads than in male gonads. cyp19a1b expression was relatively higher in female than in male brain.
Testicular type sox9 is the most upstream conserved gene in the sex determining cascade among vertebrates . Therefore, the expression pattern of this gene during zebrafish sex determination and differentiation is important. In non-mammalian vertebrates like birds and turtles, sox9 is expressed in testes and down-regulated in ovaries during gonadal differentiation similar to the expression pattern seen in mammals [30, 31]. In the present study, expression of sox9a is found in several tissues including spleen, muscle, gall bladder, brain and eye of both male and female zebrafish. However, sox9a is also expressed in gonad and gill in a male-specific manner. This is in accordance with a previous study in adult zebrafish showing different expression patterns for sox9a and sox9b in zebrafish gonads suggesting that sox9a retained its function in the testis while sox9b possibly acquired a different function in zebrafish ovary during evolution . The expression pattern of sox9a in zebrafish during sex determination and differentiation segregated the fish into two groups. In fish of the high expresser group three peaks were observed. The first small peak seen at 4 dph is early in development and might be related to skeletal development, as sox9 in mammals is also involved in this process , but it could also be related to sex determination and differentiation. The peaks at 18 and 22 dph coincide with the expected time of gonadal differentiation, i.e. ovary to testis transformation in genetic males, which is a key event in zebrafish sex differentiation. In the mammalian bipotential gonad, Sry initiates sox9 expression and translocation from cytoplasm to nucleus which induces expression of amh that is an important testis determining factor . The expression pattern of sox9a seen in this study could indicate involvement of sox9a in sex differentiation in zebrafish. However, no available data including those of the present study firmly indicates that sox9a is the sex determining gene in zebrafish, but merely suggests that it might be involved in the sex signalling cascade and gonadal sex differentiation. This is in agreement with a study in medaka where the level of testicular type sox9 (sox9a2) expression in somatic cells is equally high in both sexes at the time when dmy expression is initiated during early gonadal differentiation. However, during the period from 10–30 dph, sox9 expression continues only in the Sertoli cells in male gonads, with a marked reduction in the XX gonads .
In a study investigating amh, sox9a, sox9b and cyp19a1a expression in undifferentiated gonads of zebrafish, expression of all genes could be detected , however, in the differentiated gonads a sexually dimorphic expression pattern was found; sox9a and amh were expressed in testis whereas cyp19a1a was not. In ovaries, sox9b and cyp19a1a were expressed while amh was not . Based on the expression pattern of these genes during sex differentiation, the authors suggested that 17 dph represents a transitional stage in zebrafish gonad development and by 31 dph gonads have differentiated into testes or ovaries . This is in accordance with the present study in which distinct peaks in both cyp19a1a and sox9a expression was found at 18 dph. Furthermore, these results are in agreement with those of a recent study of the molecular mechanism of the ovary-to-testes transition which indicated that all males go through the juvenile ovary phase until approximately 21 dpf (corresponding to 19 dph). However, they differ in the extent of commitment toward femaleness during this period and can be divided into three types of males based on the intensity, onset and duration of gonadal transformation .
In a number of species including humans, mice, chickens, alligators and turtles, the dmrt1 expression is limited to the gonads, and the expression is considerably up-regulated in the developing testes compared to ovaries. The timing of this up-regulation varies between species, but generally occurs in the late sex determining or early testis-differentiation period. This characteristic dmrt1 expression pattern seen in different species indicates that this gene is specifically involved in the early formation of testes . The results from the present study correspond well with this general notion. We found peaks in expression of dmrt1 in the high expresser group prior to the expected time of ovary to testis transformation. The expression peaks of dmrt1 early in the sex determination and differentiation period could indicate that it is involved upstream in the signalling cascade that initiates sex determination and differentiation in zebrafish. In a recent model of medaka gonadal development it is shown that the dmrt1 expression is detected from 10 dph and is thereby the first gene to be differentially expressed between males and females except for DMY . This is in agreement with the results from the present study where we see a peak in Dmrt1a expression at 10 dph. However, the idea that the dmrt1 gene might be the sex-determining gene in non-mammalian vertebrates has been rejected based on the lack of DMY (Dmrt1bY) in other fish species than Oryzias latipes, including Oryzias celebensis, Oryzias mekongenesis, guppy Poecilia reticulata, pufferfish Takifugu rubripes and zebrafish Danio rerio [34, 36]. The dmrt1 gene most likely has a role downstream the sex determination event and it might be involved in testis development in teleost fish analogous to its putative role in mammalian species . This is in accordance with results from medaka where results indicated that dmrt1 regulates spermatogonial differentiation . Sexual development in sex-reversed medaka gonads indicates that DMY is not necessary for gonadal differentiation and that expression of fig α and dmrt1 correlate with the phenotypic differentiation of the gonads . The expression pattern of fig α in the present study resembles to some extent that in the hermaphrodite fish Kryptolebias marmoratus, where fig α expression was very low until 39 dpf where expression peaks and remains high until 103 dpf which was the last day measured . Furthermore, fig α showed no expression in males during gonadogenesis in mice  which might be in accordance with the results in this study where we see low expression (close to the detection level) in the low expresser group. In mice, fig α has been suggested to play a key role in preserving oocytes and normal ovarian development  and in medaka fig α is an oocyte specific marker .
In this study, the pattern of ar expression in zebrafish during sex determination and differentiation clearly segregated the fish into two groups. This is in accordance with a study in sea bass where differences in ar expression were first encountered at 150 dph and became especially marked at 250 dph (corresponding with the time of sex differentiation) with higher expression in a male dominant group compared to a female dominant group (male and female dominant groups were based on consecutive size grading) . In the current study, we find peaks in ar expression in the high expresser group at 16 and 25 dph which is shortly after a peak in expression of dmrt1 is observed and shortly before a peak in sox9a expression. These peaks in expression of the different genes around 10–18 dph coincide with the expected initiation of transformation from ovary to testis in individuals developing into males. The difference in expression profiles between the two groups at the time of sex differentiation is most likely sex-related and corresponds to differences between males and females, suggesting an important role for ar in sex differentiation. In a recent study of zebrafish gonad development, Hossain et al (2007) found ar mRNA levels at 4 weeks post fertilization (their first measurement) to be similar in males and females. However, from 5–7 wpf the ar expression was higher in the transforming testis . Furthermore, expression levels of ar mRNA at different developmental stages has been studied in hermaphrodite fish species such as the protogynous wrasse, Halichoeres trimaculatus  and the protandrous black porgy, Acanthopagrus schlegeli . In the black porgy, both the testicular tissue of the bisexual gonad and the functional testis exhibited higher ar mRNA levels than the ovarian tissue. Conversely, ar mRNA levels were lower in functional ovaries than in the ovarian part of the bisexual gonad, suggesting association of the decrease in ar mRNA levels with protandrous sex change in this fish species . The ubiquitous expression of ar in adult zebrafish of both sexes found in this study is in accordance with  and we found ar expression in all investigated tissue of adult zebrafish (gonad, brain, liver, kidney, skin, muscle, eye) of both males and females.
Also, the ovarian form of aromatase, cyp19a1a was suggested to be involved in sex determination and differentiation . cyp19a1a expression was found in undifferentiated gonads at 17 dpf which is expected based on the juvenile ovarian gonads. However, at 31 dpf the differentiated ovary showed specific expression of cyp19a1a in cells surrounding the oocytes . Furthermore, at 31 dpf the differentiated testis showed no expression of cyp19a1a . The results presented in this study correspond well with the previous studies of cyp19a1a expression. We find two peaks in expression during the sex determination and differentiation period (at 18 and 30 dph) in the high expresser group. Furthermore, we find that the expression of cyp19a1a in both the high and low expresser group is very low (almost undetectable in the low expresser group) at the time of gonadal transformation at approximately 19–21 dph. In accordance, a recent study showed cyp19a1a expression around oocytes in juvenile ovary, this expression was down-regulated and could no longer be detected when gonadal transformation was initiated .
The inter-individual grouping of the expression of cyp19a1b during sex determination and differentiation into high and low expressers observed in the present study is in accordance with previous studies [1, 27, 43]. In the present study we see one distinct peak in cyp19a1b expression at 14 dph for the high expresser group. This does not correspond completely with the peaks found in the high expresser group during sex determination and differentiation in the study by Trant et al. (2001), where cyp19a1b peaks were seen at 4–5 dpf, 13 dpf (small peak) and 25 dpf and from 34 – 40 dpf increase (with 2.5 dpf = 0 dph) . The small peak at 13 dpf is most likely corresponding to the peak we see at 14 dph. However, the general expression pattern is not similar in the two studies. The finding that individual cyp19a1b expression during sex determination and differentiation distributes into two groups indicates that the expression of cyp19a1b (the brain form) might be associated with sexual differentiation in zebrafish as previously suggested by Trant et al. (2001). In contrast, Gato-Kazeto et al. (2004) investigated high and low expressers of cyp19a2 (the brain form) and found that it did not correlate with sex in adult zebrafish  Furthermore, a recent study on cyp19a1b expression in brain of adult zebrafish found that the expression levels were in the same range in males and females . This is in accordance with the present study that indicates high brain expression of cyp19a1b in adults of both sexes. Due to the conflicting results regarding the involvement of cyp19a1b in gonadal sex differentiation, cyp19a1b is not included in the calculation of ratios between expected male and female genes in this study.
The coinciding low expression of sox9a, dmrt1, ar, fig α and cyp19a1a (in the high expresser groups) at 19–20 dph is just prior to the time of oocyte apoptosis in individuals developing into males. This could indicate that genes involved in the sex differentiation might be down-regulated at this time in development. Interestingly, it has previously been shown that the vas gene (which is a germ cell marker) was also down-regulated at 23 dpf . Furthermore, the genes expected to be highly expressed in females (fig α and cyp19a1a) both had increased expression levels from 22 dph. Also, two of the three genes expected to be highly expressed in males (ar and sox9a) have a peak in expression at 22 dph. The last gene expected to be highly expressed in males (dmrt1) does not have a significant peak in expression after the general down regulation at 19–20 dph; this could indicate that dmrt1 (which shows two significant peaks in expression before 19–20 dph, at 10 and 14 dph) is involved in gonadal sex differentiation upstream in the signalling cascade.
In conclusion, the current study investigated the expression pattern of six genes previously suggested to be involved in sex differentiation in individual zebrafish during the period of sex determination and differentiation. The three genes expected to be highest expressed in males (ar, sox9a and dmrt1) segregated into two groups of high and low expressers, respectively. These groups are putatively based on difference in genetic sex and might indicate a possible sex-related difference in expression profiles by the time of sex determination and gonadal differentiation. Thereby suggesting an important role for ar, sox9a and dmrt1 in controlling the male development process in zebrafish. Accordingly, the expression of the two genes expected to be highest expressed in females (fig α and cyp19a1a) also segregated fish into high and low expressers. Since the peaks in expression of these two genes are relatively late in the investigated period they are most likely involved in the gonadal sex differentiation of females and not in the sex determining process. The peak in expression of dmrt1, at 10 dph, suggests that dmrt1 might be upstream in the sex determining cascade, playing an important role in early gonadal sex differentiation of zebrafish males.
Trant JM, Gavasso S, Ackers J, Chung BC, Place AR: Developmental Expression of Cytochrome P450 Aromatase Genes (Cyp19a and Cyp19b) in Zebrafish Fry (Danio Rerio). Journal of Experimental Zoology. 2001, 290: 475-483. 10.1002/jez.1090.
Von Hofsten J, Olsson PE: Zebrafish Sex Determination and Differentiation: Involvement of Ftz-F1 Genes. Reproductive Biology and Endocrinology. 2005, 3: 63-73. 10.1186/1477-7827-3-63.
Wang XG, Bartfai R, Sleptsova-Freidrich I, Orban L: The Timing and Extent of 'juvenile Ovary' Phase Are Highly Variable During Zebrafish Testis Differentiation. Journal of Fish Biology. 2007, 70: 33-44. 10.1111/j.1095-8649.2007.01363.x.
Nanda I, Kondo M, Hornung U, Asakawa S, Winkler C, Shimizu A, Shan ZH, Haaf T, Shimizu N, Shima A, Schmid M, Schartl M: A Duplicated Copy of Dmrt1 in the Sex-Determining Region of the Y Chromosome of the Medaka, Oryzias Latipes. Proc Natl Acad Sci USA. 99: 11778-11783. 10.1073/pnas.182314699.
Paul-Prasanth B, Matsuda M, Lau EL, Suzuki A, Sakai F, Kobayashi T, Nagahama Y: Knock-Down of Dmy Initiates Female Pathway in the Genetic Male Medaka, Orydas Latipes. Biochemical and Biophysical Research Communications. 2006, 351: 815-819. 10.1016/j.bbrc.2006.10.095.
Orlando EF, Guillette LJ: Sexual dimorphic responses in wildlife exposed to endocrine disrupting chemicals. Environ Res. 2007, 104: 163-173. 10.1016/j.envres.2006.06.002.
Matsuda M: Sex Determination in Fish: Lessons From the Sex-Determining Gene of the Teleost Medaka, Oryzias Latipes. Development Growth & Differentiation. 2003, 45: 397-403. 10.1111/j.1440-169X.2003.00716.x.
Guo YQ, Cheng HH, Huang X, Gao S, Yu HS, Zhou RJ: Gene Structure, Multiple Alternative Splicing, and Expression in Gonads of Zebrafish Dmrt1. Biochemical and Biophysical Research Communications. 2005, 330: 950-957. 10.1016/j.bbrc.2005.03.066.
Raymond CS, Murphy MW, O'sullivan MG, Bardwell VJ, Zarkower D: Dmrt1, a Gene Related to Worm and Fly Sexual Regulators, Is Required for Mammalian Testis Differentiation. Genes & Development. 2000, 14: 2587-2595. 10.1101/gad.834100.
Maldonado LCT, Piedra AL, Mendoza NM, Valencia AM, Martinez AM, Larios HM: Expression Profiles of Dax1, Dmrt1, and Sox9 During Temperature Sex Determination in Gonads of the Sea Turtle Lepidochelys Olivacea. General and Comparative Endocrinology. 2002, 129: 20-26. 10.1016/S0016-6480(02)00511-7.
Shibata K, Takase M, Nakamura M: The Dmrt1 Expression in Sex-Reversed Gonads of Amphibians. General and Comparative Endocrinology. 2002, 127: 232-241. 10.1016/S0016-6480(02)00039-4.
Nanda I, Zend-Ajusch E, Shan Z, Grutzner F, Schartl M, Burt DW, Koehler M, Fowler VM, Goodwin G, Schneider WJ, Mizuno S, Dechant G, Haaf T, Schmid M: Conserved Synteny Between the Chicken Z Sex Chromosome and Human Chromosome 9 Includes the Male Regulatory Gene Dmrt1: a Comparative (Re)View on Avian Sex Determination. Cytogenetics and Cell Genetics. 2000, 89: 67-78. 10.1159/000015567.
Berta P, Hawkins JR, Sinclair AH, Taylor A, Griffiths BL, Goodfellow PN, Fellous M: Genetic-Evidence Equating Sry and the Testis-Determining Factor. Nature. 1990, 348: 448-450. 10.1038/348448A0.
Koopman P, Gubbay J, Vivian N, Goodfellow P, Lovellbadge R: Male Development of Chromosomally Female Mice Transgenic for Sry. Nature. 1991, 351: 117-121. 10.1038/351117a0.
Hacker A, Capel B, Goodfellow P, Lovellbadge R: Expression of Sry, the Mouse Sex-Determining Gene. Development. 1995, 121: 1603-1614.
Chiang EFL, Pai CI, Wyatt M, Yan YL, Postlethwait J, Chung BC: Two Sox9 Genes on Duplicated Zebrafish Chromosomes: Expression of Similar Transcription Activators in Distinct Sites. Developmental Biology. 2001, 231: 149-163. 10.1006/dbio.2000.0129.
Zhou RJ, Cheng HH, Zhang QY, Guo YQ, Cooper RK, Tiersch TR: Sry-Related Genes in the Genome of the Rice Field Eel (Monopterus Albus). Genetics Selection Evolution. 2002, 34: 129-137. 10.1051/gse:2001008.
Borg B: Androgens in Teleost Fishes. Comparative Biochemistry and Physiology C-Pharmacology Toxicology & Endocrinology. 1994, 109: 219-245.
Ikeuchi T, Todo T, Kobayashi T, Nagahama Y: Cdna Cloning of a Novel Androgen Receptor Subtype. Journal of Biological Chemistry. 1999, 274: 25205-25209. 10.1074/jbc.274.36.25205.
Berrevoets CA, Umar A, Brinkmann AO: Antiandrogens: Selective Androgen Receptor Modulators. Molecular and Cellular Endocrinology. 2002, 198: 97-103. 10.1016/S0303-7207(02)00373-8.
Olsson PE, Berg AH, Von Hofsten J, Grahn B, Hellqvist A, Larsson A, Karlsson J, Modig C, Borg B, Thomas P: Molecular Cloning and Characterization of a Nuclear Androgen Receptor Activated by 11-Ketotestosterone. Reproductive Biology and Endocrinology. 2005, 3: 37-54. 10.1186/1477-7827-3-37.
Blazquez M, Piferrer F: Sea Bass (Dicentrarchus Labrax) Androgen Receptor: Cdna Cloning, Tissue-Specific Expression, and Mrna Levels During Early Development and Sex Differentiation. Molecular and Cellular Endocrinology. 2005, 237: 37-48. 10.1016/j.mce.2005.04.001.
Jorgensen A, Andersen O, Bjerregaard P, Rasmussen LJ: Identification and characterisation of an androgen receptor from zebrafish Danio rerio. Comp Biochem Physiol C Toxicol Pharmacol. 2007, 146: 561-568. 10.1016/j.cbpc.2007.07.002.
Hossain MS, Larsson A, Scherbak N, Olsson PE, Orban L: Zebrafish Androgen Receptor: Isolation, Molecular, and Biochemical Characterization. Biol Reprod. 2008, 78: 361-369. 10.1095/biolreprod.107.062018.
Onichtchouk D, Aduroja K, Belting HG, Gnugge L, Driever W: Transgene Driving Gfp Expression From the Promoter of the Zona Pellucida Gene Zpc Is Expressed in Oocytes and Provides Marker for Gonad Differentiation in an Early Zebrafish. Developmental Dynamics. 2003, 228: 393-404. 10.1002/dvdy.10392.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.
Goto-Kazeto R, Kight KE, Zohar Y, Place AR, Trant JM: Localization and Expression of Aromatase mRNA in Adult Zebrafish. General and Comparative Endocrinology. 2004, 139: 72-84. 10.1016/j.ygcen.2004.07.003.
Kallivretaki E, Eggen RIC, Neuhauss SCF, Kah O, Segner H: The zebrafish, Brain-Specific, Aromatase cyp19a2 Is Neither Expressed nor Distributed in a Sexually Dimorphic Manner During sexual Differentiation. Developmental Dynamics. 2007, 236: 3155-3166. 10.1002/dvdy.21344.
Nakamoto M, Suzuki A, Matsuda M, Nagahama Y, Shibata N: Testicular Type Sox9 Is Not Involved in Sex Determination but Might Be in the Development of Testicular Structures in the Medaka, Orydas Latipes. Biochemical and Biophysical Research Communications. 2005, 333: 729-736. 10.1016/j.bbrc.2005.05.158.
Moreno-Mendoza N, Harley VR, Merchant-Larios H: Differential Expression of Sox9 in Gonads of the Sea Turtle Lepidochelys Olivacea at Male- or Female-Promoting Temperatures. Journal of Experimental Zoology. 1999, 284: 705-710. 10.1002/(SICI)1097-010X(19991101)284:6<705::AID-JEZ12>3.0.CO;2-1.
Dasilva SM, Hacker A, Harley V, Goodfellow P, Swain A, Lovellbadge R: Sox9 Expression During Gonadal Development Implies a Conserved Role for the Gene in Testis Differentiation in Mammals and Birds. Nature Genetics. 1996, 14: 62-68. 10.1038/ng0996-62.
Harley VR, Clarkson MJ, Argentaro A: The Molecular Action and Regulation of the Testis-Determining Factors, Sry (Sex-Determining Region on the Y Chromosome) and Sox9 [Sry-Related High-Mobility Group (Hmg) Box 9]. Endocrine Reviews. 2003, 24: 466-487. 10.1210/er.2002-0025.
Rodriguez-Mari A, Yan YL, Bremiller RA, Wilson C, Canestro C, Postlethwait JH: Characterization and Expression Pattern of Zebrafish Anti-Mullerian Hormone (Amh) Relative to Sox9a, Sox9b, and Cyp19ala, During Gonad Development. Gene Expression Patterns. 2005, 5: 655-667. 10.1016/j.modgep.2005.02.008.
Koopman P, Loffler KA: Sex Determination: the Fishy Tale of Dmrt1. Current Biology. 2003, 13: R177-R179. 10.1016/S0960-9822(03)00117-9.
Kurokawa H, Saito D, Nakamura S, Katoh-Fukui Y, Ohta K, Baba T, Morohashi K, Tanaka M: Germ cells are essential for sexual dimorphism in the medaka gonad. Proc Natl Acad Sci USA. 2007, 104: 16958-16963. 10.1073/pnas.0609932104.
Kondo M, Nanda I, Hornung U, Asakawa S, Shimizu N, Mitani H, Schmid M, Shima A, Schartl M: Absence of the Candidate Male Sex-Determining Gene Dmrt1b(Y) of Medaka From Other Fish Species. Current Biology. 2003, 13: 416-420. 10.1016/S0960-9822(03)00125-8.
Kobayashi T, Matsuda M, Kajiura-Kobayashi H, Suzuki A, Saito N, Nakamoto M, Shibata N, Nagahama Y: Two DM Domain Genes, Dmy and Dmrt1, Involved in Testicular Differentiation and Development in the Medaka, Oryzias Latipes. Developmental Dynamics. 2004, 231: 518-526. 10.1002/dvdy.20158.
Kanamori A, Yarnamura A, Koshiba S, Lee JS, Orlando EF, Hori H: Methyltestosterone Efficiently Induces Male Development in the Self-Fertilizing Hermaphrodite Fish, Kryptolebias Marmoratus. Genesis. 2006, 44: 495-503. 10.1002/dvg.20240.
Soyal SM, Amleh A, Dean J: Fig Alpha, a Germ Cell-Specific Transcription Factor Required for Ovarian Follicle Formation. Development. 2000, 127: 4645-4654.
Kim SJ, Ogasawara K, Park JG, Takemura A, Nakamura M: Sequence and Expression of Androgen Receptor and Estrogen Receptor Gene in the Sex Types of Protogynous Wrasse, Halichoeres Trimaculatus. General and Comparative Endocrinology. 2002, 127: 165-173. 10.1016/S0016-6480(02)00020-5.
He CL, Du JL, Lee YH, Huang YS, Nagahama Y, Chang CF: Differential Messenger Rna Transcription of Androgen Receptor and Estrogen Receptor in Gonad in Relation to the Sex Change in Protandrous Black Porgy, Acanthopagrus Schlegeli. Biology of Reproduction. 2003, 69: 455-461. 10.1095/biolreprod.102.015040.
Wang XG, Orban L: Anti-Mullerian Hormone and 11 Beta-Hydroxylase Show Reciprocal Expression to That of Aromatase in the Transforming Gonad of Zebrafish Males. Developmental Dynamics. 2007, 236: 1329-1338. 10.1002/dvdy.21129.
Andersen L, Goto-Kazeto R, Trant JM, Nash JP, Korsgaard B, Bjerregaard P: Short-Term Exposure to Low Concentrations of the Synthetic Androgen Methyltestosterone Affects Vitellogenin and Steroid Levels in Adult Male Zebrafish (Danio Rerio). Aquatic Toxicology. 2006, 76: 343-352. 10.1016/j.aquatox.2005.10.008.
Sawyer SJ, Gerstner KA, Callard GV: Real-Time Pcr Analysis of Cytochrome P450 Aromatase Expression in Zebrafish: Gene Specific Tissue Distribution, Sex Differences, Developmental Programming, and Estrogen Regulation. General and Comparative Endocrinology. 2006, 147: 108-117. 10.1016/j.ygcen.2005.12.010.
Krovel AV, Olsen LC: Sexual Dimorphic Expression Pattern of a Splice Variant of Zebrafish Vasa During Gonadal Development. Developmental Biology. 2004, 271: 190-197. 10.1016/j.ydbio.2004.04.004.
This work was supported by grant no. 272-06-0457 from the Danish National Science Research Council to O. Andersen and P. Bjerregaard.
The authors declare that they have no competing interests.
AJ conducted the RNA purification, quantitative qRT-PCR and statistical analysis. JEM setup the zebrafish and collected the fish. OA, LJR and PB devised the study, participated in the discussion of results and participated in the planning of experiments. AJ wrote the manuscript. All authors read and approved the final manuscript.