A demonstration of the H3 trimethylation ChIP-seq analysis of galline follicular mesenchymal cells and male germ cells

  • Chokeshaiusaha, Kaj (Faculty of Veterinary Medicine, Rajamangala University of Technology Tawan-OK) ;
  • Puthier, Denis (Laboratoire TAGC/INSERM U1090, Parc Scientifique de Luminy case 928) ;
  • Nguyen, Catherine (Laboratoire TAGC/INSERM U1090, Parc Scientifique de Luminy case 928) ;
  • Sananmuang, Thanida (Faculty of Veterinary Medicine, Rajamangala University of Technology Tawan-OK)
  • Received : 2017.10.09
  • Accepted : 2018.01.10
  • Published : 2018.06.01


Objective: Trimethylation of histone 3 (H3) at 4th lysine N-termini (H3K4me3) in gene promoter region was the universal marker of active genes specific to cell lineage. On the contrary, coexistence of trimethylation at 27th lysine (H3K27me3) in the same loci-the bivalent H3K4m3/H3K27me3 was known to suspend the gene transcription in germ cells, and could also be inherited to the developed stem cell. In galline species, throughout example of H3K4m3 and H3K27me3 ChIP-seq analysis was still not provided. We therefore designed and demonstrated such procedures using ChIP-seq and mRNA-seq data of chicken follicular mesenchymal cells and male germ cells. Methods: Analytical workflow was designed and provided in this study. ChIP-seq and RNA-seq datasets of follicular mesenchymal cells and male germ cells were acquired and properly preprocessed. Peak calling by Model-based analysis of ChIP-seq 2 was performed to identify H3K4m3 or H3K27me3 enriched regions ($Fold-change{\geq}2$, $FDR{\leq}0.01$) in gene promoter regions. Integrative genomics viewer was utilized for cellular retinoic acid binding protein 1 (CRABP1), growth differentiation factor 10 (GDF10), and gremlin 1 (GREM1) gene explorations. Results: The acquired results indicated that follicular mesenchymal cells and germ cells shared several unique gene promoter regions enriched with H3K4me3 (5,704 peaks) and also unique regions of bivalent H3K4m3/H3K27me3 shared between all cell types and germ cells (1,909 peaks). Subsequent observation of follicular mesenchyme-specific genes-CRABP1, GDF10, and GREM1 correctly revealed vigorous transcriptions of these genes in follicular mesenchymal cells. As expected, bivalent H3K4m3/H3K27me3 pattern was manifested in gene promoter regions of germ cells, and thus suspended their transcriptions. Conclusion: According the results, an example of chicken H3K4m3/H3K27me3 ChIP-seq data analysis was successfully demonstrated in this study. Hopefully, the provided methodology should hereby be useful for galline ChIP-seq data analysis in the future.


Supported by : Thailand Research Fund (TRF)


  1. Feil R, Khosla S. Genomic imprinting in mammals: an interplay between chromatin and DNA methylation? Trends Genet 1999;15:431-5.
  2. Lesch BJ, Silber SJ, McCarrey JR, Page DC. Parallel evolution of male germline epigenetic poising and somatic development in animals. Nat Genet 2016;48:888-94.
  3. Waclawska A, Kurpisz M. Key functional genes of spermatogenesis identified by microarray analysis. Syst Biol Reprod Med 2012;58:229-35.
  4. Guenther MG, Levine SS, Boyer LA, Jaenisch R, Young RA. A chromatin landmark and transcription initiation at most promoters in human cells. Cell 2007;130:77-88.
  5. Zhang B, Zheng H, Huang B, et al. Allelic reprogramming of the histone modification H3K4me3 in early mammalian development. Nature 2016;537(7621):553-7.
  6. Azuara V, Perry P, Sauer S, et al. Chromatin signatures of pluripotent cell lines. Nat Cell Biol 2006;8:532-8.
  7. Mikkelsen TS, Ku M, Jaffe DB, et al. Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature 2007;448(7153):553-60.
  8. Villar D, Berthelot C, Aldridge S, et al. Enhancer evolution across 20 mammalian species. Cell 2015;160:554-66.
  9. Eissenberg JC, Shilatifard A. Histone H3 lysine 4 (H3K4) methylation in development and differentiation. Dev Biol 2010;339:240-9.
  10. Landt S, Marinov G. ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res 2012;22:1813-31.
  11. Li A, Figueroa S, Jiang T-X, et al. Diverse feather shape evolution enabled by coupling anisotropic signalling modules with selforganizing branching programme. Nat Commun 2017;8: ncomms14139.
  12. Jahan S, Xu W, He S, et al. The chicken erythrocyte epigenome. Epigenetics Chromatin 2016;9:19.
  13. Luo J, Mitra A, Tian F, et al. Histone methylation analysis and pathway predictions in chickens after MDV infection. PLoS One 2012;7:e41849.
  14. Chokeshaiusaha K, Thanawongnuwech R, Puthier D, Nguyen C. Inspection of C-type lectin superfamily expression profile in chicken and mouse dendritic cells. Thai J Vet Med 2016;46: 443-53.
  15. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods 2012;9:357-9.
  16. Dobin A, Davis CA, Schlesinger F, et al. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013;29:15-21.
  17. Li H, Handsaker B, Wysoker A, et al. The sequence alignment/ map format and SAMtools. Bioinformatics 2009;25:2078-9.
  18. Ramirez F, Dundar F, Diehl S, Gruning BA, Manke T. Deep Tools: A flexible platform for exploring deep-sequencing data. Nucleic Acids Res 2014;42 (Web Server issue):W187-91.
  19. Feng J, Liu T, Qin B, Zhang Y, Liu XS. Identifying ChIP-seq enrichment using MACS. Nat Protoc 2012;7:1728-40.
  20. Vastenhouw NL, Schier AF. Bivalent histone modifications in early embryogenesis. Curr Opin Cell Biol 2012;24:374-86.
  21. Ross-Innes CS, Stark R, Teschendorff AE, et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature 2012;184(7381):389-93.