What this is
- This research investigates the relationship between of the oxytocin receptor gene (OXTR) and oxytocin needs during childbirth.
- It explores how variations in OXTR methylation correlate with the amount of oxytocin administered and the risk of .
- Findings suggest that blood-derived could serve as a biomarker for predicting risk.
Essence
- Higher OXTR in blood correlates with increased oxytocin requirements during labor and greater risk. This suggests that assessing methylation levels could help identify individuals at risk for complications.
Key takeaways
- Blood and uterine myometrium levels at several OXTR loci are significantly correlated. This finding supports the use of blood samples as a surrogate for assessing uterine methylation.
- Females with higher OXTR required more oxytocin during labor. Specifically, each 5-percentage point increase in methylation was associated with a need for 2.48 more units of oxytocin.
- Higher OXTR was linked to a greater relative risk for when oxytocin was used during labor. A 5-percentage point increase in methylation was associated with a 45% higher risk for .
Caveats
- The sample size was limited, affecting the robustness of the findings. Recruitment was impacted by the COVID-19 pandemic, leading to fewer participants than planned.
- The study's observational nature limits the ability to establish causation between OXTR methylation and oxytocin needs or .
Definitions
- DNA methylation: A biochemical process involving the addition of a methyl group to DNA, affecting gene expression without altering the DNA sequence.
- Postpartum hemorrhage: Excessive bleeding following childbirth, defined as blood loss of 1000 mL or more within 24 hours after delivery.
Simplified
Introduction
Successful physiologic parturition resulting in vaginal birth of an infant, placenta and survival of the mother depends upon effective uterine contractility1. The myometrium of the uterus consists of layers of smooth muscle that utilize several proteins to generate coordinated and rhythmic contractions that may last from hours to days throughout the process of parturition. These proteins include oxytocin receptor (OXTR), gap junction (Connexin-43) and several isoforms of prostaglandin receptor2–4. While the action of prostaglandins primarily prompt labor onset, the latter stages of labor (active and expulsive/second stage and placental delivery) are more dependent upon oxytocin, as a potent uterotonic hormone5–8.
Oxytocin produced by the placenta can stimulate some uterine contractility9,10. However, in advanced stages of labor, a positive feedback mechanism, resulting from vaginal distension and the descending fetal head/body, triggers maternal hypothalamic oxytocin pulsation, and facilitates the final stage of labor11,12. Following birth of the newborn, maternal oxytocin continues to stimulate uterine contraction, causing placental detachment from the uterine wall and expelling it from the body13–15. Contractions during the third stage of labor also minimize bleeding at the placental site as the blood vessels that were supplying the placenta close off. Clotting factors aid in uterine tone and provide stasis to the open vessels16. While some bleeding is normative or physiologic after birth17,18, if uterine contractions are delayed, infrequent or absent (uterine atony), rapid and excessive bleeding after birth can occur.
Oxytocin acts upon OXTR expressed in smooth muscle cells within the myometrium4. OXTR is also present in many regions of the central nervous system19 and other peripheral organs/tissues20. The regulation of OXTR in the myometrium has been investigated for its role in labor processes4 with the quantity and responsiveness of OXTR playing an important role in uterine contraction during labor21. Furthermore, variability in OXTR function could also influence the postpartum outcomes of uterine atony and subsequent postpartum hemorrhage22.
Postpartum hemorrhage is currently defined as blood loss accompanied by signs of hypovolemia or an accumulation of 1000 mL or greater blood loss in the first 24 h postpartum23. Earlier definitions classified postpartum hemorrhage as 500 mL after vaginal birth and 1000 mL after Cesarean birth. Hemorrhage is a leading cause of maternal death worldwide24, and though it accounts for only 11.4% of maternal mortality in the United States25, recent data indicate that postpartum hemorrhage itself, and/or the need for life-saving interventions (blood transfusions, surgical interventions) to control or treat postpartum hemorrhage, are on the rise26–29. Sequela of postpartum hemorrhage beyond the immediate hospitalization period include prolonged recovery from childbirth, delays in lactation onset or difficulty producing milk30, anemia31,32, fatigue, and some studies have linked postpartum hemorrhage with postpartum depression and post-traumatic stress33–35. Other labor/maternal characteristics associated with postpartum hemorrhage include length of labor, genital tract lacerations, fetal macrosomia, hypertension/preeclampsia, or a prior history of postpartum hemorrhage28,36. Many factors are closely associated with each other, which compounds the risk to the individual and complicates research methods. Regardless of known risk factors, up to 40% of postpartum hemorrhage occurs in individuals without risk factors, making postpartum hemorrhage often difficult to anticipate27,37–39.
Despite the important role of OXTR during parturition, little has been published about epigenetic regulation of the oxytocin receptor gene, OXTR, in connection with the normal process of parturition including postpartum uterine contractility. Located on chromosome 3, OXTR consists of three introns and four exons40. Previous research has shown that DNA methylation of a region contained within a CpG island in the human OXTR promoter (termed MT2) plays a key role in regulating OXTR transcription41. In humans and animal (prairie vole) models, our team has shown that both brain and blood-derived DNA methylation in this region is correlated with gene expression in the brain. Blood methylation accounts for approximately 18% of the variance in expression42, suggesting to some extent that blood-derived DNA methylation measurements in OXTR may be a proxy for gene expression in inaccessible tissues42,43.
Researchers examining OXTR DNA methylation (OXTRm) have found that various psychopathologies are correlated to OXTRm dysregulation44–46. Others have reported variation in social/emotional outcomes in connection with OXTRm47,48. We and others have also identified a role for OXTRm in the risk for postpartum depression44,49, preterm labor2,50, maternal antidepressant use51, and exposure to early-life adversity/ post-traumatic stress52–57. Autism and social cognition have also been linked to OXTRm58–63. Finally, promoter OXTRm has also been found to be sensitive to early caregiving experience, such that higher levels of parental care have been associated with reduced DNA methylation in both human infants64 and prairie vole pups42.
Understanding how OXTRm also relates to OXTR function and subsequent uterine contractility during parturition would have important clinical implications both physiologically and pharmacologically. Oxytocin administration is the primary pharmacological intervention used during the course of labor65,66 and is frequently used for labor induction67, rates of which have risen from 9.5% of births in 199068 to 31.4% of births in 202069. Intrapartum dose requirements for oxytocin have been found to vary by body mass index (BMI)70,71 and maternal age72, and may be influenced by differences in lipid metabolism73. While use of oxytocin for labor stimulation may be a life-saving intervention, overuse can expose the mother to risks, including postpartum hemorrhage74–76. Intrapartum oxytocin use leads to uterine oxytocin receptor desensitization (and internalization of the receptor from the myocyte membrane) and/or downregulation of OXTR21,22,77,78. These mechanisms result in diminished capacity of the uterus to contract effectively after birth and predisposes the individual to a higher likelihood of uterine atony and postpartum hemorrhage. Rates of atony-associated postpartum hemorrhage appear to be specifically rising79, which is tied to increased intrapartum oxytocin administration80 even among healthy, low-risk vaginal births81. Oxytocin also remains a first-line therapy for stimulating uterine contraction after birth for prevention and treatment of postpartum hemorrhage82.
Other than three preterm labor focused studies already mentioned, no studies have linked OXTRm to maternal oxytocin response or uterine contractility throughout parturition; however, extracting myometrial tissue would be impractical and unethical following vaginal birth. A peripheral tissue (e.g., whole blood) could be used to understand uterine OXTRm patterns if the tissue methylation levels were well-correlated with one another. Our previous work has reported significant correlations between OXTRm derived from brain and blood in prairie voles42 as well as blood and saliva in humans;61 however, no studies have reported correlation of blood and myometrial DNA methylation.
In the current study, we address this scientific/clinical gap by examining epigenetic variability in the endogenous oxytocin system. First, we measure OXTRm between matched uterine myometrium and maternal blood samples using an exploratory array method as well as targeted pyrosequencing to determine if blood can be used as an indicator of uterine OXTRm levels. Second, we test differences in OXTRm in relation to oxytocin administration during the process of parturition and indices of postpartum hemorrhage. Considering that postpartum hemorrhage is commonly attributed to uterine atony, reduced availability of OXTR might associate with greater postpartum hemorrhage risk. We therefore hypothesize that higher levels of OXTRm, reflective of less OXTR availability, might associate with greater blood loss after giving birth. In this study, we find that OXTRm in myometrium correlates with the level of DNA methylation from blood and that higher OXTRm relates to both greater oxytocin dosage needed during the birth process and more blood loss after vaginal birth. Findings from this study may further our understanding of the relationship between the endogenous oxytocin system and birth outcomes, inform clinical risk assessment, or help improve use of oxytocin administration.
Materials and methods
Correlation betweenm using matched uterine myometrium and maternal blood OXTR
Participants
Participants’ data in this study were available via a pregnancy-related tissue repository managed by Oregon Health and Science University. Potential participants were engaged by trained research staff, and after providing written informed consent, donated samples of myometrium, whole maternal blood, placenta and cord blood. The pregnancy repository was approved by the ethics committee at Oregon Health and Science University. Limited clinical information was available in the repository and no identifying information, therefore the committee found it did not constitute research with human subjects. Samples available in the repository were from participants who gave birth at two different institutions, (Oregon Health and Science University and Chelsea and Westminster Hospital, London). Material data transfer and data use agreements were established between institutions, de-identified samples were collected from (Chelsea and Westminster Hospital) and included in an existing tissue repository at (Oregon Health and Science University).
DNA obtained from whole blood and myometrium at time of planned cesarean delivery
Delivery at 37 weeks or greater was indicated due to prior cesarean, choice, or presence of chronic or gestational hypertension. In the hour prior to delivery, whole blood was collected in standard EDTA tubes and frozen in 1.5 mL aliquots at −80 °C. After birth of the newborn and placenta, a small portion of full-thickness myometrium (measuring approximately 2 cm × 2 cm in size) was excised from the hysterotomy incision and was handed off the sterile field to study staff, described by McElvey et al. (2000)83. Myometrium was rinsed in phosphate buffered saline, dissected from serosa and endometrium, if present, and flash frozen in liquid nitrogen in small sections. Uterine specimens were not collected in an RNA preservation medium. Specimens were collected from 27 participants. DNA from blood and myometrium were isolated, bisulfite converted and followed by both Illumina array processing and pyrosequencing to assess the methylation across the OXTR as well at the MT2 locus of interest (−934) (complete methods described below).
Blood-derivedm in relationship to obstetric oxytocin requirements during parturition and postpartum hemorrhage (case-control study) OXTR
Postpartum hemorrhage cases were matched to participants with normal/typical blood loss after birth. Participants were matched by induction of labor status (vs. spontaneous labor onset) and parity. The protocol for this study underwent ethics review at (Oregon Health and Science University) and was approved as human research.
Participants
Participants gave birth in the greater Portland, Oregon area including hospitals and community birth settings (home/birth center) from November 2018 through February 2020. Importantly, study recruitment and planned study visits for blood sampling were halted prematurely due to SARS-CoV-2 pandemic restrictions on in-person research activities. Those giving birth at our institution were given information about the study in the postpartum unit or at a follow-up visit based on meeting basic screening criteria in the electronic health record. Other participants responded to advertisements posted in the community and on social media.
Enrollment criteria included the following: singleton vaginal birth, 37 weeks or greater gestational age, maternal age 15–45, able to read and provide informed assent/consent in English or Spanish. Exclusions included: known coagulopathy, hemolysis, elevated liver enzymes and low platelet (HELLP) syndrome, disseminated intravascular coagulation, magnesium sulfate administered during labor, placenta accreta, vasa previa, gestational thrombocytopenia (< 80,000 platelets/mL) or severe bleeding noted from genital tract injury during birth. Participants completed baseline survey measures upon enrollment after a process of informed written consent was conducted. Recruitment goal of 100 cases and 100 matched controls was based on a mean difference between groups of a magnitude of two percentage points (SD of 5.0) with 80% power and α of 0.05. Of the 577 postpartum records screened, we approached 393 individuals and 120 met criteria and consented to participate in the study.
Cases of postpartum hemorrhage and characteristics of uterine atony
We enrolled participants as cases who had experienced a diagnostic postpartum hemorrhage, a cumulative blood loss ≥ 1000 mL in the first two hours after birth, when bleeding was attributed to lack of uterine tone. People with significant bleeding from lacerations were excluded. We also enrolled those who experienced uterine atony with blood loss totaling at least 400 mL in the first two hours after birth and had additional clinical interventions including secondary uterotonic medication to treat heavier than normal bleeding (not attributed to lacerations). These interventions may have led to avoiding an official postpartum hemorrhage diagnosis. Further rationale for this expanded definition of postpartum hemorrhage is supported by other research indicating that providers commonly underestimate the volume of blood lost during birth84,85.
This phenotyping process required extensive review of the electronic medical record for nurses’ notes and the birth attendant’s (midwife/physician) written birth narrative, in addition to viewing the medication administration records and postpartum notes. Care providers reported blood loss in the health record as either measured (e.g., blood loss was weighed or measured using graduated drapes) or by estimation methods (visual inspection of drapes and towels). We also examined post-birth hemoglobin/hematocrit values to validate the quantity of blood lost, when it was available in the chart.
Medical record abstraction and oxytocin administration quantification
The electronic medical record was abstracted for the dose, start and stop times of any oxytocin administered during labor, along with the timing and duration of titration of oxytocin doses. The dosing of oxytocin is quantified in milliunits (mU)/min; therefore, the total cumulative dosage was a product of the duration (min) of a given level of titration and the oxytocin dose (i.e., 60 min x 2 mU/min = 120 mU). Highest dose given along with duration of oxytocin, other medications given during labor and any prenatal, intrapartum or postpartum complications was also abstracted. Normal or healthy labors will not typically require augmentation with oxytocin, or oxytocin will be used for only a short time, for example, following regional analgesia (epidural) placement. In the presence of labor induction or need for labor augmentation, oxytocin regimens typically involve starting at a 1 mU/min administration rate and increasing by 1–2 mU/min every 20–40 min86. Significant institutional variation as well as situational variability in dosing exists, however.
Treatments and interventions to control atony and postpartum hemorrhage
Similarly, oxytocin given after birth for postpartum hemorrhage prophylaxis (active management of the third stage of labor) was recorded in addition to treatments for heavier bleeding including any of the following: intravenous (IV) oxytocin, intramuscular oxytocin, misoprostol, methylergonovine, carboprost tromethamine, bimanual compression, uterine massage, bladder management, manual removal of retained clots or placental/membrane fragments. Routine active management of third stage labor may consist of 10 U of oxytocin given via intramuscular injection or through the intravenous (IV) line. However, some institutions will instead utilize the existing IV solution of oxytocin being used for labor management, which may be 30 U oxytocin in 500 mL lactated ringers (however, some institutions may use 40 U in 1000 mL)66. If bleeding is ongoing, an IV solution of oxytocin would likely be ordered and administered rapidly. The medical record was examined for the dose/duration of the postpartum infusion (if any was used) as well as any intramuscular doses. Some births did not have any postpartum oxytocin administration, which may occur if a patient declined the medication and/or the practitioner assessed the patient to be low-risk for postpartum hemorrhage. In addition, the postpartum records were reviewed for IV iron therapy or blood transfusions as well as delayed postpartum hemorrhage (> 2 h after birth). In the enrollment survey, participants were asked to self-report their pregnancy, birth and medical history, which was used to confirm and validate medical record abstraction findings.
Case-control maternal blood-derived DNA sampled 6–10 weeks after birth
Timing of blood sampling was delayed until 6–10 weeks after birth to limit chimerism with any DNA/blood cells received via blood transfusion among participants with postpartum hemorrhage. Whole blood samples were obtained from participants using vacutainer EDTA tubes. Samples were stored immediately at −80 °C. Participants received remuneration in the form of a $20 gift card at the time of the study visit. A total of 91 blood samples were available for DNA extraction. Several participants (n = 6) could not have blood sampling due to pandemic restrictions and another 23 enrolled participants did not show for the study visit and could not be contacted for rescheduling. DNA isolated for pyrosequencing occurred using reagents and procedures outlined in the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). All DNA was quantified using Nanodrop. One-hundred nanograms (ng) of isolated DNA were subject to bisulfite conversion using MECOV50 kits (Thermo Fisher Scientific, Waltham, USA). Sample locations within a 96-well plate were determined using an online random number generator (random.org).
Determination of targeted CpG methylation by pyrosequencing
Following bisulfite conversion, the remaining procedures were performed in triplicate: Forty nanograms of bisulfite-converted DNA was amplified using polymerase chain reaction (PCR) with PyroMark PCR Kits (Qiagen) and 0.2 µM primers (Forward (TSL101F): 5′-TTGAGTTTTGGATTTAGATAATTAAGGATT-3′; Reverse (TSL101R): 5′-biotin-AATAAAATACCTCCCACTCCTTATTCCTAA-3′). Methylation standards (0, known, and 100% methylated) and negative controls from bisulfite conversion and PCR were included on the PCR plate. The following cycling protocol occurred simultaneously on three identical satellite thermocyclers (Bio-Rad, Hercules, USA): Step 1) 95 °C, 15 min; Step 2) 50 cycles of: 94 °C, 30 s; 56 °C, 30 s; 72 °C, 30 s; Step 3) 72 °C, 10 min; Step 4) 4 °C hold. The amplification of a 116-base-pair region on the coding strand of OXTR containing CpG −934 and −924 (hg38, chr3: 8,769,044 – 8,769,160) and no contamination in negative controls was confirmed by 2% agarose gel electrophoresis.
DNA methylation level for each sample was assessed through pyrosequencing using PyroMark Gold Q24 reagents on a PyroMark Q24 machine (Qiagen) (sequencing primer TSL101S: 5′-AGAAGTTATTTTATAATTTTT-3′). Samples with high replicate variability were identified as having a Z-score above +/−1.96. Eight samples were identified on the basis of this criterion. These 8 samples were visually inspected; when one replicate value significantly deviated from the other two values, this value was removed (n = 7). Replicate variability, assessed via mean deviation, averaged +/−1.33% for the tissue match experiment (blood = +/−1.33%, tissue = +/−1.33%) and +/−1.29% in the case-control study. Apart from samples reporting on two replicate values, reported epigenotypes are the average of three replicate values. In addition, 0% and 100% DNA methylation controls as well as a sample of known DNA methylation level, were run alongside experimental samples with expected methylation results (0% OXTRm mean = 3.57% +/− 0.93%, 100% OXTRm mean = 99.06% +/− 0.78%, and known 42.5% +/− 1.88%).
Illumina Array (myometrial tissue and matched blood only)
Isolated DNA was quantified by fluorescence using Quant-iT™ PicoGreen ® dsDNA assay kit (Invitrogen | Thermo Fisher Scientific). DNA methylation was assessed using Illumina’s Methylation EPIC BeadChip (Illumina, San Diego, USA). 500 ng of each sample was bisulfite converted using an EZ DNA Methylation Kit (Zymo), amplified, hybridized, and imaged. DNA methylation data for over 850,000 CpGs was generated per sample and preprocessed using R statistical suite (version 3.6.1).
Raw.idat files were read and preprocessed using the minfi R package87,88. The data set was preprocessed using noob for background subtraction and dye-bias normalization. All methylation values with detection P > 0.01 were set to missing. One uterine sample had excessive missing probes (809,422, 93.44%) and was excluded from further analysis. Of the remaining samples, the number of missing probes ranged from 847 to 5008 (median: 1106). Probes with > 1% missing values (n = 12,399) were removed from further analysis. Additionally, samples were checked for unusual cell mixture estimates using the using the Houseman method as implemented in minfi89,90. All samples passed these quality controls. Principal components analysis, as implemented in the shinyMethyl R package, was used to examine batch effects91. The first seven principal components were examined using plots and potential batch effects were tested using linear models. Principal component 1, which accounted for 79.47% of the total variance, was associated with tissue type (F(1, 50) = 5149, p < 2.2e-16, adjusted R2 = 0.992). Location of specimen collection site (London, UK or Portland, OR) was associated with principal components 2, 3, 6, and 7, which account for 1.57%, 0.95%, 0.72%, and 0.69% of the total variance, respectively (PC2: F(1, 50) = 4.34, p = 0.04, adjusted R2 = 0.061; PC3: F(1, 50) = 4.57, p = 0.038, adjusted R2 = 0.066; PC6: F(1, 50) = 6.48, p = 0.014, adjusted R2 = 0.097; PC7: F(1, 50) = 16.04, p = 2.06e-4, adjusted R2 = 0.228). Array was associated with PC2 (F(6, 45) = 30.97, p = 2.01e-14, adjusted R2 = 0.779). Position on the array was not associated with any of the first seven principal components. Because proportions of variance explained by associated principal components were low, no batch correction method was used. The ewastools R package was used to assess Illumina quality control metrics and call genotypes and donor IDs to ensure the identity of matched samples from the same individual92. One uterine sample failed the non-polymorphic green control, but was retained as it was just below threshold and did not appear abnormal on the other 16 metrics. All other samples passed all Illumina quality controls. Following quality controls, background subtraction and normalization were performed using the preprocess Funnorm function in the minfi package, which is recommended for between-tissue studies93. Normalized beta values were then extracted for probes that fall within OXTR and probes within 3 kb of the gene (hg19 chr3: 8,789,094-8,814,314). There are 23 such probes, but one probe (cg08449558) was excluded because it includes a SNP with minor allele frequency of 9%.
Statistics and reproducibility
DNA methylation levels between matched blood and uterine specimens were assessed with a Spearman correlation. Paired t-tests (two-sided) were used to compare mean differences in level of methylation observed in each tissue. To correct for multiple comparisons, we applied a Benjamini-Hochberg correction to the p-values from the correlation analysis.
For the case-control study, differences in sample characteristics were measured between cases and controls using parametric and non-parametric tests as appropriate (Spearman, two-sided t-tests or χ2). We also examined differences in OXTRm by baseline characteristics (e.g., age, body mass index, parity). Our primary hypothesis was that higher OXTRm would be associated with higher blood loss and postpartum hemorrhage occurrence and that this relationship would be moderated by oxytocin use during labor. We therefore tested this hypothesis using an interaction term of oxytocin use (versus no oxytocin use) in the regression models.
The intent for this study was a 1:1 match of hemorrhage cases with control (typical postpartum bleeding) participants. However, in the end, cases outnumbered controls as recruitment and study visits stopped in February of 2020. We therefore analyzed our data comparing cases to control participants in addition to using the volume of blood loss as a continuous outcome. First, we compared DNA methylation between cases and control groups with linear regression. Controlling for mean replicate deviation, we then examined this association with or without oxytocin use in labor. We then used generalized linear models with a gamma distribution to examine the outcome of the quantity of oxytocin used across the phases of birth (oxytocin quantity was significantly right skewed) (intrapartum and postpartum) with OXTRm as the primary independent variable controlling for mean deviation (across replicates). We included parity and body mass index at birth as confounders in the models as they significantly differed between cases and controls and both are associated with oxytocin use during childbirth.
In the second model, we used adjusted Poisson regression with robust standard errors to assess the relationship of OXTRm to the case or control status. The third model, a multivariable generalized linear model (GLM), assessed the total blood lost during the birth process using a gamma distribution (blood loss is skewed). The fourth and fifth models used adjusted Poisson regression with robust standard errors for the outcomes of postpartum hemorrhage (using the 500 mL definition) and again at the 1000 mL definition. For each of these models we included the interaction term of oxytocin use during labor (versus no oxytocin) and controlled for replicate variability, parity, antibiotic use during labor (to control for possibility of intrapartum infection) and body mass index at delivery. Both the traditional and updated definitions of postpartum hemorrhage were used, given the known problems with blood loss underestimation94.
Reporting summary
Further information on research design is available in thelinked to this article. Nature Portfolio Reporting Summary
Results
Uterine myometrium and maternal blood are correlated at multiple CpG sites in OXTR
Matched blood and myometrium samples from 26 participants were included in the analysis. The mean (standard deviation) age of the sample was 35.3 (5.9) years. Participants delivered via planned, pre-labor Cesarean with a mean (standard deviation) gestational age of 38.7 (0.77) weeks. The majority were multiparous (85%) with a mean (standard deviation) pre-pregnancy BMI of 27.8 (7.8) kg/m2.
We analyzed myometrium and blood-derived DNA methylation at CpG site sites across the OXTR locus including two previously reported and conserved sites, −924 and −934, located in the promoter region of OXTR (Fig. 1). Data were generated using array-based methods and targeted pyrosequencing.
Using data from the Illumina methylome array, 22 probes within OXTR and its promoter were examined for correlation between maternal blood and myometrium using Spearman’s rho and Benjamini-Hochberg correction (Fig. 1 denotes specific loci and Fig. 2 shows significant sites). Of these 22, four sites in the 3’ region of the gene were significantly correlated between tissues cg03257388 (ρ(24) = 0.66, p = 0.005), cg02192228 (ρ(24) = 0.62, p = 0.008), cg11171527 (ρ(24) = 0.55, p = 0.02) and cg15317815 (ρ(24) = 0.59, p = 0.01). In addition, two CpG sites not present on the array, but shown to partially account for gene expression and highly reported on in the literature, were assayed by pyrosequencing and found to be correlated (−934: ρ (24) = 0.65, p = 0.005) (−924: ρ(24) = 0.52, p = 0.03). Given that site −934 (chr3: 8,810,729-8,810,845) demonstrated both high statistical significance between tissues, is associated with differences in transcription of OXTR, and is a commonly studied CpG site in human adult literature43,95, we chose to focus on this site for our examination of oxytocin administration and postpartum hemorrhage outcomes in the case-control study.
gene schematic and location of CpG sites with significantly correlated methylation. OXTR Gene schematic of. Boxes represent exons and lines represent introns. Coding regions are in white and untranslated regions are in gray. The black arrow indicates the transcription start site. The black bar below the gene denotes the CpG island. The 5’ region of the gene is enlarged to show detail. Blue lines indicate CpG sites assayed by pyrosequencing and red lines indicate CpG sites assayed by Illumina Methylation EPIC 850 K array. CpG sites with significantly correlated DNA methylation in blood and myometrium are indicated by letters and Spearman’s rho is reported below the gene schematic. OXTR
DNA methylation of CpG sites inare correlated in blood and myometrium. OXTR Significantly correlated CpGs between tissues (= 26 participants with matched samples)< 0.05 after correction for multiple comparisons using the Benjamini-Hochberg method: −924 (red): ρ = 0.52,= 0.03, −934 (light blue): ρ = 0.65,= 0.005, cg02192228 (green): ρ = 0.62,= 0.001, cg03257388 (dark blue): ρ = 0.66,= 0.005, cg11171527 (apricot): ρ = 0.55,= 0.02, cg15317815 (lavender): ρ = 0.59,= 0.01.Differences in methylation levels at each CpG site between myometrium and blood, stars denote Bonferroni corrected significant differences: −934 (t(25) = 10.1,= <0.001), −924 (t(25) = 9.6,< 0.001), cg15317815 (t(24) = 3.2,= 0.004), cg02192228 (t(24) = −4.3,= 0.003), cg03257388 (t(24) = 1.12,= 0.27), cg11171527 (t(24) = 2.2,= 0.04). Whiskers denote +/−1.5 x interquartile range. A B n p p p p p p p p p p p p p
Postpartum hemorrhage cases were administered higher dosage of oxytocin used throughout labor and postpartum
Using a case-control study design, we enrolled a sample of 69 postpartum hemorrhage cases and 50 control participants. Blood samples were available for DNA isolation from 55 of the cases and 36 of the control participants. Table 1 denotes sample characteristics and differences between cases and controls. While we did match for parity during enrollment, nulliparous participants were more highly represented among cases compared to controls (57.9% vs 38.0%, p < 0.05) due to halted enrollment. The range of total blood loss was 400–2540 mL among cases and 50–600 mL for control participants. No differences were seen between cases and controls in maternal age, gestational age, and presence of antenatal complications (hypertension/diabetes). Participants self-identifying as Latin American / Hispanic ethnic background/ancestry were more frequent among the cases (n = 11,15% of cases, vs. n = 1, 2% of controls, p < 0.05). Other differences included use of antibiotics during labor and BMI, with postpartum hemorrhage cases having significantly higher BMI at delivery (31.9 vs. 29.9, p = 0.03). Finally, the total quantity, duration and maximum dosage of oxytocin needed during labor was higher among cases of postpartum hemorrhage. As expected, postpartum use of oxytocin was also significantly higher among cases, though the use of prophylactic oxytocin (active management of third stage labor) was not different between groups. The cumulative oxytocin dosage for cases was more than double that of controls (25.6 Units vs. 12.4 Units, p < 0.001).
| Case (= 69) Mean (SD)n | Control (= 50) Mean (SD)n | |
|---|---|---|
| Maternal age, years | 31.4 (4.9) | 32.4 (4.4) |
| Body Mass Index at delivery, kg/m2 | 31.9 (5.4) | 29.9 (4.2) * |
| Gestational age, weeks | 39.7 (1.3) | 39.8 (1.2) |
| Third trimester hemoglobin mg/dl | 11.8 (1.2) | 11.9 (1.2) |
| Length of first stage labor, hours | 16.4 (11.3) | 11.6 (7.8) * |
| Length of second stage labor, hours | 1.6 (1.6) | 1.3 (1.5) |
| Length of ruptured membranes, hours | 9.7 (10.9) | 6.7 (6.8) |
| Total oxytocin if administered intrapartum, Units | 9.9 (9.4) | 3.9 (4.7) † |
| Total oxytocin duration intrapartum, hours | 16.3 (9.7) | 10.5 (7.5) † |
| Highest maximum dose of oxytocin, mU/min | 13.9 (8.1) | 9.1 (6.1) † |
| Postpartum oxytocin used, Units | 21.5 (14.6) | 12.5 (9.1) † |
| Total oxytocin used intrapartum/postpartum combined, Units | 25.6 (16.2) | 12.4 (9.4) ‡ |
| Infant size, grams | 3410.7 (568.5) | 3484.6 (689.9) |
| Length of third stage of labor, minutes | 10.0 (9.5) | 9.0 (10.4) |
| Total volume blood loss, mL | 813.9 (511.9) | 277.1 (142.8) ‡ |
| Total volume blood loss, range | 400–2540 | 50–600 |
| Difference in hemoglobin (third trimester to lowest postpartum value), mean (SD) mg/dl | −1.9 (2.3) | 0.05 (1.23) ‡ |
| Case (= 69)(%)nn | Control (= 50)(%)nn | |
| Primiparous | 40 (57.9) | 19 (38.0) * |
| Self-reported ancestry/race/ethnicitya | ||
| European | 55 (79.7) | 39 (78.0) |
| Latin American | 11 (15.9) | 1 (2.0) * |
| Asian | 7 (10.1) | 6 (12.0) |
| African | 1 (1.4) | 0 (0.0) |
| Other | 4 (5.8) | 3 (6.0) |
| Bleeding during first trimester pregnancy | 7 (10.0) | 1 (2.0) |
| Prior Cesarean birth | 6 (8.8) | 4 (8.0) |
| Tobacco used prior to/during pregnancy | 4 (5.8) | 7 (14.0) |
| Gestational diabetes | 13 (18.5) | 5 (10.0) |
| Mild gestational thrombocytopenia (platelets 80,000–120,000/μl prior to delivery) | 4 (5.8) | 2 (4.0) |
| Hypertensive disorder diagnosis | 6 (8.7) | 4 (8.0) |
| Location of birth | ||
| Hospital | 68 (97.1) | 47 (95.9) |
| Community (home/birth center) | 2 (2.9) | 2 (4.0) |
| Group B streptococcus carrier | 22 (31.9) | 8 (17.0) |
| Infection/ chorioamnionitis suspected during labor | 8 (11.6) | 2 (4.0) |
| Antibiotics administered | 30 (43.5) | 9 (18.8) † |
| Epidural use in labor | 49 (71.0) | 30 (61.2) |
| Spontaneous onset labor | 34 (49.3) | 26 (52.0) |
| Oxytocin used in labor | 42 (60.9) | 28 (56.0) |
| Instrument assisted vaginal birth | 5 (7.2) | 1 (2.1) |
| Genital trauma (none/no repair needed) | 17 (24.3) | 14 (29.8) |
| Female infant | 37 (53.6) | 21 (42.0) |
| Active management of third stage labor (oxytocin) | 58 (84.1) | 41 (82.0) |
| Timing of cord clamping < 30 seconds | 8 (12.5) | 3 (6.9) |
| Received iron infusion or blood transfusion postpartum | 13 (18.8) | 1 (2.1) † |
| Available DNA specimens | 55 (80.0) | 36 (71.4) |
Maternalm was lower among cases of postpartum hemorrhage following labor without oxytocin OXTR
Average OXTRm across the 91 participants was 45.2% (SD = 6.5), ranging from 20.3% to 63.2%. Mean (SD) % OXTRm for cases was 44.9 (7.2) and 45.4 (5.4) for control participants, which was not significantly different (z = −0.02, p = 0.99, Wilcoxon rank-sum). Figure 3 illustrates the DNA methylation differences between cases and controls measured by Wilcoxon rank-sum tests when oxytocin was used in labor (n = 53) z = −1.86, p = 0.06, or when labor progressed without oxytocin stimulation (n = 38) z = 2.84, p = 0.003. Using a linear regression model and controlling for pyrosequencing replicate variability, we found OXTRm to be significantly lower in cases compared to controls (n = 22 cases and 16 controls) when labor occurred without oxytocin stimulation (4.15% lower; 95% CI −7.91 to −0.39) (Table 2). Parity, maternal age, gestational age, and maternal BMI at delivery were not associated with OXTRm in bivariate analyses (p-values > 0.05).
Cases of postpartum hemorrhage had lower DNA methylation when labor was not stimulated by oxytocin. Star denotes significant difference in total DNA methylation between cases of postpartum hemorrhage cases and controls among participants whose labor was not stimulated with oxytocin (red box plots,= 38),= 2.84,= 0.003.Difference in DNA methylation between postpartum hemorrhage cases and controls when oxytocin was used in labor (blue box plots,= 53) z = −1.86,= 0.06. Whiskers denote +/−1.5 x interquartile range. A B n z p n p
| Cases vs. Controls β(95% CI) | |
|---|---|
| mOXTR | |
| Sample without interaction (= 91)n | −0.50 (−3.27 to 2.28) |
| Oxytocin during labor (= 53)n | 2.16 (−1.76 to 6.09) |
| Without oxytocin in labor (= 38)n | −4.15 (−7.91 to −0.39) * |
| Interaction: Oxytocin in labor x Case vs. Control (= 91)n | 6.29 (0.80 to 11.79) * |
| Sample Characteristics β(95% CI) | |
| mOXTR | |
| Maternal age | 0.08 (−0.23 to 0.39) |
| Parity | 0.07 (−1.10 to 1.24) |
| Body Mass Index | 0.09 (−0.18 to 0.37) |
| Gestational Age | −0.52 (−1.64 to 0.61) |
| Infant sex | −0.17 (−2.90 to 2.56) |
| Smoking before/during pregnancy | 0.71 (−4.10 to 5.50) |
| Antepartum anemia | −0.49 (−3.73 to 2.74) |
Oxytocin administration during parturition is associated with greaterm OXTR
Cumulative oxytocin administration (intrapartum and postpartum total in Units) was correlated to higher OXTRm (ρ(87) = 0.28, p = 0.008). Among participants who needed oxytocin administered during labor, the correlation between total dosage and OXTRm was higher (ρ(50) = 0.42, p = 0.002). As the distribution of oxytocin use in labor was significantly right-skewed (Shapiro-Wilk W test for normality W = 0.87, p < 0.0001), we used a GLM with a gamma distribution to examine oxytocin use by OXTRm in lieu of a linear regression or transforming the data. We found that with an increase in OXTRm of 5 percentage points, 2.48 more units of oxytocin were required for management of parturition (stimulation of labor and treating postpartum bleeding) (β = 2.47, 95% CI 0.66-4.28, p = 0.008). We controlled for pyrosequencing replicate variability, parity, and maternal BMI at time of delivery in the model (Table 3, Fig. 4).
IncreasedDNA methylation is associated with increased oxytocin use and total postpartum blood loss only when oxytocin is administered during labor. OXTR Red dots show participants without oxytocin use in labor (= 38); blue dots denote oxytocin was used in labor (= 53). Shaded area denotes 95% CI of regression line.Higher −934m is associated with higher intra/postpartum administration of oxytocin among participants receiving oxytocin (blue dots).Postpartum bleeding volume associated with DNA methylation however, severity increased when oxytocin was used in labor. n n OXTR A B
| Total oxytocin administered (Units) | β(95% CI) |
|---|---|
| m (5% increase)OXTR | 2.47 (0.66 to 4.28) † |
| Replicate variability | −0.98 (−30.69 to 28.73) |
| Parity (total prior births) | −2.64 (−4.67 to −0.59) * |
| BMI at delivery | 0.82 (0.14 to 1.51) * |
Higherm is associated with greater relative risk for postpartum hemorrhage when oxytocin was used during labor OXTR
Examining only participants with oxytocin use during labor, we found higher OXTRm was associated with greater blood loss, ρ(51) = 0.31, p = 0.02. However, among those labors not stimulated with oxytocin, higher OXTRm correlated with lower blood loss (ρ(36) = −0.34, p = 0.04) (Fig. 4). Given that oxytocin use appeared to be an effect modifier, we used Poisson regression models and GLMs with interactions across the entire sample. In the regression models, the association between OXTRm and blood loss remained across the hemorrhage outcomes. When oxytocin was used during labor, each 5-percentage point increase in OXTRm level was associated with 45% higher relative risk of being a case compared to a control (IRR 1.45, 95% CI = 1.11–1.91). For each 5-percentage point increase in OXTRm level, 185.22 mL greater blood loss was predicted in the GLM (β = 185.22, 95% CI 70.53-299.91). Increased OXTRm was also associated with higher relative risk for postpartum bleeding of 500 mL or greater, RR 1.45 (95% CI 1.12–1.89). Finally, the combination of higher OXTRm and oxytocin use in labor was associated with a RR of 2.95 for postpartum hemorrhage (1000 mL or greater), 95% CI 1.53–5.71. All models controlled for replicate variability, parity, antibiotic use during labor (a proxy for possible intrapartum infection) and BMI at delivery (Table 4).
| Study Case (vs. Control) IRR(95%CI) | Cumulative Volume Blood Loss (mL) β(95%CI) | Postpartum Hemorrhage (500 mL) IRR(95%CI) | Postpartum Hemorrhage (1000 mL) IRR(95%CI) | |
|---|---|---|---|---|
| Main effects | ||||
| m levelOXTR | ||||
| (5% increases) without oxytocin in labor | 0.74 (0.59–0.94)† | −109.71 (−207.5 to −11.95)* | 0.71 (0.56–0.90)† | 0.56 (0.35–0.91)* |
| Oxytocin used in labor vs. none and lowmOXTR | 0.04 (0.003–0.39)* | −1567.57 (−2620.14 to −515.01)† | 0.05 (0.006–0.51)* | 0.001(2.00to 0.09)†−7 |
| Interactionm (5%) X oxytocin useOXTR | 1.45 (1.11–1.91)† | 185.22 (70.53–299.91)† | 1.45 (1.12–1.89)† | 2.95 (1.53–5.71)‡ |
Discussion
In this study, we assessed the relationship of OXTRm in pregnant uterine myometrium obtained from maternal blood and in a separate experiment tested the association of OXTRm with an individual’s need for greater uterine stimulation using oxytocin and subsequent postpartum hemorrhage. The first key finding is that pregnant uterine and blood OXTRm are significantly correlated at several sites, providing support for the use of blood OXTRm as a surrogate for uterine OXTRm. Second, we found a significant positive association between OXTRm and greater cumulative oxytocin utilization during parturition. Last, we provide evidence that the severity of postpartum blood loss was associated with higher OXTRm and moderated by oxytocin administration. To the best of our knowledge, this is the first study to examine these associations, which has relevance to understanding an important source of maternal morbidity.
Despite recognized tissue-specific DNA methylation patterns41, our finding that OXTRm derived from uterine myometrium and maternal blood are correlated is encouraging for future use of peripheral blood DNA samples in the ongoing investigation of uterine function during childbirth. Clearly, sampling uterine tissues during childbirth can only be done routinely during surgical (Cesarean) birth, which will necessarily limit the interpretation of the findings to vaginal births. Furthermore, future clinical utility of this kind of biomarker would require it be accessible prior to birth. Future work may determine if salivary OXTRm will correlate to uterine levels, as it does with blood61,62, which would allow for non-invasive measurement of this marker in uterine function. In addition, future study of the variation in OXTR expression as a result of OXTRm is needed to connect methylation to the functional changes in uterine tissue.
While our work focused only on OXTRm, the potential consequences of higher myocyte OXTRm include decreased OXTR transcription, reduced OXTR density or less effective intracellular signaling (leading to less efficient uterine contraction) leading to diminished pharmacologic responses. In other human studies, differential promoter OXTRm has been documented in samples of maternal decidua (higher with preterm labor)50, in fetal membranes (higher with preterm birth)58 and in myometrium at term without labor (lower methylation than preterm)2. One study in mice showed an association in uterine tissue of higher Oxtr promoter methylation and lower mRNA expression, though we note that the OXTR CpG site we focused on (−934) is not conserved in mouse models, which may limit translation to human interpretation96. Other studies linking higher promoter OXTRm to reduced gene transcription have been performed in brain tissue42,95. Future research is necessary to examine in vitro uterine tissue contractility in relationship to OXTRm and oxytocin administration. In concordance with our findings, other studies have shown a relationship between exogenous oxytocin and diminished physiologic responses with higher OXTRm. Among female participants given intranasal oxytocin, Chen et al. found fMRI-measured responses in visual cortex were higher in individuals with lower promoter OXTRm assessed from saliva97. Furthermore, greater promoter OXTRm derived from placental vasculature was associated with diminished vasoconstriction in response to exogenous oxytocin98. In sum, our study links the peripheral OXTRm to key obstetric outcomes as a biomarker, yet requires further investigation in terms of functional differences at the tissue level.
Prior to the hospital admission, obstetric care providers have no indication of how effectively oxytocin will stimulate uterine contractions. Current practice includes titration of the medication by 1–2 mU/min (4–6 mU/min if using a high dose approach) until uterine activity appears regular and the cervix dilates66. For some, this may take only a few hours, while others may require many hours or days. Mousa et al. (2008) reported in a sample of over 20,000 births that among individuals who developed postpartum hemorrhage, 10% of women having Cesarean birth and 21% of those with vaginal births failed to respond to uterotonic medication99. Whether epigenetic variability of OXTR might be responsible for the inability to respond to first line treatments due to lower receptor density has clear clinical relevance, given the variability in clinical oxytocin response both during and after labor. It is also well-established that higher oxytocin use/dosing is needed for people with higher BMI100, however, whether OXTRm plays a role between oxytocin needs and body mass is unknown. Our analyses did account for BMI at the time of birth, which suggest that OXTRm may be a useful predictor of oxytocin needs regardless of body characteristics. More research into DNA methylation (oxytocin sensitivity) could help explain current clinical challenges with labor inductions101 or lead to developing new personalized strategies for labor management or postpartum hemorrhage prevention/treatment. Possible clinical improvements might include using personalized dosing regimens during labor or alternative uterotonics as first line postpartum hemorrhage prophylaxis. More work will be needed to prospectively test OXTRm measured before labor starts for the utility of this biomarker in predicting oxytocin responses both during labor and after giving birth.
Future research should addressm stability over time OXTR
An important consideration when interpreting these findings is whether OXTRm is stable throughout a short time frame in parturient females (i.e., birth to 6 weeks postpartum). Given that our sampling was done after birth, we cannot claim that higher OXTRm was the mechanism responsible for greater oxytocin needs at parturition. The alternative testable hypothesis is that greater oxytocin usage led to greater OXTRm postpartum, an effect magnified by postpartum hemorrhage or higher blood loss, because far more oxytocin would be needed in this scenario. Evidence for this line of reasoning among adult humans is limited; however, oxytocin exposure in the prairie vole has been demonstrated to cause conserved promoter OXTRm and gene transcription differences102. Assessing OXTRm pre-labor as well as postpartum would aide in answering this question. If OXTRm was altered by the events and exposures occurring during labor and birth, future questions would address the impact of this change on other oxytocin-dependent outcomes in the transition to mother/parenthood103.
The sample size for this study was limited by the extenuating circumstances of the pandemic. Study staff also noted that some individuals approached about the study who experienced a much more severe postpartum hemorrhage (> 2000 mL) declined participation due to feeling traumatized by the experience and not wanting to engage with additional activities. Further work on postpartum hemorrhage should consider optimal times to engage individuals in research, given the sensitive nature of this complication. Future work should include gene transcription quantification, which was not possible due to the lack of RNA stabilizing medium at the time of Cesarean birth and institutional research operation limitations. In a larger random sample, future investigations should also consider DNAm variability across maternal age as well as differences by parity as maternal outcomes are strongly linked to these variables. Our sample was also fairly homogenous in terms of maternal race; in more diverse samples, we can examine differences in both DNAm and birth care practices by racial/ethnic demographics.
Primary strengths of this study are the definition of postpartum hemorrhage cases and the sample selection process. Cases were identified not from total blood loss alone (which is known to be typically underestimated)94 but from several variables indicating a person was bleeding heavily after birth. Furthermore, by design, the sample is comprised of individuals who experienced heavy bleeding most likely due to ineffective uterine contractility. Most large cohort or retrospective studies will use medical record diagnoses of postpartum hemorrhage or billing codes based on the estimated blood loss alone or discharge data, which may be very inaccurate104. This approach will group together multiple causes of postpartum hemorrhage into a single outcome, which limits specificity of the associations to oxytocin use or other risk factors. Another strength is the evidence showing a correlation between blood and myometrial OXTRm levels. The finding helps address concerns that blood may not be a viable tissue for understanding the methylation patterns of the uterine tissue. Finally, the pyrosequencing was performed in triplicate and multivariable analyses accounted for variability between replicates. These procedures help ensure that measured OXTRm was as reliable and valid as possible.
We present evidence that higher OXTRm is associated with greater total oxytocin requirements for labor and postpartum management as well as a significant interaction of OXTRm and oxytocin use on higher blood loss due to uterine atony. Together with the findings that maternal blood OXTRm is well correlated to uterine OXTRm, this study presents interesting findings that contribute to our understanding of variability in maternal health outcomes during parturition.
Supplementary information
Acknowledgements
We are grateful for the contributions of the participants to this research following a potentially life-threatening birth complication and during the sensitive postpartum period. We also acknowledge research support received from Brenna Park-Egan MPH, Kayla Tabari-House, the Oregon Health and Science University Women’s Health Research Unit and Pregnancy Research Team, Placenta Repository Project, especially Kierstyn Tuel, Dr. Leonardo Pereira and Dr. Monica Rincon as well as Dr. Natasha Singh and Professor Mark Johnson (Imperial College, University of London). We are grateful for the mentorship and peer support of the Oregon Building Interdisciplinary Research Careers in Women’s Health fellowship throughout this investigation and the ongoing mentorship of Dr. C. Sue Carter. Dr. Erickson was supported during the course of this work in part, by National Institutes of Health grants K12HD043488, K99NR019596-01, CTRC UL1TR002369 as well as a Sigma Theta Tau International Beta-Psi Naomi Ballard Research Award. Dr. Connelly and Dr. Krol’s work was supported in part by National Institutes of Health grant R01HD098117. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Author contributions
Conceptualization: E.N.E., L.M., J.J.C. Methodology: E.N.E., L.M., J.J.C., K.M.K., J.S.D. Investigation: E.N.E., K.M.K., J.S.D. Visualization: E.N.E., J.S.D. Supervision: E.N.E., L.M., J.J.C. Writing—original draft: E.N.E., K.M.K. Writing—review & editing: E.N.E., J.J.C., L.M., K.M.K., J.S.D.
Data availability
Source data for all figures are available as a single Excel file, Supplementary Data. Data is available upon reasonable request to corresponding author for extensions of this study; however, limited to the participants’ data who consented to use in future research during the informed consent process. 1
Code availability
No custom code was used to generate or process the data described in the manuscript.
Competing interests
The authors declare no competing interests.
Footnotes
Supplementary information
The online version contains supplementary material available at 10.1038/s43856-023-00244-6.
References
Associated Data
Supplementary Materials
Data Availability Statement
Source data for all figures are available as a single Excel file, Supplementary Data. Data is available upon reasonable request to corresponding author for extensions of this study; however, limited to the participants’ data who consented to use in future research during the informed consent process. 1
No custom code was used to generate or process the data described in the manuscript.