Skip to main content

Advertisement

  • Research article
  • Open Access

The identification of differentially expressed genes between extremes of placental efficiency in maternal line gilts on day 95 of gestation

BMC Genomics201920:254

https://doi.org/10.1186/s12864-019-5626-0

  • Received: 20 September 2018
  • Accepted: 20 March 2019
  • Published:

Abstract

Background

Placental efficiency (PE) describes the relationship between placental and fetal weights (fetal wt/placental wt). Within litters, PE can vary drastically, resulting in similarly sized pigs associated with differently sized placentas, up to a 25% weight difference. However, the mechanisms enabling the smaller placenta to grow a comparable littermate are unknown. To elucidate potential mechanisms, morphological measurements and gene expression profiles in placental and associated endometrial tissues of high PE and low PE feto-placental units were compared. Tissue samples were obtained from eight maternal line gilts during gestational day 95 ovario-hysterectomies. RNA was extracted from tissues of feto-placental units with the highest and lowest PE in each litter and sequenced.

Results

Morphological measurements, except placental weight, were not different (P > 0.05) between high and low PE. No DEG were identified in the endometrium and 214 DEG were identified in the placenta (FDR < 0.1), of which 48% were upregulated and 52% were downregulated. Gene ontology (GO) analysis revealed that a large percentage of DEG were involved in catalytic activity, binding, transporter activity, metabolism, biological regulation, and localization. Four GO terms were enriched in the upregulated genes and no terms were enriched in the downregulated genes (FDR < 0.05). Eight statistically significant correlations (P < 0.05) were identified between the morphological measurements and DEG.

Conclusion

Morphological measures between high and low PE verified comparisons were of similarly sized pigs grown on different sized placentas, and indicated that any negative effects of a reduced placental size on fetal growth were not evident by day 95. The identification of DEG in the placenta, but absence of DEG in the endometrium confirmed that the placenta responds to the fetus. The GO analyses provided evidence that extremes of PE are differentially regulated, affecting components of placental transport capacity like nutrient transport and blood flow. However, alternative GO terms were identified, indicating the complexity of the relationship between placental and fetal weights. These findings support the use of PE as a marker of placental function and provide novel insight into the genetic control of PE, but further research is required to make PE production applicable.

Keywords

  • Endometrium
  • Fetal weight
  • Gene expression
  • Pigs
  • Placenta
  • Placental efficiency
  • Placental weight

Background

Placental efficiency (PE), quantified by the ratio of fetal or birth weight to placental weight, is commonly used as a marker of placental function in humans and animals. The ratio reflects grams of fetus produced per gram of placenta [1]. In general, it is assumed that high PE values associated with averaged sized fetuses represent placentas with a greater nutrient transport capacity, while low PE values associated with growth restricted fetuses represent placentas with a reduced nutrient transport capacity or a failure to adapt.

In humans, PE (birth wt/placental wt) is reduced during pregnancy complications like fetal growth restriction, small for gestational age, gestational diabetes mellitus, and pre-eclampsia [2]. Therefore, the ratio is often used to predict abnormal fetal growth, and consequently health later in life. However, whether alterations in PE truly reflect adaptations in human placental nutrient transport capacity remains unclear. The most compelling evidence for an association between the two is in mice, with conflicting reports in humans [2]. Coan and others [3] evaluated placental nutrient transport capacity in mice with natural variations in placental size to determine if the smallest placenta in a litter of appropriately grown fetuses was the most efficient. The authors reported near term fetuses with lighter placentas were of comparable weight to fetuses with heavier placentas, and thus, PE was greater in the lightest placentas. Additionally, expression of Slc2a1, a glucose transporter gene, and Slc38a2, an amino acid transporter gene, were upregulated in the lightest placentas, providing evidence that high PE placentas adapt to meet the nutrient demands of the growing fetus.

Natural variations in PE are also apparent in pigs, a litter bearing species. These variations are not only between, but also within, breeds and even within litters [4]. Within a litter, PE can vary drastically, resulting in similarly sized pigs grown on very different placentas, with up to a 25% weight difference [5]. High PE placenta are smaller in size than low PE placenta, thus, high PE placentas occupy less space in the uterus and still grow an averaged sized littermate. While the use of PE as a selection tool to increase litter size has been debated [6], PE may provide an opportunity to optimize reproductive performance. The average litter size of U.S. production breeds has continued to increase over time and is currently 10.6 [7], but ovulation rates [8] and teat numbers [9] indicate the maximum has not been reached. At the same time, increases in litter size have resulted in lower birth weights, less uniform litters, and greater preweaning mortalities, minimizing the benefits of producing more pigs per litter. Increasing PE in these larger litters may normalize birth weights and, as a result, increase preweaning survival. In fact, Vernunft and others [10] reported on the relationship between placental size and measures of reproductive performance in modern Landrace sows. The authors concluded sows with larger litters and shorter placentas can rescue placental function. However, the compensatory mechanisms driving the growth of an adequately sized fetus on a smaller placenta are still being investigated.

In comparisons of breeds differing in PE, the increased efficiency of high PE placentas was attributed to greater vascularity [11]. Yet, variations in PE within production breed litters on day 90 of gestation could not be attributed to differences in vascular density (VD) despite increased expression of vascular endothelial growth factor and associated receptors in high PE placentas [12]. Recently, Krombeen and others [5] reported placental VD was positively related to PE on day 110 of gestation in maternal line gilts. The results of Vonnahme and Ford [12] in conjunction with Krombeen and others [5] suggest morphological adaptations, like increases in VD, may occur later in gestation (day 90 to term) to maintain fetal growth when placental size is reduced.

Krombeen and others [5] also investigated the relationship between PE and seven genes encoding glucose, amino acid, or fatty acid transporters in the placenta and adjacent endometrium of maternal line gilts on day 70, 90, and 110 of gestation. Based on conditional effects plots, variations in PE on day 70 of gestation were related to alterations in amino acid transporter expression (SLC7A7, SLC7A1, SLC3A1) in the placenta and endometrium. On day 90 of gestation, PE was positively related to placental expression of a glucose transporter (SLC2A3) and negatively related to endometrial expression of two cationic amino acid transporters (SLC7A1 and SLC7A2) and a very long chain fatty acid transporter (SLC27A1). Near term (day 110), PE was negatively related to two amino acid transporters (SLC7A7 and SLC7A1) and a glucose transporter (SLC2A3). The results of Krombeen and others [5] agree with those of Coan and others [3] and supports the use of PE as a marker of placental function.

Although the study by Krombeen and others [5] provides evidence for an association between PE and placental nutrient transport capacity, only seven genes encoding nutrient transporters were investigated. Zhou and others [13] detected 226 and 577 differentially expressed genes on gestational days 75 and 90, respectively, between two breeds with differing PE. Similarly, Kwon and others [14] identified 588 differentially expressed genes in placentas from larger litter sizes compared to smaller litter sizes. Therefore, it is likely the compensatory mechanisms enabling comparable fetal growth despite reductions in placental size (high PE) are controlled by many genes and the interactions of those genes, as seen in Zhou and others [13] and Kwon and others [14].

However, the global expression of genes in the placenta of the high and low PE feto-placental units within a litter has not been investigated. The hypothesis of this work was that extremes of PE would be associated with differentially expressed genes (DEG) that affect fetal growth via gene products that promote growth, vascularity, and/or nutrient transport. The main objective was to determine and understand the role of gene expression profiles in placental and associated endometrial tissues of high PE and low PE feto-placental units. A secondary objective was to compare fetal and utero-placental measurements of high PE and low PE units.

Results

Fetal and utero-placental measurements

Mean litter size was 10.75 and ranged from 5 to 15. There was no association between litter size and PE (r = − 0.04, P = 0.72). Table 1 contains least square means ± SE of fetal and placental measurements. Mean placental weight was lower in the high PE group compared to the low PE group (P = 0.0002, Fig. 1), but mean fetal weight was not different between high and low PE (P = 0.5914, Fig. 1). While the effect of sex was not significant for placental weight or fetal weight, there was a significant interaction between PE and sex for placental weight (P = 0.0479, Table 1). Males had lower weight placentas than females in the high PE group, but the opposite was true in the low PE group (Fig. 2).
Table 1

Least square means ± SE of fetal and utero-placental measurements of high PE and low PE units

Measurementb

Placental efficiency

P-valuea

Low

High

PE

Sex

PE × Sex

Placental wt, g

322.21 ± 24.14

172.37 ± 25.39

0.0002*

0.7257

0.0479*

Fetal wt, g

831.28 ± 99.36

782.04 ± 106.59

0.5914

0.3647

0.4386

ISL, cm

21.08 ± 3.05

22.19 ± 3.60

0.8208

0.6657

0.9710

CRL, cm

25.88 ± 1.35

25.02 ± 1.42

0.4338

0.2303

0.7022

Girth, cm

17.58 ± 1.34

18.27 ± 2.02

0.7310

0.9900

0.4043

Heart wt, g

