Risk of spontaneous preterm birth and fetal growth associates with fetal SLIT2

Abstract

Spontaneous preterm birth (SPTB) is the leading cause of neonatal death and morbidity worldwide. Both maternal and fetal genetic factors likely contribute to SPTB. We performed a genome-wide association study (GWAS) on a population of Finnish origin that included 247 infants with SPTB (gestational age [GA] < 36 weeks) and 419 term controls (GA 38–41 weeks). The strongest signal came within the gene encoding slit guidance ligand 2 (SLIT2; rs116461311, minor allele frequency 0.05, p = 1.6×10−6). Pathway analysis revealed the top-ranking pathway was axon guidance, which includes SLIT2. In 172 very preterm-born infants (GA <32 weeks), rs116461311 was clearly overrepresented (odds ratio 4.06, p = 1.55×10−7). SLIT2 variants were associated with SPTB in another European population that comprised 260 very preterm infants and 9,630 controls. To gain functional insight, we used immunohistochemistry to visualize SLIT2 and its receptor ROBO1 in placentas from spontaneous preterm and term births. Both SLIT2 and ROBO1 were located in villous and decidual trophoblasts of embryonic origin. Based on qRT-PCR, the mRNA levels of SLIT2 and ROBO1 were higher in the basal plate of SPTB placentas compared to those from term or elective preterm deliveries. In addition, in spontaneous term and preterm births, placental SLIT2 expression was correlated with variations in fetal growth. Knockdown of ROBO1 in trophoblast-derived HTR8/SVneo cells by siRNA indicated that it regulate expression of several pregnancy-specific beta-1-glycoprotein (PSG) genes and genes involved in inflammation. Our results show that the fetal SLIT2 variant and both SLIT2 and ROBO1 expression in placenta and trophoblast cells may be correlated with susceptibility to SPTB. SLIT2-ROBO1 signaling was linked with regulation of genes involved in inflammation, PSG genes, decidualization and fetal growth. We propose that this receptor-ligand couple is a component of the signaling network that promotes SPTB.

About the Authors:

Heli Tiensuu

Roles Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

Antti M. Haapalainen

Roles Formal analysis, Investigation, Methodology, Supervision, Validation, Writing – review & editing

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

ORCID logo http://orcid.org/0000-0002-9154-1552

Minna K. Karjalainen

Roles Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Visualization, Writing – original draft

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

Anu Pasanen

Roles Formal analysis, Software, Writing – review & editing

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

Johanna M. Huusko

Roles Methodology, Writing – review & editing

Affiliations PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland, Division of Human Genetics, Center for Prevention of Preterm Birth, Perinatal Institute, Cincinnati Children's Hospital Medical Center, Department of Pediatrics, University of Cincinnati College of Medicine, March of Dimes Prematurity Research Center Ohio Collaborative, Cincinnati, Ohio, United States of America

ORCID logo http://orcid.org/0000-0001-8663-9005

Riitta Marttila

Roles Data curation, Methodology, Writing – review & editing

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

Marja Ojaniemi

Roles Data curation, Methodology, Writing – review & editing

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

ORCID logo http://orcid.org/0000-0002-2660-1758

Louis J. Muglia

Roles Methodology, Writing – review & editing

Affiliation: Division of Human Genetics, Center for Prevention of Preterm Birth, Perinatal Institute, Cincinnati Children's Hospital Medical Center, Department of Pediatrics, University of Cincinnati College of Medicine, March of Dimes Prematurity Research Center Ohio Collaborative, Cincinnati, Ohio, United States of America

ORCID logo http://orcid.org/0000-0002-0301-8770

Mikko Hallman

Contributed equally to this work with: Mikko Hallman, Mika Rämet

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

Affiliation: PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland

ORCID logo http://orcid.org/0000-0002-8172-729X

Mika Rämet

Contributed equally to this work with: Mikko Hallman, Mika Rämet

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliations PEDEGO Research Unit, Medical Research Center Oulu, University of Oulu, and Department of Children and Adolescents, Oulu University Hospital, Oulu, Finland, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland

Introduction

Preterm live births that take place before 37 completed weeks of gestation and even as early as 22–24 weeks are a global problem. Up to 11.1% (15 million babies) of all births worldwide occur prematurely, and approximately 45–50% of them are idiopathic or spontaneous [1–3]. Complications caused by preterm birth are the most common cause of neonatal deaths and the largest direct cause of deaths of children <5 years of age [1,3]. The research focusing on spontaneous preterm birth (SPTB) has been complicated by etiological, pathophysiological, and genetic heterogeneities. Multiple events are associated with SPTB, either independently or in concert [4]. These include intrauterine inflammation, called chorioamnionitis, preterm premature rupture of fetal membranes (PPROM) and abnormal fetal growth relative to uterine size [5,6]. It is important to find new biomarkers for early detection of SPTB. Currently, our understanding of the early molecular pathways leading to SPTB is incomplete and there are no effective means to prevent SPTB. Knowledge of how maternal and fetal genomes contribute to the risk of SPTB could provide more personalized tools to prevent it [7].

Epidemiological studies have shown that both fetal and maternal genes affect fetal growth, birth weight, birth length, head circumference, and gestational age (GA) [8–11]. A recent study indicated that variants of the fetal and maternal genome independently affect normal variations in birth weight [12]. In addition, maternal and fetal genomes are also considered to affect the susceptibility to preterm birth and duration of pregnancy in general [13–15]. The intrauterine environment influences fetal growth, and adverse intrauterine events affect pregnancy length not only in elective preterm pregnancies but also in SPTB [11]. Genetic analysis of 244,000 Swedish births resulting in twins, full siblings, and half-siblings revealed that 13% and 21% of the variation in birth timing is explained by fetal and maternal genetic factors, respectively [9]. Overall, preterm birth is a phenotype with contributions from both, maternal and fetal genomes that may have separate contributions and together with environmental factors, interactively determine the outcome [7,13,16].

A recent study that included a population of > 40,000 women and replication cohorts of > 8000 women identified several common variants in EBF1, EEFSEC, and AGTR2 that showed associations with preterm birth at a genome-wide significance level [13]. In addition, other genome-wide association studies (GWAS) and SPTB genetic studies focused on mother [17–19] or infant [19–21] genomic signals have discovered genetic loci associated with preterm birth and gestational length. Studies that focused on fetal genomes have not revealed replicable associations between fetal genetic factors and SPTB.

Many pathways and cellular processes are reported to be associated with SPTB, including response to infection, regulation of inflammation, stress, and other immunologically mediated processes [3]. According to our current understanding, inflammatory pathways also have roles in the initiation of spontaneous term birth, as normal labor starts when there is a shift in signaling between anti-inflammatory and proinflammatory pathways in the myometrium. This shift appears to involve many chemokines such as interleukin 8 (IL8), cytokines such as IL1 and IL6, and contraction-associated proteins such as oxytocin receptor (OXTR), connexin 43 (CX43), and prostaglandin receptors [22]. Therefore, it is likely that changes in inflammation-associated pathways also contribute to preterm birth. Evidence from candidate gene studies supports the role of inflammation-related factors in SPTB. For example, polymorphisms of the genes encoding TLR4, TNF, IL1B, interferon gamma (IFNγ), IL6, and matrix metalloproteinases may be associated with increased risk of SPTB [5].

The aim of the present study was to use a GWAS to investigate fetal genetic variants that may predispose infants to SPTB in a homogeneous population of Finnish origin. A variant of slit guidance ligand 2 (SLIT2) had the most suggestive association with SPTB in the GWAS and in a genetic pathway analysis. Therefore, we characterized SPTB-associated expression of SLIT2 and its receptor ROBO1 in the placenta and subsequently conducted experiments with relevant placenta-associated cells.

Results

GWAS association signals for spontaneous preterm birth

In order to find fetal genetic factors associated with predisposition to SPTB, we analyzed polymorphisms encompassing the entire genome for associations with SPTB. After quality control, 247 infants born spontaneously preterm and 419 infants born at term remained for inclusion in the GWAS. We performed the analysis for both GA (quantitative trait) and SPTB (dichotomous setting). However, due to sample collection bias resulting in skewed GA distribution, we acknowledge that the results of the quantitative trait analysis should be interpreted with caution. Fig 1 summarizes the study workflow.

[Figure omitted. See PDF.]

Fig 1. Overview of the study workflow.

https://doi.org/10.1371/journal.pgen.1008107.g001

We detected several suggestive associations (p < 10−5) in the GWAS (Fig 2, Table 1). The two most promising regions were within the genes encoding SLIT2 (rs116461311, p = 1.6 × 10−6) and succinyl-CoA:glutarate-CoA transferase (SUGCT; rs57670997, p = 1.8 × 10−6). We also detected suggestive associations for GA as a quantitative trait (S2 Table, S1 Fig). SNP rs116461311 within the SLIT2 gene showed the most significant signal for GA (p = 3.1 × 10−7, S1 Fig). In addition to SLIT2, four regions showed suggestive signals both in the primary setting and in the GWAS of GA; these signals were within SUGCT, an intergenic region in chromosome 6 (nearest loci LOC105377949 and LOC107986634), and within the genes encoding anaplastic lymphoma receptor tyrosine kinase (ALK) and DLC1 Rho GTPase activating protein (DLC1).

[Figure omitted. See PDF.]

Fig 2. Manhattan plot of GWAS of SPTB.

Each dot represents −log10(p) value of a single SNP in association analysis. Dashed line denotes level of suggestive significance (−log10(p) > 5). SLIT2, SUGCT, and EXOSC1-ZDHHC16-MMS19-UBTD1 were among the best associating loci.

