Next Article in Journal
Stability and Robustness of Unbalanced Genetic Toggle Switches in the Presence of Scarce Resources
Next Article in Special Issue
Potential Biomarkers of Metastasizing Paragangliomas and Pheochromocytomas
Previous Article in Journal
Expression of Potential Targets for Cell-Based Therapies on Melanoma Cells
Previous Article in Special Issue
De Novo Transcriptome Profiling of Brain Tissue from the Annual Killifish Nothobranchius guentheri
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clusterin and Its Potential Regulatory microRNAs as a Part of Secretome for the Diagnosis of Abnormally Invasive Placenta: Accreta, Increta, and Percreta Cases

by
Angelika V. Timofeeva
1,*,
Ivan S. Fedorov
1,
Mariya M. Pirogova
1,
Oksana N. Vasilchenko
1,
Vitaliy V. Chagovets
1,
Larisa S. Ezhova
1,
Tatiana M. Zabelina
1,
Roman G. Shmakov
1 and
Gennadiy T. Sukhikh
1,2
1
Kulakov National Medical Research Center of Obstetrics, Gynecology, and Perinatology, Ministry of Health of Russia, Ac. Oparina 4, 117997 Moscow, Russia
2
Department of Obstetrics, Gynecology, Perinatology and Reproductology, First Moscow State Medical University Named after I.M. Sechenov, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Life 2021, 11(4), 270; https://doi.org/10.3390/life11040270
Submission received: 30 December 2020 / Revised: 17 March 2021 / Accepted: 20 March 2021 / Published: 24 March 2021
(This article belongs to the Special Issue Genomic and Transcriptomic Alterations in Cancer and Aging)

Abstract

:
Magnetic resonance imaging (MRI) and ultrasound methods used for the diagnosis of an abnormally invasive placenta (AIP) have a wide range of sensitivity (Se, 33–93%) and specificity (Sp, 71–100%) levels, which results in a high risk of unfavorable maternal and perinatal outcomes. The relevance of optimizing the diagnosis of AIP is beyond doubt. Given the epigenetic nature of trophoblast invasion, we aimed to quantitate microRNAs and proteins of their target genes that are potentially associated with AIP in blood plasma samples from 64 pregnant women at gestation weeks 30–34 by reverse transcription coupled with polymerase chain reaction (RT-PCR) and Western blotting, respectively. Statistically significant increases in the expression levels of hsa-miR-17-5p, hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, and hsa-miR-320a-3p were revealed in the groups of women with AIP (accreta, increta, percreta) relative to the group of women with scars on the uterus or to the group with placenta previa. Opposite changes in the expression level of “gene–target protein/miRNA” pairs were found for the α-subunit of the clusterin secretory form and any of the hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, hsa-miR-320a-3p, and hsa-miR-17-5p in all cases of AIP. The developed logistic regression models to diagnose AIP cases of various severity gave Se values of 88.8–100% and Sp values of 91.6–100% using a combination of hsa-miR-21-5p, hsa-miR-92a-3p, hsa-miR-320a-3p, or clusterin levels.

1. Introduction

Abnormal uterine scars after a previous cesarean section, especially in combination with the occurrence of placenta previa, are thought to serve as a basis for defective decidualization and deep trophoblast invasion [1]. Many research groups are working to understand the molecular biological mechanisms related to abnormally invasive placenta (AIP) pathogenesis and associate it with the formation of (i) a placenta with a proangiogenic phenotype (decreased expression level of antiangiogenic sFlt-1 and increased levels of angiogenic growth factors (vascular endothelial growth factor, angiopoietin-2) [2]; (ii) increased proliferation of extravillous trophoblast cells, supported by both syncytiotrophoblast (increased expression levels of epidermal growth factor receptor and c-erbB-2 oncogene) and impaired regulation of apoptosis in the extravillous trophoblast cells themselves (decreased level of INSL4 expression, inducing apoptosis) [3,4]; (iii) overexpression of matrix metalloproteinases in the placenta (MMP2, MMP9), which ensures the migration of trophoblast cells by reorganizing the extracellular matrix [5]; (iv) increased expression of pregnancy-associated plasma protein A in the syncytiotrophoblast, providing an increased local concentration of insulin-like growth factor, which controls the uptake and transport of glucose and amino acids in trophoblast cells and their invasion into the decidual tissue [6,7,8]; and (v) an abnormally aggressive epithelial–mesenchymal transition in the extravillous cells of the trophoblast, which ensures their high invasiveness and does not stop by the beginning of the second trimester but proceeds throughout pregnancy [9].
Some of the main regulators of biological processes, including those mentioned above, are small noncoding RNAs, in particular, miRNAs [10,11,12]. Regulation by miRNA molecules is based on RNA interference and occurs at the transcriptional and post-transcriptional levels, namely (i) as a part of the RNA-induced silencing complex, microRNA recognizes the sequence of the target RNA, and the endonuclease activity of Argonaute proteins leads to the cleavage of the target RNA, and (ii) protein synthesis is blocked by the interaction of the RISC complex with translation initiation factors. miRNAs may participate in the epigenetic silencing of genes by DNA methylation, either indirectly through binding to mRNAs encoding transcription factors or by directly affecting the mRNAs coding for DNA methyltransferases [13]. In turn, the activity of promoters of some miRNAs depends on the expression level of DNA methyltransferases, which form regulatory feedback loops [14]. According to the base miRWalk, miRanda, RNA22, and Targetscan, one miRNA molecule regulates the expression level of tens to hundreds of genes, and at the same time, one gene is a target for many miRNAs. Nevertheless, as regulator of mRNA, miRNA itself can be regulated by pseudogene transcripts, circRNAs, viral RNAs, and lncRNAs through sequestration mechanism [15].
In light of the above, it is clear that the pathogenesis of AIP is complex. Therefore, to diagnose this condition, it is sufficient to identify molecules that reflect components of pathogenesis and are secreted into the biological fluids of the body.
In the present study, we focused on the miRNAs that potentially regulate the epithelial–mesenchymal transition (EMT), which determines the ability of extravillous cytotrophoblast cells to invade. The EMT program is triggered by growth factors synthesized in stromal cells and the activation of various signaling pathways: transforming growth factor beta/bone morphogenetic protein, WNT/β-catenin, Hedgehog, Notch, and phosphatidylinositol 3-kinase /AKT [16]. Most prominent among these cells are macrophages and activated resident fibroblasts that accumulate in chronically inflamed tissues as a result of injury and release growth factors, chemokines, MMP-2, MMP-3, and MMP-9 [17]. This may be the case in chronic endometritis after repeated curettage of the uterine cavity, leading to morphological disorders of the decidual layer of the uterus and being a possible cause of AIP [18]. The result of the activation of these signaling pathways is the expression of the transcriptional repressors of epithelial genes, in particular E-cadherin, and as a consequence, the cell molecular phenotype changes from epithelial to mesenchymal, providing restructuring of the actin cytoskeleton, decreased intercellular adhesion contacts, and increased cell proliferation and migration. In addition, the degree of proteolysis of E-cadherin, which determines the formation of its soluble form, affects the invasiveness of cells, in particular, trophoblast cells. It has been shown that with pathological invasion of the placenta, the process of E-cadherin proteolysis is disrupted, which is probably an important molecular mechanism for controlling the invasiveness of the trophoblast during placenta accreta [19]. Thus, assessment of the soluble form of E-cadherin in the blood plasma of pregnant women could be a useful way to evaluate the risk of AIP.
Another key molecule related to the regulation of cell phenotype transformation is clusterin [20,21]. Clusterin is constitutively expressed in almost all tissues and is present in all organismal body fluids, albeit at varying levels dependent on cell type and physiological status [21]. From the initial protein precursor, different isoforms with various functions can generate: a nuclear isoform (49 kDa), the cytosol and mitochondria isoform (53 kDa), as well as glycosylated proteins that are cleavaged in the endoplasmic reticulum/Golgi apparatus and secreted in the form of α-β heterodimer (75−80 kDa) [22,23,24]. A change in the balance between the isoforms of clusterin is tightly controlled in the cell during different steps of tumor progression associated with changes in apoptotic, proliferative, and invasive processes [25,26,27].
Clusterin integrates multiple signaling pathways involved in tumorigenesis and the progression of various malignancies, among which are those involved in epithelial mesenchymal transition such as IGF-1, TGFβ, VEGF, AKT/MAPK, and Wnt signaling [21]. The expression regulation of clusterin different forms is under control of the tumor–microenvironment interactions providing local cancer growth, invasion, and metastases [28]. Similar mechanisms for controlling the level of clusterin expression can take place during the invasion of embryonic trophoblast cells into the decidual layer of the endometrium. In normal pregnancy, the complex embryo–endometrium crosstalk involving numerous cytokines, growth factors, receptors, and mediators, is necessary [29] to ensure a strictly defined degree of invasion of the trophoblast cells into the decidual layer, without affecting the deeper layers of the endometrium and the adjacent myometrium. Most likely, morphological changes in the decidua due to chronic inflammation in the uterus lead to pathological invasiveness of the placenta after embryo implantation as a result of impaired expression level regulation of clusterin.
The specificity of clusterin expression for syncytiotrophoblast cells and the endothelium of chorionic villi has been proven, and increased expression has been found in preeclampsia [30,31], the key pathogenetic mechanism of which is the dysregulation of the phenotypic transformation of trophoblast cells and their reduced invasive ability. However, there are still no data on changes in the expression and role of clusterin in the development of AIP.
In connection with the above, we aimed to quantify soluble E-cadherin and the secreted form of clusterin as well as the role of miRNAs as potential regulatory molecules in the expression levels of these compounds in the blood plasma of pregnant women with AIP of various levels of severity to assess the diagnostic value of each of them for use in clinical practice.

2. Materials and Methods

2.1. Patients

In total, 64 pregnant women aged between 18 and 45 years with cesarean section indications were enrolled in the study and comprised the five groups (Table 1). The inclusion criteria for the study were singleton pregnancy, scar on the uterus after previous cesarean section or myomectomy, placenta previa, and pathological invasion of the placenta (placenta accreta, increta, percreta). The exclusion criteria from the study were preeclampsia and acute phase or exacerbation of infectious and inflammatory diseases. Pathological invasion of the placenta was diagnosed by ultrasound and magnetic resonance imaging (MRI) studies carried out at 30–34 weeks of pregnancy, during which time peripheral blood was sampled. The pregnant women delivered their babies by cesarean section, and the diagnosis of pathological invasion of the placenta was confirmed by histological analysis of the material of the resected uterine wall with the area of pathological invasion of the placenta. Written informed consent was obtained from each patient, and the study was approved by the ethics committee of the National Medical Research Center for Obstetrics, Gynecology, and Perinatology, named after Academician V.I. Kulakov of Ministry of Healthcare of the Russian Federation (protocol No 8, approval date: 31 October 2019).

2.2. RNA Isolation from Peripheral Blood Plasma

Venous blood samples from pregnant women were collected into S-MONOVETTE tubes containing EDTA KE (Sarstedt AG&Co., Nümbrecht, Germany, cat. No 04.1915.100), centrifuged for 20 min at 300× g (4 °C) followed by plasma collection and recentrifugation for 10 min at 14,500× g. RNA was extracted from 200 µL of blood plasma using an miRNeasy Serum/Plasma Kit (Qiagen, Hilden, Germany, cat. No 217184).

2.3. RNA Isolation from the Placental Tissue

Placental tissue samples were collected for study no later than 10 min after delivery. A 5 mm-thick tissue slice was collected by passing only through the maternal surface of a placenta free of myometrium directly in the area of pathological invasion (P-area) and outside this area (N-area), and this was immediately frozen in liquid nitrogen for subsequent storage at −80 °C. Total RNA was extracted from 20–40 mg of placental tissue using an miRNeasy Micro Kit (Qiagen, Hilden, Germany, catalog No. 217084), followed by an RNeasy MinElute Cleanup Kit (Qiagen, Hilden, Germany, catalog No. 74204). The RNA concentration was measured using a Qubit fluorometer 3.0 (Life Technologies, Singapore, cat. Q33216). The sample quality of the total RNA was examined on an Agilent Bioanalyzer 2100 (Agilent, Waldbronn, Germany, cat. No G2939A) using an RNA 6000 Nano Kit (Agilent Technologies, Vilnius, Lithuania, cat. No. 5067-1511). Total RNA samples with an RNA integrity number (RIN) of at least 8 were used for further study.

2.4. Reverse Transcription and Quantitative Real-Time PCR

Seven microliters from 14 µL of total RNA column eluate (miRNeasy Serum/Plasma Kit, Qiagen, Hilden, Germany, cat. No 217184) extracted from 200 µL of blood plasma was converted into cDNA in a reaction mixture (20 µL) containing 1x Hispec buffer, 1x Nucleics mix, and miScript RT in accordance with the miScript® II RT Kit protocol (Qiagen, Germany, cat. No 218161); then, the sample volume was adjusted with deionized water to 200 µL. The synthesized cDNA (2 µL) was used as a template for real-time PCR using a forward primer specific for the studied RNA (Table 2) and the miScript SYBR Green PCR Kit (Qiagen, Germany, cat. No 218075). The following PCR conditions were used: (1) 15 min at 95 °C and (2) 40 cycles at 94 °C for 15 s, an optimized annealing temperature (52–60 °C) for 30 s and 70 °C at 30 s in a StepOnePlusTM thermocycler (Applied Biosystems, Waltham, MA, USA, cat. No 4376600). The relative expression of miRNA in the blood plasma sample was determined by the ∆Ct method using hsa-miR-382-5p (MIMAT0000737, miRBase, available online: http://www.mirbase.org/, accessed on 24 March 2021) as the reference RNA.

2.5. miRNA Deep Sequencing

cDNA libraries were synthesized using 500 ng of total RNA from the placenta using a NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (Set11 and Set2, New England Biolab®, Frankfurt am Main, Germany, cat. No E7300S, E7580S), amplified for 14 and 18 PCR cycles, respectively, and sequenced on a NextSeq 500 platform (Illumina, San Diego, CA, USA, cat. No SY-415-1001). The adapters were removed with Cutadapt. All trimmed reads shorter than 16 bp and longer than 30 bp were filtered, and only reads with a mean quality higher than 15 were retained. The remaining reads were mapped to the GRCh38.p15 human genome and miRBase v21 with the bowtie aligner [32]. Aligned reads were counted with the featureCount tool from the Subread package [33] and with the fracOverlap 0.9 option, so the whole read was forced to have a 90% intersection with sncRNA features. Differential expression analysis of the sncRNA count data was performed with the DESeq2 package [34].

2.6. Western Blotting

Eleven micrograms of the total protein of the plasma sample, the concentration of which was determined by the Biuret method, was (i) denatured at 65 °C for 5 min in a buffer containing 50 mM tris(hydroxymethyl)aminomethane (Sigma-Aldrich, St. Louis, MO, USA, cat. No T4661) hydrochloride (Tris-HCl) at a pH of 6.8, 1% sodium dodecyl sulfate (SDS) (VWR Life Science AMRESCO, Framingham, MA, USA, cat. No Am-O227-0.1), and 10% glycerol (AppliChem, Darmstadt, Germany, A4443), 238 mM 2-mercaptoethanol (VWR Life Science AMRESCO, Framingham, MA USA, cat. No Am-O482-0.1); (ii) separated in 10% SDS-PAGE in Tris-tricine buffer (100 mM Tris, 100 mM tricine (Sigma-Aldrich, St. Louis, MO, USA, cat. No T0377), 0.1% SDS), as recommended by H. Schägger [35]; (iii) transferred to a PVDF membrane (0.45 µm, Immobilon, Merck Millipore Ltd., Germany, catalog No. IPVH07850) by semi-dry transfer in a Trans-blot® SD semi-dry transfer cell (BioRad, Hercules, CA, USA, cat. No 1703940) using two buffer systems (anode buffer: 40 mM 3-cyclohexylamino-1-propanesulfonic acid (CAPS) (Sigma-Aldrich, USA, cat. No SW18805), 60 mM Tris, pH = 9.6, 15% ethanol; cathode buffer: 40 mM CAPS, 60 mM Tris, pH = 9.6, 0.1% SDS) followed by membrane blocking with 1% non-fat dry milk (M7409, Sigma-Aldrich), 0.1% Tween® 20 (Sigma-Aldrich, USA, cat. No P1379), 50 mM Tris-HCl at a pH of 7.5 and 150 mM NaCl (AppliChem Panreac ITW Companies, Darmstadt, Germany, A1371); and (iv) incubated with primary mouse monoclonal antibodies to E-cadherin at a dilution of 1:400 (5F133, Santa Cruz Biotechnology, Dallas, TX, USA, cat. No sc-71007) or primary mice monoclonal antibodies to the alpha subunit of clusterin at a dilution of 1:400 (B-5, Santa Cruz Biotechnology, Dallas, TX, USA, cat. No sc-5289), followed by incubation with goat anti-mouse secondary polyclonal antibodies conjugated with horseradish peroxidase at a dilution of 1:1000 (R&D Systems, Minneapolis, MN, USA, cat. No HAF007). Protein bands were visualized by enhanced chemiluminescence (SuperSignal™ West Femto Maximum Sensitivity Substrate, ThermoScientific, Rockford, IL, USA, catalog No. 34096) in the ChemiDoc MP gel documentation system (BioRad, Hercules, CA, USA, cat. No 12003154). Densitometric analysis of the band intensity was done using ImageLab™ Software (version 6.0 build 25 Standard Edition, BioRad Laboratories, Inc., Hercules, CA, USA). The size of the analyzed protein was determined in the ImageLabTM program based on the electrophoretic mobility of the PageRulerTM Plus Prestained Protein Ladder (ThermoScientific, Waltham, MA, USA, catalog no. 26619) added to the same SDS-PAGE as the blood plasma samples.

2.7. Statistical Analysis of the Obtained Data

For statistical processing, we used scripts written in R language [33] and RStudio [36]. The correspondence of the analyzed parameters to the normal distribution law was assessed by the Shapiro–Wilk test. When the distribution of data was different from normal, the Mann–Whitney test for paired comparison was used, and data were described as the median (Me) and quartiles Q1 and Q3 in the format Me (Q1; Q3). To reduce the instance of a false positive, Bonferroni correction for multiple testing and Holm–Bonferroni test were used as described in [37]. Bonferroni correction was applied between compared groups and among studied molecules. To identify the relationships among categorical variables, chi-square tests were performed. Since both quantitative and qualitative characteristics were analyzed, a correlation analysis was performed using Spearman’s nonparametric correlation test. The 95% confidence interval for the correlation coefficient was determined using the Fisher transformation. The value of the threshold significance level (p) was taken as equal to 0.05. If the p-value was less than 0.001, then p was indicated in the format p < 0.001.

3. Results

3.1. Histological Analysis of the Placenta

The material from the resected uterine wall obtained after cesarean section was examined by histological analysis to reveal signs of pathological invasion of the placenta. The diagnosis of placenta accrete was based on changes in the content and ratio of normal histological components in the uterine wall, in particular, the absence or thinning of decidual tissue without visualization of placenta tissue infiltration into the myometrium (Figure 1a). Placenta increta (Figure 1b) and placenta percreta (Figure 1c) were characterized by prolapse of the chorionic villi into the walls of the large veins in the myometrium and the presence of chorionic villi among the muscle bundles up to the serous layer in the case of placenta percreta (Figure 1c).

3.2. Selection of miRNAs That Potentially Regulate E-Cadherin and Clusterin

According to four electronic databases—miRWalk, miRanda, RNA22, and Targetscan—hsa-miR-320a-3p, hsa-miR-17-5p, hsa-miR-21-5p, hsa-miR-1323, hsa-miR-25-3p, hsa-miR-138-5p, hsa-miR-34a-5p, and hsa-miR-92a-3p are the common potential regulators of the expression levels of E-cadherin and clusterin; hsa-miR-30a-5p and hsa-miR-30c-5p specifically bind to clusterin mRNA; and hsa-miR-371a-5p and hsa-miR-506-3p potentially regulate E-cadherin. However, it should be noticed that even with the existence of experimentally validated interactions between miRNA and mRNA of the target gene, as in the case of hsa-miR-92a-3p, hsa-miR-25-3p, hsa-miR-138-5p, and E-cadherin, or in the case of hsa-miR-21-5p, hsa-miR-17-5p, and clusterin according to the miRTargetLink Human database, the impact of miRNA on the mRNA of the target gene is not so unambiguous. This is due to the mechanisms of the regulation of miRNA activity, among which is the sequestration mechanism [15], along with post-translational modifications that impact target gene protein activity and function, among which are methylation, acetylation, glycosylation, phosphorylation, ubiquitination, and SUMOylation [38,39,40,41,42,43,44,45].
Thus, in our study, 12 miRNAs were analyzed in 64 blood plasma samples from pregnant women from groups I–V (Table 1) by a quantitative real-time RT-PCR method using a sense primer specific for miRNA (Table 2) and a universal antisense primer from the miScript SYBR Green PCR kit. For data normalization, the use of miRNAs as normalizer seems to be more appropriate than other small noncoding RNAs such as RNU6, since this normalization strategy with the same family of RNA molecules comprises the same methods, such as extraction, reverse transcription, and PCR, and this conclusion is highlighted in the review article of Madadi S. et al. [46]. We found the expression level of hsa-miR-382-5p to be identical in all 64 analyzed samples except one (Figure 2). In this connection, the relative expression level of miRNAs in the blood plasma samples was determined by the ∆Ct method using hsa-miR-382-5p as the reference RNA.

3.3. Quantitative RT-PCR Analysis of the Selected miRNAs

The principal component analysis was used to determine the clustering of patient samples based on the RT-PCR expression data of hsa-miR-17-5p, hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, and hsa-miR-320a-3p. In general, samples were clustered separately based on the disease status of the first principal component (Figure 3): patients with AIP (III, IV, V groups) clustered away from patients without AIP (I, IIa, IIb groups), except for several outliers in some groups which were removed for further analysis.
Specifically, one patient (ID = 9) was removed from group I (n = 10), one patient (ID = 19) was removed from group IIa (n = 9), one patient (ID = 18) was removed from group IIb (n = 8), one patient (ID = 54) was removed from group III (n = 10), four patients (ID = 29, 30, 36, 44, 58) were removed from group IV (n = 20), and none were removed from group V (n = 7). These samples were excluded for particular reasons. For example, according to MRI data of one particular patient (ID = 54), there was uterine wall swelling and deformation of its contour along with expansion of the veins of the placental bed and expansion of the lacunae of the maternal part of the placenta which indirectly indicated placenta accreta, although according to histological analysis, no areas of placenta accreta were found. As for the samples from patients (IDs = 29, 30, 36, 44, 58) from group IV, despite the diagnosis of “placenta increta” identified by MRI and ultrasound methods, a placental hernia was revealed during the cesarean section, the formation of which was not so much due to the pathological invasiveness of the placenta itself but rather the thinning of the myometrium and the divergence of its fibers. That is why miRNAs implicated in the EMT did not differentiate these samples from normal ones. The cohort of patients without outliers mentioned above was reanalyzed, and the results of pairwise comparison of the AIP groups with groups I or II are presented in Table 3, Table 4, Table 5 and Table 6. To visualize the data obtained, a box diagram was plotted (Figure 4).
Five out of 12 analyzed miRNAs (miR-17-5p, miR-21-5p, miR-25-3p, miR-92a-3p, and miR-320a-3p) showed statistically significant decreases in ΔCt median values in at least one of the groups with pathological invasion of the placenta (III–V) relative either to group I (Table 3 and Table 4) or to group II (Table 5 and Table 6). It is important to note that the lower the ΔCt value, the higher the microRNA expression level. There were no statistically significant differences in the expression levels of miR-17-5p, miR-21-5p, miR-25-3p, miR-92a-3p, and miR-320a-3p between groups I and II. However, the presence of placenta previa in group II and the absence of this condition in group I strongly influenced the statistical significance of differences in the expression level of some miRNAs when compared with groups III, IV, or V. For example, when taking into account the Bonferroni correction for multiple testing or applying the Holm–Bonferroni test, changes in the expression level of miR-21-5p and miR-320a-3p were significantly higher in groups III–V in comparison with both group I and group II; in contrast, miR-17-5p and miR-92a-3p were significantly upregulated in all groups with AIP (III–V) only in comparison with group II. As for miR-25-3p, the level of expression was significantly increased only in group V when compared to groups I or II. The variability in the significance of differences in groups III–V depending on the group selected for comparison (I or II) for some miRNAs may be due to the following reasons: i) the number of samples in each of the groups and/or ii) a slight change in the expression levels of these miRNAs in compared groups and/or iii) the presence of placenta previa in groups II–V and its absence in group I. From our point of view, it is more correct to compare each of groups III–V with group II and not with group I, since placenta previa possibly arises due to the characteristics of placental tissue itself, including the expression profile of genes encoding for the proteins and miRNAs that regulate them, and this parameter should be taken into account in all compared groups. Our assumptions are supported by literature data, in particular, by increased mRNA expression levels of HMGB1 and VEGF [47] and by decreased expression of β-catenin [48] in the placenta previa group compared to that in the normal group.
To assess the relationship between the microRNA expression level and the severity of AIP, the groups of analyzed samples were ranked in the following way: I < II < III < IV < V. The Spearman’s rank correlation method revealed the following statistically significant correlations: miR-17-5p (r = −0.33, p = 0.007), miR-21-5p (r = −0.68, p < 0.001), miR-25-3p (r = −0.36, p = 0.004), miR-92a-3p (r = −0.45, p < 0.001), miR-320a-3p (r = −0.55, p < 0.001). The revealed negative correlation between ΔCt miRNA and the severity of AIP indicated a positive correlation between increased expression levels of miR-17-5p, miR-21-5p, miR-25-3p, miR-92a-3p, and miR-320a-3p in the peripheral blood plasma of pregnant women and the severity of AIP.