7.63 ± 1.06

5.92 ± 1.19

0.2071

0.7930

0.3736

Liver wt, g

21.97 ± 3.14

17.86 ± 3.38

0.1925

0.4927

0.5075

Brain wt, g

19.95 ± 0.72

19.73 ± 0.83

0.8254

0.5937

0.7535

ST wt, g

1.60 ± 0.22

1.50 ± 0.23

0.6886

0.6552

0.8463

Placental VD, %

8.19 ± 1.35

5.74 ± 1.56

0.2689

0.3511

0.4941

Endometrial VD, %

6.15 ± 1.73

7.48 ± 1.96

0.5463

0.8872

0.4179

aEffect of placental efficiency (PE), effect of sex (Sex), effect of placental efficiency by sex interaction (PE × Sex). bISL Implantation site length, CRL crown-rump length, ST semitendinosus, VD vascular density. (*) asterisk indicates P < 0.05

Fig. 1
Fig. 1

Mean placental and fetal weight. Mean placental weight and fetal weight of low PE and high PE feto-placental units on day 95 of gestation in pigs. Data presented as means ± SD. Asterisk (*) indicates P < 0.05

Fig. 2
Fig. 2

Mean placental weight by sex. Mean placental weight in low PE and high PE male and female feto-placental units on day 95 of gestation. Data presented as means ± SD

Mean implantation site length (ISL), crown-rump length (CRL), girth, heart weight, liver weight, brain weight, and semitendinosus (ST) weight were not significantly different between the high PE group and the low PE group (P = 0.8208, 0.4338, 0.7310, 0.2071, 0.1925, 0.8254, and 0.6886, respectively). There were also no significant differences in mean placental or endometrial vascular density (VD) between high PE and low PE (P = 0.2689, 0.5463, respectively). For these variables, there was no significant effect of sex and there were no significant interactions between PE and sex.

Differential gene expression and gene ontology

The number of genes expressed in the placenta and endometrium was 20,280. In total, 214 DEG (FDR < 0.1) were identified in the placenta (Fig. 3a) and 0 DEG (FDR < 0.1) were identified in the endometrium (Fig. 3b). Of the DEG in the placenta, 103 genes were upregulated (33 log fold change (log2FC) ≥ 1; 70 0 < log2FC < 1.0) and 111 genes were downregulated (49 log2FC ≤ − 1.0; 62 -1.0 < log2FC < 0). Table 2 lists a subset of the candidate genes in high PE compared to low PE placentas.
Fig. 3
Fig. 3

Gene Expression. a Volcano plot of DEG in high PE compared to low PE placental tissue on day 95 of gestation in pigs. Black dots indicate FDR > 0.10 or |log2FC| < 1. Non-black dots indicate DEG (FDR < 0.10, |log2FC| ≥ 1). Dot colors indicate log2FC range. b Volcano plot of DEG in high PE compared to low PE endometrial tissues on day 95 of gestation in pigs. Black dots indicate FDR > 0.10 or |log2FC| < 1

Table 2

Candidate genes in high PE compared to low PE placentas

Gene Symbol

Gene Name

log2FCa

FDR P-value2

Gene Ontology (MF; BP)b

SAA2

Serum Amyloid A2

4.89

0.00002

binding, molecular transducer activity; locomotion

SPHKAP

Sphingosine Kinase Type 1

2.88

0.00388

binding

DKK1

Dickkopf-1

2.80

0.00876

binding, molecular transducer activity, molecular function regulator; biological regulation

CYP4F22

Cytochrome P450 Family 4 Subfamily F Member 22

2.44

0.05077

catalytic activity; metabolic process

FBP2

Fructose-1, 6-Bisphosphatase Isoenzyme 2

2.42

0.08290

catalytic activity; metabolic process

PCDHB1

Protocadherin Beta-1 Isoform X2

2.42

0.00302

;biological adhesion

SLC45A3

Solute Carrier Family 45 Member 3

1.89

0.07267

transporter activity, metabolic process

ASIC1

Acid Sensing Ion Channel 1 Isoform X2

1.81

0.05118

transporter activity; biological regulation, localization, multicellular organismal process

CELA1

Chymotrypsin-like Elastase Family Member 1

1.55

0.04841

catalytic activity

LEP

Leptin

1.01

0.05194

growth factor activity, hormone activityc

ATP13A3

Probable Cation-Transporting ATPase 13A3

0.82

0.06425

; biological regulation

SLC4A7

Solute Carrier Family 4 Member 7

0.77

0.04207

transporter activity

CTSH

Pro-cathepsin H

0.71

0.07854

; positive regulation of angiogenesisc

SLC52A3

Solute Carrier Family 52 Member 3

0.68

0.01438

molecular transducer activity; biological regulation

SLC23A2

Solute Carrier Family 23 Member 2

0.59

0.06802

transporter activity; localization

AGR2

Anterior Gradient Protein 2 Homolog

−2.99

0.00045

; biological regulation, cellular process

EMB

Embigin

−2.06

0.01288

binding; biological adhesion, cellular process

STEAP1

Six-Transmembrane Epithelial Antigen of Prostate 1

−1.98

0.05194

catalytic activity; localization

SARDH

Sarcosine Dehydrogenase

−1.89

0.02504

catalytic activity; metabolic process

STEAP2

Six Transmembrane Epithelial Antigen of the Prostate 2

−1.72

0.07267

catalytic activity; localization

MRP4

Multidrug Resistance-Associated Protein 4-like

−1.61

0.07267

transporter activity

KCNJ2

Potassium Channel, Inwardly Rectifying Subfamily J, Member 2

−1.31

0.00875

transporter activity

ANGPT1

Angiopoietin-1

−0.75

0.07031

binding; biological regulation, biological adhesion

alog2FC log2 fold change. 2FDR P-value False discovery rate adjusted P-value, level of significance P < 0.10. bMF Molecular function, BP biological process. Genes above dotted line were among ten most upregulated or downregulated. cNo hits for Panther GO slim terms, specific GO terms used instead

The results of the functional classification analysis are presented in Fig. 4, which illustrates percentages of upregulated compared to downregulated genes designated by MF (Fig. 4a) or BP (Fig. 4b) gene ontology (GO) terms for high PE compared to low PE. The MF of catalytic activity, binding, and transporter activity had the greatest representation. More catalytic activity genes (39% vs 25.5%) and transporter activity genes (14% vs 10%) were upregulated than downregulated, while the opposite was true for binding activity genes (29% vs 37%). Additionally, a greater percentage of downregulated genes compared to upregulated genes were molecular function regulators or had transcription regulator activity (both 2% vs 10%). The BP terms with the greatest representation were metabolic process, biological regulation, and localization. Upregulated genes were involved in metabolic processes (27% vs 17%) and localization (15% vs 11%) more often than downregulated genes, while downregulated genes were involved in biological regulation more often than upregulated genes (21% vs 27%). Also, the terms developmental process (0% vs 13%) and biological adhesion (3 vs 9%) were associated with more downregulated than upregulated genes.
Fig. 4
Fig. 4

Gene Ontology Functional Classification Analysis. a Molecular function GO terms of DEG in high PE compared to low PE placentas on day 95 of gestation in pigs. b Biological Process GO terms of DEG in high PE compared to low PE placentas on day 95 of gestation in pigs. Subfigures included 160 of the 214 DEG

Table 3 contains the results of the GO enrichment analysis performed on the DEG in high PE compared low PE placentas. Four terms were significantly (FDR < 0.05) enriched in the upregulated genes and no terms were significantly enriched in the downregulated genes. The MF terms sodium-dependent multivitamin transmembrane transporter activity and nucleobase transmembrane transporter activity, the BP term nucleobase transport, and the cellular component term extracellular exosome were enriched.
Table 3

Gene Ontology Enrichment Analysis

Ontologya

Category

GO Term

No. DEG in Category

No. in Category

FDR P-value

Gene Symbolb

MF

GO:0008523

Sodium-dependent multivitamin transmembrane transporter activity

2

3

0.0486

SLC23A2, LOC102159690

MF

GO:0015205

Nucleobase transmembrane transporter activity

2

3

0.0486

SLC23A2, LOC102159690

BP

GO:0015851

Nucleobase transport

2

3

0.0486

SLC23A2, LOC102159690

CC

GO:0070062

Extracellular exosome

20

2283

0.0486

TXN, TXNDC8, MGST3, CADM4, PI16, CTSH, ABCB1, FBP2, RPL15, ITIH3, PDHB, PEBP1, PPA1, GCA, HSPE1, SPHKAP, EFHD1, ENTPD6, AK2, GLA

