Molecular cloning of doublesex genes of four cladocera (water flea) species
© Toyota et al.; licensee BioMed Central Ltd. 2013
Received: 29 September 2012
Accepted: 26 March 2013
Published: 10 April 2013
The gene doublesex (dsx) is known as a key factor regulating genetic sex determination in many organisms. We previously identified two dsx genes (DapmaDsx1 and DapmaDsx2) from a freshwater branchiopod crustacean, Daphnia magna, which are expressed in males but not in females. D. magna produces males by parthenogenesis in response to environmental cues (environmental sex determination) and we showed that DapmaDsx1 expression during embryonic stages is responsible for the male trait development. The D. magna dsx genes are thought to have arisen by a cladoceran-specific duplication; therefore, to investigate evolutionary conservation of sex specific expression of dsx genes and to further assess their functions in the environmental sex determination, we searched for dsx homologs in four closely related cladoceran species.
We identified homologs of both dsx genes from, D. pulex, D. galeata, and Ceriodaphnia dubia, yet only a single dsx gene was found from Moina macrocopa. The deduced amino acid sequences of all 9 dsx homologs contained the DM and oligomerization domains, which are characteristic for all arthropod DSX family members. Molecular phylogenetic analysis suggested that the dsx gene duplication likely occurred prior to the divergence of these cladoceran species, because that of the giant tiger prawn Penaeus monodon is rooted ancestrally to both DSX1 and DSX2 of cladocerans. Therefore, this result also suggested that M. macrocopa lost dsx2 gene secondarily. Furthermore, all dsx genes identified in this study showed male-biased expression levels, yet only half of the putative 5’ upstream regulatory elements are preserved in D. magna and D. pulex.
The all dsx genes of five cladoceran species examined had similar amino acid structure containing highly conserved DM and oligomerization domains, and exhibited sexually dimorphic expression patterns, suggesting that these genes may have similar functions for environmental sex determination in cladocerans.
KeywordsCladocera Doublesex genes Environmental sex determination Gene duplication Annotated regulatory regions
Sex determination is a fundamental developmental process, affecting the sexual differentiation of gonads, and leads to sex-specific differences in behavior, physiology and morphology. Sex-determining systems can be divided into two categories: genotypic sex determination (GSD) and environmental sex determination (ESD) [1–3]. GSD is attributed to the genetic segregation of genes, often residing on sex chromosomes that initiate alternate sex-determining developmental pathways. In contrast, ESD has repeatedly arisen during animal evolution  and is initiated by diverse environmental cues, such as temperature, photoperiod, nutrition and population density, that trigger alternative genetic signals, resulting in the regulation of male or female sex-determining genes [5, 6].
Natural selection of rare mutational variants has been suggested to mediate the transitions between GSD and ESD [7, 8]. A previous phylogenetic analysis revealed that there have been at least three independent switches from GSD to ESD in lizards, and six transitions from ESD to GSD in turtles . Moreover, previous experiments using temperature-sensitive mutations created artificially in Caenorhabditis elegans demonstrated how GSD could rapidly evolve into ESD as a consequence of a mutation in key sex determining genes . Orthologs of GSD genes such as dmrt1, sox9 and cyp19a (aromatase) are expressed in the gonads during the temperature-sensitive period in ESD of reptiles . Thus, according to the current interpretation of these data, ESD mechanisms are likely to share many sex-determining components with GSD .
Sex determination systems in insects vary considerably in key factors and regulatory mechanisms to develop sex-specific traits. The sex determination mechanism in Drosophila melanogaster is best understood. The ratio of X chromosomes to autosomes (X:A ratio) is thought to provide the initial signal for the activation of sex-lethal (sxl), a master gene of the sex determination cascade. Then, sxl is produced as the sex-specific splicing isoforms. Sxl in female acts on the pre-mRNA of transformer (tra) resulting in reproduction of functional Tra. The functional Tra in the female, in concert with Tra-2, regulates the production of female-specific doublesex (dsx) mRNA. The male-specific splice form of dsx mRNA is the default splice-variant in D. melanogaster. Dsx regulates the various sex-specific traits such as gonads. Recently, sex determination mechanisms have also been demonstrated in various insect lineages such as Diptera (Musca domestica and Ceratitis capitata), Hymenoptera (Apis mellifera and Nasonia vitripennnis) and Coleoptera (Tribolium castaneum). These studies revealed that tra and dsx are highly conserved among insects [11–14]. However, in case of Lepidoptera, Bombyx mori, tra and tra-2 are assumed not to be required for the sex-specific splicing of Bmdsx pre-mRNA, because Bmdsx has no Tra/Tra-2 binding motif. Recently, it has been revealed that binding of the BmPSI, a Bombyx homolog of P-element somatic inhibitor, to the exonic splicing suppressor sequence on expected region is involved in sex-specific splicing of Bmdsx. These data suggest that upstream genetic cascades of dsx might be diverse among insects.
The Cladocera (commonly called water fleas) is an ancient clade of branchiopod crustaceans comprising 16 or 18 family lineages [15, 16] that all reproduce by cyclical parthenogenesis involving ESD . The most well studied species are from the family Daphniidae, particularly of the genus Daphnia. Daphnia inhabit freshwater ponds and lakes on all continents and are known to switch between parthenogenetic and sexual reproduction when environmental conditions for growth and reproduction deteriorate. During normal growing conditions, populations are most often entirely composed of females. However, shortened photoperiod, a lack of food and/or increased population density all lead to the clonal production of males that are genetically identical to their sisters and mothers . First instar male juveniles are easily distinguished from the females . During maturation, daphniids undergo morphological sexual differentiation of various somatic tissues, including the armament of a first thoracic leg with the copulatory hook in males, which becomes larger during the fifth instar . Gonads develop and finally settle at both sides of the gut during embryogenesis in both sexes . The appearance of males allows sexual reproduction to occur [22, 23], when females begin producing haploid eggs requiring fertilization.
Recently, we and others found that male production occurred independently of environmental cues by treatment with exogenous juvenile hormone (JH) or its analogs [24, 25]. Exposure of D. magna to JH analogs at the stage corresponding to the environmentally-sensitive period for sex determination of a cladoceran species of the family Moinidae , produced exclusively male broods, suggesting that JH could be a key molecule for understanding mechanisms of ESD [24, 27, 28].
A doublesex (dsx) gene was originally identified in D. melanogaster as a critical and terminal transcription factor in the fly sex determining cascade. Dsx is spatially and temporally transcribed into two sex-specific splice forms conferring sexually dimorphic traits during development [29, 30]. The dsx gene contains two conserved domains: the Dsx/Mab-3 (DM) domain at the N-terminus and the oligomerization domain at the C-terminus . Genes encoding DM-domain (DM-domain genes) were discovered to play a related role in C. elegans[32, 33] and also in vertebrates [34–36]. In contrast, results from numerous studies have shown that other genes in the genetic sex-determination cascade widely diversified among species [1, 2, 37].
To understand the molecular and evolutionary relationships between GSD and ESD, we previously identified and analyzed three DM-domain genes (DMRT11E, DMRT93B and DMRT99B) from D. magna, displaying sexual dimorphic gene expression patterns in adult gonads . However, none of these DM-domain genes exhibited sexually dimorphic expression patterns during embryonic development, suggesting that they are not involved in sex determination . Two additional DM-domain genes were later found in the D. magna expressed sequence tags (ESTs) database . Therefore, we analyzed the function of these two genes from D. magna using gene manipulations that we developed . These experiments revealed that two dsx genes in D. magna were obtained by lineage-specific duplication, and then one of the paralogs, Daphnia magna dsx1 (DapmaDsx1), plays an important role in directing the major sexually dimorphic development of D. magna. In contrast, specific function of Daphnia magna dsx2 (DapmaDsx2) remains unknown. These newly identified dsx genes showed greater sequence similarity at the amino acid sequence level to known insect dsx genes than to the previously identified DM-domain containing genes in D. magna. A genome-wide study of gene functions in D. pulex suggested that lineage-specific duplicated genes are most responsive to varying environmental conditions . In the present study, we investigated the sequence and functional conservation of the two dsx genes in a broader taxonomic sampling of cladocerans by cloning dsx homologs, and determining their sex specific expression in four species representing two families and three genera. We also analyzed the structures of cloned dsx genes of D. magna and D. pulex including their putative regulatory motifs and putative transcription factor binding sites in the 5’ upstream regions of these duplicated dsx genes.
Results and discussion
Molecular cloning of doublesexgenes from cladocerans
By comparing the DapmaDsx1 sequence to orthologs from the four studied species (D. pulex, D. galeata, C. dubia and M. macrocopa), we found that DapmaDsx1 shows 88-58%, 100-95%, 88-48%, 100-78% and 96-61% sequence identities to the A, B (DM-domain), C, D (oligomerization-domain) and E domains, respectively (Figure 2A). Similarly, by comparing the DapmaDsx2 sequence with identified orthologs, we observed that DapmaDsx2 shows comparable relative ratios among each of the domains: 74-66%, 98-97%, 88-67%, 100-78% and 92-69% sequence identities to the A, B, C, D and E domains, respectively (Figure 3A). These results suggest that putative amino acids of both the DM- and oligomerization-domains are highly conserved among the Cladocera. On the other hand, amino acid similarities outside of these domains are lower, and are proportional to the evolutionary distance between each genus; Daphnia, Ceriodaphnia and Moina (Additional files 2 with 3).
The DM-domain contains zinc chelating residues, and among the insects studied to date, two highly conserved amino acid residues, threonine and glutamine (Boxed in Figures 2B, 3B), distinguish the DM-domain of DSX from DM-domains of other insect proteins . Therefore, we searched for similar highly conserved amino acid residues within the DSX DM-domains of cladocerans. Indeed, all zinc chelating residues are found to be conserved in the DM domains of DSX1 and DSX2 among the five cladoceran species (Figures 2B, 3B). Yet, although the threonine and glutamine residues were conserved in DSX2, the glutamine residue in DSX1 was substituted by arginine in all cladoceran species examined (Figures 2B, 3B). These results suggest that DSX1 in cladocerans might have gained a novel function affecting sex determination by amino acid replacement after duplication of dsx in branchiopoda lineage.
Compared to the DM-domains, more amino acid variation is observed in alignments of the oligomerization-domains (Figures 2C, 3C). Dimerization, which enhances specific DNA binding, is mediated by several residues in a non-polar interface that is conserved within oligomerization-domains. Previous study has revealed that, in DSX2 of D. magna, two of three nonpolar amino acids indispensable for formation of the nonpolar interface are substituted with the polar acidic amino acid, aspartic acid . This suggests that the daphniid DSX2 proteins are unable to dimerize and may not be functional or may have a different, unknown functions.
Sex specific expression of dsxgenes in five cladoceran species
We previously reported that both dsx genes in D. magna were transcriptionally up-regulated in males and showed no sex-specific splice isoforms. We and others also reported that exposure to JH analogs reliably produces male daphniids [24, 25]. Moreover, by using gene knock-down (RNAi) and overexpression methods in D. magna[40, 41], we discovered that DapmaDsx1 is necessary and sufficient for sex determination in D. magna, whereas the tandemly duplicated DapmaDsx2 gene does not determine sex, even though its transcription is equally sex-biased .
Annotation of dsx gene structures in D. Magna and D. Pulex
We previously cloned and described the mRNA transcripts of the DapmaDsx1 and DapmaDsx2 genes . DapmaDsx1 produces two mRNA, DapmaDsx1-α and DapmaDsx1-β, which are expressed in both sexes and differ only in their 5’ UTR, while DapmaDsx2 produces only one mRNA transcript. The D. pulex draft genome sequence was recently published , and the D. magna genome sequencing project is currently in progress by the Daphnia Genomics Consortium. The DapmaDsx1-α, DapmaDsx1-β, and DapmaDsx2 mRNA transcripts were aligned to the genome assemblies and used to annotate the dsx1 and dsx2 gene models in the D. magna and D. pulex genomes (Additional file 4).
The D. magna dsx gene cluster is located on scaffold 2190 of the D. magna genome assembly v2.4 with a ~10 Kbp intergenic region between the DapmaDsx1 (~16.1 Kbp length) and DapmaDsx2 (~1.6 Kbp length) genes (Additional file 4A). The second exon of the DapmaDsx1-β mRNA transcript fell within an assembly gap in scaffold 2190, but was located on scaffold 521. We conclude that scaffold 521 (~5.8 Kbp length) characterizes a ~3.1 Kbp gap in scaffold 2190, located in the intragenic region of DapmaDsx1.
The D. pulex dsx gene cluster is located on scaffold 32 of the D. pulex genome assembly v1.1 with a ~9.6 Kbp intergenic region between the DappuDsx1 (~21.5 Kbp length) and DappuDsx2 (~1.6 Kbp length) genes (Additional file 4B). The D. pulex v1.1 gene model predictions did not correctly identify the dsx1 and dsx2 duplicate gene cluster. Our annotations improved the dsx gene models in both D. magna and D. pulex genomes.
Potential transcriptional regulatory elements in the 5’ upstream promoter regions of the dsxgenes
To further annotate putative functional and conserved elements of the dsx genes, we searched and compared transcriptional promoter regions of the genes in D. magna and D. pulex. 1.0 Kbp upstream of the transcription start site (TSS) of dsx1-α, dsx1-β, and dsx2 were extracted as transcriptional regulatory regions. This interval spans the predicted intergenic region between adjacent loci. Promoter sequences are challenging for multiple alignment algorithms, because upstream regulatory regions are not well conserved compared to protein coding regions of genes [47, 48]. We aligned orthologous promoter sequences from D. magna and D. pulex using Pro-Coffee , an alignment algorithm specifically designed for homologous promoter regions (Additional files 5, 6, 7). The DapmaDsx1-α and DappuDsx1-α promoter alignment showed 46% sequence identity (Additional file 5), the DapmaDsx1-β and DappuDsx1-β promoter alignment showed 60% sequence identity (Additional file 6), and the DapmaDsx2 and DappuDsx2 promoter alignment showed 62% sequence identity (Additional file 7). The promoter regions of dsx1 and dsx2 have much less conservation than their respective protein coding regions, which have 84% and 83% sequence identity, respectively (92% and 84% sequence identity at synonymous sites on protein coding region, respectively).
Since we previously reported that DapmaDsx1 and DapmaDsx2 are paralogs, and that both DapmaDsx1 and DapmaDsx2 mRNAs are transcriptionally up-regulated in male D. magna, we searched for de novo conserved motifs present in all dsx promoter regions in D. magna and D. pulex (DapmaDsx1-α, DappuDsx1-α, DapmaDsx1-β, DappuDsx1-β, DapmaDsx2, and DappuDsx2), without reference to TFBS sequence databases. We identified 14 conserved motifs in the Daphnia dsx promoters (Additional file 12), which can later be functionally investigated as potential TFBSs and/or potential transcriptional promoters of Daphnia dsx. The motifs were labeled M1 through M14 and annotated onto the D. magna/D. pulex dsx promoter alignments (Figures 6B, 6D and 6F). Motifs 1 through 9 were also found in the D. melanogaster dsx upstream promoter region, supporting the conservation and potential regulatory functions of these motifs.
In order to assess whether these de novo conserved motifs are similar to any known TFBSs, we scanned the individual motif sequences for matches to TFBS PFMs and compared each motif consensus sequence to known TFBS consensus sequences. Several of the conserved motifs showed similarity to known TFBSs (Figures 7B-F). Motif 3 matches the TFBS of Mirr, a homeobox transcription factor in the Iroquois complex of transcription factors, which is predicted in all three dsx promoter TF-map alignments (Figure 7B). Motif 4 matches the TFBS of Vvl, a homeobox transcription factor which is predicted in both dsx1-β and dsx2 TF-map alignments (Figure 7C). Motif 7 matches the TFBS of Gsc/Bcd/Oc, three homeobox transcription factors with nearly identical binding sites (Figure 7D), with four of the Motif 7 sequences matching the Gsc/Bcd/Oc consensus TAATC exactly. Gsc was also predicted in the dsx2 TF-map alignment, and Oc was predicted in the dsx1-α TF-map alignment. Motif 8 matches the TFBS of Pan, a high mobility group transcription factor, which is predicted in both dsx1-β and dsx2 TF-map alignments (Figure 7E). Motif 12 also matches the TFBS of Pan, but with the reverse complement of the binding site (Figure 7F). The similarity of these de novo conserved motifs to known TFBSs that were also predicted by the TF-map alignment further supports our results describing the regulatory elements of the Daphnia dsx genes.
In summary, we identified the orthologs of DapmaDsx1 and DapmaDsx2 genes from closely related species belonging to two cladoceran families and three genera: D. pulex, D. galeata, C. dubia (Daphniidae) and M. macrocopa (Moinidae), with highly conserved DM- and oligomerization-domains. All five species examined exhibited sexually dimorphic expression pattern of dsx genes, suggesting that these genes may have similar functions for sex determination in cladocerans. Daphniids are unique animals that exhibit ESD and are; therefore, attractive for understanding the evolution of ESD. We also identified potential regulatory motifs and transcription factor binding sites in the putative promoter regions of these genes in D. magna and D. pulex. This information will facilitate future study of molecular mechanisms underlying sex-determination in cladocerans.
Daphniastrain and culture conditions
Isoclonal strains of D. magna (NIES and Belgium strains), D. galeata, C. dubia and M. macrocopa were obtained from the National Institute for Environmental Studies (NIES; Tsukuba, Japan) [25, 52]. D. pulex was obtained from Hokkaido University, Sapporo, Japan , and maintained as described previously . Briefly, culture medium was prepared using charcoal-filtered tap water and cultures of 20 individuals per liter were incubated at 21 ± 1°C under a 14-h light/10-h dark photoperiod. A 0.01-ml suspension of 4.3 × 108 cells ml−1 Chlorella (Chlorella vulgaris) was added daily to each culture. The water hardness was between 72 and 83 mg L−1, the pH between 7.0 and 7.5, and the dissolved oxygen concentration between 80 and 99%. To obtain natural male embryos, adult D. magna (Belgium clone) was reared in crowded conditions, and D. pulex was incubated at 18°C under a 10-h light/14-h dark photoperiod, and a 0.01-ml suspension of 4.3 × 108 cells ml−1 Chlorella was added every two days. To obtain male embryos of D. magna (NIES clone), D. galeata, C. dubia and M. macrocopa, (in which natural males are rarely seen) adult individuals (about 2 weeks of age) were chemically induced to produce males by treating them with a synthetic JH analog, fenoxycarb (1 μg/L) (technical grade 96.6% pure, Wako Pure Chemical Industries, Ltd., Osaka, Japan) . We confirmed the offspring sexes by the length of the first antenna  observed and photographed using a Leica MZ APO dissecting microscope (Leica, Mannheim, Germany).
Cloning of dsxgenes
The nucleotide sequences of the D. magna dsx genes were used for designing primers that amplify dsx genes in four different species. The harvested animals were homogenized using the Micro Smash MS-100R (Tomy, Tokyo, Japan). Total RNA was extracted with ISOGEN reagent according to the manufacturer’s protocol (NIPPON GENE, Tokyo, Japan). Poly (A) + RNA was isolated from purified total RNA using Fast Track (Life Technologies, Carlsbad, CA USA) and converted to cDNA using Superscript III and random primers (Life Technologies) according to the manufacturer’s protocol. cDNAs corresponding to the EST sequences were obtained by PCR amplification, and full-length cDNAs were obtained by RACE (Cap Fishing; SeeGene, Seoul, South Korea) using the oligonucleotide sequences as shown in Additional file 13.
Phylogenetic analysis of the DM-domain genes
A phylogenetic tree of DM-domain genes including newly cloned D. pulex, D. galeata, C. dubia and M. macrocopa dsx genes were constructed using amino acid sequences of DM-domain genes used in the previous study  (Additional file 14). A multiple alignment was constructed using Clustal W [45, 54] with the following settings (pairwise alignment parameters: gap opening penalty 15, gap extension penalty 6.66, identity protein weight matrix; multiple alignment parameters: delay divergent cutoff 30%, gap separation distance 4). Phylogenetic reconstruction was performed using the maximum likelihood and the neighbor-joining methods implemented in MEGA version 5 .
Two to three weeks old male and female animals of the five cladoceran species were used in quantitative-PCR (q-PCR) assays of gene expression levels. mRNAs were quantified as described previously . Animals were washed briefly and soaked in RNAlater (Life Technologies) for 10 min. Total RNA was purified and cDNA was synthesized as described above except that a random oligonucleotide was used as the primer. PCR was performed in an ABI Prism 7000 (Life Technologies) using the SYBR-Green PCR core reagents kit (Life Technologies), in the presence of appropriate primers. PCR amplifications were performed using the following conditions: 2 min at 50°C and 10 min at 95°C, followed by a total of 40 two-temperature cycles (15 s at 95°C and 1 min at 60°C).
The primers were chosen to amplify short PCR products of 150 bp; the primer sequences are listed in Additional file 15. Ribosomal protein L32 gene was used for normalization purposes [21, 55]. Data acquisition and analysis were performed by ABI Prism 7000 SDS software ver. 1.1 (Life Technologies). The baseline and threshold for the Ct (cycle threshold) were set automatically. Each gene was tested in technical triplicate samples by the relative standard curve method. In the case of D. magna and D. pulex, each experiment was performed in biological triplicate and statistical analyses were applied.
The genomic locations of DapmaDsx1-α, DapmaDsx1-β, and DapmaDsx2 mRNA transcripts were identified using BLASTN sequence similarity searches against a reference blast database of the D. magna genome assembly v2.4 scaffolds, and against a reference blast database of the D. pulex genome assembly v1.1 scaffolds. The best BLAST matches were analyzed and used to map the gene exons onto the D. magna and D. pulex scaffolds. ESTs mapped onto the genome assembly with PASA , microarray tiling path expression data, and RNA-Seq data from wFleaBase  were used as supporting evidence for exon annotations. The dsx gene annotation figures were created with AnnotationSketch  (Figure 6).
Transcription factor map alignments
Matscan  was used to search for matches to 125 JASPAR core insect TFBS matrices and 44 TRANSFAC insect TFBS matrices in each dsx promoter region in D. magna and D. pulex. A threshold of 0.85 matrix similarities was used to find TFBS matrix matches in the promoter sequences. For each promoter sequence, the collection of TFBS matrix matches is referred to as its TF-map. The TF-maps for each promoter sequence can be found in the following Additional files: DapmaDsx1-α (Additional file 16), DappuDsx1-α (Additional file 17), DapmaDsx1-β (Additional file 18), DappuDsx1-β (Additional file 19), DapmaDsx2 (Additional file 20), DappuDsx2 (Additional file 21). Meta  was then used to find the best meta-alignment of the orthologous promoters TF-maps (with parameters: a = 0.5, l = 0.1, m = 0.1).
de novoconserved promoter motifs
The 5′ 1.0 Kbp upstream region was extracted from all D. pulex gene models and used to create a background frequency model for 8 bp length motifs in Daphnia promoter regions. WeederH  was used to search for conserved regulatory motifs of length 8 present in all six dsx promoters in D. magna and D. pulex. The WeederH algorithm measures conservation based on the sequence conservation as well as the motif occurrence’s positions relative to the TSS. WeederH produces a χ2 score assessing how conserved the motif is compared to the rest of the homologous sequences. We used a χ2 score threshold of 3, discarding motifs with a χ2 score less than 3.
The D. magna and D. pulex sequence data were provided by The Center for Genomics and Bioinformatics at Indiana University, which is supported in part by the METACyt Initiative of Indiana University, funded in part through a major grant from the Lilly Endowment, Inc., and distributed via wFleaBase in collaboration with Dr. Donald Gilbert and the Daphnia Genomics Consortium https://wiki.cgb.indiana.edu/display/DGC/Home. This work was partly supported by the Japan Society for the Promotion of Science (JSPS) Research Fellowship for Young Scientists. to KT and CH, the Sasakawa Scientific Research Grant from The Japan Science Society to KT, the Saito Ho-on Kai Scientific Research Grant from The Saito Gratitude Foundation to KT, grants from the Ministry of Education, Culture, Sports, Science and Technology (TI), the Ministry of the Environment of Japan (TI), a grant of Long-Range Research Initiative by Japan Chemical Industry Association (TI) and a grant from National Institute for Basic Biology (TI).
- Marin I, Baker BS: The evolutionary dynamics of sex determination. Science. 1998, 281: 1990-1994.View ArticlePubMed
- Zarkower D: Establishing sexual dimorphism: conservation amidst diversity?. Nat Rev Genet. 2001, 2: 175-185.View ArticlePubMed
- Bull JJ: Sex determining mechanisms: an evolutionary perspective. Experientia. 1985, 41: 1285-1296. 10.1007/BF01952071.View ArticlePubMed
- Organ CL, Janes DE: Evolution of sex chromosomes in Sauropsida. Integr Comp Biol. 2008, 48: 512-519. 10.1093/icb/icn041.PubMed CentralView ArticlePubMed
- Crews D, Bull JJ: Mode and tempo in environmental sex determination in vertebrates. Semin Cell Dev Biol. 2009, 20: 251-255. 10.1016/j.semcdb.2009.02.004.View ArticlePubMed
- Korpelainen H: Sex ratios and conditions required for environmental sex determination in animals. Biol Rev Camb Philos Soc. 1990, 65: 147-184. 10.1111/j.1469-185X.1990.tb01187.x.View ArticlePubMed
- Bull JJ: Evolution of environmental sex determination from genotypic sex determination. Heredity. 1981, 47: 173-184. 10.1038/hdy.1981.74.View Article
- Bulmer MG, Bull JJ: Models of polygenic sex determination and sex ratio control. Evolution. 1982, 36: 13-26. 10.2307/2407962.View Article
- Chandler CH, Phillips PC, Janzen FJ: The evolution of sex-determining mechanisms: lessons from temperature-sensitive mutations in sex determination genes in Caenorhabditis elegans. J Evol Biol. 2009, 22: 192-200. 10.1111/j.1420-9101.2008.01639.x.View ArticlePubMed
- Shoemaker CM, Crews D: Analyzing the coordinated gene network underlying temperature-dependent sex determination in reptiles. Semin Cell Dev Biol. 2009, 20: 293-303. 10.1016/j.semcdb.2008.10.010.PubMed CentralView ArticlePubMed
- Sanchez L: Sex-determining mechanisms in insects. Int J Dev Biol. 2008, 52: 837-856. 10.1387/ijdb.072396ls.View ArticlePubMed
- Shukla JN, Nagaraju J: Doublesex: a conserved downstream gene controlled by diverse upstream regulators. J Genet. 2010, 89: 341-356. 10.1007/s12041-010-0046-6.View ArticlePubMed
- Gempe T, Beye M: Function and evolution of sex determination mechanisms, genes and pathways in insects. BioEssays. 2011, 33: 52-60. 10.1002/bies.201000043.PubMed CentralView ArticlePubMed
- Shukla JN, Palli SR: Sex determination in beetles: Production of all male progeny by parental RNAi knockdown of transformer. Sci Rep. 2012, 2: srep00602-
- Stenderup JT, Olesen J, Glenner H: Molecular phylogeny of the Branchiopoda (Crustacea)–multiple approaches suggest a ‘diplostracan’ ancestry of the Notostraca. Mol Phylogenet Evol. 2006, 41: 182-194. 10.1016/j.ympev.2006.06.006.View ArticlePubMed
- Olesen J: A phylogenetic analysis of the Conchostraca and Cladocera (Crustacea, Branchiopoda, Diplostraca). Zool J Linn Soc. 1998, 122: 491-536. 10.1111/j.1096-3642.1998.tb02161.x.View Article
- Taylor DJ, Crease TJ, Brown WM: Phylogenetic evidence for a single long-lived clade of crustacean cyclic parthenogens and its implications for the evolution of sex. Proc R Soc Lond B Biol Sci. 1999, 266: 791-797. 10.1098/rspb.1999.0707.View Article
- Kleiven OT, Larsson P, Hobæk A: Sexual reproduction in Daphnia magna requires three stimuli. Oikos. 1992, 65: 197-206. 10.2307/3545010.View Article
- Olmstead AW, LeBlanc GA: Effects of endcrine-active chemicals on the development of sex characteristics of Daphnia magna. Environ Toxicol Chem. 2000, 19: 2107-2113.View Article
- Mitchell SE: Intersex and male development in Daphnia magna. Hydrobiologia. 2001, 442: 145-156. 10.1023/A:1017564105942.View Article
- Sagawa K, Yamagata H, Shiga Y: Exploring embryonic germ line development in the water flea, Daphnia magna, by zinc-finger-containing VASA as a marker. Gene Expr Patterns. 2005, 5: 669-678. 10.1016/j.modgep.2005.02.007.View ArticlePubMed
- Stross RG, Hill JC: Diapause induction in Daphnia requires two stimuli. Science. 1965, 150: 1462-1464. 10.1126/science.150.3702.1462.View ArticlePubMed
- Hebert PD: Genotypic characteristics of cyclic parthenogens and their obligately asexual derivatives. Experientia Suppl. 1987, 55: 175-195.View ArticlePubMed
- Olmstead AW, LeBlanc GA: Juvenoid hormone methyl farnesoate is a sex determinant in the crustacean Daphnia magna. J Exp Zool. 2002, 293: 736-739. 10.1002/jez.10162.View ArticlePubMed
- Oda S, Tatarazako N, Watanabe H, Morita M, Iguchib T: Production of male neonates in four cladoceran species exposed to a juvenile hormone analog, fenoxycarb. Chemosphere. 2005, 60: 74-78. 10.1016/j.chemosphere.2004.12.080.View ArticlePubMed
- Banta AM, Brown LA: Control of sex in Cladocera. III. Localization of the critical period for control of sex. Proc Natl Acad Sci USA. 1929, 15: 71-81. 10.1073/pnas.15.2.71.PubMed CentralView ArticlePubMed
- Tatarazako N, Oda S, Watanabe H, Morita M, Iguchi T: Juvenile hormone agonists affect the occurrence of male Daphnia. Chemosphere. 2003, 53: 827-833. 10.1016/S0045-6535(03)00761-6.View ArticlePubMed
- Kato Y, Kobayashi K, Oda S, Tatarazako N, Watanabe H, Iguchi T: Sequence divergence and expression of a transformer gene in the branchiopod crustacean, Daphnia magna. Genomics. 2010, 95: 160-165. 10.1016/j.ygeno.2009.12.005.View ArticlePubMed
- Robinett CC, Vaughan AG, Knapp J-M, Baker BS: Sex and the single cell. II. There is a time and place for sex. PLoS Biol. 2010, 8: e1000365-10.1371/journal.pbio.1000365.PubMed CentralView ArticlePubMed
- Burtis KC, Baker BS: Drosophila doublesex gene controls somatic sexual differentiation by producing alternatively spliced mRNAs encoding related sex-specific polypeptides. Cell. 1989, 56: 997-1010. 10.1016/0092-8674(89)90633-8.View ArticlePubMed
- Bayrer JR, Zhang W, Weiss MA: Dimerization of doublesex is mediated by a cryptic ubiquitin-associated domain fold: implications for sex-specific gene regulation. J Biol Chem. 2005, 280: 32989-32996. 10.1074/jbc.M507990200.View ArticlePubMed
- Shen MM, Hodgkin J: mab-3, a gene required for sex-specific yolk protein expression and a male-specific lineage in C. elegans. Cell. 1988, 54: 1019-1031. 10.1016/0092-8674(88)90117-1.View ArticlePubMed
- Raymond CS, Shamu CE, Shen MM, Seifert KJ, Hirsch B, Hodgkin J, Zarkower D: Evidence for evolutionary conservation of sex-determining genes. Nature. 1998, 391: 691-695. 10.1038/35618.View ArticlePubMed
- Raymond CS, Murphy MW, O’Sullivan MG, Bardwell VJ, Zarkower D: Dmrt1, a gene related to worm and fly sexual regulators, is required for mammalian testis differentiation. Genes Dev. 2000, 14: 2587-2595. 10.1101/gad.834100.PubMed CentralView ArticlePubMed
- Matson CK, Zarkower D: Sex and the singular DM domain: insights into sexual regulation, evolution and plasticity. Nat Rev Genet. 2012, 13: 163-174.PubMed CentralPubMed
- Kopp A: Dmrt genes in the development and evolution of sexual dimorphism. Trends Genet. 2012, 28: 175-184. 10.1016/j.tig.2012.02.002.PubMed CentralView ArticlePubMed
- Williams TM, Carroll SB: Genetic and molecular insights into the development and evolution of sexual dimorphism. Nat Rev Genet. 2009, 10: 797-804.View ArticlePubMed
- Kato Y, Kobayashi K, Oda S, Colbourne JK, Tatarazako N, Watanabe H, Iguchi T: Molecular cloning and sexually dimorphic expression of DM-domain genes in Daphnia magna. Genomics. 2008, 91: 94-101. 10.1016/j.ygeno.2007.09.002.View ArticlePubMed
- Watanabe H, Tatarazako N, Oda S, Nishide H, Uchiyama I, Morita M, Iguchi T: Analysis of expressed sequence tags of the water flea Daphnia magna. Genome. 2005, 48: 606-609. 10.1139/g05-038.View ArticlePubMed
- Kato Y, Shiga Y, Kobayashi K, Tokishita S, Yamagata H, Iguchi T, Watanabe H: Development of an RNA interference method in the cladoceran crustacean Daphnia magna. Dev Genes Evol. 2011, 220: 337-345. 10.1007/s00427-011-0353-9.View ArticlePubMed
- Kato Y, Kobayashi K, Watanabe H, Iguchi T: Environmental sex determination in the branchiopod crustacean Daphnia magna: deep conservation of a Doublesex gene in the sex-determining pathway. PLoS Genet. 2011, 7: e1001345-10.1371/journal.pgen.1001345.PubMed CentralView ArticlePubMed
- Colbourne JK, Pfrender ME, Gilbert D, Thomas WK, Tucker A, Oakley TH, Tokishita S, Aerts A, Arnold GJ, Basu MK, Bauer DJ, Cáceres CE, Carmel L, Casola C, Choi J-H, Detter JC, Dong Q, Dusheyko S, Eads BD, Fröhlich T, Geiler-Samerotte KA, Gerlach D, Hatcher P, Jogdeo S, Krijgsveld J, Kriventseva EV, Kültz D, Laforsch C, Lindquist E, Lopez J: The ecoresponsive genome of Daphnia pulex. Science. 2011, 331: 555-561. 10.1126/science.1197761.PubMed CentralView ArticlePubMed
- Oliveira DCSG, Werren JH, Verhuls EC, Giebel JD, Kamping A, Beukeboom LW, Zande L: Identification and characterization of the doublesex gene of Nasonia. Insect Mol Biol. 2009, 18: 315-324. 10.1111/j.1365-2583.2009.00874.x.PubMed CentralView ArticlePubMed
- Volff J-N, Zarkower D, Bardwell VJ, Schartl M: Evolutionary dynamics of the DM domain gene family in metazoans. J Mol Evol. 2003, 57: 241-249. 10.1007/s00239-003-2470-1.View Article
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralView ArticlePubMed
- Swain TD, Taylor DJ: Structural rRNA characters support monophyly of raptorial limbs and paraphyly of limb specialization in water fleas. Proc R Soc Lond B Biol Sci. 2003, 270: 887-896. 10.1098/rspb.2002.2297.View Article
- Erb I, Bussotti G, Blanco E, Eyras E, Notredame C, lez-Vallinas JRG: Use of ChIP-Seq data for the design of a multiple promoter-alignment method. Nucleic Acids Res. 2012, 40: e52-10.1093/nar/gkr1292.PubMed CentralView ArticlePubMed
- Blanco E, Messeguer X, Smith TF, Guigo R: Transcription factor map alignment of promoter regions. PLoS Comp Biol. 2006, 2: e49-10.1371/journal.pcbi.0020049.View Article
- Bryne JC, Valen E, Tang M-HE, Marstrand T, Winther O, Piedade I, Krogh A, Lenhard B, Sandelin A: JASPAR, the open access database of transcription factor-binding profiles: new content and tools in the 2008 update. Nucleic Acids Res. 2008, 36: 102-106. 10.1093/nar/gkn449.View Article
- Matys V, Fricke E, Geffers R, Goßling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, Kloos D-U, Land S, Lewicki-Potapov B, Michael H, Munch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC(R): transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31: 374-378. 10.1093/nar/gkg108.PubMed CentralView ArticlePubMed
- Crooks GE, Hon G, Chandonia J-M, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2002, 14: 1188-1190.View Article
- Oda S, Kato Y, Watanabe H, Tatarazako N, Iguchi T: Morphological changes in Daphnia galeata induced by a crustacean terpenoid hormone and its analog. Environ Toxicol Chem. 2011, 30: 232-238. 10.1002/etc.378.View ArticlePubMed
- Hiruta C, Nishida C, Tochinai S: Abortive meiosis in the oogenesis of parthenogenetic Daphnia pulex. Chromosome Res. 2010, 18: 833-840. 10.1007/s10577-010-9159-2.View ArticlePubMed
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMed
- Kato Y, Kobayashi K, Oda S, Tatarazako N, Watanabe H, Iguchi T: Cloning and characterization of the ecdysone receptor and ultraspiracle protein from the water flea Daphnia magna. J Endocrinol. 2007, 193: 183-194. 10.1677/JOE-06-0228.View ArticlePubMed
- Haas BJ, Delcher AL, Mount SM, Wortman JR RKS, Hannick LI, Maiti R, Ronning CM, Rusch DB, Town CD, Salzberg SL, White O: Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003, 31: 5654-5666. 10.1093/nar/gkg770.PubMed CentralView ArticlePubMed
- Colbourne JK, Singan VR, Gilbert DG: wFleaBase: the Daphnia genome database. BMC Bioinformatics. 2005, 6: 45-10.1186/1471-2105-6-45.PubMed CentralView ArticlePubMed
- Steinbiss S, Gremme G, Scharfer C, Mader M, Kurtz S: AnnotationSketch: a genome annotation drawing library. Bioinformatics. 2009, 25: 533-534. 10.1093/bioinformatics/btn657.View ArticlePubMed
- Pavesi G, Zambelli F, Pesole G: WeederH: an algorithm for finding conserved regulatory motifs and regions in homologous sequences. BMC Bioinformatics. 2007, 8: 46-10.1186/1471-2105-8-46.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.