Skip to main content
  • Research article
  • Open access
  • Published:

Microarray profiling for differential gene expression in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows



The Chinese Taihu is one of the most prolific pig breeds in the world, which farrows at least five more piglets per litter than Western pig breeds partly due to a greater ovulation rate. Variation of ovulation rate maybe associated with the differences in the transcriptome of Chinese Taihu and Large White ovaries. In order to understand the molecular basis of the greater ovulation rate of Chinese Taihu sows, expression profiling experiments were conducted to identify differentially expressed genes in ovarian follicles at the preovulatory stage of a PMSG-hCG stimulated estrous cycle from 3 Chinese Taihu and 3 Large White cycling sows by using the Affymetrix Porcine Genechip™.


One hundred and thirty-three differentially expressed genes were identified between Chinese Taihu and Large White sows by using Affymetrix porcine GeneChip (p ≤ 0.05, Fold change ≥ 2 or ≤ 0.5). Gene Ontology (GO) analysis revealed that these genes belonged to the class of genes that participated in regulation of cellular process, regulation of biological process, biological regulation, developmental process, cell communication and signal transduction and so on. Significant differential expression of 6 genes including WNT10B and DKK2 in the WNT signaling pathway was detected. Real-time RT-PCR confirmed the expression pattern in seven of eight selected genes. A search of chromosomal location revealed that 92 differentially expressed transcripts located to the intervals of quantitative trait loci (QTLs) for reproduction traits. Furthermore, SNPs of two differentially expressed genes- BAX and BMPR1B were showed to be associated with litter size traits in Large White pigs and Chinese DIV line pigs (p ≤ 0.1 or p ≤ 0.05).


Our study detected many genes that showed differential expression between ovary follicles of two divergent breeds of pigs. Genes involved with regulation of cellular process, regulation of biological process, in addition to several genes not previously associated with ovarian physiology or with unknown function, were differentially expressed between two breeds. The suggestive or significant associations of BAX and BMPR1B gene with litter size indicated these genetic markers had the potentials to be used in pig industry after further validation of their genetic effects. Taken together, this study reveals many potential avenues of investigation for seeking new insights into ovarian physiology and the genetic control of reproduction.


Reproductive traits are of primary interest in livestock because they play a major role in efficiency of production. Selection for increased number of offspring has been employed in pigs with only limited success because of its low heritability and sex-limited nature [1]. Genetic characterization of litter size and its components (e.g. ovulation rate and embryo survival) will increase our understanding of the underlying physiology and could enhance genetic improvement through use of marker-assisted selection (MAS) [2]. In the past several decades, the reproductive strategy of the Chinese Taihu pigs, a breed which farrows three to five more piglets per litter than American or European pig breeds, has come under intense scrutiny [3, 4]. The greater litter size at farrowing in multiparous Chinese Taihu sows is due, in part, to a greater ovulation rate, a greater embryonic survival, a lower fertilization failure rate, and a greater uterine capacity [57].

In order to isolate the factors controlling the component traits of litter size, the differentially expressed (DE) genes were characterized during conceptus transformation, in Meishan-Landrace conceptuses and endometrial tissue when compared with conventional Landrace sows, in the porcine endometrium between pregnant and non-pregnant sows, and in the Erhualian and Large White placenta [811]. Furthermore, DE genes were identified in porcine ovarian follicles of multiparous sows on 12 d to 14 d of the estrous cycle between a line of pigs selected for an index of ovulation rate and embryo survival and its randomly selected control line [2, 12]. However, the DE genes in Chinese Taihu and Large White preovulatory follicles remain unexplored.

To develop a broader view of the gene expression in preovulatory ovary and to identify the key genes involved in ovulation, we used Affymetrix microarrays to screen the genes differentially expressed in preovulatory follicles from Chinese Taihu and Large White sows simulated by hCG at 80 h after PMSG administration. Bioinformatics analysis has revealed these DE genes were involved in important biological processes such as reproduction and the DE genes were then in silico mapped to quantitative trait loci (QTL) regions related to reproduction traits. Real-time RT-PCR was used to confirm the expression profiles of various genes. And association analyses of two DE genes (BAX and BMPR1B gene) with litter size were carried out to screen the molecular markers for litter size. This research identified candidate genes, molecular markers and pathways associated with ovulation rate, and gained a further insight into the genetic basis of the prolificacy of Chinese Taihu pigs.


Transcriptome analysis

Expression profiling experiments of Large White and Chinese Taihu ovary follicles were conducted by a commercial Affymetrix Porcine Genechip including 24,123 probe sets, which represent 23,999 transcripts and 124 controls. The transcriptome of ovarian follicles from Chinese Taihu sows was determined, and 23,921 probe sets were identified to have expression in the ovary follicles. Expression was detected for 16,066 transcripts (67.16% of all probe sets) in Large White ovarian follicles. A total of 15,477 transcripts (64.70% of all probe sets) were expressed in Chinese Taihu ovarian follicles (see additional file 1).

Differentially expressed (DE) gene analysis

The global expression profile of Chinese Taihu porcine ovaries challenged with PMSG-hCG was compared with that of the Large White sows. Probe sets whose intensities were normalized and filtered, then were subjected to identifying significantly differentially expressed genes using T-test analysis by comparing the log2 (normalized signal) of two breeds, and 1,476 transcripts were identified to be differentially expressed at the p ≤ 0.05 level (see additional file 2). Fold change (FC) is the gene expression level (normalized signal intensity) of Chinese Taihu sows compared to Large White compared to Large White sows. Taking a FC ≥ 2 or ≤ 0.5 and the p ≤ 0.05 significance level as the criteria, 133 transcripts showed differential expression. A set of 29 transcripts belonged to the up-regulated group and the other set of 104 transcripts belonged to the down-regulated group in Chinese Taihu sows (see additional file 3). Of the 133 DE transcripts, 12 were annotated pig genes, 83 transcripts could be determined for human homologous putative identities based on BLAST searches, and 38 additional probes did not have any database matches. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of 95 annotated DE gene lists were carried out by using the Database for Annotation, Visualization and Integrated Discovery (DAVID). GO annotation mapping revealed that the genes assigned GO terms for regulation of cellular process, regulation of biological process, biological regulation, developmental process, cell communication and signal transduction and so on (Figure 1). DE transcripts were involved in some important pathways including the WNT signaling pathway, p53 signaling pathway and Amyotrophic lateral sclerosis (ALS).