Significantly enriched GO terms of upregulated genes in high PE compared to low PE placentas on day 95 of gestation in pigs. Level of significance P < 0.05. aMF Molecular function, BP Biological process, CC Cellular component, GO Gene ontology, DEG Differentially expressed genes, FDR P-value False discovery rate adjusted P-value, bSLC23A2 Solute carrier family 23 member 2, LOC102159690 solute carrier family 23 member 2-like, TXN Thioredoxin, TXNDC8 Thioredoxin domain-containing protein 8, MGST3 Microsomal glutathione S-transferase 3, CADM4 Cell adhesion molecule 4, PI16 Peptidase inhibitor 16 precursor, CTSH Pro-cathepsin H, ABCB1 ATP-binding cassette subfamily B member 1 isoform X2, FBP2 Fructose-1, 6-biphosphatase isoenzyme 2, RPL15 Ribosomal protein L15, ITIH3 Inter-alpha-trypsin inhibitor heavy chain H3, PDHB Pyruvate dehydrogenase E1 component subunit beta, PEBP1 Phosphatidylethanolamine-binding protein 1, PPA1 Inorganic pyrophosphatase, GCA Grancalcin, HSPE1 10 KDa heat shock protein, mitochondrial, SPHKAP A-kinase anchoring protein, EFHD1 EF-hand domain-containing protein D1 isoform X2, ENTPD6 Ectonucleoside triphosphate diphosphohydrolase 6, AK2 Adenylate kinase 2, GLA Alpha-galactosidase A

Fetal/utero-placental measurements and differentially expressed genes

Eight significant correlations (P < 0.05) were identified between fetal and utero-placental measurements, and the 214 DEG in high PE compared to low PE placentas (Table 4), of which 5 were placental weight and 3 were PE correlations. Placental weight was positively correlated with GRINL1A complex locus 1 (GCOM1, r = 0.82, P = 0.0214), gene19656 (LOC100739517, r = 0.78, P = 0.0257), TOX high mobility group box family member 3 (TOX3, r = 0.77, P = 0.0257), and ATP-binding cassette sub-family G member 2 (ABCG2, r = 0.76, P = 0.0257), but negatively correlated with ras-related protein rab-6B (RAB6B, r = − 0.76, P = 0.0257). Placental efficiency was positively correlated with gene12188 (LOC100156118, r = 0.81, P = 0.0214), transmembrane protein 199 (TMEM199, r = 0.81, P = 0.0214), and proto-cadherin beta 1 (PCDHB1, r = 0.76, P = 0.0428).
Table 4

Significant correlations between fetal/utero-placental measurements and DEG in high PE compared to low PE placentas

Correlation Dataa

DEG Datab

Measurement* Gene

Correlation

FDR P-value

Gene Symbol

Protein

Molecular Function

log FC

FDR P-value

PW*gene1223

0.82

0.0214

GCOM1

Myocardial Zonula Adherens Protein

N/A

−0.43

0.0769

PE*gene12188

0.81

0.0214

LOC100156118

N/A

N/A

0.63

0.0484

PE*gene24847

0.81

0.0214

TMEM199

Transmembrane Protein 199

N/A

0.29

0.0421

PW*gene19656

0.78

0.0257

LOC100739517

ATP-Binding Cassette Sub-Family G Member 2

N/A

−0.61

0.0144

PW*gene13353

0.77

0.0257

TOX3

TOX High Mobility Group Box Family Member 3

Chromatin binding, Phosphoprotein binding, Protein homodimerization activity, Estrogen response element binding

−1.41

0.0183

PW*gene19659

0.76

0.0257

ABCG2

ATP-Binding Cassette Sub-Family G Member 2

ATP binding, ATPase activity coupled to transmembrane movement of substance, Cholesterol transporter activity

−0.77

0.0039

PE*gene6367

0.76

0.0428

PCDHB1

Protocadherin Beta-1 Isoform X2

Calcium ion binding

2.42

0.0030

PW*gene26368

−0.76

0.0257

RAB6B

Ras-Related Protein Rab-6B

GTP binding, GTPase activity, Myosin V binding

0.86

0.0084

aPW Placental weight, PE placental efficiency, FDR P-value false discovery rate adjusted P-value, level of significance P < 0.05. bDEG Differentially expressed genes, N/A not available, log2FC log2 fold change, FDR P-value false discovery rate adjusted P-value, level of significance P < 0.10

Discussion

Fetal and utero-placental measurements

Placental weight was lower in the high PE group compared to the low PE group, but fetal weight was not different between high and low PE. These results agree with Krombeen and others [5] and confirm high PE placentas are smaller than low PE placentas, yet each grow a littermate of comparable size. There was a significant interaction between PE and sex for placental weight, which to the authors’ knowledge has not previously been reported in PE studies.

The lack of a difference in ISL between high PE and low PE was unexpected as the basis for high PE is a smaller but more efficient placenta that occupies less room within the uterus [4, 11]. Furthermore, a previous study reported high PE ISL were shorter than low PE ISL [12]. The conflicting results suggest ISL differs when PE is used as a selection tool as opposed to observing natural variations in PE. Alternatively, the folded placental trophoblast-endometrial epithelial bilayer width may be greater in feto-placental units with reduced placental size and comparable fetal growth (high PE). Vallet and Freking [15], reported greater fold widths were associated with the smallest pig fetuses in a litter and may increase PE via a larger surface area for exchange between the maternal and fetal circulations.

In addition to ISL, CRL, girth, heart weight, liver weight, brain weight, and ST weight were not significantly different between the high PE group and the low PE group. Crown-rump length and girth are highly correlated to fetal age [16] and weight [17, 18], and can be used to predict neonatal survival [19] and post-natal growth performance [20]. Considering these fetal measurements did not differ based on PE, the survival and postnatal growth performance of pigs grown on high PE placentas should not be negatively affected. Moreover, the absence of significant differences in fetal organ and tissue weights supports data from Krombeen and others [5], indicating any negative effects of a reduced placental size on fetal growth are not evident by term.

Placental VD and endometrial VD were also not significantly different between high PE and low PE. High PE in more prolific breeds has been attributed to increased placental VD during late gestation [11]; however, the role VD plays in extremes of PE within production breeds is less clear. Vonnahme and Ford [12] reported there was no additional increase in placental or endometrial VD to account for the increased efficiency of high PE placentas on day 90 of gestation in Yorkshires, despite increased expression of a vascular growth factor and its associated receptors. Conversely, Krombeen and others [5] identified a positive relationship between placental VD and PE on day 110 of gestation in maternal line gilts. It is conceivable that increased vascular permeability or reduced placental-endometrial intercapillary distance contribute to high PE, as suggested by Vonnahme and Ford [12], and/or changes in VD occur after day 90, as suggested by Krombeen and others [5].

Differential gene expression and gene ontology

A total of 214 DEG were identified in the placenta and 0 DEG were identified in the endometrium. Since the placenta is conceptus derived, it is not surprising that a greater number of genes would be differentially expressed in the placenta than in the endometrium, which is maternal tissue. Of the genes expressed in the placenta, only 214 were differentially expressed or 1.06% of the transcriptome. The small percentage of DEG identified in this study could be attributed to the comparison of two naturally occurring states as opposed to two treatment groups, the gestational day evaluated, and/or the expression level measured.

Nonetheless, the functional classification analysis performed on the DEG in high PE compared to low PE placentas identified molecular functions (MF) and biological processes (BP) associated with the phenotype. The MF of catalytic activity, binding, and transporter activity had the greatest representation. Catalytic activity was a MF of more upregulated than downregulated genes.

Upregulated candidate genes with catalytic activity included cytochrome P450 family 4 subfamily F member 22 (CYP4F22), fructose-1, 6-bisphosphatase isoenzyme 2 (FBP2), and chymotrypsin-like elastase family member 1 (CELA1). The catalytic activity of the products of CYP4F22, recently identified as an ultra-long chain fatty acid omega hydroxylase [21], and FBP2, encoding the gluconeogenic enzyme fructose-1,-6,biphosphatase-2 [22], suggests the metabolic state of high PE and low PE placentas differs. The gene CELA1 encodes for an enzyme that degrades the protein elastin, a component of the vascular matrix. Data in mice indicate there is a positive association between CELA1 and angiogenesis [23, 24]. While the catalytic activity of CELA1 may have a role in vascularity, no differences in VD were detected between high PE and low PE placentas.

Catalytic activity was also a function of downregulated genes; six transmembrane epithelial antigen of prostate 1 (STEAP1), six transmembrane epithelial antigen prostate 2 (STEAP2), and sarcosine dehydrogenase (SARDH) were candidate genes with catalytic activity. The STEAP family of proteins function as metal reductases, enabling the transport of iron and copper across the plasma membrane, and superoxide synthases, generating superoxide [25, 26]. The gene SARDH encodes for a mitochondrial enzyme that catalyzes the conversion of sarcosine to glycine, a major amino acid involved in an array of BP [27]. The downregulation of these genes in high PE compared to low PE suggests metabolism differs by PE.

