text:
MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination
Srihari Radhakrishnan1,4,
Robert Literman2,4,
Beatriz Mizoguchi3,4 &
Nicole Valenzuela ORCID: orcid.org/0000-0003-1148-631X4
Epigenetics & Chromatin volume 10, Article number: 28 (2017) Cite this article
DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.
Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.
Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.
Epigenetic modifications are heritable changes to the DNA that do not change the nucleotide sequence. Among them, DNA methylation is a biochemical process that adds methyl groups to cytosine or adenine nucleotides. Methylated DNA alters gene expression by preventing transcription factor binding [1] or by sometimes favoring the binding of repressors [2, 3]. The regulatory role of DNA methylation is widespread across eukaryotes where it mediates development, environmental responses and disease [4,5,6,7]. In animals, the addition of methyl groups occurs on CpG dinucleotides (cytosine linked to a guanine by a phosphate group) within genes in invertebrates [8] and across genic and intergenic regions in vertebrates [9]. Importantly, changes in DNA methylation levels have been linked to the regulation of phenotypic plasticity [10,11,12]. Temperature-dependent sex determination (TSD) represents a textbook example of phenotypic plasticity (a thermal polyphenism), where individuals with identical genotypes can develop alternative phenotypes (male or female) based on environmental cues [13, 14]. Differential methylation of two genes in the sex-determining pathway has been experimentally identified in a few TSD vertebrates (a fish, a turtle and alligator) [15,16,17] and in other genes during temperature-induced sex reversal in a fish with a mixed sex-determining system (ZZ/ZW GSD and TSD) [18]. However, the extent to which TSD plasticity is mediated by DNA methylation in turtles remains unclear. As an initial step to address this question, it is necessary to characterize the level of methylation in the genome of TSD turtles and to test whether methylation patterns differ in males and females as would be expected if temperature induces sex-specific methylation profiles [18]. Additionally, if TSD is mediated by DNA methylation of the regulatory network of gonadal development, it would be expected that genes in this network would be the target of differential methylation. And if differential methylation of this network has occurred over evolutionary time, these elements should display a signature left by historic methylation at the DNA sequence level [19, 20].
In silico techniques have been used to estimate historic DNA methylation patterns in animals by measuring the normalized CpG content (or nCpG), i.e., the ratio of the CpG dinucleotide abundance observed at particular genomic regions compared to that expected at random based on the frequency of cytosines and guanines present in the genome [CpG(O/E) = CpG observed/expected] [21]. This value of nCpG is used as a proxy for DNA methylation since (a) DNA methylation is almost entirely targeted to CpG dinucleotides in animals [22], and (b) 5-methylcytosine has the tendency to undergo spontaneous deamination which converts it to thymine leaving a footprint that reflects historic methylation levels [19, 20, 23, 24]. Therefore, nCpG is inversely correlated with the extent of DNA methylation such that in hypermethylated regions (where cytosines within methylated CpGs have been converted to thymine) the nCpG is less than one. On the other hand, an nCpG ratio equal to 1 is indicative of no deviation from random expectation, while a value greater than one indicates hypomethylated regions.
The genome-wide patterns of animal DNA methylation vary within and between invertebrates and vertebrates. For instance, in multiple ant species the genomic distribution of nCpG values is unimodal and centered around 1 [25], whereas the nCpG distribution in honeybee and pea aphid genic regions is bimodal with one peak centered at 0.5 and the other at 1 [26]. Among vertebrates, the promoter regions of eutherian mammals, opossum, chicken, lizard, frog and fish also show a bimodal pattern, such that genes with lower nCpG content (LCG promoters) undergo higher methylation and are linked to tissue-specific expression, and those with higher nCpG content (HCG promoters) are hypomethylated and are linked to broad patterns of gene expression [27, 28]. In contrast, the nCpG content of promoters in platypus (and the urochordate sea squirt) is unimodal [27, 29].
Other studies across various vertebrates reported a genome-wide average CpG ratio of ~25% for mammals and birds and of ~35% for fish and amphibians using alternative approaches such as restriction enzyme assays and HPLC [19, 24, 30]. Because of earlier difficulties to study nCpG distribution in reptilian DNA sequences [31], data at this level are scarce for reptiles. Indeed, pioneering work on reptilian methylation levels focused more heavily on global levels [24, 31], and to our knowledge, only one study examined nCpG content and DNA methylation [27] and another study examined non-methylated islands [32] in anole lizards. Thus, the pattern of genome-wide nCpG distribution in turtles and in TSD vertebrates remains unknown at the DNA sequence level. This is critical given that nCpG not only reflects historic DNA methylation, but is a factor that mediates current DNA methylation levels in ways that influence transcription levels.
Here, we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle, Chrysemys picta, to identify methylated regions of interest that may be important for phenotypically plastic sexual development, which would represent candidates for further analyses in future studies. A similar approach was used to investigate the association of DNA methylation and TSD in a ZZ/ZW + TSD fish by examining DNA methylation in mature individuals [18]. We concentrate our search on genes known to regulate vertebrate primary sexual development [33], as well as genes that by their nature may help mediate TSD sex-specific development, such as genes of the epigenetic machinery [34], hormonal pathways [35] or general sensing responses [36]. Second, we examine the genome-wide nCpG distribution in the painted turtle genome [37, 38] and compare it to that of other vertebrates. We then test whether this index predicts the DNA methylation landscape in turtle gonads and also test for a signature of historic methylation on gene regulators of vertebrate sexual development. Thus, our study provides the first insight into the association between nCpG content and differential methylation in any TSD vertebrate.
Genome-wide CpG distribution is bimodal in turtle promoters, introns, exons and intergenic regions
In all genomic regions, the nCpG was much lower than the expected ratio of 1, predicting that a significant fraction of the painted turtle genome is methylated. Notably, for the exons the overall distribution of nCpG was centered at 0.35, whereas the distributions for the rest of the profiled regions were centered at 0.25. Importantly, although at first glance the profiles of genome-wide normalized CpG content (nCpG) appeared unimodal, statistical analyses fitting mixture models to the data [39] revealed higher support for bimodal distributions across gene bodies (exons plus introns), exons and introns individually, and the intergenic regions (p value of likelihood ratio test <0.00001) (Fig. 1a–d). Interestingly, bimodality was even more obvious in the nCpG profiles of promoter regions measured as 3000, 600, 300 or 150 bases upstream of exon 1 (p value of likelihood ratio test <0.00001) (Fig. 1e–h).
Distribution of normalized CpG (nCpG) content in the Chrysemys picta genome by region. a Exons only (CDS), b introns only, c exons and introns, d intergenic regions and e–h at 3000, 600, 300 and 150 bases upstream of exon 1. Fitted Gaussian density curves for the bimodal distribution along with their respective peak values are indicated in red and green
A substantial portion of the turtle genome is methylated
We used MeDIP-seq to characterize the DNA methylation landscape broadly by profiling differentially methylated regions rather than individual cytosines [40] in two pools of gonads from five male hatchlings each and two pools from five female hatchlings each (Table 1). Over 98% of the MeDIP-seq reads from the male and female hatchling gonads mapped to the C. picta genome [37]. The methylome analysis uncovered ~2.95 million methylated 500-bp windows, totaling 1.48 gigabases in size, or ~57% of the assembled genome [37], and overlapping with 17,646 genes. This corresponds to 78% of the CpG nucleotides in the genome. A total of 40% of the methylated windows fall within gene bodies which is significantly less than the 46% located within 50-kb-upstream sequences that have potential regulatory functions (permutation test p value =0.001). The remaining 14% of methylated windows fall in intergenic regions outside of gene bodies or outside sequences 50 kb upstream of all genes.
Table 1 Illumina library statistics for MeDIP sequencing of Chrysemys picta hatchling gonads
Sexually dimorphic DNA methylation varies by gene region
MeDIP-seq results from the biological replicates revealed strong differences between males and females above and beyond the differences between individual replicates (Fig. 2). Positive (methylated DNA) and negative (unmethylated DNA) controls were used during the MeDIP step and ensure that the MeDIP worked properly, such that the observed variation between the two replicates of the same sex likely reflects natural population-level variation in methylation among individuals. This was observed also using PCR of DNA of multiple males and females digested or not with a restriction enzyme sensitive to DNA methylation (Fig. 3). Our MeDIP-seq analysis revealed 5647 differentially methylated windows between the sexes. Of these, 3076 windows were hypermethylated in females (in 2414 genes) and 2571 windows in males (in 2086 genes) (Fig. 2a). The log-fold change in methylation between the sexes was highest in introns, followed by promoter sequences and finally exons (Fig. 2c). Additionally, differential methylation of exons was around half that of promoters (there were 536 differentially methylated windows in promoter sequences vs. 281 in exons). Finally, 541 differentially methylated genes contained multiple windows with contrasting sex-specific methylation, such that some windows were hypermethylated in male hatchlings and other windows within the same gene were hypermethylated in females (Fig. 2d; Additional file 1: Table S1).
a RPKM heatmap of differentially methylated genes in Chrysemys picta (rows) clustered by mean methylation level per gene. Methylation levels were scaled to [−1.5, 1.5] to indicate genes undergoing high (green) and low (red) relative methylation. b Normalized CpG content of all annotated genes (red), experimentally verified to be methylated using MeDIP-seq (yellow) and differentially methylated (purple). c Fold change in methylation (red: hypermethylated in female; green hypermethylated in male) as seen in gene bodies (exons + introns), exons only, introns only and promoters. d Examples of genes possessing multiple windows that displayed sex-specific methylation. e, f Scatterplot of normalized CpG content (nCpG) in methylated windows occurring (e) in gene bodies relative to nCpG of gene bodies and f in promoters relative to nCpG of the complete promoter sequence (~5 kb upstream). Differentially methylated windows in hatchlings are overlaid, with those hypermethylated in males indicated in blue, and those hypermethylated in females indicated in red
Validation of sexually dimorphic DNA methylation of the fezf2 gene in Chrysemys picta hatchling gonads by methylation-sensitive restriction enzyme PCR. DNA from three females (1, 2 and 3) and three males (4, 5 and 6) was digested or not with HpaII (HpaII cuts unmethylated DNA). Expected size of the PCR amplicons from PCR primers F1 plus R1, or F1 and R2 are indicated by the arrows (non-specific PCR products were also obtained). L = DNA ladder (1 kb plus). Amplification of the expected fragment in females and not males using the F1 + R2 primers confirms the hypermethylation of fezf2 in females detected using the MeDIP-seq data
Potential mediators of thermal transduction and regulators of sexual development are differentially methylated
While no particular GO terms were significantly enriched in the MeDIP-seq data after controlling for false discovery (Additional file 1: Tables S2, Additional file 1: Table S3), we detected some interesting pathways and genes with differential methylation. We focused our attention on reptilian homologs of genes that govern mammalian gonadogenesis [33, 41] and on genes involved in painted turtle gonadal development [42] whose biological functions may help convert the temperature signal into the sexual fate of TSD embryos (Additional file 1: Table S4, Additional file 1: Table S5, Additional file 1: Table S6, Additional file 1: Table S7, Additional file 1: Table S8, Additional file 1: Table S9, Additional file 1: Table S10, Additional file 1: Table S11). Full list of differentially methylated genes and GO pathways are presented in Additional file 1: Table S12, Additional file 1: Table S13, Additional file 1: Table S14). Among these were a number of kinases, androgen-/estrogen-related genes, histone- and ubiquitin-related genes, heat shock and transient potential receptor genes that displayed sexually dimorphic methylation. Interestingly, members of the Wnt signaling pathway and genes involved in transcriptional regulation tended to be hypermethylated in males relative to females, whereas genes involved in cell and neuron differentiation tended to be hypermethylated in female hatchlings. Importantly, 13 out of 53 reptilian homologs of genes involved in mammalian gonadogenesis [33, 41] were differentially methylated. For instance, genes that regulate testicular formation, some of which are highly expressed during embryogenesis at male-producing temperature in TSD turtles, including Amh, Ar, Gata4, Lhx1, Lhx9 and Sf1 [42,43,44,45,46], were significantly hypermethylated in female hatchlings. In contrast, genes important in ovarian formation in mammals, such as Wnt4 and Emx2 [47, 48], were hypermethylated in males. Differential methylation also varied by genic region, perhaps reflecting differences in the type of regulation that DNA methylation might exert in several genes in the sexual development network. For instance, while some of these genes exhibited hypermethylation in the promoter regions near the 5′ end, others were hypermethylated in their gene bodies, mostly in their intronic sequences. Among them, Wt1 exhibited three hypermethylated intronic windows in male hatchlings. Yet other genes, such as Lhx1 and Gata4, show female hypermethylation in both promoter (1 window each) and intronic sequences (two and three windows, respectively).
CpG content is a good in silico predictor of DNA methylation
There is no significant difference between the number of genes predicted to be methylated using the nCpG index and the number of genes identified as methylated by MeDIP-seq (two-sample Kolmogorov–Smirnov test, p = 0.1931), suggesting that in silico predictions from the nCpG index are fairly accurate. Specifically, 94% of all methylated gene bodies detected experimentally via MeDIP-seq had an in silico nCpG value of 0.5 or lower, thereby showing a strong association between CpG depletion and actual methylation. Further, 98.7% (16,884 genes) of the 17,104 genes with nCpG ≤ 0.5 are methylated. Only 75 out of 17,646 methylated genes show an nCpG ≥ 0.8, including 36 out of 90 methylated tRNA genes.
Differentially methylated genes suffered greater historic methylation
Interestingly, the nCpG content of the differentially methylated genes was significantly lower (mean = 0.282, range = 0–1.27) than for all methylated genes (mean = 0.306, range = 0–12) identified by MeDIP-seq (resampling test p value =0.001) (Fig. 2b), indicating that genes displaying sexually dimorphic methylation show a signature of higher historic DNA methylation. This pattern was observed in gene bodies (comprising of exons and introns) as well as promoters (Fig. 2e, f), suggesting that differential methylation of hatchling gonads was restricted to the genomic regions of greater CpG depletion.
DNA methylation appears associated with expression of some genes but not all
In the absence of hatchling transcriptomes, we leveraged gonadal transcriptomes of painted turtle embryos incubated at male-producing temperature (MPT) and female-producing temperature (FPT) that we obtained in another study [42], and uncovered potential candidates for future study when we compared the methylome signatures in hatchlings and differential transcription patterns in late-stage embryos (stage 22). Namely, some genes overexpressed in male embryos were hypermethylated in female hatchlings (58 out of 394), and some genes overexpressed in female embryos were hypermethylated in male hatchlings (40 out of 754) (Additional file 1: Table S15). Thus, indirect evidence was detected suggesting that DNA methylation of numerous genes might mediate their sexually dimorphic transcription and that this influence may not be global but a gene-by-gene effect where some genes may be affected while others are not. A direct test also revealed the lack of such global effect, as the Fisher exact test on the contingency table of genes methylated or not that were differentially expressed or not was not significant (p > 0.9999). However, we note that while intriguing, the finding of gene-by-gene effects should be taken with caution given the difference in life stage between the transcriptomic and methylomic data.
Repeat elements abound nearby methylated genes and share identical pattern of differential methylation
Because repetitive DNA sequences such as transposable elements can be subject to silencing by DNA methylation [49] which could affect nearby genes, we analyzed the repeat content of the methylome. RepeatMasker analyses revealed that around 40% of the methylome consists of repeats (which represent 9.25% of the CPI genome—Additional file 1: Table S16), with significant representation from the CR1 and HAT repeat categories (~45% of the methylome repeats). CR1 repeats were also the most abundant in the C. picta gonadal transcriptome [42] (Fig. 4a; Additional file 1: Table S16). Furthermore, the relative abundance among repeat categories as a fraction of the genome (Additional file 1: Table S16) did not differ significantly between the genome, methylome and transcriptome (pairwise Kruskal–Wallis tests; all p values >0.453, Fig. 4a), even though the absolute abundance of each category and of all repeats combined varied (e.g., 25, 9 and 1% of the genome, methylome and transcriptome consisted of repeats, respectively—Additional file 1: Table S16). Additionally, methylated repeats were significantly more concentrated around 95% of all methylated genes (resampling test p value =0.001) and relatively scarce around non-methylated genes (Table 2). Of these, HAT repeats were the most abundant. Methylated repeats were also common (although significantly less so) in the vicinity of differentially methylated genes (~80% instead of 95%; permutation test p value =0.001). Of these, DIRS repeats were the most common. Interestingly, the direction of the sex-specific methylation was identical for 70% of the differentially methylated genes and their neighboring repeats (Table 2).
Chrysemys picta repeat abundance. a Relative abundance of repeat categories in the Chrysemys picta genome [37], and their relative abundance in the hatchling gonadal methylome (this study) and embryonic gonadal transcriptome [42] as a fraction of the genome. Note that repeat abundance in the transcriptome is plotted in gray and scaled by the right-hand axis for visualization purposes. b Relative abundance of various repeat categories within the fraction of repeats present in the C. picta genome [37] versus the methylome (this study). c Regression of transcriptomic repeat abundance as a function of repeat abundance in the genome (p = 0.0001) and d in the methylome (p = 0.0007). Abundance of repeats in the transcriptome is slightly better explained by their genomic abundance (R 2 = 0.59) than by their abundance in the methylome (R 2 = 0.49)
Table 2 Overabundance of methylated repeats upstream of the transcription starting site of differentially and non-differentially methylated genes compared to non-methylated genes in Chrysemys picta hatchlings
Importantly, significantly more genes that were differentially methylated in hatchlings and differentially transcribed in stage 22 embryos were nearby methylated repeats than genes equally expressed in male and females (permutation test p value =0.001). The distribution of categories of methylated repeats did not differ between the differentially and non-differentially expressed genes. Using regression, we evaluated the effect of DNA methylation on repeat silencing, by assessing the repeat transcription level as a function of repeat abundance in the genome and as a function of repeat abundance in the methylome. In both cases, the relationship was highly significant and explained a significant proportion of variation in repeat transcription level (repeat abundance in genome: slope = 0.0239, p = 0.0001, r 2 = 0.59; repeat abundance in methylome: slope = 0.0569, p = 0.0007, r 2 = 0.49), although the variation in repeat expression itself was small (Fig. 4c, d). Further, multiple regression analyses with all variables combined did not improve the explanatory power significantly, implying that repeat transcription, repeat methylation status and repeat genomic abundance are tightly linked.
Genomic approaches are advancing our understanding of phenotypic plasticity at unprecedented rates, including the role that DNA methylation plays in mediating plastic responses to environmental inputs [10, 11, 21, 50]. Here we tested whether regulators of vertebrate sexual development (and of other sensing and epigenetic responses) are subjected to differential methylation in male and female hatchlings of a turtle with phenotypically plastic sex determination (C. picta) lacking sex chromosomes [51]. Based on these data, we identified candidate genes that may mediate TSD epigenetically. We also characterized the genome-wide nCpG distribution in the painted turtle genome, the first such analysis in reptiles and TSD vertebrates, and found that this proxy predicts reasonably well the methylation levels estimated using MeDIP-seq. Further, nCpG profiles helped assess historic methylation levels of genes regulating vertebrate sexual development. Below we highlight our most important observations and propose working hypotheses to guide further research. Our sex-specific methylomes represent an important genomic resource to aid investigations of the epigenetic regulation of environmental responses.
Although bisulfite sequencing is the gold standard for methylation studies requiring single-nucleotide resolution, MeDIP-seq is appropriate for studies seeking to profile broader patterns of DNA methylation rather than individual cytosines [40], as was our goal. MeDIP-seq provides methylation profiles at a resolution of 150–200 bp, signals tend to concentrate on CpG-rich genomic regions with high methylation levels, and the methylome information can be comprehensive [52], less sequence-biased and fairly concordant with bisulfite-sequencing data for genomes similar in size to those of turtles [37, 52, 53]. The reduction in genomic complexity afforded by MeDIP-seq is also advantageous for taxa with large genomes as is the case of the painted turtle [37, 38]. It should be noted that methylation levels are correlated with 1 kb of neighboring CpG sites [54]. Nonetheless, future bisulfite sequencing will permit evaluating methylation patterns at higher resolution, particularly in non-CpG genomic regions, which was precluded in our study. But importantly, given that methylation status may be more stable at the level of DNA domains rather than at single nucleotides [55], MeDIP-seq snapshots do contain highly valuable information as discussed below.
Bimodal distribution of nCpG turtle genomes matches most vertebrate promoters but not introns
Our data revealed that nCpG values follow a bimodal distribution at the promoter regions of genes in painted turtles (Fig. 1), consistent with most major vertebrates lineages studied ([56], Table 3), but not with platypus [29] or tunicates [27]. However, some interesting differences exist among bimodal patterns between turtles and other vertebrates. For instance, the CpG bimodality is less pronounced in turtle promoters than in human and chicken (Fig. 1 and [27]). Namely, the low- and high-CpG modes of the turtle promoter distribution overlap more extensively and are centered around lower CpG values than those of human and chicken promoters (Fig. 1 and [27]), suggesting a potentially larger role for DNA methylation as transcriptional regulator in reptiles than previously anticipated (i.e., greater historic methylation in turtle than in human and chicken). In human, high-CpG promoters are more abundant than low-CpG promoters, whereas the opposite is true in turtles (Fig. 1 and [27]). Human high-CpG promoters are hypermethylated less frequently whereas low-CpG promoters tend to be hypermethylated more often [28, 29], whereas differentially methylated genes in turtles had lower CpG content (Fig. 2b). The full implications of the observed methylation patterns in turtles are unknown. For instance, human promoters are methylated in somatic as well as germline cells and thus could be heritable [55]. Further, the CpG observed/expected ratio [CpG(O/E)] affects gene expression levels and breath in humans [57]. But further research is needed to test if the same is true in turtles. Unlike promoters, the bimodal nCpG content of turtle introns (Fig. 1) differs from other vertebrates, as it is unimodal in fish, amphibian, lizard, bird and human [27], and more strikingly bimodal in tunicates [27]. Subtle bimodality can pass undetected during qualitative evaluations (e.g., platypus [29]), unless it is explicitly tested statistically (this study and [27]). Notably, recent evidence suggests bimodal gene body methylation exists in mammals, birds and fish [58].
Table 3 Summary of exemplar studies exploring the diversity of nCpG distributions in vertebrate genomes.
The overall turtle CpG depletion agrees generally with vertebrates (where it is correlated with high DNA methylation) compared to invertebrates [19, 59, 60]. The overall distribution of nCpG values in the C. picta genome, disregarding bimodality, is higher for exons than for other regions (introns, promoters and intergenic sequences), suggesting lower exon methylation. In contrast, human exons are more highly methylated than introns [61]. CpG depletion in exons could affect transcription as in humans where it downregulates genes [62, 63]. Thus, we hypothesize that the higher nCpG content of C. picta's exons relative to other regions (i.e., lower exon CpG depletion relative to overall depletion) might be the result of natural selection to preserve gene expression given that nCpG content is much lower in turtles than in other vertebrates, and perhaps also to prevent the accumulation of CpG to TpG mutations [19] that could produce non-functional proteins [64].
nCpG content is a reasonable predictor of DNA methylation, except for tRNAs, and matches other vertebrates
Our results link DNA methylation at the resolution measured by MeDIP-seq and CpG depletion, consistent with humans [28] and insects [26]. Indeed, 94% of all the methylated genes revealed by MeDIP-seq had an nCpG ≤ 0.5, and 98.7% of all genes with nCpG ≤ 0.5 were actually methylated. Thus, CpG content is a reasonable indicator of methylation status, except for tRNAs. Namely, 90 of 182 annotated tRNAs were methylated, 40% of which (36/90) have an nCpG ≥ 0.8, indicating that CpG depletion has been suppressed historically in many tRNAs, perhaps to prevent CpG to TpG mutations [19]. This agrees with the high conservation of tRNA sequences in all domains of life [65]. Comparable methylation levels were found here (57% of the genome, including 78% of all CpG dinucleotides) as in fish, salamander, snake, birds and mammals using various approaches [24, 27, 31, 59, 66].
Turtle methylation patterns suggest cis-regulation via DNA methylation
Sequences up to 50 kb upstream of genes in painted turtle, which should have potential regulatory roles, were methylated at significantly higher levels (46%) than gene bodies (40%) in hatchling gonads. Of these upstream windows, 16% fall within the boundaries of another gene upstream of the focal gene (within an intron in 94% of those cases), suggesting that perhaps some methylated genes may serve as alternative upstream promoters for downstream genes, as in humans [67, 68]. Importantly, DNA methylation in promoters is linked to gene silencing, whereas methylation in gene bodies may modulate transcription [8, 69], but further research is needed to test these hypotheses directly in turtles.
DNA methylation varies by sex, gene and gene region
Results from our biological duplicates, which combined five males or five females each, were fairly consistent (Fig. 2), strengthening the inferences of sex-specific methylation. However, differences were also detected within sex, likely revealing genetic variation among individuals associated with DNA methylation (Fig. 3), but without obscuring some evident sex-specific methylation patterns.
Interestingly, 3% (541) of all methylated genes contained multiple differentially methylated windows, some hypermethylated in females while others in the same gene were hypermethylated in males (Fig. 2d). Introns showed the highest log-fold change in methylation between the sexes, followed by promoters and finally exons (Fig. 2c), but it is unclear whether turtle intronic methylation modulates transcription rather than silencing genes [8, 69], or whether it is linked to the transcription of antisense RNA, as observed in genes important for urogenital development in other vertebrates [70].
Because the males and females studied here were obtained at contrasting temperatures that produce a single sex (only males at 26 °C and only females at 31 °C), perhaps the sex-specific methylation in hatchlings is also temperature specific. Thermosensitive methylation of distinct windows within the same gene could potentially lead to alternative splicing of sex-specific transcripts, as in humans [71]. Alternatively, the observed differential methylation could be a consequence of sexual differentiation. Further studies are needed to test these alternative hypotheses.
Known regulators of sexual development are differentially methylated in TSD turtles
We investigated 53 genes in the mammalian urogenital regulatory network, whose reptilian homologs are differentially expressed in TSD turtles, including Wt1 [72, 73], Sf1 [74, 75], Dax1 [41, 76], Sox9 [76,77,78,79], Aromatase [17, 80], Dmrt1 [76, 78, 81], Esr1 [82, 83] and Rspo1 [79]. Of these, methylation has been studied only in aromatase and Sox9 and demonstrated to influence transcription in developing TSD reptiles (slider turtle [17], American alligator [84]) and a TSD fish [15]. Our results provide evidence (1) that many more genes in this regulatory network are differentially methylated between male and female gonads, (2) that sexually dimorphic methylation persists post-hatching and (3) that it affects some genes that are differentially expressed during gonadal development.
Genes such as Amh, Ar, Gata4, Lhx9 and Sf1 that govern testicular formation are upregulated at male-producing temperature (MPT) late in the thermosensitive period of C. picta [42, 43] and displayed hypermethylated promoters in female hatchlings (Table 4). The exception was Sf1, whose introns were differentially methylated (Table 4). Additionally, three hypermethylated intronic windows in male hatchlings were observed in Wt1, a gene important for the formation of the bipotential gonad and later testicular development [85] that is upregulated at MPT early in the thermosensitive period but not by stage 22, and is a candidate TSD master gene [43, 72, 85]. The consequences of turtle intronic methylation remain unknown for alternative splicing, the use of alternative promoters, transcription of antisense RNA or transcription modulation [8, 64, 69, 70]. Data for Amh, Ar, Gata4 and Lhx9 suggest that upregulation in one sex may result from silencing by promoter methylation in the opposite sex, which was also true for Wnt4 and Emx2 that govern ovarian formation [47, 48] and were hypermethylated in male hatchlings. These types of regulation do not appear to be ubiquitous as some other genes in the gonadal development network were differentially expressed but showed no differential methylation. Interestingly, Lhx1 and Gata4 displayed female hypermethylation in promoters (1 window each) and introns (2 and 3 windows, respectively) suggesting a potentially more complex regulation by methylation of these elements (Table 4).
Table 4 Summary of differentially methylated (FDR cutoff: 0.05) genes in hatchlings that are putatively involved in reptilian gonadogenesis
Our findings suggest that methylation marks may be stable and persist post-hatching, perhaps longer-term as in mammals [86, 87]. Yet notably, aromatase showed no differential methylation in hatchlings, counter to observations in slider turtle embryos [17]. Assuming that embryonic methylation in slider and painted turtles is similar as is aromatase transcription [43, 75], our finding in C. picta hatchlings suggests that embryonic aromatase methylation may be transient, or alternatively, that it passed undetected by unknown technical limitations. Further research will help elucidate between these hypotheses directly.
A contrast of our turtle methylomes with those from the ZZ/ZW + TSD tongue sole fish [18] revealed the overlap of 22 genes (Additional file 1: Table S17) out of 56 genes differentially transcribed or methylated in tongue sole gonads (Additional file 1: Table S8 in [18]). Fish and painted turtle differed in overall patterns of methylation and expression, non-surprisingly given that fish transcription and that in other vertebrates has diverged profoundly [43]. The only exception was Lhx9 which was upregulated at MPT in turtle embryos and was hypermethyled in both the tongue sole and painted turtle. Lastly, two other differentially methylated genes in turtle are worth mentioning, the epidermal growth factor receptor (Egfr) which has been previously implicated in sexual dimorphism in Drosophila [88], and Mafb, a gene with sexually dimorphic expression responsible for masculinization of male genitalia in mice [89].
Potential thermal transmitters are differentially methylated in TSD turtles
Sexually dimorphic methylation was present in hatchling gonads in several gene candidates that may help convert the environmental temperature into sexual development signals to establish the sexual fate in TSD taxa [13, 36, 43, 90] (Additional file 1: Table S4, Additional file 1: Table S5, Additional file 1: Table S6, Additional file 1: Table S7, Additional file 1: Table S8, Additional file 1: Table S9, Additional file 1: Table S10, Additional file 1: Table S11), rendering them targets for future study. These include several heat shock genes, transient receptor potential genes, a number of kinases, androgen-/estrogen-related genes and histone-related genes. Namely, several heat shock protein genes (e.g., Hspa4, Hspa12a and Hspa12b) were hypermethylated in males, and they are differentially expressed by temperature in alligators and could play a role in TSD [36]. Further, the cold-inducible RNA-binding protein (Cirbp), a putative key TSD gene that is upregulated at FPT in snapping turtles Chelydra serpentina (TSD) [91] and C. picta [42], showed hypermethylated exons in females (Additional file 1: Table S4), perhaps causing differential splicing as in humans [71]. Additionally, transient receptor potential genes (e.g., Trpm1, Trpm2, Trpm3, Trpm7 and Trpm8) can respond to temperature stimuli [92] and were differentially methylated. Moreover, many kinases such as members of the Mapk signaling family which helps activate Sry in mice [93, 94] were differentially methylated, along with many androgen and estrogen signaling genes. Also, genes involved in histone modification directly regulate transcription and can act in a sex-specific manner. For instance, aromatase transcription in embryonic gonads of T. scripta turtle embryos increases by demethylation [17] which depends on local histone acetylation in mammals [95]. We found hypermethylation of histone acetyltransferases (Kat2a and Kat6a) in males and of deacetylases (Hdac4, Hdac7 and Hdac8) in females, but it is unclear whether differential methylation of histone modifiers is linked to sexually dimorphic transcription in turtles.
Differentially methylated genes have undergone greater historic methylation
Intriguingly, differential methylation in C. picta hatchlings coincided with regions that showed significantly lower normalized CpG content relative to the genome-wide methylation levels, a pattern observed in both gene bodies and promoters (Fig. 2e, f) that indicates stronger historic methylation in these genomic regions. This pattern includes various genes in the sex determination/differentiation network and supports the notion that differential methylation of elements in this network has occurred over long periods during turtle evolution, leaving a footprint of historic methylation at the DNA sequence level [19, 20].
Differential repeat methylation is linked to methylation of adjacent genes
Importantly, over 95% of methylated genes (but not unmethylated genes) were nearby methylated repeats. Further, repeat methylation was associated with sexually dimorphic methylation of neighboring genes, as >70% of repeat windows nearby differentially methylated genes were methylated in the same direction as the genes (Table 2). This suggests that repeat methylation could promote DNA methylation in adjacent genes perhaps mediating their transcription in C. picta as occurs in Drosophila and human [96, 97]. Repeat methylation may vary by repeat category and developmental stage as observed in fish [98], but further research is needed to test for DNA methylation in developing turtles.
However, we found that DNA methylation targets repeat elements in general more often than expected by their overall genome abundance, perhaps as a result of repeat silencing. Namely, 10% of the C. picta genome assembly is composed of transposable elements [37] [and 25% of the genome consists overall repeats (Additional file 1: Table S16)], whereas repeats composed 40% of all methylated regions. This result is consistent with findings in rat using our same approach where the high proportion of repeats in the methylome was attributed to repeat silencing (41% repeats in the rat genome and 53% in the rat methylome) [99]. Notably, insights about the relative abundance of repeat categories vary depending on whether the entire turtle genome was considered or only the genomic fraction that was composed of repeats (Fig. 4a, b). On the one hand, relative abundance among repeat categories in the methylome followed their relative abundance in the entire genome (Fig. 4a). For instance, CR1 repeats were the most abundant in the methylome, followed by HATs, DIRS, Gypsy and Harbinger repeats, just as in the painted turtle genome [37]), and in the transcriptome [42] (Fig. 4a). On the other hand, some repeat categories were overrepresented and others underrepresented in the methylome compared to their abundance in just the fraction of the genome composed of repeats (Fig. 4b). Thus, not fully conclusive evidence was found for overall repeat silencing by DNA methylation [64] in turtles as occurs in humans [100]. Interestingly, CR1 repeats concentrate at centromeres of a few C. picta chromosomes [38] such that their high methylome abundance may reflect chromosome stabilization by DNA methylation of centromeric repeats [64].
Ours is the first genome-wide assessment of DNA methylation in reptiles and the first study of sexually dimorphic methylation levels in a purely TSD vertebrate [51]. As such, this study sheds light on the epigenetic modifications that may play a role in the development of phenotypically plastic vertebrates, complementing recent work on a fish with mixed sex determination [18]. Our MeDIP-seq data provide empirical validation of in silico predictions obtained from nCpG content never done before in reptiles and show that nCpG content is a reasonable predictor of actual methylation status at the resolution of MeDIP-seq. We found that painted turtles possess a unique pattern with nCpG values below those of other vertebrates, indicative of global and extensive historic methylation. In contrast, actual methylation levels given the turtle genome CpG content agree with those described for most vertebrates. Our data helped us identify several genes whose methylation status renders them candidates for a putative function in regulating transcription in a thermosensitive manner, a pattern that appears associated with methylation status of neighboring repeat elements. Based on this correlational pattern, we speculate that methylation of some of these elements may play a key role in mediating the sexual fate of TSD reptiles. Further research is warranted to test this hypothesis and the prevalence of DNA methylation in governing the sexual outcome in TSD turtles as it does in sex reversal of fish with mixed sex determination [18], or whether instead, most differential methylation observed in turtle hatchlings is the consequence of sexual differentiation and not its cause.
Eggs from several freshly laid clutches were obtained from a turtle farm, and equal numbers from each clutch were assigned randomly to incubators set at 26 and 31 °C that produce 100% males and 100% females, respectively, under standard incubation conditions [101]. Turtle embryos are at gastrula stage at oviposition [102]. All hatchlings were raised for 3 months in tanks with water heated to ~26 °C and fed ad libitum, to allow full gonadal differentiation prior to sexing. Individual sex was diagnosed by gonadal inspection [51, 103]. All procedures were approved by the IACUC of Iowa State University.
DNA isolation and sequencing
We extracted DNA from the gonads of 3-month-old C. picta hatchlings (ten males and ten females, and pooled in groups of five hatchlings for replication) using the Gentra Puregene DNA extraction kit (Gentra) following the manufacturer's instructions. DNA was processed by BGI Americas following the MeDIP protocol [104] where DNA was fragmented and denatured, and the methylated DNA was immunoprecipitated using the MagMEDI kit (Diagenode). A positive control (methylated DNA) and a negative control (unmethylated DNA) were used at the MeDIP step to ensure that the MeDIP worked properly. The methylated DNA was size selected (100–300 bp) and sequenced using the Illumina HiSeq paired-end protocol. We obtained between 126 million and 163 million 50-bp reads per library amounting to a total of ~564 million reads (Table 1).
Methylome construction and differential methylation analysis
We mapped the sequencing reads to the C. picta genome version 3.0.1 [37] using Bowtie2 version 2.2.5 [105]. We filtered out the unmapped reads using Samtools [106]. We used the MEDIPS package [107] in four steps: (a) to build an index for the C. picta genome and to ensure fast querying of the alignment files. (b) To model read counts under a negative binomial distribution. (c) To quantify mapped read counts per 500-bp windows in RPKM (reads per kilobase of genes per million mapped reads). Mean methylation levels were calculated for each 500-bp window by averaging the corresponding RPKM levels in the replicates. A gene was considered as methylated if it contained one or more methylated windows within its start and end coordinates. (d) To merge methylated windows and compute differential methylation by sex, while controlling for false discoveries [108] at a q-value cutoff =0.05. Any gene with a q-value <0.05 was considered as differentially methylated. Besides the positive and negative controls used during the MeDIP procedure described above, primers were also designed to validate by PCR the differential methylation identified by MeDIP-seq for the gene Fezf2 that was hypermethylated in females compared to males [109]. The primer cocktail consisted of one forward primer (F1: 5′-GGG GTG AAA AAC CAC AG-3′) and two reverse primers, one closer to the F1 (R1: 5′-CAC ACA CAA GGA GG-3′) and one further away (R2: 5′-CAG CAA CAA CTT GAT TTG G-3′) (Fig. 3). Primers F1 + R1 flank a shorter non-methylated area that should amplify in any individual, which is nested within the larger region encompassed by F1 + R2 that contains the differentially methylated area, and which amplifies only when methylation protects the DNA from digestion by a methylation-sensitive restriction enzyme (HpaII in our case). Gonadal DNA from three male and three female hatchlings was digested with HpaII for 30 min at 37 °C. PCR was carried out in 15ul reactions, using 1ul of 100 ng/ul undigested DNA (control) or digested DNA (experimental) as template, and 5ul PCR products along with 1ul 1 kb plus DNA ladder (Invitrogen) were visualized in 1% agarose gels stained with ethidium bromide. The PCR cocktail contained 1.5 mM MgCl2, 0.2 mM dNTPs, 0.4uM of each primer (F1, R1, R2), 0.4U Taq polymerase, ~100 ng DNA, 1.5ul 10X buffer and 10.5 μl water. PCR conditions included initial denaturing at 94 °C for 3 min, followed by 35 cycles of denaturing at 94 °C for 30 s, annealing at 58 °C for 30 s, and extension at 72 °C for 1 min.
We used Kolmogorov–Smirnov tests to determine whether the nCpG content distribution of all annotated genes was comparable to those identified by the MEDIPS package. To test whether the nCpG of the differentially methylated genes varied significantly from that of all the methylated genes, we performed a resampling test [110] by iteratively drawing a random subset of genes (equal to the number of differentially methylated genes) from the entire set of methylated genes. We used RepeatMasker v3.3.0 [111] to identify repeats in the C. picta genome. We used pairwise Kruskal–Wallis tests [112] to test whether repeat abundances from the genome, methylome and transcriptome came from the same population. We used Bedtools v2.17.0 [113] to compute repeats overlapping with methylated regions identified by MEDIPS. We used permutation tests to evaluate differences in the abundance of methylated DNA repeats occurring in vicinity of methylated genes versus non-methylated genes. In the absence of transcriptional data from hatchlings, we used our transcriptomic dataset from late-developing embryos (stage 22) from another study [42] to test for an association between methylation patterns in hatchlings and transcription patterns using Fisher's exact test. For this purpose, a contingency table was generated of genes that were differentially expressed (or not) and differentially methylated (or not), and a hypergeometric distribution was used to then determine the p value. A p value of ~1 would indicate no significant association between differential methylation and differential expression. We used R [114] to perform regression tests to model transcriptome abundance using genomic abundance of repeats and methylated repeats. We used the DAVID Bioinformatics knowledgebase [115] to assess the enrichment of functional categories to which the differentially methylated genes belong.
Analysis of normalized CpG content
The normalized CpG content (nCpG) is calculated as:
$${\text{nCpG}} = \frac{{\left( {\frac{cg}{l}} \right)}}{{\left( {\frac{c}{l}} \right) \times \left( {\frac{g}{l}} \right)}}$$
for a sequence of length l, where c is the number of occurrences of cytosine, g is the number of occurrences of guanine and cg is the number of times cytosine is bordered by guanine linked by a phosphate group (CpG) [21]. In theory, nCpG values can range from 0 to +infinity for an infinitely long sequence, with a value of 1 when the number of CpG dinucleotides observed is equal to the expected based on the sequence length and abundance of C and G. Values <1 denote CpG depletion from what is expected by chance, and values >1 represent overabundance of CpGs from random expectation. The methylation of CpG dinucleotides initiates the deamination of the cytosine, transforming it to thymine, thus lowering the nCpG to <1. Studies conducted in vertebrate and invertebrate animals reveal an upper limit between 2 and 2.5 for various genomic regions [21, 25,26,27, 29, 116].
We used Bedtools [113] to parse the exon, intron, promoter and intergenic coordinates from the C. picta genome [37], given the annotations in gff3 format. We used in-house perl scripts to compute the CpG contents by genomic region. For the nCpG distributions, the R package Mclust [39] was used to assess likelihood of a mixture model with one (G = 1) and two (G = 2) components, with the better fit model decided by a likelihood ratio test.
Amh :
anti-mullerian hormone
Ar :
Cirbp :
cold-inducible RNA-binding protein
CR1:
chicken repeat 1
Dax1 :
dosage-sensitive sex reversal, adrenal hypoplasia critical region, on chromosome X, gene 1
Dmrt1 :
doublesex and mab-3-related transcription factor 1
Egfr :
Emx2 :
empty spiracles homeobox 2
Esr1 :
estrogen receptor 1
Gata4 :
gata binding protein 4
GSD:
genotypic sex determination
HAT:
transposon superfamily named after Hobo from Drosophila melanogaster, Ac from maize and Tam3 from the snapdragon
Hdac x :
histone deacetylase number x
HiSeq:
high-throughput sequencing technology by Illumina
high-performance liquid chromatography
Hsp x :
heat-shock protein number x
Lhx1 :
lim homoeobox 1
Mafb :
MAF BZIP transcription factor B
Mapk :
MeDIP-seq:
methyl DNA immunoprecipitation sequencing
nCpG:
normalized CpG content
RPKM:
reads per kilobase of genes per million mapped reads
Rspo1 :
R-spondin1
Sf1 :
steroidogenic factor 1
Sox9 :
SRY (sex-determining region Y)-box 9
Sry :
sex-determining region Y
tRNA:
Trpm x :
transient receptor potential cation channel subfamily M member number x
TSD:
temperature-dependent sex determination
Wnt4 :
wingless-type MMTV integration site family, member 4
Wt1 :
Wilms tumor protein 1
Watt F, Molloy PL. Cytosine methylation prevents binding to DNA of a hela-cell transcription factor required for optimal expression of the adenovirus major late promoter. Genes Dev. 1988;2(9):1136–43.
Boyes J, Bird A. DNA methylation inhibits transcription indirectly via a methyl-CpG binding-protein. Cell. 1991;64(6):1123–34.
Hendrich B, Bird A. Identification and characterization of a family of mammalian methyl-CpG binding proteins. Mol Cell Biol. 1998;18(11):6538–47.
Varriale A. DNA methylation, epigenetics, and evolution in vertebrates: facts and challenges. Int J Evolut Biol. 2014;2014:475981.
Bergman Y, Cedar H. DNA methylation dynamics in health and disease. Nat Struct Mol Biol. 2013;20(3):274–81.
Cantone I, Fisher AG. Epigenetic programming and reprogramming during development. Nat Struct Mol Biol. 2013;20(3):282–9.
Head JA. Patterns of DNA methylation in animals: an ecotoxicological perspective. Integr Comp Biol. 2014;54(1):77–86.
Suzuki MM, Bird A. DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008;9(6):465–76.
Zemach A, et al. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010;328(5980):916–9.
Szyf M, et al. Maternal care, the epigenome and phenotypic differences in behavior. Reprod Toxicol. 2007;24(1):9–19.
Kucharski R, et al. Nutritional control of reproductive status in honeybees via DNA methylation. Science. 2008;319(5871):1827–30.
Weiner SA, et al. A survey of DNA methylation across social insect species, life stages, and castes reveals abundant and caste-associated methylation in a primitively social wasp. Naturwissenschaften. 2013;100(8):795–9.
Valenzuela N, Lance VA, editors. Temperature dependent sex determination in vertebrates. Washington, DC: Smithsonian Books; 2004.
Tree of Sex Consortium, et al. Tree of sex: a database of sexual systems. Sci Data. 2014;1:140015.
Navarro-Martin L, et al. DNA methylation of the gonadal aromatase (cyp19a) promoter is involved in temperature-dependent sex ratio shifts in the European sea bass. PLoS Genet. 2011;7(12):e1002447.
Parrott BB, et al. Gonadal DNA methylation patterning is affected by incubation temperature in the American alligator, a species undergoing temperature-dependent sex determination. Integr Comp Biol. 2014;54:E327.
Matsumoto Y, et al. Epigenetic control of gonadal aromatase cyp19a1 in temperature-dependent sex determination of red-eared slider turtles. PLoS ONE. 2013;8(6):e63599.
Shao C, et al. Epigenetic modification and inheritance in sexual reversal of fish. Genome Res. 2014;24(4):604–15.
Bird AP. DNA methylation and the frequency of CpG in animal DNA. Nucleic Acids Res. 1980;8(7):1499–504.
Coulondre C, et al. Molecular-basis of base substitution hotspots in Escherichia coli. Nature. 1978;274(5673):775–80.
Elango N, et al. DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera. Proc Natl Acad Sci USA. 2009;106(27):11206–11.
Jabbari K, Bernardi G. Cytosine methylation and CpG, TpG (CpA) and TpA frequencies. Gene. 2004;333:143–9.
Shen JC, et al. The rate of hydrolytic deamination of 5-methylcytosine in double-stranded DNA. Nucleic Acids Res. 1994;22(6):972–6.
Jabbari K, et al. Evolutionary changes in CpG and methylation levels in the genome of vertebrates. Gene. 1997;205(1–2):109–18.
Simola DF, et al. Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality. Genome Res. 2013;23(8):1235–47.
Glastad KM, et al. DNA methylation in insects: on the brink of the epigenomic era. Insect Mol Biol. 2011;20(5):553–65.
Elango N, Yi SV. DNA methylation and structural and functional bimodality of vertebrate promoters. Mol Biol Evol. 2008;25(8):1602–8.
Weber M, et al. Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007;39(4):457–66.
Yang H, et al. Relating gene expression evolution with CpG content changes. BMC Genomics. 2014;15:693.
Aissani B, Bernardi G. CpG islands—features and distribution in the genomes of vertebrates. Gene. 1991;106(2):173–83.
Varriale A, Bernardi G. DNA methylation in reptiles. Gene. 2006;385:122–7.
Long HK, et al. Epigenetic conservation at gene regulatory elements revealed by non-methylated DNA profiling in seven vertebrates. Elife. 2013;2:e00348.
Eggers S, et al. Genetic regulation of mammalian gonad development. Nat Rev Endocrinol. 2014;10(11):673–83.
Kuroki S, et al. Epigenetic regulation of mouse sex determination by the histone demethylase Jmjd1a. Science. 2013;341(6150):1106–9.
Carmi I, et al. The nuclear hormone receptor Sex-1 is an X-chromosome signal that determines nematode sex. Nature. 1998;396:168–73.
Kohno S, et al. Potential contributions of heat shock proteins to temperature-dependent sex determination in the American alligator. Sex Dev. 2010;4(1–2):73–87.
Shaffer HB, et al. The western painted turtle genome, a model for the evolution of extreme physiological adaptations in a slowly evolving lineage. Genome Biol. 2013;14(3):21–2.
Badenhorst D, et al. Physical mapping and refinement of the painted turtle genome (Chrysemys picta) inform amniote genome evolution and challenge turtle-bird chromosomal conservation. Genome Biol Evol. 2015;7(7):2038–50.
Fraley C, Raftery A. MCLUST version 3 for R: Normal mixture modeling and model-based clustering. Technical Report No. 504 2006, Department of Statistics, University of Washington.
Beck S, Rakyan VK. The methylome: approaches for global DNA methylation profiling. Trends Genet. 2008;24(5):231–7.
Valenzuela N. Evolution of the gene network underlying gonadogenesis in turtles with temperature-dependent and genotypic sex determination. Integr Comp Biol. 2008;48(4):476–85.
Radhakrishnan S, et al. Transcriptomic responses to environmental temperature by turtles with temperature-dependent and genotypic sex determination assessed by RNAseq inform the genetic architecture of embryonic gonadal development. PLoS ONE. 2017;12:e0172044.
Valenzuela N, et al. Transcriptional evolution underlying vertebrate sexual development. Dev Dyn. 2013;242(4):307–19.
Yeh SY, et al. Generation and characterization of androgen receptor knockout (ARKO) mice: an in vivo model for the study of androgen functions in selective tissues. Proc Natl Acad Sci USA. 2002;99(21):13498–503.
Oreal E, et al. Different patterns of anti-Mullerian hormone expression, as related to Dmrt1, Sf-1, Wt1, Gata-4, Wnt-4, and Lhx9 expression, in the chick differentiating gonads. Dev Dyn. 2002;225(3):221–32.
Manuylov NL, et al. Conditional ablation of Gata4 and Fog2 genes in mice reveals their distinct roles in mammalian sexual differentiation. Dev Biol. 2011;353(2):229–41.
Tripathi V, Raman R. Identification of Wnt4 as the ovary pathway gene and temporal disparity of its expression vis-a-vis testis genes in the garden lizard, Calotes versicolor. Gene (Amst). 2010;449(1–2):77–84.
Pellegrini M, et al. Emx2 developmental expression in the primordia of the reproductive and excretory systems. Anat Embryol. 1997;196(6):427–33.
Weisenberger DJ, et al. Analysis of repetitive element DNA methylation by MethyLight. Nucleic Acids Res. 2005;33(21):6823–36.
Bonasio R. The role of chromatin and epigenetics in the polyphenisms of ant castes. Brief Funct Genomics. 2014;13(3):235–45.
Valenzuela N, et al. Molecular cytogenetic search for cryptic sex chromosomes in painted turtles Chrysemys picta. Cytogenet Genome Res. 2014;144:39–46.
Clark C, et al. A comparison of the whole genome approach of MeDIP-seq to the targeted approach of the infinium humanmethylation450 Beadchip for methylome profiling. PLoS ONE. 2012;7(11):e50233.
Su Y, et al. Genome-wide DNA methylation profile of developing deciduous tooth germ in miniature pigs. BMC Genom. 2016;17(1):1–9.
Eckhardt F, et al. DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006;38(12):1378–85.
Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6–21.
Jiang N, et al. Conserved and divergent patterns of DNA methylation in higher vertebrates. Genome Biol Evol. 2014;6(11):2998–3014.
Park J, et al. What are the determinants of gene expression levels and breadths in the human genome? Hum Mol Genet. 2012;21(1):46–56.
Keller TE, et al. Evolutionary transition of promoter and gene body DNA methylation across invertebrate–vertebrate boundary. Mol Biol Evol. 2016;33(4):1019–28.
Bird AP, Taggart MH. Variable patterns of total DNA and rDNA methylation in animals. Nucleic Acids Res. 1980;8(7):1485–97.
Yi SV, Goodisman MAD. Computational approaches for understanding the evolution of DNA methylation in animals. Epigenetics. 2009;4(8):551–6.
Laurent L, et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010;20(3):320–31.
Bauer AP, et al. The impact of intragenic CpG content on gene expression. Nucleic Acids Res. 2010;38(12):3891–908.
Krinner S, et al. CpG domains downstream of TSSs promote high levels of gene expression. Nucleic Acids Res. 2014;42(6):3551–64.
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.
Widmann J, et al. Stable tRNA-based phylogenies using only 76 nucleotides. RNA. 2010;16(8):1469–77.
Tweedie S, et al. Methylation of genomes and genes at the invertebrate-vertebrate boundary. Mol Cell Biol. 1997;17(3):1469–75.
Kimura K, et al. Diversification of transcriptional modulation: large-scale identification and characterization of putative alternative promoters of human genes. Genome Res. 2006;16(1):55–65.
Maunakea AK, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466(7303):U131–253.
Huh I, et al. DNA methylation and transcriptional noise. Epigenetics Chromatin. 2013;6:9.
Dallosso AR, et al. Alternately spliced WT1 antisense transcripts interact with WT1 sense RNA and show epigenetic and splicing defects in cancer. RNA. 2007;13(12):2287–99.
Maunakea AK, et al. Intragenic DNA methylation modulates alternative splicing by recruiting MeCP2 to promote exon recognition. Cell Res. 2013;23(11):1256–69.
Valenzuela N. Relic thermosensitive gene expression in a turtle with genotypic sex determination. Evolution. 2008;62(1):234–40.
Spotila LD, et al. Sequence and expression analysis of WT1 and Sox9 in the red-eared slider turtle, Trachemys scripta. J Exp Zool. 1998;284:417–27.
Valenzuela N, et al. Comparative gene expression of steroidogenic factor 1 in Chrysemys picta and Apalone mutica turtles with temperature-dependent and genotypic sex determination. Evol Dev. 2006;8(5):424–32.
Ramsey M, et al. 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–91.
Torres Maldonado LC, et al. Expression profiles of Dax1, Dmrt1, and Sox9 during temperature sex determination in gonads of the sea turtle Lepidochelys olivacea. Gen Comp Endocrinol. 2003;129:20–6.
Barske LA, Capel B. Estrogen represses SOX9 during sex determination in the red-eared slider turtle Trachemys scripta. Dev Biol. 2010;341(1):305–14.
Valenzuela N. Multivariate expression analysis of the gene network underlying sexual development in turtle embryos with temperature-dependent and genotypic sex determination. Sex Dev. 2010;4(1–2):39–49.
Matsumoto Y, et al. Changes in gonadal gene network by exogenous ligands in temperature-dependent sex determination. J Mol Endocrinol. 2013;50(3):389–400.
Valenzuela N, Shikano T. Embryological ontogeny of Aromatase gene expression in Chrysemys picta and Apalone mutica turtles: comparative patterns within and across temperature-dependent and genotypic sex-determining mechanisms. Dev Genes Evol. 2007;217:55–62.
Kettlewell JR, et al. Temperature-dependent expression of turtle Dmrt1 prior to sexual differentiation. Genesis. 2000;26(3):174–8.
Bergeron JM, et al. Cloning and in situ hybridization analysis of estrogen receptor in the developing gonad of the red-eared slider turtle, a species with temperature-dependent sex determination. Dev Growth Differ. 1998;40(2):243–54.
Chávez B, et al. Cloning and expression of the estrogen receptor-alpha (Esr1) from the Harderian gland of the sea turtle (Lepidochelys olivacea). Gen Comp Endocrinol. 2009;162:203–9.
Parrott BB, et al. Influence of tissue, age, and environmental quality on DNA methylation in Alligator mississippiensis. Reproduction. 2014;147(4):503–13.
Wilhelm D, Englert C. The Wilms tumor suppressor WT1 regulates early gonad development by activation of Sf1. Genes Dev. 2002;16(14):1839–51.
Cedar H, Bergman Y. Linking DNA methylation and histone modification: patterns and paradigms. Nat Rev Genet. 2009;10(5):295–304.
Lande-Diner L, Cedar H. Silence of the genes—mechanisms of long-term repression. Nat Rev Genet. 2005;6(8):648–54.
Foronda D, et al. Drosophila Hox and sex-determination genes control segment elimination through EGFR and extramacrochetae activity. PLoS Genet. 2012;8(8):e1002874.
Suzuki K, et al. Sexually dimorphic expression of Mafb regulates masculinization of the embryonic urethral formation. Proc Natl Acad Sci USA. 2014;111(46):16407–12.
Morrish BC, Sinclair AH. Vertebrate sex determination: many means to an end. Reproduction. 2002;124:447–57.
Rhen T, Schroeder A. Molecular mechanisms of sex determination in reptiles. Sex Dev. 2010;4(1–2):16–28.
Dhaka A, et al. TRP ion channels and temperature sensation. Annu Rev Neurosci. 2006;29:135–61.
Bogani D, et al. Loss of mitogen-activated protein kinase kinase kinase 4 (Map3k4) reveals a requirement for MAPK signalling in mouse sex determination. PLoS Biol. 2009;7:e1000196.
Warr N, et al. Gadd45γ and Map3k4 interactions regulate mouse testis determination via p38 mapk-mediated control of Sry expression. Dev Cell. 2012;23:1020–31.
Cervoni N, Szyf M. Demethylase activity is directed by histone acetylation. J Biol Chem. 2001;276(44):40778–87.
Garrison BS, et al. Postintegrative gene silencing within the Sleeping Beauty transposition system. Mol Cell Biol. 2007;27(24):8824–33.
Cridland JM, et al. Gene expression variation in Drosophila melanogaster due to rare transposable element insertion alleles of large effect. Genetics. 2015;199(1):U85–490.
McGaughey DM, et al. Genomics of CpG methylation in developing and developed zebrafish. G3. 2014;4(5):861–9.
Sati S, et al. High resolution methylome map of rat indicates role of intragenic DNA methylation in identification of coding region. PLoS ONE. 2012;7(2):e31621.
Choi SH, et al. Changes in DNA methylation of tandem DNA repeats are different from interspersed repeats in cancer. Int J Cancer. 2009;125(3):723–9.
Valenzuela N. Egg incubation and collection of painted turtle embryos. Cold Spring Harbor Protoc. 2009;4(7):1–3.
Ewert M. Embryology of turtles. In: Gans C, editor. Biology of the reptilia, development A, vol. 14. New York: Wiley; 1985.
Ewert MA, Nelson CE. Sex determination in turtles: diverse patterns and some possible adaptive values. Copeia. 1991;1991(1):50–69.
Weber M, et al. Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat Genet. 2005;37(8):853–62.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):U354–7.
Li H, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Lienhard M, et al. MEDIPS: genome-wide differential coverage analysis of sequencing data derived from DNA enrichment experiments. Bioinformatics. 2014;30(2):284–6.
Benjamini Y, Hochberg Y. Controlling the false discovery rate—a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Caetano LC, Ramos ES. MHM assay: molecular sexing based on the sex-specific methylation pattern of the MHM region in chickens. Conserv Genet. 2008;9(4):985–7.
Crowley PH. Resampling methods for computation-intensive data-analysis in ecology and evolution. Annu Rev Ecol Syst. 1992;23:405–47.
Smit AFA et al. RepeatMasker Open-3.0. 1996–2010.
Kruskal WH, Wallis WA. Use of ranks in one-criterion variance analysis. J Am Stat Assoc. 1952;47(260):583–621.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.
R Core Development Team: R: a language and environment for statistical computing. In. Version 3.2.2 edn. Vienna: R Foundation for Statistical Computing. Available via http://cran.R-project.org; 2012.
Huang DW, et al. DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35:W169–75.
Park J, et al. Comparative analyses of DNA methylation and sequence evolution using Nasonia genomes. Mol Biol Evol. 2011;28(12):3345–54.
NV conceived of the project, guided data analyses and wrote the manuscript. SR performed all bioinformatics analyses, contributed to the interpretation of results and wrote the manuscript. RL collected samples and contributed to the interpretation of results. BM contributed to data collection. All authors read and approved the final manuscript.
We thank members of the Valenzuela, Adams and Serb Labs at Iowa State University for comments.
All authors declare that they have no competing interests.
Availability of supporting data
Raw sequencing reads are available at the Short Read Archive (SRA) database of NCBI (bioproject accession PRJNA381825, study SRP103157, accessions SRR5420558-SRR5420561).
His work was partially funded by NSF Grants MCB 1244355 and IOS 1555999 to NV.
Bioinformatics and Computational Biology Program, Iowa State University, Ames, IA, 50011, USA
Srihari Radhakrishnan
Ecology and Evolutionary Biology Program, Iowa State University, Ames, IA, 50011, USA
Robert Literman
Interdepartmental Genetics and Genomics Program, Iowa State University, Ames, IA, 50011, USA
Beatriz Mizoguchi
Department of Ecology, Evolution and Organismal Biology, Iowa State University, 251 Bessey Hall, Ames, IA, 50011, USA
, Robert Literman
, Beatriz Mizoguchi
& Nicole Valenzuela
Search for Srihari Radhakrishnan in:
Search for Robert Literman in:
Search for Beatriz Mizoguchi in:
Search for Nicole Valenzuela in:
Correspondence to Nicole Valenzuela.
13072_2017_136_MOESM1_ESM.xlsx
Additional file 1. Table S1: Genes with distinct hypermethylated sex-specific windows within the same gene. Table S2: Categories enriched (p = 0.05) in hypermethylated genes at male-producing temperature (26 °C). Table S3: Categories enriched (p = 0.05) in hypermethylated genes at female-producing temperature (31 °C). Tables S4 through S11—differentially methylated genes in C. picta hatchling methylomes along with RPKM, fold change and edgeR p value: Table S4: Heat shock genes. Table S5: Androgen-/estrogen-related genes. Table S6: Kinases. Table S7: Histone-related genes. Table S8: Ubiquitin-related genes. Table S9: Transient receptor potential genes. Table S10: Genes involved in cell proliferation. Table S11: Germline-related genes. Table S12: Number of differentially methylated genes by GO term present in testis and ovaries of C. picta hatchlings. Table S13: Comprehensive list of differentially methylated genes and associated GO pathways between testis and ovaries of C. picta hatchlings. Table S14: Statistics of all differentially methylated genes between testis and ovaries of C. picta hatchlings. Table S15: Genes differentially expressed in the C. picta transcriptome (stage 22) [42] and differentially methylated in the hatchling methylome (this study). Table S16: Relative abundance distribution of repeat categories as a fraction of the genome. Table S17: Overlapping genes between C. picta hatchling methylomes and the methylomes of tongue sole [18].
Radhakrishnan, S., Literman, R., Mizoguchi, B. et al. MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination. Epigenetics & Chromatin 10, 28 (2017). https://doi.org/10.1186/s13072-017-0136-2
Received: 08 February 2017
Sex-specific thermosensitive DNA methylation
Genome-wide normalized CpG content
MeDIP sequencing
Temperature-dependent and genotypic sex determination
Turtle gonadal embryonic development
Ecological genomics
Epigenetic modification
Sexual development
Reptile vertebrate