Figure 1
figure 1

Categories of 95 annotated transcripts in pigs or in human (FC ≥ 2 or ≤ 0.5, p ≤ 0.05) based on biological process GO term. Many categories shared the same transcripts.

Verification of gene expression pattern from microarray data by real-time RT-PCR

Eight genes (GALP, ITM2A, CTNNBIP1, WIF1, BMPR1B, CTNNAL1, LOC396850,CNN2) were selected to confirm the expression pattern by real-time RT-PCR. The results indicated that expression patterns of seven genes were consistent with the microarray (Pearson correlation coefficient ≥ 0.614, Figure 2). A contradiction expression profile was detected in CNN2 gene (Pearson correlation coefficient = - 0.766). This variance most likely came from technical differences.

Figure 2
figure 2

Validation of the microarray data by quantitative real-time RT-PCR analysis of eight representative genes. The x-axis represents the genes and the y-axis shows the relative expression (-ΔCt values for quantitative real-time RT-PCR; Log(Sample signal, 10) for microarray). The number of biological replicates is three for both assays. R represents the Pearson correlation coefficient. The significance of differences for the expression between Chinese Taihu (TH) and Large White (LW) PMSG-hCG stimulated ovary follicles was calculated using two-tailed T-test. p value: left for microarray, and right for quantitative real-time RT-PCR (qRT-PCR).

In silico mapping of DE genes to the interval of reproduction QTLs

Through BLAST of the 133 probes sets against Ensembl Sus_scrofa database (Sscrofa9, Apr 2009), 116 differentially expressed probes sets (transcripts) were in silico mapped pig chromosome covering 1-18 and X (see additional file 3). Ninety-two DE transcripts were located to porcine reproduction QTL chromosomal regions (see additional file 4). One such gene is the pig BMPR1B gene, which was localized to chromosome 8 within the interval of QTLs for ovulation rate, number of corpora lutea, litter size and prenatal survival [1317]. Another gene is Kit-ligand/stem cell factor gene, which was located in the QTLs region for ovulation rate and FSH level [18, 19].

BAX and BMPR1B as the candidate genes for litter size

BAX and BMPR1B were selected to be the candidate genes for litter size based on their biological functions on female reproduction and their positions within the intervals of QTLs. BAX T/C mutation in the intron 1 was genotyped as previously described [20], and BMPR1B C/G mutation was genotyped by PCR-Sac II-RFLP method which was established in the present study. As shown in Figure 3, the 416 bp BAX gene PCR amplicon containing T/C mutation was detectable by digestion with Ear I, resulting in allele T (416 bp) and allele C (296 and 120 bp) [20]; The 501 bp BMPR1B gene PCR amplicon containing C/G mutation was detectable by digestion with Sac II, resulting in allele C (501 bp) and allele G (321 and 180 bp). Genotyping results showed some variations in allele frequency between Chinese indigenous and Western breeds (Table 1). Briefly, the BAX allele C was fixed in Taihu, Tongcheng, Hezuo, Landrace, Duroc and White Duroc pigs; The BMPR1B allele C was fixed in Huainan, Tongcheng, Hezuo, Large White I, Landrace, Duroc, White Duroc and Pietrain pigs. Association analyses revealed that BAX CC had larger litter size than TC for both total number of piglets born (TNB) and the number of piglets born alive (NBA) of all parities in Large White pigs (p ≤ 0.05), and for NBA of both the 1st parity and all the parities in the DIV pigs (p ≤ 0.05) (see additional file 5). The additive effects of BAX gene in Large White pigs ranged from 0.41 to 1.00 (p ≤ 0.1). The BMPR1B CC was suggestively different with CG (p = 0.09) in NBA of all the parities in DIV line and the dominance effect was -0.66 piglet/litter (p = 0.06) (see additional file 5).

Figure 3
figure 3

Genotyping of single nucleotide polymorphisms in the porcine BAX and BMPR1B genes using PCR-RFLP. (A): BAX gene, Genotype CC: 296 bp + 120 bp; Genotype TT: 416 bp; Genotype TC: 416 bp + 296 bp + 120 bp. (B): BMPR1B gene, Genotype GG: 321 bp + 180 bp; Genotype CC: 501 bp; Genotype CG: 501 bp + 321 bp + 180 bp.

Table 1 Allele frequencies of BAX and BMPR1 B gene polymorphisms in 11 pig populations


It is reported that a proportionately greater increase in ovulation rate occurs in Chinese Taihu after they have experienced pregnancy than with European breeds [21]. Using microarray expression profiling to identify changes in RNA levels resulting from genetic background, our studies revealed differential expression of 133 DE genes in ovarian follicle tissues at the preovulatory phase of the estrous cycle. Several genes not previously associated with the processes of follicle maturation or ovulation, were differentially expressed in follicle pools of the Large White and Chinese Taihu sows. Moreover, 83 probes with BLASTN matches to gene sequences of known function in human and 38 probes with no database matches, were also identified as being differentially expressed.

Three DE gene pathways associated with p53 gene

WNT signaling pathway, p53 signaling pathway and Amyotrophic lateral sclerosis (ALS) were the major pathways including DE transcripts, and all of these three pathways are associated with p53 which was shown to promote apoptosis but not proliferation in ovarian granulosa cells [22, 23].

