Front. Vet. Sci. Frontiers in Veterinary Science Front. Vet. Sci. 2297-1769 Frontiers Media S.A. 10.3389/fvets.2021.760487 Veterinary Science Original Research Circulating MicroRNAs as Non-invasive Biomarkers for Canine Cushing's Syndrome Sanders Karin 1 * Veldhuizen Anouk 1 Kooistra Hans S. 1 Slob Adri 1 Timmermans-Sprang Elpetra P. M. 1 Riemers Frank M. 1 Daminet Sylvie 2 Fracassi Federico 3 van Nimwegen Sebastiaan A. 1 Meij Björn P. 1 Galac Sara 1 1Department of Clinical Sciences, Faculty of Veterinary Medicine, Utrecht University, Utrecht, Netherlands 2Small Animal Department, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium 3Department of Veterinary Medical Sciences, University of Bologna, Ozzano dell'Emilia, Bologna, Italy

Edited by: David Bruyette, Anivive Lifesciences, United States

Reviewed by: Corrin John Boyd, Murdoch University, Australia; Elena De Felice, University of Camerino, Italy

 *Correspondence: Karin Sanders k.sanders@uu.nl

This article was submitted to Comparative and Clinical Medicine, a section of the journal Frontiers in Veterinary Science

 17 11 2021 2021 8 760487 18 08 2021 21 10 2021 Copyright © 2021 Sanders, Veldhuizen, Kooistra, Slob, Timmermans-Sprang, Riemers, Daminet, Fracassi, van Nimwegen, Meij and Galac. 2021 Sanders, Veldhuizen, Kooistra, Slob, Timmermans-Sprang, Riemers, Daminet, Fracassi, van Nimwegen, Meij and Galac

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.

Canine Cushing's syndrome (hypercortisolism) can be caused by a pituitary tumor (pituitary-dependent hypercortisolism; PDH) or a cortisol-secreting adrenocortical tumor (csACT). For both cases, non-invasive biomarkers that could pre-operatively predict the risk of recurrence after surgery would greatly impact clinical decision making. The aim of this study was to determine whether circulating microRNAs (miRNAs) can be used as diagnostic (presence of PDH or csACT) and/or prognostic (disease recurrence, histological grade) non-invasive biomarkers for canine Cushing's syndrome. After a pilot study with 40 miRNAs in blood samples of healthy dogs (n = 3), dogs with PDH (n = 3) and dogs with a csACT (n = 4), we selected a total of 20 miRNAs for the definitive study. In the definitive study, these 20 miRNAs were analyzed in blood samples of healthy dogs (n = 6), dogs with PDH (n = 19, pre- and post-operative samples) and dogs with a csACT (n = 26, pre-operative samples). In dogs with PDH, six miRNAs (miR-122-5p, miR-126-5p, miR-141-3p, miR-222-3p, miR-375-3p and miR-483-3p) were differentially expressed compared to healthy dogs. Of one miRNA, miR-122-5p, the expression levels did not overlap between healthy dogs and dogs with PDH (p = 2.9x10−4), significantly decreased after hypophysectomy (p = 0.013), and were significantly higher (p = 0.017) in dogs with recurrence (n = 3) than in dogs without recurrence for at least one year after hypophysectomy (n = 7). In dogs with csACTs, two miRNAs (miR-483-3p and miR-223-3p) were differentially expressed compared to healthy dogs. Additionally, miR-141-3p was expressed significantly lower (p = 0.009) in dogs with csACTs that had a histopathological Utrecht score of ≥ 11 compared to those with a score of <11. These results indicate that circulating miRNAs have the potential to be non-invasive biomarkers in dogs with Cushing's syndrome that may contribute to clinical decision making.

 miRNA pituitary-dependent hypercortisolism Cushing's disease adrenocortical tumor miR-122 prognostic recurrence normalization Introduction

Spontaneous Cushing's syndrome, or hypercortisolism, is one of the most commonly diagnosed endocrinopathies in dogs (1). It is caused by an ACTH-secreting pituitary tumor (pituitary-dependent hypercortisolism; PDH) in ~80–85% of cases, and by a cortisol-secreting adrenocortical tumor (csACT) in ~15–20% of cases (1). Both PDH and csACT can be treated by surgically removing the causative tumor. Because surgery is not without risks and not suitable for every patient, dogs with Cushing's syndrome are often treated with the steroidogenesis inhibitor trilostane (2). Although trilostane can effectively reduce the clinical signs associated with hypercortisolism, it does not inhibit tumor growth (3).

Pituitary tumors are usually classified as adenomas, but can nonetheless compress or invade surrounding tissues (4). After hypophysectomy, recurrence of hypercortisolism occurs in ~23–27% of cases (5, 6), for which the most important predictive marker is the pituitary height to brain area value (P/B value) (6). To determine the P/B value, diagnostic imaging with computed tomography (CT) or magnetic resonance imaging (MRI) is necessary. The mRNA expression of pituitary tumor transforming gene 1 (PTTG1) was previously reported to be associated with disease-free interval after surgery, but can only be assessed post-operatively (7).

A csACT can be classified as an adrenocortical adenoma (ACA) or an adrenocortical carcinoma (ACC) (8), but making this distinction can be difficult (9). After adrenalectomy, recurrence of hypercortisolism occurs in ~30–38% of cases (9, 10), but predicting recurrence remains challenging. We have recently developed a novel histopathological scoring system, the Utrecht score, to predict the prognosis of dogs after adrenalectomy (9). In addition, we have identified three genes of which the mRNA expression in csACTs was associated with survival after adrenalectomy (11). However, these markers can only be assessed when the csACT has already been surgically removed. The only known pre-operative marker that is associated with survival after adrenalectomy is the tumor diameter, but has low predictive value (9, 12).

Having non-invasive biomarkers that can pre-operatively predict the risk of recurrence after surgery would greatly impact clinical decision making. Potential candidates for such non-invasive biomarkers are circulating microRNAs (miRNAs; miRs). MiRNAs are single-stranded, non-coding RNAs of ~20–24 nucleotides in size (13). They function as antisense RNAs which regulate their target genes post-transcriptionally, and can influence cellular differentiation, proliferation and apoptosis (13, 14). MiRNA expression patterns can be altered in multiple diseases, including cancer (15, 16) and Cushing's syndrome in humans (14, 17, 18). MiRNAs that are expressed higher in cancer are regarded as oncomiRs, while miRNAs that are expressed lower in cancer are regarded as tumor suppressor miRNAs (13). Most miRNAs are expressed intracellularly, but numerous miRNAs can also be found in biological fluids such as blood, urine, saliva and cerebrospinal fluid. These miRNAs are referred to as circulating miRNAs (19). MiRNAs that are detectable in plasma or serum samples are bound with lipid proteins or encapsulated in extracellular vesicles and are therefore resistant to RNase digestions (13). Consequently, their expression levels are remarkably stable, which makes them ideal candidates as non-invasive biomarkers for various diseases (13).

