- Research article
- Open Access
An unbiased approach to identify genes involved in development in a turtle with temperature-dependent sex determination
© Chojnowski and Braun; licensee BioMed Central Ltd. 2012
- Received: 12 September 2011
- Accepted: 15 July 2012
- Published: 15 July 2012
Many reptiles exhibit temperature-dependent sex determination (TSD). The initial cue in TSD is incubation temperature, unlike genotypic sex determination (GSD) where it is determined by the presence of specific alleles (or genetic loci). We used patterns of gene expression to identify candidates for genes with a role in TSD and other developmental processes without making a priori assumptions about the identity of these genes (ortholog-based approach). We identified genes with sexually dimorphic mRNA accumulation during the temperature sensitive period of development in the Red-eared slider turtle (Trachemys scripta), a turtle with TSD. Genes with differential mRNA accumulation in response to estrogen (estradiol-17β; E2) exposure and developmental stages were also identified.
Sequencing 767 clones from three suppression-subtractive hybridization libraries yielded a total of 581 unique sequences. Screening a macroarray with a subset of those sequences revealed a total of 26 genes that exhibited differential mRNA accumulation: 16 female biased and 10 male biased. Additional analyses revealed that C16ORF62 (an unknown gene) and MALAT1 (a long noncoding RNA) exhibited increased mRNA accumulation at the male producing temperature relative to the female producing temperature during embryonic sexual development. Finally, we identified four genes (C16ORF62, CCT3, MMP2, and NFIB) that exhibited a stage effect and five genes (C16ORF62, CCT3, MMP2, NFIB and NOTCH2) showed a response to E2 exposure.
Here we report a survey of genes identified using patterns of mRNA accumulation during embryonic development in a turtle with TSD. Many previous studies have focused on examining the turtle orthologs of genes involved in mammalian development. Although valuable, the limitations of this approach are exemplified by our identification of two genes (MALAT1 and C16ORF62) that are sexually dimorphic during embryonic development. MALAT1 is a noncoding RNA that has not been implicated in sexual differentiation in other vertebrates and C16ORF62 has an unknown function. Our results revealed genes that are candidates for having roles in turtle embryonic development, including TSD, and highlight the need to expand our search parameters beyond protein-coding genes.
- Suppression Subtraction Hybridization
- Sexual Development
- mRNA Accumulation
- Suppression Subtraction Hybridization Library
- Subtraction Library
A general overview of sexually dimorphic gene expression in turtles with TSD
late in TSP
Many reptilian taxa, including the majority of turtle species, exhibit TSD [20, 21]. Incubation temperature is the initial cue for sexual development in TSD, in contrast to genetic sex determination (GSD) that is evident in a number of vertebrate groups such as amphibians, snakes, birds, and mammals . GSD is best characterized in therian mammals and is initiated by the SRY gene, located on the Y chromosome, which causes organisms expressing the gene to develop as males . SRY orthologs have not been identified in other groups of vertebrates regardless of whether they exhibit TSD or GSD, suggesting that SRY is an innovation unique to therian mammals. In fact, only one other unique “trigger” gene for sexual development has been identified in a vertebrate taxon, a fish with GSD (medaka; see [24, 25]). It is unclear whether a trigger gene exists in organisms that exhibit TSD, since there are several models that can explain TSD. For example, TSD may reflect regulation of a trigger gene (or set of trigger genes) by incubation temperature, it may reflect the impact of temperature upon the activity of specific enzymes that have a role in signaling, or it may reflect a combination of both phenomena [9, 12, 13]. Regardless, it is clear based upon the studies in organisms with known triggers that a gene homologous to a known trigger in other organisms does not regulate TSD in turtles .
Although trigger gene(s) are not conserved, if one or more are even present, a number of genes involved in gonadal differentiation and other aspects of sexual development are conserved among vertebrates, including organisms with different sex determining systems [27, 28]. A number of orthologs of genes first identified in mammals have been identified and characterized in different vertebrate groups [29–33], including turtles (Table 1). Studies focused on orthologs of genes known to play a role in mammalian sexual development have provided valuable information, although this approach has limits for broader investigations. One setback to an ortholog approach is that the complete set of genes involved in mammalian GSD remains unknown and therefore, it is restricting in finding novel genes or pathways specific to a new system. A complementary approach is to identify candidate genes using patterns of expression rather than orthology. This raises the question of the most appropriate tissue to assay gene expression, since sexually dimorphic gene expression has been noted in multiple tissues [13, 34–36]. Thus, focusing on specific tissues may result in the omission of critical candidate genes. The use of whole embryos is an unbiased strategy for finding candidate genes for sex determination, but it also allows for the identification of genes in other pathways that occur within the same time points. Also, temperature and hormones do not exclusively affect sex determination, and therefore novel genes can also be identified that are related to temperature and the effects of hormonal exposure.
The goal of this study is to use the red-eared slider turtle (Trachemys scripta), a turtle with TSD, to identify candidate genes within the embryonic development pathway. Genes that exhibited a sexually dimorphic expression pattern during the temperature-sensitive period (TSP), which contains the critical stage for commitment to a specific sex , were targeted. The approach we used will also reveal genes that exhibit a differential expression pattern due to temperature or hormonal differences that are not directly related to sexual development. To accomplish this search we identified genes that show increased mRNA accumulation under either the male or female producing temperatures as well as genes that exhibit response to estrogen exposure. Since we also examined embryos at various times during the TSP it was also possible to identify stage effects upon mRNA accumulation during the TSP. These different searches allow for a broader investigation for not only sexually specific factors but also for general developmental genes previously unknown within the time points of the TSP. To accomplish this, we produced three subtraction libraries. Two of these libraries were enriched for genes that show higher mRNA accumulation during the TSP in one specific temperature regime (i.e., genes that exhibit greater mRNA accumulation at the female-producing temperature [31°C] than at the male-producing temperature [26°C] and vice versa). The third library was enriched for genes that show increased mRNA accumulation during the TSP in embryos produced at the male-producing temperature after treatment with exogenous estrogen, ultimately producing the female phenotype. Subsets of the cDNAs from these libraries were examined more thoroughly by macroarray hybridization and semi-quantitative PCR or quantitative real-time PCR. This approach has the potential to identify developmental genes that exhibit differential expression when red-eared slider turtles are exposed to different temperatures and hormonal conditions during the TSP without making a priori assumptions about the identity of the genes.
Suppression subtraction hybridization (SSH) libraries
SSH was used to construct libraries enriched for cDNAs that correspond to mRNAs that exhibit different levels of accumulation during the TSP. Three subtracted cDNA libraries were constructed: one enriched for mRNAs that accumulate at higher levels at the female-producing temperature than the male-producing temperature (hereafter called the “female library”); another enriched for mRNAs that accumulate at higher levels at the male-producing temperature than the female-producing temperature (hereafter called the “male library”); and a third enriched for mRNAs that accumulate at higher levels at the male-producing temperature with exogenous estradiol-17β (sufficient for sex reversal) than in similar embryos treated with the vehicle alone (hereafter called the “E2 library”). A total of 767 sequences were obtained and were previously deposited in dbEST (FG341000:FG341832). The SSH expressed sequence tags (ESTs) were processed as described [38, 39], yielding a total of 581 contigs and singletons (unigenes) after assembly using CAP3 . The results for the homology searches can be found in Additional file 1: Homology search for subtraction libraries.
SSH libraries typically contain some housekeeping genes [41–43] since it is difficult to completely eliminate genes that do not exhibit differential expression between the two experimental conditions. The quality of SSH libraries can be assessed by examining the proportion of broadly expressed genes, although the most appropriate sets of genes to view as “housekeeping” can be problematic to define precisely . However, analyses of GC-content  provide a line of evidence that the proportion of cDNAs that correspond to housekeeping gene transcripts is greatly reduced in our SSH libraries; housekeeping genes tend to have a higher GC-content than genes that exhibit lower levels of expression [45, 46] and it provides an efficient and practical method to define housekeeping genes that avoids conflicts among the available lists of housekeeping genes. We found that the GC-content of the turtle transcripts was lower than expected for other reptilian EST efforts . Thus, the SSH method did appear to enrich for genes with lower levels of mRNA accumulation despite being unable, as expected, to eliminate all housekeeping gene cDNAs.
Genes found in the SSH libraries
GeneMerge was used to test for over-represented GO (Gene Ontology) terms signifying biological processes from genes with human homologs found in all three SSH libraries  (Additional file 2: GeneMerge for biological processes). A total of 34 over-represented GO terms were significant (p < 0.05) and they represent a broad range of biological processes. A few umbrella categories that include a number of over-represented GO terms are anatomical structure morphogenesis (GO:0009653; includes face morphogenesis [GO:0060325] and skeletal system morphogenesis [GO:0048705]), cellular processing (GO:0009987; includes ribosomal small subunit biogenesis [GO:0042274], T cell differentiation in the thymus [GO:0033077], cellular membrane organization [GO:0016044], DNA packaging [GO:0006323], regulation of cell cycle [GO:0051726], negative regulation of apoptosis [GO:0043066]), and metabolic processing (GO:0044267; includes translation [GO:0006412], transcription [GO:0006350], protein folding [GO:0006457], translational initiation [GO:0006413], and translational elongation [GO:0006414]). These categories show that the genes found in the SSH libraries involve active cell differentiation and processing, as expected for mRNAs expressed in developing embryos.
Temperature responsive genes found in SSH
Library where found
Redundancy within library
Differential expression revealed by macroarray analyses
Genes that are found to be sexually dimorphic are not automatically considered to be part of TSD since differing temperatures and hormones can affect more than just sexual development. The genes that emerged from the macroarray as being differentially expressed have a mixture of biological roles in humans based on DAVID. Some genes overlap in their biological roles while others have more distinct roles. For example, 10 genes (GTPBP4, HSP90AA1, ARID4A, RAN, HBZ, SERPINA3, BRIP1, NFE2L1, CDK6, and NFIB) are involved in the regulation of metabolic processing and 6 genes (GTPBP4, BRIP1, RAN, KATNA1, CDK6, and NFIB) are involved in cell division and proliferation (5 of the 6 in the later categories are also found in the earlier category of regulation of metabolic processing). Moreover, AFP and LAPTM4A are involved in reproduction and transport, respectively. Though these genes are not directly related to human sexual development, they provide us with a clue as to which types of processes occur during turtle development under different temperatures and hormone exposures.
Semi quantitative (semiQ) PCR validation of TSD candidate genes
Only the turtle ortholog of C16ORF62, a gene of unknown function, showed evidence of sexually dimorphic expression at stage 17. The transcript of this gene showed greater accumulation in males than in females (~2.5 fold increase), the same trend that was evident in the macroarray results. It is conserved in mammals, birds, fish, insects, nematodes, and plants, its GO term is Integral to membrane, and it is found in a wide variety of adult and embryonic tissues in mammals [57, 58]. Since C16ORF62 is a gene of unknown function and it was found to be sexually dimorphic in turtle embryos during the TSP, it represents truly a novel candidate for a gene involved in turtle development, potentially including TSD; and potentially in sexual development in other vertebrates as well.
Four genes exhibited increased accumulation during stage 17 relative to stage 14 in embryos incubated at the male producing temperature (26°C). This stage effect was evident for C16ORF62, CCT3, MMP2, and NFIB (Figure 3B); the most striking is a ~8-fold increase in mRNA accumulation between stage 14 and 17 for C16ORF62. The others showed a range of relative increase in mRNA accumulation of 1.2-fold to 4.8-fold with stage progression. CCT3, MMP2, and NFIB have been implicated in gonad development (MMP2: ) or other aspects of development (NFIB and CCT3: [54, 55]) in mammals. When this information is combined with our observation that turtle orthologs exhibited increased mRNA accumulation as development proceeded from stage 14 to stage 17 (early in TSP) it is reasonable to speculate that these genes play a role in turtle development, potentially sexual development in the case of MMP2. Though NOTCH2 does not show a stage effect for turtle it has been previously seen as a developmental gene in neuronal development of a mammal and might have a more significant affect on turtle development during different stages or the stage differences are too low to identify with this study .
All five genes show a rapid (6 hours) response to E2 exposure. C16ORF62, CCT3, and MMP2 all show a downregulation of mRNA expression and NFIB and NOTCH2 show an upregulation. Four of the genes (CCT3, MMP2 NFIB, and NOTCH2) exhibited similar mRNA accumulation both 6 hours and 24 hours after E2 exposure; accumulation of the C16ORF62 mRNA almost returned to pre-exposure levels after 24 hours. MMP2 is affected by E2 based upon previous studies whereas the impact of E2 exposure upon mRNA accumulation for other genes has not been examined. Mahmoodzadeh et al.  showed that E2 inhibits MMP2 gene expression in rat fibroblasts and those results corroborate our findings for MMP2’s involvement with E2.
Expression of a long noncoding RNA (ncRNA) is sexually-dimorphic
A number of cDNAs on the macroarray (61) could not be identified using BLASTX, suggesting that they correspond either to cDNAs for which only untranslated region was included in the EST read or noncoding RNAs (ncRNA). To identify some of these cDNAs we conducted BLASTN searches and revealed that one of the cDNAs that exhibits sexual dimorphism is a ncRNA, MALAT1.
MALAT1 is a long (~7 kb) ncRNA that undergoes a cleavage that produces two RNAs, a smaller tRNA-like cytoplasmic RNA (~61nt) and a 6.7 kb RNA that localize to two different subcellular compartments, cytoplasm and nuclear speckles respectively . Characteristically, it has short blocks of high conservation across the entire transcript, especially in 3’ half of RNA, and lacks repetitive elements except for a SINE and LINE element near its 5’ end . The smaller (~61nt) transcript generated by cleavage is highly conserved across many species, including mouse, human, dog, lizard, frog, and stickleback. It has not yet been found in any of the four available sequenced bird genomes (chicken, turkey, zebrafinch, and duck). There are two possible explanations: (1) birds lost the locus or (2) it is located in a region that was consistently underrepresented in all of the available avian genome assemblies [61–64]. Both explanations are plausible since gene loss is known to be an important process during evolution  but the avian genome assemblies (like other vertebrate genome assemblies) are incomplete.
MALAT1 shows a broad distribution of expression in normal human and mouse tissues but its misregulation is correlated with the progression of cancers and it is upregulated in many human carcinomas [66–69]. More importantly for this study, MALAT1 accumulation is higher in adult mammalian ovaries than adult testes [60, 70]. However, the pattern of differential expression for MALAT1 in adult mammalian gonads is distinct from the pattern we observed using the macroarray assay, in which the mRNA accumulation appeared 1.6-fold higher in whole male turtle embryos. Since MALAT1 is a ncRNA that shows dimorphic expression in the TSP we felt it was an excellent candidate for a gene involved in TSD so we used quantitative real-time PCR (qRT-PCR) to verify the pattern of expression suggested by the macroarray.
Here we reported a survey of genes identified based upon their patterns of mRNA accumulation during embryonic development in the Red-eared slider turtle. We used a non-ortholog based strategy and identified four genes that exhibited increased mRNA expression as development proceeded from stage 14 to stage 17, a set of genes that responded to E2 exposure, and two genes (MALAT1 and C16ORF62) that show greater accumulation at the male producing temperature than at the female producing temperature. This survey focused on changes in mRNA accumulation in whole embryos. Thus, it remains possible that some or all of these genes exhibit even more strongly dimorphic expression in specific tissues (e.g., the developing gonad or brain). Moreover, the genes we identified are likely to be significant since screening for differential expression at the whole embryo level is expected to be a conservative way to examine gene expression during development.
MMP s (matrix metalloproteinases) are involved in the breakdown of extracellular matrices in physiological processes, including cancer . MMP2 is sexually dimorphic in developing male mice because it functions as a paracrine death factor in Müllerian duct regression downstream of the MIS cascade . In addition, Kim et al. 2008  found MMP2 to be sexually dimorphic and regulated by testosterone in songbirds in relation to the vocal control center during adult neurogenesis. Furthermore, estrogen affects the MMP pathway in humans by increasing MMP2 enzymatic activity . It is unclear if the increase in MMP2 is through an increase in mRNA accumulation or through other mechanisms, such as binding affinity changes. Though MMP2 was not found to be sexually dimorphic in turtles it was found to be inhibited by E2, the opposite of the regulation in mammals but potentially similar to birds . Together with the prior knowledge of its involvement in mammalian and avian development, MMP2 is a novel candidate gene for development in the turtle.
ncRNAs are believed to play a large number of biological roles (reviewed in ), but their role in development remains poorly characterized . Although there is some evidence that ncRNAs have roles in sexual development in both mammals  and birds , this is the first evidence that a ncRNA may have a role in sexual development for an organism with TSD. This hypothesis is corroborated by the fact that MALAT1 exhibits differential expression in mammalian gonads (expression is higher in adult ovaries than in testes). However, the pattern of sexual dimorphism reported for mammals is distinct from that evident in turtles (where the RNA accumulation is higher in male embryos than in female embryos). Our findings highlight the importance of examining ncRNAs when investigating vertebrate development in general and sexual development specifically.
Little is known about the genes involved in turtle development, including the processes related to TSD, regardless of whether they are protein-coding genes or ncRNAs. Much of the information available is focused on the examination of the turtle orthologs of genes involved in mammalian development. The limitations on this type of an examination are exemplified by our novel results of two genes (MALAT1 and C16ORF62) that show sexual dimorphism during TSP. These two types of genes are novel to turtle sexual development because one is a ncRNA (MALAT1), not typically found in sexual differentiation in other vertebrates, and the other (C16ORF62) is of unknown function. Our results highlight the need to diverge from focusing only on protein-coding genes when looking for developmental genes and to expand into the more diverse world of RNA in general, specifically including ncRNAs.
Incubation and experimental manipulations
Freshly laid Trachemys scripta eggs (500) were purchased from Kliebert Turtle Farms in Hammond, Louisiana in 2004 and 2006. They were kept at room temperature for less than 48 hours until they were established as viable by candling. Those viable were randomly separated equally into four experimental groups in containers with moistened vermiculite (1:1 vermiculite to water). The experimental groups were the female producing temperature of 31°C (hereafter called female), the male producing temperature of 26°C (hereafter called male), 26°C painted with exogenous estradiol-17β (E2) in 1 μg/μL in 95% ethanol (non-denatured) (phenotypically female; hereafter called E2), and 26°C painted with exogenous 95% ethanol (non-denatured) as the vehicle control (phenotypically male; hereafter called vehicle). Application of E2 and vehicle occurred at stage 14. The egg boxes were rotated daily within the incubators and we checked a random selection of eggs periodically to determine their developmental stage using the staging guidelines suggested by Yntema . The temperature was monitored daily with HOBO data loggers and in-incubator thermometers. Sex was determined for each experimental group through a visual inspection by two independent researchers at hatching (gonads are visually distinct at hatching but not before) of 10 embryos per experimental group for relevant gross anatomy. All experiments complied with the appropriate ethical guidelines.
Isolation of RNA
Whole embryos were taken between stages 17 and 20 from each experimental group and quickly frozen in liquid nitrogen and stored at −80°C. A subset of whole embryos was also collected at 0, 6, and 24 hours after E2 and vehicle application at stage 14. An average of 5 embryos was collected per stage and experimental group. Conducting the E2 time trial at stage 14 before the TSP removes E2 effects within gonadal differentiation and leaves just those from the trial. Total RNA was extracted from each embryo by homogenization in Tri-Reagent (TRIzol, Sigma USA), followed by extraction in chloroform, and precipitation in isopropanol according to Sambrook and Russell . Total RNA yield and quality were assessed with the ND-1000 Nanodrop spectrophotometer (NanoDrop Technologies, Thermo Fisher Scientific, Wilmington, DE 19810, USA), and the integrity was verified by electrophoresis on a 1% agarose gel.
Suppression subtractive hybridization (SSH)
Three libraries were selectively induced for female against male, male against female, and E2 against vehicle. Testers and drivers were made from pooled RNA from stages 17–20 (two individuals per stage) from each experimental group from 2004 (stated above). cDNA synthesis was performed with the BD SMART™ PCR cDNA Synthesis Kit (Clontech, Mountain View, CA) according to the manufacturer’s protocol. Three subtraction libraries were constructed with the Clontech PCR–SelectTM cDNA Subtraction Kit (Clontech, Mountain View, CA) according to the manufacturer’s protocol except a PEG (PolyEthylene Glycol) precipitation followed by an ethanol wash was used to purify the PCR products after cDNA synthesis instead of the column chromatography. The resulting cDNA was ligated into a pGEM-T Easy vector and transformed into E. cloni® 10 G electrocompetent cells (Lucigen) by the manufacturer’s protocol. Individual colonies were picked and stored in 96-well plates with 50% glycerol at −80°C. Plasmid inserts were purified using a modified 96-well Perfectprep® Plasmid protocol (5Prime, Gaithersburg, MD), according to Sambrook and Russell  or a TempliPhi Amplification kit (as recommended by manufacturer; GE Healthcare). Single-pass sequencing was conducted on an ABI Prism™ 3100-Avant genetic analyzer (PE Applied Biosystems) using the ABI BigDye® Terminator v.3.1 chemistry.
Analysis of SSH results
Sequences from the libraries with redundancy were aligned and all sequences (both individual and aligned) were edited in Sequencher™ 4.1 (Gene Codes Corp.). Sequences were put into FASTA format and run in GOanna from the AgBase v.2.0 database to determine the top Blast hit for each sequence and to simultaneously determine GO terms for each hit . GOanna uses BLASTX (determines gene products from sequences) therefore any sequences that did not have a hit were run through BLASTN (determines all aspects of RNA transcripts including untranslated regions and non-protein coding RNAs) on the NCBI server .
GeneMerge categorized the human homologs of genes that resulted from all three SHH libraries with over-represented GO (Gene Ontology) terms from the biological processes category given a human background set of genes . The significance cut-off was set at P < 0.05.
The genes that have human homologs from the SSH libraries were clustered into functionally related groups within the subset of biological processes by DAVID v.6.7 [48, 49]. The DAVID tool for functional annotation clustering uses GO terms and the term enrichment score was used at high stringency (based on P < 0.05).
Macroarray preparation and analyses
Three hundred and seventy four clones (including 322 known clones [those with BLAST hits] and 61 unknown clones) that were obtained from all three SSH libraries (discussed above) were spotted onto membranes (Pall Biodyne B Nylon, Nunc) using 100 nanoliter pins on a Biomek 2000 (Beckman Coulter, USA). Positive controls (Arabidopsis thaliana RCA [X14212], CAB [X56062], and RBCL [U91966]) and negative controls were also spotted onto the membranes. All samples were spotted in duplicate with four replicates per experimental group (female and male). Total RNA was collected (as stated above) from 2 embryos from stage 17 and 2 embryos from stage 19. All collected embryos from the male experimental group in 2004 were pooled as well as for the female experimental group. Pooled total RNA was mixed with control cDNAs and then reverse transcribed before labelling with α33P-dATP as described in Blum et al. . After hybridization , the membranes were rinsed and exposed to a phosphor imager and scanned using a Molecular Devices Typhoon Scanner.
Signal intensities were quantified using ImageQuant 5.1 (Amersham Biosciences UK Limited, Amersham Place Little Chalfont Buckinghamshire England) and intensity differences were calculated as described by Helbing et al. . Briefly, image data were converted to a standard 8-bit TIFF file before accounting for signal saturation. Positive controls used to standardize across arrays for hybridization efficiency were chosen based on their coefficient of variation (cv) across all 8 arrays (the cv did not exceed 0.2). Each array was normalized to the geometric mean for all positive controls chosen for that array . The non-signal background was determined using the median intensity value plus one standard deviation for the blank positions and negative controls, and the maximum value for all arrays was set as the “no signal” value. All values across all arrays equal to or below that number is considered zero. After normalizing across arrays and discarding values equal to or below the floor value, the average of the duplicates within arrays, the estimated standard deviation across replicate arrays, the median across replicate arrays, and the fold change between treatment groups were calculated. If the standard deviation across replicates was greater than or equal to 2 then those spots were not reliable for further examination. Based upon Helbing et al. , we considered a 1.5-fold difference in signal intensity relative to control treatments as sufficient in the array experiments based on the detection limitations of the cDNA array analyses. Those genes that exhibited differential expression were used for further examination.
Semi-quantitative PCR preparation and analysis
Embryos were collected and total RNA was extracted as stated above from the 31°C experimental group at stage 17, the 26°C experimental group at stage 17, and from the E2 time trial experimental group at stage 14 from 2006. Four embryos from each group were pooled and cDNA was made using Invitrogen’s Superscript III Reverse Transcription kit (Invitrogen, USA) following the manufacturer’s instructions.
Semi-quantitative and quantitative real-time PCR primers, and optimal conditions for semi-quantitative PCR
Annealing Temperature (°C)
Quantitative Real-time PCR (qRT-PCR) preparation and analysis
Five whole embryos were collected from stage 17 and 5 from stage 19 from two experimental groups, male and female, and total RNA was extracted as stated above. cDNA was generated using ImProm-II™ Reverse Transcriptase and random primers following the manufacturer’s instructions. Relative gene expression levels were quantified using an ABI StepOnePlus™ Real-time PCR cycler (STepOne™ Software v2.1) with the following cycling parameters: initial denaturing for 10 min at 95°C, followed by 40 cycles of 35 s at 95°C, 30 s at 60°C, and 30 s at 72°C. The final cycle was followed by a melting curve analysis to verify the amplification of a single product in each well. Specificities of all primer pairs were also verified by sequencing PCR products. Repeating the above procedures on RNA samples (prior to reverse transcription) verified that no products were amplified from contaminating genomic DNA. All samples were run in duplicate and included 3.75 μL of a 1:100 diluted sample, 1 μM of each primer, and 2x SYBR Green Master Mix (Applied Biosystems) in a total of 15 μL. PCR efficiencies were calculated from a gene-specific standard curve from a 10-fold dilution series. Relative transcript abundance was normalized to the expression of PP1 by using the relative standard curve method . To determine if expression differed between experimental groups a two-tailed Student’s t-test and a standard error analysis were performed. Primers used to assay gene expression were designed using Primer 3  and Amplify  (Table 3).
We are grateful to Rebecca Kimball, Gordon Burleigh, Teresa Bryan, Marta Wayne, Karen Koch and the members of the Braun-Kimball lab group for helpful discussions. This work was facilitated by National Science Foundation grants to ELB and collaborators (DEB-0228682 and DUE-0920151) as well as an American Society of Ichthyology and Herpetology grant and a Sigma Xi grant to JLC. Publication of this article was funded in part by the University of Florida Open-Access Publishing Fund.
- Francisco-Morcillo J, Hidalgo-Sánchez M, Martín-Partido G: Spatial and temporal patterns of proliferation and differentiation in the developing turtle eye. Brain Res. 2006, 1103: 32-48. 10.1016/j.brainres.2006.05.052.View ArticlePubMedGoogle Scholar
- Mehta V, Sernagor E: Early neural activity and dendritic growth in turtle retinal ganglion cells. Eur J Neurosci. 2006, 24: 773-786. 10.1111/j.1460-9568.2006.04933.x.View ArticlePubMedGoogle Scholar
- Bernier B, Bar I, Pieau C, Lambert De Rouvroit C, Goffinet A: Reelin mRNA expression during embryonic brain development in the turtle Emys orbicularis. J Comp Neurol. 1999, 413: 463-479. 10.1002/(SICI)1096-9861(19991025)413:3<463::AID-CNE8>3.0.CO;2-F.View ArticlePubMedGoogle Scholar
- Metin C, Alvarez C, Moudoux D, Vitalis T, Pieau C, Molnar Z: Conserved pattern of tangentieal neuronal migration during forebrain development. Development. 2007, 134: 2815-2827. 10.1242/dev.02869.View ArticlePubMedGoogle Scholar
- Kuratani S, Kuraku S, Nagashima H: Evolutionary developmental perspective for the origin of turtles: the folding theory for the shell based on the developmental nature of the carapacial ridge. Evol Dev. 2011, 13: 1-14. 10.1111/j.1525-142X.2010.00451.x.View ArticlePubMedGoogle Scholar
- Ramsey M, Crews D: Steroid signaling and temperature-dependent sex determination—Reviewing the evidence for early action of estrogen during ovarian determination in turtles. Semin Cell Dev Biol. 2009, 20: 283-292. 10.1016/j.semcdb.2008.10.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Fleming A, Wibbels T, Skipper JK, Crews D: Developmental expression of steroidogenic factor 1 in a turtle with temperature-dependent sex determination. Gen Comp Endocrinol. 1999, 116: 336-346. 10.1006/gcen.1999.7360.View ArticlePubMedGoogle Scholar
- Ramsey M, Shoemaker C, Crews D: Gonadal expression of Sf1 and aromatase during sex determination in the red-eared slider turtle (Trachemys scripta), a reptile with temperature-dependent sex determination. Differentiation. 2007, 75: 978-991.View ArticlePubMedGoogle Scholar
- Valenzuela N: Relic thermosensitive gene expression in a turtle with genotypic sex determination. Evolution. 2008, 62: 234-240.View ArticlePubMedGoogle Scholar
- Schmahl J, Yao HH, Pierucci-Alves F, Capel B: Colocalization of WT1 and cell proliferation reveals conserved mechanisms in temperature-dependent sex determination. Genesis. 2003, 35: 193-201. 10.1002/gene.10176.PubMed CentralView ArticlePubMedGoogle Scholar
- Spotila LD, Spotila JR, Hall SE: Sequence and expression analysis of WT1 and Sox9 in the red-eared slider turtle, Trachemys scripta. J Exp Zool. 1998, 281: 417-427. 10.1002/(SICI)1097-010X(19980801)281:5<417::AID-JEZ7>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Shoemaker CM, Queen J, Crews D: Response of candidate sex-determining genes to changes in temperature reveals their involvement in the molecular network underlying temperature-dependent sex determination. Mol Endocrinol. 2007, 21: 2750-2763. 10.1210/me.2007-0263.View ArticlePubMedGoogle Scholar
- Ramsey M, Crews D: Adrenal-kidney-gonad complex measurements may not predict gonad-specific changes in gene expression patterns during temperature-dependent sex determination in the red-eared slider turtle (Trachemys scripta elegans). J Exp Zool A Ecol Genet Physiol. 2007, 307: 463-470.View ArticlePubMedGoogle Scholar
- 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. 2007, 236: 1055-1063. 10.1002/dvdy.21096.View ArticlePubMedGoogle Scholar
- Murdock C, Wibbels T: Expression of Dmrt1 in a turtle with temperature-dependent sex determination. Cytogenet Genome Res. 2003, 101: 302-308. 10.1159/000074353.View ArticlePubMedGoogle Scholar
- Torres Maldonado LC, Landa Piedra A, Moreno Mendoza N, Marmolejo Valencia A, Meza Martinez A, Merchant Larios H: 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. 10.1016/S0016-6480(02)00511-7.View ArticlePubMedGoogle Scholar
- Murdock C, Wibbels T: Cloning and expression of aromatase in a turtle with temperature-dependent sex determination. Gen Comp Endocrinol. 2003, 130: 109-119. 10.1016/S0016-6480(02)00573-7.View ArticlePubMedGoogle Scholar
- Takada S, DiNapoli L, Capel B, Koopman P: Sox8 is expressed at similar levels in gonads of both sexes during the sex determining period in turtles. Dev Dyn. 2004, 231: 387-395. 10.1002/dvdy.20132.View ArticlePubMedGoogle Scholar
- Smith CA, Shoemaker CM, Roeszler KN, Queen J, Crews D, Sinclair AH: Cloning and expression of R-Spondin1 in different vertebrates suggests a conserved role in ovarian development. BMC Dev Biol. 2008, 8: 72-10.1186/1471-213X-8-72.PubMed CentralView ArticlePubMedGoogle Scholar
- Crews D, Bergeron JM, Bull JJ, Flores D, Tousignant A, Skipper JK, Wibbels T: Temperature-dependent sex determination in reptiles: proximate mechanisms, ultimate outcomes, and practical applications. Dev Genet. 1994, 15: 297-312. 10.1002/dvg.1020150310.View ArticlePubMedGoogle Scholar
- Janzen FJ, Phillips PC: Exploring the evolution of environmental sex determination, especially in reptiles. J Evol Biol. 2006, 19: 1775-1784. 10.1111/j.1420-9101.2006.01138.x.View ArticlePubMedGoogle Scholar
- Ezaz T, Stiglec R, Veyrunes F, Marshall Graves JA: Relationships between vertebrate ZW and XY sex chromosome systems. Curr Biol. 2006, 16: 736-743. 10.1016/j.cub.2006.08.021.View ArticleGoogle Scholar
- Sinclair AH, Berta P, Palmer MS, Hawkins JR, Griffiths BL, Smith MJ, Foster JW, Frischauf AM, Lovell-Badge R, Goodfellow PN: A gene from the human sex-determining region encodes a protein with homology to a conserved DNA-binding motif. Nature. 1990, 346: 240-244. 10.1038/346240a0.View ArticlePubMedGoogle Scholar
- Volff JN, Kondo M, Schartl M: Medaka dmY/dmrt1Y is not the universal primary sex-determining gene in fish. Trends Genet. 2003, 19: 196-199. 10.1016/S0168-9525(03)00051-9.View ArticlePubMedGoogle Scholar
- Matsuda M, Shinomiya A, Kinoshita M, Suzuki A, Kobayashi T, Paul-Prasanth B, Lau EL, Hamaguchi S, Sakaizumi M, Nagahama Y: DMY gene induces male development in genetically female (XX) medaka fish. Proc Natl Acad Sci U S A. 2007, 104: 3865-3870. 10.1073/pnas.0611707104.PubMed CentralView ArticlePubMedGoogle Scholar
- Herpin A, Braasch I, Kraeussling M, Schmidt C, Thoma EC, Nakamura S, Tanaka M, Schartl M: Transcriptional rewiring of the sex determining dmrt1 gene duplicate by transposable elements. PLoS Genet. 2010, 6: e1000844-10.1371/journal.pgen.1000844.PubMed CentralView ArticlePubMedGoogle Scholar
- Johnston CM, Barnett M, Sharpe PT: The molecular biology of temperature-dependent sex determination. Philos Trans R Soc Lond B Biol Sci. 1995, 350: 297-303. 10.1098/rstb.1995.0165. discussion 303–294View ArticlePubMedGoogle Scholar
- Western PS, Sinclair AH: Sex, genes, and heat: triggers of diversity. J Exp Zool. 2001, 290 (6): 624-631. 10.1002/jez.1113.View ArticlePubMedGoogle Scholar
- Clinton M: Sex determination and gonadal development: a bird's eye view. J Exp Zool. 1998, 281: 457-465. 10.1002/(SICI)1097-010X(19980801)281:5<457::AID-JEZ10>3.0.CO;2-6.View ArticlePubMedGoogle Scholar
- Guan G, Kobayashi T, Nagahama Y: Sexually dimorphic expression of two types of DM (Doublesex/Mab-3)-domain genes in a teleost fish, the Tilapia (Oreochromis niloticus). Biochem Biophys Res Commun. 2000, 272: 662-666. 10.1006/bbrc.2000.2840.View ArticlePubMedGoogle Scholar
- Western PS, Harry JL, Marshall Graves JA, Sinclair AH: Temperature-dependent sex determination in the American alligator: expression of SF1, WT1 and DAX1 during gonadogenesis. Gene. 2000, 241: 223-232. 10.1016/S0378-1119(99)00466-7.View ArticlePubMedGoogle Scholar
- Smith CA, Sinclair AH: Sex determination in the chicken embryo. J Exp Zool. 2001, 290: 691-699. 10.1002/jez.1119.View ArticlePubMedGoogle Scholar
- Shibata K, Takase M, Nakamura M: The Dmrt1 expression in sex-reversed gonads of amphibians. Gen Comp Endocrinol. 2002, 127: 232-241. 10.1016/S0016-6480(02)00039-4.View ArticlePubMedGoogle Scholar
- Dewing P, Shi T, Horvath S, Eric Vilain E: Sexually dimorphic gene expression in mouse brain precedes gonadal differentiation. Brain Res Mol Brain Res. 2003, 118: 82-90.View ArticlePubMedGoogle Scholar
- Lee SI, Lee WK, Shin JH, Han BK, Moon S, Cho S, Park T, Kim H, Han JY: Sexually dimorphic gene expression in the chick brain before gonadal differentiation. Poult Sci. 2009, 88: 1003-1015. 10.3382/ps.2008-00197.View ArticlePubMedGoogle Scholar
- Jänne M, Hogeveen KN, Deol HK, Hammond GL: Expression and regulation of human sex hormone-binding globulin transgenes in mice during development. Endocrinology. 1998, 140: 4166-4174.Google Scholar
- Wibbels T, Bull JJ, Crews D: Chronology and morphology of temperature-dependent sex determination. J Exp Zool. 1991, 260: 371-381. 10.1002/jez.1402600311.View ArticlePubMedGoogle Scholar
- Chojnowski JL, Franklin J, Katsu Y, Iguchi T, Guillette LJ, Kimball RT, Braun EL: Patterns of vertebrate isochore evolution revealed by comparison of expressed mammalian, avian, and crocodilian genes. J Mol Evol. 2007, 65: 259-266. 10.1007/s00239-007-9003-2.View ArticlePubMedGoogle Scholar
- Chojnowski JL, Braun EL: Turtle isochore structure is intermediate between amphibians and other amniotes. Integr Comp Biol. 2008, 48: 454-462. 10.1093/icb/icn062.View ArticlePubMedGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle Scholar
- Diatchenko L, Lau YF, Campbell AP, Chenchik A, Moqadam F, Huang B, Lukyanov S, Lukyanov K, Gurskaya N, Sverdlov ED: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Natl Acad Sci U S A. 1996, 93: 6025-6030. 10.1073/pnas.93.12.6025.PubMed CentralView ArticlePubMedGoogle Scholar
- Luo MJ, Lai MD: Identification of differentially expressed genes in normal mucosa, adenoma and adenocarcinoma of colon by SSH. World J Gastroenterol. 2001, 7: 726-731.PubMed CentralPubMedGoogle Scholar
- Chen AQ, Wang ZG, Xu ZR, Yu SD, Yang ZG: Analysis of gene expression in granulosa cells of ovine antral growing follicles using suppressive subtractive hybridization. Anim Reprod Sci. 2009, 115: 39-48. 10.1016/j.anireprosci.2008.10.022.View ArticlePubMedGoogle Scholar
- De Ferrari L, Aitken S: Mining housekeeping genes with a Naive Bayes classifier. BMC Genomics. 2006, 7: 277-291. 10.1186/1471-2164-7-277.PubMed CentralView ArticlePubMedGoogle Scholar
- Kudla G, Lipinski L, Caffin F, Helwak A, Zylicz M: High guanine and cytosine content increases mRNA levels in mammalian cells. PLoS Biol. 2006, 4: e180-10.1371/journal.pbio.0040180.PubMed CentralView ArticlePubMedGoogle Scholar
- Arhondakis S, Clay O, Bernardi G: GC level and expression of human coding sequences. Biochem Biophys Res Commun. 2008, 367: 542-545. 10.1016/j.bbrc.2007.12.155.View ArticlePubMedGoogle Scholar
- Castillo-Davis CI, Hartl DL: GeneMerge–post-genomic analysis, data mining, and hypothesis testing. Bioinformatics. 2003, 19: 891-892. 10.1093/bioinformatics/btg114.View ArticlePubMedGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticlePubMedGoogle Scholar
- Roberts LM, Visser JA, Ingraham HA: Involvement of a matrix metalloproteinase in MIS-induced cell death during urogenital development. Development. 2002, 129: 1487-1496.PubMedGoogle Scholar
- Kohno S, Katsu Y, Urushitani H, Ohta Y, Iguchi T, Guillette LJ: Potential contributions of heat shock proteins to temperature-dependent sex determination in the American alligator. Sex Dev. 2010, 4: 73-87. 10.1159/000260374.PubMed CentralView ArticlePubMedGoogle Scholar
- Picard D: Chaperoning steroid hormone action. Trends Endocrinol Metab. 2006, 17: 229-235. 10.1016/j.tem.2006.06.003.View ArticlePubMedGoogle Scholar
- Robinson LL, Sznajder NA, Riley SC, Anderson RA: Matrix metalloproteinases and tissue inhibitors of metalloproteinases in human fetal testis and ovary. Mol Hum Reprod. 2001, 7: 641-648. 10.1093/molehr/7.7.641.View ArticlePubMedGoogle Scholar
- Chaudhry AZ, Lyons GE, Gronostajski RM: Expression patterns of the four nuclear factor I genes during mouse embryogenesis indicate a potential role in development. Dev Dyn. 1997, 208: 313-325. 10.1002/(SICI)1097-0177(199703)208:3<313::AID-AJA3>3.0.CO;2-L.View ArticlePubMedGoogle Scholar
- Walkley NA, Demaine AG, Malik AN: Cloning, structure and mRNA expression of human Cctg, which encodes the chaperonin subunit CCT gamma. Biochem J. 1996, 313 (Pt 2): 381-389.PubMed CentralView ArticlePubMedGoogle Scholar
- Weinmaster G, Roberts V, Lemke G: Notch2: a second mammalian Notch gene. Development. 1992, 116: 931-941.PubMedGoogle Scholar
- Bgee: a dataBase for Gene Expression Evolution.http://bgee.unil.ch/bgee/bgee,
- Mahmoodzadeh S, Dworatzek E, Fritschka S, Pham TH, Regitz-Zagrosek V: 17beta-Estradiol inhibits matrix metalloproteinase-2 transcription via MAP kinase in fibroblasts. Cardiovasc Res. 2010, 85: 719-728. 10.1093/cvr/cvp350.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilusz JE, Freier SM, Spector DL: 3' end processing of a long nuclear-retained noncoding RNA yields a tRNA-like cytoplasmic RNA. Cell. 2008, 135: 919-932. 10.1016/j.cell.2008.10.012.PubMed CentralView ArticlePubMedGoogle Scholar
- Morisson M, Denis M, Milan D, Klopp C, Leroux S, Bardes S, Pitel F, Vignoles F, Gérus M, Fillon V, Douaud M, Vignal A: The chicken RH map: current state of progress and microchromosome mapping. Cytogenet Genome Res. 2007, 117: 14-21. 10.1159/000103160.View ArticlePubMedGoogle Scholar
- Hillier LW, Miller W, Birney E, Warren W, Hardison RC, Ponting CP, Bork P, Burt DW, Groenen MAM, Delany ME, Dodgson JB, Chinwalla AT, Cliften PF, Clifton SW, Delehaunty KD, Fronick C, Fulton RS, Graves TA, Kremitzki C, Layman D, Magrini V, McPherson JD, Miner TL, Minx P, Nash WE, Nhan MN, Nelson JO, Oddy LG, Pohl CS, Randall-Maher J: International Chicken Genome Sequencing Consortium, et al: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.View ArticleGoogle Scholar
- Dalloul RA, Long JA, Zimin AV, Aslam L, Beal K, Ann Blomberg L, Bouffard P, Burt DW, Crasta O, Crooijmans RPMA, Cooper K, Coulombe RA, De S, Delany ME, Dodgson JB, Dong JJ, Evans C, Frederickson KM, Flicek P, Florea L, Folkerts O, Groenen MAM, Harkins TT, Herrero J, Hoffmann S, Megens H-J, Jiang A, de Jong P, Kaiser P, Kim H: Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): Genome assembly and analysis. PLoS Biol. 2010, 8: e1000475-10.1371/journal.pbio.1000475.PubMed CentralView ArticlePubMedGoogle Scholar
- Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Künstner A, Searle S, White S, Vilella AJ, Fairley S, Heger A, Kong L, Ponting CP, Jarvis ED, Mello CV, Minx P, Lovell P, Velho TAF, Ferris M, Balakrishnan CN, Sinha S, Blatti C, London SE, Li Y, Lin Y, George J, Sweedler J, Southey B, Gunaratne P, Watson M: The genome of a songbird. Nature. 2010, 464: 757-762. 10.1038/nature08819.PubMed CentralView ArticlePubMedGoogle Scholar
- Braun EL: Innovation from reduction: gene loss, domain loss and sequence divergence in genome evolution. Appl Bioinformatics. 2003, 2: 13-34.PubMedGoogle Scholar
- Ji P, Diederichs S, Wang W, Boing S, Metzger R, Schneider PM, Tidow N, Brandt B, Buerger H, Bulk E: MALAT-1, a novel noncoding RNA, and thymosin beta4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene. 2003, 22: 8031-8041. 10.1038/sj.onc.1206928.View ArticlePubMedGoogle Scholar
- Guffanti A, Iacono M, Pelucchi P, Kim N, Solda G, Croft LJ, Taft RJ, Rizzi E, Askarian-Amiri M, Bonnal RJ: A transcriptional sketch of a primary human breast cancer by 454 deep sequencing. BMC Genomics. 2009, 10: 163-10.1186/1471-2164-10-163.PubMed CentralView ArticlePubMedGoogle Scholar
- Koshimizu TA, Fujiwara Y, Sakai N, Shibata K, Tsuchiya H: Oxytocin stimulates expression of a noncoding RNA tumor marker in a human neuroblastoma cell line. Life Sci. 2010, 86: 455-460. 10.1016/j.lfs.2010.02.001.View ArticlePubMedGoogle Scholar
- Guo F, Li Y, Liu Y, Wang J, Li Y, Li G: Inhibition of metastasis-associated lung adenocarcinoma transcript 1 in CaSki human cervical cancer cells suppresses cell proliferation and invasion. Acta Biochim Biophys Sin (Shanghai). 2010, 42: 224-229. 10.1093/abbs/gmq008.View ArticleGoogle Scholar
- Hutchinson JN, Ensminger AW, Clemson CM, Lynch CR, Lawrence JB, Chess A: A screen for nuclear transcripts identifies two linked noncoding RNAs associated with SC35 splicing domains. BMC Genomics. 2007, 8: 39-10.1186/1471-2164-8-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Bourboulia D, Stetler-Stevenson WG: Matrix metalloproteinases (MMPs) and tissue inhibitors of metalloproteinases (TIMPs): Positive and negative regulators in tumor cell adhesion. Semin Cancer Biol. 2010, 20: 161-168. 10.1016/j.semcancer.2010.05.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim DH, Lilliehook C, Roides B, Chen Z, Chang M, Mobashery S, Goldman SA: Testosterone-induced matrix metalloproteinase activation is a checkpoint for neuronal addition to the adult songbird brain. J Neurosci. 2008, 28: 208-216.View ArticlePubMedGoogle Scholar
- Grandas OH, Mountain DH, Kirkpatrick SS, Cassada DC, Stevens SL, Freeman MB, Goldman MH: Regulation of vascular smooth muscle cell expression and function of matrix metalloproteinases is mediated by estrogen and progesterone exposure. J Vasc Surg. 2009, 49: 185-191. 10.1016/j.jvs.2008.07.080.View ArticlePubMedGoogle Scholar
- Ponting CP, Oliver PL, Reik W: Evolution and functions of long noncoding RNAs. Cell. 2009, 136: 629-641. 10.1016/j.cell.2009.02.006.View ArticlePubMedGoogle Scholar
- McFarlane L, Wilhelm D: Non-coding RNAs in mammalian sexual development. Sex Dev. 2010, 3: 302-316.View ArticleGoogle Scholar
- Zhao D, McBride D, Nandi S, McQueen HA, McGrew MJ, Hocking PM, Lewis PD, Sang HM, Clinton M: Somatic sex identity is cell autonomous in the chicken. Nature. 2010, 464: 237-242. 10.1038/nature08852.PubMed CentralView ArticlePubMedGoogle Scholar
- Yntema C: A series of stages in the embryonic development of Chelydra serpentina. J Morph. 1968, 125: 219-252. 10.1002/jmor.1051250207.View ArticlePubMedGoogle Scholar
- Sambrook J, Russell D: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York, 3Google Scholar
- McCarthy FM, Wang N, Magee GB, Nanduri B, Lawrence ML, Camon EB, Barrell DG, Hill DP, Dolan ME, Williams WP: AgBase: a functional genomics resource for agriculture. BMC Genomics. 2006, 7: 229-10.1186/1471-2164-7-229.PubMed CentralView ArticlePubMedGoogle Scholar
- Blum JL, Prucha MS, Patel VJ, Denslow ND: Use of cDNA macroarrays and gene profiling for detection of effects of environmental toxicants. Methods Mol Biol. 2008, 410: 43-53. 10.1007/978-1-59745-548-0_4.View ArticlePubMedGoogle Scholar
- Helbing CC, Bailey CM, Ji L, Gunderson MP, Zhang F, Veldhoen N, Skirrow RC, Mu R, Lesperance M, Holcombe GW, Kosian PA, Tietge J, Korte JJ, Degitz SJ: Identification of gene expression indicators for thyroid axis disruption in aXenopus laevismetamorphosis screening assay: Part 1. Effects on the brain. Aquatic Toxicology. 2007, 82: 227-241. 10.1016/j.aquatox.2007.02.013.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515.View ArticlePubMedGoogle Scholar
- ,: Primer 3.http://frodo.wi.mit.edu,
- ImageJ. Wayne Rasband, National Institutes of Health, Bethesda, MD,http://rsb.info.nih.gov/ij,
- Larionov A, Krause A, Miller W: A standard curve based method for relative real time PCR data processing. BMC Bioinformatics. 2005, 6: 62-10.1186/1471-2105-6-62.PubMed CentralView ArticlePubMedGoogle Scholar
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.