BoTA-induced enhancement of angiogenesis in endothelial cells
To investigate the effect of BoTA treatment on the tube formation ability of human endothelial cells, we plated HUVECs onto Matrigel-coated plate and 0.5 IU or 2.0 IU of BoTA was applied. The tube formation as measured by the total number of tube loops and branching points was significantly higher in BoTA-treated HUVECs compared to saline-treated controls (Fig. 1a, 1b). More specifically, the assessment of tube formation of human uterine endothelial cells (HUtMEC) revealed that BoTA induced the angiogenic effect exhibiting dramatically thicker and higher numbers of tube formation in BoTA (0.5 IU)-treated groups compared to control (Fig. 1c, 1d). Interestingly, unlike HUVECs, HUtMECs showed no significant difference between the total numbers of tube formation at 0.5 IU and 2.0 IU of BoTA-treated condition. Furthermore, we implemented a micro-engineered 3-dimensional angiogenesis system to explore the effect of BoTA on the ability of vessel formation and sprouting using HUVECs and endometrial stromal CRL4003 cells. The device consists of 5 primary channels; 2 fluidic microchannels separate 3 hydrogel-laden microchannels from each other to facilitate the supply of fresh media through the device (Fig. 1e). This co-culture system enables paracrine interaction between HUVECs and endometrial stromal layer, mediating the angiogenic morphogenesis. The micro-engineered angiogenesis model was exposed to BoTA at 0.5 IU or 2.0 IU. These analyses revealed that BoTA treatment induced an increasing effect on angiogenesis displaying promoted angiogenic sprouting of HUVECs in BOTA-treated group compared to control (Fig. 1f), which is consistent with data from tube formation assay. This increasing angiogenic effect was further quantified through the analysis of the sprout area (Fig. 1g).
BoTA-induced elevation of endometrial receptivity markers in endometrial cells in vitro
To investigate the effect of BoTA in the endometrial cells in regards to the regulation of endometrial receptivity, BoTA was applied to endometrial epithelial (Ishikawa) and stromal (CRL4003) cells in vitro. Quantitative RT-PCR analyses revealed that BoTA treatment induces significant increases of mRNA expression levels of endometrial receptivity-related genes including ITGB3, HIF1-α, VEGFR2, VIMENTIN, IL-6, and IL-8 in both Ishikawa and CRL4003 cells (Fig. 2a, 2b) suggesting both endometrial epithelial and stromal cells are responsive to BoTA treatment. Particularly, our data show that endometrial epithelial cells are more sensitive to BoTA treatment rather than stromal cells, which are evidenced by the findings showing that Ishikawa cells were more responsive to BoTA with a short-term treatment compared to CRL4003 cells. Even though both Ishikawa and CRL4003 immortalized cells might lack normal physiological features compared to primary endometrial cells, these data might implicate that endometrial epithelial layer is the more effective target of BoTA to induce the faster response of endometrial receptivity. Induction of endometrial receptivity-related genes was validated with immunofluorescence (IF) staining analyses, displaying higher expression of HIF1α, integrin αvβ3, and LIF in BoTA-treated group compared to control (Fig. 2c). Of note, relatively high concentration (10 IU) of BoTA treatment rather led to reduction of all the expressions of endometrial receptivity-related markers (supplementary Fig. 1a, 1b). Furthermore, to evaluate the toxicity or safety of BoTA treatment we performed the cell-based in vitro assays to test if BoTA treatment induces cytotoxicity or aberrant cell proliferation in endometrial epithelial cells prior to further investigation for the clinical applicability of BoTA treatment. These analyses revealed that BoTA treatment induces impact on neither cytotoxicity nor aberrant cell proliferation in endometrial epithelial cells regardless of the concentration ranges between 0.5 IU and 20 IU (supplementary Fig. 2a, 2b). Additionally, BoTA-induced migratory effects of endometrial cells were assessed. Fully confluent cells in the condition of starvation were scratched and their migratory capacity was assessed by measuring the gap distance in between of two sides of cells (supplementary Fig. 2c, 2d). No significant differences were observed in BoTA-treated (0.5 and 2.0 IU) groups compared to controls indicating BoTA has no stimulatory effect on migratory capacity of both endometrial epithelial and stromal cells (supplementary Fig. 2e, 2f). Our observations here demonstrate that BoTA treatment at an appropriate concentration induces positive regulatory effect of endometrial receptivity in both endometrial epithelial and stromal cells in vitro implicating BoTA treatment as a potential therapeutic strategy for infertility patients with specific issues of uterine microenvironment with insufficient receptivity.
Identification of alterations in gene expression induced by intrauterine BoTA infusion in mice
To examine the in vivo effects of intrauterine application of BoTA, BoTA was infused into one side of mouse uterine horns and saline was infused into the other side of horns for the control (Fig. 3a). Mouse uterine tissues were obtained 3 or 8 days after BoTA infusion and total RNA was extracted from the whole endometrial tissues (Fig. 3b). In order to further examine the alterations of global gene expression depending on the intrauterine application of BoTA, RNA-seq data were generated from BoTA-treated versus saline-treated (control) endometrial tissues. Unsupervised hierarchical clustering analyses using a fold change cutoff of 2 and a P-value cutoff of 0.05 identified a total of 197 (88 genes were up-regulated and 109 genes were down-regulated) and 573 (266 genes were up-regulated and 307 genes were down-regulated) differentially expressed genes in BoTA-Day 3 and Day 8 groups, respectively (Fig. 3c and Supplementary Fig. 3a). To perform functional clustering of 573 differentially expressed genes, which were captured from Day 8 BoTA-treated uterine tissues, gene ontology (GO) and pathway analyses were performed by using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tools [34]. Enriched GO terms in each category and pathway including associated-gene counts, P-value, and FDR are shown in Table 1. The P-value and fold enrichment (FE)-value were calculated by Fisher’s exact test and multiple comparisons test, respectively (P < 0.05 and FE > 1.5). A total of 573 differentially expressed genes were classified according to GO terms, including biological process, BP; cellular component, CC; molecular function, MF; and KEGG pathway, KP (Fig. 3d and Table 1). These analyses demonstrate that specific BP categories, including response to steroid hormone (P = 0.023429), positive regulation of growth hormone secretion (P = 0.022572), decidualization (P = 0.002337), positive regulation of vasodilation (P = 0.012491), lung vasculature development (P = 0.001022), and cell adhesion (P = 0.006471) were enriched in BoTA-treated (Day 8) uterine tissues implicating biological processes particularly related to angiogenesis and the early phase of embryo-endometrial interaction in response to hormones. Moreover, enriched terms of the CC category include integral components of plasma membrane (P = 4.39xE− 08), apical plasma membrane (P = 8.52xE− 04), extracellular region (P = 6.03xE− 14), and cell surface (P = 0.004213). Among the enriched CC components, integrin αvβ3-integrin IGF1-IGF1R complex (P = 0.003231) is specifically reported to be strongly correlated to the early stage of embryo implantation [28]. Integrin binding (P = 0.010497), signal transducer activity (P = 0.002136), hormone activity (P = 0.018643), fibroblast growth factor binding (P = 0.023895), and extracellular matrix binding activity (P = 0.029068) were composed of enriched MF terms. Pathway enrichment analyses based on KEGG pathway (KP) analyses were performed using either Chi-square test or Fisher’s exact test, revealing that pathways including cytokine-cytokine interaction (P = 1.30xE− 04) and PI3K-AKT signaling pathway (P = 0.032839), hematopoietic cell lineage (P = 0.005628), focal adhesion (P = 0.037505), and ECM-receptor interaction (P = 0.023927) were enriched in BoTA-treated (Day 8) uterine tissues (Fig. 3d and Table 1). These analyses implicate that intrauterine BoTA treatment might induce the intracellular signaling including ERK and PI3K-AKT pathway to possibly result in promoting angiogenesis [35, 36], which is supported by enrichment of signal transducer activity classified in BP and MF categories. Moreover, intercellular interactions might be enhanced in BoTA-treated endometrium, which was evidenced by high levels of responses including cytokine-cytokine receptor interaction (P = 1.30xE− 04), chemokine-mediated signaling pathway (P = 0.009934), immune response (P = 0.014198), and oxidation-reduction process (P = 0.006635). Furthermore, in-depth clustering analyses revealed that BoTA treatment are strongly associated with the regulation of decidualization, hormonal activity, and interaction between surface adhesion molecules and ECM, which are critically correlated with endometrial receptivity, the process of embryo implantation, and maintenance of pregnancy [37].
Table 1
Whole Gene ontology (GO) and pathway analysis of differentially expressed genes of BoTA-day8 samples compared to control using DAVID tool. The cutoff for significance was set by P-value < 0.05 and Fold Enrichment (FE) > 1.5. biological process (BP), cellular component (CC), and molecular function (MF). KEGG pathway (KP) annotations.
Category
|
Term
|
Gene Count
|
P-value
|
Fold Enrichment
|
Biological process (BP)
|
GO:0070374 ~ positive regulation of
ERK1 and ERK2 cascade
|
15
|
1.77E-04
|
3.3319
|
Biological process (BP)
|
GO:0008285 ~ negative regulation of cell proliferation
|
22
|
3.94E-04
|
2.392489
|
Biological process (BP)
|
GO:0060426 ~ lung vasculature development
|
4
|
0.001022
|
18.55992
|
Biological process (BP)
|
GO:0046697 ~ decidualization
|
5
|
0.002337
|
8.699962
|
Biological process (BP)
|
GO:0007165 ~ signal transduction
|
47
|
0.002526
|
1.563913
|
Biological process (BP)
|
GO:0090131 ~ mesenchyme migration
|
3
|
0.005429
|
25.05589
|
Biological process (BP)
|
GO:0007155 ~ cell adhesion
|
22
|
0.006471
|
1.89426
|
Biological process (BP)
|
GO:0055114 ~ oxidation-reduction process
|
28
|
0.006635
|
1.729696
|
Biological process (BP)
|
GO:0070098 ~ chemokine-mediated signaling pathway
|
6
|
0.009934
|
4.555616
|
Biological process (BP)
|
GO:0045909 ~ positive regulation of vasodilation
|
5
|
0.012491
|
5.494713
|
Biological process (BP)
|
GO:0006955 ~ immune response
|
14
|
0.014198
|
2.149402
|
Biological process (BP)
|
GO:0048146 ~ positive regulation of fibroblast proliferation
|
6
|
0.019461
|
3.854752
|
Biological process (BP)
|
GO:0060124 ~ positive regulation of growth hormone secretion
|
3
|
0.022572
|
12.52794
|
Biological process (BP)
|
GO:0048545 ~ response to steroid hormone
|
4
|
0.023429
|
6.424587
|
Biological process (BP)
|
GO:0030308 ~ negative regulation of cell growth
|
8
|
0.030527
|
2.672628
|
Cellular component (CC)
|
GO:0005576 ~ extracellular region
|
94
|
6.03E-14
|
2.267361
|
Cellular component (CC)
|
GO:0005615 ~ extracellular space
|
90
|
3.68E-16
|
2.530285
|
Cellular component (CC)
|
GO:0005887 ~ integral component of plasma membrane
|
58
|
4.39E-08
|
2.178032
|
Cellular component (CC)
|
GO:0009986 ~ cell surface
|
27
|
0.004213
|
1.815047
|
Cellular component (CC)
|
GO:0016324 ~ apical plasma membrane
|
19
|
8.52E-04
|
2.449371
|
Cellular component (CC)
|
GO:0035867 ~ alphav-beta3 integrin-IGF-1-IGF1R complex
|
3
|
0.003231
|
31.7129
|
Molecular function (MF)
|
GO:0042803 ~ protein homo dimerization activity
|
31
|
0.010672
|
1.609802
|
Molecular function (MF)
|
GO:0004871 ~ signal transducer activity
|
29
|
0.002136
|
1.854542
|
Molecular function (MF)
|
GO:0005178 ~ integrin binding
|
8
|
0.010497
|
3.315154
|
Molecular function (MF)
|
GO:0005179 ~ hormone activity
|
8
|
0.018643
|
2.959959
|
Molecular function (MF)
|
GO:0017134 ~ fibroblast growth factor binding
|
4
|
0.023895
|
6.375297
|
Molecular function (MF)
|
GO:0050840 ~ extracellular matrix binding
|
4
|
0.029068
|
5.919919
|
Molecular function (MF)
|
GO:0001968 ~ fibronectin binding
|
4
|
0.031859
|
5.715783
|
KEGG pathway (KP)
|
mmu04060: Cytokine-cytokine receptor interaction
|
18
|
1.30E-04
|
2.918663
|
KEGG pathway (KP)
|
mmu04151: PI3K-Akt signaling pathway
|
16
|
0.032839
|
1.788709
|
KEGG pathway (KP)
|
mmu04510: Focal adhesion
|
11
|
0.037505
|
2.085207
|
KEGG pathway (KP)
|
mmu04640: Hematopoietic cell lineage
|
8
|
0.005628
|
3.693157
|
KEGG pathway (KP)
|
mmu04512: ECM-receptor interaction
|
7
|
0.023927
|
3.121347
|
Further analyses of enriched GO categories of Day 8 BoTA-treated endometrium compared to saline-treated group, visualized by ClueGO, revealed that BoTA-induced enrichment includes angiogenesis-related pathways such as angiogenesis, cell migration, signaling receptor binding, reproductive structure development, tube development, embryo implantation, and decidualization (Fig. 3e). Among these enrichments, 13 significantly differentially expressed genes including Ccl7, Cyr61, Itgb3, Foxc1, Clec14a, Hif3a, Gpx1, Cd34, Ccbe1, Tgfbi, Lif, Stc1, and Stc2, were specifically associated with GO terms of angiogenesis (GO:0001525) or embryo implantation (GO:0007566) (Table 2).
Table 2
Differentially expressed genes of BoTA-day8 compared to control, which are specifically associated with GO terms of angiogenesis (GO:0001525) or embryo implantation (GO:0007566).
|
Gene symbol
|
BoTA-Day8 /Control
|
Description
|
Fold change
|
P-value
|
Angiogenesis
|
Ccl7
|
7.979
|
0.029
|
chemokine (C-C motif) ligand 7
|
Cyr61
|
4.683
|
0.004
|
cysteine rich protein 61
|
Itgb3
|
2.894
|
0.002
|
integrin beta 3
|
Foxc1
|
2.803
|
0.001
|
forkhead box C1
|
Clec14a
|
2.779
|
0.034
|
C-type lectin domain family 14, member a
|
Hif3a
|
2.153
|
0.049
|
hypoxia inducible factor 3, alpha subunit
|
Gpx1
|
2.125
|
0.026
|
glutathione peroxidase 1
|
Cd34
|
2.047
|
0.038
|
CD34 antigen
|
Ccbe1
|
0.448
|
0.018
|
collagen and calcium binding EGF domains 1
|
Tgfbi
|
0.326
|
0.015
|
transforming growth factor, beta induced
|
Embryo Implantation
|
Lif
|
4.021
|
0.008
|
leukemia inhibitory factor
|
Itgb3
|
2.894
|
0.002
|
integrin beta 3
|
Stc2
|
2.522
|
0.033
|
stanniocalcin 2
|
Stc1
|
0.447
|
0.018
|
stanniocalcin 1
|
Enhanced endometrial blood vessel formation by intrauterine BoTA treatment
In order to validate our RNA-seq data, genes that are classified and strongly related to angiogenesis including Ccl7, Cry61, Itgb3, Clec14a, Hif3a, Gpx1, Cd34, Ccbe1, Tgfbi, Vegfr1, and Tie1 were subjected to quantitative QRT-PCR analyses. The validation revealed that the expression pattern examined by QRT-PCR was concordant with RNA-seq data (Fig. 4a). To further corroborate the evidence for the induction of endometrial angiogenesis with BoTA infusion in mouse uterus, our interrogation of paraffin-embedded BoTA-treated endometrial tissue sections with IF staining of CD31, a surrogate marker for blood vessel formation, revealed that BoTA treatment increased CD31 expression in mouse uterus compared to control group (Fig. 4b, supplementary Fig. 4a), which is consistent with our RNA-seq analyses. In particular, BoTA-infused endometrial tissues exhibited enhanced CD31 levels in the stromal layer especially close to the lining of epithelium. Moreover, the total number and surface area of CD31 staining were further quantified exhibiting significant and gradual increase in BoTA-treated groups (Day 3 & 8) compared to control (Fig. 4c, 4d). As a novel marker for angiogenesis of vascular endothelial cells and progenitor cells with stemness [38], CD34 expression pattern was assessed in paraffin-embedded BoTA-treated endometrial tissue sections revealing more abundantly expressed CD34 in BoTA-treated groups compared to control (Fig. 4e, 4f). These observations indicate that BoTA infusion might have an angiogenic effect to enhance the development of new capillaries from pre-existing blood vessels in the endometrium.
Improved endometrial receptivity by intrauterine BoTA treatment
Among significantly altered embryo implantation-related genes, differential expression pattern of Itgb3 and Lif upon BoTA treatment showed consistency with data from in vitro analyses shown in Fig. 2. In order to further support our RNA-seq data, genes that are classified and strongly related to embryo implantation, including Integrin β3, osteopontin (OPN), LIF, STC1, and STC2, were subjected to Immunohistochemistry or QRT-PCR analyses. Integrin β3 and its ligand OPN proteins were remarkably highly expressed in BoTA-treated endometrial tissue section compared to saline-treated control group (Fig. 5a, 5b). The validation with QRT-PCR revealed concordant pattern of expression of Lif, Stc1, and Stc2 with those from RNA-seq analyses (Fig. 5c-5e), suggesting that intrauterine infusion of BoTA might be an effective method to improve the endometrial receptivity for patients who are suffering from implantation failure with poor receptivity.
BoTA-treated uterine fertility assessment
Angiogenesis is known to participate in a wide range of the process in pregnancy including folliculogenesis, early implantation, placentation, and embryonic development [2, 39]. This led us to examine the correlation of enhanced neovascular activity and increased molecular features of endometrial receptivity induced by BoTA treatment with the outcomes of embryo implantation and pregnancy. For the assessment of BoTA effect on the stability of embryo attachment in vitro, a total of 73 day 5 mouse embryos were transferred (one embryo per well) to confluent Ishikawa cells, which were treated with BoTA (primed; 5 h of exposure to BoTA prior to embryo co-culture, non-primed; commence BoTA exposure at the time of embryo co-culture) and co-cultured for 19 h, 24 h, 28 h or 45 h. The stability of attached mouse embryo was assessed according to the 5-stage standard: score 1, floating; 2, weakly attached but detached after tapping; 3, weakly attached but stuck at the attachment site after tapping; 4, stably attached; and 5, stably attached and showed outgrowth [30, 31]. Up to 24 h of co-culture, embryos transferred onto both control and BoTA-treated cells showed no significant difference in their attachment stability. However, at the late stage of co-culture (24 h-48 h) embryos were more stably attached onto BoTA-primed treated Ishikawa cells compared to controls (Fig. 6a). Of particular note, no significant difference was observed in the stability of attached embryo onto non-primed BoTA-treated Ishikawa cells compared to that of control cells (Supplementary Fig. 5a), suggesting an effective time-point for BoTA treatment to improve the rates of embryo implantation, which should be primed for the embryo transfer onto endometrial epithelial cells. In order to verify the efficacy of BoTA in vivo, BoTA was applied to mouse uterus in a same manner as shown in Fig. 3a. After 8 days of BoTA-infusion, both sides of uterine horns were obtained to expose the embryo implantation sites 16 days after mating (Fig. 6b). These analyses revealed that the total number of implantation sites in BoTA-treated group was significantly higher than controls (P = 0.0121) (Fig. 6c, 6d). Intriguingly, increased numbers and thickness of uterine arteries, which formed the network towards each site of embryo implantation in mouse uterus, were observed at the open abdominal surgery in BoTA-administered side compared to control (Fig. 6e). This might be because of increased endometrial neomicrovessles, which might have been developed and stimulated from pre-existing uterine arteries, induced by BoTA intrauterine treatment. Importantly, no retarded embryos were detected in both groups and the weight of embryos from both groups showed no significant difference (Fig. 6f-6 h). Representative images of the implantation sites in the uterus of mice and embryos obtained from different conditions are shown in Fig. 6c and 6 f.