The MF GO term with the second greatest representation in the DEG was binding, with fewer upregulated than downregulated genes involved in binding. Serum amyloid A2 (SAA2), sphingosine kinase type 1 interacting protein (SPHKAP), and dickkopf-1 (DKK1), were the three most upregulated genes in high PE compared to low PE placentas and shared the MF of binding. These genes encode a major acute phase protein involved in the innate immune response [28], an A-kinase anchoring protein involved in second messenger intracellular signaling [29], and a glycoprotein that is an inhibitor of the Wnt signaling pathway [30], respectively. In general, it appears the binding activity of these gene products is relevant to cell signaling pathways. While the role of these genes products within the placenta requires further investigation, SAA2 and DKK1 have been implemented in lipid metabolism [31, 32] and angiogenesis [3335].

Embigin (EMB) and angiopoietin 1 (ANGPT1) were downregulated candidate genes with the MF binding. The gene product of EMB is a transmembrane glycoprotein required for the localization and function of MCT2, a plasma membrane transporter of pyruvate, lactate, and ketone bodies [36]. The downregulation of EMB in high PE placentas suggests a lack of transport of these substrates, but given that pyruvate and lactate are gluconeogenic precursors and FBP2, encoding a gluconeogenic enzyme, was upregulated, it is plausible that these substrates may be metabolized to produce glucose in high PE placentas. On the other hand, ANGPT1 belongs to a family of endothelial growth factors and is a glycoprotein that inhibits endothelial permeability [37]. Thus, the downregulation of ANGPT1 in high PE compared to low PE placentas indicates vascular permeability may be increased in high PE placentas.

Transporter activity was another MF of the DEG, with a greater percentage of upregulated than downregulated genes involved in transporter activity. The following are candidate genes for PE with transporter activity: solute carrier family 45 member 3 (SLC45A3), acid sensing ion channel subunit 1 (ASIC1), solute carrier family 4 member 7 (SLC4A7), and solute carrier family 23 member 2 (SLC23A2). The gene SLC45A3 encodes for a novel H+ sucrose symporter, suggested to also transport fructose and glucose [38, 39]. Verification of this function in the pig placenta is required, but it appears sugar transport is increased in high PE placentas. Acid sensing ion channel 1 isoform X2 (ASIC1) encodes for a proton-gated sodium ion channel that localizes to the plasma membrane and golgi apparatus. Expression is most common in neurons [40] and to the authors’ knowledge has not been reported in the pig placenta. Assuming a similar MF, upregulation of ASIC1 in high PE placentas indicates sodium ion transmembrane transport differs based on PE. The gene SLC4A7 encodes for a sodium bicarbonate (Na+: HCO3) cotransporter [41] and the gene SLC23A2 encodes for a sodium dependent ascorbate (vitamin C) co-transporter (2Na+: ascorbate) [42]. In rodent models, both transporters have been implemented in the control of vasodilation [41, 43]. Thus, the upregulation of these transporters in high PE placentas may alter placental vascular tone, but further research is required to elucidate the role of these transporters in the pig placenta.

Transporter activity was also a function of downregulated genes, like multidrug resistance associated protein 4 (MRP4) and potassium channel inwardly rectifying subfamily J member 2 (KCNJ2). The gene MRP4 encodes for an active transporter protein with a broad substrate specificity [44]. Interpretation of the significance of MRP4 downregulation in high PE placentas requires further research into the substrates of MRP4 in the pig placenta. The gene KCNJ2 encodes the inwardly rectifying potassium channel KIR2.1. Expression of KIR2.1 has been reported in the human placenta, but the tissue specific function is unknown [45, 46]. Disruption of the potassium channel in mice indicated KIR2.1 mediates vasodilation [47]. The downregulation of KCNJ2 in high PE placentas suggest potassium transport and potentially vasodilation may be altered by PE.

Interestingly, two MF terms, molecular function regulator and transcription regulator activity, were functions of mostly downregulated genes. Molecular function regulators modulate a gene products activity and are often enzyme regulators or channel regulators [48]. Downregulated genes associated with this term were mostly enzyme regulators modulating intracellular activity. Transcription regulator activity describes the function of controlling gene expression at the level of transcription [48]. Accordingly, downregulated genes with this function encoded for transcription factors. Thus, downregulated genes were involved in the control of gene expression and the activity of gene products, which is not surprising given that this study compares high PE to low PE placentas.

Functional classification of the DEG also identified BP associated with the phenotype. Metabolic process, biological regulation, and localization had the greatest representation among the DEG. The gene products of more upregulated than downregulated genes were involved in metabolic processes. This is as expected given the MF catalytic activity and transporter activity had the greatest representation among the upregulated genes, and these functions are often involved in metabolism. For instance, CYP4F22 encodes for an ultra-long chain fatty acid omega hydroxylase, an enzyme of fatty acid metabolism [21]. Similarly, the genes FBP2 and SLC45A3 encode for a gluconeogenic enzyme [22] and sucrose transporter [38, 39], respectively, both of which are involved in carbohydrate metabolism. Furthermore, the identification of metabolic processes as the BP with the greatest representation in the DEG indicates extremes of PE are related to metabolism.

Biological regulation was also a BP term of a large percentage of the DEG and is a broad term encompassing genes products that modulate part of a BP [48]. More downregulated than upregulated genes were involved in biological regulation. Biological regulation was a term of anterior gradient protein 2 (AGR2), the most downregulated gene in high PE compared to low PE placentas. The gene AGR2 encodes for a member of the protein disulfide isomerase family of endoplasmic reticulum proteins, which are essential to post-translational folding [49]. The protein has been implemented in epithelial barrier function and cell proliferation. Moreover, it has been suggested that AGR2 downregulation in sheep placentomes may serve as an adaptive placental mechanism to support fetal growth during stress by reducing the proliferative actions of AGR2 [50]. Whether AGR2 regulates a similar BP in the pig placenta remains to be determined. Other downregulated genes involved in biological regulation included ANGPT1, a regulator of vascular permeability [37], and several other genes encoding transcription factors.

While there were more downregulated genes involved in biological regulation, a significant percentage of upregulated genes were also biological regulators. Among these were probable cation-transporting ATPase 13A3 (ATP13A3) and solute carrier family 52 member 3 (SLC52A3). The gene ATP13A3 encodes a protein involved in calcium ion transmembrane transport, with evidence of polyamine transport in worms [51] and humans [52]. Polyamines perform numerous essential functions in mammalian physiology and are known regulators of placental growth and angiogenesis [53]. Thus, ATP13A3 may be a candidate gene for high PE, but additional research into the substrate specificity of ATP13A3 in the pig placenta is required. The gene SLC52A3 also encodes for a transporter, but with riboflavin (vitamin B2) specificity. Riboflavin is a regulator of metabolism via the active forms flavin mononucleotide and flavin adenine dinucleotide. Knockout of Slc52a3 in mice caused reduced riboflavin concentrations in pups resulting in death, with signs of hyperlipidemia and hypoglycemia [54]. Thus, the upregulation of SLC52A3 in high PE placentas likely regulates metabolism. Overall, the large percentage of DEG involved in biological regulation suggests extremes of PE may originate from differential regulation of several BP.

Localization was also a BP term of a large percentage of the DEG, with greater representation in upregulated than downregulated genes. The term describes the transportation or maintenance of a substance to a location [48]. The percentage of DEG involved in localization is reflective of the MF with greatest representation among the DEG (catalytic activity, binding, and transporter activity). For example, ASIC1 and SLC23A2 were upregulated and encode for a sodium ion transmembrane transporter [40] and an ascorbate transmembrane co-transporter [42], respectively. Similarly, the downregulated genes STEAP1 and STEAP2 encode for transmembrane proteins with metal reductase and superoxide synthase activity [25, 26]. Furthermore, the DEG involved in localization indicate micronutrient transport differs in extremes of PE.

The terms developmental process and biological adhesion were BP of mainly downregulated genes. The term developmental process describes BP involved in the progression of a living unit [48]. Downregulated genes with this BP mostly encoded transcription factors, suggesting extremes of PE are driven by differences in the control of development. Conversely, biological adhesion was a term of downregulated genes involved in cell adhesion. Among these were EMB and ANGPT1, encoding a transmembrane protein that localizes monocarboxylate transporters to the cell membrane [36] and an endothelial growth factor known to inhibit endothelial permeability [37], respectively. Thus, downregulated genes involved in biological adhesion may affect the nutrient transport capacity of high PE placentas.

A GO enrichment analysis was also performed on the DEG in high PE compared to low PE placentas. Only 4 enriched terms were identified, which was probably due to the limited number of DEG in the input list. The MF terms sodium-dependent multivitamin transmembrane transporter activity and nucleobase transmembrane transporter activity, and the BP term nucleobase transport were enriched, indicating vitamin transport and/or the transport of nitrogenous bases from one side of the membrane to the other are imperative to high PE. The enrichment of these terms also identified SLC23A2, mentioned previously, as a candidate gene for high PE. The fourth enriched term was the cellular component term extracellular exosome, which describes gene products that localize to vesicles that are released from cells into the extracellular region via exocytosis [48]. Exosomes are involved in cell to cell communication and modulate intercellular communication at the maternal-fetal interface in pigs [55]. Therefore, the enrichment of this term suggests exosomes play a role in the cell to cell communication required for the increased efficiency of high PE placentas. Additionally, extracellular exosome was a term of 20 upregulated genes in high PE compared to low PE placentas, validating the association with high PE.