The aim of this study was to determine whether circulating miRNAs can be used as non-invasive biomarkers for canine Cushing's syndrome. We hypothesized that we could identify circulating miRNAs of which the expression profiles are associated with the presence or absence of PDH or a csACT, the size of the pituitary tumor and whether PDH recurred after hypophysectomy, and the histopathological assessment of a csACT.

 Materials and Methods Pilot Study Animals and Samples

For the pilot study, we included blood samples of healthy dogs (n = 3), dogs with csACTs (n = 4), and dogs with PDH (n = 3). Written informed consent was obtained from the owners for the participation of their animals in this study. Serum samples were available from the healthy dogs and dogs with csACTs, while EDTA plasma samples were available from the dogs with PDH. After sample collection, serum or EDTA plasma was obtained by centrifuging whole blood for 5 min at 4,000 rpm. The serum and EDTA plasma samples were stored at −70°C until use.

The suspicion of hypercortisolism was based on the dogs' medical history and findings during physical examination. The diagnosis of PDH was made by demonstration of suppressible hypercortisolism with endocrine testing (low-dose dexamethasone suppression test or urinary corticoid to creatinine ratios (UCCR) combined with high-dose dexamethasone suppression test), or non-suppressible hypercortisolism combined with findings of symmetric bilateral adrenal enlargement (indicating bilateral hyperplasia) and the absence of unilateral adrenal enlargement (indicating an ACT) on ultrasonography or CT, together with findings of pituitary enlargement on CT or MRI. The diagnosis of a csACT was made by non-suppressible hypercortisolism demonstrated with endocrine testing, combined with the presence of an ACT found on ultrasonography or CT. All pituitary tumors were classified as adenomas based on histopathological assessment and absence of detectable metastases (4). Histopathological assessment of the ACTs confirmed their adrenocortical origin, and the histopathological Utrecht score (Ki67 proliferation index (PI) + 4 if ≥ 33% of cells have clear/vacuolated cytoplasm + 3 if necrosis is present) was used to predict the risk of recurrence (9). The dogs in the “healthy” group were regarded as healthy based on the absence of clinical signs, no abnormalities in complete blood count and blood biochemistry, and UCCRs within reference values.

 MiRNA Isolation and CDNA Synthesis

MiRNA was isolated using the miRNeasy Serum/Plasma Kit (ID: 217184, Qiagen, Venlo, the Netherlands), approved for animal and human plasma and serum samples, according to the manufacturer's instructions. The serum and EDTA plasma samples were thawed on ice. As sample input 150 μL was used, and volumes of the other components were adjusted accordingly. The miRNeasy serum/plasma kit was combined with the RNA Spike-In Kit, For RT (ID: 339390, Qiagen) for miRNA isolation quality control. The mixture of spike-ins UniSp2, UniSp4 and UniSp5 was prepared according to the manufacturer's instructions, and 1 μL of the spike-in mixture was added per sample. After isolation, the miRNA eluates were stored at −20°C until use. For subsequent cDNA synthesis, miRNA samples were thawed on ice. The cDNA was synthesized using the miRCURY LNA RT Kit (ID: 339340, Qiagen) according to the manufacturer's instructions, using 2 μL of RNA template per 10 μL total reaction volume. For cDNA synthesis quality control, 0.5 μL of a mixture of spike ins UniSp6 and cel-miR-39-3p (Qiagen) was added to the cDNA reaction. The cDNA samples were stored at 4°C when real-time quantitative polymerase chain reaction (RT-qPCR) was performed within four days after cDNA synthesis, or at −20°C when the RT-qPCR was performed at a later timepoint.

 RT-QPCR: Quality Control

For all samples, a quality control was performed to detect the presence of spike-ins UniSp2, UniSp4 and UniSp5 (to assess efficiency of miRNA isolation), of UniSp6 and cel-miR-39-3p (to assess efficiency of cDNA synthesis), and of the miRNAs miR-23a and miR-191 (to assess efficiency of endogenous miRNAs detection). To obtain miRNA-specific primers we checked the canine sequences on the miRBase database (20), and ordered miRCURY LNA miRNA PCR Assays [Qiagen; sequences and category numbers available through online data repository (21)]. The cDNA samples were diluted 1:15, and 3 μL of diluted cDNA was used per well in a total reaction volume of 10 μL. Detection of targets was performed with miRCURY LNA SYBR® Green PCR Kits (ID: 339347, Qiagen) according to the manufacturer's instructions in a MyiQ™2 Two-Color Real-Time PCR Detection System (Bio-Rad, Veenendaal, The Netherlands).

 RT-QPCR: Custom Plates

For the pilot study, miRCURY LNA miRNA Custom PCR Panels 96-well plates were ordered from Qiagen. These ready-to-use plates contained specific primer sets pre-coated in the wells. For each sample, 48 reactions were performed: UniSp3 (interplate calibrator for optimal calibration and control assessment), UniSp6 (cDNA spike-in control), one non-template control, five potential reference miRNAs, and 40 potential target miRNAs [Qiagen; sequences and category numbers available through online data repository (21)]. Detection of targets was performed with miRCURY LNA SYBR® Green PCR Kits (Qiagen) according to the manufacturer's instructions in a MyiQ™2 Two-Color Real-Time PCR Detection System (Bio-Rad). The threshold for detection of fluorescence was manually adjusted to 40 Relative Fluorescent Units to obtain the same threshold for each plate. The geometric mean of the five potential reference miRNAs was used to normalize the data. The relative expression of the target miRNAs was calculated using the 2−ΔΔCT method (22).

 Definitive Study Animals

For the definitive study, we included blood samples of healthy dogs (n = 6), pre-operative blood samples of dogs with csACTs (n = 26), and both pre- and post-operative samples of dogs with PDH (n = 19). The blood samples used in the pilot study were also included but reanalyzed in the definitive study. For the healthy dogs both EDTA plasma and serum samples were available. For the dogs with PDH only EDTA plasma samples were available, while for the dogs with csACTs only serum samples were available. Because of these differences in sample types, miRNA expression levels were not compared between dogs with PDH and with csACTs, but only between dogs with either PDH or csACTs and healthy controls of corresponding sample types. Written informed consent was obtained from the owners for the participation of their animals in this study. For all dogs with PDH, surgery and sample collection were performed at the Faculty of Veterinary Medicine in Utrecht, as was the sample collection for the healthy dogs. For dogs with csACTs, the surgery and sample collection were performed at the Utrecht University Clinic or at external clinics. In case of external clinics, the samples were sent by post to the Faculty of Veterinary Medicine in Utrecht. Samples for the dogs with PDH and csACTs were collected between 2015 and 2020, samples for the healthy dogs were collected between 2019 and 2020. The clinical data of the dogs included in this study can be found in an online data repository (21).

