Genome-wide DNA methylation analysis in Chinese Chenghua and Yorkshire pigs
BMC Genomic Data volume 22, Article number: 21 (2021)
The Chinese Chenghua pig (CHP) is a typical Chinese domestic fatty pig breed with superior meat quality characteristics, while the Yorkshire pig (YP) has the characteristics of fast growth and a high rate of lean meat. Long term natural selection and artificial selection resulted in great phenotypic differences between the two breeds, including growth, development, production performance, meat quality, and coat color. However, genome-wide DNA methylation differences between CHP and YP remain unclear.
DNA methylation data were generated for muscle tissues of CHP and YP using reduced representation bisulfite sequencing (RRBS). In this study, a total of 2,416,211 CpG sites were identified. Besides, the genome-wide DNA methylation analysis revealed 722 differentially methylated regions (DMRs) and 466 differentially methylated genes (DMGs) in pairwise CHP vs. YP comparison. Six key genomic regions (Sus scrofa chromosome (SSC)1:253.47–274.23 Mb, SSC6:148.71–169.49 Mb, SSC7:0.25–9.86 Mb, SSC12:43.06–61.49 Mb, SSC14:126.43–140.95 Mb, and SSC18:49.17–54.54 Mb) containing multiple DMRs were identified, and differences of methylation patterns in these regions may be related to phenotypic differences between CHP and YP. Based on the functional analysis of DMGs, 8 DMGs (ADCY1, AGBL4, EXOC2, FUBP3, PAPPA2, PIK3R1, MGMT and MYH8) were considered as important candidate genes associated with muscle development and meat quality traits in pigs.
This study explored the difference in meat quality between CHP and YP from the epigenetic point of view, which has important reference significance for the local pork industry and pork food processing.
Epigenetic modifications of the genome can have both short-term and long-term effects on gene expression in different environments . In turn, changes in these expression profiles have implications for multiple traits. DNA methylation was the first discovered epigenetic modification and one of the most thoroughly studied . DNA methylation predominantly occurs at the C-5 position of cytosine in cytosine and guanine dinucleotide (CpG) dinucleotides in mammals . Moreover, DNA methylation is critical for mammalian growth and development . DNA methylation is traditionally regarded as a heritable and stable silence marker, which is essential for X-inactivation , silencing of genomic elements such as transposons , and genetic imprinting . In addition, variation in DNA methylation involves in a wide range of cellular functions and pathologies , and DNA methylation also affects muscle growth and development . Recently, the role of DNA methylation dynamics on skeletal muscle development and disease has been reported .
As the main meat source and human medical research model , the pig has important research value. Long-term domestication and modern breeding have resulted in both genetic variation and epigenetic modification in different breeds in pigs. Yorkshire pig (YP) is an important commercial pig breed with a high growth rate and lean meat . Chenghua pig (CHP) is a Chinese local breed which is famous for superior meat quality . By contrast, there are significant differences in body composition, muscle, and fat content between Chinese local pigs and commercial pigs , especially between CHP and YP . Epigenetic variations, and in particular DNA methylation, might not only influence differences between individuals but also between populations . Hence, DNA methylation might contribute to phenotype variation between pig breeds.
Recently, some studies have explored methylation patterns in different pig breeds and tissues. Choi and colleagues reported the DNA methylome profiles of five different tissues . Zhang and colleagues revealed the epigenetic mechanism of hypoxic adaptation in Tibetan and Yorkshire pigs . Wang and Kadarmideen provided an epigenome-wide DNA methylation map of testis by a genome-wide DNA methylation analysis . However, few studies have investigated the different epigenetic patterns between CHP and YP.
The main objective of this study was to explore the DNA methylation differences between CHP and YP by genome-wide DNA methylation analysis and then identify key genes and candidate epigenetic biomarkers associated with these differences of meat quality traits. We identified the differentially methylation regions (DMRs) and differentially methylation genes (DMGs) of CHP and YP to determine some of the important genomic regions and key genes associated with these phenotypic differences and providing new insights into the epigenetic mechanisms underlying the differences between the two pig breeds.
Meat quality traits
Meat quality traits, including pH45min, pH24h, lightness (L*), redness (a*), and yellowness (b*), were assessed at 45 min and 24 h postmortem. Table 1 summarized the meat quality traits of the two breeds. Between the CHP and YP, the pH45min (P = 7.78e-10), pH24h (P = 1.24e-4), L*45min (P = 1.95e-4), a*45min (P = 1.06e-6), and b*45min (P = 1.33e-3) of CHP were higher than those of YP.
Summary of RRBS data
Approximately 690.32Gb raw data was generated by RRBS from 48 muscle tissue samples of CHP and YP (approximately 14.38Gb raw data per individual). After quality control, approximately 523.72Gb clean data was obtained (approximately 10.91Gb clean data per individual). Besides, approximately 65% of the reads were mapped to the porcine reference genome (Table 2). Moreover, in all individuals, the density of normalized reads mapped to the proximal and distal regions of the chromosomes was higher than that of reads mapped to other regions. Overall methylated cytosines in the CpG/CHG/CHH (whereby H can be either A, T, or C) context were 51.39%/0.96%/0.7% in CHP and 52.68%/1.04%/0.78% in YP, respectively. Besides, C methylated in an unknown context like CN or CHN (whereby N can be either A, T, G, or C) was observed to be 5.8% in CHP and 5.77 in YP. Figure 1 shows CpG- and non-CpG-methylation sites (CHG, CHH, CN, or CHN) in muscle tissue of CHP and YP.
DMRs in the two groups
Table 3 summarized the numbers of CpG sites and DMRs identified by CHP vs. YP comparison. A total of 2,416,211 CpG sites and 722 DMRs were identified by CHP vs. YP comparison (Fig. 2). Of these 2,416,211 CpG sites, the distribution of 2,416,211 CpG sites annotation within promoter, exon, intron, and intergenic regions was 28.31, 15.08, 36.01, and 20.60%, respectively. Additionally, percentages of 2,416,211 sites annotated within CpG islands, CpG island shores, and other regions were 48.93, 18.96, and 32.11%. However, of these 722 DMRs in CHP vs. YP group, 3.32% were overlapped with promoter regions, 12.33% with exons, 50.69% with introns, and 33.66% with intergenic regions (Table 3). Most DMRs were in introns regions, followed by intergenic, exons, and promoters. Furthermore, 12.19% DMRs in CHP vs. YP group were CpG island regions, 16.48% in the CpG shore region, and 71.33% in other regions. Of these DMRs, much fewer (32.55% in CHP vs. YP comparison) were hypermethylated in CHP (Table S1). Six key genomic regions were identified by the CHP vs. YP comparison (Table 4). These genomic regions contained multiple DMRs shared by the CHP vs. YP comparison and DMGs. The region on SSC1 (253.47–274.23 Mb) contained 34 DMRs and 18 DMGs. The region on SSC6 (148.71–169.49 Mb) contained 18 DMRs and 8 DMGs. The region on SSC7 (0.25–9.86 Mb) contained 26 DMRs and 6 DMGs. The region on SSC12 (43.06–61.49 Mb) contained 18 DMRs and 11 DMGs. The region on SSC14 (126.43–140.95 Mb) contained 17 DMRs and 9 DMGs. The region on SSC18 (49.17–54.54 Mb) contained 11 DMRs and 6 DMGs.
DMGs identified according to DMRs and functional annotation of DMGs
We annotated 466 DMGs from DMRs identified by comparing CHP vs. YP. Besides, 149 DMGs exhibited higher levels of DNA methylation in CHP than in YP (Table S2), while 317 DMGs exhibited lower levels of DNA methylation in CHP than in YP (Table S3).
The main GO terms enriched in 466 DMGs that were identified by CHP vs. YP comparison (Fig. 3 and Table S4) and that might be related to nucleolus (P = 9.18e-3), transcription from RNA polymerase I promoter (P = 9.94e-3), DNA-directed RNA polymerase I complex (P = 1.15e-2), and membrane (P = 1.64e-2), while the KEGG pathways included Type II diabetes mellitus (P = 3.62e-3), cAMP signaling pathway (P = 4.05e-3), and Purine metabolism (P = 8.71e-3). According to the functions of DMGs, we identified 8 DMGs (ADCY1, AGBL4, EXOC2, FUBP3, PAPPA2, PIK3R1, MGMT, and MYH8) which were possibly related to the difference in appearance, meat quality, disease resistance, reproductive capacity, and adaptability between CHP and YP (Table 5).
In this study, we found that there were differences in DNA methylation between CHP and YP. The methylation patterns of CHP may help to explain the epigenetic regulation mechanisms of traits.
Bisulfite sequencing is an ideal and practical technique for studying epigenetic modifications of different species and tissues , especially DNA methylation, which can detect the DNA methylation level at each base position of the whole genome. However, genome-wide DNA methylation sequencing with high coverage of the whole genome is required to accurately assess the methylation levels at each base position. Thus, RRBS was used in this study because of its high coverage, small data requirement, low cost, and simple operation. Compared to other studies in pigs [16, 17, 19], this study used a larger population size. Therefore, RRBS is suitable for detecting DNA methylation differences among breeds in this study.
We observed several interesting GO terms and KEGG pathways associated with muscle metabolism and development. The KEGG pathways of Type II diabetes mellitus (enriched with MAPK10, PRKCE, GCK, MTOR, PIK3R1), cAMP signaling pathway (enriched with ADCY1, ACOX3, ADCY5, ARAP3, PIK3R1, MAPK10, GRIN2B, VIPR2, and VAV2), the GO terms of skeletal muscle acetylcholine-gated channel clustering (enriched with COLQ and DNAJA3) and the cAMP-mediated signaling (enriched with ADCY1, ADCY5, and KSR1) were identified in CHP vs. YP. As a major metabolic tissue, metabolic-related pathways and GO terms, including Type II diabetes mellitus, cAMP signaling pathway, skeletal muscle acetylcholine-gated channel clustering, and cAMP-mediated signaling was enriched in this study. The results indicated that DMGs associated with these metabolic processes show significant differences between CHP and YP. The pork pH has an important relationship with muscle metabolism. In this study, the pH45min (P = 7.78e-10) and pH24h (P = 1.24e-4) of CHP were higher than those of YP. Therefore, DMGs involved in muscle metabolism were identified in CHP and YP, which suggested that the difference of pH between the two breeds may be influenced by these pathways and related genes.
cAMP signaling pathway is a crucial pathway which regulates pivotal physiologic processes including metabolism, secretion, calcium homeostasis, muscle contraction, cell fate, and gene transcription. In this study, 9 DMGs are enriched in the cAMP signaling pathway. Two of these DMGs, including ADCY1 and PIK3R1, are related to melanoma metastasis. Previous studies have shown that knockdown of ADCY1 gene leads to decreased intracellular cAMP and subsequently inhibits PKA activity, and phospho- cAMP-responsive element binding protein (CREB) and microphthalmia-associated transcription factor (MITF) levels were significantly downregulated after inactivation of PKA . Furthermore, CREB and MITF have been implicated in melanoma tumor growth and metastasis [21,22,23]. Besides, the ADCY1 gene was identified as a key candidate gene involved in melanoma metastasis . There is an important link between pigmentation and melanoma. This result suggests that ADCY1 gene may affect pigmentation through cAMP. The PI3K protein, encoded by PIK3R1 gene, is a key protein involved in the PI3K/AKT signaling pathway, which is essential for myogenic differentiation  and regulates cell survival, growth, differentiation, glucose transport, and utilization . Therefore, the high levels of methylation of ADCY1 and PIK3R1 in CHP may trigger changes in their expression, potentially leading to different meat color traits between CHP and YP.
Notably, some other key DMGs, including AGBL4, EXOC2, FUBP3, PAPPA2, MGMT, and MYH8 were found in this study. The AGBL4 gene was regarded as a candidate gene associated with the heterotic quantitative trait in beef cattle . A genome-wide association study (GWAS) suggested that one SNP (rs12210050) in EXOC2 was related to the tanning ability of Europeans . A previous study demonstrated that the FUBP3 gene was associated with the skeletal formation in Duroc population . Furthermore, the FUBP3 gene was identified as a candidate gene associated with the loin eye area in pigs . The PAPPA2 gene encodes pregnancy-associated plasma protein A2 (PAPPA2) which plays an important role in the regulation of IGF-I bioavailability . It is a metalloproteinase that can specifically clew IGFBP-3 and IGFBP-5, thereby releasing IGF-I from its ternary complex, enabling it to bind to IGF-I receptors on the cell surface, initiating growth-promoting activity . Besides, in genome-wide association analysis, PAPPA2 and its related gene, PAPPA, were identified as common genetic variants associated with adult stature in the general population . The MGMT gene is a DNA repair gene responsible for removing alkylation adducts from the O6-position of guanine in DNA. The promoter CpG island hypermethylation associated gene silencing of MGMT is involved in a wide spectrum of human cancers, including glioblastoma , gastric , colon , and ovarian . The MYH8 gene belonged to the myosin heavy chain gene family that share the common features of ATP hydrolysis, actin binding, and potential for kinetic energy transduction . Moreover, the MYH8 myosin is re-expressed during muscle regeneration and is deemed as a specific marker of regenerating fibers in the pathologic skeletal muscle [39, 40].
This study performed epigenome-wide DNA methylation analysis using RRBS data generated for muscle tissues of 48 pigs. CHP vs. YP revealed 722 DMRs and 466 DMGs based on these DMRs. Besides, 6 key genomic regions and 8 key DMGs, which might be related to phenotypic differences between CHP and YP, were identified according to the further functional analysis. Our finding may help to further understand the epigenetic mechanisms of phenotype traits and have reference significance for the local pork industry.
Animals and measurements of meat quality
Totals of 48 healthy pigs were used in this study from two pig breeds, including CHP (n = 20) and YP (n = 28). These pigs were maintained in a similar environment to avoid the effects of other confounders. There are 10 males and 10 females in the Chenghua pigs, and there are 20 males and 8 females in the Yorkshire pigs. Each population contains a certain number of males and females. In addition, a large sample size was used to reduce the influence of confounders. Animals were slaughtered at a commercial slaughterhouse when they reached the slaughter weight of 105 kg. Subsequently, muscle pH values were measured at 45 min and 24 h postmortem using a portable pH meter (model 720A; Orion Research Inc., Boston, MA, USA). Meat color, including lightness (L*), redness (a*), and yellowness (b*) were assessed at 45 min and 24 h postmortem on the longissimus dorsi muscle using a Minolta CR-300 colorimeter (Minolta Camera, Osaka, Japan). After slaughter, tissue samples from muscle were collected from each breed for DNA isolation. Tissue samples were frozen in liquid nitrogen and stored at − 20 °C until analysis. All experimental procedures and sample collection were approved by the Institutional Animal Care and Use Committee of the College of Animal Science and Technology of Sichuan Agricultural University, Sichuan, China, under permit No. DKY-B20121403.
Briefly, genomic DNA was isolated from flash frozen muscular tissue. Then, the construction of RRBS libraries and paired-end sequencing using Illumina HisSeq analyzer was performed at Novogene technology co., LTD (Beijing, China). Raw sequencing data were processed by an Illumina base-calling pipeline. Genomic DNA was digested with MspI enzyme at 37 °C for 16 h. The DNA fragments after enzyme digestion were repaired at the end, and the sequencing adapters with all cytosine methylated were attached. The inserted DNA fragments with the length ranging from 40 to 220 bp were selected for glue cutting. Then, Bisulfite conversion was carried out. After that, the unmethylated C was changed to U (after PCR amplification to T), while the methylated C remained unchanged. Finally, PCR amplification was carried out to obtain the final DNA library. Clean reads were obtained from the raw data after removing reads containing adaptor sequences, unknown, or low-quality bases. The process of quality control was carried out using Trimmomatic software . Quality control was adopted to access the high data quality by (1) removing low-quality reads using a sliding window method (SLIDINGWINDOW: 4:15); (2) removing reads including adaptor sequences (ILLUMINACLIP: adapter.fa: 2:30:7:1: true); (3) removing reads with tail quality lower than 3 or with unknown bases (TRAILING: 3).
Clean reads were aligned to the pig reference genome (Sscrofa11.1) using Bismark v0.22.1 . This progress includes three steps: genome preparation, alignment using Bowtie 2 v184.108.40.206 , and methylation extractor. Bismark methylation extractor outputs read coverage and methylation percentage of detected methylated or unmethylated reads at one genomic position. The R package methyKit v1.14.2  was used to identify DMRs (window size = 500 bp, qvalue< 0.01, methylation difference > =0.25) based on the Bismark coverage file. The R package Rldeogram v0.2.2  was used to visualize the distribution of DMR. The DMGs within DMRs were annotated using the Ensemble database. The R package genomation v1.20  was used to perform annotation of DMRs. The porcine RefSeq and CpG island database (Sscrofa11.1/susScr11) for annotations were derived from the UCSC website (http://genome.ucsc.edu/cgi-bin/hgTables).
Significant GO terms and KEGG pathways were selected after filtering with P < 0.01. R package ggplot2 v3.3.2 was used to visualize the significant GO terms and KEGG pathways for the DMGs associated with DMRs.
Reduced representation bisulfite sequencing
Differentially methylated regions
Differentially methylated genes
Sus scrofa chromosome
Kyoto Encyclopedia of Genes and Genomes
Cytosine and guanine dinucleotide
Genome-wide association study
Roudbar MA, Mohammadabadi MR, Mehrgardi AA, Abdollahi-Arpanahi R, Momen M, Morota G, et al. Integration of single nucleotide variants and whole-genome DNA methylation profiles for classification of rheumatoid arthritis cases from controls. Heredity. 2020;124(5):658–74. https://doi.org/10.1038/s41437-020-0301-4.
Holliday R, Pugh JE. DNA modification mechanisms and gene activity during development. Science. 1975;187(4173):226–32. https://doi.org/10.1126/science.1111098.
Barazandeh A, Mohammadabadi M, Ghaderi-Zefrehei M, Rafeie F, Imumorin IG. Whole genome comparative analysis of CpG islands in camelid and other mammalian genomes. Mamm Biol. 2019;98(1):73–9. https://doi.org/10.1016/j.mambio.2019.07.007.
Greenberg MVC, Bourc’his D. The diverse roles of DNA methylation in mammalian development and disease. Nat Rev Mol Cell Biol. 2019;20(10):590–607. https://doi.org/10.1038/s41580-019-0159-6.
Gartler SM, Riggs AD. Mammalian X-chromosome inactivation. Annu Rev Genet. 1983;17(1):155–90. https://doi.org/10.1146/annurev.ge.17.120183.001103.
Kawakami K, Largaespada DA, Ivics Z. Transposons as tools for functional genomics in vertebrate models. Trends Genet. 2017;33(11):784–801. https://doi.org/10.1016/j.tig.2017.07.006.
Swain JL, Stewart TA, Leder P. Parental legacy determines methylation and expression of an autosomal transgene: a molecular mechanism for parental imprinting. Cell. 1987;50(5):719–27. https://doi.org/10.1016/0092-8674(87)90330-8.
Carrió E, Suelves M. DNA methylation dynamics in muscle development and disease. Front Aging Neurosci. 2015;7:19.
Schook L, Beattie C, Beever J, Donovan S, Jamison R, Zuckermann F, Niemi S, Rothschild M, Rutherford M, Smith D. Swine in biomedical research: creating the building blocks of animal models. Anim Biotechnol. 2005;16(2):183–90. https://doi.org/10.1080/10495390500265034.
Mikhaĭlov NV, Usatov AV, Getmantseva LV, Bakoev S. Associations of PRLR/AluI gene polymorphism with reproductive, growth and meat quality traits in pigs. Tsitol Genet. 2014;48(5):60–4.
Li J, Chen D, Yu B, He J, Huang Z, Mao X, et al. The fungal community and its interaction with the concentration of short-chain fatty acids in the faeces of Chenghua, Yorkshire and Tibetan pigs. Microb Biotechnol. 2020;13(2):509–21.
Yang SL, Wang ZG, Liu B, Zhang GX, Zhao SH, Yu M, et al. Genetic variation and relationships of eighteen Chinese indigenous pig breeds. Genet Sel Evol. 2003;35(6):657–71. https://doi.org/10.1186/1297-9686-35-7-657.
Yang XR, Yu B, Mao XB, Zheng P, He J, Yu J, et al. Lean and obese pig breeds exhibit differences in prenatal gene expression profiles of muscle development. Animal. 2015;9(1):28–34. https://doi.org/10.1017/S1751731114002316.
Heyn H, Moran S, Hernando-Herraez I, Sayols S, Gomez A, Sandoval J, et al. DNA methylation contributes to natural human variation. Genome Res. 2013;23(9):1363–72.
Choi M, Lee J, Le MT, Nguyen DT, Park S, Soundrarajan N, et al. Genome-wide analysis of DNA methylation in pigs using reduced representation bisulfite sequencing. DNA Res. 2015;22(5):343–55. https://doi.org/10.1093/dnares/dsv017.
Zhang B, Ban D, Gou X, Zhang Y, Yang L, Chamba Y, et al. Genome-wide DNA methylation profiles in Tibetan and Yorkshire pigs under high-altitude hypoxia. J Anim Sci Biotechnol. 2019;10:25.
Wang X, Kadarmideen HN. An epigenome-wide DNA methylation map of testis in pigs for study of complex traits. Front Genet. 2019;10:405. https://doi.org/10.3389/fgene.2019.00405.
Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, et al. Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci U S A. 2010;107(19):8689–94. https://doi.org/10.1073/pnas.1002720107.
Wang X, Kadarmideen HN. Characterization of global DNA methylation in different gene regions reveals candidate biomarkers in pigs with high and low levels of boar taint. Vet Sci. 2020;7(2):77.
Ma M, Dai J, Tang H, Xu T, Yu S, Si L, et al. MicroRNA-23a-3p inhibits mucosal melanoma growth and progression through targeting adenylate cyclase 1 and attenuating cAMP and MAPK pathways. Theranostics. 2019;9(4):945–60. https://doi.org/10.7150/thno.30516.
Dobroff AS, Wang H, Melnikova VO, Villares GJ, Zigler M, Huang L, et al. Silencing cAMP-response element-binding protein (CREB) identifies CYR61 as a tumor suppressor gene in melanoma. J Biol Chem. 2009;284(38):26194–206. https://doi.org/10.1074/jbc.M109.019836.
Xie S, Price JE, Luca M, Jean D, Ronai Z, Bar-Eli M. Dominant-negative CREB inhibits tumor growth and metastasis of human melanoma cells. Oncogene. 1997;15(17):2069–75. https://doi.org/10.1038/sj.onc.1201358.
Hsiao JJ, Fisher DE. The roles of microphthalmia-associated transcription factor and pigmentation in melanoma. Arch Biochem Biophys. 2014;563:28–34. https://doi.org/10.1016/j.abb.2014.07.019.
Chen J, Wu F, Shi Y, Yang D, Xu M, Lai Y, et al. Identification of key candidate genes involved in melanoma metastasis. Mol Med Rep. 2019;20(2):903–14. https://doi.org/10.3892/mmr.2019.10314.
Schiaffino S, Mammucari C. Regulation of skeletal muscle growth by the IGF1-Akt/PKB pathway: insights from genetic models. Skelet Muscle. 2011;1(1):4. https://doi.org/10.1186/2044-5040-1-4.
Karadoğan AH, Arikoglu H, Göktürk F, İşçioğlu F, İpekçi SH. PIK3R1 gene polymorphisms are associated with type 2 diabetes and related features in the Turkish population. Adv Clin Exp Med. 2018;27(7):921–7. https://doi.org/10.17219/acem/68985.
Akanno EC, Chen L, Abo-Ismail MK, Crowley JJ, Wang Z, Li C, et al. Genome-wide association scan for heterotic quantitative trait loci in multi-breed and crossbred beef cattle. Genet Sel Evol. 2018;50(1):48. https://doi.org/10.1186/s12711-018-0405-y.
Nan H, Kraft P, Qureshi AA, Guo Q, Chen C, Hankinson SE, et al. Genome-wide association study of tanning phenotype in a population of European ancestry. J Investig Dermatol. 2009;129(9):2250–7. https://doi.org/10.1038/jid.2009.62.
Sato S, Uemoto Y, Kikuchi T, Egawa S, Kohira K, Saito T, et al. SNP- and haplotype-based genome-wide association studies for growth, carcass, and meat quality traits in a Duroc multigenerational population. BMC Genet. 2016;17(1):60. https://doi.org/10.1186/s12863-016-0368-3.
He Y, Ma J, Zhang F, Hou L, Chen H, Guo Y, et al. Multi-breed genome-wide association study reveals heterogeneous loci associated with loin eye area in pigs. J Appl Genet. 2016;57(4):511–8. https://doi.org/10.1007/s13353-016-0351-8.
Andrew M, Liao L, Fujimoto M, Khoury J, Hwa V, Dauber A. PAPPA2 as a therapeutic modulator of IGF-I bioavailability: in vivo and in vitro evidence. J Endocr Soc. 2018;2(7):646–56. https://doi.org/10.1210/js.2018-00106.
Overgaard MT, Boldt HB, Laursen LS, Sottrup-Jensen L, Conover CA, Oxvig C. Pregnancy-associated plasma protein-A2 (PAPP-A2), a novel insulin-like growth factor-binding protein-5 proteinase. J Biol Chem. 2001;276(24):21849–53. https://doi.org/10.1074/jbc.M102191200.
Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010;467(7317):832–8. https://doi.org/10.1038/nature09410.
Mansouri A, Hachem LD, Mansouri S, Nassiri F, Laperriere NJ, Xia D, et al. MGMT promoter methylation status testing to guide therapy for glioblastoma: refining the approach based on emerging evidence and current challenges. Neuro-oncology. 2019;21(2):167–78. https://doi.org/10.1093/neuonc/noy132.
Zhang Z, Xin S, Gao M, Cai Y. Promoter hypermethylation of MGMT gene may contribute to the pathogenesis of gastric cancer: a PRISMA-compliant meta-analysis. Medicine. 2017;96(17):e6708. https://doi.org/10.1097/MD.0000000000006708.
Farzanehfar M, Vossoughinia H, Jabini R, Tavassoli A, Saadatnia H, Khorashad AK, et al. Evaluation of methylation of MGMT (O6-methylguanine-DNA methyltransferase) gene promoter in sporadic colorectal cancer. DNA Cell Biol. 2013;32(7):371–7. https://doi.org/10.1089/dna.2012.1949.
Shilpa V, Bhagat R, Premalata CS, Pallavi VR, Ramesh G, Krishnamoorthy L. Relationship between promoter methylation & tissue expression of MGMT gene in ovarian cancer. Indian J Med Res. 2014;140(5):616–23.
Yoon SJ, Seiler SH, Kucherlapati R, Leinwand L. Organization of the human skeletal myosin heavy chain gene cluster. Proc Natl Acad Sci U S A. 1992;89(24):12078–82. https://doi.org/10.1073/pnas.89.24.12078.
Schiaffino S, Rossi AC, Smerdu V, Leinwand LA, Reggiani C. Developmental myosins: expression patterns and functional significance. Skelet Muscle. 2015;5(1):22. https://doi.org/10.1186/s13395-015-0046-6.
He YM, Gu MM. Research progress of myosin heavy chain genes in human genetic diseases. Yi chuan = Hereditas. 2017;39(10):877–87. https://doi.org/10.16288/j.yczz.17-090.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Felix K, Andrews SR. Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;11:11.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.
Fan Y, Liang Y, Deng K, Zhang Z, Zhang G, Zhang Y, et al. Analysis of DNA methylation profiles during sheep skeletal muscle development using whole-genome bisulfite sequencing. BMC Genomics. 2020;21(1):327. https://doi.org/10.1186/s12864-020-6751-5.
Smith ZD, Meissner A. DNA methylation: roles in mammalian development. Nat Rev Genet. 2013;14(3):204–20. https://doi.org/10.1038/nrg3354.
Fang X, Zhao Z, Yu H, Li G, Jiang P, Yang Y, et al. Comparative genome-wide methylation analysis of longissimus dorsi muscles between Japanese black (wagyu) and Chinese red steppes cattle. PLoS One. 2017;12(8):e0182492. https://doi.org/10.1371/journal.pone.0182492.
This study was supported by High-performance Computing Platform of Sichuan Agricultural University.
The study was supported by grants from the Sichuan Science and Technology Program (2020YFN0024), the earmarked fund for the China Agriculture Research System (CARS-35-01A), the National Key R&D Program of China (2018YFD0501204), the National Natural Science Foundation of China (C170102) and the Sichuan Innovation Team of Pig (sccxtd-2021-08). The funder XL had role in methodology of this study.
Ethics approval and consent to participate
All experimental procedures and sample collection were approved by the Institutional Animal Care and Use Committee of the College of Animal Science and Technology of Sichuan Agricultural University, Sichuan, China, under permit No. DKY-B20121403.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Summary of hypermethylation regions in CHP vs. YP.
. Summary of hypomethylation regions in CHP vs. YP.
. Summary of hypermethylation genes in CHP vs. YP.
. Summary of hypomethylation genes in CHP vs. YP.
. Summary of Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in CHP vs. YP.
About this article
Cite this article
Wang, K., Wu, P., Wang, S. et al. Genome-wide DNA methylation analysis in Chinese Chenghua and Yorkshire pigs. BMC Genom Data 22, 21 (2021). https://doi.org/10.1186/s12863-021-00977-0
- DNA methylation
- Chenghua pig
- Yorkshire pig