redpajama_set_name: RedPajamaCommonCrawl


Lock in a great price for Hotel Hegra Amsterdam Centre – rated 8.4 by recent guests!
Very traditional and cozy hotel in a fantastic location in Amsterdam. Right by the canal on Herrengracht street. I booked it last minute cause I was looking for a place to stay in Amsterdam. Room 10 is on the top floor overlooking the canal. Very nice view, perfect wake up view in the morning. The owner who is also at the front desk is super friendly and polite, he will try to do his best to make your stay as pleasant as possible. He makes a small but homey breakfast buffet in a little lounge room in the back. I did not have breakfast but probably would if I did not have plans to eat outside. Overall very cozy place but most of all extremely nice owner and staff. Highly recommended, would stay there again any time.
Property was really clean, confortable and location was good. Good for visiting this amazing city.
The bed was super comfortable and the heat worked well. Bathroom was modern and wifi was good. The room was small (as to be expected) but charming.
Amazing location, great staff, relaxed atmosphere and cozy bed !
The view and location weee great! The decor of the room was adorable!
The stairway up to room 10 was a challenge but I thought it was really cool!
If one wants to have a room with "Canal View" then in this Hotel 4 rooms are meeting this wish If one wants to have a Hotel in the best possible location, then this Hotel is perfectly meeting this criteria .
Hotel Hegra Amsterdam Centre This rating is a reflection of how the property compares to the industry standard when it comes to price, facilities and services available. It's based on a self-evaluation by the property. Use this rating to help choose your stay!
Hotel Hegra Amsterdam Center offers unique rooms which differ in décor. They are fitted with wooden floors and come with flat-screen cable TV and a private bathroom with a shower, toilet and hairdryer.
The Nine Streets shopping area is only steps away and you only have to walk 7 minutes to reach the Anne Frank House. The Westermarkt tram stop is a 2-minute walk away and offers direct access to Amsterdam Central Station.
When would you like to stay at Hotel Hegra Amsterdam Centre?
Cozy room on the 3rd floor with box-spring beds.
Public parking is available at a location nearby (reservation is not possible) and costs EUR 45 per day.
Cozy rooms with box-spring bed.
This twin room features box spring beds and canal view.
House Rules Hotel Hegra Amsterdam Centre takes special requests – add in the next step!
ATM card Debit card Hotel Hegra Amsterdam Centre accepts these cards and reserves the right to temporarily hold an amount prior to arrival.
Hotel Hegra Amsterdam Centre accepts these cards and reserves the right to temporarily hold an amount prior to arrival.
Please note that for late check-in past 18:00, guests are requested to inform the accommodations with their expected time of arrival as there is no check-in possible after 23:00.
Please note that the hotel has no elevator and that you may need to climb some steep and narrow stairs to reach your room.
A breakfast buffet is available between 08:00 - 10:00 and costs EUR 12,50 per person.
Payment is required upon check-in in cash.
Pillows for me were really big and uncomfortable. Bed was comfortable.
My wife definitely left her ear rings (sentimental value) in the room, but despite a finger tip search we couldn't find them later on. Management not very helpful or sympathetic.
There were no tea or coffee making facilities in the room.
No tea making facilities in room.
Possibly the best location in Amsterdam, facing the beautiful Herengracht canal. Staff extremely friendly and helpful (e.g. booked a brilliant restaurant for us the first night). A small hotel with narrow steep stairs. Room also fairly small and basic but a great shower and warm. Good, help yourself, breakfast. Short walk to all the sights (e.g. Rijksmuseum in 15 mins).
Everything was great, nice clean room, friendly greeting and fantastic location. Would definitely recommend.
Very close to all attractions and transport! Staff were helpful and friendly!
Just the stairs to go up is a little small to pass through but francisco helped us with the bags, so its was ok!
The guys on the reception Francisco and Paulo were so kind and helpful! The room was very clean and comfortable bed. Location is great.
steep stairs but his is Amsterdam!
The location was amazing. The staff super nice and helpful (helped with the luggages down the steep stairs, here for any recommendation, etc) The bedroom was spotless, nice and warm, with comfy pillows and a decent sized bathroom.