Diagnoses of PDH and csACTs were made as described for the pilot study, as were the confirmations of healthy dogs. Recurrence of hypercortisolism after hypophysectomy was suspected based on recurrence of clinical signs, and confirmed by UCCR values above the reference range.

 RT-QPCR

Isolation of miRNAs from the blood samples and subsequent cDNA synthesis was performed as described for the pilot study. Quality control was first performed on all samples to determine efficiency of miRNA isolation, cDNA reaction, and detection of endogenously expressed miRNAs. For the definitive study, we pre-coated 384-well plates with a mixture containing the PCR primer mix [Qiagen; sequences and category numbers available through online data repository (21)], miRCURY LNA SYBR® Green PCR Kits (Qiagen), and nuclease-free water. The plates were stored at −20°C until further use. The assessed targets in the definitive study included five potential reference miRNAs and 15 potential target miRNAs, selected based on the pilot study. Detection of miRNAs was performed in a CFX384 Touch Real-Time PCR Detection system (Bio-Rad).

 Data Analysis

To determine the pairwise variance and stability of miRNA expression, the geNorm (23) method from the SLqPCR package [v1.52.0, (24)] using R [v3.6, (25)] and RStudio [v1.3.1093, (26)] was used. Data normalization was performed by subtracting the geometric mean of the CT values of the 12 most stably expressed miRNAs from the CT values of the target miRNAs (ΔCT). The relative expression of the target miRNAs was calculated using the 2−ΔΔCT method (22). Normal distribution was assessed with the Shapiro-Wilk test. Because the data were not normally distributed, the Mann-Whitney U-test was performed to determine the significance of differences in expression between groups for independent samples, while paired sample analysis was performed with the Wilcoxon signed-rank test. The Spearman's rank correlation coefficient test was used to assess correlation between variables. P-values <0.05 were considered significant. All statistical analyses were performed with SPSS Statistics for Windows (Version 27.0, IBM Corp, Armonk, NY, USA).

The datasets generated for this study can be found in the DataverseNL data repository (21).

 Results Pilot Study Quality Control

The quality control results [available through online data repository (21)] indicated that highly (mimicked by UniSp2), moderately (UniSp4), and lowly (UniSp5) expressed miRNAs could all be efficiently isolated and detected with our protocol, and that highly (UniSp6) and moderately (cel-miR-39-3p) expressed miRNAs were efficiently transcribed to cDNA. In addition, the detection of miR-23 and miR-191 in all the samples indicated that endogenously present miRNAs were also efficiently detected.

 Target MiRNAs

Of the 40 target miRNAs assessed in the pilot study [results available through online data repository (21)], 6 were not detected in any samples (miR-34a-5p, miR-96-5p, miR-144-3p, miR-210-3p, miR-300-5p and miR-381-3p), while another 8 were not detected in most samples (miR-139-3p, miR-183-5p, miR-433-3p, miR-499a-5p, miR-455-5p, miR-483-5p, miR-497-5p, miR-499-5p). From the remaining 26 miRNAs, 15 miRNAs were selected for further analyses. The selected miRNAs showed potential differences in expression between patient groups or were potentially associated with the P/B value in the PDH group or with histopathological assessment in the csACT group.

 Definitive Study Stably Expressed MiRNAs

Five miRNAs were initially included in this study as normalization controls because we expected their expression levels to be stable: miR-23-3p, miR-191-5p, miR-222-3p, miR-423-5p and miR-425-5p. The selection of these miRNAs was based on endogenous miRNAs that are typically detected and stably expressed in human serum/plasma samples according to the Qiagen Guidelines for Profiling Biofluid miRNAs ©2019. We used geNorm software (23) to determine whether these miRNAs were indeed stably expressed in EDTA plasma samples (n = 44; 6 healthy dogs and 19 dogs with PDH, pre-op and post-op), in serum samples (n = 32; 6 healthy dogs and 26 dogs with csACT), and in the combined EDTA plasma/serum samples. The combination of miRNAs was interpreted as suitably stable when the Pairwise Variance (V-score) was lower than 0.15 (V0.15) (23, 27). However, although the V-score decreased when all 5 miRNAs were used compared to when 3 or 4 miRNAs were used, they did not reach the V0.15 (Figure 1A).

Pairwise Variance (V-score) (A,B) and Average expression stability M (C) of miRNAs included in geNorm analyses to determine the optimal number of miRNAs for data normalization. Results when including 5 (A; miR-23a-3p, miR-191-5p, miR-222-3p, miR-423-5p and miR-425-5p) and when including 16 (B,C; miR-7-5p, miR-16-5p, miR-18a-5p, miR-21-5p, miR-23a-3p, miR-122-5p, miR-126-5p, miR-191-5p, miR-194-5p, miR-222-3p, miR-223-3p, miR-320a-3p, miR-375-3p, miR-423-5p, miR-425-5p and miR-483-3p) miRNAs. Left panels (EDTA plasma) include samples from healthy dogs (n = 6) and from dogs with pituitary-dependent hypercortisolism (n = 19, pre-op and post-op samples); middle panels (serum) include samples from healthy dogs (n = 6) and from dogs with cortisol-secreting adrenocortical tumors (n = 26); right panels (all) include both EDTA plasma and serum samples (n = 76).

To determine whether the V0.15 would be reached when more miRNAs were included in the analysis, we assessed the stability of all 20 assessed miRNAs, including the 15 intended target miRNAs. Only the miRNAs that could be detected in all samples were included, which left 16 miRNAs for the geNorm analysis (5 of the originally intended reference miRNAs and 11 target miRNAs). The lowest V-score (Vmin; 0.115) in the combined EDTA plasma/serum samples group (“All”) was achieved with a combination of 12 miRNAs (Figure 1B). The most stably expressed miRNAs in both EDTA plasma and serum sample groups was miR-191-5p, while the least stably expressed miRNAs in all sample groups were miR-122-5p, miR-483-3p, and miR-375-3p (Figure 1C). Overall, both the V-scores (Figure 1B) and the Average expression stability M (Figure 1C) were lower in the EDTA plasma samples group than in the serum samples group.

For the subsequent analyses, we analyzed the expression levels of all 20 miRNAs. The geometric mean of the Vmin (12 most stably expressed miRNAs) was used as normalization control.

 MiRNAs in Pituitary-Dependent Hypercortisolism

In the EDTA plasma samples, five miRNAs were expressed significantly higher in dogs with PDH (n = 19) compared to healthy dogs (n = 6): miR-122-5p (p = 2.9x10−4; no overlap between groups), miR-141-3p (p = 0.028), miR-222-3p (p = 0.008), miR-375-3p (p = 0.001), and miR-483-3p (p = 0.009) (Figure 2A). One miRNA was expressed significantly lower in dogs with PDH compared to healthy dogs: miR-126-5p (p = 0.036) (Figure 2A).