Fetal/utero-placental measurements and differentially expressed genes

The correlation analysis performed between fetal/utero-placental measures and the DEG identified significant correlations between 8 of the DEG and placental weight or PE. GRINL1A complex locus 1 (GCOM1), gene 19,656 (LOC100739517), ATP-binding cassette sub-family G member 2 (ABCG2), and TOX high mobility group box family member 3 (TOX3) were positively correlated with placental weight and downregulated in high PE compared to low PE placentas. The MF of GCOM1 in pigs is unknown. The gene ABCG2 encodes for an active transporter that is expressed in the human placenta and transports xenobiotic compounds [56, 57]. Cholesterol activity was also a GO term of ABCG2. The protein product of TOX3 may be involved in chromatin remodeling, and the bending and unwinding of DNA [57]. Molecular function GO terms included chromatin binding, phosphoprotein binding, protein homodimerization activity, and estrogen response element binding. Conversely, ras-related protein rab-6B (RAB6B) was negatively correlated with placental weight and was upregulated in high PE compared to low PE placentas. Molecular function GO terms of RAB6B included GTP binding, GTPase activity, and myosin v binding. The protein encoded by RAB6B localizes to the golgi apparatus and may function in retrograde membrane traffic [58]. Although the function(s) of these genes within the pig placenta are largely unknown, the identification of strong correlations with placental weight, in combination with gene expression, indicates GCOM1, ABCG2, TOX3, and RAB6B may regulate the reduced placental size of high PE placentas.

Gene 12188 (LOC100156118), transmembrane protein 199 (TMEM199), and proto-cadherin beta 1 (PCDHB1) were positively correlated with PE and upregulated in high PE compared to low PE placentas. Gene 12188 encodes an uncharacterized protein in swine. The MF of TMEM199 in pigs in unknown, but the protein encoded by this gene in humans may be involved in golgi homeostasis [57]. The gene PCDHB1 was among the 10 most upregulated genes. The specific function of PCDHB1 is unknown, but PCDHB1 may be a calcium dependent cell to cell adhesion protein [57]. Given the strong positive correlations of gene 12188, TMEM199, and PCDHB1 with PE, further research is warranted to determine the specific functions these genes may have within the pig placenta.

Conclusion

Placental efficiency, quantified by the ratio of fetal weight to placental weight, was determined within maternal line gilt litters to compare expression profiles of high PE feto-placental units to low PE feto-placental units. Mean fetal weight was not significantly different between the high PE group and low PE group, but placental weight was significantly reduced in in the high PE group, verifying comparisons were of similarly sized pigs grown on different sized placentas. Likewise, the absence of significant differences in fetal measures indicated any negative effects of a reduced placental size on fetal growth were not evident by day 95 of gestation. The comparison of gene expression profiles in the placenta and adjacent endometrium of high PE and low PE feto-placental units identified 214 DEG in the placenta and no DEG in the endometrium, confirming that the placenta responds to the fetus.

Gene ontology functional classification analysis of the 103 upregulated and 111 downregulated genes identified common MF and BP. The MF with the greatest representation among the DEG were catalytic activity, binding, and transporter activity. The BP with the greatest representation among the DEG were metabolism, biological regulation, and localization. Further investigation into the candidate genes associated with these terms partially supported the hypothesis and suggested extremes of PE are differentially regulated, affecting components of placental transport capacity like nutrient transport and blood flow. Conversely, DEG with growth factor activity were minimal and alternative functions were identified, indicating the complexity of the relationship between placental and fetal weights.

Overall, the results of this study support the use of PE as a marker of placental function and provide new insights into compensatory mechanisms that enable comparable fetal growth despite a reduced placental size. In swine, PE may provide an opportunity to optimize reproductive performance by normalizing the reduced birth weights of larger litters and in turn increasing pre-weaning survival; however, further research is required to effectively incorporate PE into selection schemes. Identifying associations between the phenotype and genome may be useful. Additionally, a limitation of this study was the method used to define high and low PE. Using the highest and lowest PE value within a litter is both advantageous and restrictive as it ensures comparisons are only between the very most and least efficient units, but excludes feto-placental units with efficiencies that are closer to the average. Therefore, it is suggested that future research determine the most appropriate method of defining high and low PE. Lastly, the role extracellular exosomes play in PE and the impact of environmental effects on PE is also of interest.

Methods

Animal management

All procedures were approved by the West Virginia University Animal Care and Use Committee (WVU-ACUC; ACUC # 10–0505). Eight Camborough 23 gilts (experimental unit; N = 8), owned by the WVU Animal Science Farm (Morgantown, WV), were group housed in a hoop structure equipped with tunnel ventilation and inspected annually by WVU-ACUC. Gilts were monitored for estrous behavior beginning at 5 months of age. The first estrus was observed and recorded. Gilts were bred by artificial insemination 12 and 24 h after the onset of a second estrus (6–7 months of age) using Pig Improvement Company (PIC) 1025 pooled maternal line semen (Birchwood Genetics, West Manchester, OH). Assignments for gestational day 95 ovario-hysterectomies (113–114 average day of farrowing) were randomly assigned at the time of breeding. Following breeding, gilts remained in this structure and were group housed throughout gestation.

Surgical procedure

At least 2 days before surgery, gilts were moved to the Food Animal Research Facility at the WVU Animal Science Farm. Gilts were taken off of feed 12 h before surgery. On the morning of the day of surgery, gilts (205.63 ± 37.6 kg) were anesthetized via jugular venipuncture using ketamine (3 mg/kg) and xylazine (2 mg/kg). Atropine sulfate (0.05 mg/kg) was administered to reduce salivation and isoflurane was used to maintain anesthesia.

Gilts were placed in dorsal recumbency and a mid-ventral incision was made to expose the gravid uterus. An antimesometrial incision was then made to open the uterus and expose the feto-placental units. Two tags were attached to the umbilical cord of each feto-placental unit, identifying which uterine horn the fetus (observational unit) originated from and the location within that uterine horn. The umbilical cord was cut between the two tags to ensure a tag remained with the placenta and the fetus. All fetuses were removed, and then the uterus was removed from the dam and set aside for processing. Sodium pentobarbital (400 mg/mL, Sigma Aldrich, St. Louis, MO) was used to euthanize the gilt.

Fetal and utero-placental measurements

Fetal weight, CRL, and girth were recorded as fetuses were removed from the uterus. Fetal necropsies were performed to obtain heart weight, liver weight, brain weight, and left hindquarter ST weight. The uterus was opened along the antimesometrial side and laid flat. Boundaries of each placenta were identified and a sample (~ 6.5 cm2) of all tissue layers (placenta, endometrium, and myometrium) was collected from an area void of calcium deposits and representative of the entire placenta. Samples were placed in tissue cassettes and fixed in neutral buffered formalin for histological processing. Then each placenta was peeled away from the endometrium and weighed. At this time representative samples, as described by Krombeen and others [5], were taken from both the placenta and the adjacent endometrium, placed in 2.0 mL cryovials (filled to 1.8 mL), and snap frozen in liquid nitrogen for RNA extraction. Implantation site length for each placenta was measured in the empty uterus using avascular bands as boundaries. Placental efficiency was determined for each feto-placental unit by dividing fetal weight by placental weight.

Similar to Krombeen and others [5], tissue cross sections containing placenta, endometrium, and myometrium were fixed in formalin, dehydrated with graded ethanol and xylenes, perfused with molten paraffin, and embedded in paraffin molds. Five micrometer sections were fixed to glass slides. Two sections for each fetus were stained using periodic acid and Schiff’s reagent (Sigma Aldrich, St. Louis, MO). Along the placental-endometrial interface, two fields per section were visualized (Nikon Eclipse TE2000–5, Nikon Instruments Inc., Melville, NY) and captured (Retiga 2000R, Q Imaging, Surrey, BC, Canada; Q Capture, Quantitative Imaging Corporation, v2.90.1, Surrey, BC, Canada) for a total of four fields visualized. Northern eclipse v6.0 software (Empix Inc., North Tonawanda, NY) was used to analyze images. Placental and endometrial tissues were outlined separately, and then total number of vessels, total area of vessels, and total area selected were measured. Vascular density was determined by dividing the area of the vessels by the total area selected. Replicates were averaged.

To compare fetal and utero-placental measurements of high PE and low PE units, the feto-placental unit with the highest PE and the feto-placental unit with the lowest PE in each litter (n = 8) were selected, creating the high PE group (n = 8, PE range 3.96 to 7.84) and the low PE group (n = 8, PE range 1.92 to 3.19) used for analysis. Statistical analyses were conducted using JMP Pro version 12.2.0 (SAS Institute Inc., Cary, NC 1989–2007). A linear mixed effects model was used to analyze each dependent variable (placental wt, fetal wt, ISL, CRL, girth, heart wt, liver wt, brain wt, ST wt, placental VD, and endometrial VD), with PE, sex, and PE*sex as fixed effects, and a random effect to account for pigs nested within dam. A significance level of 0.05 was used for all statistical tests.

Differential gene expression and gene ontology

Endometrial and placental samples from the most efficient (2 female, 6 male) and least efficient (3 female, 5 male) feto-placental unit in each litter (n = 8) were processed at the Clemson University Genomics & Computational Laboratory (CU-GCL). Total RNA was extracted from the endometrial (n = 16) and placental (n = 16) samples in duplicate using a RNeasy Plus Universal Mini Kit (Qiagen, Valencia, CA) and all extractions were performed according to the manufacturer’s instructions. An aliquot of each sample was qualitatively assayed for purity using UV spectroscopy via the Nanodrop8000 (ThermoFisher Scientific, Waltham, MA) to determine the 260/280 and 260/230 ratios, respectively. RNA integrity was measured using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). All RNA purity ratios, were > 1.8 and all RNA integrity numbers (RIN) were approximately 6. Total RNA was quantitated with the Broad Range Assay in the Qubit (ThermoFisher Scientific, Waltham, MA).

