Next Article in Journal
Complications after Thermal Ablation of Hepatocellular Carcinoma and Liver Metastases: Imaging Findings
Next Article in Special Issue
Cell-Free DNA as a New Biomarker of IVF Success, Independent of Any Infertility Factor, Including Endometriosis
Previous Article in Journal
Development of “Mathematical Technology for Cytopathology,” an Image Analysis Algorithm for Pancreatic Cancer
Previous Article in Special Issue
Endometrium as Control of Endometriosis in Experimental Research: Assessment of Sample Suitability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Endometriosis Associated-miRNome Analysis of Blood Samples: A Prospective Study

1
Clinical Research Group (GRC) Paris 6: Centre Expert Endométriose (C3E), Department of Obstetrics and Reproductive Medicine, Hôpital Tenon, Sorbonne University (GRC6 C3E SU), 4 rue de la Chine, 75020 Paris, France
2
Cancer Biology and Therapeutics INSERM UMR_S_938, Centre de Recherche Saint-Antoine (CRSA), 75020 Paris, France
3
Ziwig, 19 rue Reboud, 69003 Lyon, France
4
Endometriosis Expert Center—Pays de la Loire, Department of Obstetrics and Reproductive Medicine—CHU d’Angers, 49100 Angers, France
5
Endometriosis Expert Center—Steering Committee of the EndAURA Network, Department of Obstetrics and Reproductive Medicine, Lyon South University Hospital, Lyon Civil Hospices, 69310 Pierre Bénite, France
6
Paris Brain Institute—Institut du Cerveau—ICM, Inserm U1127, CNRS UMR 7225, AP-HP—Hôpital Pitié-Salpêtrière, Sorbonne University, 75006 Paris, France
7
Genotyping and Sequencing Core Facility, iGenSeq, Institut du Cerveau et de la Moelle épinière, ICM, Hôpital Pitié-Salpêtrière, 47-83 Boulevard de l’Hôpital, 75013 Paris, France
8
Department of Obstetrics and Reproductive Medicine, Hôpital Tenon, Sorbonne University, 4 rue de la Chine, 75020 Paris, France
*
Author to whom correspondence should be addressed.
Submission received: 10 March 2022 / Revised: 19 April 2022 / Accepted: 27 April 2022 / Published: 5 May 2022
(This article belongs to the Special Issue Advances in the Diagnostics of Endometriosis)

Abstract

:
The aim of our study was to describe the bioinformatics approach to analyze miRNome with Next Generation Sequencing (NGS) of 200 plasma samples from patients with and without endometriosis. Patients were prospectively included in the ENDO-miRNA study that selected patients with pelvic pain suggestive of endometriosis. miRNA sequencing was performed using an Novaseq6000 sequencer (Illumina, San Diego, CA, USA). Small RNA-seq of 200 plasma samples yielded ~4228 M raw sequencing reads. A total of 2633 miRNAs were found differentially expressed. Among them, 8.6% (n = 229) were up- or downregulated. For these 229 miRNAs, the F1-score, sensitivity, specificity, and AUC ranged from 0–88.2%, 0–99.4%, 4.3–100%, and 41.5–68%, respectively. Utilizing the combined bioinformatic and NGS approach, a specific and broad panel of miRNAs was detected as being potentially suitable for building a blood signature of endometriosis.

1. Introduction

The mammalian genome contains sequences for RNAs coding for messenger RNA (mRNA) proteins and non-coding RNAs (ncRNA). ncRNAs represent 98% of the transcriptome [1]. The known ncRNAs are subdivided into long non-coding RNAs (lncRNAs comprising more than 200 nucleotides (nt)) and small non-coding RNAs (sncRNAs) comprising less than 50 nucleotides [2,3]. To date, miRNAs have been studied far more than their non-coding counterparts. However, during the last decades, a growing interest appears for ncRNAs because of their implication in many benign, malignant pathologies and also in neurodegenerative diseases.
Approximately 70% of studies evaluated exosomes as the source of choice for ncRNAs [4]. The RNA content in the exosomes is estimated at 40.4% mature miRNAs, 40% piwi-interacting RNAs (piRNAs), 3.7% pseudogenes, 2.4% lncRNAs, tRNAs at 2.1%, and mRNAs at 2.1% of total RNA [5]. Among sncRNAs, numerous studies have focused on the role of miRNAs, which are detectable in various body fluids, as potential biomarkers for various pathologies [6,7,8]. Currently, more than 2600 human miRNAs have been annotated [9,10,11,12]. miRNAs are single-stranded conserved sncRNAs composed of 21–25 nucleotides playing a pivotal role in gene degradation and translation by binding to their complementary messenger RNA (mRNA) [13]. The human miRNA spectrum varies according to cell type, tissue type, developmental stage, environmental factors, and health/disease state and disease stages [14,15,16,17]. Recently, numerous studies have demonstrated the values of miRNAs in various cancers and systemic disease, but rarely in the context of endometriosis [8,13,16,17].
Endometriosis, defined by the presence of endometrium-like tissue outside the uterus, affects 2–10% of the female population, i.e., around 190 million women worldwide [18]. It is well known that endometriosis is a debilitating disease associated with severe symptoms. Consequently, endometriosis negatively affects all aspects of quality of life and is considered a public health issue related to its socioeconomic impact, and treatment and clinical management costs [19,20,21]. Previous studies have evaluated the potential of circulating miRNAs as biomarkers for endometriosis [22] and association with functions and pathophysiological pathways in endometriosis [8,15,17,23] but with conflicting results. This is due firstly to pre-analytical factors such as the source of the miRNAs (serum or plasma), type of blood collection tubes (EDTA versus heparin), hemolysis, and sample processing protocols. Secondly, there are technical factors related to the method used for RNA extraction, miRNA expression analysis (microarray, qRT-PCR and next-generation sequencing (NGS) techniques), and the strategy for normalization of miRNA expression data. Finally, biological factors are also implicated: the genetic background of the study cohort, the control population (self-reported healthy versus laparoscopically proven absence of endometriosis), and the extent of endometriosis (stage I versus stage IV) [8,15,17].
Therefore, the goal of the present study was to describe using NGS and bioinformatics systematic approach the miRNome sequencing of 200 plasma samples based on the prospective data from the ENDO-miRNA study.