Circulating miRNA expression levels in EDTA plasma samples, normalized with Vmin (12 most stably expressed miRNAs), in (A) dogs with PDH (n = 19) and healthy dogs (n = 6); (B) dogs with PDH (n = 19) before (pre-op) and after (post-op) hypophysectomy; and (C) dogs with PDH with (n = 3) or without (n = 7) recurrence of hypercortisolism within a year after hypophysectomy. Circles above and below boxes indicate outliers. * p < 0.05; ** p < 0.01; *** p < 0.001. PDH, pituitary-dependent hypercortisolism.

Of all 19 dogs with PDH, post-op samples were available. The timepoint of post-op sample collection was at median 3 days after hypophysectomy (range 1–7 days). The six differentially expressed miRNAs were analyzed in the post-op EDTA plasma samples, to determine whether their expression levels either decreased or increased after surgery. Paired-sample analyses of these six miRNAs showed that the expression levels of miR-122-5p (p = 0.013) and of miR-141-3p (p = 0.035) significantly decreased after surgery (Figure 2B). Although miR-222-3p and miR-375-3p showed a tendency to decrease based on Figure 2B, these changes were not significant (p = 0.117 and p = 0.165, respectively). The expression levels of miR-126-5p and miR-483-3p did not change in the post-op samples (p = 0.647 and p = 0.520, respectively).

Association of MiRNAs With Recurrence of PDH

Because the P/B value is currently the most important predictor of recurrence after hypophysectomy, we analyzed whether expression levels of miRNAs (pre-operative samples) were correlated with the P/B value. None of the six differentially expressed miRNAs were significantly correlated with the P/B value, nor were any of the other fourteen miRNAs (Table 1).

Correlation of pre-operative circulating miRNA expression with pituitary height to brain area value in dogs with PDH.

Correlation Significance Number (n)
coefficient (r) (2-tailed; p)
miR-7-5p 0.207 0.395 19
miR-16-5p 0.044 0.857 19
miR-18a-5p −0.029 0.908 19
miR-21-5p −0.083 0.737 19
miR-23a-3p −0.201 0.409 19
miR-122-5p 0.179 0.463 19
miR-126-5p −0.309 0.198 19
miR-132-3p −0.151 0.563 17
miR-141-3p 0.145 0.554 19
miR-191-5p 0.051 0.836 19
miR-194-5p 0.193 0.428 19
miR-218-5p −0.071 0.772 19
miR-222-3p 0.175 0.474 19
miR-223-3p 0.035 0.887 19
miR-320a-3p 0.247 0.307 19
miR-375-3p 0.067 0.786 19
miR-423-5p −0.116 0.637 19
miR-425-5p 0.105 0.668 19
miR-483-3p 0.023 0.926 19
miR-503-5p 0.068 0.788 18

Data analyzed with the Spearman's rank correlation coefficient test.

Of the 19 dogs with PDH, follow-up information was available of 14 dogs. Of these dogs, 3 had recurrence of hypercortisolism (5, 7 and 12 months after hypophysectomy), 7 dogs had no recurrence of hypercortisolism and had follow-up of at least one year after surgery (median follow-up time 21 months, range 16 – 50 months), and 4 dogs had no recurrence of hypercortisolism but follow-up times of less than one year after surgery (median 5 months, range 2 – 8 months). The 3 dogs with reported recurrence were subsequently included in the “Recurrence” group, while the 7 dogs without recurrence and with follow-up time of at least one year after surgery were included in the “No recurrence” group. In assessing the expression profiles of the six differentially expressed miRNAs in the pre-operative samples, we found that miR-122-5p (p = 0.017) and miR-222-3p (p = 0.030) were expressed significantly higher in the dogs with recurrence than in the dogs without recurrence (Figure 2C).

 Differentially Expressed MiRNAs With One Reference MiRNA

Because the clinical application of miRNAs is limited when 12 miRNAs need to be assessed for data normalization, we wanted to determine whether the detected differences could still be observed when only one reference miRNA was used. As the most stably expressed miRNA in all sample groups, we chose to use the expression of miR-191-5p for data normalization. Although miR-126-5p lost its significance (p = 0.105), the other five miRNAs were still significantly higher in dogs with PDH compared to healthy dogs (Figure 3A) when normalized to only miR-191-5p. MiR-122-5p was equally significant when only one reference miRNA was used (p = 2.9x10−4) because still no overlap between sample groups was observed, but the p-values of the other four miRNAs increased (miR-141-3p: p = 0.049 vs. 0.028; miR-222-3p: p = 0.026 vs. 0.008; miR-375-3p: p = 0.003 vs. 0.001; miR-483-3p: p = 0.022 vs. 0.009).

Circulating miRNA expression levels in EDTA plasma samples, normalized with only miR-191-5p, in (A) dogs with PDH (n = 19) and healthy dogs (n = 6); and (B) dogs with PDH with (n = 3) or without (n = 7) recurrence of hypercortisolism within a year after hypophysectomy. Circles below boxes indicate outliers. * p < 0.05; ** p < 0.01; *** p < 0.001. PDH, pituitary-dependent hypercortisolism.

When comparing the expression of these miRNAs between the recurrence group and the no recurrence group, normalized with only miR-191-5p, both miR-122-5p (p = 0.017) and miR-222-3p (p = 0.017) were still significantly higher in the recurrence group (Figure 3B).

 MiRNAs in Adrenal-Dependent Hypercortisolism

In serum samples, one miRNA was expressed significantly higher in dogs with a csACT (n = 26) compared to healthy dogs (n = 6): miR-483-3p (p = 0.020), while one miRNA was expressed significantly lower in dogs with a csACT: miR-223-3p (p = 0.038) (Figure 4A).

Circulating miRNA expression levels in serum samples, normalized with Vmin (12 most stably expressed miRNAs), in (A) dogs with a csACT (n = 26) and healthy dogs (n = 6); and (B) dogs with a csACT with an Utrecht score below 11 (n = 9) compared to those with a score of 11 or more (n = 17). Circles above and below boxes indicate outliers. * p < 0.05; *** p < 0.001. csACT, cortisol-secreting adrenocortical tumor.

 Histopathological Assessment

The median Ki67 PI in the csACTs (n = 26) was 8.2 % (range 0.3–50.7 %), and the median Utrecht score was 12.6 (range 4.3 – 57.7). Although miR-132-3p and miR-141-3p showed a moderate correlation with the Utrecht score, either positively (miR-132-3p) or negatively (miR-141-3p), these correlations were not significant (Table 2). None of the other miRNAs correlated to either the Ki67 PI or Utrecht score (Table 2).

Correlation of circulating miRNA expression with Ki67 proliferation index (PI) and Utrecht score in dogs with csACTs.