The first evidence of the important role of WNT signaling in the ovary was provided by Vainio et al., which demonstrated that WNT4 mutant mice lost most of their oocyte reserve in the days prior to birth [24]. In the present study, six genes (p53, WNT10B, TCF12, DKK2, CTNNBIP1, WIF1) involved in WNT signaling pathway were differentially expressed in Large White and Chinese Taihu sows. DKK2 acts as an antagonist of WNT signaling to induce endocytosis of the WNT-Fz receptor complex [25]. CTNNBIP1 was found to negatively regulate WNT signaling via inhibition of the interaction between β-catenin and T cell factor (TCF) [26]. Wnt inhibitory factor 1 (WIF1), a secreted antagonist that can bind to WNT proteins directly and inhibit WNT signaling pathway, was also differentially expressed (fold change > 2, p = 0.07).

Pathway analysis found that two DE genes (p53, BAX) were involved in p53 signaling pathway and Amyotrophic lateral sclerosis. The BCL-2 associated X protein (BAX) is a member of the BCL-2 protein family and functions as an accelerator of apoptosis. The regulatory activity of BAX in follicular endowment likely occurs during primordial germ cell (PGC) migration, prior to PGC colonization of the gonad [27].

Hormone related genes

Animal reproduction system is regulated by a suite of hormones. The serum hormone level of Chinese Taihu pigs was compared with the European in order to detect the mechanisms of high prolificacy of Chinese Taihu pigs. Chinese Taihu ovarian follicles tended to have a higher concentration of oestradiol in follicular fluid (p < 0.06) compared to Large White hybrids [28].

In our study, very limited hormone related genes were differentially expressed. One such gene encodes the caspase 8 associated protein 2 (CASP8AP2) which differentially suppresses steroid hormone receptor - induced transcriptional activity by interfering with their association with steroid hormone receptor coactivator 2 (SRC2)/N-CoA2 and SRC3/N-CoA3 [29]. Another one encodes cysteine - rich angiogenic inducer 61 (CYR61) which was known to be responsive to estrogen or to interact with estrogen receptor [30]. GALP gene encodes a member of the galanin family of neuropeptides and has been implicated in biological processes involving the central nervous system including hypothalamic regulation of metabolism and reproduction [31].

However, the well-known hormone related genes including estrogen receptor (ESR), follicle stimulating hormone (FSH-β) beta, retinol binding protein 4 (RBP4) genes had similar expressions in the two pig breeds. There are three major gene systems associated with follicular growth and ovulation: (1) IGF-1, estrogen and FSH gene, (2) the steroid receptors as well as orphan nuclear receptors and their ligands, and (3) the Wnt/Frizzled pathways [32]. In our study, their expressions of the genes from the first two systems probably were normalized by PMSG-hCG, and the genes from the last system remained differentially expressed.

Positional candidate genes for female reproduction traits

Differentially expressed genes were in silico mapped to reproduction QTL regions. One of these genes is BMPR1B which played a pivotal role in bone formation during embryogenesis and fracture repair, and its mutation was associated with the increased ovulation rate of the Booroola sheep with FecB phenotype [3335]. The pig BMPR1B gene was mapped to 108 cM on chromosome 8 [36], where the QTL for ovulation rate at approximately 105 cM, the QTL for number of corpora lutea at approximately 101 cM and 99 cM were reported to be located [1315].

Another gene, the KIT gene encoding the transmembrane tyrosine kinase receptor KIT maps to the mouse W locus [37]. Kit-ligand/stem cell factor induces primordial follicle development and initiates folliculogenesis [38]. The pig KIT gene was located to 8p12-21, which is the same region as the QTL for ovulation rate at approximately 5 cM and for FSH level at approximately 12.4-11.1 cM [18, 19, 39].

Two DE genes mutations and their associations with litter size

Marker-assisted selection (MAS) in conjunction with traditional selection methods is most effective for the traits such as litter size, which are either expressed later in life, are sex-dependent, or are of low heritability [40]. The candidate gene approach has led to notable success in demonstrating reproduction-related genetic markers or major genes, such as ES R, the erythopoietin receptor (EPOR), the prolactin receptor (PRLR), FSH-β and RBP4[1, 5, 4143]. Differentially expressed transcripts also can be evaluated as potential candidates if they meet the following criteria: the chromosomal location to a region containing a QTL, the importance of the tissue of interest in regulating the expression of the specific quantitative trait, and the critical developing stage responsible for the differences in the specific quantitative trait [44].

Here we selected the differentially expressed genes- BAX and BMPR1B as the candidate genes for litter traits due to their chromosomal location and/or their biological function. Some mutations were genotyped in several Chinese native pig breeds and Western pig breeds, and the results showed allele BAX C and BMPR1B C were fixed in several tested pig populations including Huainan, Tongcheng, Landrace, Duroc, White Duroc and Pietrain pigs. We also noticed that in DIV line and Large White pigs allele BAX C and BMPR1B C were not fixed, but BAX TT and BMPR1B GG did not nearly appear, possibly because the homozygotes led to embryonic lethality in these specific pig populations. The association analyses of these mutations with litter traits were carried out in Large White and Chinese DIV line pigs. The BAX gene was significantly associated with the litter size in these two populations, suggesting this mutation might be a good marker for pig selection and breeding. Suggestive association (p = 0.09) was found between the BMPR1B and litter size, in some ways similar to the results of Tomas et al., which also only found suggestive associations (nominal p < 0.1) between the BMPR1B haplotypes and the number of piglets born alive and weaned when data from the first parity of Iberian × Meishan F2 sows were considered [45].


