Skip to main content
Download PDF

Detection of differentially expressed genes between Erhualian and Large White placentas on day 75 and 90 of gestation

BMC Genomics volume 10, Article number: 337 (2009) Cite this article

Abstract

Background

Placental efficiency is strongly associated with litter size, fetal weight and prenatal mortality. Together with its rapid growth during late gestation, the Large White pig breed shows a significant increase in placental size and weight, but this does not occur in the highly prolific Chinese pig breeds. To understand the molecular basis of placental development during late gestation in Chinese indigenous and Western breeds with different placental efficiency, female placental samples were collected from six pregnant Erhualian gilts at gestation day 75 (E75) and day 90 (E90) and from six pregnant Large White gilts at gestation day 75 (L75) and day 90 (L90). Two female placentas from one sow were used to extract RNA and then pooled in equal volumes. Twelve pooled samples were hybridized to the porcine Affymetrix GeneChip.

Results

A total of 226 and 577 transcripts were detected that were differentially expressed between E75 and L75 and between E90 and L90 (p < 0.01, q < 0.2), respectively. Gene Ontology (GO) analysis revealed that these genes belong to the class of genes that participate in angiogenesis and development. Real-time RT-PCR confirmed the differential expression of eight selected genes. Significant differential expression of five genes in the VEGF pathway was also detected between the breeds. A search of chromosomal location revealed that 44 differentially expressed genes located to QTL regions related to reproduction. Differential expression of six candidate imprinted genes was also confirmed. Three of the six genes (PLAGL1, DIRAS3, and SLC38A4) showed monoallelic expression in the porcine placenta.

Conclusion

Our study detected many genes that showed differential expression between placentas of two divergent breed of pigs, and confirmed the imprinting of three genes. These findings help to elucidate the genetic control of placental efficiency and improve the understanding of placental development.

Background

The ratio of birth weight to placental weight can be used as a measure of placental efficiency [1]. It is determined by many factors, including the thickness and surface area of the placenta, its vascular density, and the number and activity of transporters [2, 3]. In eutherian mammals, intrauterine growth shows a balance between fetal growth and the placental supply of nutrients and oxygen [4]. Fetal body weight is correlated positively with placental weight. In fact, birth size affects the long-term health of an individual, and is critical in determining life expectancy. In humans, smaller neonates are less likely to survive at birth and have a greater susceptibility to disease [5]. Thus, exploration of the genetic factors that regulate placental efficiency is an important research area.

The porcine placenta is chorioallantoic; it originates from the trophoblast and inner cell mass with no trophoblastic invasion of uterine vessels [6]. It is responsible for the exchange of respiratory gases, nutrients, and waste products between the maternal and the fetal systems. In fact, prenatal fetal organs do not participate in any nutrient metabolic pathways; all their metabolic demands are supplied by transplacental exchange from mother to fetus [7]. The placenta is a provisional organ, which only emerges during gestation. It secretes a variety of steroid and protein hormones that act in a paracrine manner on the endometrium and fetus to elicit changes in gene expression that support the growth and development of the conceptus.

High placental efficiency is thought to allow smaller placentas to maintain relatively larger fetuses, thereby contributing to higher uterine capacity and litter size [8, 9]. In fact, placental insufficiency is the primary mechanism through which intrauterine crowding of fetuses decreases fetal survival [10]. In the pig, the period from gestation day 75 to gestation day 90 is an important stage for placental and fetal development, during which the fetuses grow rapidly and need adequate nutrition. Vonnahme et al. [11] found that if one of two adjacent fetuses dies in the uterus, in the Meishan pig the other conceptus fails to increase its placental weight or surface area. In contrast, the adjacent conceptus of a Large White pig accelerates its placental growth. This indicates that Meishan and Large White conceptuses use different strategies in the competition for nutrients during gestation. The Meishan increases the vascular density of the placenta in order to obtain adequate nutrition, while the Large White conceptus accelerates its placental growth instead of enhancing placental efficiency [12]. Previous studies that involved translocation of embryos from Meishan to Large White sows indicated that the number of conceptuses transferred at gestation day 30 is generally reduced during late gestation as a result of the limited ability of the uterus to maintain a large number of fetuses to term [13]. It will be of great interest to understand the molecular mechanism underlying the differences in placental gene expression between Chinese indigenous and Western breeds of pig.

In this study, we investigated differences in the expression profiles of placental genes between Erhualian and Large White pigs in late gestation by use of the Affymetrix GeneChip. The Erhualian is a type of Taihu breed, and has larger litters than the Meishan. Bioinformatics analysis revealed that the differentially expressed genes are involved in important biological processes such as angiogenesis and embryonic development in utero. Real-time RT-PCR was used to confirm the differential expression of various genes. Significant differential expression of five genes in the VEGF pathway was also detected between Large White and Erhualian pigs using real-time RT-PCR. The differentially expressed genes were also mapped to QTL regions related to reproduction traits such as number born alive, fully formed piglets, number of stillborn piglets, etc. Differential expression of candidate imprinted genes was also detected, and the imprinting status of three genes was determined. Our results offer new insight towards an understanding of the molecular basis of placenta efficiency.

Results

Transcriptomes in the placenta

There are 24,123 probesets, which represent 20,201 transcripts and 198 controls, on the Affymetrix porcine GeneChip. The transcriptome analysis indicated that 7,842 probesets were expressed in the placenta (GEO accession numbers: GSM299411-GSM299422). Through BLAST sequence similarity analyses, 6087 transcripts were matched to human Refseq entries. GO annotation was used to classify these genes into groups representing different biological processes. The results revealed that these genes were involved in metabolic processes, transport and developmental process, etc. (Figure 1A). It is interesting that many genes had high and stable expression in E75, E90, L75 and L90. The strong expression of these genes in the placenta suggests that they may have important roles in placental development or nutrient transportation. We further checked the relative expression of some homologous genes in data obtained from human http://symatlas.gnf.org/SymAtlas/symform. Some genes that are specifically expressed in the placenta (EBI3, HSD3B1, HSD17B1), or are specifically expressed in a few tissues including the placenta, in other species [14], such as FOXF2, COL1A1, COL3A1, PERP, SPP1 and PAPPA, were detected to have high expression in our study.

Figure 1
figure 1

Gene Ontology (GO) molecular function classification of all genes expressed in placenta and differentially expressed genes. (A) All expressed genes in placenta. (B) Differentially expressed genes in E75 vs L75. (C) Differentially expressed genes in E90 vs L90. The x-axis shows the gene percentage of each GO category with regard to the placenta transcripts or the declared differentially expressed genes and the y-axis represents each GO category.

