Skip Navigation


Integrative and Comparative Biology Advance Access originally published online on May 13, 2008
Integrative and Comparative Biology 2008 48(4):476-485; doi:10.1093/icb/icn031
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
48/4/476    most recent
icn031v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (1)
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Valenzuela, N.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology. All rights reserved. For permissions please email: journals.permissions@oxfordjournals.org.

Evolution of the gene network underlying gonadogenesis in turtles with temperature-dependent and genotypic sex determination

Nicole Valenzuela1
Department of Ecology, Evolution and Organismal Biology, Iowa State University, 253 Bessey Hall, Ames, IA 50011, USA

Correspondence: 1Email: nvalenzu{at}iastate.edu


    Synopsis
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
The evolution of sex determination has long fascinated biologists, as it has paramount consequences for the evolution of a multitude of traits, from sex allocation to speciation and extinction. Explaining the diversity of sex-determining systems found in vertebrates (genotypic or GSD and temperature-dependent or TSD) requires a comprehensive and integrative examination from both a functional and an evolutionary perspective. Particularly revealing is the examination of the gene network that regulates gonadogenesis. Here, I review some advances in this field and propose some additional hypotheses about the composition of the gene network underlying sexual development, the functional links among some of its elements and their evolution in turtles. I focus on several pending questions about: (1) What renders TSD systems thermo-sensitive? (2) Is there one developmentally conserved or multiple TSD mechanisms? (3) Have evolutionarily derived GSD species lost all ancestral thermal-sensitivity? New data are presented on embryonic expression of Dax1 (the dosage-sensitive sex-reversal adrenal hypoplasia congenital on the X chromosome gene in the turtles Chrysemys picta (TSD) and Apalone mutica (GSD). No differential Dax1 expression was detected in C. picta at any of the stages examined, consistent with reports on two other TSD turtles and alligators. Notably, significantly higher Dax1 expression was found at 30°C than at 25°C at stage 15 in A. mutica (GSD), likely caused by Wt1's identical expression pattern previously reported. Because Sf1 is an immediate downstream target of Dax1 and its expression is not affected by temperature, it is proposed that Sf1 renders Dax1's differential signal ineffective to induce biased sex ratios in A. mutica, as previously proposed for Wt1's thermosensitive expression. Thus, it is hypothesized that Sf1 plays a major role in the lack of response of sex ratio to temperature of A. mutica, and may function as a sex-determining gene in this GSD species. These and previous data permit formulating several mechanistic hypotheses: (1) the postulation of Wt1 as a candidate thermal master switch alone, or in combination with Sf1, in the TSD turtle C. picta; (2) the proposition of Sf1 as a sex-determining gene in the GSD turtle A. mutica; and (3) the hypothesis that differing patterns of gene expression among TSD taxa reflect multiple traits from a developmental perspective. Moreover, the recent finding of relic differential Wt1 expression in A. mutica and the results for Dax1 in this species provide empirical evidence that GSD taxa can harbor thermal sensitivity at the level of gene expression, potentially co-optable during TSD evolution.


    Introduction
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
The evolutionary survival of most dioecious organisms depends on the production of males and females sufficient to enable sexual reproduction. Surprisingly, the strong selection assumed to be acting to ensure that adults of both sexes are produced in adequate numbers has not resulted in a conserved mechanism of sex determination, but instead, it has generated an amazing array of sex-determining systems among taxa (Bull 1983Go; Solari 1994Go; Valenzuela et al. 2003Go; Valenzuela and Lance 2004Go). Explaining this diversity requires a comprehensive and integrative examination from both a functional and an evolutionary perspective. Particularly revealing in this quest is the examination of the gene network that regulates gonadogenesis. Indeed, we must understand how alternative sex-determining mechanisms operate from an ecological and developmental point of view if we are to decipher why such recurrent transitions among mechanisms have taken place.

Sex determination is the process of irreversible commitment to the male or female developmental fate, which can be initially triggered by genetic factors such as those contained in sex chromosomes (GSD), or by environmental factors (ESD) such as temperature (TSD) (Bull 1983Go; Valenzuela and Lance 2004Go). Studies of model organisms have generated a wealth of information about the control of the developmental steps of cell determination and cell differentiation that transform a unicellular zygote into an adult male or female (Gilbert 2000Go; Carroll et al. 2005Go). This has revealed the existence of a delicate balance between short-range and long-range signaling molecules that are required for the formation of the bipotential gonadal pimordia, the determination to a male developmental fate or a female one, and the ensuing differentiation into testes or ovaries (Kim and Capel 2006Go). Comparative analyses have demonstrated that the gene regulatory network underlying these events contains numerous elements common to all vertebrates (reviewed in Place and Lance 2004Go) but is triggered by different initial activators in GSD and TSD systems.

Importantly, this gene network has been extensively studied in model organisms (Fig. 1), and only more recently in non-model vertebrates, such as reptiles. Consequently, our current working paradigm for reptiles is built mostly upon fragmentary data from individual genes and species under the mostly untested assumption that the functional links of this network are largely the same as those known from model vertebrates. However, given that significant evolutionary divergence appears to have accrued in the elements of the gene network of sexual development (Wilkins 2002Go; Haag and Doty 2005Go; Valenzuela and Shikano 2007Go; Valenzuela 2008aGo), relying on a chimerical model may not reflect the true nature of this network and instead, it may obscure our functional and evolutionary understanding of it. For instance, if we are to decipher the nature of the evolutionary divergence of the gene network underlying sexual development to test current hypotheses about the mode in which such evolution occurred (e.g., Wilkins 1995Go, Graham et al. 2003Go, Nuzhdin et al. 2004Go), and the forces responsible for such changes (e.g., True and Haag 2001Go; Haag and Doty 2005Go), it is essential that we compare species-specific networks in non-model species. This is particularly crucial when studying evolutionary changes among closely related taxa that exhibit a wide variation of sex-determining mechanisms, as is the case of reptiles.


Figure 1
View larger version (55K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1 Gene network and cellular interactions underlying gonadogenesis in GSD mammalian model systems, based on Koopman 2001Go with modifications according to Lovell-Badge et al. 2002Go; Wilkins 2002Go; Yao et al. 2002Go; Lasala et al. 2004Go; Place and Lance 2004Go; Salmon et al. 2005Go; Crews et al. 2006Go; Li et al. 2007Go; Yamaguchi et al. 2007Go; Wilhelm et al. 2007Go. Elements with known reptilian homologs are indicated in the dark gray boxes. The list of elements of the gene network is not exhaustive. Key to symbols is included at the bottom of the graph.

 
Reptilian homologs of mammalian and avian genes involved in gonadal differentiation have been identified (Fig. 1 gray boxes; e.g., Aromatase, Amh, Dax1, Dmrt1, Foxl2, Sox9, Sf1, Wnt4, Wt1), and indicate that common elements exist between GSD and TSD regulatory networks. Some of these genes exhibit temperature-specific expression in TSD species, concordant to their expression by sex in GSD taxa. Here I focus on four genes in particular: Wt1, Sf1, Aromatase, and Dax1.

Wt1, the Wilms’ tumor-suppressor gene, encodes a transcription factor important for urogenital development (Kent et al. 1995Go; Roberts 2005Go) and it is necessary for the establishment of the bipotential gonad and later for the maintenance of Sertoli cells and seminiferous tubules in the testes of mammals (Gao et al. 2006Go, Fig. 1) and other vertebrates (e.g., Kent et al. 1995Go; Semba et al. 1996Go; Yamamura et al. 2005Go). Sf1, the steroidogenic factor 1 (or AD4BP and NR5A1) is a nuclear orphan receptor important for sexual differentiation in vertebrates, as it is required for the formation of primary steroidogenic organs (adrenal gland and gonad) in mammals, and for the expression of steroidogenic enzymes therein (Parker and Schimmer 1997Go; Morohashi 1999Go). Aromatase encodes the enzyme that converts androgens to estrogens and is important in many estrogen-dependent functions, including sexual differentiation (reviewed by Place and Lance 2004Go; Fig. 1). Dax1 (dosage-sensitive sex-reversal [DSS], adrenal hypoplasia congenital [AHC] critical region on the X chromosome gene 1) encodes an orphan receptor that is a member of the nuclear receptor superfamily (McCabe 2007Go), and is involved in ovarian differentiation in mammals (Ramkissoon and Goodfellow 1996Go) and other vertebrates, including reptiles (Pieau et al. 1999Go; Western et al. 2000Go).

Important insights can be gained not only by comparing multiple elements of the gene network underlying sexual development, but also by contrasting patterns of gene expression across GSD and TSD reptiles from a network perspective. For instance, previous studies described significant differential expression of Sf1 and Wt1 at stage 12 in Chrysemys picta turtles (TSD) prior to the onset of the thermosensitive period (TSP) (Valenzuela et al. 2006Go; Valenzuela 2008aGo). Additionally, Wt1 expression was also affected by temperature at this early developmental stage in Apalone mutica turtles (GSD) (Valenzuela 2008aGo). Given that Dax1 is a negative regulator of the transactivating functions of Sf1, particularly those mediated by Wt1 (Iyer and McCabe 2004Go), and given that Wt1 is also a known up regulator of Dax1 (Fig. 1; Kim et al. 1999Go), it is important to examine the expression of Dax1 in the context of what we know about other genes involved in gonadogenesis in these two taxa.

Dax1 has a C-terminal portion whose structure is characteristic of a ligand-binding domain, and an N-terminal portion formed by an atypical DNA-binding domain with 3.5 amino acid repeats (McCabe 2007Go). Importantly, Dax1 has an inhibitory effect on Sf1-mediated transactivation of other genes. For instance, Dax1 obstructs the combined synergism between Wt1 and Sf1 that is required during sexual development (Fig. 1) by the formation of heterodimers. Two Dax1 spliceoforms exist at equal proportion in testis (Ho et al. 2004Go, Hossain et al. 2004Go). The most recently described alternative spliceoform is encoded by exon 1 and exon 2A, the latter of which is contained within the intron, between exons 1 and 2 (Ho et al. 2004Go; Hossain et al. 2004Go). Dax1 exhibits sexually monomorphic expression in Trachemys scripta, Lepidochelys olivacea, and Alligator mississipiensis (Western et al. 2000Go; Maldonado et al. 2002Go; Shoemaker et al. 2007bGo).

Here, I analyze Dax1 expression during the embryogenesis of C. picta and A. mutica, and formalize some hypotheses about the composition of the gene network underlying sexual development, the functional links among some of its elements and their evolution in turtles, based on these novel and previous data.


    Materials and methods
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
Collection of samples
Freshly laid eggs were collected from natural nests of C. picta (7 clutches) and A. mutica (17 clutches), and correspond to the same eggs used previously (Valenzuela et al. 2006Go). Eggs were incubated at constant temperatures that produce ~100% male (25°C) and ~100% female (30°C) C. picta hatchlings (Ewert and Nelson 1991Go). Temperature does not bias sex ratios in A. mutica (GSD) (Janzen 1993Go).

Eggs from each clutch were uniformly distributed between incubation temperatures. Eggs from each clutch were divided equally among multiple plastic containers filled 3/5 of the way with moistened vermiculite set at –150 kPa by adding 338 g of water per 300 g of vermiculite. Boxes containing only the substrate were weighed (±0.1 g), and they were weighed again after adding the eggs and the weight recorded. Eggs were partially buried in the substrate in random positions in the box. Prior to the sampling of any egg (or weekly otherwise), the boxes were reweighed and lost moisture (as determined by weight loss) was replaced. The new weight was recorded after the removal of any eggs. Egg containers were rotated at least weekly to control for the effect of any potential thermal cline in the incubators.

Embryos were collected at several developmental stages prior to and during the TSP (sensu Yntema 1968Go) equivalent for both temperature treatments and species used: Stage 9 = Before TSP, Stage 12 = Before TSP, Stage 15 = TSP Onset, Stage 19 = Middle TSP, Stage 22 = End TSP (Bull and Vogt 1981Go). Collection of embryos from each clutch was also distributed across all developmental stages sampled. Sampled embryos were placed immediately in 10 vol. of RNAlater®, stored at –20°C and subsequently at –80°C for later use, following the manufacturer's recommendations.

Cloning and QPCR
RNA was extracted from the adrenal-kidney-gonadal complex (whole embryos or trunks were used from stages 9 and 12 embryos, respectively) using QIAGEN's RNeasy Kits and DNAse-I digested to prevent DNA contamination. RNA was quantified using a NanoDrop® ND-1000 Spectrophotometer, and its quality assessed by the presence of ribosomal bands in agarose gels stained with ethidium-bromide. Individual samples were kept separate and analyzed without pooling. Total RNA (1 µg/sample) was retro-transcribed with (dT)20 primers using Superscript III (Invitrogen) following the manufacturer's protocol. For samples that yielded <1 µg total RNA in 8 µl elute volume, as much as 8 µl were used for the RT-PCR, and the total amount of RNA was recorded for standardization during the analysis of data as described below.

Degenerate primers were designed for relatively conserved regions of vertebrate Dax1 cDNA sequences as found in GeneBank (Table 1). These primers amplified a 268 bp Dax1 cDNA fragment in both species, which was cloned in a pGEM®-T Easy Vector System (Promega) and sequenced. A 343 bp β-actin cDNA fragment (Valenzuela et al. 2006Go) was used to test for the quality of the extracted RNA during QPCR, and for normalization of Dax1 gene expression. Internal primers amplified a 124 bp β-actin fragment (Valenzuela et al. 2006Go), and 98 and 79 bp Dax1 fragments in C. picta and A. mutica, respectively (Table 1), during real-time QPCR using Brilliant® SYBR® Green QPCR Master Mix in an Mx3000P thermocycler (Stratagene). ROX was used as the reference die for background correction. Standard curves were generated from pure miniprep cloned DNA obtained above by diluting it at concentrations of 5, 1, 1 x 10–2, 1 x 10–4, 1 x 10–6, 1 x 10–12, 1 x 10–14, and 1 x 10–16 ng/µl, and run in duplicate in each QPCR to ensure technical repeatability of the results. Samples from all clutches, temperatures and stages were included in each 96-well plate used for the QPCR experiments to avoid any systematic bias. QPCR conditions followed Valenzuela et al. (2006Go): 1 cycle at 95°C for 10 min; 45 cycles of 95°C for 30 s, 60°C for 1 min, 72°C for 1 min; and a dissociation-curve cycle of 95°C for 1 min, 55°C for 30 s, and 95°C for 30 s. Dissociation profiles obtained at the end of the PCR cycles were inspected to confirm amplification of a single transcript.


View this table:
[in this window]
[in a new window]

 
Table 1 Primers used for standard PCR or quantitative-real time PCR (those with qm prefix) amplification of Wt1 and β-actin cDNA fragments in Chrysemys picta and Apalone mutica

 

    Analysis of data
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
Initial template amount for each sample and each gene was calculated using standard-curve quantification via the algorithms implemented in Mx3000p v2.0 from Stratagene with background correction. Initial copy numbers were calculated from the molecular weight of the fragments for each gene, and standardized to 1 µg of initial total RNA. Expression of Wt1 was normalized using β-actin expression data, and Log2 transformed to correct for the heteroscedasticity and nonnormal distribution characteristic of amplification data. Statistical analysis was carried out on these normalized Log2-transformed data keeping individual values separate (Valenzuela et al. 2006Go). The significance of the temperature-treatment effect was determined by testing for differences in the mean values of gene expression at each sampled developmental stage via an ANOVA as implemented in Jmp 5.1.2© (SAS Institute 2004Go).


    Results
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
The fragments of Dax1 cloned from embryonic adrenal-kidney-gonadal tissue encompassed part of the 3' end of the exon 1 and part of the 5' region of the exon 2, which correspond to the C-terminal ligand-binding domain and share high similarity with other members of the nuclear hormone-receptor-gene superfamily (Guo et al. 1996Go; McCabe 2007Go). Given the location of the primers, the transcript abundance measured by QPCR does not include the recently described alternative spliceoform (Ho et al. 2004Go; Hossain et al. 2004Go). The sequence of the Dax1 fragment of C. picta's (GenBank EU433884) was 93.6% identical to A. mutica's (GenBank EU433885), 97.8% to T. scripta's, 95.0% to L. olivacea, 91.8% to alligators's, 86.9% to chicken's, 77.0% to human's, 72.4% to those of mice. Differences were also found at the amino acid level (Table 2). The amino acid sequence of the Dax1 fragment of C. picta's was 93.3% identical to A. mutica's, 95.6% to T. scripta's, 94.1% to L. olivacea, 92.2% to alligators's, 91.1% to chicken's, 74.4% to human's, 71.1% to those of mice. The amino acid sequences of C. picta and A. mutica differ by three residues from that of T. scripta (Shoemaker et al. 2007) but the changes are unique to each of these two species (Table 2).


View this table:
[in this window]
[in a new window]

 
Table 2 Amino-acid sequence of Dax1 cDNA fragments in selected species of vertebrates

 
Standard curves for β-actin and Dax1 had R2 values of 0.97 and 0.998, respectively, for C. picta, and of 0.98 and 0.991, respectively, for A. mutica. No significant differences in Dax1 expression were detected in C. picta at any of the stages examined. Notably, in A. mutica, Dax1 expression was significantly higher at 25°C than at 30°C during stage 15 (F = 5.2593, P = 0.0323; Fig. 2).


Figure 2
View larger version (23K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2 Developmental expression of Dax1 in C. picta (TSD) and A. mutica turtles (GSD). Y-axis represents Dax1 initial copy number normalized to β-actin (+SD). Boxed stages correspond to C. picta's TSP. Sample sizes at 25°C (italics) and at 30°C (bold) are presented per sampling time. Stars denote significant differential expression. Results are robust irrespective of outliers. Differences in transcript abundance at all other stages are not statistically significant and thus, no biological significance can be attributed to the crossing of the expression profiles.

 
Significance of the results is neither due to any systematic bias generated by a small number of clutches represented, nor to clutch-mates being over-represented in any particular developmental stage or temperature. For A. mutica, 17 clutches were represented, and a single egg per clutch was sampled at all temperature-by-stage combinations. For C. picta, 7 clutches were represented, and only in stages 19 and 22 did two of the total eggs sampled at each temperature come from the same clutch. Readings from these clutch-mate pairs were not more similar than from any pair of unrelated individuals sampled at the same temperature-by-stage combination. Additionally, although duplicate readings from embryonic samples were precluded because the limited amount of total RNA per sample was used to profile seven different genes (Valenzuela et al. 2006Go; Valenzuela and Shikano 2007Go; Valenzuela 2008aGo and unpublished data), duplicate readings from the standard curves were 97–99.9% identical, thus ensuring the technical reliability of the readings of transcript abundance.


    Discussion
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
An important step in explaining the outstanding lability of sex-determining mechanisms requires a better understanding of the regulation and evolution of the gene network that controls sexual development. Taxonomic groups such as reptiles that encompass a wide variety of mechanisms are particularly informative as they allow a comparative analysis of evolutionary changes that have accrued among closely related taxa with alternative sex-determining systems. Such is the case of the target species used in this study, C. picta (TSD) and A. mutica (GSD), which have been the subject of recent comparative analyses of developmental gene expression (Valenzuela et al. 2006Go; Valenzuela and Shikano 2007Go; Valenzuela 2008aGo).

In this study, expression of Dax1 was insensitive to temperature at all stages examined. Consistently, Dax1 expression during the TSP of developing gonads of T. scripta, L. olivacea, and of the gonad-adrenal-mesonephros of A. mississipiensis was not significantly different between male-producing and female-producing temperatures (Western et al. 2000Go; Maldonado et al. 2002Go; Shoemaker et al. 2007bGo). However, significant differential Dax1 expression was detected at stage 15 of development in A. mutica. The potential functional and evolutionary significance of this observation emerges when this pattern of expression is compared to that of other genes of the sexual development network in a comparative framework.

For instance, results from previous studies on C. picta and A. mutica permit the examination of several important questions. One is the still pending question of what renders TSD systems thermo-sensitive. Such a thermal master switch will be a gene(s) that activates the TSP or that enables responses during this developmental window of time. Candidates for this role would be genes whose expression occurs prior to, or exactly at, the onset of the TSP, rather than genes exhibiting differential expression within the TSP after it has already been activated. Because early Sf1 expression is required for the formation of bipotential gonads (Wilkins 2002Go, Fig. 1), Sf1 was recently proposed as one such potential master switch in the TSD turtle C. picta, based on its early significant differential expression at stage 12, prior to the onset of the TSP (Valenzuela et al. 2006Go). However, Wt1 was also found to have significant differential expression at stage 12 in C. picta (Valenzuela 2008aGo). Early Wt1 expression is also required for the formation of the bipotential gonads (Gao et al. 2006Go, Fig. 1). Because Wt1 is known to regulate Sf1 in other vertebrates (Wilhelm and Englert 2002Go), Wt1 was proposed as another candidate master switch in TSD turtles (Valenzuela 2008aGo), which, given the regulatory function of Wt1 on Sf1, makes it an even more likely candidate than Sf1 (Fig. 3). The hypothetical role of Wt1 (alone, or in combination with Sf1) as a thermal master switch is consistent with the evolutionary developmental model which proposes that the hierarchy of the genetic pathway underlying sex determination evolves from the bottom-up as a series of evolutionary steps that favor the addition of neomorphic regulators upstream of the cascade (Wilkins 1995Go; Graham et al. 2003Go).


Figure 3
View larger version (32K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3 Proposed hypothetical models of developmental gene networks regulating gonadogenesis under TSD and GSD. Models are based on patterns of gene expression from C. picta and A. mutica as described in the text. Warmer and colder refers to incubation temperature during development. Arrows indicate positive regulation or activation, and lines with end caps indicate negative regulation or repression. An X on a line indicates a signal that has been rendered ineffective. Dots within the male and female boxes denote all remaining gene and cellular interactions leading to sexual development as summarized in Fig. 1.

 
A second question is whether TSD encompasses a single or multiple mechanisms. This relates to the fundamental question of whether TSD mechanisms evolved in one developmentally conserved way, or whether multiple molecular pathways evolved to produce ecologically equivalent outcomes (e.g., different TSD molecular pathways may exist that result in the same sex-by-temperature reaction-norm, such as the pattern of male-to-female production with increasing temperature or TSDIa). The diversity of systems harbored within GSD (Wilkins 2002Go) provides evidence that such evolutionary changes have occurred in the evolution of sex-determining mechanisms, rendering the possibility for similar changes within TSD more likely. Alternatively, concordant expression of genes across species would be expected if TSD systems are pleisiomorphic and developmentally conserved. A comparative analysis of data from C. picta and A. mutica with that of other turtles reported in the literature, indicated that substantial evolutionary changes have accrued not only between TSD and GSD systems, but among TSD species as well (Valenzuela et al. 2006Go; Valenzuela and Shikano 2007Go; Valenzuela 2008aGo) in genes such as Sf1, Aromatase, and Wt1. These differences at the level of gene expression between TSD species suggest that there are multiple TSD systems that differ in their molecular underpinnings of development (sensu Abouheif and Wray 2002Go), perhaps as the result of divergence over evolutionary time (True and Haag 2001Go; Nuzhdin et al. 2004Go). Importantly, this comparative analysis also revealed that the fundamental difference between TSD and GSD mechanisms is not driven by Aromatase as previously proposed (reviewed by Pieau and Dorizzi 2004Go) because patterns of either differential or insensitive aromatase gene expression or activity during the TSP are found across TSD and GSD vertebrates (Valenzuela and Shikano 2007Go).

A third important issue is whether gene expression is totally thermo-insensitive in GSD taxa, or not. For turtles, whose ancestral state appears to be TSD (Janzen and Krenz 2004Go), this relates to the evolutionary question of whether GSD turtles have lost all thermal sensitivity in the regulation of the gene network underlying sexual development. Data from A. mutica provided evidence that this is not always the case, because this GSD turtle has retained its ancestral sensitivity to the expression of a gene involved in gonadogenesis, Wt1, the first such case ever to be reported (Valenzuela 2008aGo). The finding of differential expression in Wt1 in A. mutica is important because such a thermally sensitive gene may have been the key factor responsible for the thermal sensitivity of the TSD ancestor. Consistent with this idea, Wt1 was independently proposed as a potential thermal master switch based on gene-expression data from C. picta (Valenzuela 2008aGo). Of further relevance to the Dax1 results obtained here, it should be pointed out that such relict thermal-sensitivity must be overridden now by the master sex-determining gene in the GSD A. mutica.

At stages 12 and 15 of development, A. mutica exhibited significant and marginally significant differential Wt1 expression, respectively (Valenzuela 2008aGo). Remarkably, the expression pattern of Dax1 at those two stages mimics the pattern of Wt1 and the transcript abundance values of Wt1 and Dax1 are highly correlated with one another (r = 0.999, P = 0.0013). Because Wt1 is a known upregulator of Dax1 (Fig. 1; Kim et al. 1999Go), this result supports the functional link between Wt1 and Dax1 in A. mutica. The question follows as to why it is that the differential expression of Dax1 at stage 15 does not produce biased sex ratios in this species? The data point once more to Sf1 (Valenzuela 2008aGo). Dax1 inhibits the function of Sf1 in the transactivation of several other genes (Iyer and McCabe 2004Go). Because Sf1 is a downstream target of Dax1, and because Sf1 is insensitive to temperature in A. mutica (Valenzuela et al. 2006Go), Dax1's thermal signal is rendered ineffective by Sf1 in the same manner that Wt1's signal appears to be (Valenzuela 2008aGo, Fig. 3). Thus, the two lines of evidence (this study and that of Valenzuela 2008aGo) indicate Sf1 has an important role in controlling the lack of response to temperature of the sex ratio in A. mutica. I propose that Sf1 may function as a sex-determining gene in this GSD species. Further genomic and functional research is necessary to test this hypothesis.

In summary, it is argued here that a more comprehensive understanding of the functional composition and evolution of the gene network underlying sexual development requires a comparative approach at the level of the genes and species. Such efforts using C. picta and A. mutica have enabled the formulation of a series of mechanistic hypotheses. These include (1) the postulation of Wt1 as a candidate thermal master switch alone or in combination with Sf1 in the TSD turtle C. picta (Fig. 3); (2) the proposition of Sf1 as a sex-determining gene in the GSD turtle A. mutica (Fig. 3); and (3) the hypothesis that differing patterns of gene expression among TSD taxa reflect multiple traits from a developmental perspective. Moreover, the finding of relic differential Wt1 expression in A. mutica and the results for Dax1 in this same species provide empirical evidence that GSD taxa can harbor thermal sensitivity at the level of expression of genes involved in gonadal formation. This thermosensitive gene expression in GSD taxa represents a co-optable source for the evolution of phenotypic plasticity such as TSD, as previously proposed (Bull 1981Go, 1983Go; reviewed in Valenzuela 2004Go).

More research is needed to elucidate the forces responsible for the evolutionary divergence of gene expression patterns. Encouragingly, research on reptiles using a multigene approach are increasing (e.g., Maldonado et al. 2002Go; Rhen et al. 2007Go; Shoemaker et al. 2007aGo, bGo) and that should help in generating species-specific models of this gene network. It should be noted that more attention should be paid to other aspects of the regulatory gene network underlying sexual development. For instance, heterochronic changes (Alberch et al. 1979Go) in the timing of specification and determination between TSD and GSD mechanisms will alter the function of this network in ways that may alleviate the strong selection on the master trigger of male and female sexual development, thus facilitating the evolutionary transitions between GSD and TSD (Valenzuela 2008bGo). An integrated approach that includes those considerations, along with research on the ecological context in which alternative gene expression patterns are found in nature, should help reveal how males and females are produced, and the forces responsible for generating the existing variety of sex-determining mechanisms.


    Acknowledgments
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
I thank T. Shikano and A. LeClere for their assistance with some collection of data. This study was partially funded by the Center for Integrative Genomics at Iowa State University, and by NSF grant IOS 0743284 to Nicole Valenzuela.


    Footnotes
 
From the symposium "Reptile Genomics and Evolutionary Genetics" presented at the annual meeting of the Society for Integrative and Comparative Biology, January 2–6, 2008, at San Antonio, Texas.


    References
 Top
 Synopsis
 Introduction
 Materials and methods
 Analysis of data
 Results
 Discussion
 Acknowledgments
 References
 
Abouheif E, Wray GA. Evolution of the gene network underlying wing polyphenism in ants. Science (2002) 297::249–52.[Abstract/Free Full Text]

Alberch P, Gould SJ, Oster GF, Wake DB. Size and shape in ontogeny and phylogeny. Paleobiology (1979) 5::296–317.[Abstract]

Bull JJ. Evolution of environmental sex determination from genotypic sex determination. Heredity (1981) 47::173–84.[CrossRef][Web of Science]

Bull JJ. Evolution of sex determining mechanisms (1983) Menlo Park (CA): Benjamin/Cummings.

Bull JJ, Vogt RC. Temperature-sensitive periods of sex determination in emydid turtles. J Exp Zool (1981) 218::435–40.[CrossRef][Web of Science][Medline]

Bullejos M, Bowles J, Koopman P. Searching for missing pieces of the sex determination puzzle. J Exp Zool (2001) 290::517–22.[CrossRef][Web of Science][Medline]

Carroll SB, Grenier JK, Weatherbee SD. From DNA to diversity: molecular genetics and the evolution of animal design (2005) Malden (MA): Blackwell Publishing.

Cho S, Huang ZY, Zhang JZ. Sex-specific splicing of the honeybee doublesex gene reveals 300 million years of evolution at the bottom of the insect sex-determination pathway. Genetics (2007) 177::1733–41.[Abstract/Free Full Text]

Crews D, Lou W, Fleming A, Ogawa S. From gene networks underlying sex determination and gonadal differentiation to the development of neural networks regulating sociosexual behavior. Brain Res (2006) 1126::109–21.[CrossRef][Web of Science][Medline]

Ewert MA, Nelson CE. Sex determination in turtles: diverse patterns and some possible adaptive values. Copeia (1991) 1991::50–69.[CrossRef]

Gao F, Maiti S, Alam N, Zhang Z, Deng JM, Behringer RR, Lecureuil C, Guillou F, Huff V. The Wilms tumor gene, Wt1, is required for Sox9 expression and maintenance of tubular architecture in the developing testis. Proc Natl Acad Sci USA (2006) 103::11987–92.[Abstract/Free Full Text]

Gilbert SF. Developmental biology (2000) Sunderland (MA): Sinauer Associates, Inc.

Graham P, Penn JKM, Schedl P. Masters change, slaves remain. Bioessays (2003) 25::1–4.[CrossRef][Web of Science][Medline]

Guo WW, Burris TP, Zhang YH, Huang BL, Mason J, Copeland KC, Kupfer SR, Pagon RA, McCabe ERB. Genomic sequence of the DAX1 gene: an orphan nuclear receptor responsible for X-linked adrenal hypoplasia congenita and hypogonadotropic hypogonadism. J Clin Endocrinol Metab (1996) 81::2481–6.[Abstract]

Haag ES, Doty AV. Sex determination across evolution: connecting the dots. Plos Biol (2005) 3::21–4.[CrossRef][Web of Science]

Ho J, Zhang YH, Huang BL, McCabe ERB. NR0B1A: an alternatively spliced form of NR0B1. Mol Gen Metab (2004) 83::330–6.[CrossRef][Web of Science][Medline]

Hossain A, Li C, Saunders GF. Generation of two distinct functional Isoforms of dosage-sensitive sex reversal-adrenal hypoplasia congenita-critical region on the X chromosome gene 1 (DAX-1) by alternative splicing. Mol Endocrinol (2004) 18::1428–37.[Abstract/Free Full Text]

Iyer AK, McCabe ERB. Molecular mechanisms of DAX1 action. Mol Genet Metab (2004) 83::60–73.[CrossRef][Web of Science][Medline]

Janzen FJ. The influence of incubation temperature and family on eggs, embryos, and hatchlings of the smooth softshell turtle (Apalone mutica). Physiol Zool (1993) 66::349–73.

Janzen FJ, Krenz JG. Phylogenetics: which was first, TSD or GSD? In: Temperature dependent sex determination in vertebrates—Valenzuela N, Lance VA, eds. (2004) Washington (DC): Smithsonian Books. 121–30.

Jeyasuria P, Roosenburg WM, Place AR. Role of P-450 Aromatase in sex determination of the diamondback terrapin, Malaclemys terrapin. J Exp Zool (1994) 270::95–111.[CrossRef][Web of Science][Medline]

Kent J, Coriat AM, Sharpe PT, Hastie ND, Vanheyningen V. The evolution of Wt1 sequence and expression pattern in the vertebrates. Oncogene (1995) 11::1781–92.[Web of Science][Medline]

Khattak JZK, Torp AM, Andersen SB. A genetic linkage map of Spinacia oleracea and localization of a sex determination locus. Euphytica (2006) 148::311–8.[CrossRef][Web of Science]

Kim J, Prawitt D, Bardeesy N, Torban E, Vicaner C, Goodyer P, Zabel B, Pelletier J. The Wilms’ tumor suppressor gene (wt1) product regulates Dax-1 gene expression during gonadal differentiation. Mol Cell Biol (1999) 19::2289–99.[Abstract/Free Full Text]

Kim Y, Capel B. Balancing the bipotential gonad between alternative organ fates: a new perspective on an old problem. Dev Dyn (2006) 235::2292–300.[CrossRef][Web of Science][Medline]

Koopman P. Gonad development: signals for sex. Curr Biol (2001) 11::R481–R483.[CrossRef][Web of Science][Medline]

Lasala C, Carre-Eusebe D, Picard JY, Rey R. Subcellular and molecular mechanisms regulating anti-Mullerian hormone gene expression in mammalian and nonmammalian species. DNA Cell Biol (2004) 23::572–85.[Web of Science][Medline]

Li SR, Lu MM, Zhou DY, Hammes SR, Morrisey EE. GLP-1: a novel zinc finger protein required in somatic cells of the gonad for germ cell development. Dev Biol (2007) 301::106–16.[CrossRef][Web of Science][Medline]

Lovell-Badge R, Canning C, Sekido R. Sex-determining genes in mice: building pathways. Genetics and Biology of Sex Determination (2002) 244::4–22.[CrossRef][Web of Science]

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. Gen Comp Endocrinol (2002) 129::20–26.[CrossRef][Web of Science][Medline]

McCabe ERB. DAX1: increasing complexity in the roles of this novel nuclear receptor. Mol Cell Endocrinol (2007) 265::179–82.[CrossRef][Web of Science][Medline]

Morohashi K. Gonadal and extragonadal functions of Ad4BP/SF-1: developmental aspects. Trends Endocrinol Metab (1999) 10::169–73.[CrossRef][Web of Science][Medline]

Nuzhdin SV, Wayne ML, Harmon KL, McIntyre LM. Common pattern of evolution of gene expression level and protein sequence in Drosophila. Mol Biol Evol (2004) 21::1308–17.[Abstract/Free Full Text]

Parker KL, Schimmer BP. Steroidogenic factor 1: a key determinant of endocrine development and function. Endocrine Rev (1997) 18::361–77.[Abstract/Free Full Text]

Pieau C, Dorizzi M. Oestrogens and temperature-dependent sex determination in reptiles: all is in the gonads. J Endocrinol (2004) 181::367–77.[Abstract]

Pieau C, Dorizzi M, Richard-Mercier N. Temperature-dependent sex determination and gonadal differentiation in reptiles. Cell Mol Life Sci (1999) 55::887–900.[CrossRef][Web of Science][Medline]

Pires-DaSilva A. Evolution of the control of sexual identity in nematodes. Semin Cell Dev Biol (2007) 18::362–70.[CrossRef][Web of Science][Medline]

Place AR, Lance VA. The temperature-dependent sex determination drama – same cast, different stars. In: Temperature dependent sex determination in vertebrates—Valenzuela N, Lance VA, eds. (2004) Washington (DC): Smithsonian Books. 99–110.

Ramkissoon Y, Goodfellow P. Early steps in mammalian sex determination. Curr Opin Genet Dev (1996) 6::316–21.[CrossRef][Web of Science][Medline]

Rhen T, Metzger K, Schroeder A, Woodward R. Expression of putative sex-determining genes during the thermosensitive period of gonad development in the snapping turtle, Chelydra serpentina. Sex Dev (2007) 1::255–70.[CrossRef][Medline]

Roberts SGE. Transcriptional regulation by WT1 in development. Curr Opin Genet Dev (2005) 15::542–7.[CrossRef][Web of Science][Medline]

Salmon NA, Handyside AH, Joyce IM. Expression of Sox8, Sf1, Gata4, Wt1, Dax1, and Fog2 in the mouse ovarian follicle: Implications for the regulation of Amh expression. Mol Rep Dev (2005) 70::271–7.[CrossRef][Web of Science][Medline]

SAS Institute. IJMP (2004) Cary (NC): SAS Institute, Inc.

Semba K, SaitoUeno R, Takayama G, Kondo M. cDNA cloning and its pronephros-specific expression of the Wilms’ tumor suppressor gene, WT1, from Xenopus laevis. Gene (1996) 175::167–72.[CrossRef][Web of Science][Medline]

Shoemaker C, Ramsey M, Queen J, Crews D. Expression of Sox9, Mis, and Dmrt1 in the gonad of a species with temperature-dependent sex determination. Dev Dyn (2007a) 236::1055–63.[CrossRef][Web of Science][Medline]

Shoemaker CM, Queen J, Berkstresser K, Crews D. Response to candidate sex-determining genes to changes in temperature in the red-eared slider turtle. Dev Biol (2007b) 306::411–2.

Solari AJ. Sex chromosomes and sex determination in vertebrates (1994) Boca Raton (FL): CRC Press, Inc.

True JR, Haag ES. Developmental system drift and flexibility in evolutionary trajectories. Evol Dev (2001) 3::109–19.[CrossRef][Web of Science][Medline]

Valenzuela N. Evolution and maintenance of temperature-dependent sex determination. In: Temperature dependent sex determination in vertebrates—Valenzuela N, Lance VA, eds. (2004) Washington (DC): Smithsonian Books. 131–47.

Valenzuela N. Relic thermosensitive gene expression in genotypically-sex-determined turtles. Evolution (2008a) 62::234–40.[Web of Science][Medline]

Valenzuela N. Sexual development and the evolution of sex determination. Sexual Dev (2008b) (in press).

Valenzuela N, Adams DC, Janzen FJ. Pattern does not equal process: exactly when is sex environmentally determined? Am Nat (2003) 161::676–83.[CrossRef][Web of Science][Medline]

Valenzuela N, Lance VA. Temperature dependent sex determination in vertebrates (2004) Washington (DC): Smithsonian Books.

Valenzuela N, LeClere A, Shikano T. Comparative gene expression of steroidogenic factor 1 in Chrysemys picta and Apalone mutica turtles with temperature-dependent and genotypic sex determination. Evol Dev (2006) 8::424–32.[CrossRef][Web of Science][Medline]

Valenzuela N, Shikano T. Embryological ontogeny of Aromatase gene expression in Chrysemys picta and Apalone mutica turtles: comparative patterns within and across temperature-dependent and genotypic sex-determining mechanisms. Dev Genes Evol (2007) 217::55–62.[CrossRef][Web of Science][Medline]

Western PS, Harry JL, Marshall JA, Graves AH Sinclair. Temperature-dependent sex determination in the American alligator: expression of SF1, WT1 and DAX1 during gonadogenesis. Gene (2000) 241::223–32.[CrossRef][Web of Science][Medline]

Wilhelm D, Englert C. The Wilms tumor suppressor WT1 regulates early gonad development by activation of Sf1. Gen Dev (2002) 16::1839–51.[Abstract/Free Full Text]

Wilhelm D, Palmer S, Koopman P. Sex determination and gonadal development in mammals. Physiol Rev (2007) 87::1–28.[Abstract/Free Full Text]

Wilkins AS. Moving up the hierarchy - a hypothesis on the evolution of a genetic sex determination pathway. Bioessays (1995) 17::71–7.[CrossRef][Web of Science][Medline]

Wilkins AS. The evolution of developmental pathways (2002) Sunderland (MA): Sinauer Associates, Inc.

Yamaguchi T, Yamaguchi S, Hirai T, Kitano T. Follicle-stimulating hormone signaling and Foxl2 are involved in transcriptional regulation of aromatase gene during gonadal sex differentiation in Japanese flounder, Paralichthys olivaceus. Biochemi Biophys Res Comm (2007) 359::935–40.[CrossRef]

Yamamura Y, Aoyama S, Oshima Y, Kato T, Osawa N, Nakamura M. Molecular cloning and expression in gonad of Rana rugosa WT1 and Fgf9. Zool Sci (2005) 22::1045–50.[CrossRef][Web of Science][Medline]

Yao HHC, Whoriskey W, Capel B. Desert Hedgehog/Patched 1 signaling specifies fetal Leydig cell fate in testis organogenesis. Gen Dev (2002) 16::1433–40.[Abstract/Free Full Text]

Yntema CL. A series of stages in the embryonic development of Chelydra serpentina. J Morph (1968) 125::219–52.[CrossRef][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
CSH ProtocolsHome page
N. Valenzuela
The Painted Turtle, Chrysemys picta: A Model System for Vertebrate Evolution, Ecology, and Human Health
CSH Protocols, July 1, 2009; 2009(7): pdb.emo124 - pdb.emo124.
[Abstract] [Full Text]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
48/4/476    most recent
icn031v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (1)
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Valenzuela, N.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?