Ki67 PI Ki67 PI Utrecht score Utrecht score Number
Correlation coefficient (r) P-value Correlation coefficient (r) P-value (n)
miR-7-5p 0.027 0.896 0.015 0.942 26
miR-16-5p −0.047 0.822 0.030 0.883 26
miR-18a-5p −0.037 0.858 −0.066 0.750 26
miR-21-5p −0.167 0.414 −0.034 0.870 26
miR-23a-3p −0.024 0.909 −0.093 0.651 26
miR-122-5p −0.175 0.391 −0.063 0.759 26
miR-126-5p −0.091 0.660 −0.064 0.758 26
miR-132-3p 0.274 0.243 0.406 0.076 20
miR-141-3p −0.188 0.415 −0.406 0.068 21
miR-191-5p 0.267 0.187 0.241 0.235 26
miR-194-5p −0.127 0.538 −0.062 0.763 26
miR-218-5p −0.121 0.565 −0.156 0.456 25
miR-222-3p 0.027 0.896 0.015 0.942 26
miR-223-3p −0.042 0.840 0.031 0.881 26
miR-320a-3p 0.162 0.428 0.026 0.901 26
miR-375-3p −0.091 0.657 −0.024 0.909 26
miR-423-5p 0.166 0.418 0.015 0.941 26
miR-425-5p 0.013 0.948 0.009 0.966 26
miR-483-3p −0.105 0.610 −0.187 0.360 26
miR-503-5p −0.031 0.885 −0.029 0.894 24

Data analyzed with the Spearman's rank correlation coefficient test.

In our previously published introduction of the Utrecht score, we identified three groups with significantly different survival times after surgery when using cut-off values for the Utrecht score of <6, ≥6– <11, and ≥ 11. In the current study, only one csACT had an Utrecht score of <6, so for subsequent analyses we classified the csACTs as having an Utrecht score of <11 (n = 9) or ≥ 11 (n = 17). When comparing miRNA expression between these groups, miR-141-3p was significantly higher in in the group with an Utrecht score of <11 than in the group with an Utrecht score of ≥ 11 (p = 0.009, Figure 4B). None of the other miRNAs showed significant differences in expression between these groups (Figure 4B).

 Discussion

Circulating miRNAs have been shown to be useful diagnostic and prognostic biomarkers in several diseases and cancer types (13, 15). In this study, we identified several circulating miRNAs that are differently expressed in dogs with Cushing's syndrome. The most clearly altered miRNA was miR-122-5p, which was significantly overexpressed in dogs with PDH and did not overlap with expression in healthy dogs. In addition, miR-122-5p expression was higher in dogs with recurrence after hypophysectomy than in dogs without reported recurrence. MiR-122-5p might therefore represent a useful non-invasive biomarker for dogs with PDH.

MiR-1221 is also of interest in other diseases and cancer types, both in humans and in dogs (16, 28–31). In human cancers, the function of miR-122 seems to be dependent on its cellular context: while it has been described as an oncomiR in some cancers, including clear cell renal carcinoma (28) and colon cancer (31), it has been described as a tumor suppressor miRNA in others, including hepatocellular carcinoma (32) and bladder cancer (33). In several studies, miR-122 was found to inhibit the expression of aldolase, fructose-bisphosphate A (ALDOA) (31, 34–36). ALDOA is a glycolytic enzyme that catalyzes the reversible conversion of fructose-1,6-bisphosphate to glyceraldehyde 3-phosphate and dihydroxyacetone phosphate. Although ALDOA is an important oncogene in several cancer types, including non-small cell lung cancer (37), pancreatic cancer (38) and hepatocellular carcinoma (39), it is a tumor suppressor in others, including prostate adenocarcinoma and stomach adenocarcinoma (39). The tumorigenic activity of ALDOA therefore seems to be context-dependent (39), and could potentially explain the context-dependency of miR-122. The expression of ALDOA has not yet been studied in pituitary tumors, which would be an interesting avenue for future research.

Although miR-122 expression has a clear link with cancer, there are also other diseases in which miR-122 expression can be dysregulated. Because miR-122 is derived from hepatocytes, it was found to be overexpressed in serum of dogs with different types of liver diseases (29, 30, 40). Because hypercortisolism can induce vacuolar hepatopathy (41), the increase in miR-122 levels seen in dogs with PDH could also be related to liver damage. In future studies, it would be interesting to determine whether miR-122 levels in dogs with PDH are correlated to their ALT values, as they were in a previous study on dogs with liver diseases (40). The finding that miR-122-5p levels in dogs with PDH significantly decreased in post-op samples compared to pre-op samples, even though the time between surgery and post-op sampling was relatively short (median 3 days, range 1-7 days) and the dogs still receive high doses of glucocorticoids for the first couple of days after hypophysectomy, suggests that its expression in dogs with PDH is not related to hypercortisolism but rather to the pituitary tumorigenesis. In addition, miR-122-5p levels were not significantly higher in dogs with hypercortisolism caused by a csACT than in healthy dogs. To validate the hypothesis that miR-122-5p expression is not hypercortisolism-related, it would be interesting to compare the miRNA expression levels that were found to be altered in this study in dogs with PDH to those of dogs with iatrogenic Cushing's syndrome.

Another miRNA that was overexpressed in dogs with PDH compared to healthy dogs, as well as in dogs with PDH that had recurrence after hypophysectomy compared to dogs that did not have recurrence, is miR-222-3p. MiR-222 is overexpressed in many different human cancer types (42). A recent study that performed meta-analyses on 17 published articles on miR-222 showed that high expression of miR-222 significantly predicts poor overall survival in several human cancer types (43). MiR-222 could therefore also be a useful noninvasive biomarker in dogs with PDH.

Of the other four significantly altered miRNAs in PDH found in this study, only one was expressed significantly lower compared to healthy controls: miR-126-5p. MiR-126 has also been reported to be downregulated in other diseases in humans, including in diabetes mellitus (44). In contrast, circulating miR-126 was overexpressed in dogs with different types of neoplasia (45). Possibly, the miR-126-5p downregulation in the current study is related to the resulting endocrine syndrome, and not to the pituitary tumor. MiR-141-3p, miR-375-3p miR-483-3p were overexpressed in dogs with PDH and are also deregulated in several cancer types (46–48). For miR-375-3p, the difference in expression could also be related to the hypothalamic-pituitary-adrenal axis feedback loop: expression of miR-375 has been shown to be inhibited by corticotropin-releasing factor (CRF) (49). Since the high ACTH levels in dogs with PDH will result in low CRF levels (50), this lack of inhibition could result in higher miR-375p levels.

In dogs with csACTs, miR-483-3p was significantly overexpressed compared to healthy dogs. In humans, miR-483-5p, the reverse miRNA from the miR-483 locus, is one of the most commonly reported circulating miRNA that is overexpressed in patients with ACC compared to those with ACA (51–53). Although miR-483-3p is less studied, it is also overexpressed in ACC compared to ACA in humans (47). In the miRBase database, only the sequence for canine miR-483-3p is reported. Because of its importance in human ACTs, we nonetheless tested the primers for human miR-483-5p in our pilot study, but these did not result in a detectable product in the canine samples (21). Although the overexpression of miR-483-3p in dogs with csACTs compared to healthy dogs seems to be in line with its overexpression in human ACCs compared to ACAs, we did not see any differences in its expression when comparing csACTs with an Utrecht score of ≥ 11 with those that had an Utrecht score of <11, nor was there a correlation of miR-483-3p expression with Ki67 PI or Utrecht score. Whether this is related to a difference in classification (e.g., the fact that we only included cortisol-secreting ACTs, or a potential overrepresentation of carcinomas in the total csACT group) or a difference in mechanism is currently unknown.

MiR-141-3p was expressed significantly lower in dogs with a csACT that had an Utrecht score of ≥ 11 than in those with an Utrecht score of <11. MiR-141 has been reported as a useful biomarker in several cancer types, but was shown to be upregulated in some cancer types, including small cell lung cancer (46), and downregulated in others, including renal cell carcinoma (54). Recently, miR-141 was shown to inhibit angiogenesis both in vitro and in vivo (55). Previous research by our group showed that canine ACCs are in a more proangiogenic state compared to ACAs (56), which might be related to the decreased expression of miR-141-3p in the csACTs with an Utrecht score of ≥ 11. In future studies it would be interesting to determine whether miR-141-3p can predict the risk of recurrence after adrenalectomy.

MiRNAs that can predict malignancy or risk of recurrence before surgery would greatly impact clinical decision making. In dogs with PDH, both miR-122-5p and miR-222-3p were significantly overexpressed in dogs with recurrence after hypophysectomy compared to those without and are therefore interesting candidates as noninvasive biomarkers to predict recurrence. However, because only 3 dogs were included in the “Recurrence” group and 7 in the “No recurrence” group, these results will have to be validated in larger sample sizes. Additionally, in future studies it would be interesting to take blood samples at different times after surgery to determine the expression levels of these miRNAs at time of remission, and whether they will increase at time of recurrence before changes on diagnostic imaging are detectable. If so, these miRNAs could be highly useful early markers of recurrence.

The biggest hurdle in circulating miRNA analyses is data normalization. The five miRNAs that we originally intended to use as reference miRNAs were not stable enough, so we determined whether the stability would improve when we added all analyzed miRNAs in the geNorm analysis. Indeed, when adding more miRNAs, we were able to reach sufficiently stable results for data normalization. However, although this improves the reliability of our results in the current study, using 12 miRNAs for data normalization limits the clinical application of detecting circulating miRNAs. We therefore tested whether the differentially expressed miRNAs in dogs with PDH would retain their significance when normalized against only one miRNA. Although the p-values of some miRNAs slightly increased, in general the same results were observed as when using 12 miRNAs for normalization. Using only miR-191-5p, one of the two most stably expressed miRNAs in this study, as reference miRNA for data normalization therefore seems feasible in dogs with PDH. However, miR-191-5p can also be dysregulated in different cancer types (57, 58), which makes it risky to use only one miRNA for data normalization. Ideally, to identify miRNAs that are stably expressed in dogs with different diseases, large-scale experiments should be performed that analyze expression of all circulating miRNAs with next generation sequencing in groups of dogs with different types of diseases. This would also help to identify novel useful miRNAs.

Overall, the miRNAs seemed to be more stably expressed in the EDTA plasma samples than in the serum samples as shown in the geNorm analyses. A previous study showed that there can indeed be differences in stability in EDTA plasma compared to serum samples (59). However, in the current study, handling of samples after collection could also have influenced the results. For example, all EDTA plasma samples from dogs with PDH were collected in our University Clinic and could be processed soon after collection, while samples from the dogs with csACTs were in several cases sent in from other clinics, so it took longer before the samples were processed and frozen down. Because of the small differences in EDTA plasma compared to serum samples, we decided not to directly compare the dogs with PDH to dogs with csACTs. For future studies it would be interesting to collect blood of all dogs in the same sample types, so that direct comparisons can be made.

To conclude, we have identified several circulating miRNAs that are dysregulated in dogs with Cushing's syndrome. These miRNAs have the potential to be non-invasive biomarkers that may contribute to clinical decision making.

 Data Availability Statement

The datasets generated for this study will be made publicly available through the DataverseNL data repository.

 Ethics Statement

Ethical review and approval was not required for the animal study because no experimental interventions were carried out on the animals. Rest blood samples were collected when blood samples were taken for diagnostic purposes (dogs with PDH and csACTs) or preventive health screenings (healthy dogs). Written informed consent was obtained from the owners for the participation of their animals in this study.

 Author Contributions

KS and SG conceived and designed the experiments. KS, AV, AS, and ET-S performed the experiments. KS, AV, HK, FR, BM, and SG analyzed the data. HK, SD, FF, SN, and BM contributed materials. KS, AV, and SG wrote the paper. All authors contributed to the article and approved the final version of the manuscript for submission.

 Funding

We thank the Utrecht University Open Access Fund for partly reimbursing the Article Processing Fee.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 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.

We would like to thank Drs. Bart Sjollema, Rob Gerritsen, Giora van Straten, Jenny Buijtels and Dieneke Jongepier for contributing samples to our csACT dataset; Guy C. M. Grinwis for histopathological evaluations; all dog owners for letting us use their dogs' samples for this study; Louis C. Penning for sharing his insights and expertise; and Friends of VetMed for their support. A preprint of this manuscript has been published on biorxiv.org (https://doi.org/10.1101/2021.08.17.456598).

 References Galac S Reusch CE Kooistra HS Rijnberk A. Adrenals. In: Clinical Endocrinology of Dogs and Cats, eds. Rijnberk A. Kooistra H. S. (Schlütersche), p. 93–154. Sanders K Kooistra HS Galac S. Treating canine Cushing's syndrome: current options and future prospects. Vet J. (2018) 241:42–51. 10.1016/j.tvjl.2018.09.01430340659 Teshima T Hara Y Takekoshi S Nezu Y Harada Y Yogo T . Trilostane-induced inhibition of cortisol secretion results in reduced negative feedback at the hypothalamic-pituitary axis. Domest Anim Endocrinol. (2009) 36:32–44. 10.1016/j.domaniend.2008.10.00219041802 Sanders K Galac S Meij BP. Pituitary tumour types in dogs and cats. Vet J. (2021) 270:105623. 10.1016/j.tvjl.2021.10562333641809 Hanson JM Teske E Voorhout G Galac S Kooistra HS Meij BP. Prognostic factors for outcome after transsphenoidal hypophysectomy in dogs with pituitary-dependent hyperadrenocorticism. J Neurosurg. (2007) 107:830–40. 10.3171/JNS-07/10/083017937231 van Rijn SJ Galac S Tryfonidou MA Hesselink JW Penning LC Kooistra HS . The influence of pituitary size on outcome after transsphenoidal hypophysectomy in a large cohort of dogs with pituitary-dependent hypercortisolism. J Vet Intern Med. (2016) 30:989–95. 10.1111/jvim.1436727425149 Vastenhout N van Rijn SJ Riemers FM Tryfonidou MA Meij BP Penning LC. The mRNA expression of PTTG1 is a strong prognostic indicator for recurrence after hypophysectomy in dogs with corticotroph pituitary adenomas. Vet J. (2018) 240:19–21. 10.1016/j.tvjl.2018.08.01230268327 Labelle P Kyles AE Farver TB de Cock HE V. Indicators of malignancy of canine adrenocortical tumors: histopathology and proliferation index. Vet Pathol. (2004) 41:490–7. 10.1354/vp.41-5-49015347821 Sanders K Cirkel K Grinwis GCM Teske E van Nimwegen SA Mol JA . The utrecht score: a novel histopathological scoring system to assess the prognosis of dogs with cortisol-secreting adrenocortical tumours. Vet Comp Oncol. (2019) 17:329–37. 10.1111/vco.1247430860644 van Sluijs FJ Sjollema BE Voorhout G van den Ingh TS Rijnberk A. Results of adrenalectomy in 36 dogs with hyperadrenocorticism caused by adreno-cortical tumour. Vet Q. (1995) 17:113–6. 10.1080/01652176.1995.96945478525598 Sanders K van Staalduinen GJ Uijens MCM Mol JA Teske E Slob A . Molecular markers of prognosis in canine cortisol-secreting adrenocortical tumours. Vet Comp Oncol. (2019) 17:545–52. 10.1111/vco.1252131301217 Schwartz P Kovak JR Koprowski A Ludwig LL Monette S Bergman PJ. Evaluation of prognostic factors in the surgical treatment of adrenal gland tumors in dogs: 41 cases (1999–2005). J Am Vet Med Assoc. (2008) 232:77–84. 10.2460/javma.232.1.7718167113 Wang H Peng R Wang J Qin Z Xue L. Circulating microRNAs as potential cancer biomarkers: the advantage and disadvantage. Clin Epigenetics. (2018) 10:59. 10.1186/s13148-018-0492-129713393 Vetrivel S Zhang R Engel M Altieri B Braun L Osswald A . Circulating microRNA expression in cushing's syndrome. Front Endocrinol. (2021) 12:620012. 10.3389/fendo.2021.62001233692756 Condrat CE Thompson DC Barbu MG Bugnar OL Boboc A Cretoiu D . miRNAs as biomarkers in disease: latest findings regarding their role in diagnosis and prognosis. Cells. (2020) 9:276. 10.3390/cells902027631979244 Rupaimoole R Slack FJ. MicroRNA therapeutics: Towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. (2017) 16:203–21. 10.1038/nrd.2016.24628209991 Belaya Z Khandaeva P Nonn L Nikitin A Solodovnikov A Sitkin I . Circulating plasma microRNA to differentiate cushing's disease from ectopic ACTH syndrome. Front Endocrinol. (2020) 11:331 10.3389/fendo.2020.0033132582027 Bimpaki EI Iliopoulos D Moraitis A Stratakis CA. MicroRNA signature in massive macronodular adrenocortical disease and implications for adrenocortical tumourigenesis. Clin Endocrinol (Oxf). (2010) 72:744–51. 10.1111/j.1365-2265.2009.03725.x19849700 Beylerli O Beeraka NM Gareev I Pavlov V Yang G Liang Y . Mirnas as noninvasive biomarkers and therapeutic agents of pituitary adenomas. Int J Mol Sci. (2020) 21:1–15. 10.3390/ijms2119728733023145 Kozomara A Birgaoanu M Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. (2019) 47:D155–62. 10.1093/nar/gky114130423142 Sanders K Veldhuizen A Kooistra HS Slob A Timmermans-Sprang EPM Riemers FM . Circulating microRNAs in dogs with Cushing's syndrome. DataverseNL. (2021). 10.34894/XPIIU6 Livak KJ Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method. Methods. (2001) 25:402–8. 10.1006/meth.2001.126223525470 Vandesompele J De Preter K Pattyn I Poppe B Van Roy N De Paepe A . Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. (2002) 3:34–1. 10.1186/gb-2002-3-7-research003412184808 Kohl M. SLqPCR: functions for analysis of real-time quantitative PCR data at SIRS-Lab GmbH. (2007). Core Team R. R: a language and environment for statistical computing. (2020) Available online at: https://www.r-project.org/ RStudio Team R. RStudio: Integrated Development Environment for R. (2021) Available online at: http://www.rstudio.com/ Xiong F Cheng X Zhang C Klar RM He T. Optimizations for identifying reference genes in bone and cartilage bioengineering. BMC Biotechnol. (2021) 21:25. 10.1186/s12896-021-00685-833731065 Nie W Ni D Ma X Zhang Y Gao Y Peng C . miR-122 promotes proliferation and invasion of clear cell renal cell carcinoma by suppressing Forkhead box O3. Int J Oncol. (2019) 54:559–71. 10.3892/ijo.2019.469430720052 Dirksen K Verzijl T Grinwis GC Favier RP Penning LC Burgener IA . Use of serum MicroRNAs as biomarker for hepatobiliary diseases in dogs. J Vet Intern Med. (2016) 30:1816–23. 10.1111/jvim.1460227859748 Oosthuyzen W Ten Berg PWL Francis B Campbell S Macklin V Milne E . Sensitivity and specificity of microRNA-122 for liver disease in dogs. J Vet Intern Med. (2018) 32:1637–44. 10.1111/jvim.1525030070004 Li H Zhang X Jin Z Yin T Duan C Sun J . MiR-122 promotes the development of colon cancer by targeting ALDOA In Vitro. Technol Cancer Res Treat. (2019) 18:1–10. 10.1177/153303381987130031564215 Cheng D Deng J Zhang B He X Meng Z Li G . LncRNA HOTAIR epigenetically suppresses miR-122 expression in hepatocellular carcinoma via DNA methylation. EBioMedicine. (2018) 36:159–70. 10.1016/j.ebiom.2018.08.05530195653 Wang Y Xing QF Liu XQ Guo ZJ Li CY Sun G. MiR-122 targets VEGFC in bladder cancer to inhibit tumor growth and angiogenesis. Am J Transl Res. (2016) 8:3056–66. 27508026 Cui K Jin S Du Y Yu J Feng H Fan Q . Long noncoding RNA DIO3OS interacts with miR-122 to promote proliferation and invasion of pancreatic cancer cells through upregulating ALDOA. Cancer Cell Int. (2019) 19:202. 10.1186/s12935-019-0922-y31384177 Zeng Y Lv Y Tao L Ma J Zhang H Xu H . G6PC3, ALDOA and CS induction accompanies mir-122 downregulation in the mechanical asphyxia and can serve as hypoxia biomarkers. Oncotarget. (2016) 7:74526–36. 10.18632/oncotarget.1293127793029 Xu Z Liu G Zhang M Zhang Z Jia Y Peng L . MiR-122-5p inhibits the proliferation, invasion and growth of bile duct carcinoma cells by targeting ALDOA. Cell Physiol Biochem. (2018) 48:2596–606. 10.1159/00049270230121648 Fu H Gao H Qi X Zhao L Wu D Bai Y . Aldolase A promotes proliferation and G 1 /S transition via the EGFR/MAPK pathway in non-small cell lung cancer. Cancer Commun. (2018) 38:18. 10.1186/s40880-018-0290-329764507 Ji S Zhang B Liu J Qin Y Liang C Shi S . ALDOA functions as an oncogene in the highly metastatic pancreatic cancer. Cancer Lett. (2016) 374:127–35. 10.1016/j.canlet.2016.01.05426854714 Ma D Chen X Zhang PY Zhang H Wei LJ Hu S . Upregulation of the ALDOA/DNA-PK/p53 pathway by dietary restriction suppresses tumor growth. Oncogene. (2018) 37:1041–8. 10.1038/onc.2017.39829084207 Dirksen K Verzijl T van den Ingh TSGAM Vernooij JCM van der Laan LJW Burgener IA . Hepatocyte-derived microRNAs as sensitive serum biomarkers of hepatocellular injury in Labrador retrievers. Vet J. (2016) 211:75–81. 10.1016/j.tvjl.2016.01.01027021912 Gilor C Graves TK. Interpretation of laboratory tests for canine cushing's syndrome. Top Companion Anim Med. (2011) 26:98–108. 10.1053/j.tcam.2011.03.00121596349 Song J Ouyang Y Che J Li X Zhao Y Yang K . Potential value of miR-221/222 as diagnostic, prognostic, and therapeutic biomarkers for diseases. Front Immunol. (2017) 8:56. 10.3389/fimmu.2017.0005628261196 Ravegnini G Cargnin S Sammarini G Zanotti F Bermejo JL Hrelia P . Prognostic role of miR-221 and miR-222 expression in cancer patients: a systematic review and meta-analysis. Cancers. (2019) 11:970. 10.3390/cancers1107097031336701 Zampetaki A Kiechl S Drozdov I Willeit P Mayr U Prokopi M . Plasma MicroRNA profiling reveals loss of endothelial MiR-126 and other MicroRNAs in type 2 diabetes. Circ Res. (2010) 107:810–7. 10.1161/CIRCRESAHA.110.22635720651284 Heishima K Ichikawa Y Yoshida K Iwasaki R Sakai H Nakagawa T . Circulating microRNA-214 and−126 as potential biomarkers for canine neoplastic disease. Sci Rep. (2017) 7:2301. 10.1038/s41598-017-02607-128536479 Mao S Lu Z Zheng S Zhang H Zhang G Wang F . Exosomal miR-141 promotes tumor angiogenesis via KLF12 in small cell lung cancer. J Exp Clin Cancer Res. (2020) 39:193. 10.1186/s13046-020-01680-132958011 Decmann A Perge P Nyíro G Darvasi O Likó I Borka K . MicroRNA expression profiling in adrenal myelolipoma. J Clin Endocrinol Metab. (2018) 103:3522–30. 10.1210/jc.2018-0081729982598 Huang X Yuan T Liang M Du M Xia S Dittmar R . Exosomal miR-1290 and miR-375 as prognostic markers in castration-resistant prostate cancer. Eur Urol. (2015) 67:33–41. 10.1016/j.eururo.2014.07.03525129854 Zhang N Lin JK Chen J Liu XF Liu JL Luo HS . MicroRNA 375 mediates the signaling pathway of corticotropin-releasing factor (CRF) regulating pro-opiomelanocortin (POMC) Expression by targeting Mitogen-activated protein Kinase 8. J Biol Chem. (2013) 288:10361–73. 10.1074/jbc.M112.42550423430746 Van Wijk PA Rijnberk A Croughs RJM Voorhout G Sprang EPM Mol JA. Corticotropin-releasing hormone and adrenocorticotropic hormone concentrations in cerebrospinal fluid of dogs with pituitary-dependent hyperadrenocorticism. Endocrinology. (1992) 131:2659–62. 10.1210/endo.131.6.13328441332844 Decmann A Bancos I Khanna A Thomas MA Turai P Perge P . Comparison of plasma and urinary microRNA-483-5p for the diagnosis of adrenocortical malignancy. J Biotechnol. (2019) 297:49–53. 10.1016/j.jbiotec.2019.03.01730940435 Patel D Boufraqech M Jain M Zhang L He M Gesuwan K . MiR-34a and miR-483-5p are candidate serum biomarkers for adrenocortical tumors. Surg (United States). (2013) 154:1224–9. 10.1016/j.surg.2013.06.02224238045 Igaz P Igaz I Nagy Z Nyíro G Szabó PM Falus A . MicroRNAs in adrenal tumors: relevance for pathogenesis, diagnosis, and therapy. Cell Mol Life Sci. (2015) 72:417–28. 10.1007/s00018-014-1752-725297921 Chen X Wang X Ruan A Han W Zhao Y Lu X . MiR-141 Is a key regulator of renal cell carcinoma proliferation and metastasis by controlling EphA2 expression. Clin Cancer Res. (2014) 20:2617–30. 10.1158/1078-0432.CCR-13-322424647573 Dong H Weng C Bai R Sheng J Gao X Li L . The regulatory network of miR-141 in the inhibition of angiogenesis. Angiogenesis. (2019) 22:251–62. 10.1007/s10456-018-9654-130465119 Kool MMJ Galac S Kooistra HS Mol JA. Expression of angiogenesis-related genes in canine cortisol-secreting adrenocortical tumors. Domest Anim Endocrinol. (2014) 47:73–82. 10.1016/j.domaniend.2013.11.00424377872 Zhou L Zhang F Tong J Liu F. MiR-191-5p inhibits lung adenocarcinoma by repressing SATB1 to inhibit Wnt pathway. Mol Genet Genomic Med. (2020) 8:e1043. 10.1002/mgg3.104331724324 Ashirbekov Y Abaildayev A Omarbayeva N Botbayev D Belkozhayev A Askandirova A . Combination of circulating miR-145-5p/miR-191-5p as biomarker for breast cancer detection. PeerJ. (2020) 8:e10494. 10.7717/peerj.1049433362968 Glinge C Clauss S Boddum K Jabbari R Jabbari J Risgaard B . Stability of circulating blood-based MicroRNAs – pre-analytic methodological considerations. PLoS One. (2017) 12:e0167969. 10.1371/journal.pone.016796928151938

1When−3p or−5p miRNA was not specified in mentioned articles, only the miR-locus number will be mentioned.