We have reported differential gene expression in ovarian follicles at the preovulatory stage of a PMSG-hCG stimulated estrous cycle of two pig breeds with different ovulation rate. Some differentially expressed genes and key pathways related to biological and cellular regulation were identified. Furthermore, two differentially expressed genes- BAX and BMPR1B genes were genotyped and associated with the litter size traits in Large White and Chinese DIV pigs. The results of this study provide an opportunity to elucidate the genetic control of ovulation rate and improve our understanding of high prolificacy of Chinese Taihu pigs.


Population and experimental design, tissue collection

All animal procedures were performed according to protocols approved by the Biological Studies Animal Care and Use Committee of Hubei Province, PR China. Three multiparous Taihu cyclic sows (TH; ≥ 2 parities) from Jiangsu Changshu Xuqin Corporation and 3 multiparous Large White cyclic sows (LW; ≥ 2 parities) from Jinpin farm of Huazhong Agricultural University that exhibited normal estrous cycles were treated with 1000 IU PMSG and 500 IU hCG as previously described [46]. Given that the difference of prolificacy between Large White and Chinese Taihu pigs possibly due to the differences in the characteristics of follicle maturation [28], PMSG-hCG was used to initiate and synchronize the follicular phase. PMSG-hCG stimulated sows were killed and the ovaries at the preovulatory phase were removed immediately and kept in ice-cold physiological saline (0.9% wt/vol NaCl) during dissection of follicles. Only healthy follicles with a diameter of > 5 mm were isolated and snap-frozen in liquid nitrogen [2]. Total RNA was isolated from ovary follicles with Trizol (Invitrogen). RNA yield was quantified by spectrophotometric analysis using the convention that 1 absorbance unit at 260 nm equals 40 μg/mL RNA.

Animals from 11 different populations were used to investigate the allele frequency, including 36 Chinese Taihu pigs, 40 Tongcheng pigs, 7 Hezuo pigs, 26 Huainan pigs, 129 pigs from DIV line (the 4th Dam line of Chinese lean-type new lines), 108 Large White pigs I, 165 Large White pigs II, 16 Landrace pigs, 29 Duroc, 33 White Duroc and 34 Pietrain pigs.

The association analyses were conducted in Large White and DIV pig populations. Two hundred and seventy three Large White pigs contained 108, 165 pigs from strain I and II, which were raised in the farms owned by Huazhong Agricultural University (HAU) and Hubei Institute of Animal Science and Veterinary Medicine, respectively. Synthetic Line DIV was a result of cross of Landrace, Large White, Tongcheng or Taihu pigs, and was raised in pig farm owned by HAU. During the consecutive years (2005-2010), TNB and NBA of animals were recorded in 491 litters of DIV Line sows and in 481 litters of Large White II sows.

Microarray hybridizations and data analysis

The RNA labeling and hybridization were conducted by a commercial Affymetrix array service according to Technical Manual of GeneChip® Expression Analysis (CapitalBio Corporation, Beijing, China). Briefly, a total of 5 μg RNA was converted to double - stranded cDNA using the one-cycle cDNA Synthesis Kit (Affymetrix, Inc., Santa Clara, CA) and T7-Oligo (dT) Primer. In vitro transcription (IVT) of cRNA from cDNA was conducted using the MEGAscript ® T7 Kit (Ambion, Inc.). After the cleanup of the cDNA and cRNA using the Sample Cleanup Module (Affymetrix), the GeneChip IVT Labeling Kit (Affymetrix) was used for synthesis of Biotin-Labeled cRNA. cRNA quality and concentration was checked by UV spectrophotometer analyses and 2 μg cRNA was then checked by formaldehyde denaturing gel electrophoresis in 1.2% agarose gel. Subsequently, labeled cRNA was fractionated and hybridized with the GeneChip Porcine Genome Array according to the standard procedures provided by the manufacturer. Chips were washed and stained with a GeneChip Fluidics Station 450 (Affymetrix) using the standard fluidics protocol. The probe arrays were scanned using the Affymetrix ® GeneChip® Scanner 3000. Six microarrays were used in the experiment, corresponding to the RNAs from PMSG-hCG stimulated preovulatory ovarian follicles of three Taihu sows and three Large White sows. The Affymetrix GeneChip Porcine Genome Array probe set contained 11 pairs of perfect match (PM) and mismatch (MM) 25-mer probes. The probe-pair (PM-MM) data were used to detect the expression level of genes on the array (present call, marginal call, and absent call) by MAS 5.0 (Wilcoxon signed rank test). MAS 5.0 was used to normalize the signal values above to 500 on each array. An invariant set normalization procedure was performed to normalize the different arrays using DNA-chip analyzer (dChip). Probe sets whose intensity were above 50 on at least three arrays were remained and then subjected to identifying significantly differentially expressed genes. Identification of DE genes in Large White and Chinese Taihu was conducted using SAS (Ver.8.1, T-test) by comparing the log2 (normalized signal) of two breeds. The p value cutoff for DE genes was set at 0.05. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE23985 Genes with significant similarities to the transcripts in nr database based on BLAST, were selected for GO analyses (EASE threshold = 0.1, count threshold = 2) and pathway analysis at the DAVID Annotation summary results were obtained by inputting the gene list of interest by selecting OFFICIAL_GENE_SYMBOL as identifier.

Real-time RT-PCR to analyze gene expression profiles

Total RNAs were treated with DNase I (Ambion) and reverse transcribed by the M-MLV Reverse Transcriptase (Promega) according to the manufacturer's instructions. Quantitative real time PCR was performed on the iQ™5 Real Time PCR Detection System (Bio-Rad) using SYBR® Green Real-time PCR Master Mix (Toyobo Co., Ltd., Japan). Each 25 μL real-time RT-PCR reaction included 12.5 μL SYBR Green Real-time PCR Master Mix, 350-500 ng cDNA, 0.3 μM primers (see additional file 6). PCR conditions consisted of 1 cycle at 95°C for 3 min, followed by 45 cycles at 94°C for 20 sec, 58°C for 20 sec, and 72°C for 18 sec, with fluorescence acquisition at 74°C. cDNAs from three Taihu and three Large White ovary follicles were used as the template to detect the expression changes of the target genes. All PCRs were performed in triplicate and gene expression levels were quantified relatively to the expression of β-actin using Gene Expression Macro software (Bio-Rad, Richmond, CA, USA) by employing ΔCt value. Student's t-test was conducted to identify differentially expressed genes. Due to the negative relationship between Ct and expression level, an improved method of the previous report was used to compare the results of real time RT-PCR and microarray by plotting the -ΔCt values of real time RT-PCR versus the log of the microarray signal for each gene [47]. The Pearson correlation coefficient was calculated and used to estimate the correlation of real time RT-PCR results and microarray results.