https://doi.org/10.1371/journal.pgen.1008107.g002

[Figure omitted. See PDF.]

Table 1. Single-nucleotide polymorphisms suggestively associated with spontaneous preterm birth.

https://doi.org/10.1371/journal.pgen.1008107.t001

We further analyzed very preterm and moderate-to-late preterm SPTB infants separately against term-born controls (S3 and S4 Tables). Three regions showed suggestive signals (p < 10−5) both in the primary setting and in the analysis of very preterm birth (S3 Table). These signals were within SLIT2, in an intergenic region on chromosome 2 (nearest genes THUMPD2 and SLC8A1), and within the region encompassing the EXOSC1, ZDHHC16, MMS19, and UBTD1 genes, which encode exosome component 1; zinc finger DHHC-type containing 16; MS19 homolog, cytosolic iron-sulfur assembly component; and ubiquitin domain containing 1, respectively. The minor allele of SLIT2, SNP rs116461311, was overrepresented (OR 4.06, p = 1.55 × 10−7) in very preterm-born infants (GA < 32 weeks) compared to term-born infants. In moderate-to-late preterm infants (GA 32–36 weeks), two regions showed suggestive signals (S4 Table) that were also evident in the primary analysis: an intergenic region on chromosome 3 (nearest loci LOC105377173 and ROBO1) and within ADAMTS14, which encodes ADAM metallopeptidase with thrombospondin type 1 motif 14.

We also studied associations with SPTB within the contexts of PPROM and no PPROM. There were separate suggestive associations with SPTB-PPROM and with SPTB without PPROM (S5 and S6 Tables). For the SLIT2 region, the effects were similar for infants born after PPROM (rs116461311, OR = 3.5, p = 3.0 × 10−5) and for those born after spontaneous onset of labor with intact fetal membranes (rs116461311, OR = 3.6, p = 1.3 × 10−5).

To investigate potential maternal transmission of the minor allele of SPTB-associated rs116461311 (in SLIT2), we checked the MAF of the variant in maternal samples. The frequency of C-rs116461311 was 0.059 in SPTB mothers (n = 230) and 0.035 in mothers with term delivery (n = 378). Transmission analysis was not feasible because of the low minor allele counts.

Replication of SLIT2 association in a European population

SNPs with suggestive association signals in the region of the top GWAS gene (SLIT2) were examined for association with SPTB in a European population (n = 9890) [20]. In a European replication population of very preterm and term-born controls, rs12503652 and rs79034379, which correlate with the best GWAS SNP of the SLIT2 region (rs116461311), were associated with SPTB (Table 2). The two SNPs were in high LD with one another (r2 = 0.85, D′ = 0.95) in the Finnish individuals in 1000 Genomes phase3 data.

[Figure omitted. See PDF.]

Table 2. Two SNPs that correlated with the top GWAS SNP rs116461311 in SLIT2 showed association with SPTB in European replication data.

https://doi.org/10.1371/journal.pgen.1008107.t002

Functional annotation of the most promising genetic variants

To investigate if the SNPs within regions with the most promising signals in the GWAS (p < 10−5) have functional consequences, we screened the GTEx data to examine whether any of the suggestively associated SNPs colocalize with cis eQTLs; that is, whether they correlate with mRNA levels in the analyzed tissues. In the GTEx data, 26 of the SNPs with p < 10−5 for association with SPTB overlapped with significant eQTLs (Table 3). These SNPs were located within five different regions; the majority were in the region encompassing EXOSC1, ZDHHC16, MMS19, and UBTD1 (Tables 1, S3 and S5). These SNPs correlated with mRNA levels of the following genes in different tissues: ZDHHC16, MMS19, FRAT1, ANKRD2, UBTD1, and RRP12. Associated SNPs within the top GWAS region, SLIT2, were not associated with mRNA levels. According to the GWAS catalog, none of the suggestively associated SNPs had been significantly associated with any phenotype.

[Figure omitted. See PDF.]

Table 3. Suggestively associated SNPs that overlap with cis eQTLs.

https://doi.org/10.1371/journal.pgen.1008107.t003

We further functionally annotated the SNPs (p < 10−4) within the SLIT2 region with HaploReg, v 4.1. Some of the most promising SNPs within SLIT2 (including the best associating variant rs116461311, as well as rs60126904, rs115707845, and rs16869667) were located within regions that contain histone marks and DNase-hypersensitive sites in several tissues. Furthermore, SLIT2 SNPs (rs60126904 and rs115707845) mapped to predicted enhancers in several cell types and tissues, including neuronal cells, different cells in the brain, lung, spleen, adipose cells, colon and duodenum smooth muscle cells, fetal lung and kidney, and fetal membranes (amnion). Thus, there is some evidence of the putative regulatory effects of SLIT2 SNPs in several tissues, including fetal tissues such as placenta.

Biological pathways genetically associated with spontaneous preterm birth

To identify biological pathways associated with SPTB, we performed pathway analysis of the GWAS data to search for gene set enrichment. SPTB was associated with 16 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [corrected (S7 Table) and FDR adjusted p < 0.05, Table 4]. The most significant pathways were axon guidance (p = 8.6 × 10−10), focal adhesion (p = 6.6 × 10−7), and vascular smooth muscle contraction (p = 1.4 × 10−6). Axon guidance, the most significant pathway, included SLIT2 and the gene encoding its receptor, ROBO1. A Gene Ontology (GO) search revealed 35 GO terms associated with SPTB, with a false-discovery rate (FDR) of <0.05 (S8 Table). GO sets that included SLIT2 and ROBO1 are listed in S9 Table. The three most significant GO sets that included SLIT2 were retinal ganglion cell axon guidance, telencephalon development, and negative chemotaxis. The three most significant GO sets that included ROBO1 were telencephalon development, neuron recognition, and negative chemotaxis. These results led us to a more-detailed investigation of the roles of SLIT2 and ROBO1 in SPTB.

[Figure omitted. See PDF.]

Table 4. Pathways associated with SPTB.

https://doi.org/10.1371/journal.pgen.1008107.t004

SLIT2 and its receptor ROBO1 are localized in villous and decidual trophoblasts

We detected suggestive association signals for SNPs in a region that encompasses SLIT2 and in a region downstream of ROBO1. Protein Slit2 binds to Robo proteins specifically and with high affinity [23,24]. Therefore, we analyzed the localization of SLIT2 and its receptor ROBO1 in human placenta by immunohistochemical staining of placentas from SPTB and spontaneous term birth (STB) with anti‐human SLIT2 and ROBO1 antibodies (Fig 3). Both SLIT2 and ROBO1 localized to cytotrophoblasts, syncytiotrophoblasts, and decidual trophoblasts. In addition, we observed strong ROBO1 and faint SLIT2 staining in capillary endothelial cells. We also detected both proteins in the basal and chorionic plates of the placenta. We did not see apparent differences in staining intensities or cellular localization between placentas from SPTB and STB. This indicates that SLIT2 and ROBO1 are expressed in the placenta at the interface between mother and fetus during pregnancy.

[Figure omitted. See PDF.]

Fig 3. Placental localization of SLIT2 and ROBO1 in spontaneous preterm and term births.

Samples from 18 placentas were immunostained with anti‐human SLIT1 (A) and ROBO1 (B) antibodies. Twelve of the placentas were from SPTB deliveries, and six were from spontaneous term deliveries. Two samples from each placenta were stained; one sample from the basal plate (maternal side) and the other from the chorionic plate (fetal side). Immunostaining is indicated by filled large arrows in cytotrophoblasts, unfilled large arrows in syncytiotrophoblasts, filled small arrows in decidual trophoblasts, and filled arrowheads in capillary endothelial cells. Original magnification is 20× in all figures. Control represents isotype controls for immunostaining. Scale bar, 100 μm.

https://doi.org/10.1371/journal.pgen.1008107.g003

SLIT2 and ROBO1 mRNA level elevation in SPTB placentas

Immunohistochemistry demonstrated SLIT2 and ROBO1 in different types of trophoblasts in human placenta (Fig 3). To obtain more quantitative data about placental expression of these proteins, we analyzed SLIT2 and ROBO1 mRNA levels by qRT-PCR in samples collected from the basal and chorionic plates of placentas from SPTB (n = 23), STB (n = 23), and elective preterm birth (EPTB) (n = 34).

We first compared SLIT2 and ROBO1 expression levels between SPTB (n = 23) and STB (n = 23) placentas (Fig 4). Both SLIT2 and ROBO1 mRNA levels were higher in the basal plate of SPTB placentas (SLIT2 fold change [FC] = 1.679, SD = 0.667; ROBO1 FC = 1.387, SD = 0.670) compared to those of STB (SLIT2 p = 0.004, ROBO1 p = 0.013; Fig 4A). There were no differences in mRNA levels of SLIT2 (p = 0.173) and ROBO1 (p = 0.297) between the chorionic plates of SPTB and STB placentas (Fig 4B).

[Figure omitted. See PDF.]

Fig 4. mRNA levels of SLIT2 and ROBO1 in spontaneous preterm, spontaneous term, and elective term placentas.

Relative mRNA levels from the basal (A) or chorionic (B) plates were normalized to mRNA levels of the housekeeping gene CYC1. Differences were analyzed with nonparametric Mann–Whitney U‐test. Asterisk indicates statistically significant changes. Box and whiskers display quartiles. Band inside the box is the median, and ends of whiskers show the minimum and maximum values; square inside the box represents the mean value.

https://doi.org/10.1371/journal.pgen.1008107.g004