Not one, but two! It's a gated and fenced community you're living in. With most residents earning less than the maximum wages, it's definitely a brilliant choice to stay here. It's aesthetically practical, I must say. I have stayed in blocks of flat downtown but for me, the physical wall surrounding the community impresses security and exclusivity which I would relatively prefer.
The GRILLED CHICKEN at LIM's CAFE is to die for.
It opened just recently with loads of laundry soap and eggs promotions. In the same building sits KFC, Etisalat, Arabian Grill, and jewelry shops where you can check and buy convenient stuffs.
WESTZONE is open 24 hours.
This is huge. It's like a Philippine version of hypermarket where groceries, small food court, and assorted clothing essentials are sold. Why go out when most of your needs are here? Mini-Karama and Satwa cell phone retailers are ever present for your mobile accessories and needs. Not to mention, Du and Etisalat mobile package agents catch your attention along the walkway (with free speaker and USB). Should you have an emergency, a clinic is here for a provisional check up. Pharmacy, Exchange Center, ATM kiosks, Laundry Shop, Salons, Optical Shops, and even wellness center and flower shops are apparently within reach. Westzone is life!
Run, play, sit, (or date). Al Khail Gate's park is mostly packed with muscle-finders, cricket and basketball players, lovers, frequent grass-sitters, and rovers. For families who have kids, an hour or two in the swings and slides could make an interesting, DIY Family Day.
FRIDAY MARKET happens every month.
BUS STOPS and TAXIS are anywhere.
I could play 3-4 songs from Metro to Al Khail Gate (and vice versa). It could be lengthy for some but not at all bad for many commuters like me. There are 3 bus stops (ride F15 or F17) to and fro and taxis are wide awake. Nothing to worry when you're going home late.
Commuting may be persevering but residents having a car can also be challenging. Parking inside Al Khail Gate is free – on first come, first serve basis. Otherwise, park at your own RTA risk. All the best!
SECURITY GUARDS are doing their job.
Security is stationed on each and every building. CCTV's are everywhere. Although your window is just next to a parking lot, your snooze is rest assured safe and sound.
Before, I ever thought such place is isolated, far, and labor-full! Well, it is in a way, though. I'll put it on my terms – it's private, exclusive, and collective. People live here with their own circle, contributing and accommodating other residents from many countries. From a long and tiring day, going home to such neighborhood is pleasant and relaxing. Besides, Al Khail Gate houses my heart, it is likewise where I should rather be.
Disclaimer: The views and opinions expressed in this blog are those of the author and do not necessarily reflect the community or the position of Dubai Properties.
Next Post What to Expect in Coming to Dubai?