Full size image

Identification of differentially expressed genes

We identified 226 (p < 0.01 and q < 0.2) transcripts that were differentially expressed between E75 and L75 (see additional file 1). Among these genes, 160 transcripts were matched to human Refseq entries through BLAST sequence similarity analyses, while the other 29% of the transcripts were non-annotated. Of the differentially expressed genes, 94 transcripts were expressed at higher levels in E75, while 132 transcripts were expressed more abundantly in L75.

Between E90 and L90, 577 (p < 0.01 and q < 0.2) differentially expressed transcripts were found (see additional file 2). Among these genes, 435 transcripts were matched to human Refseq entries through BLAST sequence similarity analyses, while the other 25% of transcripts were non-annotated. In total, 277 transcripts were expressed at higher levels in E90, and 300 transcripts were expressed more abundantly in L90. A total of 102 transcripts were detected to have differential expression both in E75 vs L75 and in E90 vs L90. But only 96 transcripts had the same expression trends in pigs of the two gestational ages.

Real-time RT-PCR verification of differentially expressed genes

Thirteen genes (ALDH1A1, DIO3, DIRAS3, PLAGL1, PON2, ASCL2, WIF1, SLC38A4, VEGF, VEGFR-1, VEGFR-2, VE-cadherin and β-arrestin2) were selected to confirm the expression by the use of real-time RT-PCR. Among these genes, eight genes were used to validate the microarray data. The results indicated that expression patterns of these genes were consistent with the microarray (Pearson correlation coefficient >0.75, Figure 2A). The other five genes play important roles in placental angiogenesis and vascular permeability, so we also detected their expression by real-time RT-PCR. Four genes showed more than two-fold differences between breeds in the microarray analysis, but were not statistically significant. Gene β-arrestin2 is another important gene in the VEGF pathway but is not on the microarray. Expression of VEGFR-1 was significantly lower in E75 than L75, but significantly higher in E90 than L90. VEGFR-2 showed significantly higher expression in E90 than in L90. The expression levels of the VE-cadherin and β-arrestin2 genes were significantly lower in L90 than L75.

Figure 2
figure 2

Real time RT-PCR results of 13 genes. (A) Validation of the Microarray data by real-time RT-PCR analysis of eight representative genes. Expression levels of eight genes were detected in E75, L75, E90 and L90 by real-time RT-PCR (red) and microarray (blue). R represents the Pearson correlation coefficient. (B) real-time RT-PCR results of genes related to vascular development and vascular permeability. The x-axis represents the genes and the y-axis shows the fold change in expression. Different superscript alphanumeric characters indicate a statistically significant difference at p < 0.05.

Full size image

Cluster analysis

To gain insight into similarities at the transcriptome scale among placentas from conceptuses of two breeds and two ages, data from all the differentially expressed genes in the placentas were used in a systematic cluster analysis. The results showed that L75 and L90 were initially clustered together because their expression profiles were similar; E75 and E90 were clustered in another class (Figure 3).

Figure 3
figure 3

Hierarchical cluster of differentially expressed genes. We have performed a data adjustment (median center and normalization) in the cluster analysis. The color codes of red, white, black and dark green represented expression levels of high, average, low and absent respectively. A detailed view of the genes expression levels in clustering patterns is shown in the plot areas from A to E.

Full size image

GO analysis

GO classification on the basis of biological process was used to categorize 226 genes that were differentially expressed between E75 and L75. The results indicated that the proteins encoded by these genes are associated with metabolic process (56.5%), biological regulation (32.5%), developmental process (21.4%) and electron transport (5.2%) (Figure 1B). The GO biological process classification of 577 genes that were differentially expressed between E90 and L90 indicated that more genes encoded proteins associated with cellular process (65.5%), metabolic process (51.0%), developmental process (22.6%) and transport (16.2%) (Figure 1C).

Expression analysis of candidate imprinted genes

Imprinted genes, which are expressed in a monoallelic fashion depending on their parental origin, play important roles in mammalian fetal development, growth and behavior. They also influence the growth of the placenta. A homolog search using the available sequences of imprinted genes in the human and mouse genomes http://www.geneimprint.com/site/home identified 19 candidate imprinted genes on the porcine Affymetrix Genechip (see additional files 3 and 4). Eight candidate imprinted genes showed differential expression between the two breeds in the microarray analysis. One genes (DIRAS3) had significantly higher expression in E75 than in L75, and seven genes (PON2, PLAGL1, DCN, DIO3, NAP1L5, ASCL2, SDHD) showed significant differential expression between E90 and L90.

To further investigate the allelic expression of candidate imprinted genes in the porcine placenta, cDNA sequences of from both breeds were compared. In the PLAGL1 gene, one single nucleotide polymorphism (SNP) (T/C) was detected at position 1428 (FJ746562). In the SLC38A4 gene, an SNP (C/T) was detected at position 21 (FJ746563), and in the DIRAS3 gene, one SNP (T/G) was revealed at position 428 (NM_001044598). Restriction enzyme digestion (Taq I, Btg I, and Apek I) were used to digest the genomic DNA and cDNA of placentas. At least three heterozygotes were used to confirm the expression of imprinted genes, and the results showed that PLAGL1 and SLC38A4 were paternal expressed. Unfortunately, we could not detect the parental source of DIRAS3 gene, because that maternal genetype also was heterozygotes (Figure 4). No cSNPs were detected in other candidate imprinted genes, thus the imprinting status was not determined for the remaining genes.

Figure 4
figure 4

Imprinting status of PLAGL1 , SLC38A4 and DIRAS3 genes in porcine placentas. The results indicated that PLAGL1 and SLC38A4 genes are paternal expression in porcine placenta. DIRAS3 gene is monoallelic expression in porcine placenta.

Full size image

In silico mapping of differentially expressed genes to reproduction QTL regions

Through in silico analysis, 44 differentially expressed transcripts were found in porcine reproduction QTL chromosomal regions (see additional file 5). The traits associated with these QTL regions included ovulation rate (33 transcripts), age at puberty (8 transcripts), and uterine capacity (3 transcripts) http://www.animalgenome.org/cgi-bin/QTLdb/SS/browse (see additional file 5).

Discussion

During late gestation, Erhualian and Large White sows show different uterine responses to the rapid growth of the fetuses. For example, the placenta of the Erhualian pig stops growing, whereas Large White pigs continue to increase the surface area of the placenta [12, 15, 16]. This indicates indirectly that there may be many genes involved in the differences in placental development between Erhualian and Large White pigs (Table 1).