Integration of differentially expressed genes in porcine QTL regions

To retrieve the corresponding pig gene data of one differentially expressed EST, the probe (EST) sequence was used as a query to search homology using BLASTN against Ensembl Sus_scrofa database (Sscrofa9, Apr 2009). QTLs for reproduction traits were obtained from the PigQTLdb The differentially expressed probe sets located in these QTL intervals were identified through comparison analysis [11]. The reproduction traits included TNB, NBA, fully formed piglets, number of stillborn, ovulation rate, age at puberty and uterine capacity. We listed all Affymetrix probe sets that were located in these QTL intervals.

Detection of SNPs of differentially expressed genes and association analyses

Two differentially expressed genes- BAX and BMPR1B were selected as potential candidates for litter size due to their functions and/or chromosomal positions, and further were investigated for their effects on litter traits. PCR-RFLPs were used to identify their mutations including a C/T at intron 1 of BAX gene [20], and a C/G of BMPR1B gene (Fan B, private communication). Briefly, PCR reactions contained 1.5 mM MgCl2, 150 μM dNTPs, 0.5 μM of each PCR primers (see additional file 6), 1 μL DNA, and 1 U Taq DNA polymerase in a 25 μL volume. PCR was run as follows: 94°C for 4 min followed by 35 cycles of 94°C for 40 sec, 60°C for 40 sec, 72°C for 35 sec, and final extension of 72°C for 10 min. For PCR-RFLP profile, 8.5 μL of PCR products were digested with 5 U restriction endonuclease (TaKaRa) for 4 h at 37°C, and then separated by electrophoresis on a 1-3% agarose gel (with ethidium bromide) in 1 × TAE buffer. The relationships between differentially expressed genes genotypes and litter size traits of Large White (n = 273) and DIV (n = 129) pigs were evaluated with the general linear model (GLM) procedure of SAS version 8.0. The model: Yijkl= μ + Pi+ Sj+ Fk+ Gl+ eijkl, where Yijklis the observation of the trait, μ is the least square means, Piis the effect of i th parity (i = 1, 2, 3, 4 (parity ≥ 4 )), Sjis the effect of j th season, Fkis the effect of k th farm (k = 1, 2), Glis the effect of l th genotype (l = 1-3) and eijklis the random residual [48].



differentially expressed


fold change


gene ontology


quantitative trait loci

Chinese DIV line:

the 4th dam line of Chinese lean-type new line


galanin-like peptide


integral membrane protein


catenin, beta interacting protein 1


Wnt inhibitory factor 1


bone morphogenetic protein receptor, type IB


catenin (cadherin-associated protein) alpha-like1


calponin 2


p53 protein


BCL-2 associated X protein


v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog


estrogen receptor


prolactin receptor


follicle-stimulating hormone β


erythopoietin receptor


retinol binding protein 4


T cell factor


marker-assisted selection


Kyoto Encyclopedia of Genes and Genomes


total number of piglets born


the number of piglets born alive.


  1. Rothschild MF, Jacobson C, Vaske A, Tuggle K, Wang L, Short H, Eckardt GR, Sasaki S, Vincent A, McLaren DG: The estrogen receptor locus is associated with a major gene influencing litter size in pigs. Proc Natl Acad Sci USA. 1996, 93: 201-205. 10.1073/pnas.93.1.201.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Gladney CD, Bertani GR, Johnson RK, Pomp D: Evaluation of gene expression in pigs selected for enhanced reproduction using differential display PCR and human microarrays: I. Ovarian follicles. J Anim Sci. 2004, 82: 17-31.

    CAS  PubMed  Google Scholar 

  3. Zhang WC, Wu JS, Rempel WE: Some performance characteristics of prolific breeds of pigs in China. Livest Prod Sci. 1983, 10: 59-68. 10.1016/0301-6226(83)90007-6.

    CAS  Google Scholar 

  4. Haley CS, Lee GJ: Genetic components of litter size in Meishan and Large White pigs and their crosses. In Proc 4th World Congr Genet Appl to Livest Prod: 23-27 July 1990; Edinburgh. Edited by: Hill WG, Thompson R, Woolliams JA. 1990, Organizing Committee, 4th World Congress on Genetics Applied to Livestock Production, XV: 458-461.

    Google Scholar 

  5. Vallet JL, Freking BA, Leymaster KA, Christenson RK: Allelic variation in the erythropoietin receptor gene is associated with uterine capacity and litter size in swine. Anim Genet. 2005, 36: 97-103. 10.1111/j.1365-2052.2005.01233.x.

    CAS  PubMed  Google Scholar 

  6. Bazer FW, Thatcher WW, Martinat-Botte F, Terqui M: Conceptus development in Large White and prolific Chinese Meishan pigs. J Reprod Fertil. 1988, 84: 37-42. 10.1530/jrf.0.0840037.

    CAS  PubMed  Google Scholar 

  7. Bolet G, Martinat-Botte F, Locatelli P, Gruand J, Terqui M, Berthelot F: Components of prolificacy in hyperprolific Large White sows compared with the Meishan and Large White breeds. Genet Sel Evol. 1986, 18: 333-342. 10.1186/1297-9686-18-3-333.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Ross JW, Ashworth MD, Hurst AG, Malayer JR, Geiser RD: Analysis and characterization of differential gene expression during rapid trophoblastic elongation in the pig using suppression subtractive hybridization. Reprod Biol Endocrinol. 2003, 1: 23-10.1186/1477-7827-1-23.

    PubMed  PubMed Central  Google Scholar 

  9. Vallée M, Beaudry D, Roberge C, Matte JJ, Blouin R, Palin MF: Isolation of differentially expressed genes in conceptuses and endometrial tissue of sows in early gestation. Biol Reprod. 2003, 69: 1697-1706.

    PubMed  Google Scholar 

  10. Østrup E, Bauersachs S, Blum H, Wolf E, Hyttel P: Differential endometrial gene expression in pregnant and nonpregnant sows. Biol Reprod. 2010, 83: 277-285.

    PubMed  Google Scholar 

  11. Zhou QY, Fang MD, Huang TH, Li CC, Yu M, Zhao SH: Detection of differentially expressed genes between Erhualian and Large White placentas on day 75 and 90 of gestation. BMC Genomics. 2009, 10: 337-10.1186/1471-2164-10-337.

    PubMed  PubMed Central  Google Scholar 

  12. Caetano AR, Johnson RK, Ford JJ, Pomp D: Microarray profiling for differential gene expression in ovaries and ovarian follicles of pigs selected for increased ovulation rate. Genetics. 2004, 168: 1529-1537. 10.1534/genetics.104.029595.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Rathje TA, Rohrer GA, Johnson RK: Evidence for quantitative trait loci affecting ovulation rate in pigs. J Anim Sci. 1997, 75: 1486-1494.

    CAS  PubMed  Google Scholar 

  14. Wilkie PJ, Paszek AA, Beattie CW, Alexander LJ, Wheeler MB, Schook LB: A genomic scan of porcine reproductive traits reveals possible quantitative trait loci (QTL) for number of corpora lutea. Mamm Genome. 1999, 10: 573-578. 10.1007/s003359901047.

    CAS  PubMed  Google Scholar 

  15. Braunschweig MH, Paszek AA, Weller JI, Yang D, Hawken RJ, Wheeler MB, Schook LB, Alexander LJ: Generation and exploration of a dense genetic map in a region of a QTL affecting corpora lutea in a Meishan × Yorkshire cross. Mamm Genome. 2001, 12: 719-729. 10.1007/s003350020040.

    CAS  PubMed  Google Scholar 

  16. King AH, Jiang Z, Gibson JP, Haley CS, Archibald AL: Mapping quantitative trait loci affecting female reproductive traits on porcine chromosome 8. Biol Reprod. 2003, 68: 2172-2179. 10.1095/biolreprod.102.012955.

    CAS  PubMed  Google Scholar 

  17. Campbell EM, Nonneman D, Rohrer GA: Fine mapping a quantitative trait locus affecting ovulation rate in swine on chromosome 8. J Anim Sci. 2003, 81: 1706-1714.

    CAS  PubMed  Google Scholar 

  18. Rohrer GA, Ford JJ, Wise TH, Vallet JL, Christenson RK: Identification of quantitative trait loci affecting female reproductive traits in a multi-generation Meishan × White composite swine population. J Anim Sci. 1999, 77: 1385-1391.

    CAS  PubMed  Google Scholar 

  19. Rohrer GA, Wise TH, Lunstra DD, Ford JJ: Identification of genomic regions controlling plasma FSH concentrations in Meishan-White Composite boars. Physiol Genomics. 2001, 6: 145-151.

    CAS  PubMed  Google Scholar 

  20. Laenoi W, Knorr C, Charoensook R, Gatphayak K, Vearaslip T, ter Meuleu U, Brenig B: Molecular characterization and chromosomal assignment of porcine BAX and Tac1 genes. In Proceedings of international research on food security, natural resource management and rural development: 11-13 October 2006; Bonn. Edited by: Folkard Asch, Mathias Becker, Andreas Deininger, Prathap Pugalenthi. 2006, University of Bonn, 310-

    Google Scholar 

  21. Wilmut I, Galvin JM, Day BN, Haley CS: How does the Meishan pig achieve its prolificacy [abstract]?. J Reprod Fertil. 1991, 7: 65-

    Google Scholar 

  22. Amsterdam A, Keren-Tal I, Aharoni D, Dantes A, Land-Bracha A, Rimon E, Sasson R, Hirsh L: Steroidogenesis and apoptosis in the mammalian ovary. Steroids. 2003, 68: 861-867. 10.1016/j.steroids.2003.09.003.

    CAS  PubMed  Google Scholar 

  23. Cherian-Shaw M, Das R, Vandevoort CA, Chaffin CL: Regulation of steroidogenesis by p53 in macaque granulosa cells and H295R human adrenocortical cells. Endocrinology. 2004, 145: 5734-5744. 10.1210/en.2004-0253.

    CAS  PubMed  Google Scholar 

  24. Vainio S, Heikkila M, Kispert A, Chin N, McMahon AP: Female development in mammals is regulated by Wnt-4 signalling. Nature. 1999, 397: 405-409. 10.1038/17068.

    CAS  PubMed  Google Scholar 

  25. Witte F, Dokas J, Neuendorf F, Mundlos S, Stricker S: Comprehensive expression analysis of all Wnt genes and their major secreted antagonists during mouse limb development and cartilage differentiation. Gene Expr Patterns. 2009, 9: 215-223. 10.1016/j.gep.2008.12.009.

    CAS  PubMed  Google Scholar 

  26. Katoh M: Molecular cloning and characterization of LZIC, a novel gene encoding ICAT homologous protein with leucine zipper domain. Int J Mol Med. 2001, 8: 611-615.

    CAS  PubMed  Google Scholar 

  27. Greenfeld CR, Pepling ME, Babus JK, Furth PA, Flaw JA: BAX regulates follicular endowment in mice. Reproduction. 2007, 133: 865-876. 10.1530/REP-06-0270.

    CAS  PubMed  Google Scholar 

  28. Biggs C, Tilton JE, Craigon JC, Foxcroft GR, Ashworth CA, Hunter MG: Comparison of follicular heterogeneity and ovarian characteristics in Meishan and Large-White hybrid pigs. J Reprod Fertil. 1993, 97: 263-269. 10.1530/jrf.0.0970263.

    CAS  PubMed  Google Scholar 

  29. Kino T, Ichijo T, Chrousos GP: FLASH interacts with p160 coactivator subtypes and differentially suppresses transcriptional activity of steroid hormone receptors. J Steroid Biochem Mol Biol. 2004, 92: 357-363. 10.1016/j.jsbmb.2004.09.003.

    CAS  PubMed  Google Scholar 

  30. Wang Z, Liu BC, Lin GT, Lin CS, Lue TF, Willingham E, Baskin LS: Up- regulation of estrogen responsive genes in hypospadias: microarray analysis. J Urol. 2007, 177: 1939-1946. 10.1016/j.juro.2007.01.014.

    CAS  PubMed  Google Scholar 

  31. Kauffman AS, Buenzle J, Fraley GS, Rissman EF: Effects of galanin-like peptide (GALP) on locomotion, reproduction, and body weight in female and male mice. Horm Behav. 2005, 48: 141-151. 10.1016/j.yhbeh.2005.01.010.

    CAS  PubMed  Google Scholar 

  32. Richards JS, Russell DL, Ochsner S, Hsieh M, Doyle KH, Falender AE, Lo YK, Sharma SC: Novel signaling pathways that control ovarian follicular development, ovulation, and luteinization. Recent Prog Horm Res. 2002, 57: 195-220. 10.1210/rp.57.1.195.

    CAS  PubMed  Google Scholar 

  33. Mulsant P, Lecerf F, Fabre S, Schibler L, Monget P, Lanneluc I, Pisselet C, Riquet J, Monniaux D, Callebaut I, Cribiu E, Thimonier J, Teyssier J, Bodin L, Cognie Y, Chitour N, Elsen JM: Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc Natl Acad Sci USA. 2001, 98: 5104-5109. 10.1073/pnas.091577598.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Souza CJH, MacDougall C, Campbell BK, McNeilly AS, Baird DT: The Booroola (FecB) phenotype is associated with a mutation in the bone morphogenetic receptor type 1 B (BMPR1B) gene. J Endocrinol. 2001, 169: R1-R6. 10.1677/joe.0.169R001.

    CAS  PubMed  Google Scholar 

  35. Wilson T, Wu XY, Juengel JL, Ross IK, Lumsden JM, Lord EA, Dodds KG, Walling GA, McEwan JC, O'Connell AR, McNatty KP, Montgomery GW: Highly prolific Booroola sheep have a mutation in the intracellular kinase domain of bone morphogenetic protein IB receptor (ALK-6) that is expressed in both oocytes and granulosa cells. Biol Reprod. 2001, 64: 1225-1235. 10.1095/biolreprod64.4.1225.

    CAS  PubMed  Google Scholar 

  36. Kim JG, Song JH, Vallet JL, Rohrer GA, Johnson GA, Joyce MM, Christenson RK: Molecular characterization and expression of porcine bone morphogenetic protein receptor-IB in the uterus of cyclic and pregnant gilts. Biol Reprod. 2003, 68: 735-743. 10.1095/biolreprod.102.008433.

    CAS  PubMed  Google Scholar 

  37. Chabot B, Stephenson DA, Chapman VM, Besmer P, Bernstein A: The proto-oncogene c-kit encoding a transmembrane tyrosine kinase receptor maps to the mouse W locus. Nature. 1988, 335: 88-89. 10.1038/335088a0.

    CAS  PubMed  Google Scholar 

  38. Parrott JA, Skinner MK: Kit-ligand/stem cell factor induces primordial follicle development and initiates folliculogenesis. Endocrinology. 1999, 140: 4262-4271. 10.1210/en.140.9.4262.

    CAS  PubMed  Google Scholar 

  39. Sakurai M, Zhou JH, Ohtaki M, Itoh T, Murakami Y, Yasue H: Assignment of c-KIT gene to swine chromosome 8p12-p21 by fluorescence in situ hybridization. Mamm Genome. 1996, 7 (5): 397-10.1007/BF03035316.

    CAS  PubMed  Google Scholar 

  40. Soller M: Marker assisted selection - an overview. Anim Biotech. 1994, 5: 193-207. 10.1080/10495399409525821.

    Google Scholar 

  41. Rothschild MF, Messer L, Day A, Wales R, Short TH, Southwood OI, Plastow G: Investigation of the retinol-binding protein 4 (RBP4) gene as a candidate gene for increased litter size in pigs. Mamm Genome. 2000, 11: 75-77. 10.1007/s003350010015.

    CAS  PubMed  Google Scholar 

  42. Li N, Zhao YF, Xiao L, Zeng FJ, Chen YZ, Dai RJ, Zang JS, Shen SQ, Chen YF, Wu CX: Candidate gene analysis for identification of genetic loci controlling litter size in swine. In Proc 6th World Congr Genet Appl Livest Prod, 11-16 January 1998; Armidale, Australia. 1998, World Congress on Genetics Applied to Livestock Production: Organizing Committee, 6th World Congress on Genetics Applied to Livestock Production, 26: 183-186.

    Google Scholar 

  43. Vincent A, Evans G, Short TH, Southwood OI, Plastow GS, Tuggle CK, Rothschild MF: The prolactin receptor gene is associated with increased litter size in pigs. In Proc 6th World Congr Genet Appl Livest Prod, 11-16 January 1998; Armidale, Australia. 1998, World Congress on Genetics Applied to Livestock Production: Organizing Committee, 6th World Congress on Genetics Applied to Livestock Production, 27: 15-18.

    Google Scholar 

  44. Lee SJ, Liu J, Qi N, Guarnera RA, Lee SY, Cicila GT: Use of a panel of congenic strains to evaluate differentially expressed genes as candidate genes for blood pressure quantitative trait loci. Hypertens Res. 2003, 26: 75-87. 10.1291/hypres.26.75.

    CAS  PubMed  Google Scholar 

  45. Tomas A, Frigo E, Casellas J, Ramırez O, Ovilo C, Noguera JL, Sanchez A: An association study between polymorphisms of the porcine bone morphogenetic protein receptor type 1b (BMPR1B) and reproductive performance of Iberian X Meishan F2 sows. Anim Genet. 2006, 37: 297-298. 10.1111/j.1365-2052.2006.01456.x.

    CAS  PubMed  Google Scholar 

  46. Torner H, Brussow KP, Alm H, Ratky J, Kanitz W: Morphology of porcine cumulus-oocyte-complexes depends on the stage of preovulatory maturation. Theriogenology. 1998, 50: 39-48. 10.1016/S0093-691X(98)00111-3.

    CAS  PubMed  Google Scholar 

  47. Luo L, Ye L, Liu G, Shao G, Zheng R, Ren Z, Zuo B, Xu D, Lei M, Jiang S, Deng C, Xiong Y, Li F: Microarray-based approach identifies differentially expressed microRNAs in porcine sexually immature and mature testes. PLoS ONE. 2010, 5: e11744-10.1371/journal.pone.0011744.

    PubMed  PubMed Central  Google Scholar 

  48. Zhang YH, Mei SQ, Peng XW, Niu BY, Ren ZQ, Zuo B, Xu DQ, Lei MG, Zheng R, Jiang SW, Deng CY, Xiong YZ, Li FE: Molecular characterization and SNPs analysis of the porcine Deleted in AZoospermia Like (pDAZL) gene. Anim Reprod Sci. 2009, 112: 415-422. 10.1016/j.anireprosci.2008.05.069.

    CAS  PubMed  Google Scholar 

Download references


We thank all farmers for providing pig samples, and thank Yuhao Zhang and Guochao Shao for collecting the samples. This work was supported financially by the Trans-gene Key Project of China (2009ZX08009-148B), IFS project (B/4316-1), Natural Science Foundation of Hubei Province of China (2009CDA127) and the creative team project of Chinese education ministry (IRT-0831).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Fenge Li.

Additional information

Authors' contributions

XS carried out most of bioinformatics analysis and lab works. SM participated in the animal samples collection and statistical analysis. HT, GW, LS made substantial contributions to lab works. SJ, CD, YX participated in the experiment design and coordination. FL conceived the study and drafted the manuscript. All authors read and approved the final manuscript.

Xiaojie Sun, Shuqi Mei contributed equally to this work.

Electronic supplementary material


Additional file 1:Ovarian follicle transcriptome analyses using the Affymetrix Porcine Genechip. Data of each probe is from the three Large White sows (1, 2, 3) and three Chinese Taihu (1, 2, 3) sows. Raw data from. CEL files were converted to gene signal files and raw signals were normalized above to 500 on each array by MAS 5.0, and the dChip invariant set was used to normalize the signals of the different arrays. "P", present; "A", absent; "M", marginal; "Count (P)", the number of P flag; select while Count (P) ≥2 in two groups (Chen et al., 2009)*. Totally, 16,066 and 15,477 probesets were detected expression in the Large White and Taihu, respectively. *Chen H, Li C, Fang M, Zhu M, Li X, Zhou R, Li K, Zhao S. Understanding Haemophilus parasuis infection in porcine spleen through a transcriptomics approach. BMC Genomics. 2009, 5, 10:64. (XLS 6 MB)


Additional file 2:1476 transcripts that are differentially expressed in porcine ovary between Large White sows and Chinese Taihu sows (p ≤ 0.05). Totally 105 (row 5-109) transcripts have been annotated based on BLASTX searches. "FC", Fold change, gene expression level (normalized signal intensity) of Chinese Taihu sows compared to Large White sows, ">1" represents up regulation, "<1" represents down regulation. "p-value", significance level of differential expression for a particular gene by comparing their log 2 (normalized signal). "Gene name", top informative BLASTX hit. (XLS 806 KB)


Additional file 3:Pig unique genes (row 4-15), unique genes with BLSATX identity but no annotation (row 16-98) and unique genes with no BLASTX identity (row 99-136) (p ≤ 0.05, FC ≥ 2 or ≤ 0.5). "FC", Fold change, gene expression level (normalized signal intensity) of Chinese Taihu sows compared to Large White sows, "≥2" represents up regulation, "≤0.5" represents down regulation. "p-value", significance level of differential expression for a particular gene by comparing their log 2 (normalized signal). "Gene name", top informative BLASTX hit. (XLS 272 KB)


Additional file 4:Genes that distributed in the intervals of the known pig QTLs for reproduction traits. All the genes listed here are from 133 transcripts that are significantly up/down-regulated between Large White and Taihu sows. GenBank IDs represent annotated genes. (XLS 460 KB)


Additional file 5:Association between BAX and BMPR1B genotype and litter size traits. A and B, p ≤ 0.05; a and b, p = 0.09. ** p ≤ 0.05, *p ≤ 0.1 = 0.06. N: Number of investigated litters. (DOC 73 KB)


Additional file 6:Primer pairs designed for genes selected for validation by real time RT-PCR and SNP detection. (XLS 28 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Cite this article

Sun, X., Mei, S., Tao, H. et al. Microarray profiling for differential gene expression in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows. BMC Genomics 12, 111 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: