Profiling miRNA in salivary samples from subjects with endometriosis: a pilot study

Introduction: Endometriosis is a common, chronic, estrogen-dependent gynecological disorder characterized by the growth of endometrial-like tissue outside the uterus and frequently associated with pelvic pain and infertility. Despite its high prevalence, the molecular mechanisms underlying lesion persistence and inflammation remain poorly understood, limiting the development of reliable non-surgical diagnostic tools and targeted therapies.

A case-control study included 176 women with endometriosis and 124 controls. Salivary microRNA (miRNA) expression was analyzed using next-generation sequencing. Patients were stratified into untreated and treated groups (pharmacological or surgical therapy), while controls included healthy individuals and a technical control group of women with benign gynecological conditions to reduce potential confounding factors.

Ten salivary miRNAs were associated with disease activity. Seven (hsa-miR-130a-3p, hsa-miR-130b-3p, hsa-miR-141-3p, hsa-miR-200b-3p, hsa-miR-200c-3p, hsa-miR-203b-5p, and hsa-miR-29c-3p) had previously been linked to endometriosis, while three were associated with inflammatory or cancer-related pathways. Expression levels showed a gradient across clinical groups, with intermediate levels in treated patients, suggesting therapeutic modulation. Target analysis identified PTEN as a key regulated gene, implicating the PI3K/AKT/mTOR pathway.