3.4. Deep Sequencing of miRNAs from Placenta Tissue Samples of Women in Groups III, IV, and V

To confirm the correlations revealed for the expression levels of miR-17-5p, miR-21-5p, miR-25-3p, miR-92a-3p, and miR-320a-3p with the severity of AIP and to prove the placental origin of these miRNAs, deep sequencing of small noncoding RNAs from the placenta in the area of pathological invasion (P-area) and outside this area (N-area) was performed. Placental samples (N and P) from eight women from groups III, IV, and V were taken for sequencing, and the read numbers of miRNAs in these samples processed by the DESeq2 package [34] are presented in Table S1 for all identified miRNAs, and in Table 7 and Table 8 for miR-17-5p, miR-21-5p, miR-25-3p, miR-92a-3p, and miR-320a-3p. The fold changes in the miRNA expression level in groups IV and V relative to group III and the statistical significance of these differences were calculated. It follows from Table 7 that with increasing AIP severity, the expression levels of miR-17-5p, miR-21-5p, miR-25-3p and miR-92a-3p, miR-320a-3p in the P-areas significantly increased. Increases in the expression levels of these miRNAs in groups IV and V relative to group III were also observed in the N-areas, but these changes were comparable, and in such cases, statistically significant increases in the expression levels in both groups were observed only for miR-25-3p and miR-92a-3p (Table 8).

3.5. Quantitative Analysis of Potential Targets of miRNAs—Participants in the EMT

In order to analyze the expression levels of potential targets of miR-17-5p, miR-21-5p, miR-25-3p, miR-92a-3p, and miR-320a-3p, the contents of the soluble form of E-cadherin and the alpha subunit of the secretory form of clusterin were quantified by Western blotting in the same peripheral blood plasma samples in which the expression levels of miRNAs regulating them were analyzed.

3.5.1. E-Cadherin Quantitation in the Peripheral Blood Plasma of Pregnant Women

The E-cadherin monoclonal antibody used in the present study was raised against the extracellular domain of E-cadherin, shedding into the blood plasma as a result of proteolysis of the full-length form by any of a number of proteases including MMP3, MMP7, MMP9, plasmin, kallikrein, MT-1-MMP, ADAM10, caspases, and calpain, leading to the formation of E-cadherin fragments with different molecular weights [40,41,42,43,44].
Seven blood plasma samples from each of the groups I–V were pooled and analyzed by Western blotting. The same quantity loading of total protein into the wells of the PAAG/SDS and the efficiency of semi-dry transfer of proteins from the gel to the PVDF membrane were controlled by staining the membrane with Ponceau S dye (Figure 5a). It follows from Figure 5b that three soluble fragments of E-cadherin with molecular weights of 87.4, 60.6 and 26.4 kDa were revealed in all studied sample groups. Quantitative analysis of chemiluminescence showed no obvious differences between the analyzed groups in the content of each of the three fragments of soluble E-cadherin (the data are summarized in the table presented in the inset of the Figure 5b). Therefore, a detailed analysis of the content of soluble E-cadherin individually in each sample of the analyzed groups I-V was not carried out.

3.5.2. Quantitation of the Clusterin Secretory form in the Peripheral Blood Plasma of Pregnant Women

There are two isoforms of clusterin synthesized in the cell as a result of alternative splicing: the secretory and nuclear forms. The nuclear form of clusterin is a low-glycosylated protein and is translocated from the cytoplasm to the nucleus. The secretory form of clusterin is a glycosylated protein consisting of α and β chains linked by five disulfide bonds.
Seven blood plasma samples from each of the groups I–V were pooled and analyzed by Western blotting. The same quantity loading of total protein into the wells of the PAAG/SDS and the efficiency of semi-dry transfer of proteins from the gel to the PVDF membrane were controlled by staining the membrane with Ponceau S dye (Figure 6a). It follows from Figure 6b that the level of the α-subunit of the clusterin secretory form with molecular weight of 47.9 kDa was decreased in all samples with AIP (groups III–V) relative to group I (chemiluminescence data are presented in the inset of Figure 6b). Therefore, a detailed analysis of the content of the α-subunit of the clusterin secretory form individually in each sample of the analyzed groups I–V was carried out.
Seven blood plasma samples from each of the groups I–V were taken into analysis. To take into account the efficiency of protein transfer from PAAG/SDS to the PVDF membrane, the binding of primary and secondary antibodies, exposure in the analysis of chemiluminescence, the same sample (number 1 from group I, I-1) was applied to each PAAG/SDS to normalize all studied samples (Figure 6c).
A quantitative assessment of the alpha subunit of the secretory form of clusterin (47.9 kDa) in the peripheral blood plasma of pregnant women revealed a statistically significant decrease in the level of clusterin expression in all analyzed groups of women with AIP (III–V) relative to group I or group II (Figure 7, Table 9).
Using the Spearman rank correlation method, significant positive correlations of the content of clusterin in the peripheral blood plasma with the miRNA ΔCt values were found, namely for miR-21-5p (r = 0.74, p < 0.001), miR-25-3p (r = 0.54, p = 0.001), miR-92a-3p (r = 0.69, p < 0.001), miR-320a-3p (r = 0.64, p < 0.001), and miR-17-5p (r = 0.53, p = 0.001), which indicates that there were inverse correlations between the expression levels of these miRNAs and the concentration of clusterin. At the same time, an inverse correlation was found between the expression level of the α-subunit of clusterin and the severity of AIP (r = −0.79, p < 0.001). Inverse correlations between the expression level of miR-25-3p, miR-92a-3p, miR-320a-3p, miR-17-5p, and clusterin may not be a consequence of the direct interaction of these molecules, while the existence of experimentally validated interactions between hsa-miR-21-5p, hsa-miR-17-5p, and clusterin according to the miRTargetLink Human database have been already proven. This fact is not an obstacle for AIP diagnostics, because for this purpose, it is important to reveal the differential expression level of miR-25-3p, miR-92a-3p, miR-320a-3p, miR-17-5p, and clusterin, rather than the physical interaction of molecules.

3.6. Development of Logistic Regression Models for Each form of AIP

Based on the values of clusterin content and the expression levels of its potential regulators, miR-21-5p, miR-92a-3p, and miR-320a-3p, in the peripheral blood plasma, due to their highly significant correlations with the severity of AIP, logistic regression models were developed for calculating the probability of the presence of placenta accrete, increta, or percreta in a pregnant woman at 32–34 weeks of gestation. Based on the expression levels of miR-21-5p, miR-92a-3p, and miR-320a-3p, logistic regression models were developed for the entire cohort of patients, including pregnant women with unconfirmed MRI/ultrasound diagnoses of AIP during cesarean section and/or histological analysis (Table S2, Sheet 1: patients without AIP, n = 27; patients with placenta accreta, n = 10; patients with placenta increta, n = 20; patients with placenta percreta, n = 7), and for an optimized cohort of patients (Table S2, Sheet 2: patients without AIP, n = 24; patients with placenta accreta, n = 9; patients with placenta increta, n = 14; patients with placenta percreta, n = 7 ). It seemed interesting to find out whether the classification of samples by the presence or absence of AIP, based only on MRI/ultrasound data, affected the parameters of the developed logistic regression models. We found that the diagnostic accuracy of expression profiling of the combination of miR-21-5p, miR-92a-3p, and miR-320a-3p for the identification of placenta accreta was higher in the case of optimized cohort of patients (AUC = 1, p < 0.001; accuracy—100 %, sensitivity—100%, specificity—100%, true positive rate—1; false positive rate—0) compared to that of the case of entire cohort of patients (AUC = 0.9, p < 0.001; accuracy—88.2 %, sensitivity—90%, specificity—95.8%, true positive rate—0.9; false positive rate—0.125). Similarly, the diagnostic accuracy of expression profiling of the combination of miR-21-5p, miR-92a-3p, and miR-320a-3p for the identification of placenta increta was significantly higher in case of optimized cohort of patients (AUC = 0.981, p < 0.001; accuracy—94.9 %, sensitivity—100%, specificity—91.6%, true positive rate—1; false positive rate—0.0833) compared to that in case of entire cohort of patients (AUC = 0.848, p < 0.001; accuracy—80.9 %, sensitivity—75%, specificity—92.5%, true positive rate—0.55; false positive rate—0). Receiver operating characteristic (ROC) curves of the developed logistic regression models in cases of entire and optimized cohorts of patients are presented in Figure 8. The diagnostic accuracy of expression profiling of each of miR-21-5p, miR-92a-3p, miR-320a-3p for the identification of placenta percreta was the same for the two cohorts of patients (Table S1, Sheet 1 and Sheet 2: AUC = 1, p < 0.001; accuracy—100 %, sensitivity—100%, specificity—100%, true positive rate—1; false positive rate—0).
The α-subunit of the clusterin secretory form expression profile was also able to differentiate between patients with AIP and patients without AIP. In particular, the following parameters of the developed models were obtained: AUC = 1 (p < 0.001), accuracy—100 %, sensitivity—100%, specificity—100%, true positive rate—1, false positive rate—0, in the cases of placenta accreta and placenta increta; AUC = 0.987 (p < 0.001), accuracy—94.7%, sensitivity—100%, specificity—100%, true positive rate—1, false positive rate—0.076, in the case of placenta percreta (Table S1, Sheet 3).
The general Formula (1) of these models was as follows:
e = 1 1 + e i k 1 * x 1 k 2 * x 2 ,
where i is the intercept term; k1 and k2 are coefficients for clusterin or each of the miRNAs; and x1, x2, … are the relative expression levels of clusterin or miRNA ΔCt values.
The parameters of the developed models for optimized cohort of patients that reflect the contributions of the expression level of clusterin or miRNAs to diagnose the different forms of AIP are presented in Table 10.
It follows from Table 10 that the different forms of AIP (placenta accreta, increta, percreta) in a pregnant woman at 31–34 weeks of gestation can be diagnosed by the content of the α-subunit of the clusterin secretory form or by the expression level of miR-21-5p, miR-320a-3p, and miR-92a-3p in various combinations in the peripheral blood plasma with a sensitivity of 88.8–100% and a specificity of 91.6–100% at the corresponding cutoff level indicated in this Table 10. This type of analysis significantly increases the diagnostic accuracy of conventional ultrasound and MRI methods, the sensitivity and specificity of which vary in the range of 33–93% Se and 71–100% Sp.

4. Discussion

Research groups are trying to determine the underlying cause of AIP: (i) morphological changes in the decidua due to endometritis and/or antecedent curettage; (ii) incompetent uterine scar after a cesarean section; or (iii) the properties of placental tissue characterized by increased proliferative and reduced apoptotic activity of extravillous trophoblast cells along with abnormally aggressive EMT of these type of cells which does not stop by the end of the first trimester of pregnancy but occurs throughout the entire pregnancy [49,50]. Most likely, in our opinion, chronic inflammation in the uterus leads to pathological invasiveness of the placenta after embryo implantation. This assumption is supported by data showing that in the site of chronic inflammation due to damage, macrophages and activated resident fibroblasts secrete growth factors (TGF-β, PDGF, EGF, FGF-2), chemokines, and matrix metalloproteinases (MMP-2, -3, -9), the result of which is the activation of signaling pathways (receptor tyrosine kinase pathway, Wnt-, Notch-, TGFβ-, Hedgehog, TNFα-pathway) [17] leading to the synthesis of epithelial repressors (Snail, Slug, Smuc, Twist, E12/E47, dEF1, ZEB1, ZEB2), downregulation of epithelial genes (E-cadherin, EpCAM, Claudins, Cytokeratins), and activation of mesenchymal genes (N-cadherin, Vimentin, Fibronectin, MMPs), which causes a highly invasive cell phenotype [16,51]. In this regard, it seems reasonable to diagnose AIP using a noninvasive method of quantitatively evaluating potential participants in the EMT in trophoblast cells, in particular, clusterin and E-cadherin as well as potential regulators of their expression—hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, hsa-miR-320a-3p, and hsa-miR-17-5p—in the peripheral blood plasma of pregnant women.
Mature human E-cadherin is expressed as a glycoprotein (120 kDa) which consists of five extracellular domains, a linker region, a single transmembrane domain, and an intracellular domain. E-cadherin plays an important role in cell–cell contact formation by “homophilic adhesion” through the extracellular domains. E-cadherin can undergo intracellular and extracellular proteolytic cleavage, providing an additional mechanism to transcription repression regulation to reduce cell surface expression [52,53]. The catalytic activity of the E-cadherin-cleaving proteases triggers the extracellular release of soluble E-cadherin fragments of different molecular weights depending on the cleavage site: for MMPs, it is between aa residues 581 and 583 [54] and/or between aa residues 700 and 701 [55]; for γ-secretase, it is between aa residues 731 and 732 [55]; for caspase 3, it is between aa residues 750 and 751 [56]; and for PMN–elastase (leukocyte–elastase), it is between Val-393 and Gly-394 [57]. It is obvious that the expression levels of full-length E-cadherin and its soluble forms shed into the blood are determined by the balance between the biosynthesis and proteolytic cleavage activity of these proteases and may be crucial determinants in the adhesion and migration of cells. Extracellular release of a soluble E-cadherin fragment of about 80 kDa from the cell surface is accompanied by the simultaneous delivery of free β-catenin into the cell cytosol, which then translocates into the cell nucleus where it contributes to the modulation of gene expression [58]. In the present study, in the blood plasma of pregnant women, we found three soluble N-terminal E-cadherin fragments of 87.4, 60.6 and 26.4 kDa, and the ratio of these forms did not change in the blood plasma of pregnant women with any form of AIP in comparison with healthy pregnant women. These data are in good agreement with the absence of significant difference in circulatory soluble E-cadherin levels detected in the serum of women with normal placentation, invasive placentation, and previa alone [19], suggesting that an invasive phenotype of extravillous trophoblasts is not systemically regulated and could be the result of local processing of E-cadherin. It is known that E-cadherin has four utilized N-glycosylation sites, and removal of these sites decreases the functional E-cadherin level at the cell–cell border [59]. In addition, it has been shown that a decrease in the expression level of the membrane-bound form of E-cadherin is not sufficient for the effective regulation of EMT [60,61].
The role of another protein analyzed in this work, clusterin, in the regulation of EMT is discussed mainly in the context of carcinogenesis. Clusterin (apolipoprotein J) is involved in various physiological cellular processes; its function is determined by the structure of the molecule itself [62]. There are two isoforms of clusterin (secretory and nuclear), the transcription of which occurs from different promoter regions of the gene, forming mRNAs that differ in the sequence of the first exon and in the location of the translation initiation site. The secretory form of clusterin is a 76–80 kDa glycosylated protein consisting of α and β chains linked by five disulfide bonds. The nuclear form (55 kDa) of clusterin is a low-glycosylated protein that, after synthesis in the cytoplasm, is translocated into the nucleus. The nuclear form of clusterin performs proapoptotic functions, while the secretory form of clusterin is responsible for the survival of the cell and its proliferation. The prevalence of one form or another in the cell is finely regulated, and it remains a mystery as to what determines their synthesis. For example, it has been shown that, during oncogenic transformation, tumor cell survival is associated with increased expression of the secretory form of clusterin and significant decrease in the level of expression of the nuclear form of clusterin [63].
The expression level of clusterin is regulated at the transcriptional and post-transcriptional levels. The regulators of clusterin transcription are the p53 apoptotic cascade activator [64] and the B-MYB transcription factor, which are involved in the regulation of cell survival, proliferation, and differentiation [65]. In turn, clusterin regulates the activity of the transcription factor NF-κB, which plays important roles in cell viability, cell motility, proliferation, transformation, and inflammation [66].
The specificity of clusterin expression for syncytiotrophoblast cells and chorionic villus endothelium has been proven [30], and its increased expression was also found in preeclampsia both in the placenta [30] and in peripheral blood plasma [31], the key pathogenetic mechanism of which is dysregulation of the phenotypic transformation of trophoblast cells and a reduction in their invasive ability. Until now, there were no data on changes in the expression and role of this protein in the development of AIP. In the present study, for the first time, a significant decrease in the expression level of the secretory form of clusterin was revealed in pregnant women with placenta accreta, increta, and percreta, compared to pregnant women with a scar on the uterus and/or placenta previa.
The expression level of clusterin, like any other protein, can be regulated at the post-transcriptional level by the activity of small noncoding RNAs. In this work, a significant inverse correlation was revealed between the expression level of clusterin and its potential regulators: hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, hsa-miR-320a-3p, and hsa-miR-17-5p. In a review article by Minal Garg [16], the involvement of miR-21-5p, miR-25-3p, miR-92a-3p, miR-320a-3p, and miR-17-5p in the induction of EMT by either suppressing or activating several signaling genes is described in detail. Moreover, downregulation of miR-21-5p has been found to inhibit cancer cell proliferation and induce apoptosis by suppressing the ERK and PI3K/Akt signaling pathways, elevating full-length E-cadherin expression, reducing the expression of the mesenchymal marker vimentin as well as MMP-2 and MMP-9, and acting as an extracellular inducer of MMPs [67]. Increased expression of miR-92a-3p contributes to the EMT and metastasis in colorectal cancer through activation of the Wnt/β-catenin pathway and inhibition of mitochondrial apoptosis [68] as well as in hepatocellular carcinoma and malignant retinoblastoma through the PTEN/Akt pathway [69,70]. Blood plasma exosomes from breast cancer patients and healthy donors were found to differ in the overexpression of miR-25-3p and miR-92a-3p, and this was related to the stimulation of the EMT of nonmalignant breast cells, significant increases in the number of motile and proliferating cells, as well as the formation of capillary-like structures [71].
In addition, changes in the expression levels of these molecules in diseases associated with changes in cell adhesion and migration have been proven. Specifically, there are decreases in the expression levels of miR-92a-3p, miR-17-5p [72], miR-21-5p [73], and miR-320a-3p [74] in the placenta during preeclampsia. Thus, in AIP cases characterized by increased invasion of extravillous trophoblast cells, we identified expression profiles of hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, hsa-miR-320a-3p, and hsa-miR-17-5p and their target gene clusterin that were directly opposite to those seen in preeclampsia, the pathogenesis of which is based on reduced invasion of extravillous trophoblast cells. This underlines the importance of using these molecules for diagnostic purposes. Related to this connection, we constructed logistic regression models to diagnose AIP forms of varying levels of severity (placenta аccreta, increta, percreta) at 32–34 weeks with high specificity and sensitivity based on a quantitative evaluation of the secretory form of clusterin and its regulatory miRNAs—hsa-miR-21-5p, hsa-miR-92a-3p, and hsa-miR-320a-3p—in the peripheral blood plasma of pregnant women. The data obtained in thе present work significantly contribute to improvements in the accuracy of AIP diagnosis by instrumental methods [75] and the analysis of placental and fetal hormones (pregnancy-associated plasma protein A, human chorionic gonadotropin, alpha-fetoprotein) by cell-free placental mRNAs [49] as well as by recently discovered plasma marker miRNAs, in particular, miR-139-3p, miR-196a-5p, miR-518a-3p, and miR-671-3p, which are downregulated in placenta increta/percreta and have the potential to be used for future noninvasive prenatal placenta accreta spectrum screening [76].

5. Conclusions

1. The expression levels of hsa-miR-17-5p, hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, and hsa-miR-320a-3p were significantly increased in the peripheral blood plasma of pregnant women with AIP compared to their concentrations in the blood plasma of pregnant women with placenta previa without AIP. This correlated significantly with the severity of AIP. With an increase in the severity of AIP, the expression levels of these miRNAs increased in the placenta in the area of pathological invasion.
2. The content and ratio of soluble E-cadherin fragments in the blood plasma of pregnant women with any form of AIP did not differ from those in the blood plasma of healthy pregnant women.
3. The level of the secretory form of clusterin was significantly reduced in the peripheral blood plasma of pregnant women with AIP compared to that in the blood plasma of pregnant women without AIP. This significantly inversely correlated with the expression levels of hsa-miR-17-5p, hsa-miR-21-5p, hsa-miR-25-3p, hsa-miR-92a-3p, and hsa-miR-320a-3p and with the severity of AIP.
4. The probability of different forms of AIP (placenta accreta, increta, percreta) occurring in a pregnant woman at 31–34 weeks of gestation could be evaluated by the content of the α-subunit of the clusterin secretory form or by the expression levels of hsa-miR-21-5p, hsa-miR-320a-3p, and hsa-miR-92a-3p in various combinations in peripheral blood plasma from pregnant women with a sensitivity level of 88.8–100% and a specificity level of 91.6–100%, which increases the diagnostic accuracy of conventional ultrasound and MRI methods.
5. It is reasonable to quantify the secretory form of clusterin and the miRNAs that regulate its expression level in the first trimester of pregnancy with the aim of assessing their predictive value in determining AIP development.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life11040270/s1, Table S1: Deep sequencing of miRNAs from placenta tissue samples in N- and P-areas, Table S2: Logistic regression models parameters.

Author Contributions