Table 1 Candidate genes may have important function in placenta.
Full size table

Genes involved in vascular development and permeability pathway

Placental vascular development and permeability are important in transportation of nutrients, respiratory gases, and waste products between the maternal and fetal circulations. In many cases, to keep pace with fetal growth, uterine and umbilical blood flow increases markedly during late gestation [17–19]. Insufficient blood flow leads to nutritional deficiency, and compromises fetal growth, which affects not only the neonate but also the health and productivity of the resulting adult throughout life [20]. Thus it is important to investigate the expression of genes in the VEGF pathway in order to understand placental efficiency. The results of this study showed that the expression of these genes was significantly decreased from L75 to L90, while no significant changes were found between E75 and E90 (except in VEGFR-2). This indicates that vascular development and permeability may be stable in E75 and E90, but that the Large White showed notable variation accompanying further development of the placenta during this period.

Vascular endothelial growth factor (VEGF) is produced and secreted in the placenta; it serves as an angiogenic and permeability enhancing factor in several species, including the pig [21]. It has two specific receptors (VEGFR-1 and VEGFR-2). The interaction of VEGF and VEGFR-2 could improve the mitogenic ability of VEGF, and lead to survival, migration, and differentiation of endothelial cells, and mediation of vascular permeability [22]. Compared with VEGFR-1, VEGFR-2 plays a major role in the biological activity of VEGF [23]. A number of signal transduction pathways downstream of VEGFR-2 play important roles, especially in vasculogenesis. In the results of our real-time RT-PCR, the expression of VEGFR-2 had no significant difference between E75 and L75, but was significantly higher in E90 than in L90. This indicates that the Erhualian placenta maybe have higher blood vessel density comparing with the Large White placenta. Thus Erhualian pigs can provide sufficient nutrients despite a markedly smaller placenta. Of interest, the expression of VEGFR-1 in the E75 placenta was lower than that in the L75, but at gestation day 90 in the Large White placenta, this gene was down-regulated and showed no significant difference from that of the Erhualian. VEGFR-1 may decrease the binding of VEGF to VEGFR-2 by sequestration of VEGF, thereby preventing its binding to VEGFR-2 [24]. Of course, VEGFR-1 and -2 may also form heterodimers whose signal efficiency is equal to or even greater than that of the VEGFR-2 homodimer [25]. Even so, recent data indicate that VEGFR-1 plays a positive regulatory role in vascular permeability [26]. VE-cadherin1, a cell-cell adhesion molecule specific to the endothelium, also interacts with VEGFR-2. This complex may promote cell internalization, lead to decomposition of intercellular junctions, and promote endothelial permeability [27]. β-arrestin 2 also aids endocytosis of VE-cadherin. A knockdown experiment indicated that β-arrestin 2-/- inhibits the effect of VEGF on both endocytosis of VE-cadherin and endothelial permeability [28]. The expression of β-arrestin 2 and VE-cadherin in placentas from Erhualian pigs was lower than in Large White placentas on day 75 of gestation, but there was no significant difference on day 90.

Various genes related to vascular development were also detected in the microarray analysis. HOXA13, which is essential for placental vascular patterning and labyrinth endothelial specification, showed no differential expression between breeds at gestation day 75, but had higher expression in E90 than in L90 [29]. TSP-1, which is a negative regulator of angiogenesis, showed higher expression in L75 than in E75 [30]. HAND2 is required for vascular development and regulation of angiogenesis [31]. It was expressed at a higher level in the Erhualian than the Large White. All these results indicate that the vascular development and permeability of the placenta is different in the Erhualian and Large White. During placental development in the Large White, this difference increases. This may explain the different response of the placenta of Erhualian and Large White pigs to rapid fetal development in late gestation.

Hormone related genes

The placenta is an endocrine organ which synthesizes many kinds of hormones that regulate its development and alter maternal physiology to support gestation. It can influence and regulate the expression of the fetal hormone genes [32]. Regulation of endocrine signals in the placenta, especially in late gestation, may over-ride maternal signals, owing to the commencement of accelerated growth of the fetus. This parasecretion may lead to the death of both mother and offspring [33]. In our results, some specific hormone related genes were shown to have differential expression between the two breeds of pig. Glucocorticoids are steroid hormones that play important roles in metabolic, immunological and development, in particular in stimulation of several processes that collectively serve to increase and maintain concentrations of glucose in blood [34]. They exert their actions by binding to the complexes which is composed with hormone-receptor and glucocorticoid responsive elements complexes [35]. In our experiment, we found that the level of expression of GR (glucocorticoid receptor) was higher expression in E90 than L90. The PTGS1 gene, which is a glucocorticoid responsive element, also could regulate angiogenesis in endothelial cells [36, 37]. Its expression level was both higher in E75 and E90 than L75 and L90. MED14 (Mediator complex subunit 14), which is the GR coactivators, also showed higher expression in E90 than in L90 [38]. All these results showed that the Erhualian placenta may be more efficiently to able obtain glucose from the uterus than that of the Large White. Moreover, the level of expression of the. estrogen receptor 1, mineralocorticoid receptor, growth hormone regulated TBC protein 1, thyroid hormone receptor associated protein 2 and thyroid hormone receptor associated protein 3 all were significantly higher in the Erhualian than in the Large White at the same stage of gestation. All these differentially expressed genes are major regulators of the function of hormone gene. Therefore, we presume that the Erhualian placenta can provide an environment that supports optimal absorption of nutrients. The placenta of this breed may be able to extract nutrition from the uterus more efficiently than that of the Large White.

Candidate genes for placental development and nutrient transportation

A reduction in placental size may mediate negative regulation of maternal nutrition on the birth weight and the number of secondary muscle fibers in the offspring [39, 40]. During late gestation, the Large White placenta increases in size and surface area. In contrast, the Erhualian placenta remains unchanged in size and weight [12]. The genes that are involved with placental development are poorly understood. In our study, the expression of the gene CHD8 was higher in the Large White than in the Erhualian on gestation day 75 and 90. In the mouse, Nishiyama et al. [41] found that Chd8 could inhibit p53-dependent apoptosis through binding to p53. Chd8(-/-) mice died early during embryogenesis. It has strong expression in the placenta and other reproductive tissues. Therefore we presume that this gene may play roles in placental development. In fact, there are some analogous genes among those investigated in our study. For example, FUT8 is associated with in-utero development of the embryo [42]. It showed very high expression in the placenta and may potentially be the genes that influence placental development. The Prep1 gene regulates multiple aspects of embryonic development in the mouse through a Pbx-Meis network [43]. It was also found to have higher expression in Erhualian placentas than in those of Large White pigs. The GO molecular function classification showed that the percentage of differentially expressed genes that encoded proteins associated with developmental process was 21.4% on gestation day 75 and 22.6% on gestation day 90. All these results provide indications of potential genes that warrant further study in investigating the disparity in placental development between indigenous Chinese and exotic pig breeds.

One of the primary regulators of maternofetal transfer of nutrition is the density of transporter proteins in the placenta. In our results, we also detected some genes in the solute carrier family. SLC16A10, which is a member of a family of plasma membrane amino acid transporters, mediates the Na(+)-independent transport of aromatic amino acids across the plasma membrane [44]. It showed up-regulation in the Erhualian placenta on both gestation day 75 and gestation day 90, in comparison to that of the Large White. Slc2a12 has been designated as a new member of the glucose transporter family in the mouse [45]. SLC25A24 plays important roles in the net uptake and efflux of adenine nucleotides [46]. SLC1A1 also participates in glutamate transport [47]. These genes showed no difference in expression between E75 and L75, but showed higher expression in E90 than in L90. Differential expression was also detected in genes of other solute carrier families. The SLC20A1 gene was expressed at higher levels in L75 than in E75, but in placentas on gestation day 90, it had higher expression in E90. All these results suggest that Erhualian placentas can supply more transporters and transport nutrients more efficiently.

Imprinted genes

The genome of the fetus inherits two copies of most genes, one from the mother and the other from the father. Generally, the two single strands have similar transcription activities and equivalent function. However, imprinted genes escape this rule, and preferentially express maternally or paternally derived alleles. They have evolved in mammals because of the conflicting interests of maternal and paternal genes in relation to the transfer of nutrients from the mother to her offspring [48]. They play central roles in fetal growth through controlling the fetal demand for, and the placental supply of, maternal nutrients [49]. Detecting the imprinting status of genes, especially in the placenta, is very important. The information available is very useful in comparisons of mouse, human, cattle and pig genes because of the conservation of expression and function of imprinted genes [50]. In this study, we detected the imprinting status of three potentially imprinted genes in the porcine placenta. The SLC38A4 gene is paternal expression in the pig. This is the same as the mouse, but not cattle [51]. The Diras3 gene displays paternal-specific expression both in mice and humans [52]. We found that it is also expressed in a monoallelic fashion in the porcine placenta. The PLAGL1 gene was found to be imprinted in the porcine placenta. In the human and mouse, these genes are paternally expressed [53]. From our evaluation of the conservation of imprinting status among species, we assume that these genes have the same function and expression in the pig as in the other species. In the mouse, Diras3 inhibits growth, and transgenic expression of this gene produces individuals of small stature [54]. The Plagl1 gene encodes a growth suppressor protein that is related to developmental disorders such as growth retardation and transient neonatal diabetes mellitus (TNDM) [53, 55]. Moreover, it regulates an imprinted gene network that is critically involved in the control of embryonic growth [56]. In our study, DIRAS3 had higher expression in the placentas of E75 than L75. The expression of PLAGL1 increased from E75 to E90, but decreased from L75 to L90. The DCN gene, a strong candidate for function in placental development [57], was found to have higher expression in the Large White placenta. These expression patterns may be related to the phenotype, in that the Erhualian has a smaller fetus and placenta than the Large White pig.

Conclusion

We have reported differential gene expression in placentas at gestation day 75 and 90 from pig breeds with different placental efficiency. Cluster analysis indicated that E75 has more similarity on the transcriptome scale with E90 than with L75 and L90. This finding could indicate that from gestation day 75 to gestation day 90, Erhualian and Large White pigs show a large variance in placental development. We also identified some differentially expressed genes that are related to angiogenesis and vascular permeability, transport, fetal development and hormonal regulation. The up-regulation of genes associated with the promotion of proliferation of blood vessels and transporter proteins in Erhualian pigs may be responsible for providing sufficient nutrients to support fetal development. The results of this study provide an opportunity to elucidate the genetic control of placental efficiency and improve our understanding of placental development.

Methods

Animal and tissue preparation

All animal procedures were performed according to protocols approved by the Biological Studies Animal Care and Use Committee of Hubei Province, P. R. Six Erhualian and six Large White sows were mated with 6 boars of each corresponding breed. These boars and sows have no common grandparents, they can be treated as unrelated animals. Three sows of each breed were killed and necropsied on days 75 and 90 of gestation, respectively. The uterus was removed immediately and the occupied fetal placentas of each fetus were collected, washed briefly with PBS, and then snap frozen in liquid nitrogen for later experiments.

Porcine Affymetrix GeneChip hybridization and data analysis

Two female placentas from each sow were used to extract RNA and then the samples pooled in equal volumes. In total, 12 pools were prepared from the 12 sows, i.e. three biological replicates at each stage and for each breed were prepared for hybridization. The RNAs were sent to a commercial service for hybridization to the porcine Affymetrix GeneChip. In brief, 10 micrograms of total RNA was transcribed to first- and second-strand cDNAs according to the manufacturer's instructions (Affymetrix, Inc., Santa Clara, CA). After purification and testing on an Agilent Bioanalyser 2100 machine, the double-stranded cDNA served as a template for the in vitro transcription reaction for cRNA amplification. The cRNA was labeled with biotin using the GeneChip One-Cycle target labeling kit (Affymetrix; Expression Analysis Technical Manual). The labeled cRNA was quantified, 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). The chips were then scanned with an Affymetrix GeneChip Scanner 3000 (Affymetrix).

Data from the cel files were converted to gene signal files by MAS 5.0 (microarray analysis system 5.0, Affymetrix) default normalization. The data were normalized between slides using the quantile normalization method proposed by Bolstad et al. [58] in R software. The natural logarithm was taken for measures of expression for each chip and used for the next step of the analysis. Differentially expressed genes were identified by analyzing these normalized data using a general linear model in SAS (SAS Institute, Cary, NC) on a gene by gene basis. The model for this experiment was Yijk = μ + Bi + Dj + BDij + Eijk. Yijk is the base-e logarithm of the normalized measurement of expression level from breed i (i = 1, 2), pregnancy day j (j = 1, 2), and individual k (k = 1, 2, 3). In the equation, μ represents an overall mean value, B is the main effect for breed, D is the main effect for day of gestation, BD is the interaction effect of breed and day, and E is the stochastic error. The q values were calculated using the method of Storey and Tibshirani [59]. An estimate of the upper bound of the positive false discovery rate was represented by the largest q value in a list of genes identified to be differentially expressed. Expression on the Affymetrix GeneChip porcine genome array probeset was determined by comparing the expression signals of perfect match (PM) and mismatch (MM) 25-mer probes. In our data set, if there is more than two P (present) in the triplicate detection data for E75, L75, E90 or L90 for a probeset, then we count this gene as expressed in the placental transcriptome.

