Journal cover Journal topic
Archives Animal Breeding Archiv Tierzucht
Journal topic
AAB | Articles | Volume 62, issue 1
Arch. Anim. Breed., 62, 69-76, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
Arch. Anim. Breed., 62, 69-76, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.

Original study 25 Feb 2019

Original study | 25 Feb 2019

Hypomethylation in the promoter region of ZPBP as a potential litter size indicator in Berkshire pigs

Hypomethylation in the promoter region
Sang Mi An1, Seulgi Kwon1, Jung Hye Hwang1, Go Eun Yu1, Deok Gyeong Kang1, Da Hye Park1, Tae Wan Kim1, Hwa Chun Park2, Jeongim Ha1, and Chul Wook Kim1 Sang Mi An et al.
  • 1Swine Science and Technology Center, Gyeongnam National University of Science & Technology, Jinju, 52725, South Korea
  • 2Dasan Pig Breeding Co., Namwon, 55716, South Korea
Back to toptop

In pigs, litter size is typically defined as the total number of piglets born (TNB) or the number of piglets born alive (NBA). Increasing pig litter size is of great economic interest as a means to increase productivity. The capacity of the uterus is a critical component of litter size and may play a central role in prolificacy. In this study, we investigated litter-size-related epigenetic markers in uterine tissue from Berkshire pigs with smaller litter size groups (SLGs) and larger litter size groups (LLGs) using genome-wide bisulfite sequencing (GWBS). A total of 3269 differentially methylated regions (DMRs) were identified: 1566 were hypermethylated and 1703 hypomethylated in LLG compared to SLG. The zona pellucida binding protein (ZPBP) gene was significantly hypomethylated in the LLG promoter region, and its expression was significantly upregulated in uterine tissue. Thus, the methylation status of ZPBP gene was identified as a potential indicator of litter size. Furthermore, we verified its negative correlation with litter size traits (TNB and NBA) in whole blood samples from 172 Berkshire sows as a blood-based biomarker by a porcine methylation-specific restriction enzyme polymerase chain reaction (PMP) assay. The results suggest that the methylation status of the ZPBP gene can serve as a valuable epigenetic biomarker for hyperprolific sows.

1 Introduction
Back to toptop

In commercial pig farming, increasing litter size is of great economic interest as a means to increase production (Balcells et al., 2011; Nielsen et al., 2013; Rutherford et al., 2013). Litter size is typically defined as the total number of piglets born (TNB) or the number of piglets born alive (NBA). Litter size is controlled by many factors, such as ovulation rate (number of ovulated eggs), number of corpora lutea, fertilization rate, uterine capacity, and prenatal survival (Distl, 2007; Mesa et al., 2003). Fetal survival is primarily determined by the uterine capacity of the dam, which can be defined in terms of the relative surface area of placental endometrial attachment required to support the nutrient requirements of an individual fetus throughout gestation (Wilson et al., 1998). Therefore, uterine traits can greatly influence litter size. In addition, the low heritability of litter size (5 %–10 %) suggests that proximate environmental variables may contribute significantly to variation in it (Dube et al., 2012).

DNA methylation is among the main epigenetic mechanisms and plays significant roles in gene silencing (Newell-Price et al., 2000), tissue differentiation (Laurent et al., 2010), cellular development (Smith and Meissner, 2013), X-chromosome inactivation (Pollex and Heard, 2012), and genetic imprinting (Li et al., 1993). Importantly, DNA methylation is both stably heritable and fully reversible. DNA methylation may reflect interactions between genetic and environmental factors in the development and reproduction of pigs. In particular, when DNA methylation occurs in a gene promoter, it typically acts to repress gene transcription. Several studies have suggested a correlation between differentially methylated regions (DMRs) near promoter regions and gene expression changes (Lister et al., 2009; Meissner et al., 2008; Varley et al., 2013). Furthermore, these types of dynamic changes in DNA methylation tend to occur during embryonic development as cells differentiate or become reprogrammed (Hajkova et al., 2002; Mayer et al., 2000; Sasaki and Matsui, 2008). Many recent studies have examined the genome-wide methylation profiles of livestock phenotypes associated with disease resistance, milk production, and reproduction (Congras et al., 2014; Coster et al., 2012; Jin et al., 2014; Singh et al., 2012). DNA methylation affects the expression of many genes that are critical to reproductive traits (Calicchio et al., 2014; Messerschmidt et al., 2014; Stevenson and Prendergast, 2013). Hwang et al. (2017) recently identified DMRs and differentially expressed genes (DEGs) associated with litter size in pig placentas and suggested that the PRKG2, CLCA4, and PCK1 genes play important roles in improving litter size by increasing nutrition supply through the placenta.

