Intronic Non-CG DNA hydroxymethylation and alternative mRNA splicing in honey bees
- Pablo Cingolani†1, 2,
- Xiaoyi Cao†3,
- Radhika S Khetani†3,
- Chieh-Chun Chen3,
- Melissa Coon1,
- Alya’a Sammak1,
- Aliccia Bollig-Fischer5,
- Susan Land1,
- Yun Huang6,
- Matthew E Hudson3, 4,
- Mark D Garfinkel7,
- Sheng Zhong3,
- Gene E Robinson3, 8, 9 and
- Douglas M Ruden1, 10Email author
© Cingolani et al.; licensee BioMed Central Ltd. 2013
Received: 11 January 2013
Accepted: 22 September 2013
Published: 30 September 2013
Previous whole-genome shotgun bisulfite sequencing experiments showed that DNA cytosine methylation in the honey bee (Apis mellifera) is almost exclusively at CG dinucleotides in exons. However, the most commonly used method, bisulfite sequencing, cannot distinguish 5-methylcytosine from 5-hydroxymethylcytosine, an oxidized form of 5-methylcytosine that is catalyzed by the TET family of dioxygenases. Furthermore, some analysis software programs under-represent non-CG DNA methylation and hydryoxymethylation for a variety of reasons. Therefore, we used an unbiased analysis of bisulfite sequencing data combined with molecular and bioinformatics approaches to distinguish 5-methylcytosine from 5-hydroxymethylcytosine. By doing this, we have performed the first whole genome analyses of DNA modifications at non-CG sites in honey bees and correlated the effects of these DNA modifications on gene expression and alternative mRNA splicing.
We confirmed, using unbiased analyses of whole-genome shotgun bisulfite sequencing (BS-seq) data, with both new data and published data, the previous finding that CG DNA methylation is enriched in exons in honey bees. However, we also found evidence that cytosine methylation and hydroxymethylation at non-CG sites is enriched in introns. Using antibodies against 5-hydroxmethylcytosine, we confirmed that DNA hydroxymethylation at non-CG sites is enriched in introns. Additionally, using a new technique, Pvu-seq (which employs the enzyme PvuRts1l to digest DNA at 5-hydroxymethylcytosine sites followed by next-generation DNA sequencing), we further confirmed that hydroxymethylation is enriched in introns at non-CG sites.
Cytosine hydroxymethylation at non-CG sites might have more functional significance than previously appreciated, and in honey bees these modifications might be related to the regulation of alternative mRNA splicing by defining the locations of the introns.
KeywordsHoney Bees DNA methylation DNA hydroxymethylation Epigenetics
Methylation of DNA is increasingly appreciated as a potent way to regulate gene expression, but a comprehensive understanding of the diversity of methylation mechanisms has not yet been achieved. For example, methylation that does not occur at cytosine-guanosine dinucleotide sequences (non-CG methylation) is an underappreciated and poorly understood form of epigenetic regulation. While rare in mammalian somatic cells, non‒CG methylation occurs on 20‒40% of the cytosines in human embryonic stem cells (hESCs)[1, 2], and is thought to be involved in pluripotency. A recent comparative analysis of DNA methylation across hESC lines found that heavily methylated non‒CG sites are strongly conserved, especially within the motif TAmCAG at 3’ splice junctions, suggesting a role in splicing or alternative splicing of mRNA transcripts. Also, CTCF induced RNA polymerase II pausing was shown to link alternative mRNA splicing to DNA methylation at a CTCF binding site in genomic DNA encoding an intron. Recently it was shown that RNA interference knockdown of DNA methyltransferase 3 (Dnmt3) affects gene alternative splicing in the honey bee.
The honey bee (Apis mellifera) is emerging as a new model to study effects of methylation on genome function because, unlike Drosophila melanogaster, it possesses a fully functioning methylation system[5–7]. Three studies[6, 8, 9] have reported that honey bees have CG methylation primarily at exon coding regions, and we have confirmed these studies here. These three studies have also reported that honey bees have little, if any, non-CG methylation[6, 8, 9]. However, since all three of the previous studies’ experimental design filtered out much of the non-CG methylation ([6, 8, 9] and personal communication; see Acknowledgements), it is still an open question as to whether there are significant amounts of non-CG methylation in bees. The are several reasons for filtering out non-CG methylation: (1) non-CG methylation is much less abundant than CG methylation in mammals; (2) there are several times more non-CG sequences (i.e., CA, CT, and CC) than CG sequences, and focusing analyses on CG methylation is simpler; (3) non-CG methylation is often in less complex regions of the genome, such as introns, and is therefore difficult to map with short-read next-generation sequencing technology; and (4) since bisulfite works less well on double-stranded DNA than single-stranded DNA, less complex regions might form snap-back structures that are resistant to bisulfite conversion.
Hydroxymethylation of DNA is a newly discovered form of epigenetic regulation. It has been found recently in embryonic stem cells (ESC) and in the brains of mammals[10, 11]. In mammals, TET proteins have been shown to be dioxygenases that convert mC to 5-hydroxymethylcytosine (hmC). Honey bees have a TET ortholog[10, 12], but hmC has not yet been reported. Genome-wide wide mapping in mouse ESCs has revealed that hmC is enriched at the start sites of genes whose promoters bear histone 3 lysine 27 trimethylation (H3K27me3) and histone 3 lysine 4 trimethylation (H3K4me3) marks[13–18]. In human ESCs, this dual mark is derived from separable subpopulations of self-renewing and lineage-biased ESCs within the heterogeneous unfractionated ESC population.
The most common chemical approach to study DNA methylation is treating single-stranded DNA with bisulfite, but bisulfite cannot distinguish mC from hmC because both base-pair as cytosine after bisulfite treatment. The modification mC mostly remains in this form after bisulfite treatment, whereas hmC is converted to cytosine methylene sulfonate (CMS) after bisulfite treatment, which has the same hydrogen bond donor and acceptor configuration as cytosine for base pairing to guanine. In this paper, in addition to whole-genome shotgun bisulfite sequencing of honey bee head DNA, we developed new biochemical and bioinformatics tools to analyze non-CG methylation and hydroxymethylation. We sharpened our analyses by comparing bees endemic to North America that are derived from a mixture of European subspecies of Apis mellifera (“EHB”) with bees derived from the African subspecies Apis mellifera scutallata, introduced to South America in the last century (“AHB”). EHB and AHB are attractive candidates for comparative molecular analysis because they differ in many physiological and behavioral traits, including aggression. Differences in methylation between EHB and AHB that are reported here may be related to these physiological and behavioral differences, and could motivate further studies, beyond the scope of the present paper. We show that non-CG methylation and hydroxymethylation are both enriched in introns, and cytosine modification at splice junctions might help regulate alternative splicing and gene expression.
Results and discussion
We asked whether non-CG methylation and hydroxymethylation occurs in honey bees using unbiased approaches that do not filter out non-CG methylation. We compared AHB and EHB mainly to demonstrate the robustness of our experimental and bioinformatics approaches. To standardize our comparisons, we analyzed head methylomes only from AHB and EHB “guard” bees, a specialized group of individuals that patrol the hive entrance and respond first to a threat to their colony. Most of the tissues in the bee head consist of brain.
Whole-genome shotgun bisulfite sequencing validates that CG methylation is primarily located in the exons
Differential cytosine DNA methylation between European (EHB) and Africanized (AHB) honey bees in CG, CHG, and CHH genomic contexts (H = A, T, or C)
Methylated in EHB (%)
Methylated in AHB (%)
Symmetrically methylated (%)
54,120 (21%, EHB), 14,454 (15%, AHB)
423 (0.53%, EHB), 106 (0.43%, AHB)
54,543 (3.4%, EHB), 14,560 (2.3% AHB)
DNA methylation at CpG sites in AHB and EHB
AHB CpG methylation by genetic region
Type (alphabetical order)
EHB CpG methylation by genetic region
Type (alphabetical order)
Whole-genome shotgun bisulfite sequencing identifies non-CG modifications that are enriched in introns
As in the previous studies[6, 8, 9], we determined that CG methylation is primarily in the exons. However, a second finding from our analyses of both our data and the data from ref. is relatively high levels of non-CG modifications (i.e., either mC or hmC) in bee heads. Surprisingly, our methods detected about 5-fold more CHH modifications (H = C, A, T) than CG methylation in both AHB and EHB DNA (Table 1). As with the CG methylation, we also saw more than twice as many non-CG modifications in EHB than AHB heads. About 2.5% of the total number of CHH sequences was modified in EHB and about 1.1% in AHB (Table 1).
In order to validate our finding of high levels of non-CG modifications, we reanalyzed previously published honey bee data from ref. with our methods and again found that the amount of non-CG modifications exceeded the amount of CG methylation (Additional file2: Table S1). The differences in the amounts of non-CG modifications in our data compared with ref. might be caused by bisulfite treatment conditions (we did a single treatment and they did multiple treatments) or differences in strains of bees used in the two studies. We constructed the Illumina libraries using the identical protocol as ref. which used unmethylated oligonucleotides, followed by amplification of only bisulfite converted oligonucleotides, to ensure that only fully bisulfite converted DNA is incorporated into the libraries (see Methods).
We also validated our data analysis pipelines by using two recent independently developed bisulfite methylation analysis program (BS-Map and Bismark)[21–23] on the data from ref. and were again able to identify non-CG methylation sites (Additional file2: Table S2 and Additional file1: Figure S2). These independent algorithms also found more non-CG than CG methylation, consistent with our findings.
DNA methylation at CHH sites in AHB and EHB
AHB CHH modifications by genetic region
Type (alphabetical order)
EHB CHH modifications by genetic region
Type (alphabetical order)
DNA methylation at CHG sites in AHB and EHB
AHB CHG methylation by genetic region
Type (alphabetical order)
EHB CHG methylation by genetic region
Type (alphabetical order)
Overall, EHB had almost 2.5× more modified Cs than did AHB (1,591,851 versus 638,400). We also observed ~3-4-fold more CG methylation than in the previous three studies: we detected 253,041 methylated CGs in EHB, compared with 80,000-90,000 in the previous 3 studies (Table 2)[6, 8, 9]. This appears to be due to higher sensitivity of the analysis program used; as stated above, using our methods on data from ref. identified 334,949 methylated CGs (Additional file2: Table S1), which is even more than we identified in our data. The significance of EHB having 2.5× more modified Cs than AHB is not known.
Bees have 5-hydroxymethylcytosine
BS-Seq cannot distinguish mC from hmC because both base pair as C after BS treatment. We used highly sensitive anti-CMS antibodies to determine the levels of hmC in the heads and bodies of AHB and EHB. Consistent with the BS-Seq results, we found comparable and statistically indistinguishable levels of hmC in EHB and AHB heads (15.2 pmol/μg and 13.5 pmol/μg of genomic DNA) (Additional file1: Figure S4). For bodies, there was significantly more hmC in EHB than AHB (25.7 pmol/μg and 19.3 pmol/μg) (p < 0.05, 2-tailed t-test, Additional file1: Figure S4).
The 5-hydroxymethylcytosine in bees is enriched in introns
Pvu-seq and HMeDIP-seq results show that cytosine hydroxymethylation is enriched in introns in the honey bee genome
Type (alphabetical order)
Type (alphabetical order)
Type (alphabetical order)
Pvu-seq validates the location of hmC in introns
There are more non-CG modifications in genes with a low CG content
CHH modifications are enriched in the introns of genes that regulate transcription
Non-CG modifications might regulate alternative mRNA splicing
Based on the above numbers, only ~0.66% of the ~56,000 splice junctions were methylated in AHB (346) and ~3.3% in EHB (1627) (Figure 4a). However, the distribution was clearly not random because pathway analyses show that genes with methylated splice sites were most enriched in the GO pathway “phosphoprotein” in both AHB and EHB (FDR < 10E-9 for both; Additional file2: Table S5). Intriguingly, the GO pathway “alternative splicing” was also significantly enriched in both AHB and EHB (FDR < 0.05 for both; Additional file2: Table S5). For example, the honey bee gene that was most heavily methylated at splice junctions is GB13778, the ortholog of Drosophila dumpy, whose protein products are involved in cell adhesion in D. melanogaster; it has four methylated splice junctions in AHB and twelve methylated splice junctions in EHB at CHH and CHG sites (Figure 4b). Since dumpy has complex splicing programs in Drosophila (16 distinct mature spliced mRNAs are listed in FlyBase), and there are dozens of dumpy exons in honey bees, it is attractive to speculate that non-CG modifications at splice junctions in honey bees regulate alternative splicing at this locus and others as well.
In hESC lines, methylated non-CG sites are strongly conserved especially within the motif 5’-TAmCAG-3’ on the non-coding DNA strand at 3’ splice junctions. While both hESCs and bees have non-CG modifications at splice junctions, bees differ from hESCs in several respects. In hESCs, methylation is symmetrical at CAG sites on both the template and non-template strand only at the 3’ splice junction, which is most frequently CAG. However in bees, methylation was primarily asymmetrical at CHH sites at the +1 position on the template strand encoding splice acceptors and the -1 position on the template strand at splice donors, and very few of the CHG sites in bees were symmetrically modified (Table 2). Another difference between hESC and bees is that the 3’ splice sites are predominantly methylated in humans but both 5’ and 3’ splice sites were modified in bees (Figure 4).
Genes with more CHH modifications in EHB than AHB are enriched in behavioral response genes
Genes with significantly more CHH modifications in EHB than AHB were enriched in GO categories that are involved in neurological functions, such as “response to external stimulus”, “substrate specific channel activity”, “exocytosis”, and “neurotransmitter receptor activity” (Additional file2: Table S6). These categories were highly significant even after correcting for the higher overall CHH modifications in EHB as well as multiple testing using the false discovery rate method (Additional file2: Table S6). It is attractive to speculate that differential CHH modifications in introns might partially explain the striking behavioral differences between AHB and EHB, especially in aggression, but this is beyond the scope of the present study. Genetic studies suggest epigenetic regulation of aggression in vertebrates[34–38], and there are extensive aggression-related differences in brain gene expression between the aggressive AHB and the less-aggressive EHB. Since there also are developmental differences between AHB and EHB, in addition to differences in aggressive behavior, these two types of differences would need to be teased apart in future studies.
Gene expression positively correlates with both mC and hmC levels in the exons
We found a weak, but significant correlation between exon methylation and exon expression. This result was obtained for methylation detected by either bisulfite sequencing or Pvu-Seq (Additional file1: Figure S9). This correlation was stronger for CG methylation than non-CG.
Exons and Splice junctions with mC or hmC appear to affect alternative mRNA splicing
As stated above, several studies have suggested that DNA methylation might be involved in regulating mRNA splicing. To determine how our new hmC finding might influence our understanding of the relationship between DNA methylation and mRNA splicing, we analyzed exons and splice junctions that differ in either mC or hmC between AHB and EHB and compared this with RNA-Seq data that we generated from AHB and EHB brains. We found several examples of differential mC and hmC associated with alternative mRNA splicing, as described below.
Our findings underscore the diversity of DNA methylation mechanisms that exist. Non-CG modifications were only recently discovered in hESC and now we report them in the distantly related honey bee. We speculate that cytosine methylation at exons and splice junctions on the DNA might affect the mRNA splicing machinery. It is important to learn how these mechanisms work in relation to known regulators of splicing, such as histone acetylation, histone 3 lysine 4 methylation, histone 3 lysine 36 methylation, and histone 3 lysine 9 methylation[43, 44].
Understanding how non-CG hydroxymethylation might affect alternative splicing is an exciting new area of research. Our data are consistent with a model in which DNA is methylated at CG sites by the maintenance DNA methyltransferase, Dnmt1, at exons, and at non-CG sites by the de novo DNA methyltransferase, Dnmt3, at introns. In contrast to Dnmt1, mammalian Dnmt3 is known to methylate both CG and non-CG sites in mammalian stem cells. We speculate that the honey bee TET enzyme primarily recognizes the non-CG sites in the introns, thereby enriching DNA hydroxymethylation in the introns. We further speculate that the mRNA splicing machinery, as well as the histone modification machinery, distinguishes exons and introns by somehow recognizing the patterns of CG methylation in exons and non-CG hydroxymethylation in introns.
Our identification of GO categories related to protein phosphorylation that were enriched for genes with methylated splice junctions is consistent with a similar finding in a recent study of species-specific alternative exons. The authors present evidence that argues that alternative splicing is used to alter protein phosphorylation, which can alter protein stability, subcellular localization, activity, and other properties. Further research is needed to determine the mechanism by which splice junction methylation and hydroxymethylation affect mRNA splicing.
All sequencing was performed using Illumina Genome Analyzer GAIIx with a Paired End Cluster Generation Kit. Image analysis, base calling and sequence extraction was performed using standard Illumina Pipeline v1.6 software. We performed whole-genome shotgun sequencing of bisulfite modified DNA (BS-Seq). DNA sequencing (>20× coverage) also was performed to ensure that the C to U conversions were BS induced, rather than natural single nucleotide polymorphisms (SNPs; not presented here). BS-Seq was performed on Africanized honey bees (AHB) (Apis mellifera scutellata) and European honey bees (EHB) (a mixture of subspecies, primarily A.m. ligustica). The number of lanes sequenced was 11 (pair-end reads) for AHB and 8 for EHB, resulting in 240 million reads for AHB and 317 million reads for EHB. Reads were 76 base-pair long yielding a total of 18.2 giga bases and 24.1 giga bases for AHB and EHB respectively. We note that we used the Illumina Whole Genome Bisulfite Sequencing (WGBS) kit which first ligates non-methylated primers to the genomic DNA prior to bisulfite conversion. After bisulfite conversion, a second set of primers are used that only amplify the fully converted primers (e.g.,). We confirmed that only fully-converted primers were amplified in the BS-seq libraries (not shown).
The bioinformatics analysis was conducted by two different groups (D.M.R. and S.Z.), using different approaches and without sharing any processed data or results.
Method 1 (BS-Miner)
The reference genome was Amel2, which is ~228 million bases long. At the time we performed our first bisulfite sequencing data analysis, there were only a few BS-Seq mappers available and some of them had a tendency to filter out non-CG methylation. For that reason, we decided to create our own unbiased analysis pipeline, which we called BS-Miner. It should be noted that nowadays many more options exist and those early mappers have been greatly improved, so there is no longer need to develop ad hoc methods. Read mapping and downstream analysis was done using BS-Miner. One of our main biological questions was whether non-CG methylation was present or not so we designed our pipeline using statistical methods well known in standard base calling algorithms. As a result, we obtain better sensitivity but only on 100%, 50% or 0% methylation levels (that is: methylated, hemi-methylated or no methylation). This sensitivity comes at a cost. As expected, the method does not detect methylation as a continuous range (from 0% to 100%) like other algorithms do. This design trade-off was aligned with our research hypothesis.
BS-Miner uses either BWA[47, 48] or Bowtie for read alignment. Both alignment programs are based on the Burrows-Wheeler transformation and create SAM output format. There are other tools and methods based on similar approaches[21, 22, 52]. In this case BWA was selected as the main mapping method in order to have better alignment near insertions and deletions. BS-Miner performs methylation calls by invoking Samtools, which uses a probabilistic model. It must be noted that the BAQ model is explicitly disabled by BS-Miner, since some of its assumptions do not apply to methylation calls.
The BcfTools package was invoked to produce methylation calls in VCF format. BS-Miner was set to filter out low quality (Q < 20) methylation calls. After all mapping and filtering steps, the mean coverage was 20.8 and 27.4 for AHB and EHB respectively, which is > 2-fold more coverage than the previously most comprehensive bee study. As a final step, BS-Miner performs several statistical analyses of the methylation results, including ranking of hypo-methylated and hyper-methylated genes by means of the Wilcoxon rank test and Fisher’s exact test. Multiple testing was corrected using False Discovery Rate methodology. Some additional statistics were carried out using custom programs in R programming language Further statistics from BS-Miner as well as additional analysis are available at http://www.mcb.mcgill.ca/~pcingola/bees/.
We also performed reanalysis of our data using Bowtie and a pipeline similar to the one shown in Krueger et al.. Read trimming was performed using Trimmomatic. Quality control was performed using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and in-house software. Performing reanalysis using different pipelines (BS-Miner and Bismark) and different filtering strategies, we obtained consistent results. Removing duplicates using Bismark’s remove duplicates module, did not seem to change our results significantly.
Method 2 (BSMap)
Using standard software, BSMap, we obtained results consistent with previous studies: methylation was primarily at CG dinucleotides in exons and very little non-CG methylation was present (Additional file1: Figure S1a). There were 61,149,121 uniquely mapped cytosines (Cs) in EHB and 53,443,185 in AHB. More than 88% of EHB Cs and 83% of AHB Cs were covered by at least two sequencing reads (Additional file1: Figure S2). However, when we included the reads that have DNA methylation in a non-CpG context, we found considerable amounts of CHH and CHG methylation (Additional file1: Figure S1b). Bisulfite reads were mapped with BSMap to distinguish the methylated cytosine from unmethylated cytosines. The reads of bisulfite converted DNA were mapped to Apis mellifera genome assembly 4 after converting the Cs to Ts. Two mismatches were allowed for the alignment to be made. To reduce the potential erroneous cytosine methylation reads, the read for a cytosine was discarded if there was a mismatch event within the 2-bp surrounding context. There was also another “CHH-filter” which filtered out the entire read if the read contains three consecutive methylated CHH. The analyses are presented in two ways: with CHH-filter and without. For every gene, its 3k upstream region, 3k downstream region, and every exon or intron was divided into 30 bins and the ratio of the number of methylated cytosines over the number of all cytosines was plotted against the bin number. Bins were graphed from upstream (relative to the gene) to downstream.
Comparison of methylation in AHB and EHB
The cytosine methylation in 3k upstream of all genes was averaged to calculate "upstream" methylation levels for every gene in both AHB and EHB; the cytosine methylation levels in gene body (or exons) were averaged to calculate the methylation level of gene body (or exons). To reduce bias from low mapping coverage, genes with less than 100 cytosines covered by BS-Seq reads, in gene body and 3k upstream, were excluded.
Validating BS-Seq with MeDIP-Seq and HMeDIP-Seq analyses
Methylated DNA immunoprecipitation followed by sequencing (MeDIP-Seq or HMeDip-Seq) was performed for a total of 41.3 million 76 bp reads. Reads were aligned using BWA and SamTools. Peak-calling was performed using MACS 1.4 beta version. Genomic DNA from 3 AHB and 3 EHB heads was sheared to 300-600 bp fragments, gel purified, immunoprecipitated with antibody, ligated to the library primers, amplified, and then sequenced. The mC antibodies were mouse monoclonal (Active Motif, Inc.) and the hmC antibodies were rabbit polyclonal (Active Motif, Inc.). Immunoprecipitation was with protein G beads (Active Motif, Inc.) following the manufacturer’s protocol. Additional file1: Figure S4a shows the peak model based on the “forward reads” [those that align with the “Watson” (plus) strand] and the “reverse reads” [those that align with the “Crick” (minus) strand].
Differential methylation analyses
Counts of methylation sites per gene, transcript, intron, splice sites and other regions of interest were calculated. Gene orthologs were calculated by InParanoid[58, 59]. Enrichment was calculated using a greedy Wilcoxon rank sum method, RssGsc (http://www.rssgsc.sourceforge.net), which also performs multiple testing correction using false discovery rate. Fisher’s exact test was also applied as a secondary method, by setting a suitable threshold in the ranked list.
Differencial SNP analysis
Counts of SNPs sites per gene, transcript, intron, splice sites and other regions of interest were calculated. Quantile normalization was applied and ranked genes analyzed for enrichment using the same methods as described in the previous section.
Observed over expected (o/e) analysis
The observed over expected ratio is defined as the number of CG dinucleotides in the reference sequence, divided by the number expected under a uniform random distribution. For genes having multiple transcripts (which are a minority in amel2 reference genome), a weighted average was calculated according to each transcript’s length. This definition is consistent with.
Accession numbers for data
The next-generation DNA sequencing experiments were done by the Applied Genomics Technology Core at Wayne State University. The next-generation DNA sequencing data (RNA-seq, Pvu-seq, MeDIP-seq, BS-seq, and HMeDIP-seq) was deposited into the GEO database according to the MINSEQE standards (Minimum Information about a high-throughput SeQuencing Experiment). The GEO database accession number is GSE50990. Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/projects/geo/).
Animal use ethical use issues
Honey bees are not a regulated invertebrate. Therefore, no ethical use approval is necessary.
GEO accession number
Africanized Honey Bee
European Honey Bee
PvuRts1I digestion of 5hmC sites followed by next generation DNA sequencing
Bisulfite treatment of DNA followed by next-generation DNA sequencing
Chromatin immunoprecipitation followed by next-generation DNA sequencing
Immunoprecipitation of 5mC DNA with anti-5mC antibodies followed by next-generation DNA sequencing
Immunoprecipitation of 5hmC DNA with anti-5hmC antibodies followed by next-generation DNA sequencing
Database for Annotation, Visualization, and Integrated Discovery
Human embryonic stem cells
Observed over expected
Variable call format
Statistical analysis of microarrays
Single nucleotide polymorphism
The reference genome for Apis mellifera (honey bee)
Model based analysis for Chip-seq
Binary call format
Cytosine followed by two nucleotides that are not guanine
Cytosine followed by one nucleotide that is not a guanine and then a guanine.
This research was funded by NSF Frontiers in Biological Research Award EF25852 (B. Schatz, PI) and NIH Director's Pioneer Award 1DP10D006416 to G.E.R., NSF DBI 09-60583 (S.Z. and G.E.R.), NIH AI44432 and a grant from the California Institute of Regenerative Medicine RM1-01729 (A.R.), NIH ES012933, ES021983 and pilot funds from the Wayne State University Office of the Vice President for Research (D.M.R.). We thank Ryszard Maleszka and Anjana Rao for critically reading the manuscript prior to submission.
One sentence summary
Honey bee non-CG DNA hydroxymethylation is enriched in introns, which supplements previous findings that honey bee CG DNA methylation is enriched in exons.
- Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM: Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009, 462 (7271): 315-322. 10.1038/nature08514.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen PY, Feng S, Joo JW, Jacobsen SE, Pellegrini M: A comparative analysis of DNA methylation across human embryonic stem cell lines. Genome Biol. 2011, 12 (7): R62-10.1186/gb-2011-12-7-r62.PubMed CentralView ArticlePubMedGoogle Scholar
- Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, Oberdoerffer P, Sandberg R, Oberdoerffer S: CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. 2011, 479 (7371): 74-79. 10.1038/nature10442.View ArticlePubMedGoogle Scholar
- Li-Byarlay H, Li Y, Stroud H, Feng SH, Newman TC, Kaneda M, Hou KK, Worley KC, Elsik CG, Wickline SA: RNA interference knockdown of DNA methyltransferase 3 affects gene alternative splicing in the honey bee. Proc Natl Acad Sci USA. 2013, 110 (31): 12750-12755. 10.1073/pnas.1310735110.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Jorda M, Jones PL, Maleszka R, Ling X, Robertson HM, Mizzen CA, Peinado MA, Robinson GE: Functional CpG methylation system in a social insect. Science. 2006, 314 (5799): 645-647. 10.1126/science.1135213.View ArticlePubMedGoogle Scholar
- Lyko F, Foret S, Kucharski R, Wolf S, Falckenhayn C, Maleszka R: The honey bee epigenomes: differential methylation of brain DNA in queens and workers. PLoS Biol. 2010, 8 (11): e1000506-10.1371/journal.pbio.1000506.PubMed CentralView ArticlePubMedGoogle Scholar
- Lyko F, Maleszka R: Insects as innovative models for functional studies of DNA methylation. Trends Genet. 2011, 27 (4): 127-131. 10.1016/j.tig.2011.01.003.View ArticlePubMedGoogle Scholar
- Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-Wide Evolutionary Analysis of Eukaryotic DNA Methylation. Science. 2010, 328 (5980): 916-919. 10.1126/science.1186366.View ArticlePubMedGoogle Scholar
- Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, Hetzel J, Jain J, Strauss SH, Halpern ME: Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci USA. 2010, 107 (19): 8689-8694. 10.1073/pnas.1002720107.PubMed CentralView ArticlePubMedGoogle Scholar
- Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, Agarwal S, Iyer LM, Liu DR, Aravind L: Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009, 324 (5929): 930-935. 10.1126/science.1170116.PubMed CentralView ArticlePubMedGoogle Scholar
- Kriaucionis S, Heintz N: The nuclear DNA base 5-hydroxymethylcytosine is present in Purkinje neurons and the brain. Science. 2009, 324 (5929): 929-930. 10.1126/science.1169786.PubMed CentralView ArticlePubMedGoogle Scholar
- Iyer LM, Tahiliani M, Rao A, Aravind L: Prediction of novel families of enzymes involved in oxidative and other complex modifications of bases in nucleic acids. Cell Cycle. 2009, 8 (11): 1698-1710. 10.4161/cc.8.11.8580.PubMed CentralView ArticlePubMedGoogle Scholar
- Pastor WA, Pape UJ, Huang Y, Henderson HR, Lister R, Ko M, McLoughlin EM, Brudno Y, Mahapatra S, Kapranov P: Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells. Nature. 2011, 473 (7347): 394-397. 10.1038/nature10102.PubMed CentralView ArticlePubMedGoogle Scholar
- Booth MJ, Branco MR, Ficz G, Oxley D, Krueger F, Reik W, Balasubramanian S: Quantitative sequencing of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution. Science. 2012, 336 (6083): 934-937. 10.1126/science.1220671.View ArticlePubMedGoogle Scholar
- Yu M, Hon GC, Szulwach KE, Song CX, Zhang L, Kim A, Li X, Dai Q, Shen Y, Park B: Base-resolution analysis of 5-hydroxymethylcytosine in the Mammalian genome. Cell. 2012, 149 (6): 1368-1380. 10.1016/j.cell.2012.04.027.PubMed CentralView ArticlePubMedGoogle Scholar
- Ficz G, Branco MR, Seisenberger S, Santos F, Krueger F, Hore TA, Marques CJ, Andrews S, Reik W: Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation. Nature. 2011, 473 (7347): 398-402. 10.1038/nature10008.View ArticlePubMedGoogle Scholar
- Williams K, Christensen J, Pedersen MT, Johansen JV, Cloos PA, Rappsilber J, Helin K: TET1 and hydroxymethylcytosine in transcription and DNA methylation fidelity. Nature. 2011, 473 (7347): 343-348. 10.1038/nature10066.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu H, D'Alessio AC, Ito S, Wang Z, Cui K, Zhao K, Sun YE, Zhang Y: Genome-wide analysis of 5-hydroxymethylcytosine distribution reveals its dual function in transcriptional regulation in mouse embryonic stem cells. Genes Dev. 2011, 25 (7): 679-684. 10.1101/gad.2036011.PubMed CentralView ArticlePubMedGoogle Scholar
- Hong SH, Rampalli S, Lee JB, McNicol J, Collins T, Draper JS, Bhatia M: Cell fate potential of human pluripotent stem cells is encoded by histone modifications. Cell stem cell. 2011, 9 (1): 24-36. 10.1016/j.stem.2011.06.002.View ArticlePubMedGoogle Scholar
- Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125 (1–2): 279-284.View ArticlePubMedGoogle Scholar
- Krueger F, Andrews SR: Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011, 27 (11): 1571-1572. 10.1093/bioinformatics/btr167.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen PY, Cokus SJ, Pellegrini M: BS Seeker: precise mapping for bisulfite sequencing. BMC Bioinforma. 2010, 11: 203-10.1186/1471-2105-11-203.View ArticleGoogle Scholar
- Xi Y, Li W: BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinforma. 2009, 10: 232-10.1186/1471-2105-10-232.View ArticleGoogle Scholar
- Wang H, Guan S, Quimby A, Cohen-Karni D, Pradhan S, Wilson G, Roberts RJ, Zhu Z, Zheng Y: Comparative characterization of the PvuRts1I family of restriction enzymes and their application in mapping genomic 5-hydroxymethylcytosine. Nucleic Acids Res. 2011, 39 (21): 9294-9305. 10.1093/nar/gkr607.PubMed CentralView ArticlePubMedGoogle Scholar
- Szwagierczak A, Brachmann A, Schmidt CS, Bultmann S, Leonhardt H, Spada F: Characterization of PvuRts1I endonuclease as a tool to investigate genomic 5-hydroxymethylcytosine. Nucleic Acids Res. 2011, 39 (12): 5149-5156. 10.1093/nar/gkr118.PubMed CentralView ArticlePubMedGoogle Scholar
- Janosi L, Yonemitsu H, Hong H, Kaji A: Molecular cloning and expression of a novel hydroxymethylcytosine-specific restriction enzyme (PvuRts1I) modulated by glucosylation of DNA. J Mol Biol. 1994, 242 (1): 45-61. 10.1006/jmbi.1994.1556.View ArticlePubMedGoogle Scholar
- Kucharski R, Maleszka J, Foret S, Maleszka R: Nutritional control of reproductive status in honeybees via DNA methylation. Science. 2008, 319 (5871): 1827-1830. 10.1126/science.1153069.View ArticlePubMedGoogle Scholar
- Wang Y, Leung FC: In silico prediction of two classes of honeybee genes with CpG deficiency or CpG enrichment and sorting according to gene ontology classes. J Mol Evol. 2009, 68 (6): 700-705. 10.1007/s00239-009-9244-3.View ArticlePubMedGoogle Scholar
- Elango N, Hunt BG, Goodisman MA, Yi SV: 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-11211. 10.1073/pnas.0900301106.PubMed CentralView ArticlePubMedGoogle Scholar
- Pradhan S, Bacolla A, Wells RD, Roberts RJ: Recombinant human DNA (cytosine-5) methyltransferase. I. Expression, purification, and comparison of de novo and maintenance methylation. J Biol Chem. 1999, 274 (46): 33002-33010. 10.1074/jbc.274.46.33002.View ArticlePubMedGoogle Scholar
- Ramsahoye BH, Biniszkiewicz D, Lyko F, Clark V, Bird AP, Jaenisch R: Non-CpG methylation is prevalent in embryonic stem cells and may be mediated by DNA methyltransferase 3a. Proc Natl Acad Sci USA. 2000, 97 (10): 5237-5242. 10.1073/pnas.97.10.5237.PubMed CentralView 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 (5): 3-10.1186/gb-2003-4-5-p3.View ArticleGoogle Scholar
- Carmon A, Topbas F, Baron M, MacIntyre RJ: dumpy interacts with a large number of genes in the developing wing of Drosophila melanogaster. Fly. 2010, 4 (2): 117-127. 10.4161/fly.4.2.11953.View ArticlePubMedGoogle Scholar
- Donaldson ZR, Young LJ: Oxytocin, vasopressin, and the neurogenetics of sociality. Science. 2008, 322 (5903): 900-904. 10.1126/science.1158668.View ArticlePubMedGoogle Scholar
- Szyf M, Weaver I, Meaney M: Maternal care, the epigenome and phenotypic differences in behavior. Reprod Toxicol. 2007, 24 (1): 9-19. 10.1016/j.reprotox.2007.05.001.View ArticlePubMedGoogle Scholar
- Weaver IC, Cervoni N, Champagne FA, D'Alessio AC, Sharma S, Seckl JR, Dymov S, Szyf M, Meaney MJ: Epigenetic programming by maternal behavior [see comment]. Nat Neurosci. 2004, 7 (8): 847-854. 10.1038/nn1276.View ArticlePubMedGoogle Scholar
- Vitaro F, Brendgen M, Boivin M, Cantin S, Dionne G, Tremblay RE, Girard A, Perusse D: A monozygotic twin difference study of friends' aggression and children's adjustment problems. Child Dev. 2011, 82 (2): 617-632. 10.1111/j.1467-8624.2010.01570.x.View ArticlePubMedGoogle Scholar
- Tremblay RE: Developmental origins of disruptive behaviour problems: the 'original sin' hypothesis, epigenetics and their consequences for prevention. J Child Psychol Psychiatry. 2010, 51 (4): 341-367. 10.1111/j.1469-7610.2010.02211.x.View ArticlePubMedGoogle Scholar
- Alaux C, Sinha S, Hasadsri L, Hunt GJ, Guzman-Novoa E, DeGrandi-Hoffman G, Uribe-Rubio JL, Southey BR, Rodriguez-Zas S, Robinson GE: Honey bee aggression supports a link between gene regulation and behavioral evolution. Proc Natl Acad Sci USA. 2009, 106 (36): 15400-15405. 10.1073/pnas.0907043106.PubMed CentralView ArticlePubMedGoogle Scholar
- Schor IE, Rascovan N, Pelisch F, Allo M, Kornblihtt AR: Neuronal cell depolarization induces intragenic chromatin modifications affecting NCAM alternative splicing. Proc Natl Acad Sci USA. 2009, 106 (11): 4325-4330. 10.1073/pnas.0810666106.PubMed CentralView ArticlePubMedGoogle Scholar
- Sims RJ, Millhouse S, Chen CF, Lewis BA, Erdjument-Bromage H, Tempst P, Manley JL, Reinberg D: Recognition of trimethylated histone H3 lysine 4 facilitates the recruitment of transcription postinitiation factors and pre-mRNA splicing. Mol Cell. 2007, 28 (4): 665-676. 10.1016/j.molcel.2007.11.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Kolasinska-Zwierz P, Down T, Latorre I, Liu T, Liu XS, Ahringer J: Differential chromatin marking of introns and expressed exons by H3K36me3. Nat Genet. 2009, 41 (3): 376-381. 10.1038/ng.322.PubMed CentralView ArticlePubMedGoogle Scholar
- Allo M, Buggiano V, Fededa JP, Petrillo E, Schor I, de la Mata M, Agirre E, Plass M, Eyras E, Elela SA: Control of alternative splicing through siRNA-mediated transcriptional gene silencing. Nat Struct Mol Biol. 2009, 16 (7): 717-724. 10.1038/nsmb.1620.View ArticlePubMedGoogle Scholar
- Saint-Andre V, Batsche E, Rachez C, Muchardt C: Histone H3 lysine 9 trimethylation and HP1gamma favor inclusion of alternative exons. Nat Struct Mol Biol. 2011, 18 (3): 337-344. 10.1038/nsmb.1995.View ArticlePubMedGoogle Scholar
- Merkin J, Russell C, Chen P, Burge CB: Evolutionary dynamics of gene and isoform regulation in Mammalian tissues. Science. 2012, 338 (6114): 1593-1599. 10.1126/science.1228186.PubMed CentralView ArticlePubMedGoogle Scholar
- Pomraning KR, Smith KM, Freitag M: Genome-wide high throughput analysis of DNA methylation in eukaryotes. Methods. 2009, 47 (3): 142-150. 10.1016/j.ymeth.2008.09.022.View ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010, 26 (5): 589-595. 10.1093/bioinformatics/btp698.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): 25-10.1186/gb-2009-10-3-r25. Bowtie is open source http://bowtie.cbcb.umd.eduView ArticleGoogle Scholar
- Burrows M, Jerian C, Lampson B, Mann T: Online Data-Compression in a Log-Structured File System. Sigplan Notices. 1992, 27 (9): 2-9. 10.1145/143371.143376.View ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen PY, Cokus SJ, Pellegrini M: BS Seeker: precise mapping for bisulfite sequencing. MC Bioinform. 2010, 11: 203-Google Scholar
- Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18 (11): 1851-1858. 10.1101/gr.078212.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Dessau RB, Pipper CB: ''R"--project for statistical computing. Ugeskr Laeger. 2008, 170 (5): 328-330.PubMedGoogle Scholar
- Krueger F, Kreck B, Franke A, Andrews SR: DNA methylome analysis using short bisulfite sequencing data. Nat Methods. 2012, 9 (2): 145-151. 10.1038/nmeth.1828.View ArticlePubMedGoogle Scholar
- Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B: RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res. 2012, 40: W622-627. 10.1093/nar/gks540.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nussbaum C, Myers RM, Brown M, Li W: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9 (9): R137-10.1186/gb-2008-9-9-r137.PubMed CentralView ArticlePubMedGoogle Scholar
- Berglund AC, Sjölund E, Östlund G, Sonnhammer ELL: InParanoid 6: eukaryotic ortholog clusters with inparalogs. Nucleic Acids Res. 2008, 36 (suppl 1): D263-PubMed CentralPubMedGoogle Scholar
- O'Brien KP, Remm M, Sonnhammer ELL: Inparanoid: a comprehensive database of eukaryotic orthologs. Nucleic Acids Res. 2005, 33 (suppl 1): D476-PubMed CentralPubMedGoogle Scholar
- Foret S, Kucharski R, Pittelkow Y, Lockett GA, Maleszka R: Epigenetic regulation of the honey bee transcriptome: unravelling the nature of methylated genes. BMC Genomics. 2009, 10: 472-10.1186/1471-2164-10-472.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.