Each sample was normalized to a standard input concentration of 2 μg for sequencing library preparation. Stranded mRNA sequencing libraries were prepared manually at the CU-GCL with the TruSeq Stranded mRNA kit (Illumina, San Diego, CA) following the manufacturer’s recommended procedures. Sequencing data was collected on the HiSeq2500 (Illumina, San Diego, CA) using v4.0 chemistry and 2x125bp paired-end reads. Post sequencing, raw sequence reads were transferred to Clemson University’s Palmetto Cluster for analysis.

Sequence reads were quality validated with the FastQC software [59], followed by read preprocessing to remove adapter and primer sequences with the Trimmomatic software [60]. Processed sequence reads were aligned to the v10.2 Sus Scrofus reference genome assembly [61] with the GSNAP read alignment tool [62]. Sorted and indexed. BAM files were prepared from the. SAM output of GSNAP using Samtools [63]. Uniquely mapped read abundance per gene was determined with the featureCounts software in reversely stranded mode [64], and the count data per sample was output and transformed to tabular format.

Relative pairwise changes in gene level expression were determined with the edgeR software package [65]. Transcriptome comparisons were made using a generalized linear model and pairwise comparisons were made to compare low efficiency versus high efficiency conditions in a tissue specific manner. Differentially expressed genes were determined and filtered for significance using the FDR of 0.10 [66, 67]. Gene level fold-change values were output in tabular format and genes abounding thresholds were listed as candidate genes.

Candidate gene lists containing upregulated and downregulated genes in the placenta were functionally classified using the Panther database version 14.0 [68, 69] and AmiGO 2 version 2.5.12 was utilized to search GO term definitions [48, 70, 71]. Hierarchical functional classification was used to categorize genes according to the activity of the gene product (molecular function) and the pathway or processes the gene product functions in (biological process) [70, 71]. Gene ontology slim terms were utilized to classify gene lists according to defined terms. Percentages equal the number of genes within the input list with that MF or BP divided by the total number of MF or BP in the input list (gene hits against total number of function hits or gene hits against total number of process hits). Candidate gene lists containing upregulated and downregulated genes in the placenta were also independently tested for statistical enrichment (FDR < 0.05) with the GOSeq software tool [72].

The subset of candidate genes listed in Table 2 were associated with the MF and BP with the greatest representation in the DEG and were selected based on log2FC (10 most upregulated or downregulated genes) and/or GO terms related to nutrient transport, angiogenic activity, or growth factor activity.

Fetal/utero-placental measurements and differentially expressed genes

To identify relationships between the following measurements: fetal weight, placental weight, PE, CRL, brain weight, ST weight, ISL, placental VD, and endometrial VD, and the 214 DEG in high PE compared to low PE placentas, Pearson’s correlation coefficient (r) was estimated using the cor function in R [73]. A FDR adjustment was applied to correct for multiple comparisons using the p.adjust function in R. A significance level of 0.05 was used to identify significant correlations.

Abbreviations

ABCB1

ATP-binding cassette subfamily B member 1 isoform X2

ABCG2

ATP-binding cassette sub-family G member 2

AGR2

Anterior gradient protein 2 homolog

AK2

Adenylate kinase 2

ASIC1

Acid sensing ion channel 1 isoform X2

ATP13A3

Probable cation-transporting ATPase 13A3

CADM4

Cell adhesion molecule 4

CELA1

Chymotrypsin-like elastase family member 1

COCH

Coagulation factor c homolog

CRL: 

Crown-rump length

CTSH

Pro-cathepsin H

CYP4F22

Cytochrome P450 family 4 subfamily f member 22

DEG: 

Differentially expressed genes

DKK1

Dickkopf-1

EFHD1

EF-hand domain-containing protein D1 isoform X2

EMB

Embigin

ENTPD6

Ectonucleoside triphosphate diphosphohydrolase 6

FBP2

Fructose-1, 6-bisphosphatase isoenzyme 2

FDR: 

False discovery rate adjusted p-value

GCA

Grancalcin

GCOM1

Myocardial zonula adherens protein

GLA

Alpha-galactosidase A

GO: 

Gene ontology

HSPE1

10 KDa heat shock protein mitochondrial

ISL: 

Implantation site length

ITIH3

Inter-alpha-trypsin inhibitor heavy chain H3

IUGR: 

Intrauterine growth restriction

KCNJ2

Potassium channel inwardly rectifying subfamily J member 2

KRTAP8–1

Keratin associated protein 8–1

LEP

Leptin

LOC102159690: 

Solute carrier family 23 member 2-like

log2FC: 

Log2 fold change

MCOLN3

Mucolipin 3 isoform X2

MGST3

Microsomal glutathione S-transferase 3

MORN5

Morn repeat containing protein 5 isoform X8

MRP4 or ABCC4

Multidrug resistance-associated protein 4-like

PCDHB1

Protocadherin beta-1 isoform X2

PDHB

Pyruvate dehydrogenase E1 component subunit beta

PE: 

Placental efficiency

PEBP

Phosphatidylethanolamine-binding protein 1

PI16

Peptidase inhibitor 16 precursor

PIC: 

Pig improvement company

PPA1

Inorganic pyrophosphatase

PW: 

Placental weight

RAB6B

Ras-related protein rab-6b

RIN: 

Rna integrity number

RPL15

Ribosomal protein L15

SAA2

Serum amyloid A2

SARDH

Sarcosine dehydrogenase

SLC23A2

Solute carrier family 23 member 2

SLC45A3

Solute carrier family 45 member 3

SLC4A7

Solute carrier family 4 member 7

SLC52A3

Solute carrier family 52 member 3

SLITRK5

SLIT and NTRK-like protein 5

SPHKAP

Sphingosine kinase type 1

ST: 

Semitendinosus

STEAP1

Six-transmembrane epithelial antigen of prostate 1

STEAP2

Six transmembrane epithelial antigen of the prostate 2

TMEM199

Transmembrane protein 199

TMEM72

Transmembrane protein 72-like isoform X2

TOX3

Tox high mobility group box family member 3

TXN

Thioredoxin

TXNDC8

Thioredoxin domain-containing protein 8

VD: 

Vascular density

WVU: 

West virginia university

Declarations

Acknowledgements

The authors acknowledge Barbara Blackmon and Jeanice Troutman at Clemson University’s Genomics & Computational Lab for their assistance in the lab.

Funding

The authors acknowledge the support of a National Institute of General Medical Sciences Grant: P20GM109094.This grant supported Dr. Shankar’s salary, Dr. Noorai’s salary, provided a supplement to Dr. Saski, and paid for infrastructure at the genomics laboratory.

Availability of data and materials

All data generated or analyzed during this study are included in this published article. Raw reads are publicly available via BioProject PRJNA490325, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490325.

Authors’ contributions

SK performed RNA extractions, conducted statistical analyses of the morphological measurements under the guidance of JS, interpreted all data, and wrote the first draft of the manuscript. VS, RN, and CS analyzed the DEG data. JS provided guidance on and/or ran statistical analyses. MW and TW conceived the research and performed surgeries for data/sample collection. TW was principal investigator. All authors have read and approved the final manuscript.

Ethics approval and consent to participate