The objective of this study was to use epigenetic approaches to examine uterine tissues of Berkshire pigs with smaller and larger litter sizes, using genome-wide bisulfite sequencing (GWBS) technology. Our findings will provide useful knowledge and a clearer understanding of the reproductive phenotypes of individual pigs and could help in selecting sows with high fecundity for breeding.

Figure 1Comparison of average of TNB between SLG and LLG. The bars represent the mean ± SD; n=3.  Significantly different between both groups at p<0.05. TNB: total number of piglets born; SLG: smaller litter size group; LLG: larger litter size group.


2 Materials and methods
Back to toptop

2.1 Animals and sample preparation

All animal experiments were approved by the Gyeongnam National University of Science and Technology Institutional Animal Care and Use Committee (permit no. 2105-5). All Berkshire sows used in this study were reared under the same environmental conditions (Dasan Pig Breeding Co., Namwon, Korea) and provided with the same commercial diet and water ad libitum. To analyze DNA profiling according to litter size, animals were divided into smaller litter size groups (SLGs; average litter size ≤7) and larger litter size groups (LLGs; average litter size ≥10), and three sows were randomly selected from each group (Fig. 1). The uterus was collected immediately after slaughter and flushed with phosphate-buffered saline (PBS); the endometrium was then separated from the uterus. The collected endometrium was rapidly frozen in liquid nitrogen and stored at −80C until further analysis. To verify the usefulness of the candidate gene as a prognostic tool, whole blood samples were collected from 172 sows using BD Vacutainer K2 EDTA tubes (BD, Oxford, UK) as an anticoagulant. Blood was mixed immediately after the draw by inverting the tubes 10 times and then stored at 2–8 C until use.

2.2 DNA methylation

2.2.1 Genome-wide bisulfite sequencing for marker selection

Genomic DNA (gDNA) was isolated from three endometrium tissues per group using a DNeasy Tissue Kit (Qiagen, Valencia, CA, USA) and pooled for GWBS analyses. gDNA (≥6µg) was fragmented by ultrasonication to approximately 100–300 bps and then end-repaired, 3-end adenylated, and ligated with adapters. Fragmented DNA was bisulfite-converted using the EZ DNA Methylation-Gold Kit (Zymo Research, Orange, CA, USA). Bisulfite-converted DNA was quantified using a Quant-iT dsDNA High Sensitivity Assay Kit (Life Technologies, Rockville, MD, USA) on an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA) and used as a template for polymerase chain reaction (PCR) amplification. After quantitative PCR (qPCR) amplification, the resulting libraries were subjected to paired-end sequencing with a 100 bp read length using the Illumina HiSeq 2500 platform (Illumina, San Diego, CA, USA). The raw sequencing reads were cleaned by removing adaptor sequences, and reads with a percentage of unknown bases greater than 10 % and low-quality reads (more than 20 % of < Q20 bases) were filtered out to retain only high-quality reads. Clean reads were then mapped to the pig reference genome (Sscrofa 10.2) using Bismark software (version 0.9.0) with two allowed mismatches (Krueger and Andrews, 2011). Methylated cytosines were extracted from aligned reads using the Bismark methylation extractor with standard parameters. The methylation level of a cytosine (C) within an aligned read was determined by calculating the ratio of the number of reads containing a methylated C at the location to the number of all reads covering the location. DMRs between the two groups were predicted using CpG_MP with the default parameters (length, cytosine–guanine (CG) content, and cytosine : phosphate : guanine (CpG) ratio) (Su et al., 2013). We identified differentially methylated genes (DMGs) when a DMR and specific gene function element (e.g., a promoter) overlapped in the University of California Santa Cruz Genome Browser Database. Gene Ontology (GO) analysis was performed for gene functional annotation using DAVID Bioinformatics Resources (version 6.7;, last access: 13 April 2018).