text:
The bicycle is a curious vehicle. Its passenger is its engine.
— John Howard
Home / News / Featured Adventures / Thai Photo Pictorial
Erik Welke
Thai Photo Pictorial
Editor's Note: These photographs may look old-fashioned now, but when we first started posting photos most images were scanned in from film negatives, and the resolution on the internet was very poor. Images needed to be only a few kilobytes so that you could download them with a phone modem; now they are thousands of times bigger.
Thailand – A Photographer's Eden
By Erik Welke
Blessed with exquisite natural beauty, and possibly the friendliest population on the planet, Thailand is a joy to behold regardless of which senses you choose to employ during your visit. However, it does not yield its fruits without some amount of labor. The oppressive heat, pollution and fetid squalor of Bangkok, for example, often alienate the first-time visitor who fails to delve beneath its surface. But for those who linger and acclimate, there are rich rewards. A continuous visual feast is inevitably laid out before the eyes of anyone present to witness and, if you happen to be holding a camera, then your work is cut out for you.
River of Kings. Its murky, yellow waters plied by barges and river taxis, the Chao Phraya River leisurely snakes its way through the rapidly modernizing heart of old Bangkok. At the far right are the world famous Shangri-La and Oriental Hotels.
Monuments. The Grand Palace/Wat Phra Kaeo complex in Bangkok contains many well-preserved examples of classic Thai architecture. It was built over two hundred years ago to commemorate the origin of Bangkok as the new Thai capital. Towering in the background is the Phra Si Rattana Chedi, or Golden Stupa.
Seated Buddha. This statue of Buddha relaxes outside the Phra Mondop Library, near Wat Phra Kaeo. There are only four basic forms in which the image of Buddha is classically depicted: sitting, standing, walking or reclining. However, there are many acceptable hand and foot positions.
Golden Stupa. Entrance to the Phra Si Rattana Chedi. The gilded reliquary was erected by King Rama IV to house a piece of a bone said to be from the body of Buddha.
Water Lilies. Winged sentries guard an entrance to the Wilhan Yot, sheltering the Nak Buddha, before the tiled façade of the Auxiliary Library. The Nak Buddha and the library door panels were relocated from Ayutthaya, the previous Thai capital.
Prayer Bowls. The perpetually barefoot daughter of the author distributes 1/2 baht coins among prayer bowls, in the Temple of the Reclining Buddha.
Reclining Buddha. Wat Pho, a short walk or tuk-tuk ride from the Grand Palace, contains the largest reclining Buddha and greatest collection of Buddha images in all of Thailand, as well as the famous Institute of Thai Massage. The soles of the enormous statue's feet are adorned in mother-of-pearl with the 108 signs that are employed to recognize a living Buddha.
Stupa. Ayutthaya was the capital of Thailand (Siam) for more than 400 years until 1767 when the city was destroyed by the Burmese. This stupa is one of three at Wat Phra Si Sanphet, each in honor of a different Ayutthaya king. There are currently many bats living inside the stupa.
Clones. The enormous Chedi at Wat Yai Chai Mongkhon is surrounded by rows of seated Buddhas. At many sites in Ayutthaya, the conquering Burmese beheaded statues of Buddha of this scale. Photographing your friend with his or her head over a decapitated Buddha torso is punishable by imprisonment.
Cliffs. The Phi Phi Islands are arguably the most beautiful in all of Thailand. A long boat pilot marvels at the scenery presented, compliments of Ko Phi Phi Le, as his passenger prepares to snorkel.
Beach. This stretch of emerald water terminates on a strip of sand adjacent to the actual beach featured in the recent Leonardo Di Caprio film, "The Beach," an adaptation of the novel by British writer Alex Garland.
Fishing Village. The Muslim fishing village of Ban Ton Sai, on Ko Phi Phi Don, still holds its ground despite the growing onslaught of tourism.
Long Boats. So far, there are still no cruise ships in the serene northern bay of H-shaped Ko Phi Phi Don, but a large luxury hotel has recently been built.
PrevPreviousEgypt Photo Pictorial
NextAfrica Photo PictorialNext