2. Material and Methods

2.1. Study Population

We used data from the prospective “ENDO-miRNA” study (ClinicalTrials.gov Identifier: NCT04728152) [24]. Data collection and analysis (previously presented) were carried out under Research Protocol n° ID RCB: 2020-A03297-32 [25]. The IRB was delivered by the Comité de Protection des Personnes (C.P.P.) Sud-Ouest et Outre-Mer 1 (CPP 1-20-095 ID 10476. All patients gave informed written consent. The ENDO-miRNA study included 200 plasma samples obtained from patients with chronic pelvic pain suggestive of endometriosis. All had undergone a laparoscopic procedure (either operative or diagnostic) and/or magnetic resonance imaging (MRI) imaging proving endometriosis by the presence of endometrioma and/or deep endometriosis [26,27,28]. All laparoscopies were performed by two expert surgeons in endometriosis (ED, SB). For these patients, diagnosis of endometriosis was confirmed by histology. For the patients without laparoscopic evaluation, endometriosis was diagnosed when MRI revealed features of deep endometriosis with colorectal involvement and/or endometrioma confirmed by a multidisciplinary endometriosis committee. The study population was eventually composed of two groups: (i) endometriosis group composed of patients with endometriosis confirmed at either laparoscopy or MRI; (ii) control group without endometriosis at laparoscopy with and without other gynecological disorders. All patients included in the control group underwent a systematic laparoscopy. Among patients of the endometriosis group, 83 (54.2%) underwent an operative laparoscopy with histological confirmation of endometriosis and the remaining 70 (45.8%) had MRI confirmation [29]. The samples were collected from all the participants between January 2021 and June 2021. Statistical and miRNAs assays analysis were performed blinded to the surgical and imaging findings. The patients with endometriosis were stratified according to the revised American Society of Reproductive Medicine (rASRM) classification [30]. All patients filled online questionnaires to assess their symptoms and intensity using Visual Analogic Scale (VAS) [31].

2.2. Sample Collection

Blood samples (4 mL) were collected in EDTA tubes (BD, Franklin Lakes, NJ, USA). Plasma was isolated from whole blood within 2 h after blood sampling by two successive centrifugations at 4 °C (first at 1900× g (3000 rpm) for 10 min, followed by 13,000–14,000× g for 10 min to remove all cell debris) then aliquoted, labeled, and stored at −80 °C until analysis, as previously published [32,33,34].

2.3. RNA Sample Extraction, Preparation and Quality Control

RNA was extracted from 500 μL of plasma on a Maxwell 48® RSC automat using the Maxwell® RSC miRNA Plasma and Serum Kit (ref AS1680, Promega, USA) according to the manufacturer’s protocol. Libraries for small RNA sequencing were prepared using the QIAseq miRNA Library Kit for Illumina (Qiagen, Hilden, Germany). The resulting small RNA libraries were concentrated by ethanol precipitation and quantified using a Qubit 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). Samples were indexed in batches of 96, with a targeted sequencing depth of 17 million reads per sample. Sequencing was performed using 100 base single-end reads, using an Novaseq6000 sequencer (Illumina, San Diego, CA, USA) [35,36].

3. Bioinformatics

3.1. Raw Data Preprocessing (Raw, Filtered, Aligned Reads) and Quality Control

Sequencing reads were processed after adaptation using the bioinformatics and processing pipeline according to Potla et al., review [37,38]. FastQ files were trimmed to remove adapter sequences using Cutadapt version v.1.18 and were aligned using Bowtie version 1.1.1 to the following transcriptome databases: the human reference genome available from NCBI (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genome/guide/human/ 10 March 2022 and miRBase21) (miRNAs) using the MirDeep2 v0.1.0 package. The raw sequencing data quality was assessed using FastQC software v0.11.7 [10,14,35,39,40].

3.2. Differential Expression Analysis of miRNA

miRNA expression was quantified using miRDeep2 [41]. Differential expression tests were then conducted in DESeq2 for miRNAs with read counts in ≥1 of the samples. DESeq2 integrates methodological advances with several novel features to facilitate a more quantitative analysis of comparative RNA-seq data using shrinkage estimators for dispersion and fold change [41,42]. miRNAs were considered as differentially expressed if the absolute value of log2-fold change was >1.5 (up) and <0.5 (down) and the p value adjusted for multiple testing was <0.05 [41].

3.3. miRNome Accuracy

To evaluate the accuracy of each miRNA biomarker, sensitivity, specificity, and ROC analysis was performed, and the ROC AUC was calculated [43,44]. Additional statistical analysis was based on the Chi2 test as appropriate for categorical variables. Values of p < 0.05 were considered to denote significant differences. Data were managed with an Excel database (Microsoft, Redmond, WA, USA) and analyzed using R 2.15 software, available online (http://cran.r-project.org/, accessed on 10 March 2022).

4. Results

4.1. Description of the ENDO-miRNA Cohort

The clinical characteristics of the endometriosis and control patients are presented in Table 1. Among the 200 patients, 76.5% (n = 153) were diagnosed with endometriosis and 23.5% (n = 47) without. In the endometriosis group, 52% (80) had rASRM stages I–II and 48% (73) had stages III–IV. The control group was mainly composed of complex patients defined by patients sharing symptoms of endometriosis but without clinical or imaging features of endometriosis, and patients with other gynecologic disorders but with symptoms suggestive of endometriosis.

4.2. Global Overview of miRNA Transcriptome

Small RNA-seq of 200 plasma samples yielded ~4228 M raw sequencing reads (from ~11.7 M to ~34.98 M reads/sample). Pre-filtering and filtering steps retained 39% (~1639 M) of initial raw reads. The majority of filtered reads were of 20–23 nt length which corresponds to the range of mature miRNA sequences. Quantification of filtered reads and identification of known miRNAs yielded ~2588 M sequences to be mapped to 2633 known miRNAs from miRBase v22. The number of expressed miRNAs ranged from 666 to 1274 per sample. The distribution of expressed miRNAs in the 200 plasma samples and according to the overall composition of processed reads is shown in Figure 1A,B and Figure 2.

4.3. miRNA Expression in Patients with and without Endometriosis

A total of 2633 miRNAs were found differentially expressed in the plasma samples of patients with endometriosis compared with control patients. Among these, 8.6% (n = 229) were up- or downregulated. Respectively, 66% (152) and 34% (77) of the 229 miRNAs were up- and down regulated (Table S1). A volcano plot of the expressed miRNAs in the endometriosis patients is reported in Figure 3. Among the 152 miRNAs upregulated, only 5 (hsa-miR-29b-1-5p, hsa-miR-4748, hsa-miR-515-5p, hsa-miR-548j-5p, hsa-miR-6502-5p) had an AUC > 0.6. Among the 77 miRNAs downregulated, 2 (hsa-miR-3137, hsa-miR-3168) had an AUC > 0.6.

4.4. Diagnostic Accuracy of Regulated miRNAs

The diagnostic metrics for endometriosis for all the regulated plasma miRNAs (n = 229) are reported in Table S1. Among these 229 miRNAs, the F1-score, sensitivity, specificity, and AUC ranged from 0–88.2%, 0–99.4%, 4.3–100%, and 41.5–68%, respectively.
For AUC criteria, 96.9% (222) and 3.1% (7) had a value ranging between 41.5–59% and ≥60%, respectively.
For the F1-scores, 66.8% (n = 153) and 33.2% (n = 76) had a value ranging between 0–79%, and ≥80%, respectively.
For sensitivity, 69.8% (n = 160) and 30.1% (n = 69) had a value ranging between 0–79%, and ≥80%, respectively.
For specificity, 42% (n = 96) and 58% (n = 133) had a value ranging between 0–79%, and ≥80%, respectively.
Among the 229 regulated miRNAs, 69 had a sensitivity over 0.80 and 100 had a specificity over 0.80 but none had both sensitivity and specificity over 0.80 (Table S1).

5. Discussion

To our knowledge this is the first report which describes a sequencing and systematic bioinformatics approach for plasma miRNome of patients with endometriosis. The current study demonstrates that using an NGS technique allows the display of a specific and broad panel of miRNAs potentially suitable for statistical analysis as potential biomarker in endometriosis.
The recent literature review by Monnaka et al. of miRNA expression in endometriosis found that 30 miRNAs were deregulated in the blood; 27 in the serum, and 18 in the plasma of women with endometriosis compared with control populations. Interestingly, the authors of this critical review concluded that no particular miRNAs or miRNA combination was individually accurate enough to screen and diagnose endometriosis [17]. Similarly, Vanhie et al., reported that 42 miRNAs were deregulated in the blood samples of a biobank of patients with endometriosis, but failed to build a signature [15]. In addition, several miRNAs have been shown to be deregulated during the pathogenic process of endometriosis [8,45,46]. For example, Maged et al. have shown that serum miR-122 and miR-199a had a sensitivity of 95.6 and 100.0% and a specificity of 91.4 and 100%, respectively, for diagnosis of disease status in women [47]. All these data raise the question of the technical and methodological obstacles to identify miRNAs significantly associated with the presence of endometriosis. Among these obstacles, the technology used to investigate circulatory miRNAs appears crucial. Indeed, most previous studies evaluated miRNA profiles by microarray. Subsequently, only the most differentially expressed miRNAs in patients with endometriosis were validated by qRT-PCR which represents a major bias [8]. Moreover, fold change varies from one series to another, or is not mentioned, and AUC was not systematically reported [48,49]. In addition, few numbers of miRNA biomarkers have been studied in contrast with the large number of miRNAs associated with endometriosis assessed in the current study [17,22]. To overcome these concerns, in the current prospective study miRNAs were sequenced using NGS platforms, allowing the analysis of millions of RNA fragments. Moreover, unlike microarray, the NGS technique—recognized as one of the most efficient tools in this domain—does not require sequence specific hybridization probes [8]. Thanks to this technology, the 200 plasma samples we analyzed initially yielded ~4228 M raw sequencing reads (from ~11.7 M to ~34.98 M reads/sample) of which 39% (~1639 M) were retained after the pre-filtering and filtering steps. Among the 2633 miRNAs expressed in the plasma of patients with endometriosis compared to control patients, 229 were up- or downregulated. Of these, 30.1% had a sensitivity ≥80%, and 58% had a specificity ≥80%. The rate of miRNAs with an AUC value ≥ 60% was 3.1%. These values attest that high quality and high yields of transcriptomic miRNA information can be isolated from plasma without the need for preamplification. Another crucial goal is to evaluate the stability and the reproducibility of the miRNA reads for the 200 samples (100% of sample). In the current study, these two criteria were fulfilled: all 200 samples were used for sequencing, and bioinformatics treatment provided diagnostic accuracies according to the F1-score, sensitivity, specificity, and AUC, which ranged from 0–88.2%, 0–99.4%, 4.3–100%, and 41.5–68%, respectively.
Another issue is the difficulty of simultaneously analyzing hundreds of miRNAs taking into account the diversity of endometriosis phenotypes and the incomplete knowledge of the pathophysiology [22]. Previous studies have focused on the miRNAs involved in classic known signaling pathways of endometriosis including proliferation, apoptosis, cell differentiation, angiogenesis, inflammation, etc. [8,17,50]. While this approach is logical for a disorder with well identified pathophysiologic mechanisms, it is not suitable in the specific setting of endometriosis with incomplete knowledge of signaling pathways. In this specific issue, additional research are required to evaluate the relationship between miRNAs expression and genetic, epigenetic, and metabolomic abnormalities [51,52,53]. This suggests that there are issues concerning the methodology used to select the miRNAs, as well as the characteristics of the control group. As previously shown for cancer, to build a miRNA signature sufficiently stable to provide the same accuracy across different platforms implies the inclusion of a broad spectrum of miRNAs [14,54]. This is totally in keeping with previous studies [14,54,55,56]: it is illusory to attempt to reflect the heterogeneity of a multifactorial disorder such as endometriosis by a limited number of miRNAs. Therefore, beyond classis logistic regression, it is necessary to use both NGS and new mathematical tools such as Machine Learning (ML) and Artificial Intelligence (AI) as proved in cancer models [14,54,55,57,58].
Another issue is the discrepancy in miRNA expression from one series to another [8,15,17,23]. In addition to endogenous qRT-PCR controls, and the platforms used for miRNA analysis, these discrepancies could be related to various causes such as differences in study design, patient population, sample size, and stage of endometriosis, but mainly to the composition of the control group. As underlined by Agrawal et al., choosing an appropriate control group is challenging and crucial to study miRNA expression. To limit the risk of bias, patients with pelvic inflammatory disease or autoimmune disorders, as well as healthy women that may have asymptomatic endometriosis (unless ruled out by laparoscopic evaluation) should be excluded [59]. We specifically designed a prospective study, including 200 patients—the largest series to date—able to quantify the miRNome for (i) complex patients (women with pelvic pain suggestive of endometriosis and both negative clinical and imaging findings with and without endometriosis at diagnostic laparoscopy), (ii) women with early-stage endometriosis (stage I-II rASRM), (iii) women with advanced stages (stage III-IV rASRM), (iv) women with various phenotypes of endometriosis (peritoneal endometriosis, ovarian endometriosis also called endometrioma, deep endometriosis defined by infiltration of pelvic organs and/or pelvic anatomical structures with and without endometrioma), and (v) women with other gynecologic disorders sharing symptoms of endometriosis.
Some remaining issues for the evaluation of miRNAs in the context of endometriosis should be discussed. Two of the important determinants of miRNA analysis are the phase of the menstrual cycle and the impact of hormonal treatments [60]. Although several studies have observed differences in miRNA expression in tissues according to the menstrual phase mainly at endometrial level [8,61,62,63], no such cyclic differences were observed in the plasma of healthy women [49]. One hypothesis is that changes in miRNA expression at the endometrium level regulate gene expression locally but are insufficient to cause detectable systemic changes [8]. Moreover, previous studies report no specific impact of hormonal treatment on blood miRNA [8,15]. Another concern is that among the sncRNAs, we only focused on miRNA while it is possible that other RNAs such as siRNAs, piRNAs, and snoRNAs as well as lncRNAs could be of diagnostic and therapeutic values.

6. Future Perspectives

It is well known that endometriosis, which affects 190 million women worldwide, is a debilitating disorder with a negative impact on quality of life and fertility, and that it represents a major socioeconomic burden [64]. Our results provide further evidence that patients with endometriosis exhibit a specific panel of miRNAs potentially suitable as biomarker. Beyond the specific context of endometriosis, the methodology developed in the current study can be transposable to other benign chronic and malignant diseases.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/diagnostics12051150/s1, Table S1: Performances of the 229 differentially expressed in the cohort.

Author Contributions

Methodology and Design: S.B., P.D., F.G. and E.D.; Data collection: S.B., Y.D., A.P., L.D., M.P., C.T., F.G., P.D. and E.D.; Biologic data collection: Y.D., L.J. and D.B.; Analysis: S.S., S.B. and E.D.; Data Interpretation: S.B., S.S., Y.D., L.D., M.P., A.P., C.T. and E.D.; Manuscript writing: Y.D., S.B. and E.D.; All authors reviewed the manuscript for critical intellectual content. All authors have read and agreed to the published version of the manuscript.

Funding

Part of this work was funded by a grant from the Conseil Régional d’Ile de France (grant number EX024087) and from Ziwig, Inc.

Institutional Review Board Statement

The authors state that the data used are from the prospective ENDO-miRNA study (ClinicalTrials.gov Identifier: NCT04728152). Data collection and analysis were carried out under Research Protocol n° ID RCB: 2020-A03297-32. The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board delivered by the Comité de Protection des Personnes (C.P.P.) Sud-Ouest et Outre-Mer 1 (CPP 1-20-095 ID 10476.

Informed Consent Statement

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

Data Availability Statement

All relevant data are within the manuscript or Supplementary Files.

Acknowledgments

All authors would like to sincerely thank F. Neilson for her English revision of the manuscript.

Conflicts of Interest

S. Suisse is a former employee of Ziwig, Inc. The remaining authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript.

References

  1. Shahrouki, P.; Larsson, E. The Non-Coding Oncogene: A Case of Missing DNA Evidence? Front. Genet. 2012, 3, 170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Green, D.; Fraser, W.D.; Dalmay, T. Transfer RNA-Derived Small RNAs in the Cancer Transcriptome. Pflugers Arch. 2016, 468, 1041–1047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gibb, E.A.; Brown, C.J.; Lam, W.L. The Functional Role of Long Non-Coding RNA in Human Carcinomas. Mol. Cancer 2011, 10, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Doyle, L.M.; Wang, M.Z. Overview of Extracellular Vesicles, Their Origin, Composition, Purpose, and Methods for Exosome Isolation and Analysis. Cells 2019, 8, 727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Yuan, T.; Huang, X.; Woodcock, M.; Du, M.; Dittmar, R.; Wang, Y.; Tsai, S.; Kohli, M.; Boardman, L.; Patel, T.; et al. Plasma Extracellular RNA Profiles in Healthy and Cancer Patients. Sci. Rep. 2016, 6, 19413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bartel, D.P. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  7. Bartel, D.P. MicroRNAs: Target Recognition and Regulatory Functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [Green Version]
  8. Agrawal, S.; Tapmeier, T.; Rahmioglu, N.; Kirtley, S.; Zondervan, K.; Becker, C. The MiRNA Mirage: How Close Are We to Finding a Non-Invasive Diagnostic Biomarker in Endometriosis? A Systematic Review. Int. J. Mol. Sci. 2018, 19, 599. [Google Scholar] [CrossRef] [Green Version]
  9. Griffiths-Jones, S.; Grocock, R.J.; van Dongen, S.; Bateman, A.; Enright, A.J. MiRBase: MicroRNA Sequences, Targets and Gene Nomenclature. Nucleic Acids Res. 2006, 34, D140–D144. [Google Scholar] [CrossRef]
  10. Griffiths-Jones, S.; Saini, H.K.; van Dongen, S.; Enright, A.J. MiRBase: Tools for MicroRNA Genomics. Nucleic Acids Res. 2008, 36, D154–D158. [Google Scholar] [CrossRef] [Green Version]
  11. Griffiths-Jones, S. The MicroRNA Registry. Nucleic Acids Res. 2004, 32, D109–D111. [Google Scholar] [CrossRef] [PubMed]
  12. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. MiRBase: From MicroRNA Sequences to Function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef]
  13. Kang, J.-W.; Eun, Y.-G.; Lee, Y.-C. Diagnostic Value of Salivary MiRNA in Head and Neck Squamous Cell Cancer: Systematic Review and Meta-Analysis. Int. J. Mol. Sci. 2021, 22, 7026. [Google Scholar] [CrossRef] [PubMed]
  14. Lopez-Rincon, A.; Mendoza-Maldonado, L.; Martinez-Archundia, M.; Schönhuth, A.; Kraneveld, A.D.; Garssen, J.; Tonda, A. Machine Learning-Based Ensemble Recursive Feature Selection of Circulating MiRNAs for Cancer Tumor Classification. Cancers 2020, 12, 1785. [Google Scholar] [CrossRef] [PubMed]
  15. Vanhie, A.; O, D.; Peterse, D.; Beckers, A.; Cuéllar, A.; Fassbender, A.; Meuleman, C.; Mestdagh, P.; D’Hooghe, T. Plasma MiRNAs as Biomarkers for Endometriosis. Hum. Reprod. 2019, 34, 1650–1660. [Google Scholar] [CrossRef]
  16. Courts, C.; Madea, B. Specific Micro-RNA Signatures for the Detection of Saliva and Blood in Forensic Body-Fluid Identification. J. Forensic Sci. 2011, 56, 1464–1470. [Google Scholar] [CrossRef]
  17. Monnaka, V.U.; Hernandes, C.; Heller, D.; Podgaec, S. Overview of MiRNAs for the Non-Invasive Diagnosis of Endometriosis: Evidence, Challenges and Strategies. A Systematic Review. Einstein 2021, 19, eRW5704. [Google Scholar] [CrossRef]
  18. Zondervan, K.T.; Becker, C.M.; Missmer, S.A. Endometriosis. N. Engl. J. Med. 2020, 382, 1244–1256. [Google Scholar] [CrossRef]
  19. Zondervan, K.T.; Becker, C.M.; Koga, K.; Missmer, S.A.; Taylor, R.N.; Viganò, P. Endometriosis. Nat. Rev. Dis. Primers 2018, 4, 9. [Google Scholar] [CrossRef]
  20. As-Sanie, S.; Soliman, A.M.; Evans, K.; Erpelding, N.; Lanier, R.; Katz, N.P. Healthcare Utilization and Cost Burden among Women with Endometriosis by Opioid Prescription Status in the First Year after Diagnosis: A Retrospective Claims Database Analysis. J. Med. Econ. 2020, 23, 371–377. [Google Scholar] [CrossRef]
  21. Surrey, E.; Soliman, A.M.; Trenz, H.; Blauer-Peterson, C.; Sluis, A. Impact of Endometriosis Diagnostic Delays on Healthcare Resource Utilization and Costs. Adv. Ther. 2020, 37, 1087–1099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Dabi, Y.; Suisse, S.; Jornea, L.; Bouteiller, D.; Touboul, C.; Puchar, A.; Daraï, E.; Bendifallah, S. Clues for Improving the Pathophysiology Knowledge for Endometriosis Using Serum Micro-RNA Expression. Diagnostics 2022, 12, 175. [Google Scholar] [CrossRef] [PubMed]
  23. Moustafa, S.; Burn, M.; Mamillapalli, R.; Nematian, S.; Flores, V.; Taylor, H.S. Accurate Diagnosis of Endometriosis Using Serum MicroRNAs. Am. J. Obstet. Gynecol. 2020, 223, e1–e557. [Google Scholar] [CrossRef]
  24. Bendifallah, S. Evaluation of MiRNAs in Endometriosis. 2021. Available online: https://clinicaltrials.gov/ct2/show/NCT04728152 (accessed on 25 April 2022).
  25. Bendifallah, S.; Dabi, Y.; Suisse, S.; Jornea, L.; Bouteiller, D.; Touboul, C.; Puchar, A.; Daraï, E. MicroRNome Analysis Generates a Blood-Based Signature for Endometriosis. Sci. Rep. 2022, 12, 4051. [Google Scholar] [CrossRef] [PubMed]
  26. Ito, T.E.; Abi Khalil, E.D.; Taffel, M.; Moawad, G.N. Magnetic Resonance Imaging Correlation to Intraoperative Findings of Deeply Infiltrative Endometriosis. Fertil. Steril. 2017, 107, e11–e12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bazot, M.; Daraï, E. Diagnosis of Deep Endometriosis: Clinical Examination, Ultrasonography, Magnetic Resonance Imaging, and Other Techniques. Fertil. Steril. 2017, 108, 886–894. [Google Scholar] [CrossRef] [Green Version]
  28. Bazot, M.; Darai, E.; Hourani, R.; Thomassin, I.; Cortez, A.; Uzan, S.; Buy, J.-N. Deep Pelvic Endometriosis: MR Imaging for Diagnosis and Prediction of Extension of Disease. Radiology 2004, 232, 379–389. [Google Scholar] [CrossRef] [Green Version]
  29. Bendifallah, S.; Suisse, S.; Puchar, A.; Delbos, L.; Poilblanc, M.; Descamps, P.; Golfier, F.; Jornea, L.; Bouteiller, D.; Touboul, C.; et al. Salivary MicroRNA Signature for Diagnosis of Endometriosis. J. Clin. Med. 2022, 11, 612. [Google Scholar] [CrossRef]
  30. Canis, M.; Donnez, J.G.; Guzick, D.S.; Halme, J.K.; Rock, J.A.; Schenken, R.S.; Vernon, M.W. Revised American Society for Reproductive Medicine Classification of Endometriosis: 1996. Fertil. Steril. 1997, 67, 817–821. [Google Scholar] [CrossRef]
  31. Bendifallah, S.; Puchar, A.; Suisse, S.; Delbos, L.; Poilblanc, M.; Descamps, P.; Golfier, F.; Touboul, C.; Dabi, Y.; Daraï, E. Machine Learning Algorithms as New Screening Approach for Patients with Endometriosis. Sci. Rep. 2022, 12, 639. [Google Scholar] [CrossRef]
  32. Canlorbe, G.; Castela, M.; Bendifallah, S.; Wang, Z.; Lefevre, M.; Chabbert-Buffet, N.; Aractingi, S.; DaraÏ, E.; Méhats, C.; Ballester, M. Micro-RNA Signature of Lymphovascular Space Involvement in Type 1 Endometrial Cancer. Histol. Histopathol. 2017, 32, 941–950. [Google Scholar] [CrossRef] [PubMed]
  33. Canlorbe, G.; Wang, Z.; Laas, E.; Bendifallah, S.; Castela, M.; Lefevre, M.; Chabbert-Buffet, N.; Daraï, E.; Aractingi, S.; Méhats, C.; et al. Identification of MicroRNA Expression Profile Related to Lymph Node Status in Women with Early-Stage Grade 1-2 Endometrial Cancer. Mod. Pathol. 2016, 29, 391–401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. de Foucher, T.; Sbeih, M.; Uzan, J.; Bendifallah, S.; Lefevre, M.; Chabbert-Buffet, N.; Aractingi, S.; Uzan, C.; Abd Alsalam, I.; Mitri, R.; et al. Identification of Micro-RNA Expression Profile Related to Recurrence in Women with ESMO Low-Risk Endometrial Cancer. J. Transl. Med. 2018, 16, 131. [Google Scholar] [CrossRef] [PubMed]
  35. Gyvyte, U.; Juzenas, S.; Salteniene, V.; Kupcinskas, J.; Poskiene, L.; Kucinskas, L.; Jarmalaite, S.; Stuopelyte, K.; Steponaitiene, R.; Hemmrich-Stanisak, G.; et al. MiRNA Profiling of Gastrointestinal Stromal Tumors by Next-Generation Sequencing. Oncotarget 2017, 8, 37225–37238. [Google Scholar] [CrossRef]
  36. Gyvyte, U.; Kupcinskas, J.; Juzenas, S.; Inciuraite, R.; Poskiene, L.; Salteniene, V.; Link, A.; Fassan, M.; Franke, A.; Kupcinskas, L.; et al. Identification of Long Intergenic Non-Coding RNAs (LincRNAs) Deregulated in Gastrointestinal Stromal Tumors (GISTs). PLoS ONE 2018, 13, e0209342. [Google Scholar] [CrossRef]
  37. Potla, P.; Ali, S.A.; Kapoor, M. A Bioinformatics Approach to MicroRNA-Sequencing Analysis. Osteoarthr. Cartil. Open 2021, 3, 100131. [Google Scholar] [CrossRef]
  38. Li, Q.; Liu, G.; Bao, Y.; Wu, Y.; You, Q. Evaluation and Application of Tools for the Identification of Known MicroRNAs in Plants. Appl. Plant Sci. 2021, 9, e11414. [Google Scholar] [CrossRef]
  39. 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, R25. [Google Scholar] [CrossRef] [Green Version]
  40. Gao, L.; Zhang, L. Construction and Comprehensive Analysis of a CeRNA Network to Reveal Potential Prognostic Biomarkers for Lung Adenocarcinoma. BMC Cancer 2021, 21, 849. [Google Scholar] [CrossRef]
  41. 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] [Green Version]
  42. Bargaje, R.; Hariharan, M.; Scaria, V.; Pillai, B. Consensus MiRNA Expression Profiles Derived from Interplatform Normalization of Microarray Data. RNA 2010, 16, 16–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Harrell, F.E.; Lee, K.L.; Mark, D.B. Multivariable Prognostic Models: Issues in Developing Models, Evaluating Assumptions and Adequacy, and Measuring and Reducing Errors. Stat. Med. 1996, 15, 361–387. [Google Scholar] [CrossRef]
  44. Steyerberg, E.W.; Eijkemans, M.J.; Harrell, F.E.; Habbema, J.D. Prognostic Modelling with Logistic Regression Analysis: A Comparison of Selection and Estimation Methods in Small Data Sets. Stat. Med. 2000, 19, 1059–1079. [Google Scholar] [CrossRef]
  45. Panir, K.; Schjenken, J.E.; Robertson, S.A.; Hull, M.L. Non-Coding RNAs in Endometriosis: A Narrative Review. Hum. Reprod. Update 2018, 24, 497–515. [Google Scholar] [CrossRef] [Green Version]
  46. Ghafouri-Fard, S.; Shoorei, H.; Taheri, M. Role of Non-Coding RNAs in the Pathogenesis of Endometriosis. Front. Oncol. 2020, 10, 1370. [Google Scholar] [CrossRef]
  47. Maged, A.M.; Deeb, W.S.; El Amir, A.; Zaki, S.S.; El Sawah, H.; Al Mohamady, M.; Metwally, A.A.; Katta, M.A. Diagnostic accuracy of serum miR-122 and miR-199a in women with endometriosis. Int. J. Gynaecol. Obstet. 2018, 141, 14–19. [Google Scholar] [CrossRef]
  48. Wang, W.-T.; Zhao, Y.-N.; Han, B.-W.; Hong, S.-J.; Chen, Y.-Q. Circulating MicroRNAs Identified in a Genome-Wide Serum MicroRNA Expression Analysis as Noninvasive Biomarkers for Endometriosis. J. Clin. Endocrinol. Metab. 2013, 98, 281–289. [Google Scholar] [CrossRef] [Green Version]
  49. Rekker, K.; Saare, M.; Roost, A.M.; Kaart, T.; Sõritsa, D.; Karro, H.; Sõritsa, A.; Simón, C.; Salumets, A.; Peters, M. Circulating MiR-200-Family Micro-RNAs Have Altered Plasma Levels in Patients with Endometriosis and Vary with Blood Collection Time. Fertil. Steril. 2015, 104, 938–946.e2. [Google Scholar] [CrossRef]
  50. Fernandes, J.; Acuña, S.; Aoki, J.; Floeter-Winter, L.; Muxel, S. Long Non-Coding RNAs in the Regulation of Gene Expression: Physiology and Disease. ncRNA 2019, 5, 17. [Google Scholar] [CrossRef] [Green Version]
  51. Zubrzycka, A.; Zubrzycki, M.; Perdas, E.; Zubrzycka, M. Genetic, Epigenetic, and Steroidogenic Modulation Mechanisms in Endometriosis. J. Clin. Med. 2020, 9, 1309. [Google Scholar] [CrossRef]
  52. Angioni, S.; D’Alterio, M.N.; Coiana, A.; Anni, F.; Gessa, S.; Deiana, D. Genetic Characterization of Endometriosis Patients: Review of the Literature and a Prospective Cohort Study on a Mediterranean Population. Int. J. Mol. Sci. 2020, 21, 1765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Deiana, D.; Gessa, S.; Anardu, M.; Daniilidis, A.; Nappi, L.; D’Alterio, M.N.; Pontis, A.; Angioni, S. Genetics of Endometriosis: A Comprehensive Review. Gynecol. Endocrinol. 2019, 35, 553–558. [Google Scholar] [CrossRef] [PubMed]
  54. Lopez-Rincon, A.; Martinez-Archundia, M.; Martinez-Ruiz, G.U.; Schoenhuth, A.; Tonda, A. Automatic Discovery of 100-MiRNA Signature for Cancer Classification Using Ensemble Feature Selection. BMC Bioinform. 2019, 20, 480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Akter, S.; Xu, D.; Nagel, S.C.; Bromfield, J.J.; Pelch, K.E.; Wilshire, G.B.; Joshi, T. GenomeForest: An Ensemble Machine Learning Classifier for Endometriosis. AMIA Jt. Summits Transl. Sci. Proc. 2020, 2020, 33–42. [Google Scholar] [PubMed]
  56. Jang, H.-J.; Lee, H.-S.; Burt, B.M.; Lee, G.K.; Yoon, K.-A.; Park, Y.-Y.; Sohn, B.H.; Kim, S.B.; Kim, M.S.; Lee, J.M.; et al. Integrated Genomic Analysis of Recurrence-Associated Small Non-Coding RNAs in Oesophageal Cancer. Gut 2017, 66, 215–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Akter, S.; Xu, D.; Nagel, S.C.; Bromfield, J.J.; Pelch, K.; Wilshire, G.B.; Joshi, T. Machine Learning Classifiers for Endometriosis Using Transcriptomics and Methylomics Data. Front. Genet. 2019, 10, 766. [Google Scholar] [CrossRef] [Green Version]
  58. Kleczyk, E.J.; Peri, A.; Yadav, T.; Komera, R.; Peri, M.; Guduru, V.; Amirtharaj, S.; Huang, M. Predicting Endometriosis Onset Using Machine Learning Algorithms. 2020; in review. [Google Scholar]
  59. Zondervan, K.T.; Cardon, L.R.; Kennedy, S.H. What Makes a Good Case-Control Study? Design Issues for Complex Traits Such as Endometriosis. Hum. Reprod. 2002, 17, 1415–1423. [Google Scholar] [CrossRef] [Green Version]
  60. Kolanska, K.; Bendifallah, S.; Canlorbe, G.; Mekinian, A.; Touboul, C.; Aractingi, S.; Chabbert-Buffet, N.; Daraï, E. Role of MiRNAs in Normal Endometrium and in Endometrial Disorders: Comprehensive Review. J. Clin. Med. 2021, 10, 3457. [Google Scholar] [CrossRef]
  61. Kuokkanen, S.; Chen, B.; Ojalvo, L.; Benard, L.; Santoro, N.; Pollard, J.W. Genomic Profiling of MicroRNAs and Messenger RNAs Reveals Hormonal Regulation in MicroRNA Expression in Human Endometrium. Biol. Reprod. 2010, 82, 791–801. [Google Scholar] [CrossRef] [Green Version]
  62. Sha, A.-G.; Liu, J.-L.; Jiang, X.-M.; Ren, J.-Z.; Ma, C.-H.; Lei, W.; Su, R.-W.; Yang, Z.-M. Genome-Wide Identification of Micro-Ribonucleic Acids Associated with Human Endometrial Receptivity in Natural and Stimulated Cycles by Deep Sequencing. Fertil. Steril. 2011, 96, 150–155.e5. [Google Scholar] [CrossRef]
  63. Altmäe, S.; Martinez-Conejero, J.A.; Esteban, F.J.; Ruiz-Alonso, M.; Stavreus-Evers, A.; Horcajadas, J.A.; Salumets, A. MicroRNAs MiR-30b, MiR-30d, and MiR-494 Regulate Human Endometrial Receptivity. Reprod. Sci. 2013, 20, 308–317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Shafrir, A.L.; Farland, L.V.; Shah, D.K.; Harris, H.R.; Kvaskoff, M.; Zondervan, K.; Missmer, S.A. Risk for and Consequences of Endometriosis: A Critical Epidemiologic Review. Best Pract. Res. Clin. Obstet. Gynaecol. 2018, 51, 1–15. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (A) Distribution of expressed miRNAs in the 200 blood samples. (B) Distribution of expressed miRNAs in the samples diagnosis.
Figure 1. (A) Distribution of expressed miRNAs in the 200 blood samples. (B) Distribution of expressed miRNAs in the samples diagnosis.
Diagnostics 12 01150 g001
Figure 2. Overall composition of processed reads for plasma samples. In green: quantified reads; In Yellow: Unmapped reads; In Red: Filtered reads.
Figure 2. Overall composition of processed reads for plasma samples. In green: quantified reads; In Yellow: Unmapped reads; In Red: Filtered reads.
Diagnostics 12 01150 g002
Figure 3. Volcano plot of expressed miRNAs in plasma for endometriosis.
Figure 3. Volcano plot of expressed miRNAs in plasma for endometriosis.
Diagnostics 12 01150 g003
Table 1. Main characteristics of the patients included.
Table 1. Main characteristics of the patients included.
Controls
N = 47
Endometriosis
N = 153
p Value
Age years (mean ± SD)30.92 ± 13.7931.17 ± 10.780.19
BMI (body mass index) (mean ± SD)24.84 ± 11.1024.36 ± 8.380.52
rASRM classification
- I–II
-80 (52%)
- III–IV
-73 (48%)
Control diagnoses
- No abnormality
24 (51%)--
- Leiomyoma
1 (2%)
- Cystadenoma
5 (11%)
- Teratoma
11 (23%)
- Others gynecological disorders
6 (13%)
Dysmenorrhea100%100%
Abdominal pain outside menstruation
- Yes
21 (66%)89 (71%)0.69
Patients with pain suggesting sciatica10 (31%)70 (56%)0.02
Dyspareunia intensity at VAS (mean ± SD)4.95 ± 3.525.28 ± 3.95<0.001
Patients with lower back pain outside menstruation20 (62%)101 (81%)0.049
Intensity of pain during defecation at VAS (mean ± SD)2.84 ± 2.764.35 ± 3.47<0.001
Patient with right shoulder pain during menstruation3 (9%)26 (21%)0.21
Intensity of urinary pain during menstruation at VAS (mean ± SD)2.84 ± 2.764.35 ± 3.36<0.001
Patient with blood in the stools during menstruation4 (12%)30 (24%)0.24
Patient with blood in urine during menstruation8 (25%)21 (17%)0.41
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bendifallah, S.; Dabi, Y.; Suisse, S.; Delbos, L.; Poilblanc, M.; Descamps, P.; Golfier, F.; Jornea, L.; Bouteiller, D.; Touboul, C.; et al. Endometriosis Associated-miRNome Analysis of Blood Samples: A Prospective Study. Diagnostics 2022, 12, 1150. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics12051150

AMA Style

Bendifallah S, Dabi Y, Suisse S, Delbos L, Poilblanc M, Descamps P, Golfier F, Jornea L, Bouteiller D, Touboul C, et al. Endometriosis Associated-miRNome Analysis of Blood Samples: A Prospective Study. Diagnostics. 2022; 12(5):1150. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics12051150

Chicago/Turabian Style

Bendifallah, Sofiane, Yohann Dabi, Stéphane Suisse, Léa Delbos, Mathieu Poilblanc, Philippe Descamps, Francois Golfier, Ludmila Jornea, Delphine Bouteiller, Cyril Touboul, and et al. 2022. "Endometriosis Associated-miRNome Analysis of Blood Samples: A Prospective Study" Diagnostics 12, no. 5: 1150. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics12051150

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