Table 1Primer sequences used for PMP assay and RT-PCR.

F: forward; R: reverse.

Download Print Version | Download XLSX

2.2.2 Porcine methylation-specific restriction enzyme PCR assay

The Porcine methylation-specific restriction enzyme PCR (PMP) assay is a PCR-based methylation method that uses methylation-sensitive restriction enzymes to determine DNA methylation status. To verify the result of GWBS, whole blood samples were collected from Berkshire sows to isolate gDNA using a Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA) according to the manufacturer's instructions and digested with HpaII (NEB, Hitchin, UK) and MspI (NEB), a pair of methylation-sensitive isoschizomers that have the same recognition site (CC|GG). An undigested gDNA (5 µg) served as the negative control. Gene-specific primers were designed to flank the HpaII/MspI sites; these are described in Table 1. PCR was performed under the following conditions: 94 C for 5 min, followed by 35 cycles of 94 C for 30 s, 60 C for 30 s, and 72 C for 30 s. The products were electrophoresed on a 2 % (wv) agarose gel in 6X loading buffer (Biosesang, Seongnam, Korea).

2.3 Reverse-transcription PCR analysis

Reverse-transcription PCR (RT-PCR) analysis was performed to detect gene expression. Total RNA was extracted from three uterine tissues from each group using the TRIzol Reagent (Molecular Research Center, Cincinnati, OH, USA) and then reverse-transcribed into cDNA using Superscript II Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA). cDNA was then subjected to RT-PCR to evaluate the relative gene expression levels of the zona pellucida binding protein (ZPBP) and the gene encoding peptidylprolyl isomerase A (PPIA) (internal control), using appropriate primer pairs (Table 1). Amplification was performed using a Perkin Elmer 9700 system (Applied Biosystems, Waltham, MA, USA) under the following conditions: 95 C for 5 min; 30 cycles of 95 C for 30 s, 60 C for 30 s, 72 C for 30 s, and final elongation for 7 min at 72 C. The amplification products were separated on 2 % (wv) agarose gel and quantified using a Gel Logic model 200 imaging system (Kodak, Rochester, NY, USA).

2.4 Statistical analysis

Comparisons between groups were performed using t tests, with statistical significance determined at p<0.05. The results are expressed as means ± standard deviation (SD). Linear regression analyses were used to test the relationships between ZPBP methylation status and sow litter size traits (TNB and NBA). All statistical analyses were conducted using SPSS software (version 20.0; SPSS Inc., Chicago, IL, USA).

Table 2Summary of sequencing results and reads alignment.

Download Print Version | Download XLSX

3 Results and discussion
Back to toptop

3.1 Identification of the ZPBP gene as an epigenetic marker

GWBS was performed on gDNA from pooled uterus samples (n=3 for both SLG and LLG). In total, 1249 and 1217 million raw reads were generated in SLG and LLG, respectively. The mapped reads covered 53.08 % (SLG) and 50.48 % (LLG) of the pig genome (Table 2). A total of 3269 DMRs were discovered; 1566 were hypermethylated and 1703 were hypomethylated in LLG compared to SLG (Table 3). The DMRs were determined by considering the p value (<0.01) and false discovery rate (Q value <0.01). Among genes with DMRs, ZPBP was found to be strongly related to fecundity by GO enrichment analysis. This DMG was mainly observed at promoter regions (UP1kb; 1 kb region upstream of transcription start sites) and was hypomethylated in LLG (Table 4).

Table 3Numbers and ratio of hyper-DMRs and hypo-DMRs.

Download Print Version | Download XLSX

DNA methylation plays an important role in regulating gene expression (Wen et al., 2016). In particular, DNA methylation in gene promoters is strongly associated with gene silencing (Bell et al., 2011; Lande-Diner and Cedar, 2005; Weber et al., 2007). Recently, studies have been conducted to identify the genome-wide methylation profiles of farm animals (Hao et al., 2016; Hu et al., 2013; Xu et al., 2016; Zhang et al., 2014). Some studies have described DNA methylation for the pig uterus (Bartol et al., 2008; Franczak et al., 2017; Ko et al., 2008; Pistek et al., 2013); however, few have reported uterine genome-wide methylation patterns. GWBS, which allows unbiased genome-wide DNA methylation profiling, has been used to investigate prolificacy-related DNA methylation in unprecedented detail (Kurdyukov and Bullock, 2016). In the current study, we used GWBS to investigate the DNA methylation profiles of the genome in uterine tissues of high- and low-prolificacy pigs to explore the relationships between DNA methylation and litter size traits.