To explore the effects of mode of delivery and GA on SLIT2 and ROBO1 mRNA levels, we compared SPTB with EPTB placentas. SLIT2 and ROBO1 mRNA levels were significantly higher for SPTB (SLIT2 FC = 1.595, SD = 0.580; ROBO1 FC = 1.282, SD = 0.577) compared to EPTB (SLIT2 p = 0.005, ROBO1 p = 0.031) in basal plate samples (Fig 4A). There were no significant differences in these levels in EPTB and STB placentas (SLIT2 p = 0.216, ROBO1 p = 0.328; Fig 4A). These results suggest that higher SLIT2 and ROBO1 expression levels are associated with SPTB.

Association of SLIT2 expression with fetal growth

SLIT/ROBO are involved in many processes that involve cell migration, including axon guidance; thus, they could affect trophoblast cell invasion and decidualization. To this end, we looked at whether SLIT2 or ROBO1 expression in the basal plate of the placenta is associated with fetal growth. We compared mRNA levels with birth weight-for-GA Z-scores (weight Z-score), which included age and gender standardization of the infants. Deliveries with intrauterine growth restriction or other growth disorders were excluded. SLIT2 mRNA levels correlated with Z-scores (p = 0.023, rs = 0.351) in term and preterm fetuses delivered after spontaneous onset of labor (SPTB and STB samples together, S3 Fig). This suggests that SLIT2 expression is associated with variations in fetal growth.

siRNA-induced knockdown of ROBO1 causes changes in expression levels of inflammation related factors

To investigate potential functions of SLIT2 and ROBO1 in placental cells, we silenced SLIT2 and ROBO1 expression separately in the HTR‐8/SVneo human trophoblast cell line with small interfering RNAs (siRNAs) (Fig 5). qRT-PCR revealed that silencing percentages were 60% and 85%, for SLIT2 and ROBO1, respectively. Corresponding percentages revealed by RNA sequencing were 75% and 74% for SLIT2 and ROBO1, respectively. Next, we characterized the transcriptomes of trophoblasts in which SLIT2 or ROBO1 was silenced, as well as of cells treated with siRNA Universal Negative Control #1. Transcriptomic data analysis identified 14 upregulated (S10 Table) and 12 downregulated (S11 Table) genes after SLIT2 knockdown compared to samples treated with siRNA Universal Negative Control. The threshold was an FDR-adjusted p value of ≤0.01 and an FC of ≥2.0. By the same criteria, there were 216 upregulated (S12 Table) and 610 downregulated (S13 Table) genes after ROBO1 knockdown. KEGG pathway database analyses (S14 Table, S15 Table) identified the top pathways affected after SLIT2 and ROBO1 knockdown as inflammation-related pathways such as cytokine-cytokine receptor interaction (KEGG.ID 4060).

[Figure omitted. See PDF.]

Fig 5. mRNA expression level changes after gene knockdown.

SLIT2 (A) and ROBO1 (B) were post-transcriptionally knocked down with siRNA in human placental trophoblast continuous cell line HTR8/SVneo. mRNA levels of cells in which SLIT2 or ROBO1 were knocked down was compared with mRNA levels from siRNA negative treated control cells. Expression levels of the genes were first determined by qRT-PCR (levels normalized against housekeeping gene CYC1 mRNA levels). According to qRT-PCR, silencing percentage of SLIT2 was 60% and silencing percentage of ROBO1 was 85%. Gene expression levels also determined with high-throughput RNA sequencing. Reads per kilobase of exon per million reads mapped (RKPM) is a normalized gene count value determined in transcriptomic analysis. According to transcriptomics data, silencing percentages were 75% for SLIT2 and 74% for ROBO1. Columns represent median and SD values.

https://doi.org/10.1371/journal.pgen.1008107.g005

Far fewer genes were affected by SLIT2 knockdown than by ROBO1 knockdown, probably because transfection reagent alone upregulated SLIT2 mRNA levels up to 4-fold (p = 0.002) compared to untreated cells. Therefore, knockdown of SLIT2 expression brought SLIT2 mRNA levels close to the levels of intact cells. Transfection reagent by itself seemed to activate inflammation-related pathways. The ROBO1 expression level was not affected by transfection reagent.

ROBO1 knockdown particularly affected genes encoding membrane receptors and other membrane proteins. KEGG pathway analysis revealed that hematopoietic cell lineage (KEGG ID 4640) had the lowest p value (6.57 × 10−8) (S15 Table). ROBO1 knockdown affected 14 of the 42 genes in this pathway: KIT, IL7R, IL1R1, HLA-DRB1, IL1A, ITGB3, TFRC, ANPEP, IL1B, CD22, CD24, KITLG, CSF3, and CD14. SLIT2 (FC = 2.0, p = 0.005) was among the genes upregulated after ROBO1 knockdown (S12 Table). One of the gene families highly affected by ROBO1 knockdown was pregnancy-specific glycoproteins (PSG), a complex gene family that regulates maternal–fetal interactions [25]. Of the ten protein-coding human PSG genes, six were upregulated after ROBO1 knockdown. These data indicate that ROBO1 is an important regulator of the HTR8/SVneo cell transcriptome and suggest a role for ROBO1 in modulation of PSG gene expression. In addition, ROBO1 appears to have immunomodulatory functions in trophoblast-derived cells.

To verify our suggestive findings from the RNA sequencing data, we used qRT-PCR to analyze the effect of ROBO1 knockdown on expression levels of selected genes from a larger number of specimens. Because many inflammation‐related pathways were involved, we investigated TGFA, CXCL6, and total PSG expression levels (Table 5). TGFA was downregulated after ROBO1 knockdown, while CXCL6 was downregulated by both ROBO1 and SLIT2 knockdown. Select members of the PSG family were upregulated when ROBO1 was silenced (Table 5): PSG1, PSG2, PSG4, PSG6, PSG7, and PSG9. To verify, we measured total mRNA expression of different PSGs (PSG1, PSG2, PSG3, PSG4, PSG5, PSG6, PSG7, PSG8, PSG9, and PSG11) at the same time in one qRT-PCR reaction, as described previously [26]. We also tested the effect of knockdown on IL6, TNFA, and SRGAP3 (Table 5).

[Figure omitted. See PDF.]

Table 5. Comparison of RNA sequencing and qRT-PCR of ROBO1 and SLIT2 genes affected by knockdown.

Gene expression changes after SLIT2/ROBO1 knockdown were verified in a larger number of specimens (n = 6) by qRT-PCR. SLIT2 and ROBO1 were silenced in HTR8/SVneo commercial cell line with siRNAs. mRNA levels of SLIT2/ROBO1-silenced cells were compared with mRNA levels of cells treated with negative siRNA. In addition to genes affected by gene knockdown according to RNA sequencing data, we tested the effect of knockdown on mRNA expression of well-known cytokines IL-6 and TNF-A.

https://doi.org/10.1371/journal.pgen.1008107.t005

The results of qRT-PCR were generally in line with those of transcriptome sequencing (Table 5). Knockdown of ROBO1 downregulated mRNA expression of TGFA (S4A Fig). Knockdown of either ROBO1 or SLIT2 downregulated CXCL6, whereas PSGs were upregulated by ROBO1 knockdown, similar to the results of RNA sequencing. Inflammatory cytokines IL6 and TNFA were both downregulated after ROBO1 knockdown (S4B Fig). IL6 was also downregulated after SLIT2 knockdown. There was a trend toward downregulation of SRGAP3 by ROBO1 knockdown. Thus, the RNA sequencing and qRT-PCR results are in accordance and confirm that ROBO1 knockdown affects the expression of immune response–modifying genes in a cell culture model.

Discussion

Both maternal and fetal genetic factors contribute to SPTB. Nevertheless, variants in the fetal genome associated with SPTB predisposition are not well known. Our aim was to identify fetal genetic factors associated with SPTB. To this end, we performed a case–control GWAS study in a Finnish population that is known to be relatively genetically homogeneous. Based on these results, we identified a SLIT2 variant as a plausible factor for SPTB susceptibility. This led to further investigations to define high-risk populations and characterize SLIT2 and ROBO1 expression in the placenta and in trophoblast cells.

The association of the SLIT2 variant with SPTB was strongest in the population of very premature births (<32 weeks gestation), and this association was replicated in a European population with SPTB fetuses from 24 to 30 weeks gestation. We detected association signals for SNPs in a region that encompassed SLIT2 and for SNPs in a region downstream of ROBO1, which encodes the receptor for SLIT2. SLIT2 and ROBO1 encode proteins of the SLT2-ROBO1 signaling pathway; previous studies have indicated that this pathway is associated with different types of pregnancy complications, including preeclampsia [27,28], impaired placentation of missed and threatened miscarriage in early pregnancy [29], trophoblast invasion, and vascular remodeling during ectopic tubal pregnancy [29,30]. Our hypothesis-free GWAS study provides evidence that these genes have a role in another pregnancy complication, SPTB.

We did not detect associations that would reach the stringent level of genome-wide significance (p<10−8). This may be due to a relatively small sample size in the GWAS, which was one of the limitations in the study. However, we did identify signals below the generally used threshold of suggestive association (p <10−5). Moreover, the SPTB-associated SNPs in the SLIT2 locus showed association in an independent data set. In the future, there is a need to validate our findings in larger data sets to detect signals that may have gone undetected in the current sample size.

Previous study presented that fetal de novo mutations in genes that are involved in brain development are associated with preterm birth [31]. In line with this notion, SLIT2-ROBO1 signaling has a well-documented role in axon guidance during the development of the nervous system [32]. Therefore, it appears that at least in part the same fetal genetic factors are involved both in the onset of preterm birth and in brain development. It is known that preterm birth increases the risk of compromised brain development [33,34]. How much of the risk can be explained by shared genetic risk factors remains to be determined.

KEGG pathway analysis of GWAS data showed that the top-ranking pathway was SLIT2/ROBO1 signaling–regulated axon guidance. Previous studies have indicated that brain and placental development may share common pathways [35–37]. Although SLIT and ROBO were originally identified as axon guidance cues, they interact in many other cellular processes, including regulation of cell migration, cell death, and angiogenesis. As such, they have an essential role in the development of tissues such as lung, liver, kidney, breast, and tissues of the reproductive system [38,39].

Trophoblasts cover a large portion of the placenta and have multiple roles in the maintenance of pregnancy. Invading trophoblasts have a critical function in biogenesis of the placenta [40,41]. Later during pregnancy, decidual trophoblasts may have a role in silencing immune cells in the decidua. The lining of placental villae consists of the syncytiotrophoblast layer and cytotrophoblasts, which have roles ranging from immune protection to uptake of nutrients from maternal blood [42]. In addition to their functions in biogenesis of the placenta and maintenance of pregnancy, trophoblasts may participate in labor induction [43]. Two‐thirds of early pregnancy failures may present with reduced trophoblast invasion [29]. SLIT-ROBO signaling may play autocrine and/or paracrine roles in trophoblast functions, such as differentiation and invasion, by influencing the migration of trophoblastic cells [29,44]. Thus, SLIT2-ROBO1 signaling may be involved in the pathogenesis of pregnancy failures via their effect on trophoblastic cell functions. Indeed, our immunohistochemical experiments demonstrated that villous and decidual trophoblasts from preterm and term placentas were strongly positive for SLIT2 and ROBO1. These results are in line with those of a previous study, which demonstrated that villous syncytiotrophoblasts express high levels of SLIT2 and ROBO1. The study also found that trophoblastic endothelial cells highly coexpress multiple SLIT ligands (SLIT2, SLIT3) and ROBO receptors (ROBO1, ROBO2, and ROBO4) in full-term placenta. Thus, SLIT-ROBO signaling may also have an important role in the regulation of normal placental functions [27].

An earlier study found that levels of several SLIT/ROBO mRNAs and proteins are higher in preeclamptic placentas compared to normal controls [27]. Our data indicate that mean SLIT2 and ROBO1 mRNA levels were higher in SPTB placentas compared to placentas from spontaneous term deliveries and to placentas from elective preterm births. Consequently, increased levels of SLIT2 and ROBO1 mRNAs were associated with SPTB. In addition, SLIT2 may have a role in term spontaneous labor, as SLIT2 mRNA and protein expression are decreased in the myometrium after spontaneous term labor [45].

In addition to the associations of both the SLIT2 variant and mRNA expression of SLIT2 and ROBO1 in basal plate of placenta with SPTB, SLIT2 mRNA levels in placentas were associated with the birth weight of fetuses born after spontaneous labor. A GWAS of beef cattle identified SLIT2 as a candidate gene that affects the weight of internal organs [46]. As the fetal growth and intrauterine distention negatively associates the duration of pregnancy [47] influence of SLIT2 on fetal size is a potential mechanism that remains to be studied as a cause of SPTB.

To understand more about the function of SLIT2 and ROBO1 in trophoblast cells, we silenced their expression separately in immortalized extravillous invading trophoblasts. Altogether, 26 and 826 genes were affected by SLIT2 or ROBO1 siRNA knockdown, respectively. The low number of genes affected by SLIT2 knockdown was probably because the transfection reagent alone upregulated SLIT2 mRNA levels compared to untreated cells. Consequently, knockdown of SLIT2 expression only brought SLIT2 mRNA levels back to the levels of control cells. However, the mRNA expression level of ROBO1 was not affected by the transfection reagent. Genes affected by ROBO1 knockdown were mostly related to infection, inflammation, and immune response. These results correspond with those of previous studies, suggested that members of the SLIT and ROBO families act as regulators of the inflammatory response [45,48]. Both pro‐inflammatory [45,49,50] and anti‐inflammatory [48,51] functions have been reported. Our results from invading trophoblast cells support a proinflammatory role for SLIT2-ROBO1 signaling, since the genes downregulated by ROBO1 knockdown included proinflammatory cytokines and chemokines such as IL1A, IL1B, CXCL8, CCL2, and CXCL6. It is widely acknowledged that IL1 in particular, as well as other proinflammatory cytokines and chemokines, is associated with preterm labor [52–54]. IL1B is a primary secretory product of the inflammasome and as such is thought to have central roles in initiation of preterm labor, such as in the induction of prostaglandin synthesis. Both polymorphisms of IL1A and IL1B, as well as increased levels of IL1B, are associated with preterm birth [19,54–56]. CXCL6 is increased in amniotic fluid from patients with preterm labor complicated by intra-amniotic inflammation and from patients with SPTB without intra-amniotic infection/inflammation [57]. We propose that in trophoblast cells ROBO1 has a role in regulation of proinflammatory mediators. Genes involved in vascular formation (vasculogenesis) or development (angiogenesis) were not affected by SLIT2 or ROBO1 knockdown in trophoblasts.

The PSG family was one of the immune response–associated gene families affected by ROBO1 knockdown. PSGs include ten placental trophoblast–synthetized glycoproteins that belong to the immunoglobulin superfamily [58]. Of the six PSGs upregulated by siRNA-induced knockdown, PSG1 was ranked among the top three upregulated genes (S12 Table). PSGs are essential in the maintenance of normal pregnancy [58]; thus, altered PSG expression patterns could influence pregnancy complications. Over the years, complications such as abortion, preeclampsia, intrauterine growth retardation, fetal distress, and preterm delivery have all been linked to low PSG levels [58–63]. As ROBO1 was upregulated in SPTB placentas and knockdown of ROBO1 upregulated expression of PSG genes, we propose that ROBO1 signaling is important in downregulation of the expression of PSGs. In addition, PSG1 activates TGF-B 1 and TGFB2 [25,64,65]; TGFB1 suggestively associated with SPTB [13] and TGFB2 prevented preterm birth in experimental inflammatory stress [66].

The innate immune response and inflammation contribute to labor and delivery, particularly in preterm pregnancies [3,22,67,68]. Upregulation of proinflammatory cytokines stimulates and potentiates uterine contractions in the myometrium [69–71]. In preterm labor and delivery, it is mostly inflammatory signals that spread to the placenta, fetal membranes, and fetal compartment. It is plausible that SLIT2 and ROBO1 expressed by trophoblasts are associated with SPTB via regulation of inflammation-related factors. SLIT2-ROBO1–guided activation and propagation of inflammatory mediators throughout the fetal–maternal trophoblast interface of the uterine wall would likely influence the tissues actively involved in labor and delivery. As knockdown of ROBO1 downregulated many of the genes that encode cytokines and chemokines, it is probable that upregulation of ROBO1 in SPTB placentas compared to term placentas would also affect expression of these genes. There is both epidemiological and experimental evidence that untimely expression of cytokines and chemokines by either fetal or maternal tissues upregulates the activity of mediators, which leads to premature initiation of the parturition process [72].

In conclusion, the GWAS detected fetal association signals for SPTB and duration of pregnancy in the vicinity of SLIT2 and ROBO1. SLIT2 and ROBO1 were upregulated in SPTB placentas, and further functional studies confirmed that this signaling pathway has a role in regulation of the pathways associated with infection, inflammation, and immune response in trophoblast-derived cells. These results suggest that SLIT2 and ROBO1 play specific roles in increasing susceptibility to SPTB. SLIT2-ROBO1 signaling is associated with complications in early pregnancy and it is possible that it influences invading trophoblasts during placentation. Based on the currently available evidence, we propose that activation of SLIT2-ROBO1 expression and signaling in trophoblast cells contributes to inflammatory and immune activation, which in turn leads to early labor and preterm birth.

Materials and methods

Ethical approval

The present studies received ethical approval from the participating centers (Oulu University Hospital 79/2003, 14/2010, and 73/2013). Informed consent was obtained from study participants or their parents.

Study populations in genetic studies

Characteristics of the Finnish study populations are summarized in S1 Table. The discovery GWAS study population consisted of singleton SPTB and term infants sampled in Oulu and Tampere University Hospitals. The study subjects were recruited prospectively during 2004–2014 and retrospectively from the 1973–2003 birth diaries of Oulu University Hospital. For replication, we downloaded the summary statistics of a European population described in a recent study by Rappoport et al. [20] through ImmPort (http://www.immport.org/: SDY1205, DOI: 10.21430/M37N6PJEQT). This population includes 260 SPTB cases (139 male and 121 female infants) and 9,630 controls (4,055 males and 5,575 females). The cases were very preterm infants born between 25 and 30 weeks of gestation and were clinically defined as SPTB in 2005–2008. The control population consisted of adults, originally from the Health and Retirement Study (HRS) [73], who were matched for ethnicity with the European cases.

In the Finnish cohorts, SPTB was defined as birth occurring after spontaneous onset of labor at <36 completed weeks + 1 day of gestation. All medically indicated preterm births and deliveries that included known major risk factors were excluded. These criteria led to exclusion of preterm deliveries that involved the following conditions or characteristics: multiple gestation, preeclampsia, intrauterine growth restriction, placental abruption, polyhydramnios, fetuses with anomalies, clinical chorioamnionitis or acute septic infection in the mother, diseases in the mother that could influence timing of delivery, alcohol/narcotic use, and accidents. Term birth was defined as birth occurring at 38–41 weeks (38 wk + 0 d to 41 wk + 6 d) of gestation. The following conditions were used as exclusion criteria for the control population: multiple gestation, intrauterine growth restriction, placental abruption, polyhydramnios, fetuses with anomalies, and requirements for special care of the newborn. All control infants were from families with at least two term deliveries without any preterm deliveries in the family.

DNA sample preparation, genotyping, and SNP inclusion criteria

Umbilical cord blood, umbilical cord tissue, or saliva was obtained from the study subjects. Commercial kits were used to extract genomic DNA from blood (UltraClean Blood DNA Isolation Kit; MO BIO Laboratories, Inc., Carlsbad, CA, USA or Puregene Blood Core Kit; Qiagen, Hilden, Germany) and cord tissue (Gentra Puregene Tissue Kit, Qiagen). OraGene DNA collection kits (DNA Genotek, Ontario, Canada) were used for collecting saliva, and DNA was extracted with the prepIT-L2P kit (DNA Genotek). Genome-wide SNP genotyping was performed with the Infinium HumanCoreExome BeadChip (Illumina, San Diego, CA, USA) by the Technology Centre, Institute for Molecular Medicine Finland (FIMM), University of Helsinki.

Processing of genome-wide SNP data

Genome-wide SNP data were processed with PLINK, v. 1.9 [74]. SNPs with minor allele frequency (MAF) < 0.01, genotyping rate < 0.9, or deviation from Hardy–Weinberg equilibrium (p < 0.0001) were excluded. Individuals with > 0.1 missing genotypes were excluded. Identical by descent (IBD) clustering and multidimensional scaling (MDS) analyses were performed with a linkage disequilibrium–pruned SNP set; population outliers and close relatives (pihat > 0.2) were excluded. Prephasing of genotypes was performed with SHAPEIT2 [75], followed by statistical imputation with IMPUTE2 [76] using the 1000 Genomes Phase 3 variant set (October 2014) as the reference panel. Before association analysis, SNPs with impute info score < 0.8 or MAF < 0.05 in cases or controls were excluded. Altogether 6,778,521 SNPs or short insertions/deletions remained for analysis after these quality control steps.

Statistical analyses in genetic studies

Associations between SPTB or GA and SNPs were assessed with the frequentist test under the additive model with SNPtest, v. 2.5.2 [77]. After the primary analysis, the following subgroups of SPTB infants were assessed: (1) very preterm infants (GA 23–31 wk + 6 d), (2) moderate-to-late SPTB infants (GA 32 wk + 0 d to 36 wk + 0 d), (3) PPROM before onset of labor, and (4) no PPROM before onset of labor. To account for population substructure in the GWAS, the first two MDS dimensions were included as covariates. In the GWAS, the effect of population stratification was minimal (λ = 1.03). Gene set analysis (GSA)-SNP was used to search for gene set enrichment in pathway analysis [78]. We included only genotyped SNPs located within genes in this analysis to avoid the complicating effects of SNPs in linkage disequilibrium. R, v. 3.2.2 (https://www.r-project.org) was used to create Manhattan plots. LocusZoom [79] was used to create regional association plots.

Functional annotation of SNPs

We annotated SNPs with three approaches: (1) We used Genotype-Tissue Expression (GTEx) data to analyze whether the SNPs overlap with cis expression quantitative trait loci (eQTLs) [80]; 2) we screened whether the SNPs had been associated with any phenotypes in previous GWA studies using the GWAS catalog[81]; and (3) we assessed whether the SNPs were located within putative regulatory regions using HaploReg, v. 4.1 [82].

Placental samples for immunohistochemistry and qRT-PCR

Samples from human placenta were collected at Oulu University Hospital during 2010–2016 as described [16]. The placental samples used in immunohistochemical staining and in qRT-PCR analysis of SLIT2 and ROBO1 expression were subject to similar inclusion criteria as the samples used in GWAS. The inclusion criteria of gestational age for preterm placental samples was from 25 weeks to 36 weeks+ 6 days and 39 weeks to 41 weeks + 6 days for term samples. The same conditions (multiple gestation, intrauterine growth restriction, placental abruption, polyhydramnios, fetuses with anomalies or requirements for special care of the newborn) as in GWAS were used as exclusion criteria for term controls. Spontaneous preterm samples had almost the same exclusion criteria except the population included few cases with chorioamnionitis or oligohydramnion. The control group of elective preterm samples included cases with various pregnancy complications like IUGR or pre-eclampsia resulting in elective preterm delivery without labor.

Specifically, in total, 18 placental samples were analyzed by immunohistochemistry. Twelve samples were from SPTB deliveries (GA from 25 wk + 2 d to 35 wk + 2 d), and six were from spontaneous term deliveries (GA from 39 wk + 4 d to 41 wk +1 d). Samples from both basal and chorionic plates were included in the study.

RT-qPCR was performed with 23 placental samples from SPTB (GA from 25 wk + 2 d to 36 wk + 0 d), 34 from elective preterm birth (EPTB) (GA from 25 wk + 1 d to 36 wk + 6 d), and 23 from spontaneous term birth (STB) (GA from 39 wk + 1 d to 41 wk + 6 d).

Sample preparation and immunohistochemical staining

Localization of encoded proteins was visualized in placental tissues by immunohistochemical staining. Samples were embedded in paraffin and cut into 4-μm slices, deparaffinized, and rehydrated. Antigen retrieval was done in Tris-EDTA buffer. Endogenous peroxidase activity was blocked in blocking solution (Agilent, Santa Clara, CA, USA). Samples from the chorionic plate were incubated with mouse anti-human SLIT2 antibody (1:4000 dilution, PA5-3113; ThermoFisher Scientific, Waltham, Massachusetts, USA) or mouse anti-human ROBO1 antibody (1:2000, PA5-34931; ThermoFisher Scientific). Samples from the basal plate of the placenta were incubated in a 1:5000 dilution of mouse anti-human SLIT2 antibody and 1:1000 dilution of mouse anti-human ROBO1 antibody. Bound antibodies were detected with the Envision kit (Agilent).

Quantitative PCR

Tissue samples were homogenized, RNA was isolated with the RNeasy Mini Kit (Qiagen), and cDNA was synthetized as described previously [16]. After the RT-PCR, cDNA samples were diluted 1:2 using Rnase-free H20.

SLIT2 and ROBO1 were relatively quantified by intron spanning assays with Light-Cycler96 (Roche Diagnostics, Risch-Rotkreuz, Switzerland) and cytochrome c-1 (CYC1) as a reference gene. CYC1 was chosen as a reference gene because it is one of the most stably expressed genes in the placenta [16,83–85]. Primers and probes were: forward 5′-CTTCCAGAGACCATCACAGAAA-3′ and reverse 5′-CGTCTAAGCTTTTTATATGGTGAGAA-3′ for SLIT2 (with UPL probe #79), forward 5′-CGCAGAGAAACCTACACAGATG-3′ and reverse 5′-GGATTGGGCAGTAGGTGACT-3′ for ROBO1 (with UPL probe #31), and forward 5′-ATAAAGCGGCACAAGTGGTCA-3′ and reverse 5′-GATGGCTCTTGGGCTTGAGG-3′ for CYC1 (with UPL probe #47). Probes were from the Universal Probe library (UPL) Set (Roche Diagnostics). Each qPCR measurement was done in triplicate. Levels of SLIT2 and ROBO1 were normalized against the CYC1 level, and relative quantifications were then assessed with the ΔΔ cycle threshold method. A few randomly chosen qPCR products were also verified by agarose gel electrophoresis and Sanger sequencing.

Primers and probes genes for transforming growth factor alpha (TGFA), C-X-C motif chemokine ligand 6 (CXCL6), and SLIT-ROBO Rho GTPase activating protein 3 (SRGAP3) were: forward 5′-CCCTGGCTGTCCTTATCATC-3′ and reverse 5′-GGCACCACTCACAGTGTTTTC-3′ for TGFA (with UPL probe #74), forward 5′-CCAGAAAATTTTGGACAGTGG-3′ and reverse 5′-GGGATCTCCAGAAAACTGCTC-3′ for CXCL6 (with UPL probe #61), and forward 5′-GAAGGGCACTCGATGAGGT-3′ and reverse 5′-GCTCATGGTCTTCTCGATGTC-3′ for SRGAP3 (with UPL probe #66).

Total mRNA expression of different PSGs was measured with PCR primers: forward 5′-CCTCTCAGCCCCTCCCTG-3′ and reverse 5′-GGCAAATTGTGGACAAGTAGAAGA-3′ (with UPL probe #15), which are complementary to sequences conserved in all but two PSG transcript variants that both lack the N domain [26]. IL6 and TNFA were analyzed with primers forward 5′-GCCCAGCTATGAACTCCTTCT-3′ and reverse 5′-GCGGCTACATCTTTGGAATC-3′ for IL-6 (with UPL probe #43) and forward 5′-CAGCCTCTTCTCCTTCCGAT-3′ and reverse 5′-GCCAGAGGGCTGATTAGAGA-3′ for TNFA (with UPL probe #40).

Gene knockdown of SLIT2 and ROBO1 by transfection with small interfering RNAs (siRNAs)

Human placental trophoblast cells HTR-8/SVneo (CRL-3271™; ATCC, Manassas, Virginia, USA). were grown in RPMI-1640 culture media (Thermo Fisher Scientific) supplemented with 10% fetal bovine serum (FBS; Sigma-Aldrich, St. Louis, MO, USA) and 1× penicillin/streptomycin (Sigma-Aldrich). Cells were cultured under standard culturing conditions (37°C, 5% CO2, humidified atmosphere), and subculturing was performed with 0.05% trypsin/0.02% EDTA.

siRNAs targeting SLIT2 (s GUCAUAUCAAGAACUGUGAdTdT, a UCACAGUUCUUGAUAUGACdTdT) and ROBO1 (s CAUACCUAUGGCUACAUUUdTdT, a AAAUGUAGCCAUAGGUAUGdTdT) (Sigma-Aldrich) were reverse transfected and then forward transfected in HTR-8/SVneo cells with Lipofectamine 3000 reagent (Invitrogen, Carlsbad, CA, USA) [86]. MISSION siRNA Universal Negative Control #1 (Sigma-Aldrich) was used as a negative control for siRNA transfection and was transfected in the same manner as siRNAs targeting SLIT2 and ROBO1. In the reverse transfection, the cells (70,000 cells/well) were incubated with siRNAs at a final concentration of 30 nM. The cells were transfected again after 24 h of incubation. The second transfection was done as a forward transfection in the presence of 40 nM siRNAs. Cells were incubated with siRNAs for 24 h after the second transfection, and then fresh medium was added and cells were incubated for an additional 24 h. Cells were harvested with 1× Trypsin-EDTA (Sigma-Aldrich).

Transcriptomic analysis of SLIT2- and ROBO1-silenced HTR8/SVneo cells

Cells were disrupted with a 25 G needle and 1 ml syringe, and RNA was isolated in accordance with the manufacturer’s instructions (RNeasy Micro Kit, Qiagen). The quality of isolated RNA was checked with an Agilent 2100 Bioanalyzer system at the Biocenter Oulu Sequencing Center, Finland.

Samples containing total RNA were sent for transcriptomic analysis to the Finnish Functional Genomics Centre (FFGC; Turku, Finland), where transcriptional profiles of SLIT2 (n = 3) and ROBO1 (n = 3) -silenced cells and negative-control cells (n = 3) were detected with the Illumina HiSeq high‐throughput sequencing system. Whole cell RNA sequencing data was analyzed by the Bioinformatics Unit Core Service at the Turku Centre for Biotechnology, Finland. The transcriptomics data have been deposited in NCBI’s Gene Expression Omnibus [87] and are accessible through GEO Series accession number GSE119101 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE119101)

qRT-PCR verification of transcriptomic data from SLIT2- and ROBO1-silenced HTR8/SVneo cells

Knockdown of SLIT2 and ROBO1 with siRNAs, total RNA isolation, cDNA synthesis, and qPCR were done as described above, except 30 nM siRNAs were used in the forward transfection instead of 40 nM.

Statistical analysis of RT-qPCR results

Differences in mRNA expression levels among the phenotypes were assessed by nonparametric Mann–Whitney U-test with SPSS Statistics 20.0 (IBM Corporation).

Supporting information

S1 Fig. Manhattan plot showing the results of genome-wide association study of gestational age.

Each dot represents the–log10(p) value of a single SNP in association analysis. Blue line denotes the level of suggestive significance (–log10(p) > 5). Loci are indicated for regions with suggestive significance.

https://doi.org/10.1371/journal.pgen.1008107.s001

(TIF)

S2 Fig.

Regional association plots for SLIT2 region in GWAS of SPTB (A) and gestational age (B). Each dot represents −log10(p) value of a single SNP in association analysis. Blue line denotes level of suggestive significance (−log10(p) > 5).

https://doi.org/10.1371/journal.pgen.1008107.s002

(TIF)

S3 Fig. SLIT2 mRNA expression correlation with Z-score.

SLIT2 mRNA expression in the basal plate of the placenta compared to birthweight-for-gestational age Z-scores of spontaneously delivered infants (SPTB and STB samples together). SLIT2 mRNA levels correlated with Z-scores (p = 0.023).

https://doi.org/10.1371/journal.pgen.1008107.s003

(TIF)

S4 Fig. Gene expression changes after SLIT2 and ROBO1 silencing verified by qRT-PCR.

(A) Four genes/gene families affected by ROBO1 silencing according to RNA sequencing were verified by qRT-PCR. (B) Expression changes of inflammation-associated cytokines IL-6 and TNF-A after SLIT2 and ROBO1 silencing were also measured with qRT-PCR. SLIT2 and ROBO1 were post-transcriptionally silenced in HTR8/SVneo cell line by siRNA. mRNA levels of selected genes compared between SLIT2- or ROBO1-silenced cells and mRNA levels of untreated control cells. All mRNA levels normalized against mRNA levels of housekeeping gene CYC1. Columns represent median and SD values of the sample groups.

https://doi.org/10.1371/journal.pgen.1008107.s004

(TIF)

S1 Table. Clinical characteristics of discovery GWAS study population.

https://doi.org/10.1371/journal.pgen.1008107.s005

(DOCX)

S2 Table. Suggestive association signals in GWAS of gestational age.

https://doi.org/10.1371/journal.pgen.1008107.s006

(DOCX)

S3 Table. Suggestive association signals in GWAS of very preterm birth.

https://doi.org/10.1371/journal.pgen.1008107.s007

(DOCX)

S4 Table. Suggestive association signals in GWAS of moderate-to-late preterm birth.

https://doi.org/10.1371/journal.pgen.1008107.s008

(DOCX)

S5 Table. Suggestive association signals in GWAS of infants born after SPTB with PPROM.

https://doi.org/10.1371/journal.pgen.1008107.s009

(DOCX)

S6 Table. Suggestive association signals in GWAS of infants born after onset of preterm labor with intact fetal membranes.

https://doi.org/10.1371/journal.pgen.1008107.s010

(DOCX)

S7 Table. Pathways associated with spontaneous preterm birth.

https://doi.org/10.1371/journal.pgen.1008107.s011

(DOCX)

S8 Table. Gene ontology terms associated with SPTB.

https://doi.org/10.1371/journal.pgen.1008107.s012

(DOCX)

S9 Table. Significant gene ontology terms that include SLIT2 and ROBO1 genes.

https://doi.org/10.1371/journal.pgen.1008107.s013

(DOCX)

S10 Table. Upregulated genes after SLIT2 silencing in HTR8/SVneo cells.

SLIT2 silenced in HTR8/SVneo commercial cell line by siRNA. Transcriptome of these cells compared with transcriptome of cells treated with negative siRNA. Differentially expressed genes ranked based on FDR-adjusted p value and fold change (FC). Threshold of fold change was > 2.0, and threshold of FDR-adjusted p value was <0.05.

https://doi.org/10.1371/journal.pgen.1008107.s014

(DOCX)

S11 Table. Downregulated genes after SLIT2 silencing in HTR8/SVneo cells.

SLIT2 silenced in HTR8/SVneo commercial cell line by siRNA. Transcriptome of these cells compared with transcriptome of cells treated with negative siRNA. Differentially expressed genes ranked based on FDR-adjusted p value and fold change. Threshold of fold change was > 2.0, and threshold of FDR-adjusted p value was <0.05.

https://doi.org/10.1371/journal.pgen.1008107.s015

(DOCX)

S12 Table. Upregulated genes after ROBO1 silencing in HTR8/SVneo cells.

ROBO1 silenced in HTR8/SVneo cell line by siRNA. Transcriptome of these cells compared with transcriptome of cells treated with negative siRNA. Differentially expressed genes ranked based on FDR-adjusted p value and fold change. Threshold of fold change was > 2.0, and threshold of FDR-adjusted p value was <0.05.

https://doi.org/10.1371/journal.pgen.1008107.s016

(DOCX)

S13 Table. Downregulated genes after ROBO1 silencing in HTR8/SVneo cells.

ROBO1 silenced in HTR8/SVneo cell line by siRNA. Transcriptome of these cells compared with transcriptome of cells treated with negative siRNA. Differentially expressed genes ranked based on FDR-adjusted p value and fold change. Threshold of fold change was > 2.0, and threshold of FDR-adjusted p value was <0.01.

https://doi.org/10.1371/journal.pgen.1008107.s017

(DOCX)

S14 Table. Biological pathways affected by SLIT2 silencing in HTR8/SVneo cells.

https://doi.org/10.1371/journal.pgen.1008107.s018

(DOCX)

S15 Table. Biological pathways affected by ROBO1 silencing in HTR8/SVneo cells.

https://doi.org/10.1371/journal.pgen.1008107.s019

(DOCX)

Acknowledgments

We acknowledge the technical assistance of Maarit Haarala and sample collection by Dr. Mari Mahlman, Riitta Vikeväinen, Miia Lehto, and Sonja Eeli. We thank the Medical Bioinformatics Centre of Turku Centre for Biotechnology for the miRNA-seq data analysis. The Centre is supported by University of Turku, Åbo Akademi University, and Biocenter Finland. CSC–IT Center for Science, Finland, is acknowledged for computational resources.

Citation: Tiensuu H, Haapalainen AM, Karjalainen MK, Pasanen A, Huusko JM, Marttila R, et al. (2019) Risk of spontaneous preterm birth and fetal growth associates with fetal SLIT2. PLoS Genet 15(6): e1008107. https://doi.org/10.1371/journal.pgen.1008107

References

1. Blencowe H, Cousens S, Chou D, Oestergaard M, Say L, Moller AB, et al. Born too soon: the global epidemiology of 15 million preterm births. Reprod Health. 2013;10 Suppl 1: S2-4755-10-S1-S2. Epub 2013 Nov 15.

2. Menon R. Spontaneous preterm birth, a clinical dilemma: etiologic, pathophysiologic and genetic heterogeneities and racial disparity. Acta Obstet Gynecol Scand. 2008;87: 590–600. pmid:18568457

3. Goldenberg RL, Culhane JF, Iams JD, Romero R. Epidemiology and causes of preterm birth. Lancet. 2008;371: 75–84. pmid:18177778

4. Kemp MW. Preterm birth, intrauterine infection, and fetal inflammation. Front Immunol. 2014;5: 574. pmid:25520716

5. Agrawal V, Hirsch E. Intrauterine infection and preterm labor. Semin Fetal Neonatal Med. 2012;17: 12–19. pmid:21944863

6. Ott WJ. Intrauterine growth retardation and preterm delivery. Am J Obstet Gynecol. 1993;168: 1710–5; discussion 1715–7. pmid:8317512

7. Bezold KY, Karjalainen MK, Hallman M, Teramo K, Muglia LJ. The genomics of preterm birth: from animal models to human studies. Genome Med. 2013;5: 34. pmid:23673148

8. York TP, Strauss JF,3rd, Neale MC, Eaves LJ. Estimating fetal and maternal genetic contributions to premature birth from multiparous pregnancy histories of twins using MCMC and maximum-likelihood approaches. Twin Res Hum Genet. 2009;12: 333–342. pmid:19653833

9. York TP, Eaves LJ, Lichtenstein P, Neale MC, Svensson A, Latendresse S, et al. Fetal and maternal genes' influence on gestational age in a quantitative genetic analysis of 244,000 Swedish births. Am J Epidemiol. 2013;178: 543–550. pmid:23568591

10. Lunde A, Melve KK, Gjessing HK, Skjaerven R, Irgens LM. Genetic and environmental influences on birth weight, birth length, head circumference, and gestational age by use of population-based parent-offspring data. Am J Epidemiol. 2007;165: 734–741. pmid:17311798

11. Zhang G, Bacelis J, Lengyel C, Teramo K, Hallman M, Helgeland O, et al. Assessing the Causal Relationship of Maternal Height on Birth Size and Gestational Age at Birth: A Mendelian Randomization Analysis. PLoS Med. 2015;12: e1001865. pmid:26284790

12. Beaumont RN, Warrington NM, Cavadino A, Tyrrell J, Nodzenski M, Horikoshi M, et al. Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics. Hum Mol Genet. 2018;27: 742–756. pmid:29309628

13. Zhang G, Feenstra B, Bacelis J, Liu X, Muglia LM, Juodakis J, et al. Genetic Associations with Gestational Duration and Spontaneous Preterm Birth. N Engl J Med. 2017;377: 1156–1167. pmid:28877031

14. Hallman M, Haapalainen A, Huusko JM, Karjalainen MK, Zhang G, Muglia LJ, et al. Spontaneous premature birth as a target of genomic research. Pediatr Res. 2018.

15. Strauss JF,3rd, Romero R, Gomez-Lopez N, Haymond-Thornburg H, Modi BP, Teves ME, et al. Spontaneous preterm birth: advances toward the discovery of genetic predisposition. Am J Obstet Gynecol. 2018;218: 294–314.e2. pmid:29248470

16. Karjalainen MK, Ojaniemi M, Haapalainen AM, Mahlman M, Salminen A, Huusko JM, et al. CXCR3 Polymorphism and Expression Associate with Spontaneous Preterm Birth. J Immunol. 2015;195: 2187–2198. pmid:26209629

17. Hao K, Wang X, Niu T, Xu X, Li A, Chang W, et al. A candidate gene association study on preterm delivery: application of high-throughput genotyping technology and advanced statistical methods. Hum Mol Genet. 2004;13: 683–691. pmid:14976157

18. Uzun A, Dewan AT, Istrail S, Padbury JF. Pathway-based genetic analysis of preterm birth. Genomics. 2013;101: 163–170. pmid:23298525

19. Bacelis J, Juodakis J, Sengpiel V, Zhang G, Myhre R, Muglia LJ, et al. Literature-Informed Analysis of a Genome-Wide Association Study of Gestational Age in Norwegian Women and Children Suggests Involvement of Inflammatory Pathways. PLoS One. 2016;11: e0160335. pmid:27490719

20. Rappoport N, Toung J, Hadley D, Wong RJ, Fujioka K, Reuter J, et al. A genome-wide association study identifies only two ancestry specific variants associated with spontaneous preterm birth. Sci Rep. 2018;8: 226-017-18246-5.

21. Haataja R, Karjalainen MK, Luukkonen A, Teramo K, Puttonen H, Ojaniemi M, et al. Mapping a new spontaneous preterm birth susceptibility gene, IGF1R, using linkage, haplotype sharing, and association analysis. PLoS Genet. 2011;7: e1001293. pmid:21304894

22. Romero R, Dey SK, Fisher SJ. Preterm labor: one syndrome, many causes. Science. 2014;345: 760–765. pmid:25124429

23. Brose K, Bland KS, Wang KH, Arnott D, Henzel W, Goodman CS, et al. Slit proteins bind Robo receptors and have an evolutionarily conserved role in repulsive axon guidance. Cell. 1999;96: 795–806. pmid:10102268

24. Kidd T, Bland KS, Goodman CS. Slit is the midline repellent for the robo receptor in Drosophila. Cell. 1999;96: 785–794. pmid:10102267

25. Moore T, Dveksler GS. Pregnancy-specific glycoproteins: complex gene families regulating maternal-fetal interactions. Int J Dev Biol. 2014;58: 273–280. pmid:25023693

26. Camolotto S, Racca A, Rena V, Nores R, Patrito LC, Genti-Raimondi S, et al. Expression and transcriptional regulation of individual pregnancy-specific glycoprotein genes in differentiating trophoblast cells. Placenta. 2010;31: 312–319. pmid:20116096

27. Liao WX, Laurent LC, Agent S, Hodges J, Chen DB. Human placental expression of SLIT/ROBO signaling cues: effects of preeclampsia and hypoxia. Biol Reprod. 2012;86: 111. pmid:22262697

28. Fang M, Du H, Han B, Xia G, Shi X, Zhang F, et al. Hypoxia-inducible microRNA-218 inhibits trophoblast invasion by targeting LASP1: Implications for preeclampsia development. Int J Biochem Cell Biol. 2017;87: 95–103. pmid:28412444

29. Li P, Shi Y, Shuai H, Cai Y, Lu W, Wang G, et al. Alterted SLIT2/ROBO1 signalling is linked to impaired placentation of missed and threatened miscarriage in early pregnancy. Histopathology. 2017;71: 543–552. pmid:28485101

30. Li P, Peng H, Lu WH, Shuai HL, Zha QB, Yeung CK, et al. Role of Slit2/Robo1 in trophoblast invasion and vascular remodeling during ectopic tubal pregnancy. Placenta. 2015;36: 1087–1094. pmid:26282852

31. Li J, Oehlert J, Snyder M, Stevenson DK, Shaw GM. Fetal de novo mutations and preterm birth. PLoS Genet. 2017;13: e1006689. pmid:28388617

32. Hu H. Chemorepulsion of neuronal migration by Slit2 in the developing mammalian forebrain. Neuron. 1999;23: 703–711. pmid:10482237

33. Counsell SJ, Boardman JP. Differential brain growth in the infant born preterm: current knowledge and future developments from brain imaging. Semin Fetal Neonatal Med. 2005;10: 403–410. pmid:15993667

34. Larroque B, Ancel PY, Marret S, Marchand L, Andre M, Arnaud C, et al. Neurodevelopmental disabilities and special care of 5-year-old children born before 33 weeks of gestation (the EPIPAGE study): a longitudinal cohort study. Lancet. 2008;371: 813–820. pmid:18328928

35. Portmann-Lanz CB, Baumann MU, Mueller M, Wagner AM, Weiss S, Haller O, et al. Neurogenic characteristics of placental stem cells in preeclampsia. Am J Obstet Gynecol. 2010;203: 399.e1–399.e7.

36. Portmann-Lanz CB, Schoeberlein A, Portmann R, Mohr S, Rollini P, Sager R, et al. Turning placenta into brain: placental mesenchymal stem cells differentiate into neurons and oligodendrocytes. Am J Obstet Gynecol. 2010;202: 294.e1–294.e11.

37. Sitras V, Fenton C, Paulssen R, Vartun A, Acharya G. Differences in gene expression between first and third trimester human placenta: a microarray study. PLoS One. 2012;7: e33294. pmid:22442682

38. Hinck L. The versatile roles of "axon guidance" cues in tissue morphogenesis. Dev Cell. 2004;7: 783–793. pmid:15572123

39. Dickinson RE, Duncan WC. The SLIT-ROBO pathway: a regulator of cell function with implications for the reproductive system. Reproduction. 2010;139: 697–704. pmid:20100881

40. Gupta SK, Malhotra SS, Malik A, Verma S, Chaudhary P. Cell Signaling Pathways Involved During Invasion and Syncytialization of Trophoblast Cells. Am J Reprod Immunol. 2016;75: 361–371. pmid:26490782

41. Sharma S, Godbole G, Modi D. Decidual Control of Trophoblast Invasion. Am J Reprod Immunol. 2016;75: 341–350. pmid:26755153

42. PrabhuDas M, Bonney E, Caron K, Dey S, Erlebacher A, Fazleabas A, et al. Immune mechanisms at the maternal-fetal interface: perspectives and challenges. Nat Immunol. 2015;16: 328–334. pmid:25789673

43. Steinborn A, Gunes H, Halberstadt E. Signal for term parturition is of trophoblast and therefore of fetal origin. Prostaglandins. 1995;50: 237–252. pmid:8838236

44. Liao WX, Wing DA, Geng JG, Chen DB. Perspectives of SLIT/ROBO signaling in placental angiogenesis. Histol Histopathol. 2010;25: 1181–1190. pmid:20607660

45. Lim R, Liong S, Barker G, Lappas M. Slit2 is decreased after spontaneous labour in myometrium and regulates pro-labour mediators. J Reprod Immunol. 2014;106: 76–88. pmid:25130654

46. An B, Xia J, Chang T, Wang X, Miao J, Xu L, et al. Genome-wide association study identifies loci and candidate genes for internal organ weights in Simmental beef cattle. Physiol Genomics. 2018;50: 523–531. pmid:29676954

47. Bacelis J, Juodakis J, Adams Waldorf KM, Sengpiel V, Muglia LJ, Zhang G, et al. Uterine distention as a factor in birth timing: retrospective nationwide cohort study in Sweden. BMJ Open. 2018;8: e022929-2018-022929.

48. Lim R, Lappas M. Slit2 exerts anti-inflammatory actions in human placenta and is decreased with maternal obesity. Am J Reprod Immunol. 2015;73: 66–78. pmid:25329354

49. Tanno T, Fujiwara A, Sakaguchi K, Tanaka K, Takenaka S, Tsuyama S. Slit3 regulates cell motility through Rac/Cdc42 activation in lipopolysaccharide-stimulated macrophages. FEBS Lett. 2007;581: 1022–1026. pmid:17306799

50. Geutskens SB, Hordijk PL, van Hennik PB. The chemorepellent Slit3 promotes monocyte migration. J Immunol. 2010;185: 7691–7698. pmid:21078908

51. Zhao H, Anand AR, Ganju RK. Slit2-Robo4 pathway modulates lipopolysaccharide-induced endothelial inflammation and its expression is dysregulated during endotoxemia. J Immunol. 2014;192: 385–393. pmid:24272999

52. Velez DR, Fortunato SJ, Morgan N, Edwards TL, Lombardi SJ, Williams SM, et al. Patterns of cytokine profiles differ with pregnancy outcome and ethnicity. Hum Reprod. 2008;23: 1902–1909. pmid:18487217

53. Cappelletti M, Della Bella S, Ferrazzi E, Mavilio D, Divanovic S. Inflammation and preterm birth. J Leukoc Biol. 2016;99: 67–78. pmid:26538528

54. Nadeau-Vallee M, Obari D, Quiniou C, Lubell WD, Olson DM, Girard S, et al. A critical role of interleukin-1 in preterm labor. Cytokine Growth Factor Rev. 2016;28: 37–51. pmid:26684042

55. Genc MR, Gerber S, Nesin M, Witkin SS. Polymorphism in the interleukin-1 gene complex and spontaneous preterm delivery. Am J Obstet Gynecol. 2002;187: 157–163. pmid:12114904

56. Varner MW, Esplin MS. Current understanding of genetic factors in preterm birth. BJOG. 2005;112 Suppl 1: 28–31.

57. Mittal P, Romero R, Kusanovic JP, Edwin SS, Gotsch F, Mazaki-Tovi S, et al. CXCL6 (granulocyte chemotactic protein-2): a novel chemokine involved in the innate immune response of the amniotic cavity. Am J Reprod Immunol. 2008;60: 246–257. pmid:18782286

58. Martinez FF, Cervi L, Knubel CP, Panzetta-Dutari GM, Motran CC. The role of pregnancy-specific glycoprotein 1a (PSG1a) in regulating the innate and adaptive immune response. Am J Reprod Immunol. 2013;69: 383–394. pmid:23432664

59. Arnold LL, Doherty , Flor AW, Simon JA, Chou JY, Chan WY, et al. Pregnancy-specific glycoprotein gene expression in recurrent aborters: a potential correlation to interleukin-10 expression. Am J Reprod Immunol. 1999;41: 174–182. pmid:10326620

60. Silver RM, Heyborne KD, Leslie KK. Pregnancy specific beta 1 glycoprotein (SP-1) in maternal serum and amniotic fluid; pre-eclampsia, small for gestational age fetus and fetal distress. Placenta. 1993;14: 583–589. pmid:8290497

61. Tamsen L, Johansson SG, Axelsson O. Pregnancy-specific beta 1-glycoprotein (SP1) in serum from women with pregnancies complicated by intrauterine growth retardation. J Perinat Med. 1983;11: 19–25. pmid:6601703

62. Tamsen L, Axelsson O, Johansson SG. Serum levels of pregnancy-specific beta 1-glycoprotein (SP1) in women with pregnancies at risk. Gynecol Obstet Invest. 1983;16: 253–260. pmid:6416936

63. Pihl K, Larsen T, Laursen I, Krebs L, Christiansen M. First trimester maternal serum pregnancy-specific beta-1-glycoprotein (SP1) as a marker of adverse pregnancy outcome. Prenat Diagn. 2009;29: 1256–1261. pmid:19911417

64. Blois SM, Sulkowski G, Tirado-Gonzalez I, Warren J, Freitag N, Klapp BF, et al. Pregnancy-specific glycoprotein 1 (PSG1) activates TGF-beta and prevents dextran sodium sulfate (DSS)-induced colitis in mice. Mucosal Immunol. 2014;7: 348–358. pmid:23945545

65. Ballesteros A, Mentink-Kane MM, Warren J, Kaplan GG, Dveksler GS. Induction and activation of latent transforming growth factor-beta1 are carried out by two distinct domains of pregnancy-specific glycoprotein 1 (PSG1). J Biol Chem. 2015;290: 4422–4431. pmid:25548275

66. Bry K, Hallman M. Transforming growth factor-beta 2 prevents preterm delivery induced by interleukin-1 alpha and tumor necrosis factor-alpha in the rabbit. Am J Obstet Gynecol. 1993;168: 1318–1322. pmid:8475982

67. Bastek JA, Brown AG, Anton L, Srinivas SK, D'addio A, Elovitz MA. Biomarkers of inflammation and placental dysfunction are associated with subsequent preterm birth. J Matern Fetal Neonatal Med. 2011;24: 600–605. pmid:20822327

68. Bastek JA, Gomez LM, Elovitz MA. The role of inflammation and infection in preterm birth. Clin Perinatol. 2011;38: 385–406. pmid:21890015

69. Hertelendy F, Romero R, Molnar M, Todd H, Baldassare JJ. Cytokine-initiated signal transduction in human myometrial cells. Am J Reprod Immunol. 1993;30: 49–57. pmid:8311930

70. Molnar M, Romero R, Hertelendy F. Interleukin-1 and tumor necrosis factor stimulate arachidonic acid release and phospholipid metabolism in human myometrial cells. Am J Obstet Gynecol. 1993;169: 825–829. pmid:8238136

71. Young A, Thomson AJ, Ledingham M, Jordan F, Greer IA, Norman JE. Immunolocalization of proinflammatory cytokines in myometrium, cervix, and fetal membranes during human parturition at term. Biol Reprod. 2002;66: 445–449. pmid:11804961

72. Peltier MR. Immunology of term and preterm labor. Reprod Biol Endocrinol. 2003;1: 122-7827-1-122.

73. Sonnega A, Faul JD, Ofstedal MB, Langa KM, Phillips JW, Weir DR. Cohort Profile: the Health and Retirement Study (HRS). Int J Epidemiol. 2014;43: 576–585. pmid:24671021

74. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4: 7-015-0047-8. eCollection 2015.

75. Delaneau O, Howie B, Cox AJ, Zagury JF, Marchini J. Haplotype estimation using sequencing reads. Am J Hum Genet. 2013;93: 687–696. pmid:24094745

76. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009;5: e1000529. pmid:19543373

77. Marchini J, Howie B, Myers S, McVean G, Donnelly P. A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007;39: 906–913. pmid:17572673

78. Nam D, Kim J, Kim SY, Kim S. GSA-SNP: a general approach for gene set analysis of polymorphisms. Nucleic Acids Res. 2010;38: W749–54. pmid:20501604

79. Pruim RJ, Welch RP, Sanna S, Teslovich , Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26: 2336–2337. pmid:20634204

80. GTEx Consortium, Laboratory, Data Analysis &Coordinating Center (LDACC)-Analysis Working Group, Statistical Methods groups-Analysis Working Group, Enhancing GTEx (eGTEx) groups, NIH Common Fund, NIH/NCI, et al. Genetic effects on gene expression across human tissues.ident. 2017;550: 204–213.

81. MacArthur J, Bowler E, Cerezo M, Gil L, Hall P, Hastings E, et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res. 2017;45: D896–D901. pmid:27899670

82. Ward LD, Kellis M. HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016;44: D877–81. pmid:26657631

83. Cleal JK, Day P, Hanson MA, Lewis RM. Measurement of housekeeping genes in human placenta. Placenta. 2009;30: 1002–1003. pmid:19819546

84. Cleal JK, Day PL, Hanson MA, Lewis RM. Sex differences in the mRNA levels of housekeeping genes in human placenta. Placenta. 2010;31: 556–557. pmid:20398937

85. Drewlo S, Levytska K, Kingdom J. Revisiting the housekeeping genes of human placental development and insufficiency syndromes. Placenta. 2012;33: 952–954. pmid:23026742

86. Haapalainen AM, Karjalainen MK, Daddali R, Ohlmeier S, Anttonen J, Maatta TA, et al. Expression of CPPED1 in human trophoblasts is associated with timing of term birth. J Cell Mol Med. 2018;22: 968–981. pmid:29193784

87. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30: 207–210. pmid:11752295