Maternal Pre-Pregnancy Obesity Is Associated with Altered Placental Transcriptome

Abstract

Maternal obesity has a major impact on pregnancy outcomes. There is growing evidence that maternal obesity has a negative influence on placental development and function, thereby adversely influencing offspring programming and health outcomes. However, the molecular mechanisms underlying these processes are poorly understood. We analysed ten term placenta’s whole transcriptomes in obese (n = 5) and normal weight women (n = 5), using the Affymetrix microarray platform. Analyses of expression data were carried out using non-parametric methods. Hierarchical clustering and principal component analysis showed a clear distinction in placental transcriptome between obese and normal weight women. We identified 72 differentially regulated genes, with most being down-regulated in obesity (n = 61). Functional analyses of the targets using DAVID and IPA confirm the dysregulation of previously identified processes and pathways in the placenta from obese women, including inflammation and immune responses, lipid metabolism, cancer pathways, and angiogenesis. In addition, we detected new molecular aspects of obesity-derived effects on the placenta, involving the glucocorticoid receptor signalling pathway and dysregulation of several genes including CCL2, FSTL3, IGFBP1, MMP12, PRG2, PRL, QSOX1, SERPINE2 and TAC3. Our global gene expression profiling approach demonstrates that maternal obesity creates a unique in utero environment that impairs the placental transcriptome.

About the Authors:

Signe Altmäe

* E-mail: [email protected] (SA); [email protected] (CC)

Affiliations Department of Women’s and Children’s Health, Division of Obstetrics and Gynecology, Karolinska Institutet, Stockholm, Sweden, Centre of Excellence for Paediatric Research EURISTIKOS and Department of Paediatrics, School of Medicine, University of Granada, Granada, Spain

ORCID http://orcid.org/0000-0002-0708-1865

Maria Teresa Segura

Affiliation: Centre of Excellence for Paediatric Research EURISTIKOS and Department of Paediatrics, School of Medicine, University of Granada, Granada, Spain

Francisco J. Esteban

Affiliation: Department of Experimental Biology, University of Jaén, Jaén, Spain

Sabine Bartel

Affiliation: Division of Experimental Asthma Research, Research Center Borstel, Leibniz-Center for Medicine and Biosciences, Member of the German Center for Lung Research (DZL), Borstel, Germany

Pilar Brandi

Affiliation: Centre of Excellence for Paediatric Research EURISTIKOS and Department of Paediatrics, School of Medicine, University of Granada, Granada, Spain

Martin Irmler

Affiliation: Institute of Experimental Genetics, Helmholtz Zentrum Muenchen, Neuherberg, Germany

Johannes Beckers

Affiliations Institute of Experimental Genetics, Helmholtz Zentrum Muenchen, Neuherberg, Germany, Technische Universität München, Chair of Experimental Genetics, Freising, Germany, German Center for Diabetes Research (DZD), Neuherberg, Germany

Hans Demmelmair

Affiliation: Ludwig-Maximilians-University of Munich, Dr. Hauner Children’s Hospital, University of Munich Medical Centre, Munich, Germany

Carmen López-Sabater

Affiliation: Department of Nutrition and Bromatology, School of Pharmacy, University of Barcelona, Spain

Berthold Koletzko

Affiliation: Ludwig-Maximilians-University of Munich, Dr. Hauner Children’s Hospital, University of Munich Medical Centre, Munich, Germany

Susanne Krauss-Etschmann

Affiliations Division of Experimental Asthma Research, Research Center Borstel, Leibniz-Center for Medicine and Biosciences, Member of the German Center for Lung Research (DZL), Borstel, Germany, Institute for Experimental Medicine, Christian-Albrechts-Universitaet zu Kiel, Kiel, Germany, Comprehensive Pneumology Center, Ludwig Maximilians University Hospital and Helmholtz Zentrum München, Großhadern, Germany

Cristina Campoy

* E-mail: [email protected] (SA); [email protected] (CC)

Affiliations Centre of Excellence for Paediatric Research EURISTIKOS and Department of Paediatrics, School of Medicine, University of Granada, Granada, Spain, Biohealth Institute of Granada, Granada, Spain

Introduction

Maternal obesity has increased between 70–100% over the last decade [1]. This has had consequences for different aspects of female reproduction, as maternal adiposity is related to an increased risk of the majority of maternal and fetal complications [1]. Maternal pre-gravid obesity has been associated with an increased risk of miscarriage, gestational diabetes, pre-eclampsia, Caesarean section, instrumental vaginal delivery, birth trauma, stillbirth, and pelvic floor damage [1,2]. Infants of overweight and obese mothers are at higher risk of being born large for gestational age and/or macrosomic, and for developing obesity, cardiovascular disease, and diabetes in adulthood [2].

Being at the interface between the maternal and fetal environment, the placenta plays a central role in how maternal obesity influences the programming of health outcomes in the offspring. It has been shown that the placental structure and function is shaped by obesity already at the early developmental stages and onward [3]. Furthermore, previous studies demonstrate that obesity during pregnancy promotes a maternal environment favouring increased inflammation, lipotoxicity, and oxidative stress in the placenta [4–6], which may in turn alter maternal endothelial function [7], trophoblast invasion and differentiation [7], vascular development and function [6], and placental nutrient transport [8]. The molecular mechanisms underlying these changes are poorly understood so far and therefore further research in this field is warranted. Knowing the molecular bases of these processes would provide valuable insights into the placental development and functions and help to identify molecular mechanisms that have both immediate and long lasting effects on fetal health.

With the recent evolution of ‘omics’ techniques, a number of genome-scale transcriptional studies on normal human placenta have been performed and a unique placental transcriptome has been identified [9–14]. However, the effects of maternal obesity on the human placental transcriptome have received little attention, and only a few studies have been published [3,5,6]. Saben et al. performed RNA-sequencing on term placentas from lean and obese women and found dysregulation of genes related to lipid metabolism, angiogenesis, hormone activity and inflammation in placentas from obese women. This indicates that the obese maternal environment may adversely affect placental development and function [6]. Two other studies focused on the first trimester placental transcriptome using microarray platforms [3,5]. Lassance et al. analysed the transcriptome of trophoblast cells of first trimester placenta exposed in vitro to insulin and obesity. They concluded that maternal obesity associated with insulin resistance programs the placental transcriptome towards refractoriness to insulin with potential adverse effects for placental structure and function [3]. Another study by Saben et al. identified several genes and signalling pathways in trophoblast cells following lipotoxic challenge, providing novel insights into the possible mechanisms underlying obesity-induced placental inflammation [5]. Thus, there is only one study to date that has directly investigated how maternal obesity influences the placental transcriptome.

We therefore set out: (i) to investigate the transcriptome profile of term placentae in obese women in comparison to normal weight women in order to provide further insights into the molecular effects of maternal obesity on placental environment, and (ii) to find molecular biomarkers/targets and candidate regulatory pathways dysregulated by maternal pre-pregnancy obesity.

Materials and Methods

Study subjects

Ten women participating in the PREOBE study were carefully selected for the current study. The PREOBE study (Role of Nutrition and Maternal Genetics on the Programming of Development of Fetal Adipose Tissue) is an observational cohort study on healthy normal weight, overweight, and obese women, as well as women who developed gestational diabetes. This study was performed at the Clinical University Hospital ‘San Cecilio’ and ‘Mother-Infant’ Hospital in Granada, Spain (the study has been registered at www.ClinicalTrials.gov, NLM identifier: NCT01634464 2012) [15]. The project was approved by the Research Bioethical Committee of the University of Granada, and all women signed written informed consent after receiving the full information by a member of the research team.

Five normal weight women, whose pre-pregnancy body mass index (BMI kg/m2) was between 18–25, and five obese women, whose BMI was ≥30, were selected for this study. All women were: with European descent, healthy, with age range between 25 and 35 years, had no pregnancy complications or gestational diabetes, had normal blood pressure throughout pregnancy, were non-smoking, consumed no alcohol or drugs, had a normal pregnancy with fetal development in accordance to the pregnancy week, fetal position was cephalic, had natural term delivery, and had a normal/healthy placenta based on visual inspection. During pregnancy, routine pregnancy evaluation was performed at weeks 24, 34, and at delivery (see Table 1 for clinical characteristics). Additionally, blood was collected for fatty acid analysis during pregnancy and at delivery (the protocol for fatty acid analysis is described in our previous publication [16].

[Figure omitted. See PDF.]

Table 1. Clinical characteristics of the participating women in the study and their perinatal outcomes.

https://doi.org/10.1371/journal.pone.0169223.t001

Collection of term placenta samples

Placenta samples were collected and weighed immediately after delivery by a well-trained expert. Disc samples containing both maternal and fetal tissue were collected from identical portions of the placental plate starting from the periphery of the maternal side in order to avoid any regional variations. Next, after removal of the decidua, a 0.5 x 0.5 x 0.5 cm (200 mg) sample was excised from the middle of the radius (distance between the insertion of the umbilical cord and the periphery) of each placenta. Each sample was then rinsed twice with saline solution (NaCl 0.9%), and immediately placed into sterile 1.5 ml microtubes containing RNAlater solution (Quiagen, Venlo, The Netherlands). All samples were stored under RNase free conditions at -80°C for later analysis at the Comprehensive Pneumology Center Munich.

Total RNA isolation and microarray hybridisation

For microarray and real-time PCR analysis, total RNA was isolated from the placental samples by using the miRNeasy mini kit (Quiagen) according to the manufacturer’s instructions. RNA concentration was calculated with a Nanodrop-1000 spectrophotometer (NanoDrop Fluorometer, Thermo Scientific, Wilmington, DE, USA). Only the highest quality RNA (260/280 ratio > 1.8, no degradation as detected by RNA agarose gel) was used for microarray analysis. 70% of the RNA samples were further analysed with Bioanalyzer (A2100 Bioanalyzer, Agilent Technologies). An RNA integrity value (RIN) of ≥7.0 was considered acceptable.

Thirty ng of total RNA was amplified using the Ovation PicoSL WTA System V2 in combination with the Encore Biotin Module (NuGEN Technologies, Inc, San Carlos, CA, USA). Amplified cDNA was hybridized on Affymetrix human Gene ST 2.0 arrays containing about 50,000 probe sets (Affymetrix, Santa Clara, CA, USA). Staining (Fluidics script FS450_00002) and scanning was done according to the Affymetrix expression protocol including minor modifications as suggested in the Encore Biotin protocol (NuGEN Technologies, Inc).

Microarray data analysis

Pre-processing and differential gene expression.

Expression console (v.1.2, Affymetrix) was used for quality control and to obtain annotated normalized RMA gene-level data (standard settings including median polish and sketch-quantile normalisation).

Data analyses were performed using the R-statistical software system (Free Software Foundation, Boston, USA). Gene expression profiles were determined by comparing the obese and normal weight groups (2 by 2 comparisons) by means of the rank product non-parametric test in the Bioconductor RankProd package (www.bioconductor.org/packages/devel/bioc/html/RankProd.html). Due to the limited number of samples, a non-parametric statistical test was used as a rough filter to narrow down the list of the most relevant genes. The statistics calculated is equivalent to the geometric mean rank and it is less sensitive to outliers. Additionally, the proportion of false positives (PFP) correction was used to control errors in multiple tests, as it effectively controls the accumulation of false-positives relative to the total number of positive results [17]. A PFP of <0.05 was considered statistically significant.

Our primary microarray data are available in the public database ArrayExpress (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-4541. Two samples did not pass the stringent Model Based Quality Control Assessment of Affymetrix GeneChips, which was carried out using the corresponding affyPLM package (https://bioconductor.org/packages/release/bioc/html/affyPLM.html). Briefly, the Normalized Unscaled Standard Errors (NUSE) function was run over the data. This process accounts for differences in variability between genes based on standard error (SE) estimates obtained for each gene on each array; any given array with elevated SE relative to the other arrays is considered of lower quality. The final analyses were performed in total on 8 samples (4 in each group).

Sample clustering and principal component analysis.

In order to validate the above gene selection with a non-parametric method, a principal component analysis (PCA) and a hierarchical clustering were performed using MeV 4.2.02 software (www.tm4.org) [18]. In PCA analysis, a three-dimensional scatter plot was produced in order to visualize the differences between the sample sets based on each sample’s gene expression profile. The method used to calculate the distance was the covariance and the number of probes that differentially expressed genes. In the hierarchical clustering, the data were Z-normalized by gene and the Euclidean distance was selected as the similarity to cluster expression profiles. The method used to perform the linkage was a complete-linkage hierarchical clustering algorithm.

Enrichment analysis of the results.

Enrichment analysis of differentially regulated genes was explored by using the Database for Annotation, Visualization and Integrated Discovery (DAVID, v. 6.7) [19] and gProfiler [20], and Ingenuity Pathways Analysis (IPA®) (Ingenuity® Systems, www.ingenuity.com, release date 09.12.2014). DAVID searches blocks of functionally related genes according to different criteria such as the Gene Ontology (GO) terms as biological process, cellular component and molecular function. GO FAT search was used in order to filter the broadest terms so that they do not overshadow the more specific terms (david.abcc.ncifcrf.gov). Further enrichment analyses, including canonical pathways and networks of dysregulated genes were analysed with IPA. A multiple testing correction, false discovery rate (FDR) was applied, and FDR<5.0% and a p-value <0.05 were considered statistically significant.

Microarray validation by real-time PCR

In order to validate the microarray results, 500ng of the same RNA samples were reverse-transcribed using the Quantitect Reverse Transcription kit according to the manufacturer’s recommendations (Qiagen, Venlo, Netherlands). The obtained cDNA was diluted 10x in PCR-grade water for qRT-PCR analysis.

The genes AREG, CCL2, FSTL3, IGFBP1 and MMP12 were used for microarray validation. Specific forward and reverse primer pairs were designed for each gene and are listed in Table A in S1 File. Real-time qPCR was performed on the Roche Light Cycler II platform with the Roche Light Cycler 480 SYBR Green I Master mix (Roche, Mannheim, Germany). The final products were verified by melting curve analysis by using the provided software (Roche, Mannheim, Germany). Data are presented as dCp values (Cptargetgene-Cpreferencegene), where higher dCp values represent lower expression; while graphs are represented on reversed Y-axes. B2M and YWHAZ where used as reference genes. We tested 8 different housekeeping genes that have been used in studies of human placenta (by [21] and others) for our study: SDHA, B2M, ACTB, BRD1, KCTD2, HPRT, YWHAZ and TBP. B2M and YWHAZ genes were the most suitable for our study as they showed the lowest variation among subjects and groups. The geometric mean of B2M and YWHAZ Ct values was used for normalisation. The analysis of gene expression differences between the study groups was performed using the Mann-Whitney U-test. P-values of <0.05 were considered statistically significant.

Results

Characteristics of participating subjects

The perinatal outcomes between normal weight and obese women did not differ significantly. However, the mean weight of the placenta, as well as the newborn’s birth weight and length, head circumference, waist circumference, and BMI were somewhat bigger/higher in obese women than normal weight women (Table 1).

The analysis of plasma fatty acids at pregnancy weeks 24 and 34 demonstrated significant differences in palmitic acid, oleic acid and arachidonic acid concentrations between the study groups (fatty acid values are expressed as percentages by weight, wt %. See Table B in S1 File). At pregnancy weeks 24 and 34 palmitic acid was significantly higher in obese women when compared to normal weight women (33.7±0.6 vs. 30.5±1.0, p = 0.016; and 34.3±0.7 vs. 32.2±1.1, p = 0.032). Oleic acid was significantly lower (10.7±1.0 vs. 13.0±1.8, p = 0.032) and arachidonic acid was significantly higher (9.7±0.2 vs. 7.9±1.4, p = 0.016) in plasma of obese women when compared to normal weight women at pregnancy week 34.

Cluster analysis of microarray data

The PCA analysis demonstrated a very clear separation of placental gene expression profiles between obese and normal weight women (Fig 1A). In addition, hierarchical clustering was applied to the microarray data and a similar pattern was obtained, where two study groups clustered into two different groups (Fig 1B).

[Figure omitted. See PDF.]

Fig 1.

(A) Principal component analysis (PCA) of term placental gene expression profiles in obese (Ob_ N°) and normal weight women (N_N°). (B) Cluster analysis of dysregulated genes in term placentas from obese (Ob_ N°) vs. normal weight women (N_N°). Red represents genes with high expression levels and green represents genes with low expression levels.

https://doi.org/10.1371/journal.pone.0169223.g001

Differential gene expression analysis of term placenta from obese vs. normal weight women

A total of 72 differentially regulated genes were identified. Eleven transcripts were up-regulated and 61 genes were down-regulated in term placentas of obese vs. normal weight women (see Table 2). Interestingly, the majority of the differentially regulated genes were down-regulated (84.7%), while only 15.3% of the genes were up-regulated in the placentas of obese women.

[Figure omitted. See PDF.]

Table 2. List of dysregulated genes in term placentas from obese women vs. normal weight women.

The genes identified in previous transcriptome studies in human placentas in health and disease are indicated.

https://doi.org/10.1371/journal.pone.0169223.t002

GO analysis of differentially regulated genes revealed that a significant proportion of the genes in the placentae of obese women vs. normal weight women were involved in carbohydrate binding (11.5%, p = 0.009), and more specifically in polysaccharide (9.8%, p = 0.001) and heparin binding (8.2%, p = 0.003). A large proportion of the differentially regulated genes were located in the extracellular region (49.2%, p<0.0001).

To obtain insight into relevant biological processes we used IPA analysis, and found that differentially regulated placental gene expression in obese women involved pathways concerning crosstalk between dendritic cells and natural killer cells, hepatic fibrosis/hepatic stellate cell activation, UPD-D-xylose and UPD-D-glucuronate biosynthesis, inhibition of matrix metalloproteases, atherosclerosis signalling among several others (Fig 2).

[Figure omitted. See PDF.]

Fig 2. The Circos plot represents significantly enriched pathways associated with regulated genes in the placentae of obese women vs. normal weight women, detected using the Ingenuity Pathway Analysis library of canonical pathways.

Outside the circle dysregulated genes and IPA pathways together with biological process are indicated. All genes are down-regulated (highlighted in green), except for HLA-DRB1 and CCL2 genes that are up-regulated (highlighted in red).

https://doi.org/10.1371/journal.pone.0169223.g002

Analysis of the molecular relationships between differentially expressed genes showed two complex networks where the majority of signals were mediated through ERK, IgG, NFKB, and MAPK complexes (Fig 3A), and through TGFB, MYC and TP53 complexes (Fig 3B). The highest scoring IPA network of the molecular relationships between differentially expressed genes revealed the involvement of the molecules in top diseases and functions such as embryonic development, organismal development, and cancer, where different genes including: AREG, CCL2, FN1, HLA-DRB1, IGFBP1, IL1R2, MMP12, PRL, and TNFSF10 were intertwining (IPA score of 46; Fig 3A). The second largest network of genes were involved in cellular movement, haematological system development and function, as well as immune cell trafficking, where different genes such as those for ADAM19, AOC1, IL2RB, PLAC8, SERPINE2, SP1 and others seem to play important roles (IPA score of 23; Fig 3B).

[Figure omitted. See PDF.]

Fig 3. Shows the two highest scoring networks in IPA derived from the analysis of dysregulated genes in placentas of obese women vs. normal weight women.

(A) Genes involved in embryonic development, organismal development and cancer. (B) Dysregulated genes involved in cellular movement, haematological system development and function, and immune cell trafficking. The intensity of the node colour indicates the degree of up- (red) or down- (green) regulation of gene expression. A white node represents a gene that is not part of our dataset, but is incorporated into the network through relationships with other genes. Nodes are displayed using various shapes that represent the functional class of the gene product, and a biological relationship between two nodes is represented as a line. Detailed information about the figure symbols can be found at www.ingenuity.com.

https://doi.org/10.1371/journal.pone.0169223.g003

The IPA Upstream Regulator Analysis identified two potential upstream regulators with a significant activation score (z-score ≥|2|) among the dysregulated genes in placentas of obese vs. normal weight women. IgG was predicted to be an activated upstream regulator (z-score = 2.0, p = 0.0004), with CCL2, CRABP2, GPRC5A, LGALS3 as known target molecules. The second potential upstream regulator was ERBB2, being inhibited (z-score = -2.0, p = 0.002), with genes ADAM19, CGB8, FN1, and FSTL3 as downstream targets.

Microarray validation

Real-time PCRs using genes AREG, CCL2, FSTL3, IGFBP1 and MMP12 confirmed the array results. CCL2 was significantly up-regulated, and IGFBP1 and MMP12 were significantly down-regulated (Fig 4). AREG and FSTL3 genes expression profiles between obese and normal weight women did not reach statistical significance (p>0.05), nevertheless we observed a trend towards a similar dysregulation compared to the microarray.

[Figure omitted. See PDF.]

Fig 4. Microarray validation by real-time PCR.

The expression levels of genes AREG, CCL2, FSTL3, IGFBP1 and MMP12 from our gene array analysis in comparison with real-time PCR gene expression levels are presented. QRT-PCR data are shown as dCp values (Cptargetgene-Cpreferencegene), where higher dCp values represent lower expression and therefore, graphs are represented with reversed Y-axis. Mann-Whitney U-test *p <0.05, n = 5 samples/group.

https://doi.org/10.1371/journal.pone.0169223.g004

Discussion

We report the genome-wide transcriptome analysis of the effects of maternal obesity on the term placenta, in the absence of gestational diabetes and other complications. Our study results demonstrate that maternal pre-pregnancy obesity has adverse effects on the placental transcriptome, where previously characterised molecules and molecular pathways involved in placental development and function [14] were dysregulated. In addition, several new dysregulated molecules and signalling pathways were identified.

It is now well established that maternal obesity is associated with a pro-inflammatory milieu during pregnancy and in the placenta, [5,6]. In line with that, obese women in the present study had significantly higher plasma levels of arachidonic acid (which is involved in inflammatory processes) at pregnancy week 34. It has been demonstrated in mice that obesity during pregnancy disrupts inflammation through increased macrophage activation and elevation of cytokine gene expression, which provides potential links between placental inflammation and the programming of offspring disease by maternal obesity [29]. Furthermore, the majority of the identified dysregulated genes in our study including CCL2, HLA-DRB1, IL1R2, IL2RB, TNFSF10, FN1, and MMP12 are involved in inflammation and immune responses, and the second biggest network of the dysregulated genes being involved in immune cell trafficking. We also identified immunoglobulin G (IgG) as an upstream activator of several dysregulated genes (CCL2, CRABP2, GPRC5A, LGALS3) in placentas from obese women. IgG is the main type of antibody found in blood and extracellular fluid that controls infection in body tissue by mediating pro- and anti-inflammatory activities [30].

Lipid metabolism is another process that has shown to be altered in placentas from obese women [6,7]. Our obese women demonstrated different palmitic acid and oleic acid levels in the plasma at pregnancy weeks 24 and 34 when compared to normal weight women, supporting altered lipid metabolism. Furthermore, functional analyses of our microarray data identified lipid and cholesterol metabolism, specifically the LXR/RXR activation pathway to be dysregulated among placentas from obese women. The importance of liver X receptors (LXRs) in physiological lipid and cholesterol metabolism suggests that they are involved in the development of metabolic disorders such as hyperlipidaemia (lipotoxicity) and atherosclerosis [31]. Indeed, a recent study concluded that maternal obesity leads to a lipotoxic placental environment that is associated with decreased regulators of angiogenesis and increased markers of inflammation and oxidative stress [6].

Another interesting finding in our study was that all analyses indicated that obesity influences cancer pathways. The detection of cancer pathways among placental gene expression is not surprising, as many proliferative, invasive, and immune tolerance mechanisms that support normal pregnancy are also exploited by malignancies [12,32]. Furthermore, we identified erb-b2 receptor tyrosine kinase 2 (ERBB2) as an upstream inhibitor of several genes including ADAM19, CGB8, FN1, and FSTL3. ERBB2 is involved in the mitogen-activated protein kinase signalling pathway and its role in cancer development and evolution has been shown [33].

In line with the similarities in pregnancy and cancer pathways, we detected haematological system development and function in network analysis, and angiogenesis as important biological processes affected by the dysregulated genes in placentas from obese women. Angiogenesis is a crucial process for fetomaternal exchanges and placental development, and alterations in this are associated with different pregnancy-related pathologies [34]. This is in accordance with a previous study that has detected altered regulation of genes related to angiogenesis in placentas from obese women [6].

A new observation in our study was the dysregulation of glucocorticoid receptor signalling pathway in placentas from obese women, where expression of genes CCL2, IL1R2, PRL, and SLPI were altered. Obesity has been associated with reduced sensitivity to glucocorticoid feedback, an effect believed to be mediated via altered sensitivity to the glucocorticoid receptor [35]. Glucocorticoid receptor signalling is regulating genes controlling the development, metabolism and immune responses, and is also involved in major organ systems physiology and pathophysiology in the human body. Its important role during gestation for postnatal survival as well as during embryonic development has been shown in different animal studies [35].

Another important finding in the search of obesity-related effects on term placentas was the identification of several potential target molecules that could have important roles in understanding how the placenta affects fetal development and obesity-derived offspring’s future health problems. We detected dysregulation in several genes that have been shown to play a role in the function of normal placental development in humans, including AREG [36], CCL2 [37], FOSB [38], GKN1 [39], SERPINE2 [40], SLPI [41], and XIST [42]. In addition, several of the identified dysregulated genes, including CCL2, CGB8, FOSB, FSTL3, HSD11B, IGFBP1, PRL, RBP4, and TAC3 have been implicated in adverse pregnancy outcomes such as intrauterine growth restriction [23,25,26,28], large for gestational age [24], and recurrent miscarriage [43]. Furthermore, a set of the dysregulated genes has been associated with the risk of developing pre-eclampsia, including c12orf75 [23], DKK1 [44], FSTL3 [27], HTRA4 [45], IGFBP1 [27], LAIR2 [27], MMP12 [27], PAEP [27], and UPK1B [23]. Especially interesting and promising target molecules for obesity-derived implications in placenta could be CCL2, PRL, MMP12, TAC3, PRG2 and IGFBP1 that were the most dysregulated genes among our study group. Additionally, genes FSTL3, QSOX1 and SERPINE2 could serve as novel obesity-related biomarkers, as they have been shown to be uniquely enriched in the placenta [14], and we detected them as down-regulated among the placentas from obese women. All of these genes could serve as molecular biomarkers for potential progression towards metabolic syndrome in children that were born to overweight mothers.

A potential weakness of our and similar studies [5,10,11,22] is the small sample size. Nonetheless, cluster analyses showed a very clear separation between the groups in the present study. In addition, we have applied stringent non-parametric data analysis. Our study was strengthened by the well selected and characterised homogenous study groups, where healthy, young, non-smoking, Spanish women with no pregnancy complications and with natural delivery were enrolled.

Conclusions

The placental metabolic abnormalities resulting from the effects of maternal obesity (e.g. lipotoxicity) may be transmitted to the offspring via in utero programming and hence there could be far-reaching consequences for offspring health. Management of maternal body weight and/or manipulation of lipid metabolism using both lifestyle and pharmacological interventions may provide benefits to the obese women as well as to their offspring. In fact, it has been shown that maternal weight loss by pre-pregnancy bariatric surgery prevented transmission of obesity to offspring compared with children of obese mothers who did not undergo this surgery [46]. Further improvements/intervention programs in maternal weight control (preferably before pregnancy [47]) would improve perinatal metabolic outcomes.

Our findings provide a fundamental resource for better understanding the complex effects of maternal pre-pregnancy obesity on placental transcriptome. We have detected new molecular aspects of obesity-derived effects on placentas, where the glucocorticoid receptor signalling pathway and dysregulation of several genes including CCL2, FSTL3, IGFBP1, MMP12, PRG2, PRL, QSOX1, SERPINE2 and TAC3 might have important roles. Our study also confirms the dysregulation of previously identified important processes and pathways among placentas from obese women. This global gene expression profiling approach demonstrates and confirms that maternal obesity creates a unique in utero environment that impairs the placental transcriptome. Further elucidation of normal and aberrant placental ‘ome’ regulation will not only increase our understanding of the origins of a disease, but may also provide novel approaches for intervention.

Supporting Information

[Figure omitted. See PDF.]

S1 File.

Table A. Real-time PCR primer sequences used for microarray validation. Table B. Fatty acids in plasma of participating women throughout the pregnancy.

https://doi.org/10.1371/journal.pone.0169223.s001

(DOCX)

Acknowledgments

We thank participating women and PREOBE team. Special thanks to Anke Bettenbrock for excellent technical assistance and to Christine Delisle Nyström for the language revision.

Author Contributions

1. Conceptualization: SA MTS BK SK-E CC.

2. Data curation: SA FJE MI.

3. Formal analysis: SA FJE SB MI JB.

4. Funding acquisition: SA BK SK-E CC.

5. Investigation: SA MTS SB PB CL-S SK-E.

6. Methodology: MTS MI JB CL-S.

7. Project administration: BK SK-E CC.

8. Resources: JB HD CL-S BK SK-E CC.

9. Software: SA MI JB.

10. Supervision: SK-E CC.

11. Validation: SB HD SK-E.

12. Visualization: SA FJE.

13. Writing – original draft: SA MTS.

14. Writing – review & editing: SA MTS FJE SB PB MI JB HD CL-S BK SK-E CC.

Citation: Altmäe S, Segura MT, Esteban FJ, Bartel S, Brandi P, Irmler M, et al. (2017) Maternal Pre-Pregnancy Obesity Is Associated with Altered Placental Transcriptome. PLoS ONE 12(1): e0169223. https://doi.org/10.1371/journal.pone.0169223

References

1. Dixit A, Girling JC. Obesity and pregnancy. J Obstet Gynaecol. 2008;28: 14–23. pmid:18259892

2. Aune D, Saugstad OD, Henriksen T, Tonstad S. Maternal body mass index and the risk of fetal death, stillbirth, and infant death: a systematic review and meta-analysis. JAMA. 2014;311: 1536–46. pmid:24737366

3. Lassance L, Haghiac M, Leahy P, Basu S, Minium J, Zhou J, et al. Identification of Early Transcriptome Signatures in Placenta exposed to Insulin and Obesity. Am J Obstet Gynecol. Elsevier Ltd; 2015;

4. Malti N, Merzouk H, Merzouk S a., Loukidi B, Karaouzene N, Malti A, et al. Oxidative stress and maternal obesity: Feto-placental unit interaction. Placenta. Elsevier Ltd; 2014;35: 411–416.

5. Saben J, Zhong Y, Gomez-Acevedo H, Thakali KM, Borengasser SJ, Andres A, et al. Early growth response protein-1 mediates lipotoxicity-associated placental inflammation: role in maternal obesity. Am J Physiol Endocrinol Metab. 2013;305: E1–E14. pmid:23632636

6. Saben J, Lindsey F, Zhong Y, Thakali K, Badger TM, Andres A, et al. Maternal obesity is associated with a lipotoxic placental environment. Placenta. Elsevier Ltd; 2014;35: 171–177.

7. Jarvie E, Hauguel-de-Mouzon S, Nelson SM, Sattar N, Catalano PM, Freeman DJ. Lipotoxicity in obese pregnancy and its potential role in adverse pregnancy outcome and obesity in the offspring. Clin Sci (Lond). 2010;119: 123–129.

8. Higgins L, Greenwood SL, Wareing M, Sibley CP, Mills T a. Obesity and the placenta: A consideration of nutrient exchange mechanisms in relation to aberrant fetal growth. Placenta. Elsevier Ltd; 2011;32: 1–7.

9. Sitras V, Fenton C, Paulssen R, Vårtun Å, Acharya G. Differences in gene expression between first and third trimester human placenta: A microarray study. PLoS One. 2012;7: 1–7.

10. Sood R, Sood R, Zehnder JL, Zehnder JL, Druzin ML, Druzin ML, et al. Gene expression patterns in human placenta. Proc Natl Acad Sci U S A. 2006;103: 5478–83. pmid:16567644

11. Mikheev AM, Nabekura T, Kaddoumi A, Bammler TK, Govindarajan R, Hebert MF, et al. Profiling gene expression in human placentae of different gestational ages: an OPRU Network and UW SCOR Study. Reprod Sci. 2008;15: 866–877. pmid:19050320

12. Winn VD, Haimov-Kochman R, Paquet AC, Yang YJ, Madhusudhan MS, Gormley M, et al. Gene expression profiling of the human maternal-fetal interface reveals dramatic changes between midgestation and term. Endocrinology. 2007;148: 1059–1079. pmid:17170095

13. Sitras V, Paulssen RH, Grønaas H, Vårtun Å, Acharya G. Gene expression profile in labouring and non-labouring human placenta near term. Mol Hum Reprod. 2008;14: 61–65. pmid:18048457

14. Saben J, Zhong Y, McKelvey S, Dajani NK, Andres A, Badger TM, et al. A comprehensive analysis of the human placenta transcriptome. Placenta. Elsevier Ltd; 2014;35: 125–131.

15. Berglund SK, García-Valdés L, Torres-Espinola FJ, Segura MT, Martínez-Zaldívar C, Aguilar MJ, et al. Maternal, fetal and perinatal alterations associated with obesity, overweight and gestational diabetes: an observational cohort study (PREOBE). BMC Public Health. 2016; in press: 207.

16. Krauss-Etschmann S, Shadid R, Campoy C, Hoster E, Demmelmair H, Jimenez M, et al. Effects of fish-oil and folate supplementation of pregnant women on maternal and fetal plasma concentrations of docosahexaenoic acid and eicosapentaenoic acid: A European randomized multicenter trial. Am J Clin Nutr. 2007;85: 1392–1400. pmid:17490978

17. Fernando RL, Nettleton D, Southey BR, Dekkers JC, Rothschild MF, Soller M. Controlling the proportion of false positives in multiple dependent tests. Genetics. 2004/03/17. 2004;166: 611–619. pmid:15020448

18. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003/03/05. 2003;34: 374–378. pmid:12613259

19. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2008/11/27. 2009;37: 1–13. pmid:19033363

20. Reimand J, Arak T, Vilo J. G:Profiler—A web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res. 2011;39: 1–9.

21. Meller M, Vadachkoria S, Luthy DA, Williams MA. Evaluation of housekeeping genes in placental comparative expression studies. Placenta. 2005;26: 601–607. pmid:16085039

22. Zhang Y, Cui Y, Zhou Z, Sha J, Li Y, Liu J. Altered global gene expressions of human placentae subjected to assisted reproductive technology treatments. Placenta. 2010;31: 251–8. pmid:20116094

23. Sõber S, Reiman M, Kikas T, Rull K, Inno R, Vaas P, et al. Extensive shift in placental transcriptome profile in preeclampsia and placental origin of adverse pregnancy outcomes. Sci Rep. 2015;5: 13336. pmid:26268791

24. Ahlsson F, Akerud H, Schijven D, Olivier J, Sundstrom-Poromaa I. Gene Expression in Placentas From Nondiabetic Women Giving Birth to Large for Gestational Age Infants. Reprod Sci. 2015;

25. Madeleneau D, Buffat C, Mondon F, Grimault H, Rigourd V, Tsatsaris V, et al. Transcriptomic analysis of human placenta in intrauterine growth restriction. Pediatr Res. 2015; 1–9.

26. Okamoto A, Endo H, Kalionis B, Shinya M, Saito M, Nikaido T, et al. IGFBP1 and Follistatin-like 3 genes are significantly up-regulated in expression profiles of the IUGR placenta. Placenta. 2006;27: 317–21. pmid:16338475

27. Founds SA, Conley YP, Lyons-Weiler JF, Jeyabalan A, Hogge WA, Conrad KP. Altered global gene expression in first trimester placentas of women destined to develop preeclampsia. Placenta. 2009;30: 15–24. pmid:19027158

28. Whitehead CL, Walker SP, Ye L, Mendis S, Kaitu’u-Lino TJ, Lappas M, et al. Placental specific mRNA in the maternal circulation are globally dysregulated in pregnancies complicated by fetal growth restriction. J Clin Endocrinol Metab. 2013;98: 429–436.

29. Kim DW, Young SL, Grattan DR, Jasoni CL. Obesity during pregnancy disrupts placental morphology, cell proliferation, and inflammation in a sex-specific manner across gestation in the mouse. Biol Reprod. 2014;90: 130. pmid:24829026

30. Kaneko Y, Nimmerjahn F, Ravetch J V. Anti-inflammatory activity of immunoglobulin G resulting from Fc sialylation. Science. 2006;313: 670–673. pmid:16888140

31. Wang L, Schuster GU, Hultenby K, Zhang Q, Andersson S, Gustafsson J-A. Liver X receptors in the central nervous system: from lipid homeostasis to neuronal degeneration. Proc Natl Acad Sci U S A. 2002;99: 13878–83. pmid:12368482

32. Ferretti C, Bruni L, Dangles-Marie V, Pecking Pa, Bellet D. Molecular circuits shared by placental and cancer cells, and their implications in the proliferative, invasive and migratory capacities of trophoblasts. Hum Reprod Update. 2007;13: 121–141. pmid:17068222

33. Baselga J, Swain SM. Novel anticancer targets: revisiting ERBB2 and discovering ERBB3. Nat Rev Cancer. Nature Publishing Group; 2009;9: 463–475.

34. Alfaidy N, Hoffmann P, Boufettal H, Samouh N, Aboussaouira T, Benharouga M, et al. The multiple roles of EG-VEGF/PROK1 in normal and pathological placental angiogenesis. Biomed Res Int. 2014;2014.

35. Correia-Branco A, Keating E, Martel F. Maternal undernutrition and fetal developmental programming of obesity: the glucocorticoid connection. Reprod Sci. 2015;22: 138–45. pmid:25001018

36. Hung T-H, Chen S-F, Hsieh C-C, Hsu J-J, Li M-J, Yeh Y-L, et al. Tumor necrosis factor-alpha converting enzyme in the human placenta throughout gestation. Reprod Sci. 2008;15: 195–209. pmid:18276953

37. Amsalem H, Kwan M, Hazan A, Zhang J, Jones RL, Whittle W, et al. Identification of a novel neutrophil population: proangiogenic granulocytes in second-trimester human decidua. J Immunol. 2014;193: 3070–9. pmid:25135830

38. Renaud SJ, Kubota K, Rumi MAK, Soares MJ. The FOS transcription factor family differentially controls trophoblast migration and invasion. J Biol Chem. 2014;289: 5025–39. pmid:24379408

39. Fahlbusch FB, Ruebner M, Huebner H, Volkert G, Zuern C, Thiel F, et al. The tumor suppressor gastrokine-1 is expressed in placenta and contributes to the regulation of trophoblast migration. Placenta. 2013;34: 1027–35. pmid:23993393

40. Chern S-R, Li S-H, Chiu C-L, Chang H-H, Chen C-P, Tsuen Chen EI. Spatiotemporal expression of SERPINE2 in the human placenta and its role in extravillous trophoblast migration and invasion. Reprod Biol Endocrinol. 2011;9: 106. pmid:21806836

41. Klaffenbach D, Friedrich D, Strick R, Strissel PL, Beckmann MW, Rascher W, et al. Contribution of different placental cells to the expression and stimulation of antimicrobial proteins (AMPs). Placenta. 2011;32: 830–7. pmid:21899884

42. Moreira de Mello JC, de Araújo ESS, Stabellini R, Fraga AM, de Souza JES, Sumita DR, et al. Random X inactivation and extensive mosaicism in human placenta revealed by analysis of allele-specific gene expression along the X chromosome. PLoS One. 2010;5: e10947. pmid:20532033

43. Rull K, Nagirnaja L, Ulander V-M, Kelgo P, Margus T, Kaare M, et al. Chorionic gonadotropin beta-gene variants are associated with recurrent miscarriage in two European populations. J Clin Endocrinol Metab. 2008;93: 4697–706. pmid:18782867

44. Zhang Z, Li H, Zhang L, Jia L, Wang P. Differential expression of β-catenin and Dickkopf-1 in the third trimester placentas from normal and preeclamptic pregnancies: a comparative study. Reprod Biol Endocrinol. 2013;11: 17. pmid:23452984

45. Inagaki A, Nishizawa H, Ota S, Suzuki M, Inuzuka H, Miyamura H, et al. Upregulation of HtrA4 in the placentas of patients with severe pre-eclampsia. Placenta. 2012;33: 919–26. pmid:22964307

46. Kral JG, Biron S, Simard S, Hould F-S, Lebel S, Marceau S, et al. Large maternal weight loss from obesity surgery prevents transmission of obesity to children who were followed for 2 to 18 years. Pediatrics. 2006;118: e1644–9. pmid:17142494

47. Rönö K, Stach-Lempinen B, Klemetti MM, Kaaja RJ, Pöyhönen-Alho M, Eriksson JG, et al. Prevention of gestational diabetes through lifestyle intervention: study design and methods of a Finnish randomized controlled multicenter trial (RADIEL). BMC Pregnancy Childbirth. 2014;14: 70. pmid:24524674