Table 4Information of DNA methylation of ZPBP.

FDR: false discovery rate; S: smaller litter size group; L: larger litter size group.

Download Print Version | Download XLSX

Figure 2Methylation analysis of the ZPBP gene in whole blood samples of Berkshire sows by PMP assay. gDNAs were cut with the methylation-sensitive restriction enzymes HpaII/MspI. M: size marker; U: undigested DNA; SLG: smaller litter size group; LLG: larger litter size group.


3.2 Confirmation of the methylation status of the ZPBP gene in whole blood by PMP assay

To assess the prognostic capability of the ZPBP gene as a blood-based biomarker for increased litter size, DNA methylation patterns of the ZPBP gene were verified in whole blood samples collected from three sows, used to obtain uterus tissues in each group by PMP assay. As shown in Fig. 2, we confirmed that the ZPBP gene was hypomethylated in LLG. This result was consistent with the results of GWBS analysis, and will allow early prediction of litter sizes in prepubertal gilts from blood samples without slaughter.

3.3 Gene expression of hypomethylated ZPBP in the promoter region

Next, to investigate the gene expression of ZPBP with hypomethylated promoter, RT-PCR was performed on uterine tissues of three Berkshire sows from each group. As shown in Fig. 4, ZPBP gene expression was significantly upregulated in LLG, as expected. Promoter methylation generally impedes the binding of transcription factors and in a second stage leads to chromatin condensation, with long-term repression of gene expression (Schubeler, 2015). Therefore, this result supports the hypothesis that ZPBP overexpression in uterine tissue is due to promoter hypomethylation.

Figure 3Correlation between methylation status of ZPBP gene and litter size traits (TNB and NBA) in 172 Berkshire sows using logistic regression analysis. The band intensities of PMP assay are shown on the x axis, and the y axis shows TNB (a) and NBA (b). TNB: total number of piglets born; NBA: number of piglets born alive: a.u: arbitrary unit.


As the DNA methylation status of promoter regions could affect gene expression through changes in chromatin structure or transcription efficiency (Klose and Bird, 2006; Lorincz et al., 2004), we compared the genome-wide methylation patterns of high- and low-prolificacy pigs to identify DMGs that might affect prolificacy traits such as litter size.

Figure 4The effect on gene expression of promoter hypomethylated ZPBP in uterine tissues by RT-PCR. The bars represent the mean ± SD; n=3.  Significantly different between both groups at p<0.05. SLG: smaller litter size group; LLG: larger litter size group.


3.4 Verification of the relationship between methylation status of the ZPBP gene and litter size traits

We verified the relationship between methylation status of the ZPBP gene and litter size traits in Berkshire sows (n=172) using a PMP assay (Fig. 4). The amplified products were normalized using undigested DNA samples. To determine the relationship, linear regression analyses were performed; the methylation status of the ZPBP gene exhibited negative relationships with litter size traits. Logistic regression analyses of the results of TNB (R2=0.1648) indicated the following: TNB (y)=-6.0069×ZPBP(x)+10.879 (Fig. 3a). When NBA was included in the model, the relationship was as follows: R2=0.1322; NBA (y)=-4.6643×ZPBP(x)+9.3843 (Fig. 3b). These results confirm that ZPBP methylation status was significantly negatively associated with litter size traits (TNB and NBA) and that sows with hypomethylated ZPBP had high fecundity. Thus, the ZPBP gene can act as an epigenetic marker for early prediction of litter size in Berkshire pigs.

