- Research article
- Open Access
Sequence variants in oxytocin pathway genes and preterm birth: a candidate gene association study
BMC Medical Genetics volume 14, Article number: 77 (2013)
Preterm birth (PTB) is a complex disorder associated with significant neonatal mortality and morbidity and long-term adverse health consequences. Multiple lines of evidence suggest that genetic factors play an important role in its etiology. This study was designed to identify genetic variation associated with PTB in oxytocin pathway genes whose role in parturition is well known.
To identify common genetic variants predisposing to PTB, we genotyped 16 single nucleotide polymorphisms (SNPs) in the oxytocin (OXT), oxytocin receptor (OXTR), and leucyl/cystinyl aminopeptidase (LNPEP) genes in 651 case infants from the U.S. and one or both of their parents. In addition, we examined the role of rare genetic variation in susceptibility to PTB by conducting direct sequence analysis of OXTR in 1394 cases and 1112 controls from the U.S., Argentina, Denmark, and Finland. This study was further extended to maternal triads (maternal grandparents-mother of a case infant, N=309). We also performed in vitro analysis of selected rare OXTR missense variants to evaluate their functional importance.
Maternal genetic effect analysis of the SNP genotype data revealed four SNPs in LNPEP that show significant association with prematurity. In our case–control sequence analysis, we detected fourteen coding variants in exon 3 of OXTR, all but four of which were found in cases only. Of the fourteen variants, three were previously unreported novel rare variants. When the sequence data from the maternal triads were analyzed using the transmission disequilibrium test, two common missense SNPs (rs4686302 and rs237902) in OXTR showed suggestive association for three gestational age subgroups. In vitro functional assays showed a significant difference in ligand binding between wild-type and two mutant receptors.
Our study suggests an association between maternal common polymorphisms in LNPEP and susceptibility to PTB. Maternal OXTR missense SNPs rs4686302 and rs237902 may have gestational age-dependent effects on prematurity. Most of the OXTR rare variants identified do not appear to significantly contribute to the risk of PTB, but those shown to affect receptor function in our in vitro study warrant further investigation. Future studies with larger sample sizes are needed to confirm the findings of this study.
Preterm birth (PTB) is a major public health problem that accounts for approximately 10% of all births worldwide [1, 2]. In 2011, PTBs, occurring before 37 completed weeks of gestation, comprised 11.73% of live births in the United States . Prematurity is a leading cause of perinatal mortality and in survivors is a significant contributor to short- and long-term morbidity . The etiology of PTB is not completely understood, but multiple lines of evidence suggest that genetic factors play an important role in PTB, including a high rate of recurrence in individuals with a history of previous preterm delivery (PTD), a tendency to occur within families, and racial disparity [1, 5–7].
The oxytocin (OXT)-oxytocin receptor (OXTR) system provides promising candidate genes for studies of genetic contributions to prematurity. Its components influence a wide range of physiological, behavioral, and emotional processes in humans, but have been most extensively studied for their role in reproduction, particularly labor and parturition [8, 9]. OXT acts as an inducer of uterine contraction and the myometrium becomes increasingly sensitive to the action of OXT towards term . This increase in uterine sensitivity to OXT occurs concomitantly with an upregulation of OXTR mRNA and a dramatic increase in myometrial OXTR number, which peaks during early labor [8, 9]. The OXTR is a G protein-coupled receptor, which, upon stimulation, activates various intracellular signaling pathways such as the phosphoinositide cascade, eventually leading to uterine contraction [8, 10, 11]. The availability of OXT to the receptor is dependent on serum levels of leucyl/cystinyl aminopeptidase (LNPEP), an enzyme that hydrolyzes and inactivates the hormone [12, 13]. LNPEP, also known as PLAP, increases in maternal serum during pregnancy and plays an important role in fetal development as well as pregnancy maintenance via regulating OXT levels and activity [12, 13]. The efficacy of Atosiban, an OXT antagonist used to stop premature uterine contractions and delay PTD , provides further support for the importance of this hormonal system in prematurity.
Previous candidate gene association studies have provided evidence for genetic variation predisposing to PTB with four maternal and one fetal variants showing relatively consistent evidence of association . In the present study, we hypothesize that allelic variations in OXT, OXTR, and LNPEP contribute to the genetic predisposition to prematurity. We evaluated the association between maternal and fetal genotypes of single nucleotide polymorphisms (SNPs) in the three genes and PTD in U.S. patients. We also performed direct sequence analysis in four population groups to search for rare variants in the OXT-OXTR system (particularly in the OXTR gene) that might play a role in PTB. Selected rare OXTR variants identified in preterm cases were then further investigated to determine their functional significance.
The study population consisted of case and control patients from 4 countries, including the U.S. (4 sites; the University of Iowa Hospitals and Clinics in Iowa City, IA, Magee-Womens Hospital in Pittsburgh, PA, University of Rochester Medical Center in Rochester, NY, and Wake Forest University in Wake Forest, NC), Argentina (two centers; Instituto de Maternidad y Ginecología Nuestra Señora de las Mercedes in Tucumán and Hospital Provincial de Rosario in Rosario), Finland (the University of Helsinki in Helsinki), and Denmark (the island of Funen and the Danish National Birth Cohort (DNBC) ). All families provided signed informed consent for study enrollment in accordance with the protocols approved by research ethics committees in the U.S. (the University of Iowa Institutional Review Board (IRB), University of Pittsburgh IRB, University of Rochester Research Subjects Review Board, Wake Forest University Health Sciences IRB), Argentina (the Research Ethics Committee of Centro de Educación Médica e Investigaciones Clínicas), Finland (the Ethics Committee of the University Central Hospital, Helsinki), and Denmark (the Scientific-Ethical Committee of the Southern Danish Region and the Biomedical Research Ethics Committee of the Capital City Region of Denmark). PTD was defined as delivery before 37 completed weeks of gestation. Gestational age (GA) was determined by obstetrical assignment using the first day of the last menstrual period as well as ultrasound examination, and was confirmed by assessment at birth. Demographic and clinical data were collected by chart review and/or clinician interview. Table 1 summarizes demographic characteristics of each population studied. U.S. cases were families with either a spontaneous or medically indicated preterm birth - of note, twin births were excluded at the data analysis step. Cases from Finland (mother-infant dyads, <36 weeks’ gestation) were families with spontaneous onset of preterm singleton birth from 2003 to 2009, selected using the following exclusion criteria: elective deliveries without spontaneous onset of labor and deliveries in which either maternal (e.g. systemic infection) or fetal (e.g. malformation) disease with known predisposition to premature birth was indicated. Finnish control families were defined as those with two or more children, all of whom were delivered at 37 weeks or later. The median gestational ages for index mother-infant dyads were 34 and 40 weeks for cases and controls, respectively. Danish cases and controls were identified on Funen where the University Hospital of Odense was a major participating site for case collection, and from the DNBC . The DNBC consisted of 96,946 total births, including 5,352 preterm births, occurring between 1997 and 2002. We selected 723 case (≤36 weeks gestation) and 920 control (40 weeks gestation) mothers for our study. Individuals were considered for inclusion in this study only if there was no evidence of obstetrical induction, placental abnormalities, pre-eclampsia, congenital malformations, and multiple births. Biological samples collected from mothers were stored in the Biobank at Statens Serum Institut in Copenhagen. Epidemiologic information was collected from existing medical records and from three maternal interviews, during the first trimester, at about 30 weeks gestation and 6 months after birth.
A total of 812 premature infants born in Iowa between 1999 and 2008 with GAs between 21 and 36 weeks (mean, 31 ± 3.7 weeks) and one or both of their parents were selected for genotyping. DNA was extracted from cord blood or discarded venous blood from the infants and venous blood, buccal swab, or saliva from the parents using standard protocols. Genotyping for SNP markers was performed using TaqMan assays (Applied Biosystems, Foster City, CA, U.S.) as described in detail elsewhere . A total of 16 tagging SNPs were selected to cover haplotype blocks of the oxytocin gene (OXT; 1 SNP), the oxytocin receptor gene (OXTR; 9 SNPs), and the leucyl/cystinyl aminopeptidase gene (LNPEP; 6 SNPs). The SNPs genotyped, along with their minor allele frequencies (MAFs), are listed in Table 2. Twin births were excluded, leaving 651 singleton infants for analysis. Of the 651 infants, 565 were successfully genotyped with GAs between 23 and 36 weeks (mean 31.1 ± 3.75).
A subgroup of the genotyped Iowa patients was selected for phase 1 sequencing based on GA, including 94 infants and 95 mothers of different infants with GA between 24–32 weeks. All three coding exons of OXT and all four exons of OXTR were sequenced in both directions, as described in prior studies . Phase 2 (replication) sequencing, where we focused only on exon 3 of OXTR, was performed on cases and controls from Argentina (infants and mothers; case, N=212 and control, N=125), Denmark (mothers participating in the DNBC; case, N=723 and control, N=920), Finland (mothers; case, N=105 and control, N=95), and the U.S. (infants and mothers; case, N=719 and control, N=229). Removal of twin births as well as unsuccessfully genotyped individuals left a very small number of Argentinian and U.S. controls, making a case–control association analysis particularly for the Argentinian population impossible.
To assess the effect of maternal OXTR genotype on the risk of PTD, we performed sequence analysis on the same region of OXTR in maternal triads (each including a mother of a premature infant with GA between 23 and 36 weeks (single gestation, born after spontaneous onset of labor) and her parents) as well as in a small group of patients consisting of 71 case and 64 control mothers. These case and control mothers were all Caucasians from Iowa with a singleton spontaneous delivery between 2005 and 2010. Further exclusion criteria for cases include preeclampsia and induction of labor. Of the 71 cases (mothers of preterm infants with GA between 23–28 weeks and birth weight (BW) between 413 and 1770 grams) and the 64 controls (mothers of term infants with GA between 39–42 weeks and BW between 2050 and 4756 grams), 62 cases and 60 controls with complete data were included in the analysis. The maternal grandmother-grandfather-mother triads were created from 309 pedigrees from different populations, including Argentina (2 centers, N=187), Denmark (the island of Funen, N=13), and the U.S. (4 sites, N=109). The families were enrolled in the study between 2005 and 2010 (Argentina and U.S.) and between 2008 and 2009 (Denmark). Of note, the small sample size of the Danish group (N=13) made statistical analysis uninformative. DNA was extracted from venous blood or saliva (Argentina), saliva (Denmark), and venous blood, buccal swab, or saliva (U.S.) and subsequently used for sequencing. All sequencing was performed by Functional Biosciences, Inc., Madison, WI, U.S.
Plasmids and cells
A full-length cDNA encoding human OXTR [GenBank: AY389507] was obtained from the Missouri S&T cDNA Resource Center (Rolla, MO, U.S.). Site-directed mutagenesis was performed to generate mutant constructs possessing missense mutations of potential functional significance identified in premature cases, including P108A, W203R, and F284. The mutations were introduced individually into the wild-type OXTR sequence using a PCR-based mutagenesis strategy. The wild-type and mutant full-length cDNAs were inserted into the EcoRI and XhoI sites of expression vector pcDNA3.1 (+) (Invitrogen, Carlsbad, CA, U.S.). The sequence of all constructs was confirmed by sequencing.
African green monkey kidney COS-7 cells were maintained in Dulbecco’s Modified Eagle’s Medium supplemented with 10% newborn calf serum, 50 μg/mL gentamicin, and 10 mM HEPES, pH 7.4. Cells were plated at a density of 4 x 105 cells/well in 6-well (35-mm) plates, and transiently transfected with plasmid DNA using FuGENE® HD Transfection Reagent (Roche Applied Science, Indianapolis, IN, U.S.) according to the manufacturer’s instructions. After overnight incubation, the cells were washed and used in subsequent experiments.
Inositol phosphate assays
The accumulation of inositol phosphates was measured using [3H]inositol-labeled cells incubated in the absence or presence of 100 nM oxytocin for 1 h as described elsewhere . Results were presented as mean fold difference between stimulated and unstimulated values (in cpm/106 cells) or mean ratio of the basal level in each transfection group to that in the wild-type group.
Ligand binding assays
Two days after transfection, cells were washed twice with warm buffer (20 mM HEPES, 0.15 M NaCl, 1 mg/ml bovine serum albumin, pH 7.4), and incubated with 50 nM [3H]oxytocin (PerkinElmer, Waltham, MA, U.S.) alone (total binding) or together with 20 μM unlabeled oxytocin (nonspecific binding). After 2.5 hours of incubation at room temperature, the cells were washed twice with cold buffer and dissolved in 0.5 ml of 0.5 N NaOH overnight. After collecting the lysate, each well was washed with another 0.5 ml of 0.5 N NaOH, and the wash was combined with the lysate, followed by neutralization with 0.25 ml of 2 N HCl prior to measuring radioactivity. Specific binding was determined by calculating the difference between total and nonspecific binding and was expressed as mean cpm/106 cells.
The common SNP genotyping data were tested for genotype deviation from the Hardy–Weinberg equilibrium (HWE) using PLINK  and Haploview . We first analyzed the data for all GAs together and then performed stratified analyses. For both maternal and fetal effects, five-week sliding windows of GA were utilized, starting with 21–25 weeks and ending with 32–36 weeks. For fetal effect analysis, three-stage stratification was additionally used, resulting in three GA groups: early (between 21 and 27 weeks), middle (between 28 and 30 weeks), and late (between 31 and 36 weeks). We performed a transmission disequilibrium test (TDT) using the Family Based Association Test (FBAT) software [21, 22], which analyzes overtransmission of alleles to an affected fetus . For maternal genetic effect analysis, we used a log-linear approach to study maternally-mediated effects . This approach uses a case-parents triad as the analysis unit and it tests maternally-mediated genetic effects based on the symmetry assumption of allele counts between the mothers and the fathers in the source population . We performed likelihood ratio tests of the maternal genetic effects. The expectation-maximization algorithm was applied to fully utilize families with missing parental genotypes . We also performed a haplotype analysis for all the six SNPs in the LNPEP gene using TRIMM . The test statistics was constructed based on the vector of genotype differences between the mother and the father. Under genetic mating symmetry assumption, the difference vector has an expected value of 0 at each locus under the null hypothesis. A permutation test was used to evaluate the significance of the test statistic.
The sequencing data from the maternal triads were analyzed in a similar way except that we defined as affected those (mothers) who delivered preterm infants and used GA and BW of the infants as phenotypes of the case mothers. We performed a TDT using the FBAT program [21, 22] for each common or rare variant found in the sequenced region. The data were analyzed for all GAs together as well as for GA subgroups using five-week sliding windows of GA and three-stage (early, middle, and late) stratification as described above.
We performed chi-square and Fisher’s exact tests using SAS version 9.2 statistical program (SAS, Cary, NC, U.S.) for a 2 by 2 or 2 by 3 contingency table to analyze the sequence data on rare and common OXTR variants found in cases and controls from Argentina, Denmark, Finland, and the U.S. as well as a small group of Caucasian case and control mothers from Iowa. In the Iowa case and control study, we restricted our analysis to only 3 common polymorphisms (rs4686302, rs237902, and rs61740241) because of the absence of rare variant minor alleles in both cases and controls.
The Danish sequencing data were further analyzed using linear regression to model genetic effects of the OXTR variants identified on GA. In the Iowa case–control study, we conducted logistic regression analysis to predict case/control status and linear regression analyses to predict GA (weeks) and BW (grams) on the aforementioned three OXTR polymorphisms.
Data for inositol phosphate and ligand binding assays were analyzed using one-way ANOVA, followed by Dunnett's test with InStat software (GraphPad, La Jolla, CA, U.S.).
Analysis of genetic variation in OXTRusing publicly available reference data
The OXTR sequencing results from our study were additionally analyzed utilizing publicly available data from two databases: 1) the 1000 Genomes Project (1000GP) [28, 29], which provides a catalogue of genetic variants present in the genomes of healthy individuals from different ethnic groups, and 2) the National Heart, Lung, and Blood Institute (NHLBI) Grand Opportunity (GO) Exome Sequencing Project (ESP) , which provides a catalogue of exome variants present in the genomes of approximately 2500 European and African American individuals from cohorts sampled for heart, lung, and blood diseases . No data on these for GA is reported. We investigated the two databases in search of coding variants in exon 3 of the OXTR gene in European Americans and African Americans, and then compared the variants retrieved with those identified in our sequencing study. The variants reported in the 1000GP (February 2012 release, 20120316) [28, 29] were annotated using the SeattleSeq Annotation Server .
Effects of common polymorphisms in the OXT pathway genes on risk of PTB
We selected 16 SNPs in the OXT, OXTR, and LNPEP genes based on haplotype block structures and examined them for their association with PTB in patients from the U.S. (Iowa). All SNPs tested had a MAF greater than 10% (Table 2). Excluding one SNP (rs237889) with HWE violation (p=8 x 10-5) and one SNP (rs4686302) with insufficient informative families (n<10) left a total of 14 SNPs for analysis.
Fetal genetic effects
Figure 1 shows the p-values obtained from fetal effect analysis, using the entire set of data as well as subsets classified by GA. There was no significant SNP found either when all GAs were examined together (Figure 1A) or when the GA-stratified analyses were performed (Figure 1A and 1B).
Maternal genetic effects
Figure 2 summarizes the results from the analysis of maternally mediated genetic effects. When the analysis was performed on all GA cases (Figure 2A), three SNPs in LNPEP were found to have a p-value of less than 0.05 (rs4869315, p=0.035; rs3849749, p=0.019; rs4869317, p=0.006). These values are uncorrected for the multiple comparisons done in this study where a Bonferroni corrected p-value for significance would be <0.002 accounting for 14 SNPs and both maternal and fetal comparisons. When the analysis was stratified by GA using five-week sliding windows (Figure 2B), the same variants in LNPEP remained significant - rs4869315 for the 30–34 (p=0.019) and 32–36 (p=0.009) week GA groups, rs3849749 for the 32–36 (p=0.016) week GA group, and rs4869317 for the 31–35 (p=0.024) and 32–36 (p=0.012) week GA groups. Another LNPEP SNP, rs13175726 (p=0.035) was additionally found to be significant for the 26–30 week GA group. Haplotype analysis of all 6 SNPs in LNPEP failed to yield any significant results.
Analysis of rare variants in the OXTRgene for the association with PTB
Rare sequence variations in OXTR: phase 1 and phase 2 case–control resequencing studies
Lack of association between common polymorphisms in OXT and OXTR led us to examine whether there exist(s) rare genetic variant(s) in these two genes with potential effects on the risk of prematurity. To this end, we first sequenced all exons of OXT and OXTR in a subgroup (n=189) of the genotyped (U.S.) patients. This pilot (hereinafter referred to as “phase 1”) sequencing analysis revealed only previously reported SNPs in the OXT gene. In OXTR, we found three rare variants (V172A, L206V, W203R, Table 3) that were not present in public databases at the time of the data analysis. Of note, the mutant allele at codon 206 was always detected together with that at codon 172, indicative of a rare double-mutant haplotype, while the V172A rare allele was identified both on its own and together with the L206V rare allele. We found that all the missense variants identified were located in one high-yield amplicon of OXTR, which corresponds to exon 3, and therefore attempted to replicate these phase 1 sequencing results in more individuals.
Phase 2 (replication) sequencing was performed on preterm cases and term controls from the U.S., Argentina, Finland, and Denmark (Table 1). The OXTR variants identified and their locations in the receptor are shown in Table 3 and Figure 3, respectively. The F284L mutation (rs201783860), which was detected in one case mother from the U.S., was also transmitted to the infant. In a well-characterized cohort of case (≤36 weeks gestation, N=723) and control (40 weeks gestation, N=920) mothers selected from the DNBC Study , we found three novel variants, including A11T, V211L, and A247A (Table 3). All but four (V45L, V211L, A247A, and G252A) rare coding variants identified in phase 1 and 2 studies were found in cases only.
To gain a more complete understanding of the spectrum of genetic variation in exon 3 of OXTR and the importance of the identified OXTR rare variants in predisposition to PTB, we examined recent data from the 1000 Genomes Project (1000GP) [28, 29] and the NHLBI GO Exome Sequencing Project (NHLBI_ESP)  (see Methods for more details). A list of coding variants in exon 3 of OXTR reported in the two databases is shown in Additional file 1. We compared this list with our sequencing results and found 5 variants (V45L, V172A, L206V, G252A, and V281M) overlapping with our data in individuals of European origin. It was particularly noted that two (V172A and L206V) of the 5 variants occur at much higher frequencies in African Americans (MAF: V172A, 5.52% and L206V, 1.84%) than in European Americans (MAF: V172A, 0.01% and L206V, 0.01%). It would be interesting to explore if these variants play a role in prematurity, given that about 10% of normal individuals whose sequences are catalogued in the databases are expected to be preterm and that PTB rates are higher in African Americans than in whites . In addition, we found 4 rare variants (A11T, P108A, W203R, and F284L) that are present exclusively in our cases with PTB. Importantly, three of them were predicted to have functional effects by PolyPhen  and SIFT  (Table 4), suggesting that they are potentially etiologic.
Association analysis of all rare coding variants identified by phase 1 and 2 sequencing showed no significant results – either when the variants were tested individually or when all missense variants were combined and analyzed together (Table 5).
We further extended our analysis on Danish mothers, in particular, and performed regression analyses to more precisely define the effects of OXTR genetic variation on GA. We modeled the effects of all coding variants (both common and rare) identified by sequencing on GA treated as either a continuous (in weeks) or categorical (preterm, ≤36 weeks; and term, 40 weeks) variable. However, no results were found to be significant in both analyses (Additional file 2).
Maternal OXTR sequence variation and susceptibility to PTB
To address in more depth the importance of maternal genetic effects on the risk of PTB, we extended the results obtained from the study of Danish mothers by performing sequence analysis of OXTR in maternal triads created from 309 pedigrees from the U.S., Argentina, and Denmark (Table 1) as well as in a small group of Caucasian patients from Iowa consisting of case (n=71) and control (n=64) mothers. Each maternal triad consisted of a mother of a preterm infant with GA between 23 and 36 weeks and her parents (grandmother-grandfather-mother), and we examined the transmission of OXTR variant alleles from maternal grandparents to the mothers of affected infants using the transmission disequilibrium test (TDT).
Among 8 coding variants identified by sequencing, two (Q143R and E242Q) were previously unreported novel variants (Table 3). Our TDT analysis did not reveal any significant association either when we examined all GAs together (Figure 4) or three subgroups (early, middle, and late) stratified by GA. (Of note, due to very low frequencies of rare variant minor alleles, complete results were available only for 3 SNPs (rs4686302, rs237902, and rs61740241) with more than 10 informative families.) However, when we performed a stratified analysis for the U.S. triads (N=109) using five-week sliding windows of GA, rs4686302 and rs237902 were found to be significant for the 24–28 and 28–32 week groups and the 27–31 and 28–32 week groups, respectively (Figure 4).
The three SNPs were further studied in the Iowa case–control study (71 case and 64 control mothers). We first developed logistic regression models to predict case/control status, and found that the model containing the two polymorphisms rs4686302 and rs237902 best explains case–control status (p=0.04, 0.05, respectively, R2=0.062) (Table 6). Next, linear regression analyses were performed to predict GA and BW. As demonstrated in Table 6, the models containing (rs237902 + rs61740241) and (rs4686302 + rs61740241) were best at predicting lower GA (p=0.03, 0.04, R2=0.055) and lower BW (p=0.008, 0.01, R2=0.074), respectively. The results indicate that 2 SNPs significantly explain 5.5% of the variation in GA and 7.3% of the variation in BW with rs61740241 in a slightly more supportive role to the other SNP. The SNPs that best predict case–control status are identical to the SNPs involved as the main predictors of GA in weeks and BW in grams.
Functional analysis of OXTRmissense mutations and polymorphisms
We performed in vitro functional studies on the aforementioned three OXTR missense mutations (P108A, W203R, and F284L) that were uniquely identified in our PTB cases and predicted as deleterious by PolyPhen  and SIFT  (Table 4). We first tested the ligand-binding properties of the mutant receptors, using transiently transfected COS-7 cells. To this end, the level of [3H]OXT binding was assessed in empty vector-, wild-type or mutant OXTR-expressing cells in the presence or absence of excess unlabeled OXT. The specific binding of [3H]OXT to intact transfected cells was defined as the difference between total and non-specific binding. We initially observed a considerable decrease in specific binding in cells transfected with the mutant plasmids, compared to those transfected with the wild-type plasmid. Because this could be a reflection of lower expression of the mutant receptor proteins, we attempted to increase their expression levels using twice more plasmid DNA (4 μg per well). However, the mean binding levels (in cpm/million cells) in each transfection group hardly varied. Based on this observation, we combined data from all experiments performed, regardless of the amount of transfected DNA. As shown in Figure 5, two of the three mutations (P108A and W203R) tested exhibited a statistically significant reduction in [3H]OXT binding. For the W203R mutation, there was nearly complete abolishment of specific binding. Our results suggest that these mutations may cause either reduced receptor expression or ligand affinity, leading to impaired ligand binding.
To gain insight into whether the missense mutations affect receptor signaling, OXT-stimulated inositol phosphate (InsP) production by wild-type and mutant OXTRs was evaluated. Figure 6A summarizes the results obtained. The accumulation of [3H]InsPs increased more than twofold in response to oxytocin in wild-type OXTR-expressing cells, while this increase was not observed in mutant plasmid- and empty vector-transfected cells. In addition, basal (unstimulated) levels of inositol phosphates were almost indistinguishable between wild-type OXTR and the three mutant receptors (Figure 6B), suggesting that the mutant receptors are unlikely to be constitutively active. Taken together, these data indicate that the examined mutations may have no effect on the coupling of the receptor to the second messenger signaling system.
Pregnancy and parturition involve an intricate biochemical and molecular interplay between mother and fetus . The study of PTB has been hampered by incomplete knowledge of these complex, time-varying processes, and an investigation of static makers like sequence variants may serve as one way to predict the risk of PTB.
In this study, we used a candidate gene approach to better understand genetic factors contributing to PTB. Previous studies have identified some polymorphisms suggestively associated with PTB, but with a focus mostly on immune-related and inflammatory genes such as interleukin 1 receptor antagonist (IL1RN), interferon gamma (IFNG), and coagulation factor II (F2) . The present study was undertaken to explore the contribution of genetic variations in another functionally relevant pathway, the oxytocin pathway to the etiology of PTB, given the important role of constituents of the pathway in labor and parturition.
Our analysis of fourteen SNPs in OXT, OXTR, and LNPEP revealed no evidence of fetal genetic effects on PTB. However, the analysis of maternal genetic effects performed all GA cases identified three SNPs in LNPEP that correlate with prematurity when not corrected for multiple comparisons. This finding is in line with the results of recent studies demonstrating that the maternal genome has a major impact on the risk of prematurity [35, 36] and can be considered hypothesis generating for further investigations in the role of LNPEP in preterm labor (PTL).
Ishii et al. reported that in LNPEP knockout mice, an increased sensitivity to OXT is observed and pregnancy duration is significantly shortened, providing functional evidence supporting our genetic data. It was noted that the same three LNPEP SNPs showing suggestive association in the analysis of all GAs remained significant in the GA-stratified analysis, further supporting the importance of these SNPs. The activity of LNPEP is known to gradually increase during late pregnancy and reach a very high level at 11 days before the onset of labor . Interestingly, in our five-week sliding window analysis for maternal effects, the three SNPs were found to be significant in late GA groups (30–34, 31–35, and 32–36 week groups) although there was an additional significant SNP (rs13175726) detected for the 26–30 week group. Pending replication with larger samples, we postulate that these polymorphisms, all of which are non-coding SNPs, may serve as regulatory elements for the LNPEP gene rather than having a direct functional effect.
We screened the OXTR gene for rare genetic variation by resequencing and identified both previously reported and unreported rare variants in exon 3, which is highly conserved across species (Additional file 3). The variants identified were not found to be significantly associated with PTB, but these results should be interpreted cautiously because of the small sample size and the paucity of sequencing data on controls, especially those from Argentina and the U.S.
Our analysis comparing the rare variants identified by sequencing with those catalogued in the public databases demonstrates that some of the rare variants may be of etiologic importance. The observation that three (P108A, W203R, and F284L) of the four rare missense variants not found in either the 1000GP [28, 29] or the NHLBI_ESP  data were predicted to affect the function of the receptor protein particularly merits further attention. These results suggest that rare OXTR coding variants appear to make some contribution to the genetic risk of prematurity, but larger studies will be needed to measure their combined as well as individual impact more accurately and determine the diagnostic or therapeutic relevance of the findings.
We studied the role of maternal OXTR genotype in PTB in more depth using maternal grandparent-mother triads. Transmission disequilibrium analysis revealed no significant transmission distortion of the OXTR coding variants examined when the analysis was performed on all GA cases or stratified into 3 GA groups. However, two missense SNPs (rs4686302 and rs237902) showed nominal associations in the U.S. triads for the sliding window GA subgroups primarily involving 24–32 weeks of gestation. These results suggest GA-specific effects of maternal genotypes of these OXTR SNPs on susceptibility to PTB.
We characterized the effects of three mutations of potential functional significance with in vitro assays. Our radioligand binding assay revealed that P108A and W203R mutant-transfected cells display significantly diminished levels of [3H]OXT specific binding, compared to wild-type-transfected cells. Although it should be further clarified if the mutations affect OXTR protein expression, ligand binding affinity, or both, our results indicate the importance of these residues in receptor function and support previous findings that the first and second extracellular loops of OXTR, where P108 and W203 residues, respectively, are located, are important for agonist binding and selectivity . Despite their important biological effects, it is yet unclear how the mutations may contribute to the disease as it is counterintuitive that attenuated OXT signaling leads to PTB. However, it is possible that the variants may act indirectly to stimulate uterine contractions. Given that the regulation of OXTR function is dependent on steroids such as estradiol, progesterone, and its metabolites [39–41], it is probable that the mutations may affect receptor-hormone interaction in a way that interferes with the timely onset of labor.
Another mutation, F284L, also exhibited a reduced level of specific binding, but the result did not quite reach statistical significance. An earlier mutagenesis study has shown that when two transmembrane residues Y209 and F284 are mutated to F and Y residues, respectively, arginine vasopressin (AVP), a partial agonist of the OXTR, becomes a full agonist [8, 42]. The finding indicates that these residues may be critical for modulation of the receptor response to AVP. Given the uterotonic action of AVP , an investigation into the sensitivity of the F284 as well as the other two mutants to other agonists like AVP could possibly provide an explanation about the pathogenic effects of those mutations.
Our study has certain limitations. First of all, the sample size of our association study may not be large enough to provide sufficient statistical power for the analysis of rare variants with small effects. Second, the functional assays were not carried out in a cellular context relevant to address mechanisms of PTB. The functional role of the identified OXTR mutations should be determined using physiologically more relevant cell types such as uterine myometrial cells. In addition, we only focused on the InsP pathway to examine the downstream effects of the OXTR mutations. We did not observe any significant difference in OXT-induced InsP production between wild type and mutant receptor-expressing cells, and it is possible that their effects are mediated through other signal transduction pathway(s). Therefore, it would be important to examine the impact of the mutations on different biologic pathways like the MAP kinase pathway .
This study provides preliminary evidence for an association between maternal common polymorphisms in the LNPEP gene and the risk of PTB. In addition, we found that maternal OXTR missense SNPs rs4686302 and rs237902 may have gestational age-dependent effects on prematurity. These observations support that sequence variants in maternal genes may have a primary role in susceptibility to PTB. Most of the OXTR rare variants identified do not appear to significantly contribute to the risk of PTB. However, the importance of the non-synonymous rare variants shown to cause altered ligand binding warrants further investigation. Future genetic studies with larger sample sizes are needed to confirm the findings of the association analyses performed in the present study. Additional functional studies may provide more detailed insight into the precise physiological effects and clinical significance of the genetic variants identified and ultimately, accelerate translation of such findings into genetic-based diagnostic and therapeutic applications.
Single nucleotide polymorphism database
Danish National Birth Cohort
Extracellular N-terminal domain
Family Based Association Test
Coagulation factor II (thrombin)
4-(−2- hydroxyethyl)-1-piperazineethanesulfonic acid
Interleukin 1 receptor antagonist
Institutional Review Board
Minor allele frequency
The National Heart, Lung, and Blood Institute Grand Opportunity Exome Sequencing Project
Placental alkaline phosphatase
Sorting intolerant from tolerant
Single nucleotide polymorphism
Transmission disequilibrium test
The 1000 Genomes Project
Dizon-Townson DS: Preterm labour and delivery: a genetic predisposition. Paediatr Perinat Epidemiol. 2001, 15: 57-62. 10.1046/j.1365-3016.2001.00008.x.
Beck S, Wojdyla D, Say L, Pilar Betran A, Merialdi M, Harris Requejo J, Rubens C, Menon R, Van Look P: The worldwide incidence of preterm birth: a systematic review of maternal mortality and morbidity. Bull World Health Organ. 2010, 88: 31-38. 10.2471/BLT.08.062554.
Martin JA, Hamilton BE, Ventura SJ, Osterman MJ, Mathews TJ: Births: final Data for 2011. Natl Vital Stat Rep. 2013, 62 (1): 1-90.
Gotsch F, Romero R, Erez O, Vaisbuch E, Kusanovic JP, Mazaki-Tovi S, Kim SK, Hassan S, Yeo L: The preterm parturition syndrome and its implications for understanding the biology, risk assessment, diagnosis, treatment and prevention of preterm birth. J Matern Fetal Neonatal Med. 2009, 22 (Suppl 2): 5-23.
Crider KS, Whitehead N, Buus RM: Genetic variation associated with preterm birth: a huge review. Genet Med. 2005, 7 (9): 593-604. 10.1097/01.gim.0000187223.69947.db.
Adams KM, Eschenbach DA: The genetic contribution towards preterm delivery. Semin Fetal Neonatal Med. 2004, 9 (6): 445-452. 10.1016/j.siny.2004.04.001.
Muglia LJ, Katz M: The enigma of spontaneous preterm birth. N Eng J Med. 2010, 362 (6): 529-535. 10.1056/NEJMra0904308.
Gimpl G, Fahrenholz F: The oxytocin receptor system: structure, function, and regulation. Physiol Rev. 2001, 81 (2): 629-683.
Blanks AM, Shmygol A, Thornton S: Regulation of oxytocin receptors and oxytocin receptor signaling. Semin Reprod Med. 2007, 25 (1): 52-59. 10.1055/s-2006-956775.
Trzidou V: Biochemical and endocrinological preparation for parturition. Best Pract Res Clin Obstet Gynaecol. 2007, 21: 729-756. 10.1016/j.bpobgyn.2007.05.001.
Vrachnis N, Malamas FM, Sifakis S, Deligeoroglou E, Iliodromiti Z: The oxytocin-oxytocin receptor system and its antagonists as tocolytic agents. International journal of endocrinology. 2011, 2011: 350546-
Lampelo S, Vanhaperttula T: Fractionation and characterization of cystine aminopeptidase (oxytocinase) and arylamidase of human serum during pregnancy. J Reprod Fertil. 1980, 58 (1): 225-235. 10.1530/jrf.0.0580225.
Nomura S, Ito T, Yamamoto E, Sumigama S, Iwase A, Okada M, Shibata K, Ando H, Ino K, Kikkawa F, et al: Gene regulation and physiological function of placental leucine arninopeptidase/oxytocinase during pregnancy. Biochimica Et Biophysica Acta-Proteins and Proteomics. 2005, 1751 (1): 19-25. 10.1016/j.bbapap.2005.04.006.
Tsatsaris V, Carbonne B, Cabrol D: Atosiban for preterm labour. Drugs. 2004, 64 (4): 375-382. 10.2165/00003495-200464040-00003.
Dolan SM, Hollegaard MV, Merialdi M, Betran AP, Allen T, Abelow C, Nace J, Lin BK, Khoury MJ, Ioannidis JP, et al: Synopsis of preterm birth genetic association studies: the preterm birth genetics knowledge base (PTBGene). Public Health Genomics. 2010, 13 (7–8): 514-523.
Olsen J, Melbye M, Olsen SF, Sorensen TI, Aaby P, Andersen AM, Taxbol D, Hansen KD, Juhl M, Schow TB, et al: The Danish national birth cohort–its background, structure and aim. Scand J Public Health. 2001, 29 (4): 300-307.
Ehn NL, Cooper ME, Orr K, Shi M, Johnson MK, Caprau D, Dagle J, Steffen K, Johnson K, Marazita ML, et al: Evaluation of fetal and maternal genetic variation in the progesterone receptor gene for contributions to preterm birth. Pediatr Res. 2007, 62 (5): 630-635. 10.1203/PDR.0b013e3181567bfc.
Ascoli M, Pignataro OP, Segaloff DL: The inositol phosphate/diacylglycerol pathway in MA-10 leydig tumor cells. Activation by arginine vasopressin and lack of effect of epidermal growth factor and human choriogonadotropin. J Biol Chem. 1989, 264 (12): 6674-6681.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, Bender D, Maller J, Sklar P, de Bakker P, Daly M, et al: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Rabinowitz D, Laird N: A unified approach to adjusting association tests for population admixture with arbitrary pedigree structure and arbitrary missing marker information. Hum Hered. 2000, 50 (4): 211-223. 10.1159/000022918.
Laird NM, Horvath S, Xu X: Implementing a unified approach to family-based tests of association. Genet Epidemiol. 2000, 19: S36-S42. 10.1002/1098-2272(2000)19:1+<::AID-GEPI6>3.0.CO;2-M.
Spielman RS, McGinnis RE, Ewens WJ: Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet. 1993, 52 (3): 506-516.
Weinberg CR, Wilcox AJ: Distinguishing the effects of maternal and offspring genes through studies of 'case-parent triads' " and "a new method for estimating the risk ratio in studies using case-parental control design. Am J Epidemiol. 1999, 150 (4): 428-429.
Schaid DJ, Sommer SS: Genotype relative risks: methods for design and analysis of candidate-gene association studies. Am J Hum Genet. 1993, 53 (5): 1114-1126.
Weinberg CR: Allowing for missing parents in genetic studies of case-parent triads. Am J Hum Genet. 1999, 64 (4): 1186-1193. 10.1086/302337.
Shi M, Umbach DM, Vermeulen SH, Weinberg CR: Making the most of case-mother/control-mother studies. Am J Epidemiol. 2008, 168 (5): 541-547. 10.1093/aje/kwn149.
The 1000 Genomes Project Consortium: A map of human genome variation from population-scale sequencing. Nature. 2010, 467 (7319): 1061-1073. 10.1038/nature09534.
The 1000 genomes project. http://www.1000genomes.org/,
Exome variant server, NHLBI exome sequencing project (ESP). Seattle, WA, (URL: http://evs.gs.washington.edu/EVS/) [date (February, 2013) accessed]
The SeattleSeq annotation server. [http://gvs.gs.washington.edu/SeattleSeqAnnotation/]
Dolan SM: Genetic and environmental contributions to racial disparities in preterm birth. Mt Sinai J Med. 2010, 77 (2): 160-165. 10.1002/msj.20169.
Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic Acids Research. 2002, 30 (17): 3894-3900. 10.1093/nar/gkf493.
Ng PC, Henikoff S: Accounting for human polymorphisms predicted to affect protein function. Genome Research. 2002, 12 (3): 436-446. 10.1101/gr.212802.
Boyd HA, Poulsen G, Wohlfahrt J, Murray JC, Feenstra B, Melbye M: Maternal contributions to preterm delivery. American Journal of Epidemiology. 2009, 170 (11): 1358-1364. 10.1093/aje/kwp324.
Svensson AC, Sandin S, Cnattingius S, Reilly M, Pawitan Y, Hultman CM, Lichtenstein P: Maternal effects for preterm birth: a genetic epidemiologic study of 630,000 families. American Journal of Epidemiology. 2009, 170 (11): 1365-1372. 10.1093/aje/kwp328.
Ishii M, Naruse K, Hattori A, Tsujimoto M, Ishiura S, Numaguchi Y, Murohara T, Kobayashi H, Mizutani S: Oxytocin hypersensitivity in pregnant P-LAP deficient mice. Life Sciences. 2009, 84 (19–20): 668-672.
Postina R, Kojro E, Fahrenholz F: Separate agonist and peptide antagonist binding sites of the oxytocin receptor defined by their transfer into the V-2 vasopressin receptor. Journal of Biological Chemistry. 1996, 271 (49): 31593-31601. 10.1074/jbc.271.49.31593.
Grazzini E, Guillon G, Mouillac B, Zingg HH: Inhibition of oxytocin receptor function by direct binding of progesterone. Nature. 1998, 392 (6675): 509-512. 10.1038/33176.
Gellersen B, Fernandes MS, Brosens JJ: Non-genomic progesterone actions in female reproduction. Hum Reprod Update. 2009, 15 (1): 119-138.
Dunlap KA, Stormshak F: Nongenomic inhibition of oxytocin binding by progesterone in the ovine uterus. Biol Reprod. 2004, 70 (1): 65-69.
Chini B, Mouillac B, Balestre MN, TrumppKallmeyer S, Hoflack J, Hibert M, Andriolo M, Pupier S, Jard S, Barberis C: Two aromatic residues regulate the response of the human oxytocin receptor to the partial agonist arginine vasopressin. Febs Letters. 1996, 397 (2–3): 201-206.
Chan WY, Wo NC, Manning M: The role of oxytocin receptors and vasopressin V-1a receptors in uterine contractions in rats: implications for tocolytic therapy with oxytocin antagonists. American Journal of Obstetrics and Gynecology. 1996, 175 (5): 1331-1335. 10.1016/S0002-9378(96)70050-9.
Tahara A, Tsukada J, Tomura Y, Wada K, Kusayama T, Ishii N, Yatsu T, Uchida W, Tanaka A: Pharmacologic characterization of the oxytocin receptor in human uterine smooth muscle cells. British Journal of Pharmacology. 2000, 129 (1): 131-139. 10.1038/sj.bjp.0702996.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/14/77/prepub
We thank, in particular, the families that participated in this study and we are ever grateful for the clinical care provided by our colleagues that allowed those families to be enrolled. The research nurses Susan Berends, Gretchen Cress, Karen Johnson, Laura Knosp, Nancy Krutzfield, Silvina Argañaraz, Mirta Gladys Leguizamón, Mercedes López, Marta Padilla, and Azucena Singh were essential in contacting, enrolling and sampling families. Kristin Orr, Cara Zimmerman, Keegan Kelsey, Marla K. Johnson, and Tamara Busch were integral in sample organization or testing. We acknowledge Elizabeth Leslie for her help in obtaining the 1000GP data as well as her input on the analysis of the reference sequence data. Susie McConnell, Dan Benton, Erin Brothers-Smith, and Nancy Davin were essential in administration. The authors would like to thank the NHLBI GO Exome Sequencing Project and its ongoing studies which produced and provided exome variant calls for comparison: the Lung GO Sequencing Project (HL-102923), the WHI Sequencing Project (HL-102924), the Broad GO Sequencing Project (HL-102925), the Seattle GO Sequencing Project (HL-102926) and the Heart GO Sequencing Project (HL-103010). This work was supported by the National Institutes of Health [R01 HD052953 and HD57192 to J.C.M., Z01 ES040007 and Z01 ES045002 to M.S.]; and the March of Dimes Foundation [FY2006-575 and FY2008-260 to J.C.M., 21-FY07-597 to L.J.M.].
The authors declare that they have no competing interests.
JK performed SNP genotyping, sequencing, comparative sequence analysis, and functional assays (including data analysis), and drafted and revised the manuscript and prepared figures and tables. KJS carried out SNP genotyping and sequencing, in part, and drafted part of the manuscript. MEC conducted statistical analyses of SNP genotype (fetal genetic effects) and sequencing (Danish cohort, maternal triads, and Iowa cases and controls) data, and assisted in study design. MA was involved in the analysis of functional assay data and assisted in study design. AMM participated in the management and preparation of genotyping and sequencing samples and data for analysis. KKR participated in statistical analysis of Danish sequencing data. ELM performed maternal triad sequencing. LR was involved in the management and sequencing of Danish cohort samples. KLS helped with sample organization. VC, EG, and CS coordinated patient recruitment, phenotyping, and sample collection, and assisted in study design. MS conducted statistical analysis of SNP genotype (maternal genetic effects) data. MH, JP, KAT, and LJM provided Finnish DNA samples. BF, FG, HAB, and MM provided Danish DNA samples and patient information, and were involved in study design and discussion of statistical analysis of Danish sequencing data. MLM participated in study design and discussion of statistical analysis of genotype and sequencing data. JMD assisted in study design and reviewed the draft of the manuscript. JCM is the principal investigator and obtained funding, established key collaborations, provided study design and supervision, and oversaw manuscript preparation. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: List of rare missense variants in exon 3 of OXTR catalogued in public databases (1000GP, NHLBI_ESP). (PDF 78 KB)
Additional file 2: Table S2: Results of regression analyses of OXTR coding variants predicting gestational age (GA) in Danish case and control mothers. (PDF 69 KB)
Additional file 3: Table S3: Alignment of the amino acid sequence of OXTR from different species at mutant residues of all coding variants identified in the present study. (PDF 63 KB)
Authors’ original submitted files for images
About this article
Cite this article
Kim, J., Stirling, K.J., Cooper, M.E. et al. Sequence variants in oxytocin pathway genes and preterm birth: a candidate gene association study. BMC Med Genet 14, 77 (2013). https://doi.org/10.1186/1471-2350-14-77
- Preterm birth
- Genetic association analysis
- Oxytocin pathway
- Single nucleotide polymorphism
- Rare variant