Original Research ARTICLE
The effects of antenatal depression and antidepressant treatment on placental gene expression
- 1Department of Women's and Children's Health, Uppsala University, Uppsala, Sweden
- 2Department of Behavioral Physiology, University of Groningen, Groningen, Netherlands
- 3Department of Medicine, Centre for Gender Medicine, Karolinska Institute, Stockholm, Sweden
The effects of antenatal depression and antidepressant treatment during pregnancy on both mother and child are vigorously studied, but the underlying biology for these effects is largely unknown. The placenta plays a crucial role in the growth and development of the fetus. We performed a gene expression study on the fetal side of the placenta to investigate gene expression patterns in mothers with antenatal depression and in mothers using antidepressant treatment during pregnancy. Placental samples from mothers with normal pregnancies, from mothers with antenatal depression, and from mothers using antidepressants were collected. We performed a pilot microarray study to investigate alterations in the gene expression and selected several genes from the microarray for biological validation with qPCR in a larger sample. In mothers with antenatal depression 108 genes were differentially expressed, whereas 109 genes were differentially expressed in those using antidepressants. Validation of the microarray revealed more robust gene expression differences in the seven genes picked for confirmation in antidepressant-treated women than in depressed women. Among the genes that were validated ROCK2 and C12orf39 were differentially expressed in both depressed and antidepressant-treated women, whereas ROCK1, GCC2, KTN1, and DNM1L were only differentially expressed in the antidepressant-treated women. In conclusion, antenatal depression and antidepressant exposure during pregnancy are associated with altered gene expression in the placenta. Findings on those genes picked for validation were more robust among antidepressant-treated women than in depressed women, possibly due to the fact that depression is a multifactorial condition with varying degrees of endocrine disruption. It remains to be established whether the alterations found in the gene expression of the placenta are found in the fetus as well.
Unfortunately pregnancy is not a lighthearted period for all women. About 10% of pregnant women in economically developed countries and up to 25% of pregnant women in poorer countries develop symptoms of depression, such as fatigue, trouble sleeping, sense of sadness or hopelessness, during pregnancy (O'Keane and Marsh, 2007). The 5th edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-V) also acknowledged the peripartum onset of depression (American Psychiatric Association, 2013). Antenatal depression is not only affecting the mother's well-being but also affects the unborn child and has been associated with child internalizing difficulties (Barker et al., 2011), attention problems (Van Batenburg-Eddes et al., 2012), and violent behavior during adolescence (Hay et al., 2010). Moreover, the risk of developing depression during adolescence (Pawlby et al., 2009) or adulthood (Pearson et al., 2013) is higher. Although the genetic setup of the mother, the hormonal/reproductive history, current stressors, and life experiences are known risk factors (Miller and LaRusso, 2011), the underlying biological mechanisms of antenatal depression and especially its influence on the developing child remain largely unknown. So far one of the mainly suggested biological mechanism underlying the effects of antenatal depression is the activation of the HPA-axis (reviewed in: Field, 2011; Olivier et al., 2014; Waters et al., 2014). In addition to increased cortisol levels Field et al. (2004) also reported on reduced serotonin and dopamine levels in urine samples of depressed pregnant women. Alterations in cortisol/HPA axis responses or in catecholamines/serotonin may explain some effects in the offspring, however the intrauterine environment is directly passed into the embryo-fetal epigenetic programming. For this reason exposure to antenatal depression in utero may also increase the risk for adverse outcome in the offspring via epigenetic alterations (Babenko et al., 2015). Several treatments for antenatal depression are available, including antidepressant treatment. Antidepressants pass the placenta and are found in the amniotic fluid (Hostetter et al., 2000; Loughhead et al., 2006). Although the exact effect on the offspring is unknown, the use of antidepressants during pregnancy has increased during the last decades. From 1998 to 2005 a 300% increase in antidepressant use during pregnancy was reported (Alwan et al., 2011) and this number is still increasing. Selective serotonin reuptake inhibitors (SSRIs) are the most frequently used antidepressants during pregnancy (Andrade et al., 2008), and are generally considered safe (Gentile, 2005). However, epidemiologic studies have found associations between SSRI use and neurodevelopmental disorders, e.g., autism (Croen et al., 2011; Rai et al., 2013), and attention-deficit hyperactivity disorder (Clements et al., 2014). Cohort studies are underway for the study of SSRI use during pregnancy and the neurodevelopmental disorders in the offspring (Malm et al., 2012). Although some SSRI effects during pregnancy have been reported in the offspring (see: Olivier et al., 2013; Bourke et al., 2014) there is still a great need to investigate the molecular mechanisms involved in antenatal depression that may be altered by antidepressants. By unraveling these pathways we generate more insight into the effects of antenatal depression and antidepressant treatment on the developing child, which ultimately helps in future decisions of using antidepressants during pregnancy.
The placenta plays a pivotal role in supporting fetal growth and development and is a crucial regulator of maternal-fetal interactions and fetal brain development (Hsiao and Patterson, 2012). In fact, placental serotonin synthesis directly modulate fetal brain development (Bonnin et al., 2011). As the placenta carries important information about the pregnancy, investigation of the placenta provides valuable insight to the molecular mechanisms that may have both immediate and long lasting effects on fetal health. This study was designed as a hypothesis-generating study and investigated the impact of antenatal depression and antidepressant treatment during pregnancy on the gene expression in the fetal placenta. Placental samples from mothers with normal pregnancies, from mothers with antenatal depression, and from mothers using antidepressants were collected from the “Biology, Affect, Stress, Imaging and Cognition in Pregnancy and the Puerperium” (BASIC) project. In a first pilot microarray experiment 108 genes were differentially expressed in antenatal depressed mothers whereas 109 genes were differentially expressed in those using antidepressants. Of these genes, seven were chosen for biological validation in a larger sample.
Materials and Methods
This study was carried out at the Department of Women's and Children's health, Uppsala University Hospital, as part of an ongoing longitudinal study on antenatal and postpartum depression: the Biology, Affect, Stress, Imaging and Cognition in Pregnancy and the Puerperium (BASIC) project. The BASIC project started in 2010 and aims to include a study population of 5000 pregnant women in the Uppsala County. Women attending the routine ultrasound examination (gestational week 16–17) at Uppsala University Hospital are approached for participation, enabling a population-based sampling. Upon informed consent, women fill out web-based questionnaires in gestational week 17 and 32 including questions on physical and socio-demographic characteristics, medical, psychiatric, gynecologic and obstetric history variables, lifestyle, medication parameters, and the Swedish version of the Edinburgh Postnatal Depression Scale (EPDS). Information concerning the maternal depression, SSRI use, delivery and neonatal outcome were retrieved from the medical records. Placental biopsies are collected at delivery.
For the specific aim of this sub-study, inclusion criteria were women of Caucasian origin, normal pregnancies and deliveries and healthy offspring (no diagnoses and no admittance to neonatal care). Exclusion criteria were smoking or alcohol use during pregnancy, any daily use of prescribed drugs during pregnancy, any other chronic condition or disease, gestational age <35 weeks, and maternal age <18 or >42 years. Women on antidepressants used their treatment during the entire pregnancy in clinically relevant doses (low-dose use was excluded). The study was approved by the Regional Ethics Committee, Uppsala, Sweden, and performed in accordance with relevant guidelines and regulations.
Study population for micro-array analysis
Women with pregnancies complicated by ongoing depression (n = 5), SSRI treatment (n = 5) and women with normal pregnancies (n = 10) were selected from the BASIC biobank. Depressed women had medical records confirming major depression and ongoing treatment for their depression in terms of psychotherapy. In the SSRI group sertraline (n = 3), fluoxetine (n = 1) and escitalopram (n = 1) was used. Women on SSRIs displayed significantly lower depression scores than depressed women, i.e., the two groups were not readily comparable (Table 1). Hence, the exposure of ongoing depression or SSRI treatment, respectively, were compared against two control groups. Controls were matched with respective depressed or SSRI-treated women by age (±2 years), BMI (±one unit) and gestational length (±1 week) on an individual level. The control group consisted of women with no history and no current symptoms/diagnoses of mood or anxiety disorders, and their EPDS scores at gestational week 17 and 32 were 6 or lower.
Study population for validation of the microarray
The samples described above were extended to 24 women with pregnancies complicated by ongoing depression, 29 antidepressant-treated women and 31 women with normal pregnancies. The depressed group included the five microarray cases and 18 women with EPDS >12 in gestational week 17 and 32, or an EPDS score >14 on at least one time point (n = 1). The average EPDS score for all depressed women was >15 in gestational week 17 and 32. The EPDS questionnaire is validated for use in both pregnant and postpartum women (Gibson et al., 2009), and has been validated for the Swedish setting (Rubertsson et al., 2011). The EPDS contains ten items (rated on a scale from 0 to 3), based on the past 7 days. While its sensitivity is relatively low, a cut-off score of >12 points during pregnancy has a specificity of 98–99% for major depression (Bergvink et al., 2011). Thirteen of the depressed women were also evaluated by Mini International Neuropsychiatric Interview, which confirmed that all but one (she had social phobia) had major depressive disorder during pregnancy. In the antidepressant-treated group, women used sertraline (n = 11), fluoxetine (n = 8), citalopram/escitalopram (n = 7), venlafaxine (n = 2) and clomipramine (n = 1). Treatment had been initiated by primary care physicians as well as by psychiatrists.
Sample Collection, Processing, and Storing
Placental tissues (containing both maternal and fetal side) were obtained after delivery, rinsed carefully in sterile phosphate-buffered saline to wash off maternal and fetal blood, and frozen on dry ice within 60 min of delivery and stored at −70°C until further use. Each placenta was individually processed as a single biological replicate in the microarray and validation study.
A biopsy was taken with a 3 mm cube from the fetal side of the placenta. Total RNA was isolated using miRNeasy mini kit (Qiagen, Hilden, Germany). Tissue was lysed with QIAzol reagent (Qiagen) using a rotor-stator homogenizer (up to 33.000 rpm; Ingenieursbűro CAT M Zipper Gmbh, type x120, Staufen, Germany) and chloroform (Sigma Aldrich, St. Louis, MO, USA) was added for phase-separation. The rest of the procedure was performed as described in manufactures protocol.
A biopsy was taken from the fetal side of the placenta with a 3 mm cube. Total RNA was isolated using RNeasy mini (Qiagen, Hilden, Germany). Tissue was lysed with QIAzol reagent (Qiagen) using TissueLyser (20Hz, 2 × 5 min) with stainless steel beads (Qiagen) and chloroform was added for phase-separation. The rest of the procedure was performed as described in manufactures protocol.
For both studies RNA concentration was measured with ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, Delaware, USA) and RNA quality was evaluated using the Agilent 2100 Bioanalyzer system (Agilent Technologies Inc, Palo Alto, California, USA).
Microarray Expression Analysis
250 nanograms of total RNA from each sample were used to generate amplified and biotinylated sense-strand cDNA from the entire expressed genome according to the Ambion WT Expression Kit (P/N 4425209 Rev B 05/2009) and Affymetrix GeneChip® WT Terminal Labeling and Hybridization User Manual (P/N 702808 Rev. 4, Affymetrix Inc., Santa Clara, CA). GeneChip® ST Arrays (GeneChip® XXX Gene 1.0 ST Array) were hybridized for 16 h in a 45°C incubator, rotated at 60 rpm. According to the GeneChip® Expression Wash, Stain and Scan Manual (PN 702731 Rev 3, Affymetrix Inc., Santa Clara, CA) the arrays were then washed and stained using the Fluidics Station 450 and finally scanned using the GeneChip® Scanner 3000 7G.
Microarray Data Analysis
The raw data was normalized in the free software Expression Console provided by Affymetrix (http://www.affymetrix.com) using the robust multi-array average (RMA) method first suggested by Li and Wong (2001) and Irizarry et al. (2003). Subsequent analysis of gene expression data was carried out in the freely available statistical computing language R (http://www.r-project.org) using packages available from the Bioconductor project (www.bioconductor.org). In order to search for differentially expressed genes between the depressed and controls, and the SSRI and the control groups an empirical Bayes moderated t-test with robust regression was applied (Smyth, 2004), using the “limma” package (Smyth, 2005). To address the problem of multiple testing, p-values were adjusted according to Benjamini and Hochberg (1995). The Genesis software, version 1.7.1 (http://genome.tugraz.at/), was used to produce hierarchical clustering and to visualize differentially expressed genes by heat maps (Sturn et al., 2002). The expression data were further analyzed using ingenuity pathway analysis (IPA) in order to determine significantly deregulated genes and pathways (Ingenuity® Systems, Mountain View, CA, USA; www.ingenuity.com). IPA reveals Top genes, which are genes with the largest normalized enrichment scores, and IPA computes a score for each network according to the fit of that network to the user-defined set of Focus Genes. The score, derived from a p-value, indicates the likelihood of the association between the set of focus genes (Bonferroni-corrected significance threshold of P < 0.05 and a fold change of 0.5) and a given pathway. The smaller the p-value, the less likely that the association is found due to random chance. In general P-values below 0.05 indicate a non-random significant association. The p-value is calculated using the right-tailed Fisher Exact Test (for more details see www.ingenuity.com).
cDNA was synthesized using SuperScriptIII reversed transcriptase (Invitrogen, Carlsbad, California, USA) according to manufacturer's protocol. Briefly, 250 ng of total RNA was used to reverse transcribe using the random primer to prepare 20 μl of cDNA.
Real-Time Quantitative Reverse Transcriptase Polymerase Chain Reaction Analysis
The validity of the microarray results was tested via quantitative real-time PCR (qRT-PCR) employing the StepOne Plus qPCR machine (Applied Biosystems, Life Technologies, Carlsbad, California, USA). For validation we selected seven genes [NEXN (Hs00332124_m1), GCC2 (Hs00206083_m1), ROCK 1 (Hs01127699_m1), ROCK2 (Hs00178154_m1), DNM1 (Hs00247147_m1), KTN1 (Hs00192160_m1), and C12ORF39 (Hs00228976_m1)] which showed a fold change (FC) > 0.5 in the microarray. GAPDH (Hs99999905_m1) and β-actin (4326315E) were selected as reference genes for normalization. cDNA of the samples was used for quantification. TaqMan Gene expression Assay primers, probes and gene expression master mix (all Applied Biosystems, Life Technologies, Carlsbad, CA, USA) were used to run the qRT-PCR according to manufacturer's instructions. Mean plate efficiencies were calculated by LinReg.
qPCR Data Analysis
All samples were performed in triplicates and averaged for further calculations. Mean normalized expression (MNE) based on the ratio between Ct-values of target and reference genes and the efficiency of the PCR reactions, was calculated as a measure of target gene transcription, as described previously (Muller et al., 2002; Helmestam et al., 2012). Data are presented as log2 MNE to illustrate the difference between the groups.
Clinical characteristics of women in the microarray study and the validation study were compared by means of One-Way ANOVA, followed by a Bonferroni post hoc test when appropriate. Differences in expression between subsets in the validation study (qPCR) were calculated using a univariate ANOVA with age, BMI, parity, and week of delivery as covariates. Data were analyzed using the SPSS 20.0 software. Level of significance was set at P < 0.05. Data are presented as mean ± S.E.M.
The demographic characteristics of the women in the microarray and validation studies are presented in Tables 1, 2. Microarray study: A significant group difference was found for EPDS scores at gestational week 17 [F(3, 18) = 8.9, P < 0.01] and gestational week 32 [F(3, 17) = 22.8, P < 0.001]. Depressed women had significantly higher EPDS scores in gestational week 17 and 32 compared with controls (P < 0.01), whereas women on SSRIs did not differ from controls at any time-point. In gestational week 32, depressed women also had significantly higher EPDS scores than SSRI-treated women, Table 1. Validation study: All women on SSRI treatment reported previous anxiety and/or depression at the first antenatal booking, whereas 41% of depressed cases had no previous psychiatric history. Women using antidepressants had a significantly shorter gestational length than controls (P < 0.01). As expected, significant group differences were found in EPDS scores at gestational weeks 17 [F(2, 74) = 69.6, P < 0.001] and 32 [F(2, 74) = 74.3, P < 0.001]. In gestational week 17 and 32, depressed women had significantly higher EPDS scores than antidepressant-treated women (P < 0.001 and P < 0.001, respectively), whom in turn had higher EPDS scores than controls (P < 0.001 and P < 0.001, respectively), Table 2. No other parameters differed between groups.
Differentially Expressed Genes and Pathways Between the Depressed and Control Placentas
At a Bonferroni-corrected significance threshold (P < 0.05) and a log2-fold change of 0.5 or higher we found 108 genes differentially expressed between the depressed women and their respective controls; 100 were down-regulated and 8 up-regulated, see Table 3. The raw microarray data is found as an excel file in the Supplementary Data. The ingenuity pathway analysis (IPA) revealed 17 differentially expressed top genes; 7 were up-regulated and 10 down-regulated. Top up- and down-regulated molecules are summarized in Table 4. We then clustered placentas according to their gene expression profiles for the 17 genes that displayed differential expression (see Figure 1A). In order to determine the biological relevance, analysis with IPA was performed, focusing on genes that differed in expression between placentas from depressed women and controls. As shown in Table 5, we identified five gene networks that were significantly enriched, classified as follows: (I) DNA Replication, Recombination, and Repair, Cellular Assembly and Organization, Cell Cycle with an IPA score of 27; (II) Cell Cycle, Cancer, Connective Tissue Disorders with a IPA score of 18; (III) Cellular movement, Hematological System Development and Function, Immune Cell Trafficking with an IPA score of 16; (IV) Cardiovascular System Development and Function, Organismal Development, Visual System Development and Function with an IPA score of 14; and (V) Molecular Transport, RNA Trafficking, Connective Tissue Disorders with an IPA score of 9 (Table 5). Further, the significant canonical pathways identified by IPA (P < 0.05) are shown in Table 6, along with the included genes and p-values. Pathways included Actin Nucleation by ARP-WASP Complex, RhoA Signaling, VEGF Signaling, Protein Kinase A Signaling, and 1D-myo-inositol Hexakisphosphate Biosynthesis V [from Ins(1,3,4)P3].
Table 3. Significantly up and down regulated genes in the control vs. depressed groups (− = reduction in expression levels).
Table 4. Significantly up- and down-regulated top molecules in the control vs. the depressed group (− = reduction in expression levels).
Figure 1. Visualization of differentially expressed genes using hierarchical clustering of genes in depressed (A) and SSRI-treated (B) vs. control fetal placentas.
Table 5. Enriched ingenuity pathway analysis (IPA) categories including differentially expressed genes in the depressed group.
Differentially Expressed Genes and Pathways Between SSRI-Treated and Control Placentas
Similarly, we found 109 genes to be differentially expressed between the SSRI-treated women and their respective controls at a Bonferroni-corrected significance (P < 0.05) threshold with a fold change of 0.5 or higher. 82 genes were down-regulated and 27 up-regulated, see Table 7. The raw microarray data is found as an excel file in the Supplementary Data. IPA analysis revealed 20 differentially expressed top genes, of which 10 were up- and 10 down-regulated (see Table 8). We then clustered placentas according to their gene expression profiles for the 20 genes that displayed differential expression (see Figure 1B). With use of IPA we focused on genes that differed in expression between placentas from antidepressant-treated women and controls. As shown in Table 9, we identified five gene networks that were significantly enriched. Of biological relevance were: (I) Infectious Disease, Cellular Assembly and Organization, Cellular Function and Maintenance with an IPA score of 13; (II) Cellular Growth and Proliferation, Inflammatory Response, Lipid Metabolism with an IPA score of 11; (III) Cell Death and Survival, Inflammatory Response, Cellular Movement with an IPA score of 9; (IV) Cell death and Survival, Liver Necrosis/Cell Death, Hematological System Development and Function with an IPA score of 8; and (V) Cardiovascular Disease, Skeletal and Muscular Disorders, Cardiovascular System Development and Function with an IPA score of 2. Further, the significant canonical pathways identified by IPA (P < 0.05) are shown in Table 10, along with the included genes and p-values. Pathways included Ephrin A Signaling, RhoA Signaling, PEDF Signaling, Breast Cancer Regulation by Stathmin1, and Signaling by Rho Family GTPases.
Table 7. Significantly up and down regulated genes in the control vs. SSRI groups (− = reduction in expression levels).
Table 8. Significantly up- and down-regulated top molecules in the control vs. SSRI group (- = reduction in expression levels).
Table 9. Enriched ingenuity pathway analysis (IPA) categories including differentially expressed genes in the SSRI group.
Commonly Altered Genes in Depressed and SSRI-Treated Groups
Of the 108 genes that were differentially expressed between the depressed and the control cases, and the 109 genes that were differentially expressed between the SSRI-treated and the control cases, only 20 genes were overlapping. These genes are displayed in Table 11.
Table 11. Genes commonly altered in depressed and SSRI groups compared with controls. (− = reduction in expression levels).
Validation of Microarray Data Using qPCR Analysis
For validation of the microarray results we selected seven genes that were detected in top up- or down-regulated genes, in the pathway analysis or in the canonical pathway analysis in placental tissue of both depressed and antidepressant-treated women. ROCK1 and ROCK2 are both involved in the actin nucleation by ARP-WASP complex, RhoA signaling, VEGF Signaling and protein kinase A signaling in the canonical pathway analysis of depressed women. Moreover they also appeared in the canonical pathway analysis of ephrin A Signaling, RhoA signaling, PEDF signaling, breast cancer regulation by stathmin 1 and signaling by Rho family GTPases of antidepressant-treated women. In addition, ROCK1 and GCC2 belonged to the top down-regulated genes in depressed as well as antidepressant-treated women and ROCK2 is a top down-regulated gene in the depressed women as well. KTN1 is involved in molecular transport, RNA trafficking, connective tissue disorders pathway of depressed women, but also in infectious disease, cellular assembly and organization, cellular function and maintenance pathway of antidepressant-treated women. Moreover, KTN1 is involved in the RhoA signaling of the canonical pathway analysis in both depressed and antidepressant women. DNM1L is involved in cardiovascular system development and function, organismal development, visual system development and function of the pathway analysis in depressed women but also in cell death and survival, inflammatory response, cellular movement pathway analysis of antidepressant-treated women. NEXN was the top down-regulated gene in antidepressant-treated women and also turned out to be involved in cardiovascular disease, skeletal and muscular disorders, cardiovascular system development and function pathway in antidepressant-treated women. Finally, C12orf39 was chosen as it appeared to be the top up-regulated gene in antidepressant-treated women.
In the placenta of depressed women (Figure 2A) a significant down-regulation compared to placenta of controls was found for the C12orf39 gene [F(1, 45) = 3.83, p = 0.05] and a tendency for down-regulation of the ROCK2 gene was found [F(1, 44) = 3.03, p = 0.08]. No other genes were differentially expressed between control and depressed placentas (NEXN [F(1, 45) = 0.57], GCC2 [F(1, 45) = 0.52], ROCK1 [F(1, 45) = 2.43], DNM1L [F(1, 45) = 2.53] and KTN1 [F(1, 44) = 1.57]). When the placental gene expression was compared between antidepressant-treated women and controls (Figure 2B), we found a significant down-regulation of ROCK1 [F(1, 47) = 4.26, P < 0.05], ROCK2 [F(1, 46) = 9.48, P < 0.01], GCC2 [F(1, 47) = 3.78, p = 0.05], KTN1 [F(1, 46) = 6.31, P < 0.05], DNM1L [F(1, 47) = 6.40, P < 0.05], and a tendency for down-regulation of C12orf39 [F(1, 47) = 3.39, P = 0.07]. The gene expression of NEXN [F(1, 47) = 1.28, ns] did not differ between antidepressant-treated women and controls.
Figure 2. Differences in gene expression (qPCR) between subsets in the validation study of the microarray data. Log2 mean normalized expression (MNE) is shown for the ROCK1, ROCK2, GCC2, KTN1, DNM1L, NEXN, and C12orf39 genes in depressed (A) and SSRI women (B). (A) *P = 0.05, #P = 0.08, (B) *P ≤ 0.05, **P < 0.01, #P = 0.07.
We performed a gene expression study in the fetal placenta of depressed women and antidepressant-treated women, and compared them with the gene expression of the placentas from women with normal pregnancies. We found that antenatal depression and antidepressant exposure during pregnancy has an influence on the gene expression of the placenta. In the microarray 108 genes were differentially expressed in women with antenatal depression, while 109 genes were differentially expressed in antidepressant-treated women. Only 20 genes were overlapping between depressed women and women on antidepressant treatment. Among the genes we chose for validation, only 2 were validated with qPCR for depressed women. In the antidepressant-treated women, 6 genes were validated, indicating a more robust effect in alterations of these genes due to antidepressant treatment during pregnancy.
Antenatal depression is a relatively heterogeneous condition with different causes (primary or secondary to somatic disease) and differential degree of endocrine disturbances. Furthermore, women with antenatal depression may also differ between the pilot microarray and the validation study as to depression severity and duration of the depressive episode. These factors may have precluded the possibility to confirm the microarray findings, and it is a major limitation that not all women in the validation part of the study were diagnosed by a structured psychiatric interview. As depression per se has effects on perinatal outcomes such as birth weight and gestational age (Chambers et al., 1996), suggesting that placental function is altered, further studies in more homogeneous women (and with larger sample sizes) of depression are warranted. In addition, women on antidepressant treatment during pregnancy possibly may have had a more severe depression, since continued treatment apparently was needed, and that this effect is reflected in the gene expression pattern. However, it is also possible that alterations in placental gene expression already occur because of the antenatal depression, but become (more) apparent, when antidepressants are used. Although several alterations in the gene expression of the fetal placenta were found, it remains to establish if these alterations are found in the fetus as well.
When the microarray was validated with a larger sample-size, ROCK2 was down-regulated in both depressed and SSRI-treated women. ROCK1 and ROCK2 are part of the Rho-associated coiled-coil kinase family (Nakagawa et al., 1996; Amano et al., 2000) and are downstream effectors of RhoA-GTP. Rho-ROCK signaling pathways are involved in the regulation of actin cytoskeleton, cell migration and proliferation (Schofield and Bernard, 2013). In mice, ROCK1 is highly expressed in the lung, liver, spleen, kidney and testis, whereas ROCK2 is most abundant in the brain and heart (Nakagawa et al., 1996; Di et al., 2000; Wei et al., 2001). The role of the Rho/ROCK family in cardiovascular diseases has been extensively studied (Shi and Wei, 2013). Cardiac malformations (Diav-Citrin et al., 2008), including pulmonary hypertension (Chambers et al., 1996, 2006; Källén and Olausson, 2008; Kieler et al., 2012), have also been reported in the SSRI-exposed offspring. In the present study we found that ROCK2 was down-regulated in antidepressant-exposed placentas (and to lesser extent in placentas of depressed women). Moreover, IPA analysis revealed that SSRI treatment affects the “Cardiovascular Disease, Skeletal and Muscular Disorders, Cardiovascular System Development and Function” network. Furthermore, it seems that the use of SSRIs intensifies the alterations in ROCK2 expression compared to depression. Due to lowered ROCK2 expression found in the fetal placenta of antidepressant-treated women, it is tempting to speculate that a normal expression of ROCK2 in placenta is important for a normal function of the cardiovascular system in the fetus.
Interestingly, we also found that NEXN was down-regulated in antidepressant-treated women. NEXN is a Z-disk gene which is associated with dilated cardiomyopathy (Hassel et al., 2009). Hence, further studies are necessary to investigate the effects of down-regulated placental NEXN and ROCK2 on the development of the cardiovascular system in the fetus.
Besides the role of the Rho kinase pathway in cardiovascular diseases, a role has also been proposed for the modulation of the placental vasculature. Although expression of ROCK1 and ROCK2 was not different, a higher RhoA mRNA expression was found in placentae from women who suffered from preeclampsia compared with placentae from those that were normotensive (Friel et al., 2008). Interestingly, the use of antidepressants have been associated with an increased risk for preeclampsia (Palmsten et al., 2012). In addition, normal ROCK1 and ROCK2 activity is required for normal inner cell mass morphogenesis, which is of importance for successful fetal development (Laeno et al., 2013). These data indicate that antidepressant use, mainly SSRI, during pregnancy may influence pregnancy complications and fetal development, and that ROCK1 and ROCK 2 may be involved with these processes.
Another gene that was down-regulated in depressed women and to a lesser extent in antidepressant-treated women is the C12orf39. C12orf39 is mainly expressed in the placenta and brain, suggesting that C12orf39 may function in these active secretory tissues (Wan et al., 2010). With regard to its function, C12orf39 is mainly extracellular and located in the villous trophoblasts. Trophoblasts are important in exchanges between the fetus and the mother and possess endocrine activity, releasing hormones that are important in the homoeostasis of pregnancy (Lunghi et al., 2007). In addition, trophoblasts are involved with the secretion of placental growth hormone and are related to the development of the placenta (Zeck et al., 2008). Together, these findings suggest that C12orf39 is implicated in the regulation of placenta development by its role in the biological functions of the trophoblasts and that antenatal depression during pregnancy may influence this development. Of interest is the fact that in antidepressant-treated women the effect is no longer significant, which may indicate that antidepressants may restore the effects of the depression-induced down-regulation of C12of39.
GCC2 is a peripheral membrane protein localized to the trans-Golgi network (Luke et al., 2003) which is involved in the maintenance of Golgi structure or transport vesicle tethering (Brown et al., 2011). More research is needed to investigate the effects of altered gene expression of GCC2 in the placenta on the developing fetus.
Similarly, KTN1 encodes the full kinectin which is found in the endoplasmic reticulum and is responsible for the transport of vesicles along microtubules. KTN1 is mainly expressed in the brain, liver, ovarian, and hematopoietic cells (Tran et al., 2002; Bai et al., 2006). Of interest is that kinectin can interact with RhoA (Hotta et al., 1996), and that the RhoA signaling pathway is affected in the depressed and antidepressant-treated women. As described before normal RhoA signaling is important to prevent pregnancy complications. Possible consequences for the fetus due to the down-regulation of KTN1 in the placenta remains to be investigated.
The last gene that was validated was DNM1L, which is a GTPase regulating the mitochondrial fission. In mice it was shown that ablation of the DNM1L gene induced defects in trophoblast giant cells and cardiomyocytes (Wakabayashi et al., 2009). Moreover, brain-specific DNM1L ablation caused developmental defects in the cerebellum (Wakabayashi et al., 2009). Although these results were found in knockout mice it is tempting to speculate that the down-regulation of the DNM1L found in the antidepressant-treated women might have an effect on embryonic and brain development as well. However, again, this effect remains to be established.
Despite the strengths of our study, such as the longitudinal nature of the study and the information on the state of the mothers mood at multiple time points, some limitations need to be discussed. First, we investigated alterations in gene expression of the fetal side of the placenta due to antenatal depression and antidepressant treatment. The results of this study may give us an indication on altered pathways in the placenta. However, the placenta is a separate organ and is not part of the fetus itself, therefore findings need to be replicated in the developing fetus. In humans this is not easily feasible therefore experiments are ongoing in a rodent study. Second, antenatal depression is a relatively heterogeneous condition and outcome of diagnoses and treatment plans (before the women entered the study) were diagnosed by different doctors which may have biased the outcome. As a result dosages and types of medication may not have been appropriate for the diagnosed depression. Nevertheless, all women did undergo the EPDS screening providing comparable data between the groups concerning the mood state at different time points during (and after) pregnancy. Third, in the validation study we included different types of antidepressants, although they were mainly SSRIs, this may have influenced the outcome of the gene expression in the SSRI treatment group.
In conclusion, gene expression in the fetal placenta is altered by antenatal depression and SSRI treatment. As more placental genes alterations were validated in a larger subset of SSRI-treated women compared to those with antenatal depression we conclude that for these subset of genes, the effects of SSRI-intake during pregnancy are more robust. It remains to be established how these differentially affected genes influence the development of the child, and whether these differences are found in the fetus as well.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by grants from the Swedish Research Council (Grant No. K2014-54X-20642-07-4), the Marianne and Marcus Wallenberg Foundation (2010:0031), KI fonder research (2013fobi37758), and the Swedish Society of Medicine (SLS-303881 & SLS-384001).
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fncel.2014.00465/abstract
Alwan, S., Reefhuis, J., Rasmussen, S. A., and Friedman, J. M. (2011). Patterns of antidepressant medication use among pregnant women in a United States population. J. Clin. Pharmacol. 51, 264–270. doi: 10.1177/0091270010373928
Andrade, S. E., Raebel, M. A., Brown, J., Lane, K., Livingston, J., Boudreau, D., et al. (2008). Use of antidepressant medications during pregnancy: a multisite study. Am. J. Obstet. Gynecol. 198, 194–195. doi: 10.1016/j.ajog.2007.07.036
Babenko, O., Kovalchuk, I., and Metz, G. A. (2015). Stress-induced perinatal and transgenerational epigenetic programming of brain development and mental health. Neurosci. Biobehav. Rev. 48C, 70–91. doi: 10.1016/j.neubiorev.2014.11.013
Bai, J. Z., Mon, Y., and Krissansen, G. W. (2006). Kinectin participates in microtubule-dependent hormone secretion in pancreatic islet beta-cells. Cell Biol. Int. 30, 885–894. doi: 10.1016/j.cellbi.2006.06.008
Barker, E. D., Jaffee, S. R., Uher, R., and Maughan, B. (2011). The contribution of prenatal and postnatal maternal anxiety and depression to child maladjustment. Depress. Anxiety 28, 696–702. doi: 10.1002/da.20856
Bergvink, V., Kooistra, L., Lambregtse-van den Berg, M. P., Wijnen, H., Bunevicius, R., van Baar, A., et al. (2011). Validation of the Edinburgh depression scale during pregnancy. J. Psychosom. Res. 70, 385–389. doi: 10.1016/j.jpsychores.2010.07.008
Brown, F. C., Schindelhaim, C. H., and Pfeffer, S. R. (2011). GCC185 plays independent roles in Golgi structure maintenance and AP-1-mediated vesicle tethering. J. Cell Biol. 194, 779–787. doi: 10.1083/jcb.201104019
Chambers, C. D., Hernandez-Diaz, S., Van Marter, L. J., Werler, M. M., Louik, C., Jones, K. L., et al. (2006). Selective serotonin-reuptake inhibitors and risk of persistent pulmonary hypertension of the newborn. N. Engl. J. Med. 354, 579–587. doi: 10.1056/NEJMoa052744
Chambers, C. D., Johnson, K. A., Dick, L. M., Felix, R. J., and Jones, K. L. (1996). Birth outcomes in pregnant women taking fluoxetine. N. Engl. J. Med. 335, 1010–1015. doi: 10.1056/NEJM199610033351402
Clements, C. C., Castro, V. M., Blumenthal, S. R., Rosenfield, H. R., Murphy, S. N., Fava, M., et al. (2014). Prenatal antidepressant exposure is associated with risk for attention-deficit hyperactivity disorder but not autism spectrum disorder in a large health system. Mol. Psychiatry. doi: 10.1038/mp.2014.90 [Epub ahead of print].
Croen, L. A., Grether, J. K., Yoshida, C. K., Odouli, R., and Hendrick, V. (2011). Antidepressant use during pregnancy and childhood autism spectrum disorders. Arch. Gen. Psychiatry 68, 1104–1112. doi: 10.1001/archgenpsychiatry.2011.73
Di, C. F., Imarisio, S., Hirsch, E., Broccoli, V., Bulfone, A., Migheli, A., et al. (2000). Defective neurogenesis in citron kinase knockout mice by altered cytokinesis and massive apoptosis. Neuron 28, 115–127. doi: 10.1016/S0896-6273(00)00090-8
Diav-Citrin, O., Shechtman, S., Weinbaum, D., Arnon, J., di Gianantonio, E., Clementi, M., et al. (2008). Paroxetine and fluoxetine in pregnancy: a multicenter, prospective, controlled study. Reprod. Toxicol. 20, 459. doi: 10.1111/j.1365-2125.2008.03261.x
Field, T., Diego, M., Dieter, J., Hernandez-Reif, M., Schanberg, S., Kuhn, C., et al. (2004). Prenatal depression effects on the fetus and the newborn. Infant Behav. Dev. 27, 216–229. doi: 10.1016/j.infbeh.2003.09.010
Friel, A. M., Hynes, P. G., Sexton, D. J., Smith, T. J., and Morrison, J. J. (2008). Expression levels of mRNA for Rho A/Rho kinase and its role in isoprostane-induced vasoconstriction of human placental and maternal vessels. Reprod. Sci. 15, 179–188. doi: 10.1177/1933719107310306
Gibson, J., McKenzie-McHarg, K., Shakespeare, J., Price, J., and Gray, R. (2009). A systematic review of studies validating the Edinburgh Postnatal Depression Scale in antepartum and postpartum women. Acta Psychiatr. Scand. 119, 350–364. doi: 10.1111/j.1600-0447.2009.01363.x
Hassel, D., Dahme, T., Erdmann, J., Meder, B., Huge, A., Stoll, M., et al. (2009). Nexilin mutations destabilize cardiac Z-disks and lead to dilated cardiomyopathy. Nat. Med. 15, 1281–1288. doi: 10.1038/nm.2037
Hay, D. F., Pawlby, S., Waters, C. S., Perra, O., and Sharp, D. (2010). Mothers' antenatal depression and their children's antisocial outcomes. Child Dev. 81, 149–165. doi: 10.1111/j.1467-8624.2009.01386.x
Helmestam, M., Andersson, H., Stavreus-Evers, A., Brittebo, E., and Olovsson, M. (2012). Tamoxifen modulates cell migration and expression of angiogenesis-related genes in human endometrial endothelial cells. Am. J. Pathol. 180, 2527–2535. doi: 10.1016/j.ajpath.2012.02.026
Hostetter, A., Stowe, Z. N., Strader, J. R. Jr., McLaughlin, E., and Llewellyn, A. (2000). Dose of selective serotonin uptake inhibitors across pregnancy: clinical implications. Depress. Anxiety 11, 51–57. doi: 10.1002/(SICI)1520-6394(2000)11:2<51::AID-DA1>3.0.CO;2-R
Hotta, K., Tanaka, K., Mino, A., Kohno, H., and Takai, Y. (1996). Interaction of the Rho family small G proteins with kinectin, an anchoring protein of kinesin motor. Biochem. Biophys. Res. Commun. 225, 69–74. doi: 10.1006/bbrc.1996.1132
Irizarry, R. A., Hobbs, B., Collin, F., Beazer-Barclay, Y. D., Antonellis, K. J., Scherf, U., et al. (2003). Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4, 249–264. doi: 10.1093/biostatistics/4.2.249
Källén, B., and Olausson, P. O. (2008). Maternal use of selective serotonin re-uptake inhibitors and persistent pulmonary hypertension of the newborn. Pharmacoepidemiol. Drug Saf. 17, 801–806. doi: 10.1002/pds.1570
Kieler, H., Artama, M., Engeland, A., Ericsson, O., Furu, K., Gissler, M., et al. (2012). Selective serotonin reuptake inhibitors during pregnancy and risk of persistent pulmonary hypertension in the newborn: population based cohort study from the five Nordic countries. BMJ 344:d8012. doi: 10.1136/bmj.d8012
Laeno, A. M., Tamashiro, D. A., and Alarcon, V. B. (2013). Rho-associated kinase activity is required for proper morphogenesis of the inner cell mass in the mouse blastocyst. Biol. Reprod. 89, 122. doi: 10.1095/biolreprod.113.109470
Loughhead, A. M., Fisher, A. D., Newport, D. J., Ritchie, J. C., Owens, M. J., DeVane, C. L., et al. (2006). Antidepressants in amniotic fluid: another route of fetal exposure. Am. J. Psychiatry 163, 145–147. doi: 10.1176/appi.ajp.163.1.145
Luke, M. R., Kjer-Nielsen, L., Brown, D. L., Stow, J. L., and Gleeson, P. A. (2003). GRIP domain-mediated targeting of two new coiled-coil proteins, GCC88 and GCC185, to subcompartments of the trans-Golgi network. J. Biol. Chem. 278, 4216–4226. doi: 10.1074/jbc.M210387200
Malm, H., Artama, M., Brown, A. S., Gissler, M., Gyllenberg, D., Hinkka-Yli-Salomaki, S., et al. (2012). Infant and childhood neurodevelopmental outcomes following prenatal exposure to selective serotonin reuptake inhibitors: overview and design of a Finnish Register-Based Study (FinESSI). BMC Psychiatry 12:217. doi: 10.1186/1471-244X-12-217
Nakagawa, O., Fujisawa, K., Ishizaki, T., Saito, Y., Nakao, K., and Narumiya, S. (1996). ROCK-I and ROCK-II, two isoforms of Rho-associated coiled-coil forming protein serine/threonine kinase in mice. FEBS Lett. 392, 189–193. doi: 10.1016/0014-5793(96)00811-3
Olivier, J. D., Åkerud, H., Kaihola, H., Pawluski, J. L., Skalkidou, A., Hogberg, U., et al. (2013). The effects of maternal depression and maternal selective serotonin reuptake inhibitor exposure on offspring. Front. Cell Neurosci. 7:73. doi: 10.3389/fncel.2013.00073
Olivier, J. D., Akerud, H., and Sundström-Poromaa, I. (2014). Antenatal depression and antidepressants during pregnancy: unraveling the complex interactions for the offspring. doi: 10.1016/j.ejphar.2014.07.049 [Epub ahead of print].
Palmsten, K., Setoguchi, S., Margulis, A. V., Patrick, A. R., and Hernandez-Diaz, S. (2012). Elevated risk of preeclampsia in pregnant women with depression: depression or antidepressants? Am. J. Epidemiol. 175, 988–997. doi: 10.1093/aje/kwr394
Pawlby, S., Hay, D. F., Sharp, D., Waters, C. S., and O'Keane, V. (2009). Antenatal depression predicts depression in adolescent offspring: prospective longitudinal community-based study. J. Affect. Disord. 113, 236–243. doi: 10.1016/j.jad.2008.05.018
Pearson, R. M., Evans, J., Kounali, D., Lewis, G., Heron, J., Ramchandani, P. G., et al. (2013). Maternal depression during pregnancy and the postnatal period: risks and possible mechanisms for offspring depression at age 18 years. JAMA Psychiatry 70, 1312–1319. doi: 10.1001/jamapsychiatry.2013.2163
Rai, D., Lee, B. K., Dalman, C., Golding, J., Lewis, G., and Magnusson, C. (2013). Parental depression, maternal antidepressant use during pregnancy, and risk of autism spectrum disorders: population based case-control study. BMJ 346:f2059. doi: 10.1136/bmj.f2059
Rubertsson, C., Borjesson, K., Berglund, A., Josefsson, A., and Sydsjo, G. (2011). The Swedish validation of Edinburgh Postnatal Depression Scale (EPDS) during pregnancy. Nord. J. Psychiatry 65, 414–418. doi: 10.3109/08039488.2011.590606
Smyth, G. K. (2005). “Limma:linear models for microarray data,” in Bioinformatics and Computational Biology Solutions using R and Bioconductor, eds R. Gentlemen, V. Carey, S. Dudoit, R. Irizarry, and W. Huber (NewYork, NY: Springer), 397–420.
Van Batenburg-Eddes, T., Brion, M. J., Henrichs, J., Jaddoe, V. W., Hofman, A., Verhulst, F. C., et al. (2012). Parental depressive and anxiety symptoms during pregnancy and attention problems in children: a cross-cohort consistency study. J. Child Psychol. Psychiatry 54, 591–600. doi: 10.1111/jcpp.12023
Wakabayashi, J., Zhang, Z., Wakabayashi, N., Tamura, Y., Fukaya, M., Kensler, T. W., et al. (2009). The dynamin-related GTPase Drp1 is required for embryonic and brain development in mice. J. Cell Biol. 186, 805–816. doi: 10.1083/jcb.200903065
Wan, B., Wang, X. R., Zhou, Y. B., Zhang, X., Huo, K., and Han, Z. G. (2010). C12ORF39, a novel secreted protein with a typical amidation processing signal. Biosci. Rep. 30, 1–10. doi: 10.1042/BSR20080156
Waters, C. S., Hay, D. F., Simmonds, J. R., and van Goozen, S. H. (2014). Antenatal depression and children's developmental outcomes: potential mechanisms and treatment options. Eur. Child Adolesc. Psychiatry 23, 957–971. doi: 10.1007/s00787-014-0582-3
Zeck, W., Widberg, C., Maylin, E., Desoye, G., Lang, U., McIntyre, D., et al. (2008). Regulation of placental growth hormone secretion in a human trophoblast model–the effects of hormones and adipokines. Pediatr. Res. 63, 353–357. doi: 10.1203/01.pdr.0000304935.19183.07
Keywords: antenatal depression, antidepressants, fetal, placenta, gene expression, microarray
Citation: Olivier JDA, Åkerud H, Skalkidou A, Kaihola H and Sundström-Poromaa I (2015) The effects of antenatal depression and antidepressant treatment on placental gene expression. Front. Cell. Neurosci. 8:465. doi: 10.3389/fncel.2014.00465
Received: 30 October 2014; Accepted: 20 December 2014;
Published online: 13 January 2015.
Edited by:Shawn Hayley, Carleton University, Canada
Copyright © 2015 Olivier, Åkerud, Skalkidou, Kaihola and Sundström-Poromaa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jocelien D. A. Olivier, Department of Behavioral Physiology, Center for Behaviour and Neurosciences, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, Netherlands e-mail: firstname.lastname@example.org