Zona pellucida (ZP) is a filamentous matrix of well-structured and glycosylated glycoproteins surrounding the oocyte that acts as a morphological criterion for oocyte selection. This matrix is formed of three proteins encoded by three different genes: ZP1, ZP2, and ZP3 (Pokkyla et al., 2011; Wassarman, 2008). ZPBP, which localizes to the acrosomal membrane and likely interacts with multiple acrosomal matrix proteins, was named for its function, binding to the oocyte ZP protein following the acrosome reaction (Yu et al., 2009). ZPBP mainly acts in acrosome compaction and sperm morphogenesis during spermiogenesis. Most research on ZPBP has focused on its location in the acrosome of the sperm and its function in sperm–oocyte interactions during fertilization (Mori et al., 1993, 1995; Yu et al., 2006). However, Campbell et al. (2006) reported that murine ZPBP was a luminal epithelium-specific gene with 20-fold or higher expression in the uterine luminal epithelium than in the stroma–glandular epithelium. Implantation is essential for the establishment of normal pregnancy and is initiated by a physical interaction between the trophoblast and the apical surface of the luminal epithelium, followed rapidly by adhesion and then by penetration through the luminal epithelium to the underlying stroma, which responds by decidualization (Abrahamsohn and Zorn, 1993). In other words, the uterine luminal epithelium plays a critical role in implantation. ZPBP is a serine protease (serine protease 38) with ZP-binding properties that were initially identified in porcine epididymal sperm (Mori et al., 1993). Serine proteases are characterized by the presence of serine as a nucleophilic amino acid at the active site of the enzyme (Hedstrom, 2002). Some serine proteases are detectable in the uterus and are involved in female reproduction, especially in oocyte development, ovulation, implantation, and decidualization (Diao et al., 2013; Nie et al., 2005). Therefore, we believe that ZPBP might play a role in female reproduction, including implantation.

4 Conclusions
Back to toptop

This study reports the DNA methylation patterns of porcine uterine tissues, which are associated with litter size. We identified DMRs and detected hypomethylation of the ZPBP gene in its promoter region in LLG; its expression was upregulated in uterine tissues. We also verified that the methylation status of the ZPBP gene was significantly negatively associated with litter size traits in the larger pig population. Our results demonstrate that this gene can be used as a biomarker for hyperprolific sows and will likely contribute to improving reproductive capacity.

Data availability
Back to toptop
Data availability. 

The data sets are available upon request from the corresponding author.

Author contributions
Back to toptop
Author contributions. 

SMA curated data and wrote this paper. SK, JHH and GEY performed the experiments. DGK, DHP and TWK analyzed the data. HCP provided the resources. JH and CWK supervised the research project.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

Research was supported by grants from the Basic Science Research Program (no. 2017R1A6A3A11035414) and Priority Research Centers Program (no. 2011-0022965) through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology and the Export Promotion Technology Development Program (no. 313012-05) of the Ministry of Food, Agriculture, Forestry and Fisheries, Republic of Korea.

Edited by: Steffen Maak
Reviewed by: three anonymous referees

Back to toptop

Abrahamsohn, P. A. and Zorn, T. M.: Implantation and decidualization in rodents, J. Exp. Zool., 266, 603–628, 1993. 

Balcells, I., Castello, A., Mercade, A., Noguera, J. L., Fernandez-Rodriguez, A., Sanchez, A., and Tomas, A.: Analysis of porcine MUC4 gene as a candidate gene for prolificacy QTL on SSC13 in an Iberian x Meishan F2 population, BMC Genet., 12, 1–6,, 2011. 

Bartol, F. F., Wiley, A. A., and Bagnell, C. A.: Epigenetic programming of porcine endometrial function and the lactocrine hypothesis, Reprod. Domest. Anim., 43, 273–279,, 2008. 

Bell, J. T., Pai, A. A., Pickrell, J. K., Gaffney, D. J., Pique-Regi, R., Degner, J. F., Gilad, Y., and Pritchard, J. K.: DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines, Genome Biol., 12,, 2011. 

Calicchio, R., Doridot, L., Miralles, F., Mehats, C., and Vaiman, D.: DNA methylation, an epigenetic mode of gene expression regulation in reproductive science, Curr. Pharm. Design, 20, 1726–1750, 2014. 

Campbell, E. A., O'Hara, L., Catalano, R. D., Sharkey, A. M., Freeman, T. C., and Johnson, M. H.: Temporal expression profiling of the uterine luminal epithelium of the pseudo-pregnant mouse suggests receptivity to the fertilized egg is associated with complex transcriptional changes, Hum. Reprod., 21, 2495–2513, 2006.  