Majorca, sometimes spelled Mallorca is a wonderful Balearic island off the coast of Valencia, Spain. Every year it is visited by thousands of tourists who want to enjoy the delights of the island and sample its culture. Thanks to its origins the culture is diverse and there are plenty of amazing views to savor, sunshine to soak up, and things to do and sights to see. This is made all the easier as regardless of whether you're using a Majorca taxi services firm or relying on public Majorca transport you can see everything Majorca has to offer quite easily.
Majorca is well served by taxi firms and one phone call and you can find your way traveling to Cap Formentor in the north of the island or Magaluf in the south. There are plenty of reliable taxi companies. The drivers tend to be locals with good knowledge of Majorca's quieter roads so they can use shortcuts to get you to your destination sooner. Many travel with ROIG who are both trusted and respected.
Busses are frequent on Majorca and if you don't mind longer journeys then why not take a bus to Pollenca and other delights of the island. Busses tend to go everywhere and as such journey times are longer. Nonetheless, with good planning, you can catch a bus and spend a day in a quaint town or village that you never knew existed.
The busses do not run 24 hours a day, however, so if you want to have a bit of a late night you're going to need to phone a taxi cab to get home.
As well as taxi services and busses you may want to hire a car to see the island. Hire cars firms are quite plentiful so getting behind the wheel and whizzing off to see all corners of Majorca is easily accomplished. The downside, of course, is that you're the driver so you may get lost without good map reading skills or good digital tools.
As the driver, this means you cannot enjoy the fine drinks you will find in Majorca's bars and clubs. This can be disappointing given that the nightlife culture on Majorca is strong. That said, getting around by hire car is a viable option, maybe allow a little longer for journeys unless you know the island well.
Majorca is well known for its cycling. Thanks to its comprehensive smooth road network it is easy to cycle from one point on Majorca to the other. Even if your fitness is not the best and you are staying centrally, there will be a point of interest such as a quaint town or village within easy reach.
Cycling is a great way to experience the island as you simply don't miss any of it. The island is blessed with orange and olive orchards that line the slopes of the Tramuntana Mountain Range. Cycling makes you feel part of this land.
Majorca is predominantly Spanish but it's Arabic and Roman roots are visible too. In fact, walking through some of the narrow streets of the towns and villages gives you a feeling that you are walking back in time.
There are plenty relics of the past such as the Old Town of Alcudia. Here you will see city walls built by the Romans and get a flavor of walking through narrow streets such as those of Morroco.
Take a trip across the UNESCO protected Tramuntana Mountain Range and you are in for a treat. People live on the mountain range and so visiting villages and towns is a must. Each one has its own flavor and many have great places to eat and take coffee. Some have vibrant markets to savor.
Two standout gems of the mountain range are Puig Major the highest reachable point on the island and the town of Soller.
Amazing views and equally amazing selfies are a must at the Puig Major summit. Views of the ocean and the island will take your breath away.
While the town of Soller provides a contrasting experience. Here the beach features golden soft sand and gorgeous clear water. It has rocks to explore. One of the great aspects is to reach it you descend down Sa Calobra which is a hairpin road caved into the mountains.
You can also reach it by ferry. Either way, the journey is worth the effort and a great beach day await you.
Alternatively, you can take a taxi out to Cap Formentor and savor the 'meeting point of the wind.' You can marvel at the lighthouse and enjoy views of the ocean and the whole northern side of the peninsula for a great sense of scale.
Wherever you go on the island Majorca transport has you covered. Using cabs, busses, and hire cars you can reach any place and savor it. Taxis are arguably the best transport as all you have to do is sit in the back and take it all it in. Perfect for any holiday. Why not explore these Game of Thrones filming locations while you are in Spain? Looking for more of my travel tips? Check them out here.


New Listing on RidgeClassifieds.com!!! Its fast, easy to use and best of all it's FREE to list!!!
2001 Chevy Extreme. Good condition, low mileage, custom paint, AC, stereo, automatic, full moon hubcaps.
Located in Lake Wales. Call Billy at 407-272-1622.


Part of a major £1bn urban regeneration project, Kidbrooke Village is now a well-established and highly desireable community of four distinct neighbourhoods. The Square lies at the very heart of the development and will offer residential, commercial and lifestyle amenities, plus the most convenient transport links of the entire village.
A new focal point for the community, Birch House offers the very definition of contemporary living on a grand scale. With stunning views over the Capital, it is the unique shape and style that first draw the eye. Inspired by the birch trees that line the paths of Kidbrooke Village, Birch House is the centrepiece of the area.