These findings support salivary miRNAs as promising non-invasive biomarkers for the diagnosis and monitoring of endometriosis and provide further insight into the molecular pathways involved in disease pathogenesis. 1 Introduction Endometriosis is a common, chronic, and estrogen-dependent gynecological condition defined by the presence of endometrial-like tissue (glands and stroma) in ectopic locations (). The lesions are primarily found in the pelvic peritoneum, ovaries, and rectovaginal septum (). In light of these widespread effects, endometriosis should be considered a public health issue rather than merely a disease of individuals. Endometriosis can be categorized into three well-recognized phenotypes: superficial peritoneal lesions (SUP), ovarian endometriomas (OMA), and deep infiltrating endometriosis (DIE). OMA are cystic masses derived from ectopic endometrial tissue within the ovary, while the most severe phenotype, DIE, is defined as sub-peritoneal lesions that penetrate tissue deeper than 5 mm under the peritoneal surface, potentially infiltrating organs surrounding the uterus, such as bladder or intestine (; ), leading to its characterization as an “abdominal-pelvic multifocal disease” (). Furthermore, the presence of OMA is often an indicator of more severe associated DIE (Anonymous, 1985; Anonymous, 1997; ) and independently contributes to pain (Vercellini et al., 2014a; Naftalin et al., 2014). The disease is characterized by a wide variability in presentation and symptoms, which include dysmenorrhea, dyspareunia, chronic pelvic pain, irregular uterine bleeding, and/or infertility (). Notably, infertility affects 30%–50% of women diagnosed with endometriosis (). The most well-accepted pathophysiological hypothesis for endometriosis is based on retrograde menstruation, where menses transport viable endometrial fragments through the fallopian tubes to the peritoneal cavity for implantation (; Vercellini et al., 2014b). Genetically, endometriosis is considered a complex trait that exhibits familial aggregation, with up to a six-fold increased risk for first-degree relatives of patients (Stefansson et al., 2002; Treloa et al., 1999; Rahmioglu et al., 2012; Treloar et al., 2002), underscoring the complex, multi-factorial nature of the disease. The natural course of endometriosis is characterized by cyclic bleeding () and repeated tissue injury and repair at the site of the ectopic lesions, resulting in local inflammation (; ; Zhang et al., 2016). Endometriosis-related pain is mediated by several mechanisms. Ectopic endometrial cells outside the uterus stimulate immune cell infiltration, which secretes inflammatory mediators (such as cytokines and nerve growth factor), fostering a proinflammatory peritoneal microenvironment (). The close topographical relationship between endometriotic foci and nerves (; ; ) as well as structural changes within the central pain system (; ; ; Zhang et al., 2010) all contribute to pelvic pain (). Endometriosis-related infertility is a clear association, although the diagnosis is not synonymous with infertility. Adverse effects on fecundity include chronic inflammation of the peritoneal fluid (; ), disruption of ovarian function, pelvic adhesions, and reduced frequency of sexual intercourse due to dyspareunia (Raffi et al., 2012). Crucially, the infertility associated with OMA may be linked to the negative effect of surgical excision on the ovarian reserve (Raffi et al., 2012; Streuli et al., 2012; ) rather than the OMA lesions themselves (Santulli et al., 2016; ). Despite decades of research, the pathogenesis of endometriosis remains controversial and poorly understood. Although various factors contribute to its etiology, genetic predisposition accounts for only approximately 7% of cases (). Furthermore, the combination of several factors including the non-specific nature of the symptoms and inadequate public awareness contributes to a substantial diagnostic delay (). As a result, the typical delay period between the onset of symptoms and diagnosis has been estimated to be as long as 8 years, underlying the need for new biomarkers allowing an early and accurate diagnosis. In recent years, microRNAs (miRNAs or miRs) have emerged as powerful post-transcriptional regulators in disease pathogenesis. MiRNAs are small, endogenous non-coding RNA molecules (18–25 nucleotides) that regulate gene expression by binding to target messenger RNAs (mRNAs), leading to the repression of protein synthesis or degradation of the targeted mRNA (). A significant body of research has demonstrated a dysregulated miRNA profile in the eutopic and ectopic endometrium of women with endometriosis (; ), suggesting their critical involvement in altered endometrial biology, the inflammatory cascade, and the fibrotic changes characteristic of the disease. Specifically, miRNAs are implicated in regulating key pathways that contribute to the disease, such as impaired steroid biosynthesis, increased invasive potential associated with neoangiogenesis, and the proinflammatory profile observed in endometrial tissue (). Evidence further suggests that miRNAs can be released from damaged or apoptotic cells (), potentially constituting a specific circulating miRNA signature reflective of the disease burden. However, their precise role as gene drivers in the overall pathogenesis is still being investigated. To date, most studies investigating miRNA-based biomarkers for endometriosis have focused on blood-derived samples, including plasma and serum (Vanhie et al., 2019; ). These studies have identified several candidate miRNAs associated with disease presence and severity, highlighting the potential of circulating miRNAs as non-invasive diagnostic tools (Vanhie et al., 2024; Suryawanshi et al., 2013; Wang et al., 2013; ). However, blood-based miRNA analysis presents some limitations, including the invasive nature of sample collection, variability associated with menstrual cycle phase, systemic inflammatory conditions, haemolysis, and other physiological factors that may influence circulating miRNA profiles. These limitations have contributed to the ongoing search for alternative biological matrices suitable for biomarker discovery. Among these, saliva has emerged as an attractive source of diagnostic biomarkers. Saliva collection is simple, safe, non-invasive, and readily repeatable, making it particularly suitable for longitudinal monitoring. In addition, salivary miRNAs exhibit remarkable stability because they are protected within extracellular vesicles or associated with RNA-binding proteins. Accumulating evidence further suggests that salivary miRNAs can reflect systemic physiological and pathological processes rather than exclusively local oral conditions, supporting their use as biomarkers for diseases affecting distant organs. Based on these considerations, in the present study we investigated the miRNAome in non-invasive samples, represented by saliva obtained by healthy donors (without any gynecological issues) and compared it with that of patients with a confirmed diagnosis of endometriotic lesions or diagnosed with benign gynecological conditions, such adenomas or myomas with the final aim of finding a panel of miRNA with a non-invasive diagnostic potential. We also examined the possible modulation of the miRNAome as a consequence of a therapeutic intervention. Once the cluster with diagnostic potential was identified, we further investigated the biological role of the miRNA forming the cluster in order to understand its biological significance and specificity in the context of endometriosis, in the future perspective of facilitating the diagnosis of such debilitating condition affecting women worldwide. 2 Materials and methods 2.1 Ethic statement All procedures in this study were conducted in accordance with the Helsinki Declaration. The study was approved by the Regional Ethics Committee for Clinical Trials of the Region of Tuscany Area Vasta Sud EST (ID23576), by Lombardy Regional Ethics Committee 3 (ID 4165_S_N). The study has also been submitted and registered on ClinicalTrials.gov (ID protocol number: NCT06100471 registration date 10-13-2023). All women enrolled in the study were properly informed and signed an informed consent. 2.2 Study cohort and sample collection The case-control study was conducted considering two patients’ cohorts composed by300 patients. The case group consisted of 176 patients diagnosed with endometriosis, while the control group included a total of 124 subjects belonging to two different sub-groups: healthy (no symptoms) or patients with other gynecological conditions (Figure 1). FIGURE 1 Endometriosis was diagnosed either clinically, supported by imaging techniques, or through histological confirmation following surgical intervention. Among 300 collected samples, 228/300 (76%) were eligible for miRNA expression profiling analysis, while 72/300 (24%) were classified inconclusive for eligibility criteria, or not suitable considering pre-analytical criteria in sampling or timeline shipping. Saliva samples were collected by using OMNIgene kit (Danagen-Bioted, S.L., Badalona, Barcellona-Spain) following the manufacturer’s instructions. Salivary sample were obtained by subjects avoiding eating, drinking, or smoking at least 30 min before collection, stored at 4 °C, and shipped within a maximum of 5 days. In order to improve the interpretability of biomarker data also reducing potential confounding factors, endometriosis cases were further stratified based on their therapeutic status at the time of sample collection. Specifically, patients were further categorized into: Treated cases: patients with a diagnosis of endometriosis undergoing pharmacological therapy (e.g., hormonal treatment) or who had recently undergone surgical intervention. Untreated cases: patients with a first diagnosis of endometriosis not receiving at the time of the diagnosis any pharmacological treatment or surgical procedures. The control group was also stratified as follow: Technical controls: patients with benign gynecological conditions (e.g., ovarian cysts, fibroids) but without a diagnosis of endometriosis. Healthy controls: patients with no gynecological issues, referring to a visit as yearly routine control. 2.3 miRNA extraction and quality control Saliva samples were immediately processed or stored at 4 °C until use. All processes were conducted in an RNase-free environment. Both RNA extraction and miRNA-enrichment were performed by using 1 mL of salivary sample with the miRNeasy Mini Kit and MyQIAcube, (Qiagen, Germantown, MD, USA) according to the manufacturer’s instructions. Briefly, 1 mL of saliva was centrifuged at 3000 RCF for 20 min, the supernatant was removed, and the pellet was resuspended in 700 µL of QIAzol Lysis Reagent (Qiagen). Then, 140 µL of chloroform (Sigma-Aldrich Chemie GmbH, Taufkirchen, Germany) were added followed by a centrifugation at 1200 RCF for 15 min. The supernatants were then collected and used for the subsequent purification of the miRNA-enriched fractions separated from the larger RNAs, by using the miRNeasy Mini Kit and MyQIAcube instrument, according to the manufacturer’s instructions. 2.4 Library preparation and next generation sequencing (NGS) In order to assure performance in downstream applications, following purification, miRNA underwent to an internal Quality Check by using a Qubit Fluorometer (ThermoFisher Scientific, Waltham, MA, USA), with samples having a concentration lower than 15 ng/μL being defined as not compliant for the sequencing analysis. An equal quantity of miRNA (100 ng) was used for each sample analyzed. Compliant miRNAs were subjected to library preparation before sequencing (QIAseq miRNA Library Kit, Qiagen). Sequencing metrics were: 72 bp read lenght and sequencing depth of 10 milion reads, single end. The analysis was performed by using a NextSeq 550 Sequencer with NextSeq 500/550 v2.5 Kits (Illumina, San Diego, CA, USA). 2.5 Raw data preprocessing and quality controls Demultiplexing was performed using the bcl2fastq (v2.20.0.422) software (Illumina), and FASTQ files’ quality was visually inspected using FastQC (v0.11.8) (). Reads were filtered using UMI-tools (http://umi-tools.rtfd.io/) and SeqIO python package (https://pypi.org/project/seqio/), to exclude sequences shorter than 16 bp or with missing or incomplete UMI (UMI sequence shorter than 10 bp). Alignment was performed via the bowtie2 software (v2.3.4.1) using miRBase (v22.1) as a database, permitting 1 mismatch and keeping only the best alignment. Eventually, the count matrix was obtained by collapsing reads with UMI tools (with--method cluster). Only samples with at least 10000 reads mapping to the miRNAome were considered for further analysis. 2.6 Statistical analysis A two-steps statistical analysis was carried out, with a first explorative differential expression analysis and a more focused cluster analysis. 2.6.1 Differential expression analysis of miRNAs The expression level of miRNAs was first estimated using UMI tools, and then differential expression (DE) was performed using R (3.6.3). Samples were labeled based on their clinical classification (Supplementary Table S1), and case group was tested against the control group. The DE model only considered miRNA expression level as the dependent variable, with no additional covariate. Results of the DE were obtained by using the DESeq2 package (v1.26.0) (). Data normalization was performed by using DESeq2. Receiver operating characteristic (ROC) curve analysis was performed to assess the ability of the Composite Risk Index (CRI) to discriminate between cases and controls. Diagnostic performance was evaluated by calculating the area under the ROC curve (AUC) with corresponding sensitivity and specificity values. The optimal cutoff value was determined using the maximum Youden index (Supplementary Figure S1). 2.6.2 miRNA prioritization strategy Candidate miRNAs were selected through a multi-step prioritization process integrating both experimental and literature-derived evidence. First, differential expression analysis was performed on the entire salivary miRNA dataset using DESeq2. Subsequently, miRNAs were filtered according to statistical significance (p-value 10). The resulting candidate list was then compared with previously published studies investigating miRNA dysregulation in endometriosis. Particular attention was given to miRNAs involved in biological processes known to contribute to endometriosis pathogenesis, including inflammation, angiogenesis, epithelial–mesenchymal transition (EMT), fibrosis, cellular proliferation, and hormone-responsive signalling pathways. The final panel was selected by integrating statistical significance, expression abundance, consistency with previous literature, and biological relevance based on experimentally validated target genes and pathway analyses. This approach aimed to identify a biologically coherent miRNA signature rather than individual biomarkers showing isolated differential expression. 2.6.3 miRNA cluster analysis As results from the DE a cluster of 10 miRNAs was selected and evaluated against published data. The representativity of this cluster has been evaluated on the samples that previously underwent DE analysis by considering the share of sample reads which mapped to any of the miRNAs in the cluster, producing a cluster representativity index CRI for the sample. We used finer classification labels in this cluster analysis. Samples were divided based on the classification criteria described in Section 2.2 and the role of symptoms and treatment in miRNAome expression was investigated. For each of the groups, the median value of the CRIs was evaluated (Supplementary Material 1), showing a first clue of clustering behavior. 2.7 Gene ontology analyses Gene ontology (GO) analyses were performed on the genes targeted by miRNAs. Genes targeted by the single miRNAs in the cluster were extracted using miRTarBase (). GO were performed on the genes targeted by at least 5 and at least 4 miRNA using g: Profiler (Raudvere et al., 2019) version e113_eg59_p19_6be52918, accessed on January 24th, 2026. GO analysis regarded biological processes, molecular functions, and cellular components (GO:BP, GO:MF, GO:CC). Statistical domain scope was set to “Annotaded genes only” and significance threshold was set to Set Counts and Sizes (g:SCS threshold) with a threshold at 0.05. 2.8 Informed consent Written informed consent was obtained from all participants, including consent for the publication of anonymized data. 3 Results The first aim of the present study was to identify a non-invasive biomarker panel. A specific cluster of 10 salivary miRNAs, which are highly relevant to the pathophysiology of endometriosis and/or its associated processes, such as inflammation or precancerous lesions, was identified by literature review and differential expression analysis. In terms of the statistics, 7 of the 10 miRNAs (hsa-miR-130a-3p, hsa-miR-130b-3p, hsa-miR-141-3p, hsa-miR-203b-5p, hsa-miR-203a-3p, hsa-miR-29c-3p, hsa-miR-6809-5p) are part of a larger set of 86 miRNAs (47 upregulated, 39 downregulated) showing both significant differences in expression (p-value 10). Complete results of the Differential Expression procedure are available at Supplementary Material 1, where the relevant miRNAs are highlighted. In literature, 7 of the 10 miRNAs, namely hsa-miR-130a-3p, hsa-miR-130b-3p, hsa-miR-141-3p, hsa-miR-200b-3p, hsa-miR-200c-3p, hsa-miR-203b-5p, and hsa-miR-29c-3p, are known for their direct correlation with endometriosis, while remaining 3 miRNAs exhibit a clear association with inflammation and/or cancer progression (Table 1). TABLE 1 | miRNA | Main correlation | Molecular targets | References | |---|---|---|---| | hsa-miR-6809-5p | Inflammation and Cancer progression | FLOT1 | Yang et al. (2025) | | hsa-miR-200b-3p | Endometriosis | ZEB1, ZEB2, DMD, NFASC, SNAP25, HDAC4 | Ohlsson Teague et al. (2009), Zhang et al. (2021), | | hsa-miR-130a-3p | Endometriosis | SNAP25 | , | | hsa-miR-130b-3p | Endometriosis | SNAP25, SPIN90 | , , | | hsa-miR-29c-3p | Endometriosis | CDK6, BCCIP, TCF7L1, TCF7L2, PTEN, COL4A1, E‐Cadherin, and N‐Cadherin | Wentges et al. (2025), Yadav et al. (2026) | | hsa-miR-203b-5p | Endometriosis | VEGFA | Wang et al. (2022) | | hsa-miR-141-3p | Endometriosis | TGFβ2, ZEB1, ZEB2, PPP1R12A and PPP1R12B | Ohlsson Teague et al. (2009), Yang et al. (2021), | | hsa-miR-200c-3p | Endometriosis | NFASC, SNAP25, MALAT1 | , , Senturk et al. (2026), Xia et al. (2026) | | hsa-miR-27a-3p | Inflammation | IGF1 | | | hsa-miR-203a-3p | Cancer progression | SLUG, SOCS3/JNK/c-Jun, JAK1, STAT1 | Pashaei et al. (2022), Zhang et al. (2020), | Cluster of salivary miRNAs correlated to the pathophysiology of endometriosis and/or its associated processes. Given the hypothesis that endometriosis shares inflammatory and proliferative pathways with cancer, the presence of miRNAs such as hsa-miR-203b-5p (associated with inflammation) and hsa-miR-200c-3p (associated with cancer progression) in the cluster reinforces the biological validity of the panel as a potential disease signature. For classification purposes, the expression value of the cluster was normalized against the total miRNA reads per sample. As reported in Figure 2, a significant difference was observed when comparing the mean cluster expression values between patients with a first endometriosis diagnosis or under clinical management and controls represented by healthy subjects or patients with no diagnosis of endometriosis, but with other gynecological issues. FIGURE 2 The analysis of the 10-miRNA cluster had a 1.6-fold increase compared to the control group. We evaluated the discriminative ability of this index by assuming its use as a classification score. Under this assumption, the resulting classifier would achieve an AUC of 0.696 (Supplementary Figure S1). Selecting the optimal threshold (0.0273), the corresponding sensitivity and specificity would be 0.8533 and 0.5273, respectively. Despite saliva is recognized as a reliable sample matrix in bioanalytics (), its complexity is intrinsically associated with high inter-individual variability and significant biological noise (Resz et al., 2025), reason why a more rigorous and detailed clinical stratification analysis was performed. This methodological approach had the crucial dual objective of isolating the specific impact of the underlying pathology from the effect due to ongoing drug therapies, while simultaneously validating the robustness of the observed differential miRNA cluster expression pattern. To achieve this finer analysis, the two cohorts, controls and patients, were subsequently disaggregated into clinically more homogeneous sub-clusters represented by healthy controls (patients without symptoms or history of endometriosis), technical controls (patients with no diagnosis of endometriosis, but with other gynecological issues), treated patients (undergoing therapy or recent surgery), and naive patients (first endometriosis diagnosis). The results from this stratification (Figure 3) strongly confirmed the initial finding, highlighting a coherent and gradual progression of the mean miRNA cluster representation going from healthy subjects to naive patients, enlarging the perspective for this cluster to become not only a promising diagnostic tool but also a longitudinal instrument for assessing the efficacy of endometriosis treatment. FIGURE 3 As clearly depicted in Figure 3, no significant differences were observed when comparing healthy controls (0.027) to technical controls (0.028), having almost completely overlapping values. Of note, the treated group showed a mean value (0.033) that was similar to that of both control groups, but significantly lower (p < 0.001) compared to that of naive patients (0.054), showing the efficacy of the treatment (therapy or surgery) in negatively modulating the considered 10-miRNA cluster. While the increase in the cluster representation index for the cluster of genes in cases compared to the controls is clear, the different contribution of the differential expression of the single miRNAs was worth exploring. The log fold change (LFC) of differential gene expression for the single miRNAs in the cluster is shown in Figure 4. FIGURE 4 Among the miRNAs in the cluster, hsa-miR-6809-5p has the biggest differential gene expression of the 10 miRNAs in the cluster between case an CT with a LFC of 1.24 (significant with a p-val<0.01). Other miRNAs have shown a significant differential gene expression: hsa-miR-130b-3p with a LFC of 0.53 (p-val<0.001); hsa-miR-203b-5p with a LFC of 0.54 (p-val<0.05); hsa-miR-141-3p with a LFC of 0.68 (p-val<0.01). The other miRNAs in the cluster have not shown a significant differential gene expression but all of them show an increasing trend in the endometriosis group compared to controls. To improve the understanding of the role of the 10 miRNAs in the cluster and their involvement in the pathology of endometriosis, genes that are their target were extracted using miRTarBase, an annotated database based on experimentally validated miRNA-target interactions. Target genes were then ranked according to the number of miRNAs they interacted with. The genes that were common target of at least 5 miRNAs are shown in Figure 5, together with the miRNAs they are targeted by. FIGURE 5 Interestingly, the Phosphatase and Tensin homolog (PTEN) is the most central gene, targeted by 7/10 miRNAs. PTEN is an oncosuppresor gene, which targets the mTOR/PI3K/AKT pathway and that is downregulated in tumours (). Moreover, what is worth mentioning is that the miRNA hsa-miR203b-5p, which has shown a significant differential gene expression (Figure 4) does not share any of the most common target with other miRNAs, this allowing to hypothesize a different pattern of regulations it is responsible for. On the other hand, hsa-miR-200c-3p, which shows the slightest differential gene expression (Figure 4) with a LFC of 0.01, is very central and shares almost all the most common targets with other miRNAs, suggesting a role in regulating converging pathways altered in the pathology. Given the inhibitory role of miRNAs in eucaryotic organisms, all the genes targeted by the increased expressed miRNA have to be considered downregulated in the pathology (). The biological processes and molecular functions altered in relationship to the 9 genes that are common target of at least 5 miRNAs was then investigated, and results are shown in Figure 6. FIGURE 6 Two main processes emerged from the GO analysis performed, both known to be involved in endometriosis: regulation of epithelial to mesenchymal transition (EMT) () and the negative regulation of cell-substrate junctions and focal adhesion (). This shows that the miRNAs emerged from the cluster analysis performed on the miRNAome are targeting pathways strongly related to the pathology. To broaden the understanding of the regulatory network altered by the genes in the cluster, the GO analysis was extended to genes that were target of at least 4 miRNAs, this extending the analysis from 9 to 34 genes, for which the GO is shown in Figure 7. Among these, JUN, SMAD4, VEGFA are some of the most known. FIGURE 7 In this case, together with already present EMT regulation, biological processes related to embryonal development (e.g., embryo development ending in birth or egg hatching, tube development, outflow tract morphogenesis) and molecular functions related to the regulation of transcription (e.g., RNA polymerase II cis-regulatory region sequence-specific DNA binding transcription coregulator binding) emerge as significantly relevant, opening new path to study the pathophysiology of endometriosis. The results together give strong molecular support to the relevance of the cluster identified which deserves further studies as a non-invasive tool for the diagnosis of endometriosis pathology. 4 Discussion 4.1 Salivary miRNAs as non-invasive biomarkers of endometriosis Endometriosis represents a common chronic gynecological disorder affecting approximately 5%–10% of women of reproductive age; notably, its prevalence may rise to nearly 50% in selected populations, largely due to substantial underdiagnosis. This diagnostic gap is partly due to the absence of reliable, non-invasive diagnostic tools. Currently, the diagnostic “gold standard” for endometriosis remains an invasive laparoscopic surgical procedure often associated to risks of infections and bleeding typically performed in an advanced stage of the pathology (; Taylor et al., 2018). Consequently, there is a critical need for accurate, accessible, and non-invasive diagnostic approaches that would enable earlier detection, facilitate timely initiation of medical therapy, and allow longitudinal monitoring of disease recurrence without reliance on repeated laparoscopic interventions. In this context, several efforts have been directed toward the development of non-invasive diagnostic strategies for endometriosis (; Sadeghzadeh Oskouei et al., 2025), and among these the analysis of miRNA expression represents a promising approach for detecting molecular changes associated with pathogenesis and progression of several diseases, including endometriosis (). MiRNAs are present both intracellularly and in circulation and have been detected in multiple body fluids with remarkable stability (Weber et al., 2010; ). Several studies have proved the role of miRNA as potential biomarkers in the plasma of endometriosis patients, allowing to identify multi-marker panels strictly associated to the progression of this disease (Vanhie et al., 2024; Zafari et al., 2022; Rekker et al., 2015; Papari et al., 2020). Numerous miRNAs, including members of the miR-200 family and miR-17-5p, have shown strong correlations with laparoscopic biopsy-derived findings, while additional candidates such as miR-451a, let-7b, miR-20a-5p, and miR-3613-5p have been reported in serum (Zafari et al., 2022; Rekker et al., 2015; Papari et al., 2020; Nothnick et al., 2017; Yu et al., 2025). However, despite the consistency of these results, the diagnostic value of miRNAs remains limited as a consequence of the variability related to the menstrual phase during which samples are collected that may influence miRNA expression, or possible alterations of miRNA levels during proliferative and secretory phases (; Petracco et al., 2019; ). In this context, saliva-based miRNA profiling represents an attractive and, more importantly, non-invasive alternative, offering easy sample collection and potentially reduced sensitivity to systemic hormonal fluctuations. Despite the high endogenous RNase activity of saliva, salivary miRNAs have been shown to exhibit remarkable stability and reproducibility as molecular biomarkers (; ; ; Setti et al., 2020; ; ). This stability is primarily attributable to their encapsulation within extracellular vesicles, including exosomes and microvesicles, as well as their association with RNA-binding proteins such as Argonaute complexes, which protect miRNAs from enzymatic degradation. As a result, salivary miRNAs remain detectable under physiological conditions and can reliably reflect molecular alterations occurring in distant tissues. Indeed, salivary miRNA signatures have been successfully investigated as non-invasive biomarkers in a variety of pathological conditions, including cancer, inflammatory disorders, neurodegenerative diseases, and gynaecological diseases (; ). Together with the ease, safety, and non-invasive nature of sample collection, these characteristics support the use of saliva as a reliable biological matrix for assessing miRNA expression patterns associated with endometriosis pathogenesis and related inflammatory and cancer-associated pathways (; ; ; ; ). 4.2 Biological relevance of the identified miRNA signature In our work we identified a cluster of ten salivary miRNAs strongly associated with endometriosis-related pathways (Table 1). Specifically, seven miRNAs (hsa-miR-130a-3p, hsa-miR-130b-3p, hsa-miR-141-3p, hsa-miR-200b-3p, hsa-miR-200c-3p, hsa-miR-203b-5p, and hsa-miR-29c-3p) have been linked directly to endometriosis, while the remaining three are primarily associated with inflammatory processes and cancer-related pathways. Several miRNAs previously reported in association with endometriosis, including miR-17-5p, miR-451a, and let-7b-5p, were evaluated during the literature review process. However, these candidates were not retained in the final panel because they did not meet all predefined prioritization criteria within the present salivary NGS dataset. To date, only a limited number of studies have investigated salivary miRNAs in endometriosis; for instance, an untargeted high-throughput analysis identified a 109-miRNA salivary signature that lacked independent validation (). In contrast, targeted approaches have identified hsa-miR-135a as significantly upregulated in both saliva and plasma of women with endometriosis (Perricos et al., 2022), supporting its systemic relevance. Consistent with these findings, the inclusion of inflammation- and cancer-associated miRNAs supports the biological relevance of the identified cluster as an endometriosis-associated molecular signature. After normalization, the 10-miRNA cluster showed a significant 1.6-fold increase in patients with endometriosis compared with controls (Figure 2), aligning with mechanistic evidence that key miRNAs implicated in endometriosis, such as hsa-miR-135a, directly target HOXA10 and contribute to impaired endometrial receptivity, dysregulated hormonal signaling, and ectopic lesion development (Petracco et al., 2011). To further validate our findings, the cohort was stratified into healthy controls, technical controls, treated patients, and treatment-naive patients, revealing a progressive increase in miRNA cluster expression from healthy subjects to naive endometriosis patients. Notably, treated patients showed expression levels comparable to controls and significantly lower than naive patients, indicating that medical or surgical treatment effectively modulates this miRNA signature (Figure 3). We investigated the specific contribution of each miRNA belonging to the identified cluster associated with endometriosis (Figure 4) and despite the evident overall increase of the miRNA cluster in endometriosis patients, analysis of individual components revealed heterogeneous contributions. hsa-miR-6809-5p showed the strongest differential expression between cases and controls, followed by significant upregulation of hsa-miR-130b-3p, hsa-miR-203b-5p, and hsa-miR-141-3p. While the remaining miRNAs did not reach statistical significance individually, all exhibited a consistent upward trend in endometriosis patients, supporting the cumulative discriminatory power of the cluster. Our results are in line with the previous data described by Filigheddu et al., in which our salivary 10-miRNA cluster overlap with differential miRNA depicted in both eutopic and ectopic endometrial tissues derived from women affected by endometriosis (). In particular, it has been described a dysregulation of members of the miR-200 family, including hsa-miR-200c-3p and hsa-miR-141-3p, in ectopic lesions strictly correlated with pathways involved in epithelial–mesenchymal transition (EMT), cell proliferation, and hormonal signaling. The overlap of tissue and salivary miRNAs suggests that salivary profiles mirror tissue-level molecular alterations in endometriosis, supporting their value as non-invasive biomarkers. 4.3 miRNA-target networks and pathogenic pathways To better investigate the biological role of the 10-miRNA cluster and, mainly, their involvement in the pathological processes undergone endometriosis, experimentally validated target genes have been identified using miRTarBase and ranked based on the number of interacting miRNAs. PTEN emerged as the most central target, being regulated by 7 of the 10 miRNAs (Figure 5), highlighting convergence on the PI3K/AKT/mTOR pathway, which is known to be dysregulated in endometriosis (). In fact, PTEN-mediated regulation of PI3K/AKT signaling is well established in endometrial cancer, where PTEN loss drives disease progression by promoting cellular proliferation and resistance to TGFβ-induced apoptosis. In this context, the miR-424 (322)∼503 cluster has been identified as a key regulator of these processes (Vidal-Sabanés et al., 2025). Included in our 10 miRNAs cluster, hsa-miR-203b-5p has shown a distinct regulatory profile, targeting unique genes not shared with other cluster members, while hsa-miR-200c-3p, even though minimal differential expression, has occupied a central position by sharing most common targets, suggesting a role in coordinating convergent pathogenic pathways (Figure 5). The specific crosstalk between miRNA dysregulation in endometriosis and other chronic proliferative disorders is represented by the hsa-miR-203b-5p and hsa-miR-200c-3p interconnection. Hsa-miR-200c-3p is a key regulator of EMT and cellular plasticity, processes that are crucial for both tumor invasion and endometriotic lesion establishment. In contrast, hsa-miR-203b-5p is primarily associated with inflammatory signaling and epithelial differentiation, suggesting a complementary role in modulating immune responses and tissue remodeling. The concurrent dysregulation of these miRNAs supports a model in which endometriosis shares cancer-like regulatory networks involving EMT, inflammation, and aberrant cell survival, while maintaining disease-specific molecular features. Our results are corroborated by the study of Góźdź et al. (), which demonstrated that EMT-related genes and miRNAs, including miR-200 family, are constitutively expressed in eutopic endometrium and are mostly regulated by the menstrual cycle rather than by disease status, with no evidence of overt EMT activation in eutopic tissue from endometriosis patients. Based on these concepts, the presence of miR-200c-3p within our 10-miRNA cluster suggests that salivary miRNAs more likely reflect systemic or ectopic lesion-derived molecular alterations rather than eutopic endometrial EMT changes, supporting the idea that salivary miRNA profiles reflect signals from ectopic lesions and inflammatory microenvironments, reinforcing, once again, their utility as non-invasive biomarkers tool. Our findings related to the functional enrichment analysis of genes further support the disruption of key regulatory networks in endometriosis, identifying at least five miRNAs (Figure 7) involved in two major biological processes as EMT and negative regulation of cell-substrate junctions and focal adhesion, both of which exert a central role in the endometriosis pathophysiology. Extending the analysis to genes targeted by multiple miRNAs, we identified SMAD4, together with JUN and VEGFA, as central nodes within miRNA–mRNA interaction networks. In fact, in line with our findings, it has been suggested that miR-542-3p-mediated repression of the BMP7–SMAD4–CDH1 axis may contribute to the acquisition of mesenchymal-like features in ectopic endometrial lesions. Notably, the concomitant downregulation of BMP7, SMAD4, and CDH1 in ectopic tissue and the augmented miR-542-3p expression level suggest that post-transcriptional regulation may play a role in destabilizing epithelial integrity and promoting EMT-related processes (Zubrzycka et al., 2023). Concerning the JUN pathway, which is primarily implicated in endometriosis through the AP-1 transcription factor complex, evidence indicates that it plays a key role in promoting invasive behavior. JUN/AP-1 activity has been associated with the regulation of EMT-related genes, adhesion molecules, and extracellular matrix remodelers, and experimental knockdown of JUN has been shown to suppress invasion in highly invasive endometriotic cells (Wilson et al., 2022). Emerging evidence suggests that miRNA dysregulation may contribute to altered JUN/AP-1 signaling in endometriosis. Several disease-associated miRNAs are predicted or validated regulators of JUN or its upstream activators, and integrative miRNA–mRNA analyses frequently identify JUN as a hub gene within networks linked to inflammation and angiogenesis, supporting a role for post-transcriptional regulation of AP-1-related pathways in disease pathogenesis (Ohlsson Teague et al., 2009; ; ; ). Finally, a recent multi-fluid miRNA profiling study further supports the involvement of systemic miRNA dysregulation in endometriosis (). Distinct miRNA signatures were identified across serum, saliva, and vaginal mucus, yet pathway analyses converged on inflammation, remodeling, angiogenesis, and TGF-β-related signaling, emphasizing network-level regulation rather than single biomarkers (; ; ; Suárez and Sessa, 2009; Sun et al., 2018; ). This is consistent with our identification of a coordinately upregulated salivary miRNA cluster and the convergence of multiple miRNAs on key targets such as PTEN, JUN, and SMAD4, supporting a role for post-transcriptional control of central signaling nodes in endometriosis pathogenesis. 4.4 Clinical implications and future perspectives Collectively, our study identifies a salivary miRNA signature associated with endometriosis, supporting the use of saliva as a feasible and truly non-invasive biological matrix for disease investigation. The coordinated upregulation of a miRNA cluster converging on key signaling nodes, including PTEN, JUN, and SMAD4, highlights the relevance of miRNA-mediated regulatory networks in pathways related to inflammation, angiogenesis, tissue remodeling, and EMT. Notably, our findings suggest that systemic miRNA dysregulation characterizes endometriosis and that salivary miRNAs may reflect underlying pathogenic mechanisms regulated by treatment. However, a possible limitation of the present pilot study is the lack of independent qRT-PCR validation of the 10-miRNA signature identified by NGS. Future studies will validate the expression of the selected salivary miRNAs by qRT-PCR, including direct comparison between healthy controls and treatment-naïve endometriosis patients, strengthening the robustness and reproducibility of our results. Although further validation is required, these results indicate that miRNA-based signatures detected in saliva could represent promising predictive biomarkers for endometriosis, offering potential clinical utility for early, non-invasive diagnosis and disease monitoring. Statements Data availability statement The data presented in the study are deposited in the Figshare repository, available at https://doi.org/10.6084/m9.figshare.32968895. Ethics statement The studies involving humans were approved by all procedures in this study were conducted in accordance with the Helsinki Declaration. The study was approved by the Regional Ethics Committee for Clinical Trials of the Region of Tuscany Area Vasta Sud EST (ID23576), by Lombardy Regional Ethics Committee 3 (ID 4165_S_N). The study has also been submitted and registered on ClinicalTrials.gov (ID protocol number: NCT06100471 registration date 10-13-2023). All women enrolled in the study were properly informed and signed an informed consent. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article. Author contributions GrC: Project administration, Conceptualization, Writing – review and editing, Resources, Writing – original draft, Funding acquisition, Visualization, Formal Analysis. FB: Methodology, Writing – review and editing, Investigation. CP: Investigation, Writing – review and editing. LF: Writing – review and editing. NM: Methodology, Investigation, Writing – review and editing. RG: Writing – review and editing, Formal Analysis. FD: Writing – review and editing, Formal Analysis. AG: Formal Analysis, Investigation, Writing – review and editing. LD: Visualization, Writing – original draft, Writing – review and editing, Formal Analysis. LB: Formal Analysis, Writing – review and editing, Writing – original draft. FM: Writing – review and editing. DD: Writing – review and editing. GL: Writing – review and editing. SC: Writing – review and editing. GM: Writing – review and editing. SLu: Writing – review and editing. GD: Writing – review and editing. EG: Writing – review and editing. SLa: Writing – review and editing. MM: Writing – review and editing. AL: Writing – review and editing. VP: Writing – review and editing. SF: Writing – review and editing. GiC: Writing – original draft, Formal Analysis, Visualization, Supervision, Writing – review and editing, Conceptualization. Funding The author(s) declared that financial support was not received for this work and/or its publication. Acknowledgments The authors thank the participants who donated saliva samples for this study, the nurses, and lab technician involved in and sample collection, and preparation of the biological samples. The authors also thank the Next Generation Sequencing (NGS) Platform (Eurofins, Genoma) for miRNA library preparation and sequencing, every single professional healthcare that helps with patient recruitment and selection. A special thanks is deserved to the memory of Prof. Massimo Luerti, who was one of the pioneers of the study related to miRNA and endometriosis. Conflict of interest Authors GrC, FB, CP, LF, NM, RG, and FD were employed by company Eurofins Genoma Group. Author GM was employed by the company GVM S. Catherine of Siena Clinic. Author AL was employed by the company Hera Center. The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author GC declared that he was an editorial board member of Frontiers at the time of submission. This had no impact on the peer review process and the final decision. Generative AI statement The author(s) declared that generative AI was not used in the creation of this manuscript. Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us. Publisher’s note All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher. Supplementary material The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2026.1869937/full#supplementary-material SUPPLEMENTARY FIGURES 1Receiver operating characteristic (ROC) curve for CRI-based discrimination between cases and controls. The ROC analysis yielded an area under the curve (AUC) of 0.696, indicating moderate discriminatory performance. The optimal cutoff value (CRI = 0.03) was associated with a sensitivity of 85% and a specificity of 53%. The diagonal line represents the performance of a random classifier. SUPPLEMENTARY TABLES 1Differential expression analysis of circulating miRNAs in cases versus controls. Differential expression analysis was performed to identify miRNAs significantly altered between cases and controls. For each miRNA, the table reports mean normalized counts, log fold change (LFC), raw p-value, and multiple-testing adjusted p-value (false discovery rate, FDR). Positive LFC values indicate higher expression in cases relative to controls, whereas negative LFC values indicate lower expression in cases. The “DE” column indicates whether the miRNA was classified as differentially expressed according to the predefined statistical significance criteria.

1 AhnS.KwonA.HuhY. H.RheeS.SongW. K. (2022). Tumor-derived miR-130b-3p induces cancer-associated fibroblast activation by targeting SPIN90 in luminal A breast cancer. Oncogenesis11, 47. 10.1038/s41389-022-00422-6 2 Anonymous. (1985). Revised american fertility society classification of endometriosis: 1985. Fertil. Steril.43, 351–352. 10.1016/s0015-0282(16)48430-x 3 Anonymous. (1997). Revised American society for reproductive medicine classification of endometriosis: 1996. Fertil. Steril.67, 817–821. 10.1016/s0015-0282(97)81391-x 4 AlimiY.IwanagaJ.LoukasM.TubbsR. S. (2018). The clinical anatomy of endometriosis: a review. Cureus10, e3361. 10.7759/cureus.3361 5 AnafV.SimonP.El NakadiI.FaytI.BuxantF.SimonartT.et al (2000). Relationship between endometriotic foci and nerves in rectovaginal endometriotic nodules. Hum. Reprod.15, 1744–1750. 10.1093/humrep/15.8.1744 6 AnafV.ChapronC.El NakadiI.De MoorV.SimonartT.NoelJ. C. (2006). Pain, mast cells, and nerves in peritoneal, ovarian, and deep infiltrating endometriosis. Fertil. Steril.86, 1336–1343. 10.1016/j.fertnstert.2006.03.057 7 AndrewsS. (2010). Fastqc: a quality control tool for high throughput sequence data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc.Fastqc/(accessed on July 11, 2018). 8 As-SanieS.HarrisR. E.NapadowV.KimJ.NeshewatG.KairysA.et al (2012). Changes in regional gray matter volume in women with chronic pelvic pain: a voxel-based morphometry study. Pain153, 1006–1014. 10.1016/j.pain.2012.01.032 9 AveryJ. C.DeslandesA.FregerS. M.LeonardiM.LoG.CarneiroG.et al (2024). Noninvasive diagnostic imaging for endometriosis part 1: a systematic review of recent developments in ultrasound, combination imaging, and artificial intelligence. Fertil. Steril.121, 164–188. 10.1016/j.fertnstert.2023.12.008 10 BaxterE.NairS.WestZ.SalomonC.ObermairA. (2025). Extracellular vesicles as biomarkers for endometrial cancer - a systematic review. Transl. Oncol.62, 102543. 10.1016/j.tranon.2025.102543 11 BendifallahS.SuisseS.PucharA.DelbosL.PoilblancM.DescampsP.et al (2022). Salivary MicroRNA signature for diagnosis of endometriosis. J. Clin. Med.11, 612. 10.3390/jcm11030612 12 BendifallahS.DabiY.SuisseS.DelbosL.SpiersA.PoilblancM.et al (2023). Validation of a salivary miRNA signature of endometriosis - interim data. NEJM Evid.2, EVIDoa2200282. 10.1056/evidoa2200282 13 BendifallahS.DabiY.SuisseS.IlicJ.DelbosL.PoilblancM.et al (2024). Saliva-based microRNA diagnostic signature for the superficial peritoneal endometriosis phenotype. Eur. J. Obstet. Gynecol. Reprod. Biol.297, 187–196. 10.1016/j.ejogrb.2024.04.020 14 BendifallahS.RomanH.SuisseS.SpiersA.PetitE.DelbosL.et al (2025). Validation of a saliva Micro-RNA signature for endometriosis. NEJM Evid.4, EVIDoa2400195. 10.1056/evidoa2400195 15 BjorkmanS.TaylorH. S. (2019). Erratum: MicroRNAs in endometriosis biological function and emerging biomarker candidates. Biol. Reprod.101, 1179. 10.1093/biolre/ioz219 16 BonavinaG.TaylorH. S. (2022). Endometriosis-associated infertility: from pathophysiology to tailored treatment. Front. Endocrinol. (Lausanne)13, 1020827. 10.3389/fendo.2022.1020827 17 BricouA.BattR. E.ChapronC. (2008). Peritoneal fluid flow influences anatomical distribution of endometriotic lesions: why sampson seems to be right. Eur. J. Obstet. Gynecol. Reprod. Biol.138, 127–134. 10.1016/j.ejogrb.2008.01.014 18 BurneyR. O.GiudiceL. C. (2012). Pathogenesis and pathophysiology of endometriosis. Fertil. Steril.98, 511–519. 10.1016/j.fertnstert.2012.06.029 19 CarlinoC.RippoM. R.LazzariniR.MonsurròV.MorroneS.AngeliniS.et al (2018). Differential microRNA expression between decidual and peripheral blood natural killer cells in early pregnancy. Hum. Reprod.33, 2184–2195. 10.1093/humrep/dey323 20 ChapronC.FauconnierA.VieiraM.BarakatH.DoussetB.PansiniV.et al (2003). Anatomical distribution of deeply infiltrating endometriosis: surgical implications and proposition for a classification. Hum. Reprod.18, 157–161. 10.1093/humrep/deg009 21 ChapronC.BourretA.ChopinN.DoussetB.LeconteM.Amsellem-OuazanaD.et al (2010). Surgery for bladder endometriosis: long-term results and concomitant management of associated posterior deep lesions. Hum. Reprod.25, 884–889. 10.1093/humrep/deq017 22 ChapronC.SantulliP.de ZieglerD.NoelJ. C.AnafV.StreuliI.et al (2012). Ovarian endometrioma: severe pelvic pain is associated with deeply infiltrating endometriosis. Hum. Reprod.27, 702–711. 10.1093/humrep/der462 23 ChoS.MutluL.GrechukhinaO.TaylorH. S. (2015). Circulating microRNAs as potential biomarkers for endometriosis. Fertil. Steril.103, 1252–1260.e1251. 10.1016/j.fertnstert.2015.02.013 24 ChoiH. J.ChungT. W.ChoiH. J.HanJ. H.ChoiJ. H.KimC. H.et al (2018). Increased α2-6 sialylation of endometrial cells contributes to the development of endometriosis. Exp. Mol. Med.50, 1–12. 10.1038/s12276-018-0167-1 25 CuiS.YuS.HuangH. Y.LinY. C.HuangY.ZhangB.et al (2025). miRTarBase 2025: updates to the collection of experimentally validated microRNA-target interactions. Nucleic Acids Res.53, D147–d156. 10.1093/nar/gkae1072 26 DabiY.SuisseS.PucharA.DelbosL.PoilblancM.DescampsP.et al (2023). Endometriosis-associated infertility diagnosis based on saliva microRNA signatures. Reprod. Biomed. Online46, 138–149. 10.1016/j.rbmo.2022.09.019 27 Di PietroC.CarusoS.BattagliaR.Iraci SareriM.La FerlitaA.StrinoF.et al (2018). MiR-27a-3p and miR-124-3p, upregulated in endometrium and serum from women affected by chronic endometritis, are new potential molecular markers of endometrial receptivity. Am. J. Reprod. Immunol.80, e12858. 10.1111/aji.12858 28 DunselmanG. A.VermeulenN.BeckerC.Calhaz-JorgeC.D'HoogheT.De BieB.et al (2014). ESHRE guideline: management of women with endometriosis. Hum. Reprod.29, 400–412. 10.1093/humrep/det457 29 FaurC. I.RotaruH.OsanC.JurjA.RomanR. C.MoldovanM.et al (2021). Salivary exosomal microRNAs as biomarkers for head and neck cancer detection-a literature review. Maxillofac. Plast. Reconstr. Surg.43, 19. 10.1186/s40902-021-00303-9 30 FehlmannT.LudwigN.BackesC.MeeseE.KellerA. (2016). Distribution of microRNA biomarker candidates in solid tissues and body fluids. RNA Biol.13, 1084–1088. 10.1080/15476286.2016.1234658 31 FerreroS.ScalaC.RaccaA.CalanniL.RemorgidaV.VenturiniP. L.et al (2015). Second surgery for recurrent unilateral endometriomas and impact on ovarian reserve: a case-control study. Fertil. Steril.103, 1236–1243. 10.1016/j.fertnstert.2015.01.032 32 FilighedduN.GregnaninI.PorporatoP. E.SuricoD.PeregoB.GalliL.et al (2010). Differential expression of microRNAs between eutopic and ectopic endometrium in ovarian endometriosis. J. Biomed. Biotechnol.2010, 369549. 10.1155/2010/369549 33 GaiC.CamussiF.BroccolettiR.GambinoA.CabrasM.MolinaroL.et al (2018). Salivary extracellular vesicle-associated miRNAs as potential biomarkers in oral squamous cell carcinoma. BMC Cancer18, 439. 10.1186/s12885-018-4364-z 34 GalloA.TandonM.AlevizosI.IlleiG. G. (2012). The majority of microRNAs detectable in serum and saliva is concentrated in exosomes. PLoS One7, e30679. 10.1371/journal.pone.0030679 35 GebertL. F. R.MacRaeI. J. (2019). Regulation of microRNA function in animals. Nat. Rev. Mol. Cell Biol.20, 21–37. 10.1038/s41580-018-0045-7 36 GóźdźA.ŻeberkiewiczM.JaniukI.HycA.IwanA.ZwierzchowskaA.et al (2025). The influence of menstrual cycle and endometriosis on endometrial expression of epithelial-to-mesenchymal transition (EMT)-Related genes. Reprod. Biol. Endocrinol.23, 144. 10.1186/s12958-025-01486-w 37 GröschlM. (2017). Saliva: a reliable sample matrix in bioanalytics. Bioanalysis9, 655–668. 10.4155/bio-2017-0010 38 GuoL.ZhangY.ZhangL.HuangF.LiJ.WangS. (2016). MicroRNAs, TGF-β signaling, and the inflammatory microenvironment in cancer. Tumour Biol.37, 115–125. 10.1007/s13277-015-4374-2 39 HashmiA. S.GuptaN. D.KhanS.AliS. A.KulsumO.VelloneV. (2026). Diagnostic and prognostic potential of salivary microRNA in oral and head and neck squamous cell carcinomas: a systematic review. Int. J. Oral Maxillofac. Surg.55, 621–629. 10.1016/j.ijom.2026.01.004 40 HillM.TranN. (2021). miRNA interplay: mechanisms and consequences in cancer. Dis. Model Mech.14, dmm047662. 10.1242/dmm.047662 41 HoffmanD. (2015). Central and peripheral pain generators in women with chronic pelvic pain: patient centered assessment and treatment. Curr. Rheumatol. Rev.11, 146–166. 10.2174/1573397111666150619094524 42 HosseiniradH.JeongJ. W.BarrierB. F. (2025). Insights into the molecular mechanisms and signaling pathways of epithelial to mesenchymal transition (EMT) in the pathophysiology of endometriosis. Int. J. Mol. Sci.26, 7460. 10.3390/ijms26157460 43 HsuC. Y.HsiehT. H.TsaiC. F.TsaiH. P.ChenH. S.ChangY.et al (2014). miRNA-199a-5p regulates VEGFA in endometrial mesenchymal stem cells and contributes to the pathogenesis of endometriosis. J. Pathol.232, 330–343. 10.1002/path.4295 44 HufnagelD.LiF.CosarE.KrikunG.TaylorH. S. (2015). The role of stem cells in the etiology and pathophysiology of endometriosis. Semin. Reprod. Med.33, 333–340. 10.1055/s-0035-1564609 45 JalouliM.RahmanM. A.NahdiS.HarrathA. H. (2026). Epigenetic landscape of female infertility: an integrated bioinformatics perspective on DNA methylation, MicroRNAs, and gene regulatory networks across PCOS, endometriosis, and diminished ovarian reserve. Int. J. Mol. Sci.27, 1785. 10.3390/ijms27041785 46 JiaoY.ZhangT.ZhangC.JiH.TongX.XiaR.et al (2021). Exosomal miR-30d-5p of neutrophils induces M1 macrophage polarization and primes macrophage pyroptosis in sepsis-related acute lung injury. Crit. Care25, 356. 10.1186/s13054-021-03775-3 47 Kaczor-UrbanowiczK. E.Martin Carreras-PresasC.AroK.TuM.Garcia-GodoyF.WongD. T. (2017). Saliva diagnostics - current views and directions. Exp. Biol. Med. (Maywood)242, 459–472. 10.1177/1535370216681550 48 KhanK. N.KitajimaM.FujishitaA.HirakiK.MatsumotoA.NakashimaM.et al (2013). Pelvic pain in women with ovarian endometrioma is mostly associated with coexisting peritoneal lesions. Hum. Reprod.28, 109–118. 10.1093/humrep/des364 49 kumarP.KumarA.KumarA.KumarH. (2025). Essential role of hsa-miR-203a-3p in type I interferons immune homeostasis during influenza and NDV infection. bioRxiv. 10.1101/2025.04.01.646562 50 Laux-BiehlmannA.d'HoogheT.ZollnerT. M. (2015). Menstruation pulls the trigger for inflammation and pain in endometriosis. Trends Pharmacol. Sci.36, 270–276. 10.1016/j.tips.2015.03.004 51 LeonovaA.TurpinV. E.AgarwalS. K.LeonardiM.FosterW. G. (2021). A critical appraisal of the circulating levels of differentially expressed microRNA in endometriosis†. Biol. Reprod.105, 1075–1085. 10.1093/biolre/ioab134 52 LeyendeckerG.BilgicyildirimA.InackerM.StalfT.HuppertP.MallG.et al (2015). Adenomyosis and endometriosis. Re-visiting their association and further insights into the mechanisms of auto-traumatisation. An MRI study. Arch. Gynecol. Obstet.291, 917–932. 10.1007/s00404-014-3437-8 53 LiX.CuiD.LiuY. (2025). Combination of circular RNA-miRNA-mRNA expression profiles and bioinformatic analysis in ovarian endometriosis. Genomics117, 111112. 10.1016/j.ygeno.2025.111112 54 LiangZ.ChenY.ZhaoY.XuC.ZhangA.ZhangQ.et al (2017). miR-200c suppresses endometriosis by targeting MALAT1 in vitro and in vivo. Stem Cell Res. Ther.8, 251. 10.1186/s13287-017-0706-z 55 LiuY.FanL.JinL.LuC.LiT.ZhangZ.et al (2022). Integrated bioinformatic analysis of dysregulated microRNA-mRNA co-expression network in ovarian endometriosis. Acta Obstet. Gynecol. Scand.101, 1074–1084. 10.1111/aogs.14430 56 LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biology15, 550. 10.1186/s13059-014-0550-8 57 LyuS.LiQ.GuZ.YanH.TangX.DaiY.et al (2026). Identification of candidate microRNA biomarkers of endometriosis in different bodily fluids. Sci. Rep.16, 6218. 10.1038/s41598-026-37277-5 58 MadanesD.BilotasM. A.BastónJ. I.SinglaJ. J.MeresmanG. F.BarañaoR. I.et al (2020). PI3K/AKT pathway is altered in the endometriosis patient's endometrium and presents differences according to severity stage. Gynecol. Endocrinol.36, 436–440. 10.1080/09513590.2019.1680627 59 MaphuthaJ.TwilleyD.LallN. (2024). The role of the PTEN tumor suppressor gene and its anti-angiogenic activity in melanoma and other cancers. Molecules29, 721. 10.3390/molecules29030721 60 MarquardtR. M.TranD. N.LesseyB. A.RahmanM. S.JeongJ. W. (2023). Epigenetic dysregulation in endometriosis: implications for pathophysiology and therapeutics. Endocr. Rev.44, 1074–1095. 10.1210/endrev/bnad020 61 McKinnonB. D.BertschiD.BersingerN. A.MuellerM. D. (2015). Inflammation and nerve fiber interaction in endometriotic pain. Trends Endocrinol. Metab.26, 1–10. 10.1016/j.tem.2014.10.003 62 MenniK.FacchettiL.CabassaP. (2016). Extragenital endometriosis: assessment with MR imaging. A pictorial review. Br. J. Radiol.89, 20150672. 10.1259/bjr.20150672 63 MichaelA.BajracharyaS. D.YuenP. S.ZhouH.StarR. A.IlleiG. G.et al (2010). Exosomes from human saliva as a source of microRNA biomarkers. Oral Dis.16, 34–38. 10.1111/j.1601-0825.2009.01604.x 64 MiscianinovV.MartelloA.RoseL.ParishE.CathcartB.MitićT.et al (2018). MicroRNA-148b targets the TGF-β pathway to regulate angiogenesis and endothelial-to-mesenchymal transition during skin wound healing. Mol. Ther.26, 1996–2007. 10.1016/j.ymthe.2018.05.002 65 MoustafaS.BurnM.MamillapalliR.NematianS.FloresV.TaylorH. S. (2020). Accurate diagnosis of endometriosis using serum microRNAs. Am. J. Obstet. Gynecol.223, 557.e551. 10.1016/j.ajog.2020.02.050 66 MuziiL.BianchiA.BellatiF.CristiE.PerniceM.ZulloM. A.et al (2007). Histologic analysis of endometriomas: what the surgeon needs to know. Fertil. Steril.87, 362–366. 10.1016/j.fertnstert.2006.06.055 67 NaftalinJ.HooW.PatemanK.MavrelosD.FooX.JurkovicD. (2014). Is adenomyosis associated with menorrhagia?Hum. Reprod.29, 473–479. 10.1093/humrep/det451 68 NothnickW. B.FalconeT.JoshiN.FazleabasA. T.GrahamA. (2017). Serum miR-451a levels are significantly elevated in women with endometriosis and recapitulated in baboons (Papio anubis) with experimentally-induced disease. Reprod. Sci.24, 1195–1202. 10.1177/1933719116681519 69 Ohlsson TeagueE. M.Van der HoekK. H.Van der HoekM. B.PerryN.WagaarachchiP.RobertsonS. A.et al (2009). MicroRNA-regulated pathways associated with endometriosis. Mol. Endocrinol.23, 265–275. 10.1210/me.2008-0387 70 PapariE.NoruziniaM.KashaniL.FosterW. G. (2020). Identification of candidate microRNA markers of endometriosis with the use of next-generation sequencing and quantitative real-time polymerase chain reaction. Fertil. Steril.113, 1232–1241. 10.1016/j.fertnstert.2020.01.026 71 PashaeiM.MashayekhiF.ZahiriZ.SalehiZ. (2022). miR-203a-3p, ABL1 and TP63 gene expression is altered in the endometrium of women with endometriosis. Gynecol. Endocrinol.38, 603–607. 10.1080/09513590.2022.2076830 72 PerricosA.ProestlingK.HussleinH.KuesselL.HudsonQ. J.WenzlR.et al (2022). Hsa-mir-135a shows potential as A putative diagnostic biomarker in saliva and plasma for endometriosis. Biomolecules, 12(8), 1144. 10.3390/biom12081144 73 PetraccoR.GrechukhinaO.PopkhadzeS.MassasaE.ZhouY.TaylorH. S. (2011). MicroRNA 135 regulates HOXA10 expression in endometriosis. J. Clin. Endocrinol. Metab.96, E1925–E1933. 10.1210/jc.2011-1231 74 PetraccoR.DiasA. C. O.TaylorH.PetraccoÁ.BadalottiM.MichelonJ. D. R.et al (2019). Evaluation of miR-135a/b expression in endometriosis lesions. Biomed. Rep.11, 181–187. 10.3892/br.2019.1237 75 RaffiF.MetwallyM.AmerS. (2012). The impact of excision of ovarian endometrioma on ovarian reserve: a systematic review and meta-analysis. J. Clin. Endocrinol. Metab.97, 3146–3154. 10.1210/jc.2012-1558 76 RahmiogluN.MissmerS. A.MontgomeryG. W.ZondervanK. T. (2012). Insights into assessing the genetics of endometriosis. Curr. Obstet. Gynecol. Rep.1, 124–137. 10.1007/s13669-012-0016-5 77 RaudvereU.KolbergL.KuzminI.ArakT.AdlerP.PetersonH.et al (2019). g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res.47, W191–w198. 10.1093/nar/gkz369 78 RekkerK.SaareM.RoostA. M.KaartT.SõritsaD.KarroH.et al (2015). Circulating miR-200-family micro-RNAs have altered plasma levels in patients with endometriosis and vary with blood collection time. Fertil. Steril.104, 938–946.e932. 10.1016/j.fertnstert.2015.06.029 79 ResztakM.CzyrskiA.SobiakJ. (2025). Saliva as a matrix for therapeutic drug monitoring and disease biomarkers in children and adolescents. Pharmacol. Rep.77, 921–961. 10.1007/s43440-025-00732-7 80 Sadeghzadeh OskoueiB.AsadiZ.Jahanban EsfahlanR. (2025). Non-invasive blood tests for earlier diagnosis and treatment of endometriosis. J. Reprod. Immunol.169, 104521. 10.1016/j.jri.2025.104521 81 SantulliP.LamauM. C.MarcellinL.GayetV.MarzoukP.BorgheseB.et al (2016). Endometriosis-related infertility: ovarian endometrioma per se is not associated with presentation for infertility. Hum. Reprod.31, 1765–1775. 10.1093/humrep/dew093 82 SenturkM.Gumusoglu-AcarE.TopuzS.GunelT. (2026). Effect of the hsa-microRNA 200c-3p/ZEB1 loop in epithelial-mesenchymal transition of ovarian cancer cells. Oncol. Lett.31, 131–138. 10.3892/ol.2026.15484 83 SettiG.PezziM. E.VianiM. V.PertinhezT. A.CassiD.MagnoniC.et al (2020). Salivary MicroRNA for diagnosis of cancer and systemic diseases: a systematic review. Int. J. Mol. Sci.21, 907. 10.3390/ijms21030907 84 StefanssonH.GeirssonR. T.SteinthorsdottirV.JonssonH.ManolescuA.KongA.et al (2002). Genetic factors contribute to the risk of developing endometriosis. Hum. Reprod.17, 555–559. 10.1093/humrep/17.3.555 85 StreuliI.de ZieglerD.GayetV.SantulliP.BijaouiG.de MouzonJ.et al (2012). In women with endometriosis Anti-Mullerian hormone levels are decreased only in those with previous endometrioma surgery. Hum. Reprod.27, 3294–3303. 10.1093/humrep/des274 86 SuárezY.SessaW. C. (2009). MicroRNAs as novel regulators of angiogenesis. Circ. Res.104, 442–454. 10.1161/CIRCRESAHA.108.191270 87 SunL. L.LiW. D.LeiF. R.LiX. Q. (2018). The regulatory role of microRNAs in angiogenesis-related diseases. J. Cell Mol. Med.22, 4568–4587. 10.1111/jcmm.13700 88 SuryawanshiS.VladA. M.LinH. M.Mantia-SmaldoneG.LaskeyR.LeeM.et al (2013). Plasma microRNAs as novel biomarkers for endometriosis and endometriosis-associated ovarian cancer. Clin. Cancer Res.19, 1213–1224. 10.1158/1078-0432.CCR-12-2726 89 TaylorH. S.AdamsonG. D.DiamondM. P.GoldsteinS. R.HorneA. W.MissmerS. A.et al (2018). An evidence-based approach to assessing surgical versus clinical diagnosis of symptomatic endometriosis. Int. J. Gynaecol. Obstet.142, 131–142. 10.1002/ijgo.12521 90 TreloarS. A.O'ConnorD. T.O'ConnorV. M.MartinN. G. (1999). Genet. Influ. Endometr. Aust. twin sample. [email protected]. Fertil Steril71, 701–710.10.1016/s0015-0282(98)00540-8 91 TreloarS.HadfieldR.MontgomeryG.LambertA.WicksJ.BarlowD. H.et al (2002). The International Endogene Study: a collection of families for genetic research in endometriosis. Fertil. Steril.78, 679–685. 10.1016/s0015-0282(02)03341-1 92 VanhieA.OD.PeterseD.BeckersA.CuéllarA.FassbenderA.et al (2019). Plasma miRNAs as biomarkers for endometriosis. Hum. Reprod.34, 1650–1660. 10.1093/humrep/dez116 93 VanhieA.CaronE.VermeerschE.OD.TomassettiC.MeulemanC.et al (2024). Circulating microRNAs as non-invasive biomarkers in endometriosis Diagnosis-A systematic review. Biomedicines12, 888. 10.3390/biomedicines12040888 94 VercelliniP.ConsonniD.DridiD.BraccoB.FrattaruoloM. P.SomiglianaE. (2014a). Uterine adenomyosis and in vitro fertilization outcome: a systematic review and meta-analysis. Hum. Reprod.29, 964–977. 10.1093/humrep/deu041 95 VercelliniP.ViganoP.SomiglianaE.FedeleL. (2014b). Endometriosis: pathogenesis and treatment. Nat. Rev. Endocrinol.10, 261–275. 10.1038/nrendo.2013.255 96 Vidal-SabanésM.BonifaciN.NavaridasR.EgeaJ.EncinasM.Rodriguez-BarruecoR.et al (2025). Endometrial cancer progression driven by PTEN-deficiency requires miR-424 (322)∼ 503. Cell Death and Dis.16, 705. 10.1038/s41419-025-08022-z 97 WangW. T.ZhaoY. N.HanB. W.HongS. J.ChenY. Q. (2013). Circulating microRNAs identified in a genome-wide serum microRNA expression analysis as noninvasive biomarkers for endometriosis. J. Clin. Endocrinol. Metab.98, 281–289. 10.1210/jc.2012-2415 98 WangL. G.SongH. Y.ZhangK. P.LiuL. M.JiangY. T.LiuS. (2022). LINC00536 promotes hepatocellular carcinoma progression via the miR-203b-5p/VEGFA axis. Neoplasma69, 136–144. 10.4149/neo_2021_210815n1158 99 WeberJ. A.BaxterD. H.ZhangS.HuangD. Y.HuangK. H.LeeM. J.et al (2010). The microRNA spectrum in 12 body fluids. Clin. Chem.56, 1733–1741. 10.1373/clinchem.2010.147405 100 WentgesT. H.El-ShorafaH. M.BeckmannJ.GabrielM.PoutanenM.GreveB.et al (2025). Combined targeting of TCF7L1/2, PTEN, CDK6, and BCCIP by microRNA miR-29c-3p is associated with reduced invasion and proliferation of endometriotic cells. Reprod. Med. Biol.24, e12645. 10.1002/rmb2.12645 101 WilsonM. R.ReskeJ. J.ChandlerR. L. (2022). AP-1 subunit JUNB promotes invasive phenotypes in endometriosis. Reprod. Sci.29, 3266–3277. 10.1007/s43032-022-00974-3 102 XiaB.LiuY.LiJ.JiangS. (2026). Mechanisms of KLF10 in regulating proliferation of endometriotic stromal cells in endometriosis. Gynecol. Obstet. Invest.91, 210–219. 10.1159/000546836 103 YadavV.SivaramA.VyasR. (2026). Exploring molecular signature and prognostic biomarkers in ovarian cancer: insights from Late-Stage, recurrent, and metastatic tumors. Biotechnol. Appl. Biochem. 10.1002/bab.70134 104 YangL. J.GaoL.GuoY. N.LiangZ. Q.LiD. M.TangY. L.et al (2021). Upregulation of microRNA miR-141-3p and its prospective targets in endometrial carcinoma: a comprehensive study. Bioengineered12, 2941–2956. 10.1080/21655979.2021.1943111 105 YangP. W.MaS. P.ChenX. J.YangF. M.WangS. M.WangQ.et al (2025). Luteolin suppresses cell migration and invasion via targeting miR-6809-5p/FLOT1/FAK and eliciting EMT in hepatocellular carcinoma. Transl. Oncol.61, 102511. 10.1016/j.tranon.2025.102511 106 YuY.WongW. H.HuY.ShenY.XuX.LuS.et al (2025). Translation of miRNA blood-based discovery to molecular testing for clinical diagnosis of endometriosis. Npj Women's Health3, 67. 10.1038/s44294-025-00116-5 107 ZafariN.BahramyA.Majidi ZolbinM.Emadi AllahyariS.FaraziE.HassannejadZ.et al (2022). microRNAs as novel diagnostic biomarkers in endometriosis patients: a systematic review and meta-analysis. Expert Rev. Mol. Diagn22, 479–495. 10.1080/14737159.2021.1960508 108 ZhangX.YaoH.HuangX.LuB.XuH.ZhouC. (2010). Nerve fibres in ovarian endometriotic lesions in women with ovarian endometriosis. Hum. Reprod.25, 392–397. 10.1093/humrep/dep427 109 ZhangQ.DuanJ.OlsonM.FazleabasA.GuoS. W. (2016). Cellular changes consistent with epithelial-mesenchymal transition and fibroblast-to-myofibroblast transdifferentiation in the progression of experimental endometriosis in baboons. Reprod. Sci.23, 1409–1421. 10.1177/1933719116641763 110 ZhangY.ZhengY.ZhuG. (2020). MiR-203a-3p targets PTEN to promote hepatocyte proliferation by regulating PI3K/Akt pathway in BRL-3A cells. Biosci. Biotechnol. Biochem.84, 725–733. 10.1080/09168451.2019.1694860 111 ZhangF.ChengN.DuJ.ZhangH.ZhangC. (2021). MicroRNA-200b-3p promotes endothelial cell apoptosis by targeting HDAC4 in atherosclerosis. BMC Cardiovasc Disord.21, 172. 10.1186/s12872-021-01980-0 112 ZubrzyckaA.Migdalska-SękM.JędrzejczykS.Brzeziańska-LasotaE. (2023). Assessment of BMP7, SMAD4, and CDH1 expression profile and regulatory miRNA-542-3p in eutopic and ectopic endometrium of women with endometriosis. Int. J. Mol. Sci.24, 6637. 10.3390/ijms24076637 Summary

endometriosis, inflammation, miRNA, non-invasive biomarkers, oxidative stress, saliva Citation Calugi G, Blandino F, Proserpio C, Forlani L, Meschino N, Giannico R, Dalla Pozza F, Graziani A, Di Pietro L, Buccarello L, Martire FG, Dealberti D, Lanzo G, Cosma S, Marchino G, Luisi S, Damiani GR, Garavaglia E, Landi S, Marziali M, Liprino A, Pontello VV, Fracchioli S and Caruso G (2026) Profiling miRNA in salivary samples from subjects with endometriosis: a pilot study. Front. Mol. Biosci. 13:1869937. doi: 10.3389/fmolb.2026.1869937 Received 30 April 2026 Revised 16 June 2026 Accepted 24 June 2026 Published 22 July 2026 Volume 13 - 2026 Edited by Daniel X. Zhang, Hong Kong Metropolitan University, Hong Kong SAR, China Reviewed by Laura Duran-Lozano, Vall d’Hebron Institute of Oncology (VHIO), Spain Sanjit Roy, Glycomontra, Inc, United States Updates Copyright © 2026 Calugi, Blandino, Proserpio, Forlani, Meschino, Giannico, Dalla Pozza, Graziani, Di Pietro, Buccarello, Martire, Dealberti, Lanzo, Cosma, Marchino, Luisi, Damiani, Garavaglia, Landi, Marziali, Liprino, Pontello, Fracchioli and Caruso. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms. *Correspondence: Giuseppe Caruso, [email protected] Disclaimer All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.