Congras, A., Yerle-Bouissou, M., Pinton, A., Vignoles, F., Liaubet, L., Ferchaud, S., and Acloque, H.: Sperm DNA methylation analysis in swine reveals conserved and species-specific methylation patterns and highlights an altered methylation at the GNAS locus in infertile boars, Biol. Reprod., 91, 137, 1–14, 2014. 

Coster, A., Madsen, O., Heuven, H. C., Dibbits, B., Groenen, M. A., van Arendonk, J. A., and Bovenhuis, H.: The imprinted gene DIO3 is a candidate gene for litter size in pigs, Plos One, 7, e31825,, 2012. 

Diao, H., Xiao, S., Li, R., Zhao, F., and Ye, X.: Distinct spatiotemporal expression of serine proteases Prss23 and Prss35 in periimplantation mouse uterus and dispensable function of Prss35 in fertility, Plos One, 8, e56757,, 2013. 

Distl, O.: Mechanisms of regulation of litter size in pigs on the genome level, Reprod. Domest. Anim., 42, 10–16, 2007. 

Dube, B., Mulugeta, S. D., and Dzama, K.: Estimation of genetic and phenotypic parameters for sow productivity traits in South African Large White pigs, S. Afr. J. Anim. Sci., 42, 389–397, 2012. 

Franczak, A., Zglejc, K., Waszkiewicz, E., Wojciechowicz, B., Martyniak, M., Sobotka, W., Okrasa, S., and Kotwica, G.: Periconceptional undernutrition affects in utero methyltransferase expression and steroid hormone concentrations in uterine flushings and blood plasma during the peri-implantation period in domestic pigs, Reproduc. Fert. Develop., 29, 1499–1508, 2017. 

Hajkova, P., Erhardt, S., Lane, N., Haaf, T., El-Maarri, O., Reik, W., Walter, J., and Surani, M. A.: Epigenetic reprogramming in mouse primordial germ cells, Mech. Develop., 117, 15–23, 2002. 

Hao, Y., Cui, Y., and Gu, X.: Genome-wide DNA methylation profiles changes associated with constant heat stress in pigs as measured by bisulfite sequencing, Sci. Rep.-UK, 6, 27507,, 2016. 

Hedstrom, L.: Serine protease mechanism and specificity, Chemical reviews, 102, 4501–4524, 2002. 

Hu, Y., Xu, H., Li, Z., Zheng, X., Jia, X., Nie, Q., and Zhang, X.: Comparison of the genome-wide DNA methylation profiles between fast-growing and slow-growing broilers, Plos One, 8, e56411,, 2013. 

Hwang, J. H., An, S. M., Kwon, S., Park, D. H., Kim, T. W., Kang, D. G., Yu, G. E., Kim, I. S., Park, H. C., Ha, J., and Kim, C. W.: DNA methylation patterns and gene expression associated with litter size in Berkshire pig placenta, Plos One, 12, e0184539,, 2017. 

Jin, L., Jiang, Z., Xia, Y., Lou, P., Chen, L., Wang, H., Bai, L., Xie, Y., Liu, Y., Li, W., Zhong, B., Shen, J., Jiang, A., Zhu, L., Wang, J., Li, X., and Li, M.: Genome-wide DNA methylation changes in skeletal muscle between young and middle-aged pigs, BMC Genomics, 15, 653,, 2014. 

Klose, R. J. and Bird, A. P.: Genomic DNA methylation: the mark and its mediators, Trends Biochem. Sci., 31, 89–97, 2006. 

Ko, Y.-G., Cha, B.-H., Hwang, S., Im, G.-S., Yang, B.-C., Kim, M.-J., Cho, J.-H., and Seong, H.-H.: Altered DNA Methylation of Repetitive Sequences in Cloned Porcine Fetus, The Journal of Reproduction and Development Supplement, 578–578, 2008.  

Krueger, F. and Andrews, S. R.: Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications, Bioinformatics, 27, 1571–1572, 2011. 

Kurdyukov, S. and Bullock, M.: DNA Methylation Analysis: Choosing the Right Method, Biology, 5, 3,, 2016. 

Lande-Diner, L. and Cedar, H.: Silence of the genes–mechanisms of long-term repression, Nat. Rev. Genet., 6, 648–654, 2005. 

Laurent, L., Wong, E., Li, G., Huynh, T., Tsirigos, A., Ong, C. T., Low, H. M., Kin Sung, K. W., Rigoutsos, I., Loring, J., and Wei, C. L.: Dynamic changes in the human methylome during differentiation, Genome Res., 20, 320–331, 2010. 

Li, E., Beard, C., and Jaenisch, R.: Role for DNA methylation in genomic imprinting, Nature, 366, 362–365, 1993. 

Lister, R., Pelizzola, M., Dowen, R. H., Hawkins, R. D., Hon, G., Tonti-Filippini, J., Nery, J. R., Lee, L., Ye, Z., Ngo, Q. M., Edsall, L., Antosiewicz-Bourget, J., Stewart, R., Ruotti, V., Millar, A. H., Thomson, J. A., Ren, B., and Ecker, J. R.: Human DNA methylomes at base resolution show widespread epigenomic differences, Nature, 462, 315–322, 2009. 

Lorincz, M. C., Dickerson, D. R., Schmitt, M., and Groudine, M.: Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells, Nat. Struct. Mol. Biol., 11, 1068–1075, 2004. 

Mayer, W., Niveleau, A., Walter, J., Fundele, R., and Haaf, T.: Demethylation of the zygotic paternal genome, Nature, 403, 501–502, 2000. 

Meissner, A., Mikkelsen, T. S., Gu, H., Wernig, M., Hanna, J., Sivachenko, A., Zhang, X., Bernstein, B. E., Nusbaum, C., Jaffe, D. B., Gnirke, A., Jaenisch, R., and Lander, E. S.: Genome-scale DNA methylation maps of pluripotent and differentiated cells, Nature, 454, 766–770, 2008. 

Mesa, H., Safranski, T. J., Johnson, R. K., and Lamberson, W. R.: Correlated response in placental efficiency in swine selected for an index of components of lifter size, J. Anim. Sci., 81, 74–79, 2003. 

Messerschmidt, D. M., Knowles, B. B., and Solter, D.: DNA methylation dynamics during epigenetic reprogramming in the germline and preimplantation embryos, Gene. Dev., 28, 812–828, 2014. 

Mori, E., Baba, T., Iwamatsu, A., and Mori, T.: Purification and characterization of a 38-kDa protein, sp38, with zona pellucida-binding property from porcine epididymal sperm, Biochem. Biophys. Res. Co., 196, 196–202, 1993. 

Mori, E., Kashiwabara, S., Baba, T., Inagaki, Y., and Mori, T.: Amino acid sequences of porcine Sp38 and proacrosin required for binding to the zona pellucida, Dev. Biol., 168, 575–583, 1995. 

Newell-Price, J., Clark, A. J., and King, P.: DNA methylation and silencing of gene expression, Trends in endocrinology and metabolism: TEM, 11, 142–148, 2000. 

Nie, G., Li, Y., and Salamonsen, L. A.: Serine protease HtrA1 is developmentally regulated in trophoblast and uterine decidual cells during placental formation in the mouse, Dev. Dynam., 233, 1102–1109, 2005. 

Nielsen, B., Su, G., Lund, M. S., and Madsen, P.: Selection for increased number of piglets at d 5 after farrowing has increased litter size and reduced piglet mortality, J. Anim. Sci., 91, 2575–2582, 2013. 

Pistek, V. L., Fürst, R. W., Kliem, H., Bauersachs, S., Meyer, H. H. D., and Ulbrich, S. E.: HOXA10 mRNA expression and promoter DNA methylation in female pig offspring after in utero estradiol-17β exposure, J. Steroid Biochem., 138, 435–444, 2013. 

Pokkyla, R. M., Lakkakorpi, J. T., Nuojua-Huttunen, S. H., and Tapanainen, J. S.: Sequence variations in human ZP genes as potential modifiers of zona pellucida architecture, Fertil. Steril., 95, 2669–2672, 2011. 

Pollex, T. and Heard, E.: Recent advances in X-chromosome inactivation research, Curr. Opin. Cell Biol., 24, 825–832, 2012. 

Rutherford, K., Baxter, E., D'Eath, R., Turner, S., Arnott, G., Roehe, R., Ask, B., Sandøe, P., Moustsen, V., and Thorup, F.: The welfare implications of large litter size in the domestic pig I: biological factors, Anim. Welfare, 22, 199–218, 2013. 

Sasaki, H. and Matsui, Y.: Epigenetic events in mammalian germ-cell development: reprogramming and beyond, Nat. Rev. Genet., 9, 129–140, 2008. 

Schubeler, D.: Function and information content of DNA methylation, Nature, 517, 321–326, 2015. 

Singh, K., Molenaar, A. J., Swanson, K. M., Gudex, B., Arias, J. A., Erdman, R. A., and Stelwagen, K.: Epigenetics: a possible role in acute and transgenerational regulation of dairy cow milk production, Animal, 6, 375–381, 2012. 

Smith, Z. D. and Meissner, A.: DNA methylation: roles in mammalian development, Nat. Rev. Genet., 14, 204–220, 2013. 

Stevenson, T. J. and Prendergast, B. J.: Reversible DNA methylation regulates seasonal photoperiodic time measurement, P. Natl. Acad. Sci. USA, 110, 16651–16656, 2013. 

Su, J., Yan, H., Wei, Y., Liu, H., Liu, H., Wang, F., Lv, J., Wu, Q., and Zhang, Y.: CpG_MPs: identification of CpG methylation patterns of genomic regions from high-throughput bisulfite sequencing data, Nucleic Acids Res., 41, e4,, 2013. 

Varley, K. E., Gertz, J., Bowling, K. M., Parker, S. L., Reddy, T. E., Pauli-Behn, F., Cross, M. K., Williams, B. A., Stamatoyannopoulos, J. A., Crawford, G. E., Absher, D. M., Wold, B. J., and Myers, R. M.: Dynamic DNA methylation across diverse human cell lines and tissues, Genome Res., 23, 555–567, 2013. 

Wassarman, P. M.: Zona pellucida glycoproteins, J. Biol. Chem., 283, 24285–24289, 2008. 

Weber, M., Hellmann, I., Stadler, M. B., Ramos, L., Paabo, S., Rebhan, M., and Schubeler, D.: Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome, Nat. Genet., 39, 457–466, 2007. 

Wen, Y., Chen, F., Zhang, Q., Zhuang, Y., and Li, Z.: Detection of differentially methylated regions in whole genome bisulfite sequencing data using local Getis-Ord statistics, Bioinformatics, 32, 3396–3404, 2016. 

Wilson, M. E., Biensen, N. J., Youngs, C. R., and Ford, S. P.: Development of Meishan and Yorkshire littermate conceptuses in either a Meishan or Yorkshire uterine environment to day 90 of gestation and to term, Biol. Reprod., 58, 905–910, 1998. 

Xu, J., Bao, X., Peng, Z., Wang, L., Du, L., Niu, W., and Sun, Y.: Comprehensive analysis of genome-wide DNA methylation across human polycystic ovary syndrome ovary granulosa cell, Oncotarget, 7, 27899–27909, 2016. 

Yu, Y., Xu, W., Yi, Y. J., Sutovsky, P., and Oko, R.: The extracellular protein coat of the inner acrosomal membrane is involved in zona pellucida binding and penetration during fertilization: characterization of its most prominent polypeptide (IAM38), Dev. Biol., 290, 32–43, 2006. 

Yu, Y., Vanhorne, J., and Oko, R.: The origin and assembly of a zona pellucida binding protein, IAM38, during spermiogenesis, Microsc. Res. Techniq., 72, 558–565, 2009.  

Zhang, D., Cong, J., Shen, H., Wu, Q., and Wu, X.: Genome-wide identification of aberrantly methylated promoters in ovarian tissue of prenatally androgenized rats, Fertil. Steril., 102, 1458–1467, 2014. 

Short summary
Litter size is an important economic factor in pig production. This study was conducted to investigate litter-size-related epigenetic markers in uterine tissue from Berkshire pigs. The zona pellucida binding protein (ZPBP) gene was significantly hypomethylated in the larger litter size groups by DNA profiling analysis, and this gene can be used as a valuable epigenetic biomarker for hyperprolific sows.
Litter size is an important economic factor in pig production. This study was conducted to...