All procedures were approved by the West Virginia University Animal Care and Use Committee (ACUC # 10–0505), providing written consent for the use of the gilts in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Animal and Veterinary Science, Clemson University, Clemson, SC 29634, USA
(2)
Center for Human Genetics, Clemson University, Greenwood, SC 29646, USA
(3)
Genomics and Bioinformatics Facility, Clemson University, Clemson, SC 29634, USA
(4)
Department of Plant and Environmental Sciences, Clemson University, Clemson, SC 29634, USA
(5)
Department of Statistics, Colorado State University, Fort Collins, CO 80523, USA
(6)
Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV 26506, USA

References

  1. Wilson ME, Ford SP. Comparative aspects of placental efficiency. Reprod Suppl. 2001;58:223–32.PubMedGoogle Scholar
  2. Hayward CE, Lean S, Sibley CP, Jones RL, Wareing M, Greenwood SL, Dilworth MR. Placental adaptation: what can we learn from birthweight: placental weight ratio? Front Physiol. 2016;7:28.PubMedPubMed CentralView ArticleGoogle Scholar
  3. 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 Physiol. 2008;586(18):4567–76.PubMedPubMed CentralView ArticleGoogle Scholar
  4. Wilson M, Biensen N, Youngs C, Ford S. 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(4):905–10.PubMedView ArticleGoogle Scholar
  5. Krombeen SK, Bridges WC, Wilson ME, Wilmoth TA. Factors contributing to the variation in placental efficiency on days 70, 90, and 110 of gestation in gilts. J Anim Sci. 2019;97(1):359–73.PubMedView ArticleGoogle Scholar
  6. Vallet JL, McNeel AK, Johnson GA, Bazer FW. Triennial reproduction symposium: limitations in uterine and conceptus physiology that lead to fetal losses. J Anim Sci. 2013;91(7):3030–40.PubMedView ArticleGoogle Scholar
  7. National Agricultural Statistics Service: Quarterly hogs and pigs. http://usda.mannlib.cornell.edu/MannUsda/viewDocumentinfo.do?-documentID=1086. Accessed 20 Sept 2018.
  8. Town S, Patterson J, Pereira C, Gourley G, Foxcroft G. Embryonic and fetal development in a commercial dam-line genotype. Anim Reprod Sci. 2005;85(3–4):301–16.PubMedView ArticleGoogle Scholar
  9. Vonnahme K, Wilson M, Foxcroft G, Ford S. Impacts on conceptus survival in a commercial swine herd. J Anim Sci. 2002;80(3):553–9.PubMedView ArticleGoogle Scholar
  10. Vernunft A, Maass M, Brüssow K. Placental characteristics of german landrace sows and their relationships to different fertility parameters. Czech J Anim Sci. 2018;63(9):339–46.View ArticleGoogle Scholar
  11. Biensen N, Wilson M, Ford S. 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–76.PubMedView ArticleGoogle Scholar
  12. Vonnahme K, Ford S. Placental vascular endothelial growth factor receptor system mRNA expression in pigs selected for placental efficiency. J Physiol. 2004;554(1):194–201.PubMedView ArticleGoogle Scholar
  13. Zhou Q, Fang M, Huang T, Li C, Yu M, Zhao S. Detection of differentially expressed genes between erhualian and large white placentas on day 75 and 90 of gestation. BMC Genomics. 2009;10(1):337.PubMedPubMed CentralView ArticleGoogle Scholar
  14. Kwon SG, Hwang JH, Park DH, Kim TW, Kang DG, Kang KH, Kim I, Park HC, Na C, Ha J, Kim WK. Identification of differentially expressed genes associated with litter size in berkshire pig placenta. PLoS One. 2016;11(4):e0153311.PubMedPubMed CentralView ArticleGoogle Scholar
  15. Vallet J, Freking B. Differences in placental structure during gestation associated with large and small pig fetuses. J Anim Sci. 2007;85(12):3267–75.PubMedView ArticleGoogle Scholar
  16. Ullrey D, Sprague J, Becker D, Miller E. Growth of the swine fetus. J Anim Sci. 1965;24(3):711–7.PubMedView ArticleGoogle Scholar
  17. Jang Y, Ma Y, Lindemann M. Intrauterine position affects fetal weight and crown-rump length throughout gestation. J Anim Sci. 2014;92(10):4400–6.PubMedView ArticleGoogle Scholar
  18. Groesbeck CN, Lawrence KR, Young MG, Goodband RD, DeRouchey JM, Tokach MD, Nelssen JL, Dritz SS. Using heart girth to determine weight in finishing pigs. Kans Agric Exp Station Res Ser. 2002;0(10):166–8.Google Scholar
  19. Baxter EM, Jarvis S, D’eath RB, Ross DW, Robson SK, Farish M, Nevison IM, Lawrence AB, Edwards SA. Investigating the behavioural and physiological indicators of neonatal survival in pigs. Theriogenology. 2008;69(6):773–83.PubMedView ArticleGoogle Scholar
  20. Douglas S, Edwards S, Kyriazakis I. Are all piglets born lightweight alike? Morphological measurements as predictors of postnatal performance. J Anim Sci. 2016;94(8):3510–8.PubMedView ArticleGoogle Scholar
  21. Ohno Y, Nakamichi S, Ohkuni A, Kamiyama N, Naoe A, Tsujimura H, Yokose U, Sugiura K, Ishikawa J, Akiyama M, Kihara A. Essential role of the cytochrome P450 CYP4F22 in the production of acylceramide, the key lipid for skin permeability barrier formation. Proc Natl Acad Sci U S A. 2015;112(25):7707–12.PubMedPubMed CentralView ArticleGoogle Scholar
  22. Li H, Wang J, Xu H, Xing R, Pan Y, Li W, Cui J, Zhang H, Lu Y. Decreased fructose-1, 6-bisphosphatase-2 expression promotes glycolysis and growth in gastric cancer cells. Mol Cancer. 2013;12(1):110.PubMedPubMed CentralView ArticleGoogle Scholar
  23. Liu S, Young SM, Varisco BM. Dynamic expression of chymotrypsin-like elastase 1 over the course of murine lung development. Am J Phys Lung Cell Mol Phys. 2014;306(12):L1104–16.Google Scholar
  24. Varisco BM, Liu S. Chymotrypsin like elastase-1 (CELA1) regulates pulmonary vascular morphogenesis [abstract]. Am J Respir Crit Care Med. 2013;187:A5923.Google Scholar
  25. Kim K, Mitra S, Wu G, Berka V, Song J, Yu Y, Poget S, Wang D, Tsai A, Zhou M. Six-transmembrane epithelial antigen of prostate 1 (STEAP1) has a single b heme and is capable of reducing metal ion complexes and oxygen. Biochemistry. 2016;55(48):6673–84.PubMedView ArticleGoogle Scholar
  26. Ohgami RS, Campagna DR, McDonald A, Fleming MD. The steap proteins are metalloreductases. Blood. 2006;108(4):1388–94.PubMedPubMed CentralView ArticleGoogle Scholar
  27. Porter DH, Cook RJ, Wagner C. Enzymatic properties of dimethylglycine dehydrogenase and sarcosine dehydrogenase from rat liver. Arch Biochem Biophys. 1985;243(2):396–407.PubMedView ArticleGoogle Scholar
  28. Sandri S, Borbely AU, Fernandes I, de Oliveira EM, Knebel FH, Ruano R, Zugaib M, Filippin-Monteiro F, Bevilacqua E, Campa A. Serum amyloid a in the placenta and its role in trophoblast invasion. PLoS One. 2014;9(3):e90881.PubMedPubMed CentralView ArticleGoogle Scholar
  29. Kovanich D, van der Heyden MG, Aye TT, van Veen TB, Heck AR, Scholten A. Sphingosine kinase interacting protein is an A-kinase anchoring protein specific for type I cAMP-dependent protein kinase. Chembiochem. 2010;11(7):963–71.PubMedView ArticleGoogle Scholar
  30. Peng S, Li J, Miao C, Jia L, Hu Z, Zhao P, Li J, Zhang Y, Chen Q, Duan E. Dickkopf-1 secreted by decidual cells promotes trophoblast cell invasion during murine placentation. Reproduction. 2008;135(3):367–75.PubMedView ArticleGoogle Scholar
  31. Poitou C, Divoux A, Faty A, Tordjman J, Hugol D, Aissat A, Keophiphath M, Henegar C, Commans S, Clément K. Role of serum amyloid a in adipocyte-macrophage cross talk and adipocyte cholesterol efflux. J Clin Endocrinol Metab. 2009;94(5):1810–7.PubMedView ArticleGoogle Scholar
  32. Strakovsky RS, Pan Y. A decrease in DKK1, a WNT inhibitor, contributes to placental lipid accumulation in an obesity-prone rat model. Biol Reprod. 2012;86(3):1–11.View ArticleGoogle Scholar
  33. Mullan RH, Bresnihan B, Golden-Mason L, Markham T, O'hara R, FitzGerald O, Veale DJ, Fearon U. Acute-phase serum amyloid a stimulation of angiogenesis, leukocyte recruitment, and matrix degradation in rheumatoid arthritis through an NF-κB–dependent signal transduction pathway. Arthritis Rheum. 2006;54(1):105–14.PubMedView ArticleGoogle Scholar
  34. Chen Y, Hu Y, Zhou T, Zhou KK, Mott R, Wu M, Boulton M, Lyons TJ, Gao G, Ma J. Activation of the wnt pathway plays a pathogenic role in diabetic retinopathy in humans and animal models. Am J Pathol. 2009;175(6):2676–85.PubMedPubMed CentralView ArticleGoogle Scholar
  35. Smadja DM, d'Audigier C, Weiswald LB, Badoual C, Dangles-Marie V, Mauge L, Evrard S, Laurendeau I, Lallemand F, Germain S, Grelac F, Dizier B, Vidaud M, Bieche I, Gaussem P. The wnt antagonist dickkopf-1 increases endothelial progenitor cell angiogenic potential. Arterioscler Thromb Vasc Biol. 2010;30(12):2544–52.PubMedView ArticleGoogle Scholar
  36. Halestrap AP. The SLC16 gene family–structure, role and regulation in health and disease. Mol Asp Med. 2013;34(2–3):337–49.View ArticleGoogle Scholar
  37. Ong SP, Ng ML, Chu JJH. Differential regulation of angiopoietin 1 and angiopoietin 2 during dengue virus infection of human umbilical vein endothelial cells: implications for endothelial hyperpermeability. Med Microbiol Immunol (Berl). 2013;202(6):437–52.View ArticleGoogle Scholar
  38. Bartolke R, Heinisch JJ, Wieczorek H, Vitavska O. Proton-associated sucrose transport of mammalian solute carrier family 45: an analysis in saccharomyces cerevisiae. Biochem J. 2014;464(2):193–201.PubMedView ArticleGoogle Scholar
  39. Vitavska O, Edemir B, Wieczorek H. Putative role of the H+/sucrose symporter SLC45A3 as an osmolyte transporter in the kidney. Pflugers Arch. 2016;468(8):1353–62.PubMedView ArticleGoogle Scholar
  40. Wemmie JA, Price MP, Welsh MJ. Acid-sensing ion channels: advances, questions and therapeutic opportunities. Trends Neurosci. 2006;29(10):578–86.PubMedView ArticleGoogle Scholar
  41. Boedtkjer E, Praetorius J, Matchkov VV, Stankevicius E, Mogensen S, Fuchtbauer AC, Simonsen U, Fuchtbauer EM, Aalkjaer C. Disruption of na+,HCO(3)(−) cotransporter NBCn1 (slc4a7) inhibits NO-mediated vasorelaxation, smooth muscle ca(2)(+) sensitivity, and hypertension development in mice. Circulation. 2011;124(17):1819–29.PubMedView ArticleGoogle Scholar
  42. Takanaga H, Mackenzie B, Hediger MA. Sodium-dependent ascorbic acid transporter family SLC23. Pflugers Arch. 2004;447(5):677–82.PubMedView ArticleGoogle Scholar
  43. Richter H, Camm E, Modi B, Naeem F, Cross C, Cindrova-Davies T, Spasic-Boskovic O, Dunster C, Mudway I, Kelly F. Ascorbate prevents placental oxidative stress and enhances birth weight in hypoxic pregnancy in rats. J Physiol. 2012;590(6):1377–87.PubMedPubMed CentralView ArticleGoogle Scholar
  44. Russel FG, Koenderink JB, Masereeuw R. Multidrug resistance protein 4 (MRP4/ABCC4): a versatile efflux transporter for drugs and signalling molecules. Trends Pharmacol Sci. 2008;29(4):200–7.PubMedView ArticleGoogle Scholar
  45. Mylona P. Clarson l, Greenwood SL, Sibley CP. Expression of the Kir2.1 (inwardly rectifying potassium channel) gene in the human placenta and in cultured cytotrophoblast cells at different stages of differentiation. Mol Hum Reprod. 1998;4(2):195–200.PubMedView ArticleGoogle Scholar
  46. Wareing M, Greenwood S. Potassium channels in the human fetoplacental vasculature. Placenta. 2011;32:S203–6.PubMedView ArticleGoogle Scholar
  47. Zaritsky JJ, Eckman DM, Wellman GC, Nelson MT, Schwarz TL. Targeted disruption of Kir2.1 and Kir2.2 genes reveals the essential role of the inwardly rectifying K(+) current in K(+)-mediated vasodilation. Circ Res. 2000;87(2):160–6.PubMedView ArticleGoogle Scholar
  48. Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S. AmiGO hub, web presence working group. AmiGO: online access to ontology and annotation data. Bioinformatics. 2009;25(2):288–9.PubMedView ArticleGoogle Scholar
  49. Galligan JJ, Petersen DR. The human protein disulfide isomerase gene family. Hum Genomics. 2012;6(1):6.PubMedPubMed CentralView ArticleGoogle Scholar
  50. Del Curto HM. Anterior gradient homolog 2 and its potential roles in ovine pregnancy. Master's thesis, Texas A & M University. 2014. Available electronically from http://hdl.handle.net/1969.1/153372.
  51. Heinick A, Urban K, Roth S, Spies D, Nunes F. Phanstiel O,4th, Liebau E, Luersen K. Caenorhabditis elegans P5B-type ATPase CATP-5 operates in polyamine transport and is crucial for norspermidine-mediated suppression of RNA interference. FASEB J. 2010;24(1):206–17.PubMedPubMed CentralView ArticleGoogle Scholar
  52. Madan M, Patel A, Skruber K, Geerts D, Altomare DA, Phanstiel O IV. ATP13A3 and caveolin-1 as potential biomarkers for difluoromethylornithine-based therapies in pancreatic cancers. Am J Cancer Res. 2016;6(6):1231.PubMedPubMed CentralGoogle Scholar
  53. Wu G, Bazer FW, Hu J, Johnson GA, Spencer TE. Polyamine synthesis from proline in the developing porcine placenta. Biol Reprod. 2005;72(4):842–50.PubMedView ArticleGoogle Scholar
  54. Yoshimatsu H, Yonezawa A, Yamanishi K, Yao Y, Sugano K, Nakagawa S, Imai S, Omura T, Nakagawa T, Yano I. Disruption of Slc52a3 gene causes neonatal lethality with riboflavin deficiency in mice. Sci Rep. 2016;6:27557.PubMedPubMed CentralView ArticleGoogle Scholar
  55. Bidarimath M, Khalaj K, Kridli RT, Kan FW, Koti M, Tayade C. Extracellular vesicle mediated intercellular communication at the porcine maternal-fetal interface: a new paradigm for conceptus-endometrial cross-talk. Sci Rep. 2017;7:40476.PubMedPubMed CentralView ArticleGoogle Scholar
  56. Gedeon C, Anger G, Piquette-Miller M, Koren G. Breast cancer resistance protein: mediating the trans-placental transfer of glyburide across the human placenta. Placenta. 2008;29(1):39–43.PubMedView ArticleGoogle Scholar
  57. Stelzer G, Rosen R, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Stein T, Nudel R, Lieder I, Mazor Y, Kaplan S, Dahary D, Warshawsky D, Guan-Golan Y, Kohn A, Rappaport N, Safran M, Lancet D. The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. 2016;54(1) 1.30.1,1.30.33.Google Scholar
  58. Opdam FJ, Echard A, Croes HJ, van den Hurk JA, van de Vorstenbosch RA, Ginsel LA, Goud B, Fransen JA. The small GTPase Rab6B, a novel Rab6 subfamily member, is cell-type specifically expressed and localised to the golgi apparatus. J Cell Sci. 2000;113(15):2725–35.PubMedGoogle Scholar
  59. Babraham Bioinformatics: FastQC. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 20 Aug 2018.
  60. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics. 2014;30(15):2114–20.PubMedPubMed CentralView ArticleGoogle Scholar
  61. Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491(7424):393–8.PubMedPubMed CentralView ArticleGoogle Scholar
  62. Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005;21(9):1859–75.PubMedPubMed CentralView ArticleGoogle Scholar
  63. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.PubMedPubMed CentralView ArticleGoogle Scholar
  64. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2013;30(7):923–30.PubMedView ArticleGoogle Scholar
  65. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.PubMedView ArticleGoogle Scholar
  66. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57(1):289–300.Google Scholar
  67. Park BS, Mori M. Balancing false discovery and false negative rates in selection of differentially expressed genes in microarrays. Open Access Bioinformatics. 2010;2:1–9.View ArticleGoogle Scholar
  68. Mi H, Muruganujan A, Casagrande JT, Thomas PD. Large-scale gene function analysis with the PANTHER classification system. Nat Protoc. 2013;8:1551–66.PubMedView ArticleGoogle Scholar
  69. Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, Thomas PD. PANTHER version 11: expanded annotation data from gene ontology and reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2017;45(D1):D183–9.PubMedPubMed CentralView ArticleGoogle Scholar
  70. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.PubMedPubMed CentralView ArticleGoogle Scholar
  71. The Gene Ontology Consortium. Expansion of the gene ontology knowledgebase and resources. Nucleic Acids Res. 2017;45(D1):D331–8.View ArticleGoogle Scholar
  72. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.PubMedPubMed CentralView ArticleGoogle Scholar
  73. R: The R project for Statistical Computing. https://www.R-project.org/. Accessed 20 Aug 2018.

Copyright

© The Author(s). 2019

Advertisement