Endometrial gene expression profiling in pregnant Meishan and Yorkshire pigs on day 12 of gestation
BMC Genomics volume 15, Article number: 156 (2014)
Litter size in pigs is a major factor affecting the profitability in the pig industry. The peri-implantation window in pigs is characterized by the coordinated interactions between the maternal uterine endometrium and the rapidly elongating conceptuses and represents a period of time during which a large percentage of the developing conceptuses are lost. However, the gene expression and regulatory networks in the endometrium contributing to the establishment of the maternal: placental interface remain poorly understood.
We characterized the endometrial gene expression profile during the peri-implantation stage of development by comparing two breeds that demonstrate very different reproductive efficiencies. We employed the porcine Affymetrix GeneChip® to assay the transcriptomic profiles of genes expressed in the uterine endometrium obtained from Meishan and Yorkshire gilts (n = 4 for each breed) on day 12 of gestation (M12 and Y12, respectively). Total of 17,076 probesets were identified as "present" in at least two arrays. A mixed model-based statistical analysis predicted a total of 2,656 (q < 0.1) transcripts as differentially expressed between Meishan and Yorkshire pigs. Eighteen differentially expressed transcripts of interest were validated by quantitative real-time PCR. Gene ontology (GO) annotation revealed that the known functions of the differentially expressed genes were involved in a series of important biological processes relevant to early pregnancy establishment in the pig.
The results identified endometrial gene expression profiles of two breeds differing in litter size and identified candidate genes that are related to known physiological pathways related to reproductive prolificacy. These findings provide a deeper understanding of molecular pathways differing between two breeds at the critical peri-implantation stage of pregnancy, which can be utilized to better understand the events contributing to pregnancy establishment in the pig.
The pig (Susscrofa) has been used as model animal in many fields, including biological , agricultural, and biomedical studies [2–4]. Pig litter size, which can be impacted by multiple factors including ovulation rate, embryonic survival, uterine capacity and fetal survival, is one of the most important economic traits affecting production efficiency in the pig industry . Improvements in litter size across the swine industry in recent decades has occurred through different selection schemes such as phenotypic selection, family index selection, BLUP-based selection, component trait selection, and hyperprolific selection [5, 6]. However, despite current improvements, the low heritability of litter size has limited the genetic selection progress in substantially improving this trait . Improvement of litter size is ideally suited for application of molecular breeding approaches such as marker-assisted selection (MAS) and whole-genome selection, where molecular markers associated with superior trait values may accelerate genetic improvement. In addition, numerous studies have attempted to detect candidate genes underlying the molecular mechanism of sow prolificacy. Although some candidate genes for sow prolificacy have been reported , true causal genes remain scant. A deeper understanding of the physiological and molecular mechanisms associated with reproductive efficiency is needed to further the biological underpinnings that contribute to controlling litter size in pigs.
It is now well established that prenatal mortality can be a major determinant of litter size in pigs [5, 9, 10]. When early conceptuses undergo rapid differentiation and expansion of their trophoblastic membranes during the first few weeks of gestation, abnormalities of embryogenesis that produce embryonic mortality can occur . The extent and timing of embryonic and fetal losses during gestation has been well documented by previous studies, and it is estimated that about 25-45% of fertilized ova do not survive through gestation . The majority of embryonic loss has been demonstrated to occur before day 20 of gestation, primarily around day 11–12 . Early events in conceptus development, including rapid trophoblastic elongation and successful establishment of the fetal: maternal interface, play a key part in determining litter size. The conceptus itself undergoes a dramatic transcriptomic reorganization during trophoblastic elongation and during initial interactions with the uterine endometrium . Additionally, the copious amounts of estrogen synthesized and released by the elongating conceptus serves as the maternal recognition of pregnancy signal around day 12 of gestation . The timing of endometrial exposure to estrogen is critical as specific molecular cues associated with uterine receptivity occur simultaneously with conceptus estrogen synthesis and abnormal exposure to estrogen prior to day 12 results in compromised pregnancy establishment [13–15].
Studies have revealed that there are very different patterns of gene expression in uterine endometrium between successful and abortive embryo implantation, as well as between high and low litter size breeds [16, 17]. It is now clear that the molecular events happening in the maternal endometrium during the early stages play a crucial role in determining the success of embryonic survival, and the early-stage endometrium provides an ideal window to investigate the molecular mechanisms that may be associated early conceptus loss and sow reproductive prolificacy. Many genes involved in angiogenesis and immune responses are expressed at the endometrium-fetal surface . Interestingly, polymorphisms in genes expressed during embryonic implantation, such as in EphA4 and TGF-β, have been associated with litter size [19, 20]. However, what is lacking is a systematic comparison investigating global gene expression patterns in the uterine endometrium during the peri-implantation stage of pregnancy between pig breeds with established differences in litter size.
Microarrays have been widely used to detect differentially expressed genes hypothesized to contribute to controlling litter size. Some interesting genes, such as steroidogenic acute regulatory protein and interleukin-1beta (IL1B), were found to be expressed in during early development using multiple global expression profiling approaches [21–24]. Differentially expressed genes have been detected in endometrium between pregnant and non-pregnant sows  and between endometrium of normal and developmentally delayed fetuses . Furthermore, genes involved in several important pathways regulating fetal development have been identified comparing gestational day 75 and 90 placentae from Erhualian and Large White pig breeds .
Differences in reproductive efficiency exist between different pig breeds. The prolific Meishan pigs have higher embryo survival rates during early gestation,which may be due to the fact that the Meishan conceptus exhibits a reduced trophectoderm mitotic rate during the preimplantation period, elongates from fewer cells and remains smaller throughout gestation compared with other low litter size breeds . To further the comparison between Meishan and domestic swine females, reciprocal embryo transfers were performed and revealed that litter size is controlled primarily by the recipient females rather than the sire or embryos [29–31]. These data suggest that developing a more comprehensive understanding of the maternal uterine microenvironment influencing embryonic survival and implantation during the establishment of pregnancy could be of significance towards understanding maternal mechanisms contributing to successful pregnancy establishment.
In the present investigation we conducted a genome-wide screen to characterize gene expression of endometrial tissue on day 12 of gestation, during maternal recognition of pregnancy. We used the Affymetrix GeneChip® to analyze endometrial gene expression in Meishan and Yorkshire, which have phylogenetically distinct genomes and also have striking prolificacy differences with respect to litter size. We analyzed these endometrial transcriptomic data through a comparative annotation approach to gain a deeper molecular insight into genes controlling the prolific nature of pigs. A number of genes and transcripts were detected to be differentially expressed in endometrial tissue between breeds. The data and our analysis in this investigation contribute toward a deeper understanding of the molecular events associated with reproductive success during early pregnancy establishment by comparing the molecular events between two separate breeds with different levels of prolificacy.
Transcriptome analysis in early pregnancy endometrium
The Affymetrix Porcine GeneChip® was used for transcriptomic analyses of endometrium from Yorkshire and Meishan pigs early in pregnancy with four biological replicates for each pig breed. To compare similar stages of endometrial development during gestation, we collected tissue from pregnant Meishan and Yorkshire gilts bred at their second estrous cycle to boars of the same breed. At the time of endometrial tissue collection, we determined the morphological stage conceptuses present with the uteri. We only processed endometrial tissue for this comparison when all embryos were at the filamentous stage. If the signal of a transcript from at least two of the four biological replicates were called “present” by the Affymetrix algorithm, we considered this transcript expressed in the endometrium. Using this criterion, 14,951 transcripts were identified in Meishan day 12 (M12) endometrium, and 14,558 transcripts were detected in Yorkshire day 12 (Y12) endometrium.
Identification of differentially expressed genes
Following statistical analysis, we identified 2656 (q < 0.1) transcripts that were differentially expressed between M12 and Y12 (Table 1). Among these transcripts, 2264 were annotated as known genes; thus about 15% of the transcripts were unable to be annotated. Among the differentially expressed genes, 1180 were significantly more abundant in M12 endometrium, while 1084 genes were more abundant in Y12.
Cluster analysis and GO annotations
To gain insight into similarities between the transcriptome of the two breeds, data from all the differentially expressed genes in the endometrium were used in a systematic cluster analysis (Figure 1). Although there is substantial variation in expression within the replicates of one breed, the heat map suggests clearly that the selected differential gene set associates the biological samples correctly into two groups each representing one breed. The quality of this association enabled further analysis to identify potential pathways and functions impacted in endometrium using the genes detected at statistically different abundance between these two breeds.
The Database for Annotation, Visualization and Integrated Discovery (DAVID 6.7) was utilized to annotate biological terms and processes preferentially represented in Meishan and Yorkshire endometrium during the peri-implantation stage of development. We used the gene lists of q < 0.1 for DAVID annotation, and tested using the Annotation Chart function in DAVID separately for enrichment of genes with greater abundance in Yorkshire endometrium compared to Meishan endometrium, as well as those with greater abundance in Meishan endometrium compared to Yorkshire endometrium. Terms enriched utilizing these gene lists are provided as Additional file 1: Table S1 and Additional file 2: Table S2.
The five most significant terms in each of the most enriched in the annotation clusters expressed at a higher level in Yorkshire tissue were membrane-enclosed lumen, RNA-splicing, RNA binding, nucleus, and blood vessel development (Table 2). The functions of these enriched genes are closely associated with the nucleus, transcriptional and RNA processing functions in cells. Other functions represented in additional significantly enriched clusters in this list of differentially expressed (DE) genes include apoptosis, cell-cell tight junction, and other functions (Additional file 2: Table S2A). It may be that the Yorkshire endometrium is more transcriptionally responsive to during maternal: fetal crosstalk following exposure to the maternal recognition of pregnancy signals. It is also interesting that genes associated with angiogenesis appear to be more prominently expressed in Yorkshire endometrium. These genes include transcription factors, such as JUN, Foxo1, and angiogenesis genes, such as VEGFC, KDR.
Alternatively, the five most significant gene ontology terms enriched in the genes expressed more abundantly in Meishan endometrium were translational elongation, ribosome biogenesis, protein folding, CHORD-containing proteins and SGT1(CS) domain, and intracellular organelle lumen (Table 3). In Meishan endometrium, enriched functions seem to focus on post-transcriptional processes, especially protein synthesis and maturation (Additional file 2: Table S2b). Genes such as ribosome proteins (RPL10, RPL11, RPL12, RPL14, RPL17, RPL18, RPL 35A, RPL36A, RPS2, RPS3, MRPS36) and translation elongation factors (EEF1A1, EEF1A2, EEF1G, EEF2) are more highly expressed in Meishan endometrium.
Validations by qRT-PCR
To validate the microarray results, we chose 18 differentially expressed genes identified following the microarray analysis which have also been reported in the literature as relevant to the physiology of implantation. The qRT-PCR results are shown in Figure 2. Of the 18 genes assayed, 14 confirmed statistically significant differential expression in the same direction as the microarray data, and the remaining four genes showed the same direction of expression between microarray and qRT-PCR analyses despite not being statistically different.
Regulatory network analysis
Identification of small molecules or proteins with a known regulatory link in the biological literature to multiple highly differentially expressed genes is of interest. To identify such sub-networks, a Sub-Network Enrichment Analysis (SNEA) was performed on the gene lists representing those greater in the Yorkshire breed and those greater in the Meishan breed. In all, 37 proteins and 10 small molecules were identified following SNEA with a p-value less than 0.05 (Additional file 3: Table S3). The most significantly enriched protein hub was regulated by RORA, retinoic acid receptor-related orphan receptor alpha (RORA; also known as nuclear receptor subfamily 2, group F, member 2 (NR2F2)) being linked to five genes that were significantly more abundant in endometrium on day 12 from Yorkshire compared to Meishan (Additional file 3: Table S3; Figure 3A). Another retinoic acid receptor, NR2F1, was also highly ranked on the protein list in Yorkshire, due to a link with seven differentially expressed genes (Additional file 3: Table S3; Figure 3A).
Using the list representing DE genes with greater abundance in day 12 endometrium of Meishan, the SNEA identified 59 proteins and 13 small molecules (P < 0.05) as potential regulatory hubs (Additional file 3: Table S3).Of interest, leptin was linked with 38 genes involved in multiple pathways, such as solute carrier family members (SLC25A25, SLC9A1), growth factors (SST, MSTN), and immune system (MYD88, IL12B, CD63); all genes with greater abundance in Meishan endometrium. Also of interest, GHRL was linked with 13 genes more abundantly expressed in Meishan (Additional file 3: Table S3; Figure 3B).
Identifying the endometrial mechanisms and pathways contributing to reproductive success in pigs has long been a research interest [8, 32]. The majority of prenatal mortality in pigs occurs before gestation day 20, particularly during the conceptus trophoblastic elongation and attachment to the uterine endometrium, beginning around day 12 of gestation . Establishing the interface between the developing conceptus and the maternal endometrium involves a number of complex signaling networks. It is well known that Meishan pigs farrow 3–5 more piglets than Yorkshire . A study focused on physiological differences between Meishan and Yorkshire embryos showed that in vitro cultured Meishan embryos had lower mitotic rate and fewer cells in day 12 conceptuses in vivo. Morphological comparison studies revealed that variation of conceptuses diameter was smaller for Meishan pigs, indicating that Meishan conceptuses potentially develop more uniformly between days 8 and 14 of gestation [35, 36]. Meishan and Yorkshire embryos transferred to Meishan uterine environment developed slower and secreted less estrogen, compared to being transferred to a Yorkshire recipient, suggesting that the recipient’s genotype was the major factor controlling embryonic growth and estrogen secretion . Comparison of gene expression during physiologically important events across different breeds could enable the identification of novel pathways contributing to reproductive success in pigs. This approach has previously been utilized with the Meishan breed characterizing the ovarian gene expression networks of F2 lberian x Meishan sows .
The advantage of microarray analysis is the simultaneous measurement of the expression patterns of large number of genes , thus enabling classification of molecular characteristics of phenotypes and generating hypotheses predicting phenotypic outcomes . In the pig, a number of high throughput platforms have been established for global expression profiling analysis [11, 15, 40–42]. In this study, we conducted a genome-wide transcriptomic analysis on uterine endometrium collected on day 12 of gestation to identify and better understand the molecular basis of differences in conceptus survival between Yorkshire and Meishan pigs. Collectively, we identified the presence of approximately 17,000 expressed transcripts in the uterine endometrium on day 12 of gestation and following statistic alanalysis identified 2656 transcripts with differentially abundance of which 1031 demonstrated 2-fold or greater difference. Following validation of the microarray results with a sub-set of 16 genes assayed with qRT-PCR, we utilized these lists to identify pathways and gene ontology terms associated with each breed by using those gene lists representing expression greater in one breed over another.
Interpretation of gene annotation analysis: selected pathways of interest
Hormone genes and transcriptional activity
Oxytocin (OXT) decreases uterine receptivity in rodents, as rats treated with the OXT antagonist, atosiban, before embryo transfer had improved implantation rates [43, 44]. It is of interest that we observed that OXT abundance was greater in Yorkshire compared to Meishan as assessed by both microarray and qRT-PCR. Given that mRNA expression is greater in Yorkshire, it is possible that OXT mediated biological mechanisms may exist that contribute to and/or impact conceptus attachment and interaction with the uterine endometrium.
Prostaglandin production is an important component contributing to the successful establishment of pregnancy in the pig. Expression of prostaglandin endoperoxide synthase 1 (PTGS1) and PTGS2 are substantially increased during implantation . Our microarray results showed that PTGS1 was more abundant in Meishan endometrium (1.7-fold difference, q value = 0.27) whereas PTGS2 was more abundant in Yorkshire endometrium (2.6-fold difference, q value = 0.15), which is consistent with the results in Erhualian endometrium compared to large white . These results indicated that higher PTGS1 may have an important role in embryo survival and attachment in Chinese breeds.
As shown in Additional file 2: Table S2, 157 genes expressed more abundantly in Yorkshire were clustered to the gene ontology (GO) term membrane enclosed lumen, 52 genes to GO term RNA processing, and 40 genes to GO term transcription factor binding pathway following DAVID analysis. Transcription factor transcripts more abundant in Yorkshire, such as HOXA10 and FOXO1, are known to be regulated by steroid hormones, including estrogen, progesterone and testosterone [47–49]. The potential that increased abundance of steroid-hormone dependent factors may be related to the greater quantity of steroid hormone secretion from Yorkshire embryos [34, 50].
Blood vessel development
Successful apposition and attachment of pig conceptuses followed by the early formation of a placental interface requires a maternal vascular support. We observed that genes detected at a greater abundance in Yorkshire pigs are enriched for GO terms related to blood vessel development. For example, cadherin 5, type 2 vascular endothelium (CDH5) is not only an essential adhesion molecule, but also contributes to angiogenesis which regulates cell proliferation and modulates VEGF receptors . The genes with greater abundance in Meishan, such as PHKA2, PRKAG3, PPP1CA, PYGL, PHKG1, CPS1 were clustered by DAVID analyses into the pathways of glycogen metabolism and glucan metabolism process.
Cell adhesion and tight junction genes
Cell-to-cell communication through connexin protein function has been shown in porcine endometrium to be lower than in species with invasive implantation . Thus endometrial cell-to-cell interaction may be involved in limiting trophoblast invasiveness, as porcine implantation is superficial, or in developing specific channels for this implantation type. Proteins such as matrix metalloproteinase (MMP) are endopeptidases that play an important role in the transient invasive property of trophoblastic cells during implantation. Expression of MMP7 was lower in human endometrium in women with repeated implantation failure compared to fertile women, indicating a potential role in human embryonic implantation, an invasive process . In pigs, our qRT-PCR and microarray results show that MMP7 expression is greater in Yorkshire endometrium on day 12 of pregnancy compared to Meishan. Integrin beta 3 (ITGB3) protein has also been detected at porcine embryonic attachment sites on day 24 . ITGB3 has been reported to be involved with signal transduction and formation of cytoplasmic focal adhesions, which was interpreted as supportive for implantation . ITGB3 transcript was in greater abundance in Yorkshire endometrium than Meishan. We also observed other genes expressed more abundantly in Yorkshire compared to Meishan that are annotated as members of tight junction pathways, such as Claudin (CLDN) family members (CLDN1, CLDN5, CLDN8, CLDN10, CLDN11) and tight junction protein 2(TJP2). Tight junction molecules, such as CLDN1, improved invasion activity in cancer cell lines through activation of MT1-MMP-1 and MMP-2 , and peak expression of claudin-4 RNA has been temporally associated with the implantation window in humans . These data may suggest that development of membrane channels for specific ions, a major function of claudins, may be important to facilitate successful conceptus attachment in Yorkshire pigs.
Immune related pathways
We also investigated the involvement of immune related pathways in DE genes between these breeds based on the importance of immunity in early maternal: fetal interactions. Successful pregnancy requires appropriate spatio-temporal regulation of genes contributing to the immune response during the peri-implantation period. Although reports exist that aborted placentas expressed less Th2 cytokines and greater Th1 cytokines than successful pregnancies [56, 57], recent literature indicated that uterine receptivity is actually promoted by a Th1 inflammatory response, whereas the Th2-humoral response facilitated pregnancy maintenance . We observed that genes expressed more abundantly in Meishan endometrium compared to Yorkshire are annotated as functioning in immune response pathways, such as CD28, which may indicate a positive role of CD28 during implantation. CD28 expression has been significant positively correlated with Th1 cytokines expression (IL2 and INFγ) . While reports on mice indicated that suppression of CD28 and other cytokines by cyclosporine A improved pregnancy rate , the exact involvement of CD28 in abortion was not clear. In this context, it is interesting that three peptidyl-prolylisomerases (PPIA, PPID and PPIG) are more abundantly expressed in Meishan (Additional file 2: Table S2B). Proteins in the PPI family are known to bind and potentially mediate cyclosporine A function, which in turn has been shown to improve pregnancy rate when administered at the implantation window stage in mice . Thus the PPI proteins may bind to natural cyclosporine analogs impacting the immunological regulation of conceptus attachment.
While not sufficiently enriched at the level of immune annotation clustering, there were also other immune-related genes that are more highly expressed in Yorkshire in the DAVID analysis. Among these, MAPS2 may function by activating complement component factors  and alpha-2-macroglobin (A2M) is involved in immune cell migration in pregnancy . Additionally, transferrin (TF) may be involved in removal of allergens from serum  and it is expressed differently in the pregnant pig compared the non-pregnant pig , indicating the potential role during pregnancy. Higher expression of immunity associated genes such as MAPS2, A2M and TF in Meishan compared to Yorkshire, but not the major cytokines, such as IL2 and IL10, may be suggestive that the difference between these two breeds was not as large as that between females with normal versus aborted pregnancies. However, the SNEA does indicate that genes expressed more abundantly in Meishan are linked to cytokines IL4 and IL6 (Additional file 3: Table S3).
Regulatory network analysis
RORA (or NR2F2) is a transcription factor belonging to the steroid hormone receptor superfamily  and has been demonstrated to be expressed in human endometrium . This transcription factor, as well as NR2F1, has been suggested to play an important role in mediating the cellular response to retinoic acid . RORA is shown to activate hypoxia-inducible factor 1 alpha (HIF-1α)  and HIF-1 is an important mediator of estrogen-induced VEGF expression in the uterus . The function of these networks is concordant with the conclusion from our DAVID analysis suggesting greater activation of angiogenesis at this stage in the Yorkshire breed.
Leptin is a small peptide synthesized by white adipose tissue that regulates body weight homeostasis. Leptin transcript abundance was significantly higher in the endometrium of fertile women as compared to women with recurrent implantation failure  and thus higher Meishan expression of Leptin in our study is consistent with the prolificacy seen in this breed. Ghrelin is thought as a novel paracrine/autocrine factor expressed in endometrium and placenta that is involved in crosstalk between endometrium and embryo during embryo implantation , which may indicate that such intra-tissue communication is critical for implantation success.
Conceptus and fetal loss during early pregnancy establishment can be a critical determinant of litter size. However, identification of molecular pathways that may contribute to successful pregnancy establishment and the formation of a fetal: maternal interface has not yet been investigated in breeds with established differences in litter size performance. In this study, the Affymetrix GeneChip® was used to compare gene expression profiles and identify associated pathways which may be differentially regulated in the uterine endometrium of Meishan and Yorkshire gilts during day 12 of gestation. The results have described the molecular portraits of uterine endometrium during the opening of the implantation window revealing differences between the Meishan and Yorkshire breeds. Our results revealed several classifications of genes participating in regulation of endometrium function during early pregnancy. Future analysis will require an exhaustive look into the protein repertoire of the uterine endometrium of the two breeds to strengthen our understanding of the relationship between the transcriptome and proteome in the uterine endometrium.
All animal use was performed under Iowa State University Committee on Animal Care approved protocols. Heat-check positive Yorkshire and Meishan gilts were bred twice to boars of the respective breed, 24 h apart, via natural service at their second or third estrus. Reproductive tracts were recovered immediately after slaughter at ISU using IACUC approved protocols. The gestation date was calculated from the first mating day. Endometrium was collected from 10 cm sample near junction with other uterine horn at gestation day 12. Following removal of the reproductive tract, the uterus was opened by cutting down the anti-mesometrial border and uterine endometrium was collected from all areas of the lumen and snap-frozen in liquid nitrogen. If no elongated conceptuses were observed in a bred female, tissue samples from that female were not used in the study. Uterine endometrial samples from a total of four Yorkshire and four Meishan gilts were chosen using these criteria.
Total RNA was isolated from 200 mg of tissue using the QIAGEN RNeasy kit (Cat no.714104). The DNA was removed by DNase I digestion and RNeasy mini elute kit cleanup as recommended by QIAGEN. PCR assay without reverse transcription was used to confirm that the RNA samples were absent for genomic DNA. The quantity and quality of the RNA were determined using Agilent 2100 Bioanalyzer(Agilent Technologies, Santa Clara, CA) and Nanodrop 2000 (Thermo Scientific, Wilmington, DE). RNA samples with RIN number lower than 7 or yield less than 3 μg was identified as low quality and excluded from the experiment.
Microarray hybridization and statistical analyses
The porcine genome microarrays were purchased from Affymetrix (Cat. no. 900623; Santa Clara, CA). The Affymetrix GeneChip® Porcine Genome Array contains 23,937 probesets which represents 20,201 genes and 190 controls. The sequence information for this array was selected from public data sources including UniGene Build 28, GenBank® mRNAs. The RNA labeling, porcine gene chip hybridization, washing and signal detection were done at the GeneChip Facility, Iowa State University, Ames, IA according to the manufacturer’s instructions.
Data from .cel files were converted to gene signal files by GCOS1.2. Natural log was taken for all the raw gene signals for normalization. The SAS ANOVA analysis was used to identify differential gene expression (SAS Institute, Cary, NC). A mixed model including random effects for slide and animal was fit to the expression data for each gene using GLM procedure in SAS. The p-values were obtained for the Meishan and Yorkshire comparison. The p-values for each test were converted to q-values for false discovery rate estimation .
RNA samples from the eight animals which were used for the Affymetrix experiment were analyzed by quantitative reverse-transcription PCR. The reverse transcription was performed using SuperScript II Reverse Transcriptase and Oligo(dT) primer according to the manufacturer’s instructions (Invitrogen, Carlsbad, CA). All reverse transcription reactions were run along with no-template controls. The no-template controls gave non-detectable signals in all samples, confirming the high specificity of the assays. Real-time PCR was performed using Stratagene Brilliant Kit (La Jolla,CA). The sequences for probes and primers have been previously published . All reactions were run in duplicate and data were normalized using the house keeping gene RPL32. The significance level was set to 0.05.
Differentially expressed gene list analysis
An unsupervised learning procedure in Gene Cluster 2 software was used for the cluster analysis and used normalized LSMean values of differentially expressed genes. Hierarchical cluster analysis was done by using Cluster and TreeView software (http://rana.lbl.gov/EisenSoftware.htm). The p-value cut off used for differential expressed genes in the cluster analysis is p < 0.005 (q < 0.05).
Pathway studio analysis of common pathways represented by DE genes in a group
Pathway Studio 9.0 (Ariadne Genomics, Rockville, MD, US). This software was used to find common pathways determined by differentially expressed genes and subsequently gain a better insight in Yorkshire versus Meishan day 12 endometrium expression profiles. Pathway Studio is a text-mining tool that scans manuscripts and abstracts from multiple biomedical web resources to establish known relationships. We used the differentially expressed list at q < 0.01 and performed a sub-network enrichment analysis (SNEA) on the up- and down-regulated gene lists. Our analyses focused on proteins and small molecules.
Availability of the supporting data
All Affymetrix data has been submitted to GEO database with accession number GSE46332 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=kpifaowqnjczdut&acc=GSE46332).
Meishan day 12
Yorkshire day 12 endometrium
Retinoic acid receptor-related orphan receptor alpha
Nuclear receptor subfamily 2, group F, member 2
Prostaglandin endoperoxide synthase 1
Vascular endothelial growth factor
Cadherin 5, type 2 vascular endothelium
Integrin beta 3
Tight junction protein 2
Sub-network enrichment analysis.
Spencer TE, Bazer FW: Uterine and placental factors regulating conceptus growth in domestic animals. J Anim Sci. 2004, 82 E-Suppl: E4-13.
Lunney JK: Advances in swine biomedical model genomics. Int J Biol Sci. 2007, 3 (3): 179-184.
Ross JW, de Castro JP F, Zhao J, Samuel M, Walters E, Rios C, Bray-Ward P, Jones BW, Marc RE, Wang W, et al: Generation of an inbred miniature pig model of retinitis pigmentosa. Invest Ophthalmol Vis Sci. 2012, 53 (1): 501-507. 10.1167/iovs.11-8784.
Yang CX, Ross JW: Genetic Modification of Domestic Animals for Agriculture and Biomedical Applications. Biomedical Science, Engineering and Technology. Edited by: Ghista DN. 2012, Rijeka, Croatia: InTech Publishing, 697-726.
Spotter A, Distl O: Genetic approaches to the improvement of fertility traits in the pig. Vet J. 2006, 172 (2): 234-247. 10.1016/j.tvjl.2005.11.013.
Legault C, Gruand J: Amélioration de la prolificité des truies par la création d’une lignée "hyperprolifique" et l’usage de l’insemination artificielle: principe et résultats experimentaux préliminaires. J Rech Porcine Fr. 1976, 8: 201-206.
Rothschild M, Bidanel J: Biology and genetics of reproduction. The Genetics of the pig. Edited by: Rothschild MF, Ruvinsky A. 1998, Cambridge, UK: CAB international, University Press, 313-343.
Rothschild M, Jacobson C, Vaske D, Tuggle C, Wang L, Short T, Eckardt G, Sasaki S, Vincent A, McLaren D, et al: The estrogen receptor locus is associated with a major gene influencing litter size in pigs. Proc Natl Acad Sci USA. 1996, 93 (1): 201-205. 10.1073/pnas.93.1.201.
Pope WF: Embryonic mortality in swine. Embryonic Mortality in Domestic Species. Edited by: Zavy MT, Geisert RD. 1994, Boca Raton: CRC press, 53-77.
Yelich JV, Pomp D, Geisert RD: Detection of transcripts for retinoic acid receptors, retinol-binding protein, and transforming growth factors during rapid trophoblastic elongation in the porcine conceptus. Biol Reprod. 1997, 57 (2): 286-294. 10.1095/biolreprod57.2.286.
Ross JW, Ashworth MD, Stein DR, Couture OP, Tuggle CK, Geisert RD: Identification of differential gene expression during porcine conceptus rapid trophoblastic elongation and attachment to uterine luminal epithelium. Physiol Genomics. 2009, 36 (3): 140-148.
Geisert RD, Ross JW, Ashworth MD, White FJ, Johnson GA, DeSilva U: Maternal recognition of pregnancy signal or endocrine disruptor: the two faces of oestrogen during establishment of pregnancy in the pig. Soc Reprod Fertil Suppl. 2006, 62: 131-145.
Ashworth MD, Ross JW, Stein DR, Allen DT, Spicer LJ, Geisert RD: Endocrine disruption of uterine insulin-like growth factor expression in the pregnant gilt. Reproduction. 2005, 130 (4): 545-551. 10.1530/rep.1.00821.
Joyce MM, Burghardt RC, Geisert RD, Burghardt JR, Hooper RN, Ross JW, Ashworth MD, Johnson GA: Pig conceptuses secrete estrogen and interferons to differentially regulate uterine STAT1 in a temporal and cell type-specific manner. Endocrinol. 2007, 148 (9): 4420-4431. 10.1210/en.2007-0505.
Ross JW, Ashworth MD, White FJ, Johnson GA, Ayoubi PJ, DeSilva U, Whitworth KM, Prather RS, Geisert RD: Premature estrogen exposure alters endometrial gene expression to disrupt pregnancy in the pig. Endocrinol. 2007, 148 (10): 4761-4773. 10.1210/en.2007-0599.
Vallee 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 (5): 1697-1706. 10.1095/biolreprod.103.019307.
Tayade C, Black GP, Fang Y, Croy BA: Differential gene expression in endometrium, endometrial lymphocytes, and trophoblasts during successful and abortive embryo implantation. J Immunol. 2006, 176 (1): 148-156.
Wessels JM, Linton NF, Croy BA, Tayade C: A review of molecular contrasts between arresting and viable porcine attachment sites. Am J Reprod Immunol. 2007, 58 (6): 470-480. 10.1111/j.1600-0897.2007.00534.x.
Fu Y, Fu J, Wang A: Association of EphA4 polymorphism with swine reproductive traits and mRNA expression of EphA4 during embryo implantation. Mol Biol Rep. 2012, 39 (3): 2689-2696. 10.1007/s11033-011-1023-8.
Wu YP, Wang AG, Li N, Fu JL, Zhao XB: Association with TGF-beta1 gene polymorphisms and reproductive performance of large white pig. Reprod Domest Anim. 2010, 45 (6): 1028-1032. 10.1111/j.1439-0531.2009.01480.x.
Lee SH, Zhao SH, Recknor JC, Nettleton D, Orley S, Kang SK, Lee BC, Hwang WS, Tuggle CK: Transcriptional profiling using a novel cDNA array identifies differential gene expression during porcine embryo elongation. Mol Reprod Dev. 2005, 71 (2): 129-139. 10.1002/mrd.20291.
Tuggle CK, Green JA, Fitzsimmons C, Woods R, Prather RS, Malchenko S, Soares BM, Kucaba T, Crouch K, Smith C, et al: EST-based gene discovery in pig: virtual expression patterns and comparative mapping to human. Mamm Genome. 2003, 14 (8): 565-579. 10.1007/s00335-002-2263-7.
Ross JW, Ashworth MD, Hurst AG, Malayer JR, Geisert 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.
Ross JW, Malayer JR, Ritchey JW, Geisert RD: Characterization of the interleukin-1beta system during porcine trophoblastic elongation and early placental attachment. Biol Reprod. 2003, 69 (4): 1251-1259. 10.1095/biolreprod.103.015842.
Ostrup E, Bauersachs S, Blum H, Wolf E, Hyttel P: Differential endometrial gene expression in pregnant and nonpregnant sows. Biol Reprod. 2010, 83 (2): 277-285. 10.1095/biolreprod.109.082321.
le Blomberg A, Schreier L, Li RW: Characteristics of peri-implantation porcine concepti population and maternal milieu influence the transcriptome profile. Mol Reprod Dev. 2010, 77 (11): 978-989. 10.1002/mrd.21253.
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.
Ford SP: Embryonic and fetal development in different genotypes in pigs. J Reprod Fertil Suppl. 1997, 52: 165-176.
Haley CS, Lee GJ: Genetic basis of prolificacy in Meishan pigs. J Reprod Fertil Suppl. 1993, 48: 247-259.
Haley CS, Lee GJ, Ritchie M: Comparative reproductive performance in Meishan and large white pigs and their crosses. Anim Sci. 1995, 60: 259-267. 10.1017/S1357729800008420.
Biensen NJ, Wilson ME, Ford SP: The impact of either a Meishan or Yorkshire uterus on Meishan or Yorkshire fetal and placental development to days 70, 90, and 110 of gestation. J Anim Sci. 1998, 76 (8): 2169-2176.
Stinckens A, Mathur P, Janssens S, Bruggeman V, Onagbesan OM, Schroyen M, Spincemaille G, Decuypere E, Georges M, Buys N: Indirect effect of IGF2 intron3 g.3072G>A mutation on prolificacy in sows. Anim Genet. 2010, 41 (5): 493-498. 10.1111/j.1365-2052.2010.02040.x.
Youngs CR, Ford SP, McGinnis LK, Anderson LH: Investigations into the control of litter size in swine: I. Comparative studies on in vitro development of Meishan and Yorkshire preimplantation embryos. J Anim Sci. 1993, 71 (6): 1561-1565.
Wilson ME, Ford SP: Differences in trophectoderm mitotic rate and P450 17alpha-hydroxylase expression between late preimplantation Meishan and Yorkshire conceptuses. Biol Reprod. 1997, 56 (2): 380-385. 10.1095/biolreprod56.2.380.
Bazer FW, Thatcher WW, Martinat-Botte F, Terqui M: Conceptus development in large white and prolific Chinese Meishan pigs. J Reprod Fertil. 1988, 84 (1): 37-42. 10.1530/jrf.0.0840037.
Biensen NJ, Wilson ME, Ford SP: The impacts of uterine environment and fetal genotype on conceptus size and placental vascularity during late gestation in pigs. J Anim Sci. 1999, 77 (4): 954-959.
Fernandez-Rodriguez A, Munoz M, Fernandez A, Pena RN, Tomas A, Noguera JL, Ovilo C, Fernandez AI: Differential gene expression in ovaries of pregnant pigs with high and low prolificacy levels and identification of candidate genes for litter size. Biol Reprod. 2011, 84 (2): 299-307. 10.1095/biolreprod.110.085589.
Lockhart DJ, Dong H, Byrne MC, Follettie MT, Gallo MV, Chee MS, Mittmann M, Wang C, Kobayashi M, Horton H, et al: Expression monitoring by hybridization to high-density oligonucleotide arrays. Nat Biotechnol. 1996, 14 (13): 1675-1680. 10.1038/nbt1296-1675.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.
Zhao SH, Recknor J, Lunney JK, Nettleton D, Kuhar D, Orley S, Tuggle CK: Validation of a first-generation long-oligonucleotide microarray for transcriptional profiling in the pig. Genomics. 2005, 86 (5): 618-625. 10.1016/j.ygeno.2005.08.001.
Couture O, Callenberg K, Koul N, Pandit S, Younes R, Hu ZL, Dekkers J, Reecy J, Honavar V, Tuggle C: ANEXdb: an integrated animal ANnotation and microarray EXpression database. Mamm Genome. 2009, 20 (11–12): 768-777.
Tuggle CK, Wang Y, Couture O: Advances in swine transcriptomics. Int J Biol Sci. 2007, 3 (3): 132-152.
Pierzynski P, Reinheimer TM, Kuczynski W: Oxytocin antagonists may improve infertility treatment. Fertil Steril. 2007, 88 (1): 213-e219-222
Moraloglu O, Tonguc E, Var T, Zeyrek T, Batioglu S: Treatment with oxytocin antagonists before embryo transfer may increase implantation rates after IVF. Reprod Biomed Online. 2010, 21 (3): 338-343. 10.1016/j.rbmo.2010.04.009.
Ashworth MD, Ross JW, Hu J, White FJ, Stein DR, Desilva U, Johnson GA, Spencer TE, Geisert RD: Expression of porcine endometrial prostaglandin synthase during the estrous cycle and early pregnancy, and following endocrine disruption of pregnancy. Biol Reprod. 2006, 74 (6): 1007-1015. 10.1095/biolreprod.105.046557.
Zhang H, Wang S, Liu M, Zhang A, Wu Z, Zhang Z, Li J: Differential gene expression in the endometrium on gestation day 12 provides insight into sow prolificacy. BMC Genomics. 2013, 14: 45-10.1186/1471-2164-14-45.
Wu D, Song D, Li X, Yu M, Li C, Zhao S: Molecular characterization and identification of the E2/P4 response element in the porcine HOXA10 gene. Mol Cell Biochem. 2013, 374 (1–2): 213-222.
Blitek A, Morawska E, Kiewisz J, Ziecik AJ: Effect of conceptus secretions on HOXA10 and PTGS2 gene expression, and PGE2 release in co-cultured luminal epithelial and stromal cells of the porcine endometrium at the time of early implantation. Theriogenology. 2011, 76 (5): 954-966. 10.1016/j.theriogenology.2011.05.002.
Kyo S, Sakaguchi J, Kiyono T, Shimizu Y, Maida Y, Mizumoto Y, Mori N, Nakamura M, Takakura M, Miyake K, et al: Forkhead transcription factor FOXO1 is a direct target of progestin to inhibit endometrial epithelial cell growth. Clin Cancer Res. 2011, 17 (3): 525-537. 10.1158/1078-0432.CCR-10-1287.
Wilson ME, Ford SP: Effect of estradiol-17beta administration during the time of conceptus elongation on placental size at term in Meishan pigs. J Anim Sci. 2000, 78 (4): 1047-1052.
Vestweber D: VE-cadherin: the major endothelial adhesion molecule controlling cellular junctions and blood vessel formation. Arterioscler Thromb Vasc Biol. 2008, 28 (2): 223-232.
Tapia A, Gangi LM, Zegers-Hochschild F, Balmaceda J, Pommer R, Trejo L, Pacheco IM, Salvatierra AM, Henriquez S, Quezada M, et al: Differences in the endometrial transcript profile during the receptive period between women who were refractory to implantation and those who achieved pregnancy. Hum Reprod. 2008, 23 (2): 340-351.
Massuto DA, Kneese EC, Johnson GA, Burghardt RC, Hooper RN, Ing NH, Jaeger LA: Transforming growth factor beta (TGFB) signaling is activated during porcine implantation: proposed role for latency-associated peptide interactions with integrins at the conceptus-maternal interface. Reproduction. 2010, 139 (2): 465-478. 10.1530/REP-09-0447.
Oku N, Sasabe E, Ueta E, Yamamoto T, Osaki T: Tight junction protein claudin-1 enhances the invasive activity of oral squamous cell carcinoma cells by promoting cleavage of laminin-5 gamma2 chain via matrix metalloproteinase (MMP)-2 and membrane-type MMP-1. Cancer Res. 2006, 66 (10): 5251-5257. 10.1158/0008-5472.CAN-05-4478.
Riesewijk A, Martin J, van Os R, Horcajadas JA, Polman J, Pellicer A, Mosselman S, Simon C: Gene expression profiling of human endometrial receptivity on days LH + 2 versus LH + 7 by microarray technology. Mol Hum Reprod. 2003, 9 (5): 253-264. 10.1093/molehr/gag037.
Raghupathy R: Pregnancy: success and failure within the Th1/Th2/Th3 paradigm. Semin Immunol. 2001, 13 (4): 219-227. 10.1006/smim.2001.0316.
Makhseed M, Raghupathy R, Azizieh F, Omu A, Al-Shamali E, Ashkanani L: Th1 and Th2 cytokine profiles in recurrent aborters with successful pregnancy and with subsequent abortions. Hum Reprod. 2001, 16 (10): 2219-2226. 10.1093/humrep/16.10.2219.
Granot I, Gnainsky Y, Dekel N: Endometrial inflammation and effect on implantation improvement and pregnancy outcome. Reproduction. 2012, 144 (6): 661-668. 10.1530/REP-12-0217.
Jin LP, Fan DX, Zhang T, Guo PF, Li DJ: The costimulatory signal upregulation is associated with Th1 bias at the maternal-fetal interface in human miscarriage. Am J Reprod Immunol. 2011, 66 (4): 270-278. 10.1111/j.1600-0897.2011.00997.x.
Zhou WH, Dong L, Du MR, Zhu XY, Li DJ: Cyclosporin A improves murine pregnancy outcome in abortion-prone matings: involvement of CD80/86 and CD28/CTLA-4. Reproduction. 2008, 135 (3): 385-395. 10.1530/REP-07-0063.
Du MR, Dong L, Zhou WH, Yan FT, Li DJ: Cyclosporin a improves pregnancy outcome by promoting functions of trophoblasts and inducing maternal tolerance to the allogeneic fetus in abortion-prone matings in the mouse. Biol Reprod. 2007, 76 (5): 906-914. 10.1095/biolreprod.106.056648.
Matsushita M, Thiel S, Jensenius JC, Terai I, Fujita T: Proteolytic activities of two types of mannose-binding lectin-associated serine protease. J Immunol. 2000, 165 (5): 2637-2642.
Tayade C, Esadeg S, Fang Y, Croy BA: Functions of alpha 2 macroglobulins in pregnancy. Mol Cell Endocrinol. 2005, 245 (1–2): 60-66.
Sass-Kuhn SP, Moqbel R, Mackay JA, Cromwell O, Kay AB: Human granulocyte/pollen-binding protein. Recognition and identification as transferrin. J Clin Invest. 1984, 73 (1): 202-210. 10.1172/JCI111192.
Chae JI, Kim J, Lee SG, Jeon YJ, Kim DW, Soh Y, Seo KS, Lee HK, Choi NJ, Ryu J, et al: Proteomic analysis of pregnancy-related proteins from pig uterus endometrium during pregnancy. Proteome Sci. 2011, 9: 41-10.1186/1477-5956-9-41.
Jetten AM: Recent advances in the mechanisms of action and physiological functions of the retinoid-related orphan receptors (RORs). Curr Drug Targets Inflamm Allergy. 2004, 3 (4): 395-412. 10.2174/1568010042634497.
Zenri F, Hiroi H, Momoeda M, Tsutsumi R, Hosokawa Y, Koizumi M, Nakae H, Osuga Y, Yano T, Taketani Y: Expression of retinoic acid-related orphan receptor alpha and its responsive genes in human endometrium regulated by cholesterol sulfate. J Steroid Biochem Mol Biol. 2012, 128 (1–2): 21-28.
Laursen KB, Mongan NP, Zhuang Y, Ng MM, Benoit YD, Gudas LJ: Polycomb recruitment attenuates retinoic acid-induced transcription of the bivalent NR2F1 gene. Nucleic Acids Res. 2013, 41 (13): 6430-6443. 10.1093/nar/gkt367.
Kim EJ, Yoo YG, Yang WK, Lim YS, Na TY, Lee IK, Lee MO: Transcriptional activation of HIF-1 by RORalpha and its role in hypoxia signaling. Arterioscler Thromb Vasc Biol. 2008, 28 (10): 1796-1802. 10.1161/ATVBAHA.108.171546.
Koos RD, Kazi AA, Roberson MS, Jones JM: New insight into the transcriptional regulation of vascular endothelial growth factor expression in the endometrium by estrogen and relaxin. Ann N Y Acad Sci. 2005, 1041: 233-247. 10.1196/annals.1282.037.
Dos Santos E, Serazin V, Morvan C, Torre A, Wainer R, de Mazancourt P, Dieudonne MN: Adiponectin and leptin systems in human endometrium during window of implantation. Fertil Steril. 2012, 97 (3): 771-778. 10.1016/j.fertnstert.2011.12.042. e771
Tanaka K, Minoura H, Isobe T, Yonaha H, Kawato H, Wang DF, Yoshida T, Kojima M, Kangawa K, Toyoda N: Ghrelin is involved in the decidualization of human endometrial stromal cells. J Clin Endocrinol Metab. 2003, 88 (5): 2335-2340. 10.1210/jc.2002-021024.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.
Dawson K, Rodriguez RL, Malyj W: Sample phenotype clusters in high-density oligonucleotide microarray data sets are revealed using Isomap, a nonlinear algorithm. BMC Bioinformatics. 2005, 6: 195-10.1186/1471-2105-6-195.
This work was supported by USDA NRI 2004-35205-14202 and the National Transgenic Project of China (grant no. 2011ZX08006-005)
The authors have declared no competing interests.
TG performed the DAVID analyses, wrote the majority of the first draft of the manuscript and edited the final version. MJZ wrote part of the draft and assisted TG in making heatmap for microarray differentiation genes. MS ran the Pathway Studio analysis and provided interpretation and text for that section. LQ and DN performed the statistical analysis of the microarray data. DK and JKL performed the Q-PCR analysis. SHZ, JWR and CKT collected tissues and assisted in interpretation of the data and writing the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Gu, T., Zhu, Mj., Schroyen, M. et al. Endometrial gene expression profiling in pregnant Meishan and Yorkshire pigs on day 12 of gestation. BMC Genomics 15, 156 (2014). https://doi.org/10.1186/1471-2164-15-156