Cluster analysis and annotation of gene ontology Hierarchical cluster analysis of differentially expressed genes in the placentas was conducted using Gene Cluster 3.0 and tree view software (Stanford University, 2002). The Database for Annotation, Visualization and Integrated Discovery (DAVID 2.0 and 2.1 beta) provides a comprehensive set of tools to summarize gene annotation data visually ([60]; http://david.abcc.ncifcrf.gov/). Functional annotations were pursued for the differentially expressed genes selected. The TC accession numbers were first updated from TIGR 5.0 to TIGR 11.0 and the corresponding Human Gene IDs were pulled out so that the DAVID analysis software could be interrogated. The data were then analyzed for each of the above comparisons using DAVID 2.0 and 2.1 beta.

Real-time RT-PCR confirmation of differentially expressed genes

Real-time RT-PCR was used to verify the differential expression of eight genes that were detected by the Affymetrix GeneChip. The primers used are listed in additional file 6. Each real-time RT-PCR reaction (in 25 μL) contained 2×SYBR Green Realtime PCR Master Mix, 0.4 μM primers, and 0.5 μL of template cDNA. The cycling conditions consisted of an initial, single cycle of 5 min at 95°C, followed by 40 cycles of 30 sec at 95°C, 30 sec at 60°C, 15 sec at 72°C, and fluorescence acquisition at 83°C for 1 sec. The cDNA was synthesized using reverse transcriptase (Promega), oligo(dT) and random primers with 5 μg RNA treated with DNase I (Ambion, Austin, Texas) from the same placental samples as those used in the microarray. The PCR amplifications were performed in triplicate for each sample. The gene expression levels were quantified relative to the expression of RPL32 using Gene Expression Macro software (Bio-Rad, Richmond, CA), by employing an optimized comparative Ct (ΔΔCt) value method. Dissociation curves were generated to ensure that a single amplicon was produced for each gene. The differences in gene expression levels between groups were compared using the Student's t-test. A p value < 0.05 was considered significant. The Pearson correlation coefficient was calculated and used to estimate the correlation of real time RT-PCR results and microarray results.

Detection of imprinting

The differentially expressed candidate imprinted genes were further investigated for their imprinting status in the placental samples collected in this study. The unigene mRNA and expressed sequence tags (EST) were assembled to produce a contig, and then used to design primers using Primer 5.0. Sequencing of the PCR products from Erhualian and Large White pigs was used to detect SNPs in the cDNA region. Restriction fragment length polymorphisms (RFLP) were used to confirm the SNPs. The PCR products were separated in 2% agarose gels containing 0.5 g/ml ethidium bromide after digestion by restriction enzymes. All genomic DNA obtained from 132 piglets was used to detect heterozygous animals. Total RNA from the placenta of a heterozygous fetus was treated with the TURBO DNA-free kit (Ambion, Austin, USA) and was reverse transcribed to cDNA. The primer pairs shown in Table 1 were employed to amplify the genomic DNA and cDNA from the same heterozygous samples. The amplicons were digested with restriction enzymes. Sequencing was further conducted to confirm the genotype. The epigenetic status was determined by comparing the allelic expression of genomic DNA and cDNA from the same samples.

Integration of differentially expressed genes in porcine QTL regions

The Pig Quantitative Trait Loci (QTL) database http://www.animalgenome.org/cgi-bin/QTLdb/SS/browse has gathered all pig QTL data published during the past 10 years. The reproductive traits include total number born, total number born alive, fully formed piglets, number of stillborn, ovulation rate, age at puberty and uterine capacity. We downloaded all Affymetrix probesets that were located in these QTL regions. The differentially expressed probesets located in these regions were identified through comparison analysis.

Abbreviations

APPA:

pregnancy-associated plasma protein A

ALDH1A1:

aldehyde dehydrogenase family 1, subfamily A1

ASCL2:

achaete-scute complex homolog 2

CHD8:

chromodomain helicase DNA binding protein 8

COL1A1:

collagen, type I, alpha 1

COL3A1:

collagen, type III, alpha 1

DCN:

decorin

DIO3:

deiodinase, iodothyronine type III

DIRAS3:

DIRAS family, GTP-binding RAS-like 3

EBI3:

Epstein-Barr virus induced gene 3

ESR1:

estrogen receptor 1

EST:

expressed sequence tags

FOXF2:

forkhead box F2

FUT8:

fucosyltransferase 8

GO:

Gene Ontology

GR:

glucocorticoid receptor

HAND2:

heart and neural crest derivatives expressed 2

HOXA13:

homeo box A13

HSD17B1:

hydroxysteroid (17-beta) dehydrogenase 1

HSD3B1:

hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 1

IHH:

Indian hedgehog homolog

NAP1L5:

nucleosome assembly protein 1-like 5

PCR:

polymerase chain reaction

PERP:

PERP, TP53 apoptosis effector

PKNOX1:

PBX/knotted 1 homeobox 1

PLAGL1:

pleiomorphic adenoma gene-like 1

PON2:

paraoxonase 2

PREP1:

PBX/knotted 1 homeobox 1

PTGS1:

prostaglandin-endoperoxide synthase 1

QTL:

quantitative trait locus

SDHD:

succinate dehydrogenase complex, subunit D, integral membrane protein

SLC16A10:

solute carrier family 6, member 10

SLC1A1:

solute carrier family 1, member 1

SLC20A1:

solute carrier family 20, member 1

SLC25A24:

solute carrier family 25, member 4

Slc2a12:

solute carrier family 2, member 12

SLC38A4:

solute carrier family 38, member 4

SNP:

single nucleotide polymorphisms

SPARC:

secreted protein, acidic, cysteine-rich

SPP1:

secreted phosphoprotein 1

TSP-1:

thrombospondin 1

VEGF:

vascular endothelial growth factor

VEGFR-1:

soluble vascular endothelial growth factor receptor-1

VEGFR-2:

soluble vascular endothelial growth factor receptor-2

WIF1:

WNT inhibitory factor 1.

References

  1. Molteni RA, Stys SJ, Battaglia FC: Relationship of fetal and placental weight in human beings: fetal/placental weight ratios at various gestational ages and birth weight distributions. J Reprod Med. 1978, 21: 327-34.

    CAS  PubMed  Google Scholar 

  2. Carter AM: Factors affecting gas transfer across the placenta and the oxygen supply to the fetus. J Dev Physio. 1989, 12: 305e22-

    Google Scholar 

  3. Angiolini E, Fowden A, Coan P, Sandovici I, Smith P, Dean W, Burton G, Tycko B, Reik W, Sibley C, Constância M: Regulation of placental efficiency for nutrient transport by imprinted genes. Placenta. 2006, S98-102. 10.1016/j.placenta.2005.12.008. Suppl A

  4. Salafia CM, Zhang J, Miller RK, Charles AK, Shrout P, Sun W: Placental growth patterns affect birth weight for given placental weight. Birth Defects Res A Clin Mol Terato. 2007, 79: 281-8. 10.1002/bdra.20345.

    Article  CAS  Google Scholar 

  5. Coan PM, Angiolini E, Sandovici I, Burton GJ, Constância M, Fowden AL: Adaptations in placental nutrient transfer capacity to meet fetal growth demands depend on placental size in mice. J Physio. 2008, 586: 4567-76. 10.1113/jphysiol.2008.156133.

    Article  CAS  Google Scholar 

  6. Kaufmann P, Burton G: Anatomy and Genesis of the Placenta. The Physiology of Reproduction. 1999, Raven Press, NY

    Google Scholar 

  7. Wu G, Bazer FW, Cudd TA, Meininger CJ, Spencer TE: Maternal nutrition and fetal development. J Nutr. 2004, 134: 2169-72.

    CAS  PubMed  Google Scholar 

  8. 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: 2169-76.

    CAS  PubMed  Google Scholar 

  9. Wilson ME, Biensen NJ, Ford SP: Novel insight into the control of litter size in pigs, using placental efficiency as a selection tool. J Anim Sci. 1999, 77: 1654-8.

    CAS  PubMed  Google Scholar 

  10. Foxcroft GR, Dixon WT, Novak S, Putman CT, Town SC, Vinsky MD: The biological basis for prenatal programming of postnatal performance in pigs. J Anim Sci. 2006, 84: E105-12.

    PubMed  Google Scholar 

  11. Vonnahme KA, Wilson ME, Ford SP: onceptus competition for uterine space: different strategies exhibited by the Meishan and Yorkshire pig. J Anim Sci. 2002, 80: C1311-6.

    Google Scholar 

  12. Wilson ME, Biensen NJ, Youngs CR, Ford SP: Development of Meishan and Yorkshire littermate conceptuses in either a Meishan or Yorkshire uterine environment to day 90 of gestation and to term. Biol Reprod. 1998, 58: 905-10. 10.1095/biolreprod58.4.905.

    Article  CAS  PubMed  Google Scholar 

  13. Vallet JL, Christenson RK: Uterine space affects placental protein secretion in swine. Biol Reprod. 1993, 48: 575-84. 10.1095/biolreprod48.3.575.

    Article  CAS  PubMed  Google Scholar 

  14. 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-25. 10.1016/j.ygeno.2005.08.001.

    Article  CAS  PubMed  Google Scholar 

  15. Rothschild M, Jacobson C, Vaske D, Tuggle C, Wang L, Short T, Eckardt G, Sasaki S, Vincent A, McLaren D, Southwood O, Steen van der H, Mileham A, Plastow G: The estrogen receptor locus is associated with a major gene influencing litter size in pigs. Proc Natl Acad Sci USA. 1996, 93: 201-5. 10.1073/pnas.93.1.201.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Vonnahme KA, Wilson ME, Foxcroft GR, Ford SP: Impacts on conceptus survival in a commercial swine herd. J Anim Sci. 2002, 80: 553-9.

    CAS  PubMed  Google Scholar 

  17. Reynolds LP, Redmer DA: Utero-placental vascular development and placental function. J Anim Sci. 1995, 73: 1839-51.

    CAS  PubMed  Google Scholar 

  18. Konje JC, Howarth ES, Kaufmann P, Taylor DJ: Longitudinal quantification of uterine artery blood volume flow changes during gestation in pregnancies complicated by intrauterine growth restriction. BJOG. 2003, 110: 301-5. 10.1046/j.1471-0528.2003.t01-1-02163.x.

    Article  PubMed  Google Scholar 

  19. Reynolds LP, Caton JS, Redmer DA, Grazul-Bilska AT, Vonnahme KA, Borowicz PP, Luther JS, Wallace JM, Wu G, Spencer TE: Evidence for altered placental blood flow and vascularity in compromised pregnancies. J Physiol. 2006, 572: 51-8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Zhu MJ, Ford SP, Means WJ, Hess BW, Nathanielsz PW, Du M: Maternal nutrient restriction affects properties of skeletal muscle in offspring. J Physiol. 2006, 575: 241-50. 10.1113/jphysiol.2006.112110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Clauss M: Functions of the VEGF receptor-1 (FLT-1) in the vasculature. Trends Cardiovasc Med. 1998, 8: 241-5. 10.1016/S1050-1738(98)00015-2.

    Article  CAS  PubMed  Google Scholar 

  22. Matsumoto T, Claesson-Welsh L: VEGF receptor signal transduction. Sci STKE. 2001, 112: RE21-

    Google Scholar 

  23. Holmes K, Roberts OL, Thomas AM, Cross MJ: Vascular endothelial growth factor receptor-2: structure, function, intracellular signalling and therapeutic inhibition. Cell Signal. 2007, 19: 2003-12. 10.1016/j.cellsig.2007.05.013.

    Article  CAS  PubMed  Google Scholar 

  24. Shibuya M: Vascular endothelial growth factor receptor-1 (VEGFR-1/Flt-1): a dual regulator for angiogenesis. Angiogenesis. 2006, 9: 225-30. 10.1007/s10456-006-9055-8.

    Article  CAS  PubMed  Google Scholar 

  25. Olsson AK, Dimberg A, Kreuger J, Claesson-Welsh L: VEGF receptor signalling – in control of vascular function. Nat Rev Mol Cell Biol. 2006, 7: 359-71. 10.1038/nrm1911.

    Article  CAS  PubMed  Google Scholar 

  26. Shibuya M, Claesson-Welsh L: Signal transduction by VEGF receptors in regulation of angiogenesis and lymphangiogenesis. Exp Cell Res. 2006, 312: 549-60. 10.1016/j.yexcr.2005.11.012.

    Article  CAS  PubMed  Google Scholar 

  27. Lampugnani MG, Orsenigo F, Gagliani MC, Tacchetti C, Dejana E: Vascular endothelial cadherin controls VEGFR-2 internalization and signaling from intracellular compartments. J Cell Biol. 2006, 174: 593-604. 10.1083/jcb.200602080.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Gavard J, Gutkind JS: VEGF controls endothelial-cell permeability by promoting the beta-arrestin-dependent endocytosis of VE-cadherin. Nat Cell Biol. 2006, 8: 1223-34. 10.1038/ncb1486.

    Article  CAS  PubMed  Google Scholar 

  29. Shaut CA, Keene DR, Sorensen LK, Li DY, Stadler HS: HOXA13 Is essential for placental vascular patterning and labyrinth endothelial specification. PLoS Genet. 2008, 4: e1000073-10.1371/journal.pgen.1000073.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Zhang YW, Su Y, Volpert OV, Woude Vande GF: Hepatocyte growth factor/scatter factor mediates angiogenesis through positive VEGF and negative thrombospondin 1 regulation. Proc Natl Acad Sci USA. 2003, 100: 12718-23. 10.1073/pnas.2135113100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Morikawa Y, Cserjesi P: Extra-embryonic vasculature development is regulated by the transcription factor HAND1. Development. 2004, 131: 2195-204. 10.1242/dev.01091.

    Article  CAS  PubMed  Google Scholar 

  32. Hiden U, Maier A, Bilban M, Ghaffari-Tabrizi N, Wadsack C, Lang I, Dohr G, Desoye G: Insulin control of placental gene expression shifts from mother to foetus over the course of pregnancy. Diabetologia. 2006, 49: 123-31. 10.1007/s00125-005-0054-x.

    Article  CAS  PubMed  Google Scholar 

  33. Gootwine E: Placental hormones and fetal-placental development. Anim Reprod. 2004, 82–83: 551-66. 10.1016/j.anireprosci.2004.04.008.

    Article  Google Scholar 

  34. Mandl M, Ghaffari-Tabrizi N, Haas J, Nöhammer G, Desoye G: Differential glucocorticoid effects on proliferation and invasion of human trophoblast cell lines. Reproduction. 2006, 132: 159-67. 10.1530/rep.1.00976.

    Article  CAS  PubMed  Google Scholar 

  35. Yamamoto KR: Multilayered control of intracellular receptor function. Harvey Lect. 1995, 91: 1-19.

    PubMed  Google Scholar 

  36. Rheinlaender C, Weber SC, Sarioglu N, Strauss E, Obladen M, Koehne P: Changing expression of cyclooxygenases and prostaglandin receptor EP4 during development of the human ductus arteriosus. Pediatr Res. 2006, 60: 270-5. 10.1203/01.pdr.0000233066.28496.7c.

    Article  CAS  PubMed  Google Scholar 

  37. Ioannou IM, Tsawdaroglou N, Sekeris CE: Presence of glucocorticoid responsive elements in the mitochondrial genome. Anticancer Res. 1988, 8: 1405-9.

    CAS  PubMed  Google Scholar 

  38. Chen W, Rogatsky I, Garabedian MJ: MED14 and MED1 differentially regulate target-specific gene activation by the glucocorticoid receptor. Mol Endocrinol. 2006, 20: 560-72. 10.1210/me.2005-0318.

    Article  CAS  PubMed  Google Scholar 

  39. Michael K, Ward BS, Moore WM: Relationship of fetal to placental size: the pig model. Eur J Obstet Gynecol Reprod Biol. 1983, 16: 53-62. 10.1016/0028-2243(83)90220-4.

    Article  CAS  PubMed  Google Scholar 

  40. Dwyer CM, Stickland NC: Supplementation of a restricted maternal diet with protein or carbohydrate alone prevents a reduction in fetal muscle fibre number in the guinea-pig. Br J Nutr. 1994, 72: 173-80. 10.1079/BJN19940022.

    Article  CAS  PubMed  Google Scholar 

  41. Nishiyama M, Oshikawa K, Tsukada Y, Nakagawa T, Iemura S, Natsume T, Fan Y, Kikuchi A, Skoultchi AI, Nakayama KI: CHD8 suppresses p53-mediated apoptosis through histone H1 recruitment during early embryogenesis. Nat Cell Biol. 2009, 11: 172-82. 10.1038/ncb1831.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Roos C, Kolmer M, Mattila P, Renkonen R: Composition of Drosophila melanogaster proteome involved in fucosylated glycan metabolism. J Biol Chem. 2002, 277: 3168-75. 10.1074/jbc.M107927200.

    Article  CAS  PubMed  Google Scholar 

  43. Ferretti E, Villaescusa JC, Di Rosa P, Fernandez-Diaz LC, Longobardi E, Mazzieri R, Miccio A, Micali N, Selleri L, Ferrari G, Blasi F: Hypomorphic mutation of the TALE gene Prep1 (pKnox1) causes a major reduction of Pbx and Meis proteins and a pleiotropic embryonic phenotype. Mol Cell Biol. 2006, 26: 5650-62. 10.1128/MCB.00313-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Halestrap AP, Meredith D: The SLC16 gene family-from monocarboxylate transporters (MCTs) to aromatic amino acid transporters and beyond. Pflugers Arch. 2004, 447: 619-28. 10.1007/s00424-003-1067-2.

    Article  CAS  PubMed  Google Scholar 

  45. Rogers S, Chandler JD, Clarke AL, Petrou S, Best JD: Glucose transporter GLUT12-functional characterization in Xenopus laevis oocytes. Biochem Biophys Res Commun. 2003, 308: 422-6. 10.1016/S0006-291X(03)01417-7.

    Article  CAS  PubMed  Google Scholar 

  46. Fiermonte G, De Leonardis F, Todisco S, Palmieri L, Lasorsa FM, Palmieri F: Identification of the mitochondrial ATP-Mg/Pi transporter. Bacterial expression, reconstitution, functional characterization, and tissue distribution. J Biol Chem. 2004, 279: 30722-30. 10.1074/jbc.M400445200.

    Article  CAS  PubMed  Google Scholar 

  47. Rainesalo S, Keränen T, Saransaari P, Honkaniemi J: GABA and glutamate transporters are expressed in human platelets. Brain Res Mol Brain Res. 2005, 141: 161-5. 10.1016/j.molbrainres.2005.08.013.

    Article  CAS  PubMed  Google Scholar 

  48. Coan PM, Burton GJ, Ferguson-Smith AC: Imprinted genes in the placenta – a review. Placenta. 2005, S10-20. 10.1016/j.placenta.2004.12.009. Suppl A

  49. Kaneko-Ishino T, Kohda T, Ishino F: The regulation and biological significance of genomic imprinting in mammals. J Biochem. 2003, 133: 699-711. 10.1093/jb/mvg090.

    Article  CAS  PubMed  Google Scholar 

  50. Zhou QY, Huang JN, Xiong YZ, Zhao SH: Imprinting analyses of the porcine GATM and PEG10 genes in placentas on days 75 and 90 of gestation. Genes Genet Syst. 2007, 82: 265-9. 10.1266/ggs.82.265.

    Article  CAS  PubMed  Google Scholar 

  51. Zaitoun I, Khatib H: Assessment of genomic imprinting of SLC38A4, NNAT, NAP1L5, and H19 in cattle. BMC Gene. 2006, 7: 49-10.1186/1471-2156-7-49.

    Article  Google Scholar 

  52. Xu F, Xia W, Luo RZ, Peng H, Zhao S, Dai J, Long Y, Zou L, Le W, Liu J, Parlow AF, Hung MC, Bast RC, Yu Y: The human ARHI tumor suppressor gene inhibits lactation and growth in transgenic mice. Cancer Res. 2000, 60: 4913-20.

    CAS  PubMed  Google Scholar 

  53. Arima T, Yamasaki K, John RM, Kato K, Sakumi K, Nakabeppu Y, Wake N, Kono T: The human HYMAI/PLAGL1 differentially methylated region acts as an imprint control region in mice. Genomics. 2006, 88: 650-8. 10.1016/j.ygeno.2006.07.005.

    Article  CAS  PubMed  Google Scholar 

  54. Temple IK, James RS, Crolla JA, Sitch FL, Jacobs PA, Howell WM, Betts P, Baum JD, Shield JP: An imprinted gene(s) for diabetes. Nat Genet. 1995, 9: 110-2. 10.1038/ng0295-110.

    Article  CAS  PubMed  Google Scholar 

  55. Kas K, Voz M, Roijer E, Astrom AK, Meyen E, Stenman G, Ven Van de W: Promoter swapping between the genes for a novel zinc finger protein and bcatenin in pleiomorphic adenomas with t(3;8)(p21;q12) translocations. Nat Genet. 1997, 15: 170-4. 10.1038/ng0297-170.

    Article  CAS  PubMed  Google Scholar 

  56. Varrault A, Gueydan C, Delalbre A, Bellmann A, Houssami S, Aknin C, Severac D, Chotard L, Kahli M, Le Digarcher A, Pavlidis P, Journot L: Zac1 regulates an imprinted gene network critically involved in the control of embryonic growth. Dev Cell. 2006, 11: 598-9. 10.1016/j.devcel.2006.09.003.

    Article  Google Scholar 

  57. Mizuno Y, Sotomaru Y, Katsuzawa Y, Kono T, Meguro M, Oshimura M, Kawai J, Tomaru Y, Kiyosawa H, Nikaido I, Amanuma H, Hayashizaki Y, Okazaki Y: Asb4, Ata3, and Dcn are novel imprinted genes identified by high-throughput screening using RIKEN cDNA microarray. Biochem Biophys Res Commun. 2002, 290: 1499-505. 10.1006/bbrc.2002.6370.

    Article  CAS  PubMed  Google Scholar 

  58. Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.

    Article  CAS  PubMed  Google Scholar 

  59. Storey JD, Tibshirani R: Statistical significance for genome wide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-5. 10.1073/pnas.1530509100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This research was supported by grants from the National Natural Science Foundation of China (30571328), the Key Project of the National Basic Research and Developmental Plan (2006CB102105) of China, the Program for New Century Excellent Talents in University (NCET-05-0669), the important national science & technology specific projects (2009ZX08009-153B) and the creative team grant of Education Ministry of China (IRT0831). We thank Dr Wu Zhenfang for help with sample collection.

Author information

Authors and Affiliations

  1. Key Laboratory of Agricultural Animal Genetics, Breeding, and Reproduction of Ministry of Education & Key Laboratory of Swine Genetics and Breeding of Ministry of Agriculture, Huazhong Agricultural University, Wuhan, 430070, PR China

    Quan-Yong Zhou, Ming-Di Fang, Ting-Hua Huang, Chang-Chun Li, Mei Yu & Shu-Hong Zhao

  2. Institute of Animal Husbandry and Veterinary, Jiangxi Academy of Agricultural Science, Nanchang, 330200, PR China

    Quan-Yong Zhou

Authors
  1. Quan-Yong Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ming-Di Fang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ting-Hua Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Chang-Chun Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mei Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Shu-Hong Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Shu-Hong Zhao.

Additional information

Authors' contributions

SHZ conceived the research. CCL and MY supervised the research; all authors carried out the research. MDF and THH analyzed the data. QYZ analyzed the data, carried out the real time RT-PCR, and wrote the first version of the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

12864_2008_2221_MOESM1_ESM.xls

Additional file 1: Comparison of normalized expression levels between E75 and L75. The data provided represent the normalized expression levels of E75 and L75. Statistical analysis indicated that 226 transcripts were differentially expressed between E75 and L75 (p < 0.01, q < 0.2). (XLS 4 MB)

12864_2008_2221_MOESM2_ESM.xls

Additional file 2: Comparison of normalized expression levels between E90 and L90. The data provided represent the normalized expression levels of E90 and L90. Statistical analysis indicated that 577 transcripts were differentially expressed between E90 and L90 (p < 0.01, q < 0.2). (XLS 4 MB)

12864_2008_2221_MOESM3_ESM.xls

Additional file 3: Comparison of normalized expression levels of candidate imprinted genes between E75 and L75. The data listed the normalized expression levels of imprinted genes in E75 and L75. (XLS 28 KB)

12864_2008_2221_MOESM4_ESM.xls

Additional file 4: Comparison of normalized expression levels of candidate imprinted genes between E90 and L90. The data listed the normalized expression levels of imprinted genes in E90 and L90. (XLS 28 KB)

12864_2008_2221_MOESM5_ESM.xls

Additional file 5: Differentially expressed genes located in reproduction QTL regions. The file listed all differentially expressed genes located in reproduction QTL regions. (XLS 20 KB)

12864_2008_2221_MOESM6_ESM.xls

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

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zhou, QY., Fang, MD., Huang, TH. et al. Detection of differentially expressed genes between Erhualian and Large White placentas on day 75 and 90 of gestation. BMC Genomics 10, 337 (2009). https://doi.org/10.1186/1471-2164-10-337

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-337

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us