Genome-wide analysis of DNA methylation patterns in horse
- Ja-Rang Lee†1,
- Chang Pyo Hong†2,
- Jae-Woo Moon†2,
- Yi-Deun Jung1,
- Dae-Soo Kim3,
- Tae-Hyung Kim2,
- Jeong-An Gim1,
- Jin-Han Bae1,
- Yuri Choi1,
- Jungwoo Eo1,
- Yun-Jeong Kwon1,
- Sanghoon Song2,
- Junsu Ko2,
- Young Mok Yang4,
- Hak-Kyo Lee5,
- Kyung-Do Park5,
- Kung Ahn2,
- Kyoung-Tag Do5,
- Hong-Seok Ha6,
- Kyudong Han7,
- Joo Mi Yi8,
- Hee-Jae Cha9,
- Byung-Wook Cho1,
- Jong Bhak2Email author and
- Heui-Soo Kim1Email author
© Lee et al.; licensee BioMed Central Ltd. 2014
Received: 23 August 2013
Accepted: 7 July 2014
Published: 15 July 2014
DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known.
We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density.
We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.
KeywordsThoroughbred horse Jeju horse Genome-wide DNA methylation Differential methylated region (DMR) MeDIP-seq
DNA methylation is a stably inherited epigenetic modification in eukaryotes. The regulation and characteristics of the DNA methylation still remain enigmatic, although the importance of it has been demonstrated in many biological processes such as gene expression regulation, genomic imprinting, X chromosome inactivation, maintenance of genomic stability by transposon silencing. It has also been implicated in the development of diseases such as cancer [1–7]. DNA methylation is also essential for the proper differentiation and development of mammalian tissues [8, 9]. For instance, the knockout of genes encoding the DNA-methyltransferase (DNMT) enzymes, which are responsible for de novo methylation of DNA, results in embryonic lethality in mice [10, 11]. In mammals, methycytosine is observed mostly at CpG dinucleotides, except for the CpGs in CpG islands . DNA methylation is unevenly distributed in genomes: the intergenic regions, and repetitive elements are usually hypermethylated, while the 5′ and 3′ flanking regions of genes are relatively hypomethylated compared with the intragenic regions [13–15]. Recently, whole genome methylation has been extensively examined in mammalian species [16, 17] due to advanced sequencing technologies.
Previous studies have revealed the patterns of global DNA methylation in a single or few tissues across species [18–23], or in multiple tissues or developmental stages in a single organism [8, 18, 24–28]. The DNA methylation pattern is generally conserved, and through comparative analyses of DNA methylation across mammalian species, it has been suggested to play a role in tissue-specific gene regulation . When tissue-specific differentially methylated regions (T-DMRs) in human and mouse tissues including heart, colon, kidney, testis, spleen, and muscle were compared, they could be distinguished clearly according to the corresponding tissues based on their methylation status . It is probable that there are a large number of potentially important functional differences in methylation levels across species. In primates, relative tissue methylation levels generally differ among species . However, there is insufficient evidence indicating that methylation differences exist at subspecies or breeds level.
Thoroughbred horse (TH) is a horse breed that has been manipulated by humans for improved speed, agility, and endurance in England. THs have been selected for racing ability. Thus the genetic traits related to athletic performance against TH have been extensively studied, including genotyping and transcriptome analysis [29–35]. Jeju horse (JH; a natural monument No. 347) is an indigenous Korean horse, is physically a small and rugged pony . They have been raised for meet source, farm labor, riding, and racing in Jeju Island, South Korea. Detailed genetic characterization of JH is thought to be crucial for the conservation of and for effective breeding strategies of this indigenous animal. Thus, many studies have been performed to analyses phylogenetic relationships, and discovering genetic marker [37–39]. However, until now, there have been no studies associate the traits of Jeju horse with epigenetic patterns. With the advent of next-generation sequencing (NGS) and genome-wide association studies, some studies were performed using NGS and microarray technology in thoroughbred horses [35, 40, 41]. These studies concentrated only on gene expression and genetic markers of athletic ability during and after exercise. Methylation analyses in animals exhibiting racing traits have not yet been reported. Many previous studies suggested that exercise induces methylation changes [42, 43], and athletic ability is closely associated with methylation [44, 45]. The regulation of methylation profiling related to exercise genes is important for exercising horses. Therefore, identifying methylation profiles related to exercise ability will be invaluable in studying athletic traits in racing horses. Nonetheless, there are no studies about the influence of methylation on the racing ability of TH let alone JH while the traits governing the economics of horse racing, such as the racing ability, speed, disease resistance, and recovery ability, are of important resource in the horse industry.
Here, we report the data and analyses of genome-wide DNA methylation patterns in the skeletal muscle, heart, lung, and cerebrum of TH and JH, and tissue-specific DNA methylation differences between the two horse breeds produced by methyl-DNA immunoprecipitation sequencing (MeDIP-seq).
Global methylation analysis of thoroughbred and Jeju horses
We profiled the global DNA methylation status of physicality-associated organs (skeletal muscle, heart, lung, and cerebrum) of TH and JH using MeDIP sequencing. About 21 - 24 million raw reads from each samples were sequenced resulting in on average 820 K/mm2 of cluster density, producing about 1.05 - 1.2 Gbp. After low-quality data filtration, about 81.8% - 87.5% reads, assessed as clean data, were analyzed and mapped (Additional file 1: Table S1). On average, 17.5 and 16.0 million unique mapped reads were obtained from the four tissues of TH and JH, respectively, with a high-quality read lignment against the horse reference genome (Additional file 1: Table S1).
Peak distribution in different components of the thoroughbred horse and the Jeju horse
Total peak number
Upstream 2 kb
Downstream 2 kb
Summary of methylated CGIs in the different tissues of the thoroughbred and Jeju horses
Upstream 2 kb
Downstream 2 kb
Total methylated CGIs
Differential DNA methylation in thoroughbred and Jeju horses
MeDIP-seq data validation
Analysis of functional categories of DMR-containing genes
Tissue-specific DMRs were identified by k-mean clustering in the methylation regions in the four tissues. Several genes containing DMRs were clustered, and were divided into 11 clusters (Figure 5C). The k-mean clustering of 1188 genes revealed differential methylation in each tissue (P = 0.0005 ~ 0.00051 for skeletal muscle-related clusters (MC1-MC3), P = 0.0005 for heart-related clusters (HC1-HC3), P = 0.00128 for lung-related clusters (LC1 and LC2), and P = 0.017 ~ 0.015 for cerebrum-related clusters (CC1-CC3)). Clusters of tissue-specific DMRs were located upstream of the TSS, which is 5 kb upstream of genes in the skeletal muscle and cerebrum. However, in heart and lung tissues, each cluster of DMRs was evenly spread over the region upstream and downstream of the TSS site. In heart tissue, tissue-specific DMR clusters were detected in several genes, while in the lung tissue, tissue-specific DMR clusters were detected in 85 genes.
We report the analyses and data generated by methyl-DNA immunoprecipitation sequencing to provide the genome-wide DNA methylation patterns in skeletal muscle, heart, lung, and cerebrum tissues of TH and JH. In the horse genome, gene body regions showed a higher methylation density than the intergenic regions. Also the repetitive elements had a high methylation density while CpG islands showed a low methylation density. These patterns revealed in this study were similar to those previously reported in other species, from plants to mammalians [13, 17, 50].
The promoter and 5'-UTR regions play an important role in the regulation of gene expression and they have been reported to be hypomethylated . In the case of the gene body region, except for the 5'-UTR, DNA methylation contributed to chromatin structure alteration and regulation of the transcription elongation efficiency . We report an increased level of methylation in the CDS, intron, and 3'-UTR regions in TH and JH, these results are similar to those from previously reported animal studies [22, 28]. Repeat elements occupied about 30–50% of the mammalian genome; among these, LINE elements were predominantly interspersed. In the horse genome, LINE elements were also the most predominantly interspersed repeat elements . Repeat elements are usually associated with genomic instability through structural changes such as transposition, translocation, and recombination [54, 55]. To maintain genomic stability, DNA methylation functions as a silencing mechanism for repeat elements . Thus, a major proportion of genomic methylation has been reported to occur in repeat elements, which is supported by our data. We found that DNA methylation was predominantly seen in LINE elements, consistent with findings from previous animal studies . Additionally, SINE and LTR elements were hypermethylated in the horse genome, similar to the results in other animal studies [47, 57]. Methylation of these elements is known to be a crucial factor in the maintenance of genomic stability through the suppression of transcription, transposition, and recombination . Thus, hypermethylation of repeat elements in the horse genome might play an essential role, as a defense mechanism to maintain genomic stability in the presence of active repeat elements. CpG islands have been universally reported to be regions of gene regulation via methylcytosine, possibly through the mechanism of transcriptional repression. These regions in the mammalian genome are known to be generally demethylated, in spite of having a high GC content . Intragenic and intergenic methylated CpG islands affect functional gene expression through the regulation of promoter activity, and intergenic methylated CpG islands play a crucial role in the regulation of alternative promoters and splicing . In this study, we found that about 10.73% and 10.52% of the CpG islands were methylated in TH and JH genomes, respectively, which is similar to that observed in the human genome (about 6–8%) . Further analysis of the density of methylated CpG islands in intragenic regions showed a higher methylation level in exons (11.06 ± 1.78) than in introns (1.28 ± 0.28) in the horse genome. These results were consistent with the findings in humans and rats [8, 17]. Taken together, we provide a comprehensive data and information of the whole methylome in horse, They can enable researchers to perform in depth analyses of the roles played by DNA methylation in horses and probably in other mammals.
DNA methylation is one of the main epigenetic modification mechanisms; thus, the study of DMRs within tissues or individual organisms is important. In several studies, various levels of DNA methylation could regulate tissue-specific transcription and may be important during development and differentiation . Thus, the analysis of DMRs among tissues is essential in understanding tissue specific gene expression. In particular, methylation analysis between breeds in a well-known subspecies can provides invaluable information on the evolutionary divergence and evidence for useful traits. We successfully identified differentially methylated regions within four tissues in two horse breeds. Similar results have been reported in pig tissues from various breeds [59, 60] that can be compared. Differential methylation patterns were observed in seven tissues (muscle, heart, liver, spleen, lung, kidney, and stomach) from Laiwu, a specific pig breed . In addition, the level of methylation in the liver tissue genome of other breeds of pigs (such as Berkshire, Duroc, and Landrace) also differed . Distribution patterns of DNA methylation are generally conserved among these three pig breeds, but some DMRs were detected in the coding genes and repetitive element regions in liver tissue. In this study, we also observed that distribution of DNA methylation in the two breeds showed generally conserved pattern but, some DMRs were detected a high density in the gene body, including the coding regions and introns. Gene body methylation, especially intronic DNA methylation, may be associated with alternative splicing . Thus, these results suggest that methylation has important effects on gene transcription in individual breeds. Furthermore, in the repeat region, the density of DMRs was dominant. Thus, the high density of DMRs in repeat regions could also induce differences in transcript variation and expression. In summary, differences in DNA methylation patterns and the density of DMRs in the four tissues of individual breeds may play a crucial role in the process of development and the corresponding gene expression.
Gene containing DMRs in the tissues of TH showed high representation in the categories of ATP binding and cytoskeletal protein binding. ATP binding functions play a role during exercise, as they affect ATPase activity. ATPase activity-induced ATP lysis subsequently caused intermediate molecular interactions using the energy of ATP lysis . In TH, these functions may play important roles and the dominant expression of these gene categories is required. In particular, during exercise, ATP binding could induce muscle contraction . After ATP binds the myosin head, muscle contractions are initiated due to the detachment of myosin from actin filaments . DMRs in the tissues of TH are also overrepresented in cytoskeletal protein binding. Generally, the cytoskeleton plays important roles in both intracellular transport and cellular division . In eukaryotic cells, the cytoskeleton can be classified into three types: microfilaments, intermediate filaments, and microtubules . Muscle activity-related units such as actin, keratin, and tubulin are included in the cytoskeleton. Thus, genes having DMRs could influence their binding and activities, thus differentially affecting the exercise ability in TH. These functional categories of genes with DMRs suggest that DNA methylation has an important effect on the regulation of genes categorized as being involved in ATP and cytoskeletal binding. Thus, these differences in methylation status in the tissues of TH and JH may indicate differences in their gene expression levels and exercise ability characteristics.
We have generated, for the first time, DNA methylomes for TH and JH. We provide the DNA methylation landscape and differentially methylated genomic regions in these horse species, indicating that DMRs represent comprehensive maps of the DNA methylome in TH and JH. These DNA methylome maps could be useful for further studies of epigenetic gene regulation in various horse breeds. The epigenetic system existing in the horse genome lays the foundation for studying the involvement of epigenetic modifications in horse domestication and improvement and provides a more systemic analysis of DNA methylation.
The animal protocol used in this study has been reviewed by the Pusan National University-Institutional Animal Care and Use Committee (PNU-IACUC) on their ethical procedures and scientific care, and it has been approved (Approval Number PNU-2013-0411).
Genomic DNA extraction
The healthy thoroughbred (retired racing horse, Korea Racing Authority registered number: 016222; 5 years old; a castrated horse) and Jeju horses (tested Jeju native horse breed registered number: P06071M1; 6 years old; male) were sacrificed in compliance with the international guidelines for experimental animals, and the tissues were separated and stored at -80°C. The use of these samples was approved by the National Institute of Subtropical Agriculture in Jeju Island, South Korea. Genomic DNA was isolated from 4 tissue samples from each healthy horse (skeletal muscle, heart, lung, and cerebrum) using the DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s manual for MeDIP-Seq and bisulfite-treatment experiments. DNA concentration and quality were estimated by UV spectrophotometry on a NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). For quality control, we selected only those DNA samples in which the A260/A280 ratio range was 1.6 to 2.2, the A260/A230 ratio was >1.6, and the main band was identified by agarose gel electrophoresis.
Methyl-DNA immunoprecipitation sequencing
Eight genomic DNA samples with 1 μgof genomic DNA starting material (DNA concentration of 0.1 μg/μl) were sonicated to produce DNA fragments ranging from 100 to 500 bp. After DNA end-repair and the generation of 3'-dA overhangs using the Paired-End DNA Sample Prep Kit (Illumina, San Diego, CA, USA), the DNA samples were ligated to Illumina sequencing adaptors. The fragments were denatured and then immunoprecipitated using magnetic methylated DNA immunoprecipitation kit including a 1:10 diluted antibody mix (0.3 ul antibody, 0.6 ul buffer A, and 2.10 ul distilled water) following the manufacturer’s recommendation (Diagenode, Delville, NJ, USA). The immunoprecipitated DNA was quantified by quantitative real-time PCR (qPCR). DNA fragments between 200 and 300 bp were excised from the gel and purified using a gel extraction kit (Qiagen). The products were quantified with a Quant-iTTM dsDNA High Sensitivity Assay Kit (Invitrogen, Carlsbad, CA, USA) on an Agilent 2100 Analyzer (Agilent Technologies, Santa Clara, CA, USA). After qPCR analysis, DNA libraries were subjected to paired-end sequencing with a 50-bp read length using the Illumina HiSeq 2000 platform (Illumina). After the completion of a sequencing run, raw image files were processed by Illumina Real-Time Analysis (RTA) for image analysis and base calling. Sequencing reads have been submitted to the NCBI Short Read Archive (SRA) under an SRA accession no.SRP041333.
Raw sequence data were first processed to filter out adapters and low-quality reads with the follow criteria; (1) N’s per read ≥ 10%, (2) average of quality score (QS) per read < 20, (3) number of nucleotides with < QS 20 per read ≥ 5%, and (4) having called the same bases in paired-end reads. The filtered data were then aligned to the horse reference genome (EquCab2) using the SOAPaligner (version 2.21) with mismatches of no more than 2 bp . Uniquely mapped reads were retained for further analyses. To identify genomic regions that are enriched in a pool of specifically immunoprecipitated DNA fragments, genome-wide peak scanning was carried out by MACS (version 1.4.2) with a cutoff of P-value of 1 × 10-4 to exclude false positive peaks or noises . In addition, an option of ‘--mfold’ to select the regions with MFOLD range of high-confidence enrichment ration against background to build model was used with lower limit 10 and upper limit 30. The distribution of peaks in different regions of the horse genome in each sample, including the promoter, 5'-untranslated region (UTR), 3'-UTR, exons, introns, intergenic regions, CpG islands (CGIs), and repeats, was analyzed. Methylated peaks corresponding to different genomic regions were selected by mapping at least 50% of the peak on a particular genomic region. In particular, CGI can be defined by 3 criteria: length greater than 200 bp, ≥50% GC content, ≥0.6 of CpG observed/expected ≥0.6. The methylation densities in the different regions of the genome were also compared.
To identify candidate differentially methylated regions (DMRs) in any 2 samples, their peaks were merged, and the number of reads within those peaks were assessed with chi-square and FDR statistics (P < 0.05). DMRs with a greater than 2-fold difference in read numbers were finally selected and classified as hyper- or hypo-methylated regions. All DMR-containing genes were used for subsequent gene ontology (GO) enrichment analyses using the DAVID Functional Annotation Tool with P < 0.05 . Moreover, co-existing DMRs within genes among different tissues were plotted and centered at a transcription start site (TSS) using seqMINER with the k-mean clustering method .
Bisulfite sequencing (BSP)
Three pairs of primers (Additional file 1: Table S3) were designed with MethPrimer tool (http://www.urogene.org/cgi-bin/methprimer/methprimer.cgi), including one pair for the validation of relatively high methylated region, one pair for relatively moderate methylated region, and one pair for differentially methylated regions between TH and JH. Bisulfite modification of 1 μg of genomic DNA was performed using the Imprint® DNA Modification kit by standard methods (SIGMA). The bisulfite-treated DNA was amplified by PCR with BSP specific primer pair. After a hot start, PCRs were carried out for 40 cycles of 94°C for 40 sec, 50-55°C for 40 sec, and 72°C for 40 sec. PCR products were separated on a 1.5% agarose gel, purified with the LaboPass gel extraction kit (COSMO GENETCH) and cloned into the pGEM-T-easy vector (Promega). The cloned DNA was isolated using the Plasmid DNA mini-prep kit (GeneAll). Positive clones were randomly collected for sequencing at COSMO GENETCH company (Seoul, Korea).
Availability of supporting data
All sequencing reads from this study have been submitted to the NCBI Sequence Read Archive; SRA (http://www.ncbi.nlm.nih.gov/sra/) under accession no. SRP041333.
This work was supported by a grant from the Next generation BioGreen 21 Program (No. PJ0081062011), Rural Development Administration, Republic of Korea.
- Sasaki H, Allen ND, Surani MA: DNA methylation and genomic imprinting in mammals. EXS. 1993, 64: 469-486.PubMedGoogle Scholar
- Courtier B, Heard E, Avner P: Xce haplotypes show modified methylation in a region of the active X chromosome lying 3' to Xist. Proc Natl Acad Sci U S A. 1995, 92 (8): 3531-3535.PubMed CentralPubMedView ArticleGoogle Scholar
- Siegfried Z, Eden S, Mendelsohn M, Feng X, Tsuberi BZ, Cedar H: DNA methylation represses transcription in vivo. Nat Genet. 1999, 22 (2): 203-206.PubMedView ArticleGoogle Scholar
- Bird A: DNA methylation patterns and epigenetic memory. Genes Dev. 2002, 16 (1): 6-21.PubMedView ArticleGoogle Scholar
- Robertson KD: DNA methylation and human disease. Nat Rev Genet. 2005, 6 (8): 597-610.PubMedView ArticleGoogle Scholar
- Conerly M, Grady WM: Insights into the role of DNA methylation in disease through the use of mouse models. Dis Model Mech. 2010, 3 (5–6): 290-297.PubMed CentralPubMedView ArticleGoogle Scholar
- Kulis M, Esteller M: DNA methylation and cancer. Adv Genet. 2010, 70: 27-56.PubMedView ArticleGoogle Scholar
- Illingworth R, Kerr A, Desousa D, Jorgensen H, Ellis P, Stalker J, Jackson D, Clee C, Plumb R, Rogers J, Humphray S, Cox T, Langford C, Bird A: A novel CpG island set identifies tissue-specific methylation at developmental gene loci. PLoS Biol. 2008, 6 (1): e22-PubMed CentralPubMedView ArticleGoogle Scholar
- Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33 (Suppl): 245-254.PubMedView ArticleGoogle Scholar
- Li E, Bestor TH, Jaenisch R: Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell. 1992, 69 (6): 915-926.PubMedView ArticleGoogle Scholar
- Okano M, Bell DW, Haber DA, Li E: DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell. 1999, 99 (3): 247-257.PubMedView ArticleGoogle Scholar
- Weber M, Davies JJ, Wittig D, Oakeley EJ, Haase M, Lam WL, Schubeler D: Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat Genet. 2005, 37 (8): 853-862.PubMedView ArticleGoogle Scholar
- Zhang X, Yazaki J, Sundaresan A, Cokus S, Chan SW, Chen H, Henderson IR, Shinn P, Pellegrini M, Jacobsen SE, Ecker JR: Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis. Cell. 2006, 126 (6): 1189-1201.PubMedView ArticleGoogle Scholar
- Zilberman D, Gehring M, Tran RK, Ballinger T, Henikoff S: Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription. Nat Genet. 2007, 39 (1): 61-69.PubMedView ArticleGoogle Scholar
- Gehring M, Bubb KL, Henikoff S: Extensive demethylation of repetitive elements during seed development underlies gene imprinting. Science. 2009, 324 (5933): 1447-1451.PubMed CentralPubMedView ArticleGoogle Scholar
- Ruike Y, Imanaka Y, Sato F, Shimizu K, Tsujimoto G: Genome-wide analysis of aberrant methylation in human breast cancer cells using methyl-DNA immunoprecipitation combined with high-throughput sequencing. BMC Genomics. 2010, 11: 137-PubMed CentralPubMedView ArticleGoogle Scholar
- Sati S, Tanwar VS, Kumar KA, Patowary A, Jain V, Ghosh S, Ahmad S, Singh M, Reddy SU, Chandak GR, Raghunath M, Sivasubbu S, Chakraborty K, Scaria V, Sengupta S: High resolution methylome map of rat indicates role of intragenic DNA methylation in identification of coding region. PLoS One. 2012, 7 (2): e31621-PubMed CentralPubMedView ArticleGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41 (2): 178-186.PubMed CentralPubMedView ArticleGoogle Scholar
- Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, Hetzel J, Jain J, Strauss SH, Halpern ME, Ukomadu C, Sadler KC, Pradhan S, Pellegrini M, Jacobsen SE: Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci U S A. 2010, 107 (19): 8689-8694.PubMed CentralPubMedView ArticleGoogle Scholar
- Enard W, Fassbender A, Model F, Adorjan P, Paabo S, Olek A: Differences in DNA methylation patterns between humans and chimpanzees. Curr Biol. 2004, 14 (4): R148-R149.PubMedView ArticleGoogle Scholar
- Gama-Sosa MA, Midgett RM, Slagel VA, Githens S, Kuo KC, Gehrke CW, Ehrlich M: Tissue-specific differences in DNA methylation in various mammals. Biochim Biophys Acta. 1983, 740 (2): 212-219.PubMedView ArticleGoogle Scholar
- Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010, 328 (5980): 916-919.PubMedView ArticleGoogle Scholar
- Igarashi J, Muroi S, Kawashima H, Wang X, Shinojima Y, Kitamura E, Oinuma T, Nemoto N, Song F, Ghosh S, Held WA, Nagase H: Quantitative analysis of human tissue-specific differences in methylation. Biochem Biophys Res Commun. 2008, 376 (4): 658-664.PubMed CentralPubMedView ArticleGoogle Scholar
- Weber M, Hellmann I, Stadler MB, Ramos L, Paabo S, Rebhan M, Schubeler D: Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007, 39 (4): 457-466.PubMedView ArticleGoogle Scholar
- Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Gräf S, Tomazou EM, Backdahl L, Johnson N, Herberth M, Howe KL, Jackson DK, Miretti MM, Fiegler H, Marioni JC, Birney E, Hubbard TJ, Carter NP, Tavaré S, Beck S: An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Res. 2008, 18 (9): 1518-1529.PubMed CentralPubMedView ArticleGoogle Scholar
- Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, Burton J, Cox TV, Davies R, Down TA, Haefliger C, Horton R, Howe K, Jackson DK, Kunde J, Koenig C, Liddle J, Niblett D, Otto T, Pettett R, Seemann S, Thompson C, West T, Rogers J, Olek A, Berlin K, Beck S: DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006, 38 (12): 1378-1385.PubMed CentralPubMedView ArticleGoogle Scholar
- Kitamura E, Igarashi J, Morohashi A, Hida N, Oinuma T, Nemoto N, Song F, Ghosh S, Held WA, Yoshida-Noro C, Nagase H: Analysis of tissue-specific differentially methylated regions (TDMs) in humans. Genomics. 2007, 89 (3): 326-337.PubMed CentralPubMedView ArticleGoogle Scholar
- Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai SL, Arepalli S, Dillman A, Rafferty IP, Troncoso J, Johnson R, Zielke HR, Ferrucci L, Longo DL, Cookson MR, Singleton AB: Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010, 6 (5): e1000952-PubMed CentralPubMedView ArticleGoogle Scholar
- Hill EW, Gu J, Eivers SS, Fonseca RG, McGivney BA, Govindarajan P, Orr N, Katz LM, MacHugh DE: A sequence polymorphism in MSTN predicts sprinting ability and racing stamina in thoroughbred horses. PLoS One. 2010, 5 (1): e8645-PubMed CentralPubMedView ArticleGoogle Scholar
- McGivney BA, Browne JA, Fonseca RG, Katz LM, Machugh DE, Whiston R, Hill EW: MSTN genotypes in Thoroughbred horses influence skeletal muscle gene expression and racetrack performance. Anim Genet. 2012, 43 (6): 810-812.PubMedView ArticleGoogle Scholar
- Bower MA, McGivney BA, Campana MG, Gu J, Andersson LS, Barrett E, Davis CR, Mikko S, Stock F, Voronkova V, Bradley DG, Fahey AG, Lindgren G, MacHugh DE, Sulimova G, Hill EW: The genetic origin and history of speed in the Thoroughbred racehorse. Nat Commun. 2012, 3: 643-PubMedView ArticleGoogle Scholar
- Webbon P: Harnessing the genetic toolbox for the benefit of the racing Thoroughbred. Equine Vet J. 2012, 44 (1): 8-12.PubMedView ArticleGoogle Scholar
- Hill EW, Gu J, McGivney BA, MacHugh DE: Targets of selection in the Thoroughbred genome contain exercise-relevant gene SNPs associated with elite racecourse performance. Anim Genet. 2010, 41 (Suppl 2): 56-63.PubMedView ArticleGoogle Scholar
- Hill EW, McGivney BA, Gu J, Whiston R, Machugh DE: A genome-wide SNP-association study confirms a sequence variant (g.66493737C > T) in the equine myostatin (MSTN) gene as the most powerful predictor of optimum racing distance for Thoroughbred racehorses. BMC Genomics. 2010, 11: 552-PubMed CentralPubMedView ArticleGoogle Scholar
- Park KD, Park J, Ko J, Kim BC, Kim HS, Ahn K, Do KT, Choi H, Kim HM, Song S, Lee S, Jho S, Kong HS, Yang YM, Jhun BH, Kim C, Kim TH, Hwang S, Bhak J, Lee HK, Cho BW: Whole transcriptome analyses of six thoroughbred horses before and after exercise using RNA-Seq. BMC Genomics. 2012, 13: 473-PubMed CentralPubMedView ArticleGoogle Scholar
- Cho BW, Lee KW, Kang HS, Kim SK, Shin TS, Kim YG: Application of polymerase chain reaction with short oligonucletide primers of arbitrary sequence for the genetic analysis of Cheju native horse. J Agr Tech Dev Inst. 2001, 5: 109-114.Google Scholar
- Cho GJ: Genetic Relationship and Characteristics Using microsatellite. J Life Sci. 2007, 17 (5): 699-705.View ArticleGoogle Scholar
- Kim KI, Yang YH, Lee SS, Park C, Ma R, Bouzat JL, Lewin HA: Phylogenetic relationships of Cheju horses to other horse breeds as determined by mtDNA D-loop sequence polymorphism. Anim Genet. 1999, 30 (2): 102-108.PubMedView ArticleGoogle Scholar
- Shin JA, Yang YH, Kim HS, Yun YM, Lee KK: Genetic polymorphism of the serum proteins of horses in Jeju. J Vet Sci. 2002, 3 (4): 255-263.PubMedGoogle Scholar
- Schroder W, Klostermann A, Stock KF, Distl O: A genome-wide association study for quantitative trait loci of show-jumping in Hanoverian warmblood horses. Anim Genet. 2012, 43 (4): 392-400.PubMedView ArticleGoogle Scholar
- Corbin LJ, Blott SC, Swinburne JE, Sibbons C, Fox-Clipsham LY, Helwegen M, Parkin TD, Newton JR, Bramlage LR, McIlwraith CW, Bishop SC, Woolliams JA, Vaudin M: A genome-wide association study of osteochondritis dissecans in the Thoroughbred. Mamm Genome. 2012, 23 (3–4): 294-303.PubMedView ArticleGoogle Scholar
- Barres R, Yan J, Egan B, Treebak JT, Rasmussen M, Fritz T, Caidahl K, Krook A, O'Gorman DJ, Zierath JR: Acute exercise remodels promoter methylation in human skeletal muscle. Cell Metab. 2012, 15 (3): 405-411.PubMedView ArticleGoogle Scholar
- Gomez-Pinilla F, Zhuang Y, Feng J, Ying Z, Fan G: Exercise impacts brain-derived neurotrophic factor plasticity by engaging mechanisms of epigenetic regulation. Eur J Neurosci. 2011, 33 (3): 383-390.PubMed CentralPubMedView ArticleGoogle Scholar
- Brutsaert TD, Parra EJ: What makes a champion? Explaining variation in human athletic performance. Respir Physiol Neurobiol. 2006, 151 (2–3): 109-123.PubMedView ArticleGoogle Scholar
- Terruzzi I, Senesi P, Montesano A, La Torre A, Alberti G, Benedini S, Caumo A, Fermo I, Luzi L: Genetic polymorphisms of the enzymes involved in DNA methylation and synthesis in elite athletes. Physiol Genomics. 2011, 43 (16): 965-973.PubMedView ArticleGoogle Scholar
- Hu Y, Xu H, Li Z, Zheng X, Jia X, Nie Q, Zhang X: Comparison of the genome-wide DNA methylation profiles between fast-growing and slow-growing broilers. PLoS One. 2013, 8 (2): e56411-PubMed CentralPubMedView ArticleGoogle Scholar
- Jones PA: Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012, 13 (7): 484-492.PubMedView ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.PubMedView ArticleGoogle Scholar
- Booth FW, Chakravarthy MV, Spangenburg EE: Exercise and gene expression: physiological regulation of the human genome through physical activity. J Physiol. 2002, 543 (Pt 2): 399-411.PubMed CentralPubMedView ArticleGoogle Scholar
- Li Q, Li N, Hu X, Li J, Du Z, Chen L, Yin G, Duan J, Zhang H, Zhao Y, Wang J, Li N: Genome-wide mapping of DNA methylation in chicken. PLoS One. 2011, 6 (5): e19428-PubMed CentralPubMedView ArticleGoogle Scholar
- Klose RJ, Bird AP: Genomic DNA methylation: the mark and its mediators. Trends Biochem Sci. 2006, 31 (2): 89-97.PubMedView ArticleGoogle Scholar
- Lorincz MC, Dickerson DR, Schmitt M, Groudine M: Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells. Nat Struct Mol Biol. 2004, 11 (11): 1068-1075.PubMedView ArticleGoogle Scholar
- Adelson DL, Raison JM, Garber M, Edgar RC: Interspersed repeats in the horse (Equus caballus); spatial correlations highlight conserved chromosomal domains. Anim Genet. 2010, 41 (Suppl 2): 91-99.PubMedView ArticleGoogle Scholar
- Portela A, Esteller M: Epigenetic modifications and human disease. Nat Biotechnol. 2010, 28 (10): 1057-1068.PubMedView ArticleGoogle Scholar
- Gordenin DA, Lobachev KS, Degtyareva NP, Malkova AL, Perkins E, Resnick MA: Inverted DNA repeats: a source of eukaryotic genomic instability. Mol Cell Biol. 1993, 13 (9): 5315-5322.PubMed CentralPubMedView ArticleGoogle Scholar
- Walsh CP, Chaillet JR, Bestor TH: Transcription of IAP endogenous retroviruses is constrained by cytosine methylation. Nat Genet. 1998, 20 (2): 116-117.PubMedView ArticleGoogle Scholar
- Gao F, Luo Y, Li S, Li J, Lin L, Nielsen AL, Sørensen CB, Vajta G, Wang J, Zhang X, Du Y, Yang H, Bolund L: Comparison of gene expression and genome-wide DNA methylation profiling between phenotypically normal cloned pigs and conventionally bred controls. PLoS One. 2011, 6 (10): e25901-PubMed CentralPubMedView ArticleGoogle Scholar
- Grunau C, Hindermann W, Rosenthal A: Large-scale methylation analysis of human genomic DNA reveals tissue-specific differences between the methylation profiles of genes and pseudogenes. Hum Mol Genet. 2000, 9 (18): 2651-2663.PubMedView ArticleGoogle Scholar
- Bang WY, Kim SW, Kwon SG, Hwang JH, Kim TW, Ko MS, Cho IC, Joo YK, Cho KK, Jeong JY, Kim CW: Swine liver methylomes of Berkshire, Duroc and Landrace breeds by MeDIPS. Anim Genet. 2013, 44 (4): 463-466.PubMedView ArticleGoogle Scholar
- Yang C, Zhang M, Niu W, Yang R, Zhang Y, Qiu Z, Sun B, Zhao Z: Analysis of DNA methylation in various swine tissues. PLoS One. 2011, 6 (1): e16229-PubMed CentralPubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Hirano M, Anderson DE, Erickson HP, Hirano T: Bimodal activation of SMC ATPase by intra- and inter-molecular interactions. EMBO J. 2001, 20 (12): 3238-3250.PubMed CentralPubMedView ArticleGoogle Scholar
- Karki S, Holzbaur EL: Cytoplasmic dynein and dynactin in cell division and intracellular transport. Curr Opin Cell Biol. 1999, 11 (1): 45-53.PubMedView ArticleGoogle Scholar
- Fletcher DA, Mullins RD: Cell mechanics and the cytoskeleton. Nature. 2010, 463 (7280): 485-492.PubMed CentralPubMedView ArticleGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967.PubMedView ArticleGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9 (9): R137-PubMed CentralPubMedView ArticleGoogle Scholar
- Ye T, Krebs AR, Choukrallah MA, Keime C, Plewniak F, Davidson I, Tora L: seqMINER: an integrated ChIP-seq data interpretation platform. Nucleic Acids Res. 2011, 39 (6): e35-PubMed CentralPubMedView ArticleGoogle 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 credited.