- Research article
- Open Access
Gene expression analysis at the onset of sex differentiation in turbot (Scophthalmus maximus)
© Robledo et al. 2015
- Received: 19 June 2015
- Accepted: 23 October 2015
- Published: 18 November 2015
Controlling sex ratios is essential for the aquaculture industry, especially in those species with sex dimorphism for relevant productive traits, hence the importance of knowing how the sexual phenotype is established in fish. Turbot, a very important fish for the aquaculture industry in Europe, shows one of the largest sexual growth dimorphisms amongst marine cultured species, being all-female stocks a desirable goal for the industry. Although important knowledge has been achieved on the genetic basis of sex determination (SD) in this species, the master SD gene remains unknown and precise information on gene expression at the critical stage of sex differentiation is lacking. In the present work, we examined the expression profiles of 29 relevant genes related to sex differentiation, from 60 up to 135 days post fertilization (dpf), when gonads are differentiating. We also considered the influence of three temperature regimes on sex differentiation.
The first sex-related differences in molecular markers could be observed at 90 days post fertilization (dpf) and so we have called that time the onset of sex differentiation. Three genes were the first to show differential expression between males and females and also allowed us to sex turbot accurately at the onset of sex differentiation (90 dpf): cyp19a1a, amh and vasa. The expression of genes related to primordial germ cells (vasa, gsdf, tdrd1) started to increase between 75–90 dpf and vasa and tdrd1 later presented higher expression in females (90-105 dpf). Two genes placed on the SD region of turbot (sox2, fxr1) did not show any expression pattern suggestive of a sex determining function. We also detected changes in the expression levels of several genes (ctnnb1, cyp11a, dmrt2 or sox6) depending on culture temperature.
Our results enabled us to identify the first sex-associated genetic cues (cyp19a1a, vasa and amh) at the initial stages of gonad development in turbot (90 dpf) and to accurately sex turbot at this age, establishing the correspondence between gene expression profiles and histological sex. Furthermore, we profiled several genes involved in sex differentiation and found specific temperature effects on their expression.
- Sex ratio
Sex was thought to have arisen in a single evolutive event in the last common ancestor of all eukaryotes, since sexual reproduction is almost universal and exclusive of this group . Considering its consequences over the lifespan of an organism and its influence on population demography, it is thought that the sex determination (SD) mechanism should be under strong selection forces . However, sex can be established by many different and fast-evolving mechanisms [3, 4], indicating that SD triggers have emerged several times throughout evolution . Within vertebrates, different sex determining systems have been described. In therian mammals, with a XX/XY chromosome system, sex depends on the presence of the Sry gene, a paralogue of SOX3, on the Y chromosome , while in birds (chicken) with a ZZ/ZW chromosome system, the DMRT1 gene with a double dosage is required for testis development . Also, in Xenopus laevis the DM-W gene, a paralogue of DMRT1, is responsible for SD . When dmY, belonging to the DM family of transcription factors like DMRT1 and DM-W, was found to be the sex determining gene (SDG) of the fish Oryzias latipes, a biased and recurrent recruitment of specific SDGs or families throughout evolution was suggested . However, later findings in fish do not seem to support this hypothesis.
Fish, with around 30,000 species , is the most diverse group of vertebrates and its study has broadened our knowledge on SD. Fish diversity is also reflected by the variety of reproductive strategies: unisexuality, different types of hermaphroditism and gonochorism; and also by the diversity of SD mechanisms . In the last years, an important effort has been made in order to identify the SDG in several model and aquaculture fish species. Different productive traits are sex-associated in farm fish such as growth rate, color, taste and flesh quality; hence, the interest of the industry in producing monosex stocks . Detailed information at gene level is available for only a limited number of fish species. Five different master SDGs have been identified so far: dmY /dmrt1by in Oryzias latipes and in O. curvinotus , gsdf in O. luzonensis , amhy in Odontesthes hatchery , amhr2 in Takifugu rubripes, T. pardalis and T. poecilonotus , and sdY in salmonid family . Recently, a distant cis-regulatory element of Sox3 necessary for male determination in O. dancena, a species with a XX/XY SD system, has also been identified , and dmrt1 has been suggested as the SDG in Cynoglossus semilaevis . However, little information is available, not only on the SDG, but also on the initial molecular pathways related to sexual differentiation.
Traditionally, SD has been related to the switching mechanism of a hierarchical genetic network that causes the activation of downstream genes involved in gonad differentiation (GD) leading to the differentiation of testes or ovaries . Thus, concerning whether the first difference between sexes is the expression of a gene/s or the strength of an environmental factor, SD can be genetic or environmental, although both ways can coexist [12, 21]. In the classical view of SD and GD, the downstream genes of the cascade were assumed to be highly conserved, and only the genes at the top of the cascade would change by gene duplication (and by the recruitment of a downstream gene) or by allelic diversification, establishing a new SD mechanism . Nowadays, the conservation of the downstream cascade has been questioned [22, 23] and a new view, which considers sex as a threshold phenotype in which both genetic and environmental factors can act alone or in combination and, importantly, at different times during the period of GD is gaining support [4, 24]. In this new view, SD encompasses not only the initial trigger, be it genetic, environmental or both, but also the whole GD process, and different factors such as cell proliferation and hormone levels would be involved in determining a threshold which would give rise to a testis or an ovary, thus fitting to a threshold quantitative trait [4, 12].
Turbot is one of the most important species cultured in Europe, being Galicia (North-West of Spain) the main production region since the eighties. Production and quality of farmed fish rely on a deep knowledge of biological functions, especially those related to reproduction, growth and disease resistance. In this context, the production of monosex stocks to exploit sex-associated dimorphisms related to productive traits, especially growth and sexual maturation, has been a sought in finfish aquaculture by different approaches [11, 25]. Turbot shows one of the strongest sexual growth dimorphisms amongst marine species and females can reach up to 50 % bigger size than males , thus industry is interested in the production of all-female populations. In the last years, an important effort has been devoted to understanding SD and GD in this species. Analysis on mitotic and meiotic chromosomes revealed the absence of an heteromorphic sex chromosome pair related to sex [27, 28]. The major SD region was located on the linkage group (LG) 5 at 2.6 cM of Sma-USC30 marker (R2 = 86.1 %) , but other minor sex-related quantitative trait locus (QTLs) were detected at LG6, LG8 and LG21 . In that study, a ZZ/ZW chromosome system was established in accordance with the sex ratios of progenies obtained from hormonally sex-reversed parents . Temperature also showed a minor influence on sex ratios in this species . Close to the sex-associated marker Sma-USC30 several candidate genes were identified (sox2, dnajc19, fxr1, atp11b, fkbp2 and dlg1), but eventually discarded because no association to sex was detected at the species level, so the SD gene remains unidentified in turbot . Considering the lack of information on the SD mechanism of turbot, we determined the expression patterns of a suite of 29 genes shown to be involved in GD in other species at the initial critical stages of sex differentiation using a large amount of fish and sampling times in turbot. In a previous study, we analyzed reproduction in this species through the use of oligo-microarrays, spanning a larger age period but in a lower number of samples and finding several genes involved in ovary or testis development (Ribas et al. submitted). In this work, our aim was to study gene expression at a very specific time point coinciding with the onset of sex differentiation. We also evaluated the effect of temperature to clarify whether it has a major role in turbot SD and its possible interaction with genetic factors. Our results enabled us to establish the correspondence between gene expression profiles and histological sex and to identify the first sex-associated genetic cues at the initial stages of gonad development in turbot.
Sampling and sexing
Length comparison between temperatures for each age group
Mean length ± SD (cm)
Percentage (%) / 18 °C
P value/18 °C
3.14 ± 0.27
3.44 ± 0.26
3.72 ± 0.18
3.63 ± 0.43
4.13 ± 0.45
4.73 ± 0.38
4.95 ± 0.44
5.90 ± 0.42
6.02 ± 0.52
7.21 ± 0.68
9.76 ± 0.51
10.75 ± 0.47
8.33 ± 0.37
11.26 ± 0.70
11.85 ± 0.68
9.57 ± 0.33
12.34 ± 0.36
13.28 ± 0.71
Co-localization of targeted genes with sex-related QTLs
Global expression patterns
PCA analysis on the 180 samples (Fig. 3b) revealed that histologically sexed males and females (≥105 dpf; referred to as males and females in the figure), could be in most cases clearly discriminated by their differential expression. However, younger individuals (60–90 dpf; only genetically sexed and referred to as genetic males and females in the figure), appeared fully overlapped in the PCA, likely because they are still undifferentiated. A 66 % prediction ellipse for each group is shown in Fig. 3b indicating that if new individuals were added to our analysis from a certain group, 66 % of them would expect to be placed inside the corresponding ellipse. Some of the genetically sexed individuals (60, 75 and 90 dpf) are also found in the ellipses drawn for phenotypic males or females which, also considering the previous heatmap results, indicates that turbot GD might start before 105 dpf. Interestingly, about half a dozen genetic males were included in the female circle. The arrows indicate the weight of each gene on the two first principal components. Among the analyzed genes there is a large group seemingly contributing to female differentiation (e.g., cyp19a1a, sox19, tdrd1, dact1), while the presumed male-related genes are fewer and not so markedly pointing towards male differentiation (sox9, amh, sox2, hh1). Also, as suggested in the heatmap, some genes are clearly related to undifferentiated individuals (sox6, fxr1, wnt4, hsp27, ptges).
Sex differences were first observed at 90 dpf, when cyp19a1a (FCF/M = 2.0) and vasa (FCF/M = 1.2) (Fig. 4) which were more expressed in females, and amh (FCM/F = 1.1) in males (Fig. 4). These three genes also increased their expression from 75 to 90 dpf, in both sexes (~5 cm), which is also observed for foxl2, tdrd1, gsdf (Fig. 4) and sox19 and rdh3 (Fig. 5). Gsdf expression increased dramatically even earlier, from 60 to 75 dpf (FC75/60 = 4.3, ~2.5 cm). Foxl2 did not show significant differences in expression in females until 120 dpf, mainly because its expression decreased slowly in males (Fig. 4). Tdrd1 pattern of expression resembled that of vasa, being more highly expressed in females at 105 dpf and onwards. Other genes increased their expression at 105 dpf (~9 cm length) in females, while their expression remained low or even decreased in males. Fold change values of females vs. males at 105 dpf were above 1 for sox19 (FCF/M = 3.8), dnmt1 (FCF/M = 3.1), dact1 (FCF/M = 2.4), rdh3 (FCF/M = 1.8), ctnnb1 (FCF/M = 1.3), sf1 (FCF/M = 1.1) and piwi2 (FCF/M = 2.3) (Fig. 5) and these differences increased from 105 to 135 dpf, when the FCs were between 2 and 7 for the seven genes. Fxr1 was also up-regulated in females at 105 dpf (FCF/M = 1.2) and showed a 47 % increase at 135 dpf (Fig. 6), but presented higher expression levels in undifferentiated individuals, irrespective of sex. A similar pattern was also observed for wnt4, dmrt2 and zar1 (Fig. 6), genes which at some point during sex differentiation, 105–135 dpf (above 10 cm length), showed a higher expression in females, but its expression decreased from 75 to 90 dpf (~5 cm, Fig. 6).
In contrast to that observed in females, in our study there were very few genes in the subset assayed whose expression were higher in males and, even in these cases, the differences between males and females were low. One of them was sox9a for which sex differences increased from 105 dpf (FCM/F = 1.2, ~11 cm) up to 135 dpf (FCM/F = 2.0) (Fig. 6). Another sox family gene, sox8, highly expressed in undifferentiated individuals in both sexes, was up-regulated in males at 105 dpf (FCM/F = 0.8), but the difference remained constant at 135 dpf (Fig. 6). Two additional genes, fshb and cyp11a, showed mean expression values slightly higher in males, but not significant (Additional file 3).
Some genes like ptges3, hh1, hsp27 or T4_30483 did not show sex differences (Additional file 3). Among these, a gene of the sox family, sox17, showed some groups of outliers whose expression was not explained either by sex or by length/age. Other two genes of this family, sox2 and sox6, did not present any clear expression pattern along development or by sex, and androgen receptor 1, ar1 did not show sex dimorphic expression either.
Temperature effects on gene expression
Finally, since 90 dpf is the first age where we find sex dimorphic expression and can be considered the onset of sex differentiation, we decided to test for temperature differences at this age which were independent of both length and sex, since most genes did not show dimorphic gene expression at this age. A single gene showed significant differences between temperatures at 90 dpf: ctnnb1 (Additional file 5). Ctnnb1 highest expression was observed at 18 °C, while its lowest expression was at 23 °C.
Recently, sex determination and differentiation have begun to be perceived as a single modular process rather than a cascade, where sex behaves in many instances as a threshold-like character [33, 34]. Under this view, gonad fate depends on several factors acting coordinately, including many genes and possibly also environmental variables. In this sense, and given the interest in obtaining monosex stocks in many finfish aquaculture species, understanding the different gene patterns during early sex differentiation, supposedly the moment where gonad fate can be more easily controlled or altered, is of great importance to manipulate sex determination in order to control sex ratios. In this study we have analyzed the expression of 29 genes previously connected to sex differentiation in other species, some of them studied in fish regarding sex differentiation for the first time. A total of 21 genes were found to show dimorphic expression at some point during early sex differentiation in turbot. The influence of temperature was also assessed, finding differences between temperatures for 10 genes. This study has broadened our knowledge of gene expression patterns during early sex determination in turbot in particular and in fish in general.
Early sex differentiation and primodial germ cells
Although morphological gonad differences between sexes were not detected before 105 dpf, the first molecular signs of sex differentiation were observed between 75 and 90 dpf (5–6 cm length), characterized by an expression increase of gsdf and tdrd1, but also by the expression increase and differential expression of cyp19a1a, amh and vasa, which allowed discriminating males (high amh expression) and females (high cyp19a1a and vasa expression) at 90 dpf. Vasa, tdrd1 and gsdf are genes related to primordial germ cell development. Tdrd1 proteins were detected in the primordial germ cells of zebrafish (Danio rerio) at 4 dpf and were involved in both oocyte and sperm development . Tdrd1 was found to interact with vasa, which is also a highly specific marker of germ cells [36–40] required for their development  and conserved along several invertebrate and vertebrate taxa . While vasa and tdrd1 are specific germ cell markers, gsdf is a teleost-specific member of the TGF-β superfamily which has been reported to be expressed in the somatic cells surrounding the primordial germ cells in rainbow trout and promoting their proliferation . Gsdf has shown higher levels of expression in testis in zebrafish, three-spot wrasse (Halichoeres trimaculatus) and coelacanth (Latimeria menadoensis) [44–46]. Furthermore, gsdf has been proposed as the male sex determining gene in Anoploma fimbria , and a copy named gsdf Y has been found to be the sex determinant gene in Oryzias luzonensis, with a higher expression in males 10 days after hatching . Gsdf has also been found to be directly up-regulated by dmy, male SDG of Oryzias latipes , and sox3, male SDG of Oryzias dancena  during the first stages of GD. Although gsdf does not seem to have such a male-like function in turbot, since it did not present sex dimorphic expression, it seems to be important for gonad development in both sexes since its expression greatly increased between 75 and 90 dpf. Even more, this is the only gene which showed a significant expression increase from 60 to 75 dpf in our study, consistent with a function as germ cell inductor since its expression increase precedes that of vasa or tdrd1.
The expression patterns of vasa, tdrd1 and gsdf in turbot suggested that primordial germ cells start proliferating between 75 and 90 dpf in both sexes, either faster in females or suffering a certain delay in males, as suggested by vasa and tdrd1 expression levels. Another germ-line specific gene, piwi2, was investigated in our study, but it showed a different expression pattern and its induction was limited to females and delayed until 105 dpf. In zebrafish, piwi2 has been found to play a crucial role in meiosis  and perhaps its different pattern of expression in turbot may be related to the onset of meiosis in female germ cells.
The amount of primordial cells is recognized as one of the initial differences between male and female gonads in some fish species like zebrafish , medaka  and stickleback (Gasterosteus aculeatus) , although not in others such as loach (Misgurnus anguillicaudatus)  and goldfish (Carassius auratus) . Further, germ cell proliferation has been found to be associated with SD in several fish species. In medaka, germ cell proliferation is inhibited in males when the sex determining gene, dmY, is expressed at the SD stage before testis differentiation . When dmy is not active in XY embryos, germ cells proliferate and enter meiosis as in XX embryos. Surprisingly, amh does not present a sex dimorphic expression in medaka during GD . Although the Müllerian ducts are not present in modern teleosts , amh orthologs have been described in fish and characterized as gonad specific key factors for male sex differentiation [55, 57, 58]. Yet, even if amh does not present sex dimorphic expression in medaka, it has been found to control germ cell proliferation in this species, and mutations on its receptor (amhrII) lead to excessive proliferation of germ cells which caused male-to-female sex reversal , although phenotypic female XY gonads still expressed dmY. Furthermore, if these amhrII mutants are depleted of germ cells, testis development takes place . Thus, amh seems to be a repressor of germ cell proliferation in medaka and to play a major role in sex determination. This is also the case in fugu (Fugu rubripes), where a single SNP in the coding region of amhrII is likely responsible for SD . This SNP encodes a protein with a reduced function and is fixed in females, which are not sensitive to amh . Also, in the Patagonian pejerrey, a copy of amh, amhY, has been found to be the sex determining gene and its action has been suggested to regulate germ cell proliferation, being upstream to the autosomal amh and relegating the former to a function in testicular maturation and/or spermatogenesis .
Interestingly, in the female gonad of zebrafish the maintenance of cyp19a1a expression, but not its activation, has been related to the presence of the primordial germ cells , pointing towards a model where amh is responsible for the control of germ cell proliferation while germ cells aid to maintain cyp19a1a expression. Amh and cyp19a1a are genes involved in sex differentiation across all vertebrate taxa [35, 61], and have been reported as male and female like genes, respectively, in several fish species [62, 63]. In our experiment, cyp19a1a showed higher expression in females and amh in males at 90 dpf. A threshold expression of amh could be controlling sexual fate. If amh does not reach the required expression level, primordial germ cells will proliferate and maintain cyp19a1a levels while amh levels decrease. On the contrary, if amh expression reaches a certain threshold, germ cells stop proliferating and cyp19a1a expression decreases. Although further data is required, this hypothesis seems to be consistent with the findings in turbot and other fish species, and could be a common mechanism controlling the balance between male and female gonad differentiation in fish.
Female sex differentiation
Turbot GD is in progress at 105–135 dpf, either towards males or females, and sex was easily identified by cyp19a1a expression alone at these developmental stages. Several female up-regulated genes were detected during this period of GD (foxl2, vasa, tdrd1, sox19, dnmt1, dact1, rdh3), while male gene expression pattern was very similar to undifferentiated fish for the assayed genes, excluding a few classical male-like genes (amh, sox9a, sox8). Foxl2 is a transcription factor that activates cyp19a1a transcription  and both genes are strictly co-expressed in mammals . The expression pattern of Foxl2 was similar to that previously described for cyp19a1a, vasa, tdrd1 and gsdf. Its expression in turbot is consistent with an activation of cyp19a1a; however, the later decrease of cyp19a1a while foxl2 expression is still high in males suggests other roles for foxl2 not related to cyp19a1a activation at early stages of development. Foxl2 expression has also been described in the male gonad of tilapia (Oerochromis niloticus), southern catfish (Silurus meridionalis) and goldfish [54, 64, 66]. Also, cells with cyp19a1a expression without foxl2 expression have been reported in medaka . Although the authors did not exclude that those cells had earlier foxl2 expression, at least foxl2 did not seem to be essential for cyp19a1a expression maintenance . Foxl2 seems to be involved in the initiation of the female gonad development cascade through the activation of cyp19a1a expression in fish, as previously suggested , however foxl2 up-regulation alone is not sufficient for the activation of cyp19a1a and it might have additional roles important for gonad development both in males and females.
Ctnnb1 is another female-like gene with a critical function in female differentiation in mammals, antagonizing sox9 and blocking testis development, thus promoting ovarian development . Ctnnb1 is the key downstream component of the canonical wnt signaling pathway and our results are in agreement with a conserved female function of this pathway, which has also been shown in zebrafish  and rainbow trout . However, in turbot, the wnt pathway seems to be at least partially independent of wnt4. The highest expression of this gene was found at undifferentiated stages, consistent with a role in early gonad development in both sexes also observed in mammals, although with a different function since in mammals it is responsible for the development of Müllerian ducts . Wnt4 expression dropped at 90 dpf irrespective of sex, and later was found scarcely up-regulated in females. Wnt4 is a key female gene in mammals which antagonizes fgf9 and down-regulates sox9 expression . Wnt4 not so clear female-pattern during GD in turbot is likely related to the absence of fgf9 in teleosts . No dimorphic wnt4 expression has been observed in rainbow trout (Oncorhynchus mykiss) , zebrafish  or the more distant Rana rugosa . Wnt4 does not show a conserved function in female SD throughout evolution, and the results in our study suggest that it is not involved in the expression of ctnnb1 in the female gonad development and therefore other WNT proteins could be responsible for activating the wnt pathway, which seems to have a conserved female prominent function.
Six genes, dnmt1, rdh3, sox19, dact1, ctnnb1 and sf1 showed a similar increasing expression pattern in females from 105 to 135 dpf. These genes showed high pair-wise correlation values and also the highest negative correlations with amh amongst all the assayed genes. Hence, amh down-regulation may be required for the activation of genes involved in female gonad development. The coactivation of these genes in female differentiation does not seem to be related to a specific pathway considering their functional diversity, but with the activation of several concomitant pathways at the beginning of ovarian development, and their up-regulation is possibly connected to the down-regulation of amh.
Besides amh, we found two up-regulated genes in males from 105 dpf onwards: sox9a and sox8. Furthermore, the heatmap and network analyses suggested that cyp11a and fshb are related to male development as well, but they did not display a dimorphic expression pattern at any stage. Sox9a is an essential player in sex differentiation and its male-like nature seems to be rather conserved throughout evolution, but its relevance seems to be variable [75, 76]. In mammals, this gene is directly activated by sry and is responsible and sufficient for fating the male gonad , also activating amh transcription . In birds, sox9 is co-expressed with amh and induced by the SD gene dmrt1 . However, in medaka, sox9 is not required for testis development . In turbot, sox9a dimorphic expression is found later than that of amh, and its expression is more stable along the assayed stages, suggesting a less important role in GD in this species.
Genes in the main SD region
Fxr1 is located in the turbot genome very close to Sma-USCE30 , the highest sex-associated marker in this species, within the main SD region at LG5 , thus representing a potential SD candidate gene. In this experiment, fxr1 is highly expressed before the first GD signs and precedes the expression of cyp19a1a and amh, although at this time its expression is not sexually dimorphic. Fxr1 is a RNA-binding protein and is an autosomal paralogue of fmrp (fragil X mental retardation 1), important for normal female reproductive function and cognition development in humans. Further, fxr1 It has been related to female gametogenesis in pigs  and Xenopus laevis . In turbot, fxr1 position in the functional network, connecting male and female clusters, and its higher expression at low temperatures, suggest that further analysis on this gene in younger individuals looking for a putative dimorphic pattern between sexes or other type of dimorphic cue are desirable. On the contrary sox2, also located in the main SD region of turbot , did not present a dimorphic expression pattern and also showed a relatively steady expression along all the assayed stages, suggesting no role of this gene in turbot SD.
A higher proportion of phenotypic females were obtained at 15 and 18 °C, however sex ratio differences between temperatures were not significant. This effect was not due to genetic segregation distortion, since genotyping of the SD marker strongly suggested that some genetic males developed as females. Although proper male-to-female sex reversal at low temperatures has not been described in fish , it has been described in reptiles, with genetic males developing as females [84, 85]. The influence of temperature on sex ratio in turbot was previously reported , but it was family-dependent and not always in the same way: in this study, two families presented a higher proportion of females than expected at 23 °C while another family presented more females at 15 °C. In any case, if temperature effects exist in turbot, they seem to be limited and family dependent, however further work should be done to evaluate this issue. Another possible explanation for genetic males developing as phenotypic females is the existence of a family-specific genetic interaction between a secondary SD QTL and the main SD QTL in LG5, which could be responsible for genetic males developing as phenotypic females. Independently of the mechanism, families with higher proportions of females could be very interesting for aquaculture, since it would lead to higher growth rates given the important sex size dimorphism of turbot.
Nevertheless, temperature effects on gene expression are very important to explain possible sex ratio shifts and other effects on sex differentiation. We found temperature effects on the expression of some of the genes assayed: amh, ctnnb1, cyp11a, dmrt2, piwi2, sf1, sox2, sox6, sox9a and sox17. Among these, only sox2 temperature effects were sex independent. This gene showed higher expression both at 15 and 23 °C, which may indicate some kind of stress response. Sox2 is a transcription factor regulating several genes and it is also involved in the maintenance of stem-cell identity .
Among the genes which presented sex dependent temperature effects, ctnnb1, piwi2, sf1 and sox6 showed higher expression at low temperatures in males. Among these, ctnnb1 is remarkable because it not only occupies a downstream position in the wnt pathway , which needs to be up-regulated for the development of an ovary in zebrafish , but also because it showed temperature-dependent expression differences in our study at 90 dpf, the onset of SD. Genes found upstream in the wnt signalling pathway are likely regulated by temperature and responsible for this increase in ctnnb1 expression in turbot. Consistently with our results, elevated ctnnb1 expression has been reported connected to low temperatures in rats  and tilapia , albeit in other tissues. A recent study in oyster found a biased sex ratio towards females related to higher ctnnb1 expression at lower rearing temperatures . Thus, ctnnb1 and the wnt pathway are good candidates for future studies aimed at investigating temperature effects on sex ratios in turbot and possibly other fish.
Other two genes showed the opposite pattern in females, with higher expression at higher temperatures: dmrt2 and sox17. Dmrt2 was upregulated during gonad development and also expressed in germ cells in the swamp eel (Monopterus albus) , while in Rana rugosa it was expressed in the developing gonad during SD without any dimorphic pattern, suggesting a function both in testicular and ovarian differentiation . Sox17 has been associated with ovarian development in Dicentrarchus labrax , although it did not present dimorphic expression in turbot and so, apparently, it is not related to female differentiation in this species. Yet, this gene seemed to have a peak of expression at 90 dpf in both sexes, so it could have some function in early gonad development in turbot for both males and females. In the swamp eel, sox17 was also expressed both in testis and ovary .
Temperature effects are gene and sex specific. Given the labile nature of SD in fish related to specific morphogenetic thresholds, several genes could be responsible for sex ratio shifts. As seen in this study, several genes involved in sex show expression differences due to temperature, and so these genes are potential candidates for sex ratio alterations.
Turbot sex differentiation is ongoing at 90 dpf and sex can be distinguished by the expression levels of three genes when fish are 5–6 cm length: cyp19a1a, amh and vasa; while later females are easily discriminated by the expression of cyp19a1a. The first molecular signs of sex differentiation are the dimorphic expression of these three genes and an increase in the expression of vasa, gsdf and tdrd1, connected with primordial germ cells, suggesting their proliferation from 75 to 90 dpf and an important role in sex differentiation. The primary sex determining gene of turbot remains unknown, since neither sox2 nor fxr1, genes located in the main SD region of turbot, showed expression patterns clearly consistent with this role. Our data suggest that female development has more complex machinery and is strongly regulated, suggesting the involvement of both methylation and splicing mechanisms. Furthermore, we have observed that temperature affects the expression of several genes and suggest that the Wnt/β-catenin pathway could be a likely candidate to explain possible temperature-induced sex ratio shifts.
The present study, analyzing the expression pattern of many genes related to sex differentiation, has revealed that turbot sex differentiation is a complex process with many factors involved. These results are more compatible with a view of sex determination as a network where the activation or repression of several genes can affect gonad fate. This view of sex determination as a threshold character could help us to understand temperature effects during sex differentiation.
Rearing conditions and sampling
Turbot fertilized eggs were obtained at the IEO (Instituto Oceanográfico de Vigo, Spain) from a stock of wild spawners obtained at Celeiro (Galicia, North-West of Spain) and adapted to captivity. Turbot eggs were incubated in a 140 L cylindrical-conical tank at a temperature of 13–14 °C until hatching (5 days post fertilization, dpf). Hatched larvae were transferred to a 6 fiberglass cylindrical tanks (500 L) with an initial density of 30 larvae/L. The tanks were divided in three groups (2 tanks/group) and temperature was gradually adjusted to 15 °C (cold group), 18 °C (ambient group) and 23 °C (warm group). Water temperature was monitored 6–8 times per day. After metamorphosis the fish were transferred to 6 flat-bottom tanks and the temperatures maintained until the end of the experiment. Gonad samples were taken at 60, 75, 90, 105, 120 and 135 dpf. This period was chosen based on a preliminary analysis on the expression GD key genes (cyp19a1a, amh, sox9a, vasa, foxl2) every five days from 5 dpf, which did not show any expression changes before 90 dpf. We started our sampling 30 days earlier trying not to miss the expression of a possible sex determining gene. Ten fish per temperature and age were sampled and gonads dissected as accurately as possible considering the size of the fish. A total of 180 samples were used in this study: 6 ages × 3 temperatures/age × 10 fish for each age-temperature combination. In fish of 105 dpf and older, gonads were hemi dissected, one used for quantitative PCR (qPCR) and the other one for histological sexing. Animals were treated according to the Directive 2010/63/UE of the European Parliament and of the Council of 22 September 2010 on the protection of animals used for experimentation and other scientific purposes. All experimental protocols were approved by the Institutional Animal Care and Use Committee of the University of Santiago de Compostela (Spain).
Sexing by histology and molecular markers
Samples for histological analysis were kept in 4 % paraformaldehyde buffer overnight, rinsed with PBS the next day and stored in 70 % ethanol until further analysis. Samples were dehydrated and embedded in paraffin, cut at 7 μm thick and stained with hematoxylin-eosin for the determination of phenotypic sex. Additionally, all samples were genetically sexed using the SmaUSC-E30 marker, which demonstrated a ~98 % accuracy for the identification of genetic sex in turbot . To establish the association between sex and alleles at this marker, parents and grandparents of each family were genotyped, and the expected genotypes of male and female offspring obtained following .
RNA isolation and cDNA synthesis
Upon dissection, samples for qPCR were immediately embedded in RNAlater for preservation (Qiagen, Valencia, CA). Total RNA was extracted by homogenization in TRIzol (Invitrogen, Paisley, UK) following the manufacturer’s protocol. Extracted RNA was treated with RNase-free Recombinant DNaseI (Roche Diagnostics, Mannheim, DE) and RNA concentration was assessed by spectrophotometry and its quality checked using an Agilent 2100 bionalyzer (Agilent Technologies, Santa Clara, US). RNA (1.2 μg) was reverse transcribed by random primers using Affinity Script Multiple Temperature cDNA Synthesis Kit (Agilent Technologies) following the manufacturer’s protocol and then diluted 1:2 with nuclease-free water.
The 29 target genes were selected by: i) their importance for GD in other fish species; ii) previous data from our group in turbot [82, 95]; and iii) previous results from Ribas et al. (submitted) (Additional file 2). The 29 target genes are: amh, ar1, ctnnb1, cyp11a, cyp19a1a, dact1, dmrt2, dnmt1, fig-h, foxl2, fshb, fxr1, gsdf, hh1, hsp27, lhx8, piwi2, ptges3, rdh3, sf1, sox2, sox6, sox8, sox9a, sox17, sox19, tdrd1, vasa, wnt4 and zar1. (GenBank NCBI database accession numbers available in Additional file 6).
qPCR was performed on a Stratagene Mx3005P thermocycler (Agilent Technologies) using Brilliant III Ultra-Fast SYBR Green QPCR Master Mix in a final volume of 12.5 μL following the manufacturer’s protocol with 1 μL of cDNA per reaction. Specific primers for targeted genes were designed using Primer3  from sequences obtained from the turbot EST database enriched with sex differentiation-related organs (gonad and brain ). When possible, primers were designed spanning different exons (Additional file 6). Primer concentration was 300 nM and each sample was run in duplicate. The cycling parameters were: 50 °C for 2 min, 95 °C for 10 min, followed by 40 cycles of amplification at 95 °C for 15 s and 60 °C for 1 min. After amplification, a dissociation step was performed raising the temperature from 65 to 95 °C to create a melting curve and ensure the presence of a single amplification product. Specificity for each primer pair was also confirmed by PCR product sequencing. In every PCR plate, non-template controls were included to confirm the absence of contamination. In addition, the same three samples were run in triplicate in every plate in order to correct inter-assay variation. qPCR data were obtained by the MxPro software (Agilent Technologies) and quantification cycle values (Cq) calculated for each replicate and then averaged to obtain the final Cq value. Three reference genes (ubq, rps4, rpl17) were used for normalization and LinRegPCR software  was used for efficiency determination following the recommendations in . qPCR was performed in all the 180 samples for every gene. Samples with missing Cq values or inconsistencies between replicates (Cq difference > 1 cycle) were removed. Raw Cq values were transformed to the final fold difference values (FD) following the equations present in . Briefly, Cq values were normalized using the reference genes, efficiency corrected, log transformed and finally mean centered to obtain mean centered fold change values which were used for statistical analysis. All expression data is in Additional file 7, presented together with all the data collected for every sample.
Statistical analyses were performed using R (version 3.0.2) . Pearson correlations for the heatmap were obtained using the “cor” function. Principal component analysis (PCA) was computed by the “prcomp” function. Length and gene expression differences between sexes and ages were checked by Mann–Whitney tests (P < 0.05) since our data mostly did not conform to a normal distribution. Discriminant analysis was performed using the “lda” function on the “MASS” package . A Chi-square test was performed to assess if sex ratio difference between temperatures were significative. Multiple regression (p < 0.05) was used to assess temperature effects on gene expression, introducing temperature and length in the model. The gvlma function of the gvlma R package was used to check if our dataset met the assumptions of the multiple regression. Furthermore, we performed two additional tests for every temperature significant effect on gene expression: i) a moderation analysis, to check if length was modulated by the temperature, a temperature-length interaction term was added to our model checking if the new model improved the previous one; and ii) a mediation analysis by Sobel test, to explore if the detected temperature effect on gene expression is partially or fully explained by size differences between individuals.
Co-localization of targeted genes with sex-related QTLs
Several SD related QTLs were previously reported in turbot , and therefore, we considered relevant to establish the mapping position of the targeted genes regarding these QTLs in the last turbot map . For this, we established the correspondence between the turbot linkage groups and the scaffolds of the recently sequenced turbot genome (Figueras et al., unpublished) using the mapped markers and their sequences. Target gene sequences were located in the turbot genome using local blast  and then placed in the linkage map using the correspondence between linkage groups and scaffolds as far as accurately depending on the availability of markers in the vicinity.
Weighted correlation network analysis
Weighted correlation network analysis was performed in R (version 3.0.2)  using the WGCNA package  following the author’s tutorial. Co-expression networks were built for our genes and Cytoscape 3.0.2 was used to visualize the network . This allowed us to obtain information about the functional relationships between the target genes.
The data set supporting the results of this article (expression values for every gene in each sample and the data collected for every sample) is included in the Additional file 7.
The study was supported by a project from the Spanish Ministerio de Ciencia e Innovación (AGL2010-22326-C02-01 and CDS2007-0002). Diego Robledo was supported by a FPU fellowship from the Ministerio de Educación Cultura y Deporte of the Spanish Government. Laia Ribas was supported by Aquagenomics and EpigenAqua postdoctoral contract. We thank Lucía Insua, Susana Sánchez and Sonia Gómez for technical assistance.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Javaux EJ, Knoll AH, Walter MR. Morphological and ecological complexity in early eukaryotic ecosystems. Nature. 2001;412:66–9.PubMedView ArticleGoogle Scholar
- van Doorn GS. Evolutionary transitions between sex-determining mechanisms: a review of theory. Sex Dev. 2014;8(1-3):7–19.PubMedView ArticleGoogle Scholar
- Bull JJ. Evolution of sex determining mechanisms. Menlo Park: Benjamin/Cummings Publishing Company; 1983.Google Scholar
- Uller T, Helanterä. From the origin of sex-determining factors to the evolution of sex-determining systems. Q Rev Biol. 2011;86(3):163–80.PubMedView ArticleGoogle Scholar
- Graves JAM. Weird animal genomes and the evolution of vertebrate sex and sex chromosomes. Annu Rev Genet. 2008;42:565–86.View ArticleGoogle Scholar
- Sinclair AH, Berta P, Palmer MS, Hawkins JR, Griffiths BL, Smith MJ, et al. A gene from the human sex-determining region encodes a protein with homology to a conserved DNA-binding motif. Nature. 1990;346:240–4.PubMedView ArticleGoogle Scholar
- Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, et al. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature. 2009;461(7261):267–71.PubMedView ArticleGoogle Scholar
- Yoshimoto S, Okada E, Umemoto H, Tamura K, Uno Y, Nishida-Umehara C, et al. A W-linked DM-domain gene, DM-W, participates in primary ovary development in Xenopus laevis. Proc Nat Acad Sci U S A. 2008;105(7):2469–74.View ArticleGoogle Scholar
- Graves JAM, Peichel CL. Are homologies in vertebrate sex due to shared ancestry or limited options? Genome Biology. 2010;11(4):205.View ArticleGoogle Scholar
- Nelson JS. Fishes of the world. 4th ed. New York: Wiley-Interscience; 2006.Google Scholar
- Devlin RH, Nagahama Y. Sex determination and sex differentiation in fish: an overview of genetic, physiological, and environmental influences. Aquaculture. 2002;208(3-4):191–364.View ArticleGoogle Scholar
- Martínez P, Viñas AM, Sánchez L, Díaz N, Ribas L, Piferrer F. Genetic architecture of sex determination in fish: applications to sex ratio control in aquaculture. Frontiers in Genetics. 2014;5:340.PubMed CentralPubMedGoogle Scholar
- Matsuda M, Nagahama Y, Shinomiya A, Sato T, Matsuda C, Kobayashi T, et al. DMY is a Y-specific DM-domain gene required for male development in the medaka fish. Nature. 2002;417(6888):559–63.PubMedView ArticleGoogle Scholar
- Myosho T, Otake H, Masuyama H, Matsuda M, Kuroki Y, Fujiyama A, et al. Tracing the emergence of a novel sex-determining gene in medaka, Oryzias luzonensis. Genetics. 2012;191:163–70.PubMed CentralPubMedView ArticleGoogle Scholar
- Hattori RS, Murai Y, Oura M, Masuda S, Majhi SK, Sakamoto T, et al. A Y-linked anti-Müllerian hormone duplication takes over a critical role in sex determination. Proc Natl Acad Sci U S A. 2012;109:2955–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Kamiya T, Kai W, Tasumi S, Oka A, Matsunaga T, Mizuno N, et al. A Trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu). PLoS Genetics. 2012;8(7):e1002798.PubMed CentralPubMedView ArticleGoogle Scholar
- Yano A, Guyomard R, Nicol B, Jouanno E, Quillet E, Klopp C, et al. An immune-related gene evolved into the master sex-determining gene in rainbow trout, Oncorhynchus mykiss. Curr Biol. 2012;22(15):1423–8.PubMedView ArticleGoogle Scholar
- Takehana Y, Matsuda M, Myosho T, Suster ML, Kawakami K, Shin-I T, et al. Co-option of Sox3 as the male-determining factor on the Y chromosome in the fish Oryzias dancena. Nat commun. 2014;5:4157.PubMedView ArticleGoogle Scholar
- Chen S, Zhang G, Shao C, Huang Q, Liu G, Zhang P, et al. Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. Nat Genet. 2014;46(3):253–60.PubMedView ArticleGoogle Scholar
- Schartl M. Sex chromosome evolution in non-mammalian vertebrates. Curr Opin Genet Dev. 2004;14(6):634–41.PubMedView ArticleGoogle Scholar
- Penman DJ, Piferrer F. Fish gonadogenesis. Part I: genetic and environmental mechanisms of sex determination. Rev Fish Sci. 2008;16:14–32.View ArticleGoogle Scholar
- Böhne A, Heule C, Boileau N, Salzburger W. Expression and sequence evolution of aromatase cyp19a1 and other sexual development genes in East African cichlid fishes. Mol Biol Evol. 2013;30(10):2268–85.PubMed CentralPubMedView ArticleGoogle Scholar
- Herpin A, Schindler D, Kraiss A, Hornung U, Winkler C, Schartl M. Inhibition of primordial germ cell proliferation by the medaka male determining gene Dmrt1bY. BMC Dev Biol. 2007;7:99.PubMed CentralPubMedView ArticleGoogle Scholar
- Heule C, Salzburger W, Böhne A. Genetics of sexual development: an evolutionary playground for fish. Genetics. 2014;196(3):579–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Piferrer F. Endocrine sex control strategies for the feminization of teleost fish. Aquaculture. 2001;197:229–81.View ArticleGoogle Scholar
- Imsland AK, Folkvord A, Grung GL, Stefansson SO, Taranger GL. Sexual dimorphism in growth and maturation of turbot, Scophthalmus maximus (Rafinesque, 1810). Aquaculture Research. 1997;28:101–14.View ArticleGoogle Scholar
- Bouza C, Sanchez L, Martinez P. Karyotipic characterization of turbot (Scophthalmus maximus) with conventional fluorochrome and restriction endonuclease-banding techniques. Mar Biol. 1994;120:609–13.View ArticleGoogle Scholar
- Cuñado N, Terrones J, Sánchez L, Martínez P, Santos JL. Sex-dependent synaptic behavior in triploid turbot, Scophthalmus maximus (Pisces: Scophathalmidae). Heredity. 2001;89:460–4.View ArticleGoogle Scholar
- Hermida M, Bouza C, Fernández C, Sciara AA, Rodríguez-Ramilo ST, Fernández J, et al. Compilation of mapping resources in turbot (Scophthalmus maximus): A new integrated consensus genetic map. Aquaculture. 2013;414-415:19–25.View ArticleGoogle Scholar
- Martínez P, Bouza C, Hermida M, Fernandez J, Toro MA, Vera M, et al. Identification of the major sex-determining region of turbot (Scophthalmus maximus). Genetics. 2009;183:1443–52.PubMed CentralPubMedView ArticleGoogle Scholar
- Haffray P, Lebegue E, Jeu S, Guennoc M, Guiguen Y, Baroiller JF, et al. Genetic determination and temperature effects on turbot Scophthalmus maximus sex differentiation: An investigation using steroid sex-inverted males and females. Aquaculture. 2009;294(1-2):30–6.View ArticleGoogle Scholar
- Taboada X, Hermida M, Pardo BG, Vera M, Piferrer F, Viñas A, et al. Fine mapping and evolution of the major sex determining region in turbot (Scophthalmus maximus). G3. 2014;4(10):1871–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Mittwoch U. Sex is a threshold dichotomy mimicking a single gene effect. Trends Genet. 2006;22:96–100.PubMedView ArticleGoogle Scholar
- Quinn AE, Sarre SD, Ezaz T, Graves JAM, Georges A. Evolutionary transitions between mechanisms of sex determination in vertebrates. Biol Lett. 2011;7:443–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Cutting A, Chue J, Smith CA. Just how conserved is vertebrate sex determination. Dev Dyn. 2013;242:380–7.PubMedView ArticleGoogle Scholar
- Huang H, Houwing S, Kaaij LJT, Meppelink A, Redl S, Gauci S, et al. Tdrd1 acts as a molecular scaffold for Piwi proteins and piRNA targets in zebrafish. EMBO. 2011;30:3298–308.View ArticleGoogle Scholar
- Lasko PF, Ashburner M. The product of the Drosophila gene vasa is very similar to eukaryotic initiation factor-4A. Nature. 1988;335:611–7.PubMedView ArticleGoogle Scholar
- Hay B, Jan LY, Jan YN. A protein component of Drosophila polar granules is encoded by vasa and has extensive sequence similarity to ATP-dependent helicases. Cell. 1988;55:577–87.PubMedView ArticleGoogle Scholar
- Komiya T, Itoh K, Ikenishi K, Furusawa M. Isolation and characterization of a novel gene of the DEAD box protein family which is specifically expressed in germ cells of Xenopus laevis. Dev Biol. 1994;162:354–63.PubMedView ArticleGoogle Scholar
- Castrillon DH, Quade BJ, Wang TY, Quigley C, Crum CP. The human VASA gene is specifically expressed in the germ cell lineage. Proc Natl Acad Sci USA. 2000;97:9585–90.PubMed CentralPubMedView ArticleGoogle Scholar
- Johnstone O, Lasko P. Interaction with elF5B is essential for Vasa function during development. Development. 2004;131:4167–78.PubMedView ArticleGoogle Scholar
- Ephrussi A, Lehmann R. Induction of germ cell formation by oskar. Nature. 1992;358:387–92.PubMedView ArticleGoogle Scholar
- Sawatari E, Shikina S, Takeuchi T, Yoshizaki G. A novel transforming growth factor-β superfamily member expressed in gonadal somatic cells enhances primordial germ cell and spermatogonial proliferation in rainbow trout (Oncorhynchus mykiss). Dev Biol. 2007;301:266–75.PubMedView ArticleGoogle Scholar
- Gautier A, Sohm F, Joly JS, Le Gac F, Lareyre JJ. The proximal promoter region of the zebrafish gsdf gene is sufficient to mimic the spatio-temporal expression pattern of the endogenous gene in Sertoli and granulose cells. Biol Reprod. 2011;85:1240–51.PubMedView ArticleGoogle Scholar
- Horiguchi R, Nozu R, Hirai T, Kobayashi Y, Nagahama Y, Nakamura M. Characterization of gonadal soma-derived factor expression during sex change in the protogynous wrasse, Halichoeres trimaculatus. Dev Dyn. 2013;242(4):388–99.PubMedView ArticleGoogle Scholar
- Forconi M, Canapa A, Barucca M, Biscotti MA, Capriglione T, Buonocore F, et al. Characterization of sex determination and sex differentiation genes in Latimeria. PLoS One. 2013;8(4):e56006.PubMed CentralPubMedView ArticleGoogle Scholar
- Rondeau EB, Messmer AM, Sanderson DS, Jantzen SG, von Schalburg KR, Minkley DR, et al. Genomics of sablefish (Anoploma fimbria): expressed genes, mitochondrial phylogeny, linkage map and identification of a putative sex gene. BMC Genomics. 2013;14:452.PubMed CentralPubMedView ArticleGoogle Scholar
- Shibata Y, Paul-Prasanth B, Suzuki A, Usami T, Nakamoto M, Matsuda M, et al. Expression of gonadal soma derived factor (GSDF) is spatially and temporally correlated with early testicular differentiation in medaka. Gene Expr Patterns. 2010;10:283–9.PubMedView ArticleGoogle Scholar
- Houwing S, Berezikov E, Ketting RF. Zili is required for germ cell differentiation and meiosis in zebrafish. EMBO J. 2008;27:2702–11.PubMed CentralPubMedView ArticleGoogle Scholar
- Siegfried KR, Nüsslein-Volhard C. Germ line control of female sex determination in zebrafish. Dev Biol. 2008;324:277–87.PubMedView ArticleGoogle Scholar
- Kurokawa H, Saito D, Nakamura S, Katoh-Fukui Y, Ohta K, Baba T, et al. Germ cells are essential for sexual dimorphism in the medaka gonad. Proc Natl Acad Sci USA. 2007;104:16958–63.PubMed CentralPubMedView ArticleGoogle Scholar
- Lewis ZR, McClellan MC, Postlethwait JH, Cresko WA, Kaplan RH. Female-specific increase in primordial germ cells marks sex differentiation in threespine stickleback (Gasterosteus aculeatus). J Morphol. 2008;269:909–21.PubMedView ArticleGoogle Scholar
- Fujimoto T, Nishimura T, Goto-Kazeto R, Kawakami Y, Yamaha E, Arai K. Sexual dimorphism of gonadal structure and gene expression in germ cell-deficient loach, a teleost fish. Proc Natl Acad Sci USA. 2010;107:17211–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Goto R, Saito T, Takeda T, Fujimoto T, Takagi M, Arai K, et al. Germ cells are not the primary factor for sexual fate determination in goldfish. Dev Biol. 2012;370:98–109.PubMedView ArticleGoogle Scholar
- Klüver N, Pfennig F, Pala I, Storch K, Schlieder M, Froschauer A, et al. Differential expression of anti-Müllerian hormone receptor type II (amhrII) in the teleost medaka. Dev Dyn. 2007;236:271–81.PubMedView ArticleGoogle Scholar
- Nagahama Y. The functional morphology of teleost gonads. In: Hoar WS, DJ Randall, Donaldson EM, editors. Fish Physiology vol 9 part A. New York: Academic Press; 1983. p. 223–264.Google Scholar
- Miura T, Miura C, Konda Y, Yamauchi K. Spermatogenesis-preventing substance in Japanese eel. Development. 2002;129:2689–97.PubMedGoogle Scholar
- Halm S, Rocha A, Miura T, Prat F, Zanuy S. Anti-Müllerian hormone (AMH/AMH) in the European sea bass: its gene structure, regulatory elements, and the expression of alternatively-spliced isoforms. Gene. 2007;388:148–58.PubMedView ArticleGoogle Scholar
- Morinaga C, Saito D, Nakamura S, Sasaki T, Asakawa S, Shimizu N, et al. The hotei mutation of medaka in the anti-Müllerian hormone receptor causes the dysregulation of germ cell and sexual development. Proc Natl Acad Sci USA. 2007;104:9691–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Nakamura S, Watakabe I, Nishimura T, Picard JY, Toyoda A, Taniguchi Y, et al. Hyperproliferation of mitotically active germ cells due to defective anti-Müllerian hormone signaling mediates sex reversal in medaka. Development. 2012;139:2283–7.PubMedView ArticleGoogle Scholar
- Valenzuela N, Neuwald JL, Literman R. Transcriptional evolution underlying vertebrate sexual development. Dev Dyn. 2013;242:307–19.PubMedView ArticleGoogle Scholar
- Callard GV, Tchoudakova AV, Kishida M, Wood E. Differential tissue distribution, developmental programming, estrogen regulation and promoter characteristics of cyp19 genes in teleost fish. J Steroid Biochem Mol Biol. 2001;79:305–14.PubMedView ArticleGoogle Scholar
- Johnsen H, Tveiten H, Torgersen JS, Andersen Ø. Divergent and sex-dimorphic expression of the paralogs of the Sox9-Amh- cyp19a1 regulatory cascade in developing and adult atlantic cod (Gadus morhua L.). Mol Reprod Dev. 2013;80:358–70.PubMedView ArticleGoogle Scholar
- Wang DF, Kobayashi T, Zhou LY, Paul-Prasanth B, Ijiri S, Sakai F, et al. Foxl2 up-regulates aromatase gene transcription in a female-specific manner by binding to the promoter as well as interacting with Ad4 binding protein/steroidogenic factor 1. Mol Endocrinol. 2007;21:712–25.PubMedView ArticleGoogle Scholar
- Pannetier M, Fabre S, Batista F, Kocer A, Renault L, Jolivet G, et al. FOXL2 activates P450 aromatase gene transcription: towards a better characterization of the early steps of mammalian ovarian development. J Mol Endrocrinol. 2006;36:399–413.View ArticleGoogle Scholar
- Liu Z, Wu F, Jiao B, Zhang X, Hu C, Huang B, et al. Molecular cloning of doublesex and mab-3-related transcription factor 1, forkhead transcription factor gene 2, and two types of cytochrome P450 aromatase in Southern catfish and their possible roles in sex differentiation. J endocrinol. 2007;194:223–41.PubMedView ArticleGoogle Scholar
- Herpin A, Adolfi MC, Nicol B, Hinzmann M, Schmidt C, Klughammer J, et al. Divergent expression regulation of gonad development genes in medaka shows incomplete conservation of the downstream regulatory network of vertebrate sex determination. Mol Biol Evol. 2013;30(10):2328–46.PubMed CentralPubMedView ArticleGoogle Scholar
- Maatouk DM, DiNapoli L, Taketo MM, Capel B. Investigating the role of beta-catenin in sex determination. Biol Reprod. 2008;78:577.View ArticleGoogle Scholar
- Sreenivasan R, Jiang J, Wang X, Bártfai R, Kwan HY, Chrostoffels A, et al. Gonad differentiation in zebrafish is regulated by the canonical wnt signaling pathway. Biol Reprod. 2014;90:45.PubMedView ArticleGoogle Scholar
- Nicol B, Guerin A, Fostier A, Guiguen Y. Ovary-predominant wnt4 expression during gonadal differentiation is not conserved in the rainbow trout (Oncorhynchus mykiss). Mol Reprod Dev. 2013;79:51–63.View ArticleGoogle Scholar
- Vainio S, Heikkila M, Kispert A, Chin N, McMahon AP. Female development in mammals is regulated by Wnt-4 signalling. Nature. 1999;397:405–9.PubMedView ArticleGoogle Scholar
- Kim Y, Kobayashi A, Sekido R, DiNapoli L, Brennan J, Chaboissier MC, et al. Fgf9 and Wnt4 act as antagonistic signals to regulate mammalian sex determination. PLoS Biology. 2006;4(6):e187.PubMed CentralPubMedView ArticleGoogle Scholar
- Nicol B, Yano A, Jouanno E, Guerin A, Fostier A, Guiguen Y. Follistatin is an early player in rainbow trout ovarian differentiation and is both localized with aromatase and regulated by the Wnt pathway. Sex Dev. 2013;7:267–76.PubMedGoogle Scholar
- Oshima Y, Hayashi T, Tokonuga S, Nakamura M. Wnt4 expression in the differentiating gonad of the frog Rana rugosa. Zoolog Sci. 2005;22:689–93.PubMedView ArticleGoogle Scholar
- Smith CA, Sinclair AH. Sex determination: insights from the chicken. Bioessays. 2004;26(2):120–32.PubMedView ArticleGoogle Scholar
- Shoemaker CM, Crews D. Analyzing the coordinated gene network underlying temperature-dependent sex determination in reptiles. Semin Cell Dev Biol. 2009;20:293–303.PubMed CentralPubMedView ArticleGoogle Scholar
- Ramkissoon Y, Goodfellow P. Early steps in mammalian sex determination. Curr Opin Genet Dev. 1996;6:316–21.PubMedView ArticleGoogle Scholar
- Knower K, Kelly S, Harley V. Turning on the male – SRY, SOX9 and sex determination in mammals. Cytogenet Genome Res. 2003;101:185–98.PubMedView ArticleGoogle Scholar
- Lambeth LS, Raymond CS, Roeszler KN, Kuroiwa A, Nakata T, Zarkower D, et al. Over-expression of DMRT1 induces the male pathway in embryonic chicken gonads. Developmental Biology. 2014;389:160–72.PubMed CentralPubMedView ArticleGoogle Scholar
- Yang CX, Wright EC, Ross JW. Expression of RNA-binding proteins DND1 and FXR1 in the porcine ovary, and during oocyte maturation and early embryo development. Mol Reprod Dev. 2012;79:541–52.PubMedView ArticleGoogle Scholar
- Mortensen RD, Serra M, Steitz JA, Vasudevan S. Posttranscriptional activation of gene expression in Xenopus laevis oocytes by microRNA-protein complexes (microRNPs). Proc Natl Acad Sci U S A. 2011;108:8281–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Viñas A, Taboada X, Vale L, Robledo D, Hermida M, Vera M, et al. Mapping of DNA sex-specific markers and genes related to sex differentiation in turbot (Scophthalmus maximus). Mar Biotechnol. 2012;14(5):655–63.PubMedView ArticleGoogle Scholar
- Ospina-Álvarez N, Piferrer. Temperature-dependent sex determination in fish revisited: prevalence, a single sex ratio response pattern, and possible effects of climate change. PLoS One. 2008;3(7):e2837.PubMed CentralPubMedView ArticleGoogle Scholar
- Radder RS, Quinn AE, Georges A, Sarre SD, Shine R. Genetic evidence for co-occurrence of chromosomal and thermal sex-determining systems in a lizard. Biol Lett. 2008;4:176–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Holleley CE, O’Meally D, Sarre SD, Graves JAM, Ezaz T, Matsubara K, et al. Sex reversal triggers the rapid transition from genetic to temperature-dependent sex. Nature. 2015;523:79–82.PubMedView ArticleGoogle Scholar
- Andreu-Agullo C, Marin T, Thompson CB, Lai EC. Ars2 maintains neural stem-cell identity through direct transcriptional activation of Sox2. Nature. 2012;481:195–8.PubMed CentralView ArticleGoogle Scholar
- Capel B. R-spondin1 tips the balance in sex determination. Nat Genet. 2006;38(11):1233–4.PubMedView ArticleGoogle Scholar
- Zhang H, Ren C, Gao X, Takahashi T, Sapolsky RM, Steinberg GK, et al. Hypothermia blocks β-catenin degradation after focal ischemia in rats. Brain Res. 2008;10:182–7.View ArticleGoogle Scholar
- Tsai CL, Wang LH, Shiue YL, Chao TY. Influence of temperature on the ontogenic expression of neural development-related genes from developing tilapia brain expressed sequence tags. Mar Biotechnol. 2007;9:243–61.PubMedView ArticleGoogle Scholar
- Santerre C, Sourdaine P, Marc N, Mingant C, Robert R, Martinez A. Oyster sex determination is influenced by temperature – first clues in spat during first gonadic differentiation and gametogenesis. Comp Biochem Physiol. 2013;165:61–9.View ArticleGoogle Scholar
- Sheng Y, Chen B, Zhang L, Luo M, Cheng H, Zhou R. Identification of Dmrt genes and their up-regulation during gonad transformation in the swamp eel (Monopterus albus). Mol Biol Rep. 2014;41:1237–45.PubMedView ArticleGoogle Scholar
- Matsushita Y, Oshima Y, Nakamura M. Expression of DMRT genes in the gonads of Rana rugosa during sex determination. Zoolog Sci. 2007;24:95–9.PubMedView ArticleGoogle Scholar
- Navarro-Martín L, Galay-Burgos M, Sweeney G, Piferrer F. Different sox17 transcripts during sex differentiation in sea bass, Dicentrarchus labrax. Mol Cell Endocrinol. 2009;299:240–51.PubMedView ArticleGoogle Scholar
- Wang R, Cheng H, Xia L, Guo Y, Huang X, Zhou R. Molecular cloning and expression of Sox17 in gonads during sex reversal in the rice field eel, a teleost fish with a characteristic of natural sex transformation. Biochem Biophys Res Commun. 2003;303:452–7.PubMedView ArticleGoogle Scholar
- Taboada X, Robledo D, Del Palacio L, Rodeiro A, Felip A, Martínez P, et al. Comparative expression analysis in mature gonads, liver and brain of turbot (Scophthalmus maximus) by cDNA-AFLPS. Gene. 2012;492(1):250–61.PubMedView ArticleGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3 - new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.PubMed CentralPubMedView ArticleGoogle Scholar
- Ribas L, Pardo BG, Fernández C, Alvarez-Dios JA, Gómez-Tato A, Quiroga MI, et al. A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in turbot (Scophthalmus maximus). BMC Genomics. 2013;14:180.PubMed CentralPubMedView ArticleGoogle Scholar
- Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009;37:e45.PubMed CentralPubMedView ArticleGoogle Scholar
- Robledo D, Hernández-Urcera J, Cal RM, Pardo BG, Sánchez L, Martínez P, et al. Analysis of qPCR reference gene stability determination methods and a practical approach for efficiency calculation on a turbot (Scophthalmus maximus) gonad dataset. BMC Genomics. 2014;15:648.PubMed CentralPubMedView ArticleGoogle Scholar
- Kubista M, Sindelka R, Tichopad A, Bergkvist A, Lindh D, Forootan A. The prime technique. Real-time PCR data analysis. GIT Laboratory Journal. 2007;9-10:33–5.Google Scholar
- Development Core Team R. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2008.Google Scholar
- Venables WN, Ripley BD. Modern applied statistics with S. 4th ed. New York: Springer; 2002.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.PubMedView ArticleGoogle Scholar
- Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.PubMed CentralPubMedView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.PubMed CentralPubMedView ArticleGoogle Scholar