Gene expression associated with unfavorable vaginal bleeding in women using the etonogestrel subdermal contraceptive implant: a prospective study

To evaluate gene expression associated with unfavorable vaginal bleeding in users of the Etonogestrel (ENG) contraceptive implant. Prospective study involving 100 women who intended to use the ENG implant. Exclusion criteria included abnormal uterine bleeding, inability to attend a 1-year follow-up, and implant removal for reasons unrelated to vaginal bleeding or loss of follow-up. We obtained endometrial biopsies before implant placement and assessed the expression of 20 selected genes. Users maintained a uterine bleeding diary for 12 months post-implant placement. For statistical analysis, we categorized women into those with or without favorable vaginal bleeding at 3 and 12 months. Women with lower CXCL1 expression had a 6.8-fold increased risk of unfavorable vaginal bleeding at 3 months (OR 6.8, 95% CI 2.21–20.79, p < 0.001), while those with higher BCL6 and BMP6 expression had 6- and 5.1-fold increased risks, respectively. By the 12-month follow-up, women with lower CXCL1 expression had a 5.37-fold increased risk of unfavorable vaginal bleeding (OR 5.37, 95% CI 1.63–17.73, p = 0.006). Women with CXCL1 expression < 0.0675, BCL6 > 0.65, and BMP6 > 3.4 had a higher likelihood of experiencing unfavorable vaginal bleeding at 3 months, and CXCL1 < 0.158 at 12 months. Users of ENG contraceptive implants with elevated BCL6 and BMP6 expression exhibited a higher risk of breakthrough bleeding at the 3-month follow-up. Conversely, reduced CXCL1 expression was associated with an elevated risk of bleeding at both the 3 and 12-month follow-ups.


Participants
We included a convenience sample of 100 women aged between 18 and 45 years who expressed a desire to use the ENG subdermal implant (Implanon NXT ® ; Organon, Oss, The Netherlands) for contraception and displayed no contraindications to the method 17 .Women with a history of abnormal uterine bleeding (AUB) 18 , including endometriosis, endometrial polyps, uterine fibroids (leiomyomas), and adenomyosis were excluded, women in use of hormonal method of contraception, as well as those unable to attend a 1-year follow-up.Participants who had the implant removed for reasons unrelated to unfavorable vaginal bleeding within one year of placement were also discontinued from the study.

Procedures
All participants underwent a comprehensive gynecological examination immediately before ENG implant insertion, during which an endometrial biopsy was obtained using exclusively a Pipelle of Cornier ® in the first menstrual phase (until 7 days of the cycle).These specimens were then stored at − 80 °C until analysis.Following the gynecological examination, the ENG implant was placed.Participants were instructed to maintain a vaginal bleeding diary throughout 12 months after device placement.Instances of vaginal bleeding were characterized as any bloody vaginal discharge necessitating the use of pads or tampons, while spotting was described as any bloody vaginal discharge that did not require protective measures 19 .Regular telephone contacts occurred every three months, with a comprehensive face-to-face assessment at the 12-month mark to collect and assess all menstrual diaries.

RNA extraction
For the targeted genes, RT-qPCR employed pre-designed and validated hydrolysis probes obtained from Thermo Fisher Scientific™.Serial dilutions were used to create a standard curve to determine the optimal RNA quantity for each qPCR reaction, which was determined to be 25 ng.Three reference genes were utilized: ACTB, GAPDH, and PRDM4.The first two genes, widely recommended by the manufacturer as qPCR controls, were supplemented by PRDM4, chosen based on its stability in endometrial samples after comprehensive testing 29 .
The reverse transcription process employed the SuperScript™ IV VILO™ Master Mix with ezDNase™ Enzyme (Cat# 11766050, Thermo Fisher Scientific).Each sample consisted of 2.5 µg of the RNA reaction mix, strictly following the manufacturer's guidelines.This process included DNAse enzyme treatment at 37 °C for 2 min to ensure the absence of DNA contamination in reverse transcription.The cDNA was generated through steps including primer annealing, reverse transcription, and enzyme inactivation.The resulting cDNA was promptly utilized in qPCR reactions or stored at − 80 °C.
For qPCR, the TaqMan™ Fast Advanced Master Mix (Cat# 4444965, Thermo Fisher Scientific) was employed, with each 20-µl reaction containing 25 ng cDNA.Following the manufacturer's hydrolysis probe protocol, samples were run in duplicate on 96-well plates using the 7500 Real-Time PCR System (Cat# 4377354, Thermo Fisher Scientific).The ensuing analysis utilized the Relative Quantification application available on the Thermo Fisher ConnectTM website (version 1.1).Replicates with more than a 0.5 cycle variation in Cqs were excluded, and interplate variation was mitigated by employing a consistent known sample for calibration.

Statistical analysis
Participants were categorized into groups experiencing either favorable or unfavorable vaginal bleeding at 3 and 12 months.The concept of favorable bleeding encompassed reduced menstrual flow (amenorrhea, infrequent bleeding), while unfavorable bleeding was characterized by increased or unpredictable flow (defined as frequent, prolonged, and irregular bleeding) as delineated by Hankel and Belsey.To explore factors linked to unfavorable bleeding at 3 and 12 months, both simple and multiple logistic regression analyses were conducted, with variable selection employing a stepwise criterion.Statistical significance was set at 5%.Mann-Whitney analysis was utilized to associate genes with bleeding patterns at 3 and 12 months.In instances where genes demonstrated a trend (significant level up to 1%) or statistically significant difference, ROC curves were applied to determine predictive gene expression thresholds.
We calculated the sample power considering proportion values in each group setting the significance level alpha at 5% (Type I error), and the size of groups in the current sample (n = 59 favorable and n = 37 unfavorable in 3 months of follow up and n = 50 favorable and n = 40 unfavorable in 12 months follow up).Calculation of the sample size considering an 80% power and a 5% significance level.

Results
Figure 1 presents details regarding participant selection.A total of 100 ENG implants were inserted.During the follow-up period, contact was lost with 4 participants, resulting in a total sample of 96 women available at the 3rd month of analysis.Between the 3rd and 12th months, an additional 4 participants were lost to follow-up: 2 removed the device for personal reasons and were consequently excluded, while the other 2 did not return for their appointments.This led to a final count of 92 women available for the 12-month follow-up analysis.Within the 3-to 12-month period, a total of 10 women chose to remove the device due to unfavorable bleeding.
At the 3-month follow-up, uterine bleeding profiles were categorized as favorable in 59 women (61.5%) and unfavorable in 37 women (38.5%).The mean age of the participants was 27.8 (± 6.4) years, with no significant statistical distinctions observed between the groups experiencing favorable and unfavorable bleeding (Table 1).There were no disparities noted in terms of ethnicity, years of schooling, BMI (kg/m 2 ), parity, age at menarche, menstrual period duration, or prior cycle length before the study initiation.
After the results, the sample power was calculated, with the genes CXCL1 < 0.0675 having a power of 0.695, BMP6 > 3.4 with a power of 0.760, and BCL2 > 0.65 with a power of 0.746 to distinguish between women with favorable and unfavorable bleeding at the 3-month follow-up.For the gene CXCL2, the power was 0.963 to differentiate women with favorable and unfavorable bleeding at the 12-month follow-up.

Gene expression analysis
Regarding gene expression, BCL6 exhibited higher expression (p = 0.035), whereas CXCL1 showed lower expression among women experiencing unfavorable bleeding during the third-month follow-up (p = 0.05).Additionally, BMP6 displayed a trend toward elevated expression among women with unfavorable bleeding (p = 0.063) (Table 2).Multiple regression analyses revealed that women with decreased CXCL1 expression had a 6.8-fold higher likelihood of experiencing unfavorable vaginal bleeding within 3 months (OR 6.8, 95% CI 2.21-20.79,p < 0.001).Similarly, women exhibiting higher expression of BCL6 and BMP6 faced sixfold and 5.1-fold increased odds of unfavorable vaginal bleeding, respectively (Table 3).

Receiver operating characteristic (ROC) analysis
Figure 2 presents the results of Receiver Operating Characteristic (ROC) curve analysis aimed at establishing cutoff points for significant gene values as predictors of vaginal bleeding.A BCL6 value exceeding 0.65 demonstrated a substantial area under the curve, indicating a heightened likelihood of vaginal bleeding at 3 months.Conversely, CXCL1 values below 0.0675 were associated with an increased probability of vaginal bleeding at the 12-month mark (p = 0.049).

Discussion
Our study uncovered a significant association between gene expression patterns and the occurrence of breakthrough bleeding in women using the ENG implant.Breakthrough bleeding, a common concern in contraceptive users, has been inadequately understood in terms of its underlying mechanisms.Our research aimed to shed light on this phenomenon by investigating the expression of specific genes related to endometrial stability and bleeding patterns 25 .
At the 3-month follow-up, we observed that ENG-implant users with elevated BCL6 gene expression had a substantially increased likelihood of experiencing breakthrough bleeding.The BCL6 gene, located on chromosome 3q27.3,plays a crucial role in B cell maturation and has a profound impact on the immune system 30 .Notably, its upregulation has been associated with conditions such as endometriosis, infertility, and pre-eclampsia, www.nature.com/scientificreports/which involve vascular complications 31,32 .Therefore, heightened BCL6 expression may provide insights into the bleeding experienced by ENG-implant users.
In addition to BCL6, the levels two other genes, CXCL1 (with expression < 0.0675, p < 0.001) and BMP6 (with expression > 3.4, p = 0.006), showed potential connections to unfavorable bleeding patterns at the 3-month follow-up.A similar trend was observed for CXCL1 expression (< 0.158, p = 0.006) at the 12-month follow-up.The BMP6 gene, situated on chromosome 6p24.3,encodes a growth factor that binds to TGF-β receptors.Research emphasizes its importance in endometriosis, infertility, and its influence on other genes like GDNF, impacting progesterone availability, a known factor in breakthrough bleeding 33,34 .BMP6 expression also correlates with dysmenorrhea in young women, reflecting responses to menstrual cycle-related inflammation processes 35 .
Furthermore, the CXCL1 gene exerts a profound influence on inflammation, immune responses, and tumor progression 36 .It is closely linked with vascular endothelial growth factors (VEGFs) and mitogen-activated protein kinases (MAPKs) in decidual angiogenesis and arteriogenesis processes 35,37 .Elevated VEGF levels, associated with microvascular density in the endometrium, contribute to bleeding in progestogen users 38 .A previous study demonstrated that downregulation of CXCL1 can lead to adenomyosis development due to its interaction with STAT3 39 .These findings underscore the significance of BMP6 and CXCL1 expressions in uterine tissue pathology, potentially contributing to unfavorable bleeding associated with the contraceptive method under evaluation.
Our study also revealed that the simultaneous involvement of CXCL1, BCL6, and BMP6 genes amplified the likelihood of experiencing breakthrough bleeding.Furthermore, our findings identified specific threshold values: a BCL6 value exceeding 0.65 at the 3-month follow-up and CXCL1 values below 0.0675 at 12 months correlated with breakthrough bleeding occurrences.These findings, to the best of our knowledge, lack direct comparisons in existing literature.
A notable strength of our study lies in addressing the dearth of literature exploring the connection between genes and breakthrough bleeding.Enhancing our understanding of the physiopathology behind this adverse effect linked to contraceptive methods containing progestogens holds promise for improved guidance for users and the development of targeted treatment strategies.However, a limitation of our study is the absence of endometrial samples after 12 months of use due to endometrial atrophy resulting from the method.While clinical signs of endometritis were not observed, confirmatory tests were not conducted.Future studies focusing on cellular processes and integrating different metabolic pathways, such as gene expression and tissue metabolomics, will be crucial for a deeper understanding of the physiopathology of breakthrough bleeding.Future studies with suppression of BCL6 with Gonadotropin-Releasing Hormone analog for 2 months prior to ENG implant insertion and reviewing the percentage of unfavorable bleeding to support your findings can also be interesting.www.nature.com/scientificreports/ In conclusion, our study highlights the association between heightened BCL6 and BMP6 gene expressions and the increased likelihood of unfavorable vaginal at the 3-month follow-up in ENG-implant users.Reduced CXCL1 expression is also linked to bleeding occurrences, with significantly higher odds at both 3 and 12 months.These findings contribute to our understanding of the mechanisms behind breakthrough bleeding in ENG-implant users and have implications for further research and clinical management.Table 2. Gene expression according to favorable bleeding and follow up of ENG implant insertion.Genes: BCL6 transcription repressor, BCL6; bone morphogenetic protein 6, BMP6; complement C3, C3; C-C motif chemokine ligand 2, CCL2; C-C motif chemokine ligand 3, CCL3; C-C motif chemokine ligand 4, CCL4 ; C-C motif chemokine receptor 1, CCR1; CD40 molecule, CD40; C-X-C motif chemokine ligand 1, CXCL1; C-X-C motif chemokine ligand 9, CXCL9; C-X-C motif chemokine ligand 10, CXCL10; C-X-C motif chemokine ligand 12, CXCL12; interleukin 15, IL15; interleukin 17A, IL17A; matrix metallopeptidase 2, MMP2; matrix metallopeptidase 19, MMP19; spleen associated tyrosine kinase, SYK; TIMP metallopeptidase inhibitor 1, TIMP1; TIMP metallopeptidase inhibitor 2, TIMP2; TNF receptor superfamily member 11b, TNFRSF11B; actin beta, ACTB; glyceraldehyde-3-phosphate dehydrogenase, GAPDH; PR/SET domain 4, PRDM4.

of follow up before 12 Figure 1 .
Figure 1.Flow chart of the women included in the study.

Table 1 .
Characteristics of the women according to menstrual bleeding and ENG subdermal insertion time.BMI body mass index.*Exato de Fischer; **X 2 ; ***Mann-Whitney test; # 4 loss of follow up before 3 months; Total (

Table 3 .
Gene expression and unfavorable vaginal bleeding regression analysis (3 months follow up).