Conceptualization, A.V.T. and R.G.S.; methodology, A.V.T.; validation, Patient consent was waived due to reason A.V.T., I.S.F., and V.V.C.; investigation, A.V.T. and I.S.F.; clinical database, M.M.P., O.N.V. and T.M.Z.; data curation, A.V.T. and R.G.S.; writing—original draft preparation, A.V.T.; writing—review and editing, A.V.T.; visualization, I.S.F. and L.S.E.; supervision, R.G.S.; project administration, G.T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the State assignment of the Ministry of Healthcare of the Russian Federation [State Registration No. AAAA-A19-119021490133-6].

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of the National Medical Research Center for Obstetrics, Gynecology, and Perinatology, named after Academician V.I. Kulakov of Ministry of Healthcare of the Russian Federation (protocol No 8, approval date: 31 October 2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The NGS data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB43546 (https://www.ebi.ac.uk/ena/browser/view/PRJEB43546), accessed on 9 March 2021.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. Silver, R.M. (Ed.) Placenta Accreta Syndrome; CRC Press: Boca Raton, FL, USA, 2017; ISBN 9781498745963. [Google Scholar]
  2. Tseng, J.J.; Chou, M.M. Differential expression of growth-, angiogenesis- and invasion-related factors in the development of placenta accreta. Taiwan. J. Obstet. Gynecol. 2006, 45, 100–106. [Google Scholar] [CrossRef]
  3. Tseng, J.J.; Hsu, S.L.; Wen, M.C.; Ho, E.S.C.; Chou, M.M. Expression of epidermal growth factor receptor and c-erbB-2 oncoprotein in trophoblast populations of placenta accreta. Am. J. Obstet. Gynecol. 2004, 191, 2106–2113. [Google Scholar] [CrossRef] [PubMed]
  4. Goh, W.; Yamamoto, S.Y.; Thompson, K.S.; Bryant-Greenwood, G.D. Relaxin, its rfeceptor (RXFP1), and insulin-like peptide 4 expression through gestation and in placenta accreta. Reprod. Sci. 2013, 20, 968–980. [Google Scholar] [CrossRef] [PubMed]
  5. Soyama, H.; Miyamoto, M.; Ishibashi, H.; Iwahashi, H.; Matsuura, H.; Kakimoto, S.; Suzuki, R.; Sakamoto, T.; Hada, T.; Takano, M. Placenta previa may acquire invasive nature by factors associated with epithelial-mesenchymal transition and matrix metalloproteinases. J. Obstet. Gynaecol. Res. 2020, 1–8. [Google Scholar] [CrossRef]
  6. Desai, N.; Krantz, D.; Roman, A.; Fleischer, A.; Boulis, S.; Rochelson, B. Elevated first trimester PAPP-A is associated with increased risk of placenta accreta. Prenat. Diagn. 2014, 34, 159–162. [Google Scholar] [CrossRef]
  7. Laursen, L.S.; Kjaer-Sorensen, K.; Andersen, M.H.; Oxvig, C. Regulation of insulin-like growth factor (IGF) bioactivity by sequential proteolytic cleavage of IGF binding protein-4 and -5. Mol. Endocrinol. 2007, 21, 1246–1257. [Google Scholar] [CrossRef] [PubMed]
  8. Giudice, L.C.; Conover, C.A.; Bale, L.; Faessen, G.H.; Ilg, K.; Sun, I.; Imani, B.; Suen, L.F.; Irwin, J.C.; Christiansen, M.; et al. Identification and regulation of the IGFBP-4 protease and its physiological inhibitor in human trophoblasts and endometrial stroma: Evidence for paracrine regulation of IGF-II bioavailability in the placental bed during human implantation. J. Clin. Endocrinol. Metab. 2002, 87, 2359–2366. [Google Scholar] [CrossRef] [PubMed]
  9. Jordan, N.V.; Johnson, G.L.; Abell, A.N. Tracking the intermediate stages of epithelial-mesenchymal transition in epithelial stem cells and cancer. Cell Cycle 2011, 10, 2865–2873. [Google Scholar] [CrossRef]
  10. Barbosa, C.C.; Calhoun, S.H.; Wieden, H. Non-coding RNAs: What are we missing? Biochem. Cell Biol. 2020, 98, 23–30. [Google Scholar] [CrossRef]
  11. Khan, S.; Ayub, H.; Khan, T.; Wahid, F. MicroRNA biogenesis, gene silencing mechanisms and role in breast, ovarian and prostate cancer. Biochimie 2019, 167, 12–24. [Google Scholar] [CrossRef]
  12. Lee, H.J. Additional stories of microRNAs. Exp. Biol. Med. 2014, 239, 1275–1279. [Google Scholar] [CrossRef]
  13. Denis, H.; Ndlovu, M.N.; Fuks, F. Regulation of mammalian DNA methyltransferases: A route to new mechanisms. EMBO Rep. 2011, 12, 647–656. [Google Scholar] [CrossRef]
  14. Arif, K.M.T.; Elliott, E.K.; Haupt, L.M.; Griffiths, L.R. Regulatory Mechanisms of Epigenetic miRNA Relationships in Human Cancer and Potential as Therapeutic Targets. Cancers 2020, 12, 2922. [Google Scholar] [CrossRef]
  15. Tang, L.; Chen, H.-Y.; Hao, N.-B.; Tang, B.; Guo, H.; Yong, X.; Dong, H.; Yang, S.-M. microRNA inhibitors: Natural and artificial sequestration of microRNA. Cancer Lett. 2017, 407, 139–147. [Google Scholar] [CrossRef]
  16. Garg, M. Targeting microRNAs in epithelial-to-mesenchymal transition-induced cancer stem cells: Therapeutic approaches in cancer. Expert Opin. Ther. Targets 2015, 19, 285–297. [Google Scholar] [CrossRef]
  17. Kalluri, R.; Weinberg, R.A. The basics of epithelial-mesenchymal transition. J. Clin. Investig. 2009, 119, 1420–1428. [Google Scholar] [CrossRef] [PubMed]
  18. Jauniaux, E.; Burton, G.J. Pathophysiology of placenta accreta spectrum disorders: A review of current findings. Clin. Obstet. Gynecol. 2018, 61, 743–754. [Google Scholar] [CrossRef] [PubMed]
  19. Duzyj, C.M.; Buhimschi, I.A.; Motawea, H.; Laky, C.A.; Cozzini, G.; Zhao, G.; Funai, E.F.; Buhimschi, C.S. The invasive phenotype of placenta accreta extravillous trophoblasts associates with loss of E-cadherin. Placenta 2015, 36, 645–651. [Google Scholar] [CrossRef]
  20. Peng, M.; Deng, J.; Zhou, S.; Tao, T.; Su, Q.; Yang, X.; Yang, X. The role of Clusterin in cancer metastasis. Cancer Manag. Res. 2019, 11, 2405–2414. [Google Scholar] [CrossRef] [PubMed]
  21. Wilson, M.R.; Zoubeidi, A. Clusterin as a therapeutic target. Expert Opin. Ther. Targets 2017, 21, 201–213. [Google Scholar] [CrossRef] [PubMed]
  22. Rodríguez-Rivera, C.; Garcia, M.M.; Molina-Álvarez, M.; González-Martín, C.; Goicoechea, C. Clusterin: Always protecting. Synthesis, function and potential issues. Biomed. Pharmacother. 2021, 134, 111174. [Google Scholar] [CrossRef] [PubMed]
  23. Satapathy, S.; Wilson, M.R. The Dual Roles of Clusterin in Extracellular and Intracellular Proteostasis. Trends Biochem. Sci. 2021. [Google Scholar] [CrossRef] [PubMed]
  24. Rohne, P.; Prochnow, H.; Koch-Brandt, C. The CLU-files: Disentanglement of a mystery. Biomol. Concepts 2016, 7, 1–15. [Google Scholar] [CrossRef] [PubMed]
  25. Rizzi, F.; Coletta, M.; Bettuzzi, S. Clusterin (CLU): From One Gene and Two Transcripts to Many Proteins, 1st ed.; Elsevier Inc.: Philadelphia, PA, USA, 2009; Volume 104. [Google Scholar]
  26. Pucci, S.; Bettuzzi, S. The Shifting Balance between CLU Forms during Tumor Progression, 1st ed.; Elsevier Inc.: Philadelphia, PA, USA, 2009; Volume 104. [Google Scholar]
  27. Tellez, T.; Garcia-Aranda, M.; Redondo, M. The Role of Clusterin in Carcinogenesis and its Potential Utility as Therapeutic Target. Curr. Med. Chem. 2016, 23, 4297–4308. [Google Scholar] [CrossRef] [PubMed]
  28. Pucci, S.; Mazzarelli, P.; Nucci, C.; Ricci, F.; Spagnoli, L.G. CLU “In and Out”: Looking for a Link. In Advances in Cancer Research; Academic Press: New York, NY, USA, 2009; Volume 105, pp. 93–113. ISBN 0065-230X. [Google Scholar]
  29. Massimiani, M.; Lacconi, V.; La Civita, F.; Ticconi, C.; Rago, R.; Campagnolo, L. Molecular Signaling Regulating Endometrium-Blastocyst Crosstalk. Int. J. Mol. Sci. 2019, 21, 23. [Google Scholar] [CrossRef]
  30. Shin, J.K.; Han, K.A.; Kang, M.Y.; Kim, Y.S.; Park, J.K.; Choi, W.J.; Lee, S.A.; Lee, J.H.; Choi, W.S.; Paik, W.Y. Expression of clusterin in normal and preeclamptic placentas. J. Obstet. Gynaecol. Res. 2008, 34, 473–479. [Google Scholar] [CrossRef]
  31. Watanabe, H.; Hamada, H.; Yamada, N.; Sohda, S.; Yamakawa-Kobayashi, K.; Yoshikawa, H.; Arinami, T. Proteome analysis reveals elevated serum levels of clusterin in patients with preeclampsia. Proteomics 2004, 4, 537–543. [Google Scholar] [CrossRef]
  32. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, 10. [Google Scholar] [CrossRef]
  33. Team, R.C. A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. Available online: https://www.r-project.org (accessed on 10 March 2021).
  34. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  35. Schägger, H. Tricine–SDS-PAGE. Nat. Protoc. 2006, 1, 16–23. [Google Scholar] [CrossRef]
  36. Team, R. RStudio: Integrated Development for R. RStudio. Available online: http://www.rstudio.com/ (accessed on 23 March 2021).
  37. McLaughlin, M.J.; Sainani, K.L. Bonferroni, Holm, and Hochberg corrections: Fun names, serious changes to p values. Phys. Med. Rehabil. 2014, 6, 544–546. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, R.; Wang, G. Protein Modification and Autophagy Activation. Adv. Exp. Med. Biol. 2019, 1206, 237–259. [Google Scholar] [CrossRef] [PubMed]
  39. Yang, X.; Qian, K. Protein O-GlcNAcylation: Emerging mechanisms and functions. Nat. Rev. Mol. Cell Biol. 2017, 18, 452–465. [Google Scholar] [CrossRef]
  40. Tang, M.K.S.; Yue, P.Y.K.; Ip, P.P.; Huang, R.L.; Lai, H.C.; Cheung, A.N.Y.; Tse, K.Y.; Ngan, H.Y.S.; Wong, A.S.T. Soluble E-cadherin promotes tumor angiogenesis and localizes to exosome surface. Nat. Commun. 2018, 9, 1–15. [Google Scholar] [CrossRef] [PubMed]
  41. Tang, M.K.S.; Ip, P.P.; Wong, A.S.T. New insights into the role of soluble E-cadherin in tumor angiogenesis. Cell Stress 2018, 2, 236–238. [Google Scholar] [CrossRef] [PubMed]
  42. Theodoraki, M.-N.; Yerneni, S.S.; Brunner, C.; Theodorakis, J.; Hoffmann, T.K.; Whiteside, T.L. Plasma-derived Exosomes Reverse Epithelial-to-Mesenchymal Transition after Photodynamic Therapy of Patients with Head and Neck Cancer. Oncoscience 2018, 5, 75. [Google Scholar] [CrossRef]
  43. David, J.M.; Rajasekaran, A.K. Dishonorable discharge: The oncogenic roles of cleaved E-cadherin fragments. Cancer Res. 2012, 72, 2917–2923. [Google Scholar] [CrossRef]
  44. Christou, N.; Perraud, A.; Blondy, S.; Jauberteau, M.-O.; Battu, S.; Mathonnet, M. The extracellular domain of E cadherin linked to invasiveness in colorectal cancer: A new resistance and relapses monitoring serum-bio marker? J. Cancer Res. Clin. Oncol. 2017, 143, 1177–1190. [Google Scholar] [CrossRef]
  45. French, M.E.; Koehler, C.F.; Hunter, T. Emerging functions of branched ubiquitin chains. Cell Discov. 2021, 7, 6. [Google Scholar] [CrossRef]
  46. Madadi, S.; Schwarzenbach, H.; Lorenzen, J.; Soleimani, M. MicroRNA expression studies: Challenge of selecting reliable reference controls for data normalization. Cell. Mol. Life Sci. 2019, 76, 3497–3514. [Google Scholar] [CrossRef]
  47. Xie, H.A.N.; Qiao, P.; Lu, Y.I.; Li, Y.; Tang, Y. Increased expression of high mobility group box protein 1 and vascular endothelial growth factor in placenta previa. Mol. Med. Rep. 2017, 9051–9059. [Google Scholar] [CrossRef]
  48. Han, Q.; Zheng, L.; Liu, Z. Expression of b-catenin in human trophoblast and its role in placenta accreta and placenta previa. J. Int. Med. Res. 2019. [Google Scholar] [CrossRef]
  49. Bartels, H.C.; Postle, J.D.; Downey, P.; Brennan, D.J. Placenta Accreta Spectrum: A Review of Pathology, Molecular Biology, and Biomarkers. Dis. Markers 2018, 2018, 1507674. [Google Scholar] [CrossRef] [PubMed]
  50. Das, V.; Bhattacharya, S.; Chikkaputtaiah, C.; Hazra, S.; Pal, M. The basics of epithelial–mesenchymal transition (EMT): A study from a structure, dynamics, and functional perspective. J. Cell. Physiol. 2019, 234, 14535–14555. [Google Scholar] [CrossRef] [PubMed]
  51. Greening, D.W.; Gopal, S.K.; Mathias, R.A.; Liu, L.; Sheng, J.; Zhu, H.-J.; Simpson, R.J. Emerging roles of exosomes during epithelial-mesenchymal transition and cancer progression. Semin. Cell Dev. Biol. 2015, 40, 60–71. [Google Scholar] [CrossRef] [PubMed]
  52. Devaux, C.A.; Mezouar, S.; Mege, J.-L. The E-Cadherin Cleavage Associated to Pathogenic Bacteria Infections Can Favor Bacterial Invasion and Transmigration, Dysregulation of the Immune Response and Cancer Induction in Humans. Front. Microbiol. 2019, 10, 2598. [Google Scholar] [CrossRef] [PubMed]
  53. Grabowska, M.M.; Day, M.L. Soluble E-cadherin: More than a symptom of disease. Front. Biosci. Landmark Ed. 2012, 17, 1948–1964. [Google Scholar] [CrossRef] [PubMed]
  54. Ito, K.; Okamoto, I.; Araki, N.; Kawano, Y.; Nakao, M.; Fujiyama, S.; Tomita, K.; Mimori, T.; Saya, H. Calcium influx triggers the sequential proteolysis of extracellular and cytoplasmic domains of E-cadherin, leading to loss of beta-catenin from cell-cell contacts. Oncogene 1999, 18, 7080–7090. [Google Scholar] [CrossRef]
  55. Marambaud, P.; Shioi, J.; Serban, G.; Georgakopoulos, A.; Sarner, S.; Nagy, V.; Baki, L.; Wen, P.; Efthimiopoulos, S.; Shao, Z.; et al. A presenilin-1/gamma-secretase cleavage releases the E-cadherin intracellular domain and regulates disassembly of adherens junctions. EMBO J. 2002, 21, 1948–1956. [Google Scholar] [CrossRef]
  56. Steinhusen, U.; Weiske, J.; Badock, V.; Tauber, R.; Bommert, K.; Huber, O. Cleavage and shedding of E-cadherin after induction of apoptosis. J. Biol. Chem. 2001, 276, 4972–4980. [Google Scholar] [CrossRef]
  57. Mayerle, J.; Schnekenburger, J.; Krüger, B.; Kellermann, J.; Ruthenbürger, M.; Weiss, F.U.; Nalli, A.; Domschke, W.; Lerch, M.M. Extracellular cleavage of E-cadherin by leukocyte elastase during acute experimental pancreatitis in rats. Gastroenterology 2005, 129, 1251–1267. [Google Scholar] [CrossRef]
  58. Loh, C.-Y.; Chai, J.Y.; Tang, T.F.; Wong, W.F.; Sethi, G.; Shanmugam, M.K.; Chong, P.P.; Looi, C.Y. The E-Cadherin and N-Cadherin Switch in Epithelial-to-Mesenchymal Transition: Signaling, Therapeutic Implications, and Challenges. Cells 2019, 8, 1118. [Google Scholar] [CrossRef]
  59. Hall, M.K.; Weidner, D.A.; Dayal, S.; Schwalbe, R.A. Cell surface N-glycans influence the level of functional E-cadherin at the cell-cell border. FEBS Open Bio 2014, 4, 892–897. [Google Scholar] [CrossRef]
  60. Chen, A.; Beetham, H.; Black, M.A.; Priya, R.; Telford, B.J.; Guest, J.; Wiggins, G.A.R.; Godwin, T.D.; Yap, A.S.; Guilford, P.J. E-cadherin loss alters cytoskeletal organization and adhesion in non-malignant breast cells but is insufficient to induce an epithelial-mesenchymal transition. BMC Cancer 2014, 14, 552. [Google Scholar] [CrossRef]
  61. Hollestelle, A.; Peeters, J.K.; Smid, M.; Timmermans, M.; Verhoog, L.C.; Westenend, P.J.; Heine, A.A.J.; Chan, A.; Sieuwerts, A.M.; Wiemer, E.A.C.; et al. Loss of E-cadherin is not a necessity for epithelial to mesenchymal transition in human breast cancer. Breast Cancer Res. Treat. 2013, 138, 47–57. [Google Scholar] [CrossRef] [PubMed]
  62. Shannan, B.; Seifert, M.; Leskov, K.; Willis, J.; Boothman, D.; Tilgen, W.; Reichrath, J. Challenge and promise: Roles for clusterin in pathogenesis, progression and therapy of cancer. Cell Death Differ. 2006, 13, 12–19. [Google Scholar] [CrossRef] [PubMed]
  63. Pucci, S.; Bonanno, E.; Pichiorri, F.; Angeloni, C.; Spagnoli, L.G. Modulation of different clusterin isoforms in human colon tumorigenesis. Oncogene 2004, 23, 2298–2304. [Google Scholar] [CrossRef] [PubMed]
  64. Criswell, T.; Klokov, D.; Beman, M.; Lavik, J.P.; Boothman, D.A. Repression of IR-inducible clusterin expression by the p53 tumor suppressor protein. Cancer Biol. Ther. 2003, 2, 372–380. [Google Scholar] [CrossRef] [PubMed]
  65. Cervellera, M.; Raschella, G.; Santilli, G.; Tanno, B.; Ventura, A.; Mancini, C.; Sevignani, C.; Calabretta, B.; Sala, A. Direct transactivation of the anti-apoptotic gene apolipoprotein J (clusterin) by B-MYB. J. Biol. Chem. 2000, 275, 21055–21060. [Google Scholar] [CrossRef]
  66. Santilli, G.; Aronow, B.J.; Sala, A. Essential requirement of apolipoprotein J (clusterin) signaling for IkappaB expression and regulation of NF-kappaB activity. J. Biol. Chem. 2003, 278, 38214–38219. [Google Scholar] [CrossRef]
  67. Shen, K.-H.; Liao, A.C.-H.; Hung, J.-H.; Lee, W.-J.; Hu, K.-C.; Lin, P.-T.; Liao, R.-F.; Chen, P.-S. α-Solanine inhibits invasion of human prostate cancer cell by suppressing epithelial-mesenchymal transition and MMPs expression. Molecules 2014, 19, 11896–11914. [Google Scholar] [CrossRef] [PubMed]
  68. Hu, J.L.; Wang, W.; Lan, X.L.; Zeng, Z.C.; Liang, Y.S.; Yan, Y.R.; Song, F.Y.; Wang, F.F.; Zhu, X.H.; Liao, W.J.; et al. CAFs secreted exosomes promote metastasis and chemotherapy resistance by enhancing cell stemness and epithelial-mesenchymal transition in colorectal cancer. Mol. Cancer 2019, 18, 91. [Google Scholar] [CrossRef]
  69. Yang, B.; Feng, X.; Liu, H.; Tong, R.; Wu, J.; Li, C.; Yu, H.; Chen, Y.; Cheng, Q.; Chen, J.; et al. High-metastatic cancer cells derived exosomal miR92a-3p promotes epithelial-mesenchymal transition and metastasis of low-metastatic cancer cells by regulating PTEN/Akt pathway in hepatocellular carcinoma. Oncogene 2020, 39, 6529–6543. [Google Scholar] [CrossRef] [PubMed]
  70. Wan, W.; Wan, W.; Long, Y.; Li, Q.; Jin, X.; Wan, G.; Zhang, F.; Lv, Y.; Zheng, G.; Li, Z.; et al. MiR-25-3p promotes malignant phenotypes of retinoblastoma by regulating PTEN/Akt pathway. Biomed. Pharmacother. 2019, 118, 109111. [Google Scholar] [CrossRef]
  71. Konoshenko, M.; Sagaradze, G.; Orlova, E.; Shtam, T.; Proskura, K.; Kamyshinsky, R.; Yunusova, N.; Alexandrova, A.; Efimenko, A.; Tamkovich, S. Total Blood Exosomes in Breast Cancer: Potential Role in Crucial Steps of Tumorigenesis. Int. J. Mol. Sci. 2020, 21, 7341. [Google Scholar] [CrossRef] [PubMed]
  72. Xu, P.; Zhao, Y.; Liu, M.; Wang, Y.; Wang, H.; Li, Y.; Zhu, X.; Yao, Y.; Wang, H.; Qiao, J.; et al. Variations of MicroRNAs in Human Placentas and Plasma From Preeclamptic Pregnancy. Hypertension 2014, 63, 1276–1284. [Google Scholar] [CrossRef]
  73. Choi, S.-Y.; Yun, J.; Lee, O.-J.; Han, H.-S.; Yeo, M.-K.; Lee, M.-A.; Suh, K.-S. MicroRNA expression profiles in placenta with severe preeclampsia using a PNA-based microarray. Placenta 2013, 34, 799–804. [Google Scholar] [CrossRef]
  74. Xie, N.; Jia, Z.; Li, L. miR-320a upregulation contributes to the development of preeclampsia by inhibiting the growth and invasion of trophoblast cells by targeting interleukin 4. Mol. Med. Rep. 2019. [Google Scholar] [CrossRef]
  75. Jauniaux, E.; Kingdom, J.C.; Silver, R.M. A comparison of recent guidelines in the diagnosis and management of placenta accreta spectrum disorders. Best Pract. Res. Clin. Obstet. Gynaecol. 2020. [Google Scholar] [CrossRef]
  76. Chen, S.; Pang, D.; Li, Y.; Zhou, J.; Liu, Y.; Yang, S.; Liang, K.; Yu, B. Serum miRNA biomarker discovery for placenta accreta spectrum. Placenta 2020, 101, 215–220. [Google Scholar] [CrossRef]
Figure 1. Histological analysis of the placenta with pathological invasion: placenta accreta (a), placenta increta (b), placenta percreta (c). In the images of histological sections (ac), stained with hematoxylin and eosin, “V” denotes chorionic villi and ”M” denotes myometrium. In image (a) with 100× magnification, the mononuclear intermediate trophoblast is marked with a black arrow, and the fibrinoid is marked with a red arrow; in image (b) with 100× magnification and in image (c) with 40× magnification, the arrow indicates the lumen of the myometrium blood vessel.
Figure 1. Histological analysis of the placenta with pathological invasion: placenta accreta (a), placenta increta (b), placenta percreta (c). In the images of histological sections (ac), stained with hematoxylin and eosin, “V” denotes chorionic villi and ”M” denotes myometrium. In image (a) with 100× magnification, the mononuclear intermediate trophoblast is marked with a black arrow, and the fibrinoid is marked with a red arrow; in image (b) with 100× magnification and in image (c) with 40× magnification, the arrow indicates the lumen of the myometrium blood vessel.
Life 11 00270 g001
Figure 2. Amplification curves of hsa-miR-382-5p in 64 analyzed samples of blood plasma with indication of a threshold value, the point at which fluorescence reaches values above baseline levels.
Figure 2. Amplification curves of hsa-miR-382-5p in 64 analyzed samples of blood plasma with indication of a threshold value, the point at which fluorescence reaches values above baseline levels.
Life 11 00270 g002
Figure 3. Principal component analysis (PCA) plot based on the miRNA dataset.
Figure 3. Principal component analysis (PCA) plot based on the miRNA dataset.
Life 11 00270 g003
Figure 4. Quantitative RT-PCR data on the miRNA expression level (ΔCt values) in the peripheral blood plasma samples from pregnant women of groups I–V.
Figure 4. Quantitative RT-PCR data on the miRNA expression level (ΔCt values) in the peripheral blood plasma samples from pregnant women of groups I–V.
Life 11 00270 g004
Figure 5. Western blot analysis of soluble E-cadherin in peripheral blood plasma pooled samples from groups I–V. Blot was stained with Ponceau S dye (a) and with antibody against the extracellular domain of E-cadherin (b). Мolecular weights (kDa) of the protein ladder and fragments of soluble E-cadherin are indicated on the blot.
Figure 5. Western blot analysis of soluble E-cadherin in peripheral blood plasma pooled samples from groups I–V. Blot was stained with Ponceau S dye (a) and with antibody against the extracellular domain of E-cadherin (b). Мolecular weights (kDa) of the protein ladder and fragments of soluble E-cadherin are indicated on the blot.
Life 11 00270 g005
Figure 6. Western blot analysis of the α-subunit of the clusterin secretory form in peripheral blood plasma samples from groups I–V. Blot with pooled samples from groups I-V was stained with Ponceau S dye (a) and with antibody against the α-subunit of the clusterin secretory form (b). Blots with individual samples from each group I–V were stained with antibody against the α-subunit of the clusterin secretory form (c), wherein arabic numerals, written with a hyphen next to the group number, indicate the sample numbers. Мolecular weights (kDa) of the protein ladder and clusterin are indicated on the blots.
Figure 6. Western blot analysis of the α-subunit of the clusterin secretory form in peripheral blood plasma samples from groups I–V. Blot with pooled samples from groups I-V was stained with Ponceau S dye (a) and with antibody against the α-subunit of the clusterin secretory form (b). Blots with individual samples from each group I–V were stained with antibody against the α-subunit of the clusterin secretory form (c), wherein arabic numerals, written with a hyphen next to the group number, indicate the sample numbers. Мolecular weights (kDa) of the protein ladder and clusterin are indicated on the blots.
Life 11 00270 g006
Figure 7. The expression level of the α-subunit of the clusterin secretory form in the peripheral blood plasma from women of groups I–V analyzed by Western blotting. Data are presented as Me ± standard deviation (stdev).
Figure 7. The expression level of the α-subunit of the clusterin secretory form in the peripheral blood plasma from women of groups I–V analyzed by Western blotting. Data are presented as Me ± standard deviation (stdev).
Life 11 00270 g007
Figure 8. Receiver operating characteristic (ROC) curves of the logistic regression models in case of entire cohort of patients (a,c) and optimized cohort of patients (b,d) for placenta accreta (a,b) and placenta increta (c,d).
Figure 8. Receiver operating characteristic (ROC) curves of the logistic regression models in case of entire cohort of patients (a,c) and optimized cohort of patients (b,d) for placenta accreta (a,b) and placenta increta (c,d).
Life 11 00270 g008aLife 11 00270 g008b
Table 1. Clinical characteristics of the pregnant women.
Table 1. Clinical characteristics of the pregnant women.
Group of Pregnant Women 1AIPScar on the UterusPlacenta PreviaAge 2Number of Pregnancies in History 2Number of Previous Cesarean Sections 2
I (n = 10)absyesno36.5 (33; 39)3 (3; 3.8)1 (1; 1.8)
IIa (n = 9)absnoyes36 (36; 39)3 (1; 3)0 (0; 0) **
IIb (n = 8)absyesyes37.5 (34; 38.3)3 (2.8; 6.3)2 (1; 2.3)
III (n = 10)accretayesyes34 (33; 35.8)4 (3.3; 5)2 (2; 3) *
IV (n = 20)incretayesyes33 (29.5; 36.3)4 (3; 5)2 (1; 2)
V (n = 7)percretayesyes33 (32; 33)4 (2.5; 4)2 (1.5; 2)
1 Roman numerals indicate the group number; the size of the group is indicated in brackets. 2 Data are presented as the median (Me) and quartiles Q1 and Q3 in the format: Me (Q1; Q3). * The value of statistical significance of differences when compared to group I is less than 0.05. ** The value of statistical significance of differences when compared to group I is less than 0.0001.
Table 2. miRNA parameters.
Table 2. miRNA parameters.
miRNAmiRNA Gene TargetmiRNA Accession Number (miRBase)Nucleotide Sequence of Sense Primer for PCR, 5′-3′PCR Primers Annealing Temperature, ℃
miR-320a-3pCDH, CLUMIMAT0000510aaaagctgggttgagagggcga59.5
miR-17-5pCDH, CLUMIMAT0000070caaagtgcttacagtgcaggtag55
miR-21-5pCDH, CLUMIMAT0000076tagcttatcagactgatgttga54.6
miR-1323CDH, CLUMIMAT0005795tcaaaactgaggggcattttct51
miR-25-3pCDH, CLUMIMAT0000081cattgcacttgtctcggtctga56
miR-138-5pCDH, CLUMIMAT0000430agctggtgttgtgaatcaggccg54.6
miR-34a-5pCDH, CLUMIMAT0000255tggcagtgtcttagctggttgt51.3
miR-92a-3pCDH, CLUMIMAT0000092tattgcacttgtcccggcctgt60
miR-30a-5pCLUMIMAT0000087tgtaaacatcctcgactggaag54.6
miR-30c-5pCLUMIMAT0000244tgtaaacatcctacactctcagc49.1
miR-371a-5pCDHMIMAT0004687actcaaactgtgggggcact54.6
miR-506-3pCDHMIMAT0002878taaggcacccttctgagtaga51.3
miR-382-5pendogenous control for PCRMIMAT0000737gaagttgttcgtggtggattcg49.1
Table 3. Pairwise comparison of groups II–V to group I by the miRNA ΔCt value for an optimized cohort of patients.
Table 3. Pairwise comparison of groups II–V to group I by the miRNA ΔCt value for an optimized cohort of patients.
miRNAGroupМеQ1Q3Mann–Whitney Test, p-ValueBonferroni Adjustments for 20 Tests (5 miRNAs & 4 Comparisons), p-Value
miR-17-5pI10.956.2311.04
miR-17-5pIIa + IIb11.0710.8711.190.0551631
miR-17-5pIII86.5110.850.1359111
miR-17-5pIV7.025.4410.750.0408580.81716
miR-17-5pV9.675.249.750.0548951
miR-21-5pI10.897.2311.03
miR-21-5pIIa + IIb8.557.411.060.7258521
miR-21-5pIII5.24.425.93<0.00010.0198
miR-21-5pIV5.394.465.89<0.0001<0.0001
miR-21-5pV0.670.40.990.0001750.0035
miR-25-3pI7.286.338.53
miR-25-3pIIa + IIb8.777.1110.870.3175471
miR-25-3pIII5.945.066.350.0770051
miR-25-3pIV6.344.719.450.2895031
miR-25-3pV1.650.912.090.0001750.0035
miR-320a-3pI43.394.49
miR-320a-3pIIa + IIb3.913.664.660.8152531
miR-320a-3pIII2.752.292.940.0004940.00988
miR-320a-3pIV2.842.63.22<0.00010.00092
miR-320a-3pV1.010.921.690.0001750.0035
miR-92a-3pI3.743.714.04
miR-92a-3pIIa + IIb4.173.754.40.2152281
miR-92a-3pIII2.822.6130.0001650.0033
miR-92a-3pIV2.982.623.640.0251470.50294
miR-92a-3pV−0.01−0.580.480.0001750.0035
Data are presented as the median (Me) of the ΔCt values and quartiles Q1 and Q3 and the statistical significance of the differences p while applying Mann–Whitney test with Bonferroni correction for multiple testing.
Table 4. Holm–Bonferroni test for pairwise comparison of groups II–V to group I by the miRNA ΔCt value from Table 3.
Table 4. Holm–Bonferroni test for pairwise comparison of groups II–V to group I by the miRNA ΔCt value from Table 3.
miRNAI vs. II, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni Test, p-ValueSignificance 1
miR-17-5p0.0551610.010000
miR-92a-3p0.2152320.012500
miR-25-3p0.3175530.016670
miR-21-5p0.7258540.025000
miR-320a-3p0.8152550.050000
miRNAI vs. III, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni Test, p-ValueSignificance 1
miR-21-5p0.000110.010001
miR-92a-3p0.000220.012501
miR-320a-3p0.000530.016671
miR-25-3p0.077040.025000
miR-17-5p0.135950.050000
miRNAI vs. IV, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni Test, p-ValueSignificance 1
miR-21-5p0.0001010.010001
miR-320a-3p0.0001020.012501
miR-92a-3p0.0251530.016670
miR-17-5p0.0408640.025000
miR-25-3p0.2895050.050000
miRNAI vs. V, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni Test, p-ValueSignificance 1
miR-21-5p0.0001810.010001
miR-92a-3p0.0001820.012501
miR-320a-3p0.0001830.016671
miR-25-3p0.0001840.025001
miR-17-5p0.0549050.050000
1 “1” means statistically significant difference, “0” means insignificant difference.
Table 5. Pairwise comparison of groups III–V to group II by miRNA ΔCt values for the optimized cohort of patients.
Table 5. Pairwise comparison of groups III–V to group II by miRNA ΔCt values for the optimized cohort of patients.
miRNAGroupМеQ1Q3Mann–Whitney Test, p-ValueBonferroni Adjustments for 15 Tests (5 miRNAs & 3 Comparisons)
miR-17-5pIIa + IIb11.0710.8711.19
miR-17-5pIII86.5110.850.001670.02499
miR-17-5pIV7.025.4410.750.000140.0021
miR-17-5pV9.675.249.750.002890.043275
miR-21-5pIIa + IIb8.557.411.06
miR-21-5pIII5.24.425.930.000150.00222
miR-21-5pIV5.394.465.89<0.0001<0.0001
miR-21-5pV0.670.40.99<0.0001<0.0001
miR-25-3pIIa + IIb8.777.1110.87
miR-25-3pIII5.945.066.350.008260.12383
miR-25-3pIV6.344.719.450.089201
miR-25-3pV1.650.912.09<0.00010.00117
miR-320a-3pIIa + IIb3.913.664.66
miR-320a-3pIII2.752.292.940.000290.00441
miR-320a-3pIV2.842.63.22<0.00010.00396
miR-320a-3pV1.010.921.69<0.0001<0.0001
miR-92a-3pIIa + IIb4.1710.8711.19
miR-92a-3pIII2.826.5110.850.000150.00222
miR-92a-3pIV11.075.4410.750.003160.04733
miR-92a-3pV85.249.75<0.00010.00093
Data are presented as the median (Me) of the ΔCt values and quartiles Q1 and Q3 and the statistical significance of the differences p while applying Mann–Whitney test with Bonferroni correction for multiple testing.
Table 6. Holm–Bonferroni test for pairwise comparison of groups II–V to group I by the miRNA ΔCt value from Table 5.
Table 6. Holm–Bonferroni test for pairwise comparison of groups II–V to group I by the miRNA ΔCt value from Table 5.
miRNAII vs. III, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni Test, p-ValueSignificance
miR-21-5p0.0001510.011
miR-92a-3p0.0001520.01251
miR-320a-3p0.0002930.016671
miR-25-3p0.0082640.0251
miR-17-5p0.0016750.051
miRNAII vs. IV, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni Test, p-ValueSignificance
miR-21-5p<0.000110.011
miR-320a-3p<0.000120.01251
miR-17-5p0.0001430.016671
miR-92a-3p0.0031640.0251
miR-25-3p0.0892050.050
miRNAII vs. V, Mann–Whitney Test, p-ValueOrderHolm–Bonferroni test, p-ValueSignificance
miR-21-5p<0.000110.011
miR-320a-3p<0.000120.01251
miR-92a-3p<0.000130.016671
miR-25-3p<0.000140.0251
miR-17-5p0.0028950.051
“1” means a statistically significant difference, “0” means an insignificant difference.
Table 7. Deep sequencing data on the miRNA expression level in the placenta in the area of pathological invasion (P-area).
Table 7. Deep sequencing data on the miRNA expression level in the placenta in the area of pathological invasion (P-area).
miRNA 1III_P1III_P2III_P3IV_P4IV_P5IV_P8V_P6V_P7FC (IV/III)FC (V/III)p (IV vs. III)p (V vs. III)
miR-17-5p25525532826810001287173415673.96.50.026<0.001
miR-21-5p18,63511,6908478652212,44916,98023,70325,6531.12.10.8510.030
miR-25-3p298819541830128047,99182,55675,469108,39324.647.0<0.001<0.001
miR-92a-3p19,47019,21534,32023,702453,638590,837857,646911,57923.345.4<0.001<0.001
miR-320a-3p10,18311,13411,30914,03319,97821,48939,84933,0921.83.30.023<0.001
1 III, IV, and V are the group numbers; P1–P8 are the numbers of placenta samples from the P-area; FC is the fold change in the expression level as the ratio of the medians of the miRNA read numbers in the compared groups; p—statistical significance of differences.
Table 8. Deep sequencing data on miRNA expression level in the placenta in the area outside of pathological invasion (N-area).
Table 8. Deep sequencing data on miRNA expression level in the placenta in the area outside of pathological invasion (N-area).
miRNA 1III_N1III_N2III_N3IV_N4IV_N5IV_N8V_N6V_N7FC (IV/III)FC (V/III)p (IV vs. III)p (V vs. III)
miR-17-5p30121821220874111227673893.42.70.1040.048
miR-21-5p82906505626511,06617,85615,61726,97871422.42.60.1580.064
miR-25-3p245117161222121140,05484,35049,601294023.315.3<0.001<0.001
miR-92a-3p15,61218,05327,16822,394288,269722,332647,51314,94016.018.3<0.0010.002
miR-320a-3p6485861110,78115,28730,89322,91829,83155132.72.10.0750.191
1 III, IV, and V are the group numbers; P1–P8 are the numbers of placental samples from the N-area; FC is the fold change in the expression level as the ratio of the medians of the miRNA read numbers in the compared groups; p—statistical significance of differences.
Table 9. Pairwise comparison of groups III–V to group I or group II by the expression level of the α-subunit of the clusterin secretory form in the peripheral blood plasma.
Table 9. Pairwise comparison of groups III–V to group I or group II by the expression level of the α-subunit of the clusterin secretory form in the peripheral blood plasma.
IIIIIIIVV
Me1.137271.553840.705640.622910.54514
stdev0.210910.699330.112780.091990.21387
p * 0.06561<0.001<0.001<0.001
p ** 0.004050.002220.00167
* statistical significance of differences between groups II–V with group I. ** statistical significance of differences between groups III–V with group II.
Table 10. Parameters of the logistic regression models.
Table 10. Parameters of the logistic regression models.
BiomarkerAUCSpSeCutoffiК
Placenta accreta
“Clusterin”1110.51832.89−2232.65
“miR-21-5p +
miR-92a-3p +
miR-320a-3р”
1110.54526.11−141.53
−868.31
−244.95
“miR-21-5p +
miR-92a-3p”
0.995110.26264.43−3.32
−13.08
“miR-92a-3p +
miR-320a-3p”
0.9860.95810.27522.9−4.02
−3.25
“miR-21-5p +
miR-320a-3р”
0.95810.8880.42511.39−0.66
−2.19
Placenta increta
“Clusterin”1110.5287.91−372.61
“miR-21-5p +
miR-320a-3p”
0.9810.95810.26316.70−1.25
−2.41
“miR-21-5p +
miR-92a-3p +
miR-320a-3р”
0.9810.91610.27716.96−1.24
−0.16
−2.34
“miR-21-5p”0.95810.9330.7312.35−1.89
“miR-21-5p +
miR-92a-3p”
0.95810.9330.54515.59−1.68
−1.25
Placenta percreta
“miR-320a-3p”1110.5828.41−298.44
AUC—area under the curve; Sp—specificity; Se—sensitivity; cutoff—intersection point at which there is a balance between sensitivity and specificity; i and K—coefficients used in Equation (1).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Timofeeva, A.V.; Fedorov, I.S.; Pirogova, M.M.; Vasilchenko, O.N.; Chagovets, V.V.; Ezhova, L.S.; Zabelina, T.M.; Shmakov, R.G.; Sukhikh, G.T. Clusterin and Its Potential Regulatory microRNAs as a Part of Secretome for the Diagnosis of Abnormally Invasive Placenta: Accreta, Increta, and Percreta Cases. Life 2021, 11, 270. https://doi.org/10.3390/life11040270

AMA Style

Timofeeva AV, Fedorov IS, Pirogova MM, Vasilchenko ON, Chagovets VV, Ezhova LS, Zabelina TM, Shmakov RG, Sukhikh GT. Clusterin and Its Potential Regulatory microRNAs as a Part of Secretome for the Diagnosis of Abnormally Invasive Placenta: Accreta, Increta, and Percreta Cases. Life. 2021; 11(4):270. https://doi.org/10.3390/life11040270

Chicago/Turabian Style

Timofeeva, Angelika V., Ivan S. Fedorov, Mariya M. Pirogova, Oksana N. Vasilchenko, Vitaliy V. Chagovets, Larisa S. Ezhova, Tatiana M. Zabelina, Roman G. Shmakov, and Gennadiy T. Sukhikh. 2021. "Clusterin and Its Potential Regulatory microRNAs as a Part of Secretome for the Diagnosis of Abnormally Invasive Placenta: Accreta, Increta, and Percreta Cases" Life 11, no. 4: 270. https://doi.org/10.3390/life11040270

APA Style

Timofeeva, A. V., Fedorov, I. S., Pirogova, M. M., Vasilchenko, O. N., Chagovets, V. V., Ezhova, L. S., Zabelina, T. M., Shmakov, R. G., & Sukhikh, G. T. (2021). Clusterin and Its Potential Regulatory microRNAs as a Part of Secretome for the Diagnosis of Abnormally Invasive Placenta: Accreta, Increta, and Percreta Cases. Life, 11(4), 270. https://doi.org/10.3390/life11040270

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop