- 1Department of Molecular Biosciences, The Wenner-Gren Institute, Stockholm University, Stockholm, Sweden
- 2Centre for Paleogenetics, Stockholm University, Stockholm, Sweden
- 3Department of Veterinary Clinical Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
- 4Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Introduction: Non-infectious inflammatory diseases of the central nervous system in dogs, such as steroid responsive meningitis-arteritis (SRMA) and meningoencephalitis of unknown origin (MUO), represent a common clinical challenge that needs extensive and multimodal work-up to reach a presumptive diagnosis. Both diseases are presumably caused by dysregulations of the immune system, but further research is needed in order to understand the molecular mechanisms behind each disease and to optimize treatment.
Methods: By next-generation sequencing and subsequent quantitative real-time PCR (qPCR) verification, we designed a prospective case–control pilot study to analyze the small RNA profiles of cerebrospinal fluid from dogs suffering from MUO (N = 5), dogs suffering from SRMA (N = 8), and healthy dogs (N = 5) presented for elective euthanasia used as the Control group.
Results: Our results showed an overall enrichment in Y-RNA fragments across all samples, followed by microRNAs (miRNAs) and ribosomal RNAs as the major findings. Additional traces of short RNA reads mapped to long non-coding RNAs and protein-coding genes were also found. From the detected canine miRNAs, miR-21, miR-486, miR-148a, miR-99a, miR-191 and miR-92a were among the most abundant. Dogs with SRMA showed higher differences in miRNA abundance than dogs with MUO when compared to healthy dogs, and miR-142-3p was consistently detected as differentially upregulated in both diseases, although at a low concentration. Moreover, miR-405-5p and miR-503-5p showed different profiles between SRMA and MUO dogs. Subsequent qPCR analyses confirmed miR-142-5p, miR-191-5p and miR-92a-3p as significantly upregulated miRNAs in dogs with SRMA and/or MUO.
Discussion: Cerebrospinal fluid is a challenging biological material to use for profiling miRNAs due to the low content of circulating RNAs. Despite this, we could confirm several miRNAs being differentially abundant when comparing healthy dogs and dogs with MUO and SRMA, respectively. The results of this study indicate a potential role of miRNAs in the underlying molecular mechanisms of these diseases and establish the basis for further studies.
1. Introduction
Meningoencephalitis of unknown origin (MUO) and steroid responsive meningitis-arteritis (SRMA) in dogs are both common non-infectious inflammatory diseases, affecting the meninges and/or the central nervous system (CNS) (1, 2). For both diseases, the underlying etiology is presumed to be an immune system dysregulation, but other than that, they are different disease entities, with MUO having a graver prognosis than SRMA (3). MUO is an umbrella term covering several encephalitides that can currently only be differentiated by histopathology (4). Clinical phenotypes for the different types of encephalitides overlap greatly from a clinical perspective, and studies examining MUO often do not differentiate dogs into pathological subgroups, hence MUO is used as a common term to cover very similar clinical phenotypes (1, 5, 6). MUO is usually diagnosed based on pathological changes on neuroimaging, i.e., magnetic resonance imaging (MRI) or computed tomography (CT), and cerebrospinal fluid (CSF) analysis. However, in some patients, both neuroimaging and CSF analysis are normal, making an ante-mortem diagnosis extremely challenging (5). For all subtypes, an early diagnosis and initiation of treatment is, however, of great importance, as early initiation of treatment has been shown to affect survival (7). Treatment is based on general immunosuppressive therapy, usually involving high dose corticosteroids and other immunomodulatory drugs in combination, often introducing a number of unwanted side effects, but as the specific underlying immunopathological mechanisms are largely unknown, targeted treatment toward these is not an option at present (6). A better understanding of the disease mechanisms is therefore warranted. SRMA, in contrast to MUO, is a systemic inflammatory disease affecting the meninges and meningeal arteries (8, 9), and usually responds well to treatment with corticosteroids (10). Dogs are commonly presented with cervical hyperesthesia, anorexia, reluctance to move and systemic inflammation (10). In the absence of specific serum biomarkers, the mainstay in diagnosing SRMA is analysis of CSF, which requires an invasive sampling procedure under general anesthesia (2, 11).
MicroRNAs (miRNAs) are a group of small non-coding RNAs that can be found in body fluids such as serum, CSF and saliva (12, 13). MiRNAs are known to play a post-transcriptional regulatory role in molecular pathways involved in multiple pathological processes, including non-infectious inflammatory CNS diseases. In humans, a relationship between the expression of some miRNAs in inflammation and neurodegenerative diseases such as Alzheimer’s and Parkinson’s diseases has been convincingly demonstrated (12, 14, 15). However, miRNA research in dogs is sparse so far. Previous studies have investigated selected miRNAs by probe-based hypothesis-driven methods, i.e., NanoString or microarray platforms, in a small number of diseases, including neurological and cardiac pathologies (16, 17). Research on miRNA abundance profiles in CSF in dogs with inflammatory neurological disease is also limited (16), and has proven to be challenging (18). Despite this, miRNAs remain of great interest, as they appear more stable in circulation than traditional biomarkers (19–22), qualifying their use and possible implementation in clinical research (18, 19). More importantly, the identification of changes in the abundance of certain circulating miRNAs in CSF from dogs suffering from MUO or SRMA has the potential to improve our knowledge of the mechanisms involved in each disease, hence ameliorating the treatment of affected patients.
This pilot study aimed to investigate CSF-miRNA abundance profiles by hypothesis-free methods, i.e., small RNA sequencing, and subsequent verification by qPCR, in order to shed light on the molecular mechanisms that characterize SRMA and MUO, and to potentially reveal novel targets for treatment.
2. Materials and methods
2.1. Animals
The study was performed as a prospective case–control study. CSF was collected from a total of 15 dogs allocated to three different groups: (i) MUO group (N = 5), (ii) SRMA group (N = 5), and iii) Control group (N = 5). Three additional dogs suffering from SRMA were included in the qPCR verification analyses.
For the MUO and SRMA groups, dogs with characteristic clinical and paraclinical findings as deemed by the responsible clinician at the Neurology service at the University Hospital for Companion Animals (Copenhagen, Denmark) were included. For the dogs with MUO, a pleocytosis in CSF dominated by mononuclear cells following general recommendations for a tentative diagnosis (23) was required for inclusion, supported by diagnostic imaging or necropsy when available. Besides the characteristic clinical presentation, dogs included in the SRMA group were required to have a pleocytosis in CSF, defined as >5 nucleated cells/μL dominated by neutrophilic granulocytes or monocytes as generally advised (23), as well as a systemic inflammatory response, defined as a C-reactive protein (CRP) above normal reference range (> 25 mg/L). In the SRMA and MUO groups, a positive response to treatment with immunosuppressive drugs was supportive of the diagnosis. Dogs were excluded from the two disease groups if they suffered from non-related inflammatory or systemic diseases, if a CSF tap was contraindicated and/or if they had been treated with nonsteroidal anti-inflammatory drugs or corticosteroids for >7 days prior to presentation. For the healthy Control group, dogs presented for elective euthanasia with no clinical neurological or systemic disease were included, provided they had a CSF analysis within standard reference values. This was defined as a total nucleated cell count (TNCC) ≤ 5 cells/μL, and a protein count <30 mg/dL for cerebellomedullary samples (23–25). In all groups, dogs were excluded if there were any signs of blood contamination on visual inspection of CSF or if the amount of red blood cells (RBC) compromised the cytological evaluation and the TNCC, defined as a maximum of 8,480 RBC/μL as previously reported for canine CSF (26). A comprehensive list of all dogs analyzed in the present study, as well as clinical results and follow-up from the analysis of their sampled CSF with CRP and TNCC measurements (in SRMA and MUO groups) is available at Table 1.
2.2. Sample collection and handling
Cerebrospinal fluid was collected lege artis from the cerebellomedullary cistern for all analyzed groups. In the two disease groups (MUO and SRMA), surplus CSF from the clinical work-up was included. In the event of death or euthanasia before CSF collection could be performed, CSF was collected within 30 min of euthanasia and included for routine analysis and miRNA investigations. In all healthy dogs, CSF was collected postmortem within 30 min of euthanasia. The time limit of 30 min from death was set to avoid postmortem changes to the CSF composition, as also applied in previous human studies (27).
Standard analysis of CSF for clinical verification was performed within 1 h of collection. The analysis included macroscopic inspection of color and turbidity, microscopic investigation of cell morphology, and analysis of protein level and nucleated cell content. The remaining CSF was centrifuged at 2000 g for 15 min at 4°C (Multifuge 1 S-R, Heraeus). Aliquots of the supernatant were then transferred to sterile cryotubes at a maximum volume of 200 μL/tube and stored at −80°C until RNA purification.
2.3. Purification of RNA and sequencing
RNA was purified from 200 μL of stored CSF samples using the miRNeasy Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The quality and quantity of the purified RNA was evaluated using the Nanodrop 1,000 spectrophotometer (ThermoFisher Scientific). Six μL of RNA from each sample (N = 15) were shipped to the Center for Genomic Regulation (CRG), Genomic Service at Barcelona, Spain, for sequencing. Small RNA-seq libraries were prepared with the NEBNext® Small RNA Library Prep kit (New England Biolabs) and sequenced on an Illumina HiSeq 2,500 system to generate 50 bp single-end reads.
2.4. Small RNA-seq pre-processing, mapping, and quantification
Raw FASTQ files were processed for sequencing adapter trimming with the Cutadapt v3.2 software (28) and allowing a maximum of 10% error rate for adapter identification, as well as a minimum read length after trimming of at least 18 nucleotides (nt). Quality check filtering was performed with the fastp v0.12.4 software tool (29), allowing a qualified PHRED score per nucleotide ≥30. Reads fulfilling quality check requirements were then mapped against canine precursor miRNA sequences (pre-miRNAs, N = 502) belonging to the domestic dog genome assembly (CanFam3.1) according to the miRBase v22.1 database (30). Alignment was performed with the Bowtie tool v1.3.0 (31) using the following specifications tailored for small RNA sequences: i) no mismatches allowed, (ii) reporting the best alignment with high sensitivity within the stacked multimapping repertoire, and (iii) removing any reverse-complement match (−v 0−k 1−y --best -no-rc). Quantification of miRNA abundance was assessed independently for each 5p and 3p mature miRNA ends from each pre-miRNA detected.
In addition, we aimed to characterize the presence of RNA fragments belonging to transcripts derived from loci categories other than miRNAs (i.e., mRNAs, ribosomal RNAs, long non-coding RNAs or Y-RNAs). For this purpose, quality-filtered reads were mapped to the whole-genome assembly of the domestic dog (ROS_Cfam_1.0) by using the Bowtie v1.3. 0 aligner as implemented for miRNA-targeted mapping. We allowed up to a maximum of 1 mismatch with a mapping seed of 18 nt, equal to the minimum read length allowed (−n 1−l 18−k 1−y --best). Quantification was then performed with the featureCounts v2.0.3 tool (32) and focused on the exonic fraction of the successfully assigned reads. Gene annotations were retrieved from the latest Ensembl release available when performing the analyses (Cfam_1.0 v.105) (33).
2.5. Differential abundance
The significance of differences in miRNA abundance for Control vs. MUO dogs, Control vs. SRMA dogs, and MUO vs. SRMA dogs was assessed by using the edgeR tool (34). MiRNAs with extremely low abundance were tagged with the filterByExpr R function and removed from differential abundance analyses. Abundance-filtered raw counts were then normalized for library depth with the trimmed mean of M-values normalization (TMM) method (35). Statistical significance of mean abundance differences was tested with a quasi-likelihood F-test (34). Multiple hypothesis testing correction was implemented with the false discovery rate (FDR) method (36). Significant differential abundance was set at an absolute fold-change (FC) value >2, i.e., |log2FC| > 1 in the log2 scale, and FDR corrected p-value (q-value) < 0.05. Relevant miRNAs showing significant differential abundance and, whenever possible, high abundance, were selected for qPCR assessment. A comprehensive pipeline of bioinformatics approaches for RNA pre-processing, mapping, quantification, and differential abundance analyses is shown in Supplementary Figure S1.
2.6. Validation of candidate miRNAs using qPCR
Three additional dogs were recruited for the SRMA group after the sequencing part of the project was finished and were included for qPCR verification. RNA extraction for these three additional samples was done as described above. No additional RNA was available for dog M2 from the MUO group after sequencing, hence it was excluded from qPCR analyses. Two μL of RNA from each sample were used for cDNA synthesis according to methods reported previously (37). For each RNA sample, two replicates of cDNA were done when enough RNA was available (for dogs M1, M3 and M4 from the MUO group, as well as for dogs S1, S2, S4 and S5 from the SRMA group, only 1 replicate could be done). Negative technical controls (NTC) without adding cDNA to the reaction mix, as well as control samples with cDNA synthesized but not including poly(A) polymerase (noPAP) were also included. The cDNAs were finally diluted eight times and used for qPCR procedure.
A panel of 18 miRNAs were included in the verification phase. Specific forward and reverse primers were designed for each miRNA selected according to the miRprimer software (38) following protocols described previously (39). Primer sequences are shown in Supplementary Table S1. qPCR was done in a Mx 3005P machine (Agilent) using QuantiFast SYBR Green Master Mix (Qiagen) with the following cycling conditions: Five minutes at 95°C, 40 cycles of 10 s at 95°C and 30 s at 60°C, followed by melting curve analysis to ensure specific amplification. Data was manually curated and assays with no Cq (quantification cycle) values or values >33 cycles were excluded. NTC and nonPAP control samples were visually inspected and considered valid if showing no amplification or amplification with Cq at least >5 cycles away from CSF samples. Data pre-processing was done using GenEx Pro v.6 software (MultiID, Sweden). Briefly, data was normalized to the most stable miRNA (cfa-let-7a) according to the NormFinder algorithm (40). Technical repeats were then averaged. Subsequently, relative quantities were calculated by scaling to the lowest abundant sample for each assay. Finally, data was log2 transformed to run statistical analyses. A t-test with Welch correction was performed for the following comparisons: Control vs. MUO, Control vs. SRMA and MUO vs. SRMA. Multiple testing correction was implemented by using the FDR approach (36). MiRNAs with FC > 2 and q-value <0.05 were considered as significantly deregulated by qPCR technique. Raw data from qPCR and subsequent intermediate processing steps are shown in Supplementary Table S2.
2.7. Tissue clustering
A miRNA tissue atlas from the domestic dog was used to infer whether the miRNA abundance profiles obtained from CSF in our study resembled that from any of the available canine reference tissues and their homogeneity. In this way, an overall equivalent miRNA abundance among the CSF samples would result in an adjacent clustering of CSF miRNA profiles among each other. Moreover, similar miRNA abundance patterns to other tissues would be reflected in their proximity to the corresponding samples from the dog miRNA atlas, thus indicating the putative tissue origin of the miRNAs present in CSF. Small RNA-seq samples for building the dog tissue atlas were retrieved (41, 42) and mapped against the reference domestic dog miRNA annotation according to miRGeneDB 2.1 (43). To allow full equivalence, small RNA-seq data from CSF obtained in this study was also mapped to dog primary miRNAs (pri-miRNAs) available at miRGeneDB 2.1 (43). The miRNA quantification of each considered reference dog tissue (hypothalamus, cerebellum, brain, sciatic nerve, bone marrow, pancreas, lung, liver kidney, plasma, heart, and skeletal muscle), as well as of CSF samples, was normalized within tissue by transforming individual counts to counts-per-million (CPM) values with respect to the total reads successfully mapped to the whole canine genome (ROS_Cfam_1.0). The umap function from the umap v0.2.7.0 R package1 was then implemented to determine a dimensionality reduction of the reference dog tissue atlas (n_neighbors = 15, metric = “pearson,” spread = 15, random_state = 123) using the uniform manifold approximation and projection (UMAP) algorithm (44). Finally, the CSF miRNA profiles were projected onto the previously learned UMAP embedding.
2.8. MiRNA target prediction
The TargetScan v8.0 (45) with H. sapiens reference database was employed to predict putative highly conserved mRNA targets for three of the miRNAs (cfa-miR-142-5p, cfa-miR-191-5p and cfa-miR-92a-3p) detected as significantly upregulated in the CSF pathological groups (MUO and/or SRMA) by both small RNA-seq and qPCR techniques. Since the seed region (2nd to 8th nt from the 5′ end of the mature miRNA sequence) of the selected miRNAs is highly conserved among human, domestic dog and many other mammals, we prioritized highly conserved mRNA targets as predicted by the TargetScan algorithm. In this way, we sorted the predicted mRNA targets for each of the three miRNAs by their aggregate PCT score, which represents an approximation to the probability of a given miRNA-mRNA interaction to have highly conserved functionality across mammals (46). Only predicted mRNA targets with PCT score > 0.4 were used for further analyses.
2.9. Pathway enrichment
The putative mRNA targets predicted for cfa-miR-142-5p, cfa-miR-191-5p and cfa-miR-92a-3p miRNAs according to TargetScan v8.0 (45) were then subjected to pathway enrichment analysis using the ClueGO v.2.5.7 plug-in application (47) within the Cytoscape v.3.6 (48) software. We implemented the KEGG pathway human database (v.2020) as reference, jointly with a right-sided hypergeometric test for gene enrichment calculation. Pathways were considered as significantly enriched if they were represented by at least five predicted targeted mRNA genes by the set of selected miRNAs, and a q-value <0.05 after multiple testing correction with the false discovery rate (FDR) approach (36).
3. Results
3.1. Dog population and clinical manifestation
The current pilot study included three different groups of dogs (MUO, N = 5; SRMA, N = 5 for sequencing and N = 8 for subsequent qPCR analysis; and healthy Control dogs, N = 5). A variety of breeds and weights were represented in the included cohort of dogs (Table 1), with dogs in the MUO group having a lower average weight compared to dogs suffering from SRMA, as well as healthy Control dogs. This was not a surprising finding, as MUO is known to affect small sized breeds in particular (5). The CRP levels for SRMA dogs were 20-fold higher than in MUO dogs (Supplementary Figure S2), while TNCC was highly variable in SRMA dogs, with S4 having ~14-fold higher TNCC than the average of the remaining SRMA dogs (Table 1). When disregarding S4 sample from the SRMA group, considered as an extreme outlier for the TNCC phenotype, SRMA dogs had a median TNCC value of 133 TNCC/μL, while MUO dogs showed a median of 299 TNCC/μL (Supplementary Figure S2). Moreover, SRMA dogs had an increased standard deviation for TNCC (s = 592.46) compared to MUO dogs (s = 516.739, Supplementary Figure S2).
3.2. Small RNA-seq analysis
An average of 15.7 million reads were successfully generated for each CSF sample after sequencing, except for one of the healthy Control samples (C3), which failed during library preparation and was hence excluded from further analyses based on small RNA-seq data (MUO, N = 5; SRMA, N = 5; and Control, N = 4). From the total number of raw RNA reads produced per sample, an average of 90.43% were successfully trimmed and passed the minimum length filtering criterion (at least 18 nt after adapter trimming), and 96.85% of these successfully passed quality-check filtering.
Genome-wide mapping to the domestic dog assembly (ROS_Canfam_1.0) produced a variable number of successful alignments, ranging between 28.35 and 82.29%. From these, between 3.9 and 46.5% of the aligned reads mapped to annotated genes in the dog assembly. Roughly half of the reads mapped to annotated loci belonged to Y-RNAs, followed by miRNAs, ribosomal RNAs (rRNAs), long non-coding RNAs (lncRNAs) and traces of mRNA transcripts from protein-coding genes (Figure 1A). A comprehensive list and abundance profiles of all quantified loci genome-wide in each defined CSF group is shown in Supplementary Table S3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144084/fvets-10-1144084-HTML/image_m/fvets-10-1144084-g001.jpg)
Figure 1. (A) Pie charts depicting the proportion of mapped reads from small RNA-seq data to different types of annotated loci in healthy dogs (Control), dogs with meningoencephalitis of unknown origin (MUO) and dogs with steroid responsive meningitis-arteritis (SRMA). Among the most abundant loci represented are long non-coding RNAs (lncRNAs), miRNAs, protein-coding genes, ribosomal RNAs (rRNAs) and Y-RNAs. (B) Volcano plots depicting miRNAs significantly upregulated (log2FC > 1; q-value <0.05, in blue) or downregulated (log2FC < −1; q-value <0.05, in red) for dogs with MUO or SRMA with respect to healthy Control dogs according to small RNA-seq data. For the MUO vs. SRMA contrast, upregulated miRNAs (in blue) were considered as those showing increased abundance in the SRMA group with respect to the MUO group, and vice versa. In grey are miRNAs showing nominal p-value <0.05. In black are miRNAs showing nominal p-value >0.05.
For reads mapped to dog miRNAs according to the miRBase reference (N = 502), between 1.29 million and 29 thousand reads were successfully assigned. Principal component analysis (PCA) based on miRNA abundance for each of the 14 CSF samples included initially revealed that sample S4 (belonging to the SRMA pathological group) showed a clear separation over principal component 1 (92.45% of explained variance, Supplementary Figure S3A), while the other samples clustered more closely. Upon analyzing the overall clustering of the PCA without including sample S4, which was previously considered as an extreme outlier for TNCC (Supplementary Figure S3B), we decided to remove this sample for further analyses, retaining 4 Control, 4 SRMA and 5 MUO CSF samples for differential abundance assessment.
3.3. miRNA differential abundance
Differential abundance analyses were performed by using CSF from healthy Control samples as reference; hence, any obtained upregulation would imply an overexpression of the given miRNA in either SRMA or MUO group with respect to the Control group, and vice versa. Control samples (N = 4) were compared separately with SRMA (N = 4) and MUO samples (N = 5). Additionally, SRMA and MUO samples were compared independently to assess the presence of discriminative differentially abundant miRNAs between both diseases. In this way, MUO samples (N = 5) were considered as reference, meaning that any given miRNA upregulation would imply an overexpression of the miRNA in the SRMA group with respect to the MUO group, and vice versa.
Only one very lowly abundant miRNA (cfa-miR-142-3p) was detected as significantly upregulated in MUO dogs compared to healthy Control dogs (|log2FC| > 1, q-value <0.05; Supplementary Table S4, in grey, and Figure 1B, in blue), while a total of 36 miRNAs were significantly differentially abundant in Control vs. SRMA (26 upregulated and 10 downregulated in SRMA dogs compared to healthy dogs, Supplementary Table S4, in grey, and Figure 1B). Of these, 19 (52.77%) showed average abundance levels above 100 read counts. Regarding the contrast between MUO and SRMA samples, three lowly abundant miRNAs were detected as differentially upregulated in SRMA dogs (cfa-miR-450b-5p, cfa-miR-450a-5p and cfa-miR-503-5p, Supplementary Table S4, in grey, and Figure 1B, in blue). A complete list of differential abundance analyses for all analyzed miRNAs by small RNA-seq is available in Supplementary Table S4.
After qPCR verification, two significantly upregulated miRNAs (cfa-miR-191-5p and cfa-miR-142-5p) were detected in MUO dogs with respect to Control healthy dogs (|FC| > 2, q-value <0.05, Table 2 and Figure 2A). Their abundance change agreed with that obtained in small RNA-seq analyses, although using sequencing technique they were only detected as significant at nominal p-value <0.05 (Supplementary Table S4). Interestingly, three miRNAs (cfa-miR-191-5p, cfa-miR-142-5p and cfa-miR-92a-3p) were also detected as differentially upregulated by qPCR when comparing SRMA dogs with Control healthy dogs (Table 2; Figure 2A), from which two of them (cfa-miR-191-5p and cfa-miR-142-5p) were also significantly upregulated in the MUO group by qPCR, and all three were significantly differentially abundant according to small RNA-seq data in the Control vs. SRMA contrast (Supplementary Table S4, in grey). No significant differences were found when comparing MUO with SRMA dogs by qPCR. A complete list of differential abundance analyses for all analyzed miRNAs by qPCR is available in Supplementary Table S5.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144084/fvets-10-1144084-HTML/image_m/fvets-10-1144084-t002.jpg)
Table 2. Differentially abundant miRNAs (|FC| > 2; q-value <0.05) detected by qPCR when comparing Control vs. MUO dogs and Control vs. SRMA dogs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144084/fvets-10-1144084-HTML/image_m/fvets-10-1144084-g002.jpg)
Figure 2. (A) Barplots depicting qPCR log2 transformed relative quantities and their significance levels for cfa-miR-142-5p, cfa-miR-191-5p and cfa-miR-92a-3p miRNAs measured in cerebrospinal fluid (CSF) from healthy dogs (Control), dogs with meningoencephalitis of unknown origin (MUO) and dogs with steroid responsive meningitis-arteritis (SRMA). (B) UMAP plot depicting sample clustering of a collection of tissues from a canine miRNA expression atlas including miRNA profiles of hypothalamus, cerebellum, brain, sciatic nerve, bone marrow, pancreas, lung, liver kidney, plasma, heart, and skeletal muscle. The sample clustering for miRNA profiles from CSF samples based on small RNA-seq data of healthy dogs (Control, N = 4), dogs with MUO (N = 5) and dogs with SRMA (N = 4) was then predicted using the learned UMAP embedding using the initial dog miRNA expression atlas.
3.4. Pathway enrichment analysis
Among the three queried miRNAs that showed significant abundance difference between healthy and disease states in CSF by small RNA-seq and qPCR (cfa-miR-142-5p, cfa-miR-191-5p and cfa-miR-92a-3p), two of them (cfa-miR-142-5p and cfa-miR-92a-3p) gathered significantly enriched pathways after mRNA target prediction (Supplementary Table S6). A good proportion of the pathways highlighted as significant were related to immune response, such as TGF-β signaling pathway (KEGG:04350), MAPK signaling pathway (KEGG:04010), bacterial infection (KEGG:05135, KEGG:05100), leukocyte migration (KEGG:04670), FoxO signaling pathway (KEGG:04068), cAMP signaling pathway (KEGG:04024) or Hippo signaling pathway (KEGG:04392). All predicted highly conserved and putatively targeted mRNAs used for enrichment analyses are shown in Supplementary Tables S7-S9.
3.5. Tissue clustering
We analyzed the putative tissue/cell origin of the miRNAs quantified in the dog CSF samples by using a reference dog miRNA atlas with a total of 12 different tissues (i.e., hypothalamus, cerebellum, brain, sciatic nerve, bone marrow, pancreas, lung, liver, kidney, plasma, heart and skeletal muscle) and our small RNA sequencing data. After tissue clustering based on miRNA expression profiles, the projection of CSF miRNA profiles using the UMAP algorithm resulted in an overall positioning of CSF samples as an independent cluster compared to other tissues, although two of them (M5 and S1, belonging to MUO and SRMA group, respectively), were located adjacent to dog plasma (Figure 2B).
4. Discussion
In the present study, we have investigated the profiles of cell-free small non-coding RNAs in canine CSF, with a specific focus on miRNAs, and included three experimental groups of dogs: a Control group consisting of healthy dogs, and two disease groups with dogs suffering from either MUO or SRMA. We first used small RNA sequencing on CSF samples, which provides higher sensitivity to quantify low abundant transcripts than qPCR, followed by subsequent verification of the best discriminating and abundant candidates by qPCR. Our goal in this pilot study was to shed light on the underlying immune mechanisms of SRMA and MUO.
To our knowledge, only a few studies have profiled miRNAs in CSF of dogs with MUO (49) or with MUO and SRMA (18). By using solely qPCR technique, Gaitero and coworkers identified higher levels of expression of miR-21 and miR-181c in the CSF from MUO dogs compared to dogs with non-inflammatory neurological diseases (49). Furthermore, they found a positive correlation between CSF cellularity and CSF expression of miR-21 in the MUO group (49). An additional study profiled miRNAs by qPCR using CSF of dogs with neoplastic, inflammatory, and degenerative disorders affecting the CNS (50). The authors found a significantly higher abundance of miR-10b-5p in the neoplastic group compared to other groups, but no correlation between miRNA expression and CSF cellularity nor CSF protein content. Moreover, in a previous report of our group using serum and CSF samples, we investigated the abundance of several miRNAs in dogs with MUO (N = 7) and SRMA (N = 6) using qPCR, where we managed to confirm the presence of several miRNAs, including miR-21, but struggled to reproduce consistent results in CSF (18). In the current study, we were only able to detect cfa-miR-142-3p as significantly differentially abundant in MUO dogs compared to healthy Control dogs using small RNA sequencing data. Moreover, qPCR verification highlighted cfa-miR-142-5p and cfa-miR-191-5p as significantly upregulated in MUO dogs, but both were only detected as significant at the nominal p-value with small RNA-seq data. However, cfa-miR-191-5p was detected as highly abundant according to small RNA sequencing, nearly 6-fold more abundant than cfa-miR-142-5p in MUO and SRMA dogs and might be a more interesting miRNA to be considered. For the dogs in the SRMA group, we were able to identify statistically significant levels of a total of 36 miRNAs using small RNA sequencing of CSF (see Results), and 3 out of the 12 miRNAs (cfa-miR-142-5p, cfa-miR-191-5p and cfa-miR-92a-3p) analyzed by qPCR were also significantly upregulated in SRMA dogs when compared to the Control group.
Focusing on the three differentially abundant miRNAs found in this study (cfa-miR-142-5p, cfa-miR-191-5p and cfa-miR-92a-3p), miR-142-5p has emerged as one of the most critical miRNAs during development, homeostasis, and disease, with important functions in infection and inflammation (51). Also, significant changes in its abundance have been found in plasma of Alzheimer’s disease compared to healthy subjects (52, 53). Moreover, the nuclear factor erythroid 2-like 2 (NRF2) is an experimentally validated target of miR-142-5p, and it is a relevant regulator involved in the response to injury and inflammation (53). Talebi and coworkers (54) showed that the miR-142 mature miRNA transcripts (−3p and -5p) in brain tissue of a mouse model for multiple sclerosis could target mRNAs encoding proteins involved in cytokine signaling and T cell differentiation, and confirmed that SOCS1 mRNA, a negative regulator of cytokine signaling, is also a direct target of miR-142a-5p.
Regarding miR-191-5p, its expression has been shown to be stable in human CSF from the first year of life (55), it is abundantly expressed in the brain, and it has been associated with cell differentiation, proliferation, and metastasis (56). It has been found upregulated in plasma of patients with prolonged disorders of consciousness after severe brain damage (57). Furthermore, it has been proposed as a promising serum biomarker in Alzheimer’s disease (58).
Finally, miR-92a-3p belongs to the miR-17-92 cluster and it has been identified as a biomarker for many human cancers, acting as a tumor-promoting regulator through the PI3K/AKT/mTOR signaling pathway (59). Additionally, it has been found downregulated in plasma from Alzheimer’s disease patients (60).
Pathway enrichment analyses of putative predicted mRNA targets for cfa-miR-142-5p and cfa-miR-92a-3p highlighted several signaling pathways closely associated with immune response regulation, such as TGF-β (61), MAPK (62), FoxO (63), cAMP (64) or Hippo (65), as well as pathways related to bacterial infection and leukocyte migration. These results might indicate a link between an active immune response elicited by different types of leukocytes in the CSF. This might be explained by migration of leukocytes from the blood stream to CSF due to inflammation or, alternatively, blood contamination related to CSF collection. Nevertheless, blood contamination, if present, seems to have not affected the results of the present study as only two samples belonging to MUO and SRMA disease groups (M5 and S1) clustered adjacent to blood plasma samples. Overall, miRNA abundance in CSF samples did not resemble any of the tissues present in the miRNA expression atlas employed according to tissue clustering. Since no reference immune cell populations were present in our clustering analyses, it is difficult to infer a cellular origin of the profiled miRNAs in CSF. However, the close relationship with immune response predicted for the mRNA targets of cfa-miR-142-5p and cfa-miR-92a-3p may further confirm the hypothesis of both MUO and SRMA as immune-mediated diseases.
When we attempted to detect specific miRNAs differentiating MUO from SRMA dogs in the sequencing data, only 3 miRNAs were highlighted (cfa-miR-450a-5p, cfa-miR-450b-5p and cfa-miR-503-5p), but none of them were successfully profiled by qPCR, probably due to their very low abundance as detected by small RNA-seq. Interestingly, all these three miRNAs are reported to be involved in cancer development (66–68).
MiRNA studies in CSF have not previously been performed in dogs with SRMA apart from a few cases included in a population of various CNS diseases (50), and a previous work from our team where we were unsuccessful in profiling CSF by qPCR alone (18). Hence, similar canine studies for comparison are not available. In contrast to MUO, SRMA is characterized by systemic inflammatory changes, and several of the conventional biomarkers of systemic inflammatory conditions, e.g., CRP, fibrinogen, serum amyloid A, and IgA, are also elevated in dogs with SRMA (69–72). Moreover, SRMA is a well-characterized disease with well-described pathological findings (49). In this study, the SRMA group did not meet the same challenges as in the MUO group, since the study population consisted of a more homogeneous group of dogs, all of which fitted the characteristic description for the disease, both clinically and para-clinically. This might explain the higher significance in terms of abundance change detected in dogs suffering from SRMA when compared to the Control group, as opposed to the limited results obtained for dogs suffering from MUO. All dogs had a signalment resembling what has previously been reported regarding age, size, and breed (72, 73), and all showed cervical pain, as well as signs of systemic inflammation, including an increase in CRP above normal reference range (74). Furthermore, CSF analysis showed a pleocytosis dominated by neutrophilic granulocytes in all included dogs, as is classically seen in SRMA (23).
Several clinical challenges arose when recruiting healthy dogs, dogs suffering from SRMA, and especially when recruiting dogs suffering from MUO for CSF collection:
(i) The CSF of the Control group was collected upon euthanasia for ethical reasons, as opposed to most samples in the MUO and SRMA groups where CSF was taken as part of the diagnostic work-up. This might have biased the findings of the present study. Nevertheless, all CSF samples were collected within 30 min of euthanasia, in order to avoid postmortem changes to the CSF composition as suggested in human studies (27). These studies indicated, that if death is rapid, the endothelial lining of the cerebral blood vessels remains intact for hours, meaning that the total nucleated cell count would remain at the same level as at the time of death. No significant rise in protein count was observed within the first 24 h from death either. This, together with the high stability demonstrated by circulating miRNAs, provide robustness to the reliability of results reported in the current study.
(ii) The clinical diagnosis of MUO covers a rather heterogeneous group of non-infectious encephalitides, that, due to similar clinical presentation, can only be differentiated using histopathology. The main subtypes include necrotizing meningoencephalitis (NME), necrotizing leukoencephalitis (NLE) and granulomatous meningoencephalitis (GME), each characterized by specific histopathological findings, which presumably mirror specific underlying immune mechanisms (75). Ideally, future miRNA studies should be performed on dogs with a histopathologically confirmed MUO subtype (NME, NLE, GME). This would possibly offer a study of well-defined subgroups of MUO with characteristic changes in miRNA profiling. However, this approach would introduce several challenges, and a selection bias toward the more severe cases, as histopathology is usually only performed postmortem. Dogs that respond to treatment and survive would therefore not be included in a study where a histopathological diagnosis is an inclusion criterion.
(iii) The biological material used in the present study was CSF, since CSF is known to be the vehicle for intracerebral transport of biologically active substances (24), as well as having a low cell composition in both healthy and pathological states. The low cell composition leads to lesser “transcriptional noise” in sequencing that might be present in bodily fluids with a higher cell and protein content, such as serum. On the other hand, the levels of some relevant miRNAs in CSF may be too low to be confirmed by qPCR as it happened in the present study. Moreover, CSF biomarkers are not ideal for testing in a clinical setting, as CSF collection is an invasive procedure that is only performed under generalized anesthesia. For dogs suffering from SRMA, a tentative diagnosis is usually based on results of CSF analysis when seen in the context of clinical signs, signalment and other paraclinical results, i.e., signs of systemic inflammation (71). A biomarker that can only be identified in CSF is therefore less appealing. For miRNAs to be used as a biomarker in the clinical setting, identification of specific miRNAs in other easily accessible bodily fluids, such as serum or urine, would therefore be preferable. If that was possible, miRNA-profiling could have the potential to replace CSF analysis in dogs where SRMA is the most likely diagnosis.
(iv) Sampling of the CSF in small dogs (as representing the typical case of MUO) can be challenging because only a small amount of CSF can be safely collected. After performing small RNA sequencing, very little amounts of RNA isolated from CSF were left for qPCR verification. Sampling a bigger volume of CSF might be challenging, as the traditional recommendations for CSF collection involve a maximum of 1 mL/5 kg bodyweight (25). For four out of the five dogs included in the MUO group, the weight was below 5 kg. Unfortunately, a standard CSF analysis takes approximately 250 μL, and a minimum of 200 μL are needed for ensuring enough RNA concentration to perform RNA sequencing, leaving very little research material, if any, in dogs smaller than 5 kg.
The RNA isolation method employed in the present study allows purification of cell-free circulating small RNA biotypes including those inside extracellular vesicles (EV) and those attached to protein complexes (Ago2 or HDL). MiRNAs and other small ncRNAs are detected in body fluids at different concentrations depending on the investigated fluid and the experimental group (healthy versus disease). Many of these circulating small RNA biotypes other than miRNAs seem to be derived from larger RNAs as a result of processing events or degradation. In this study, a high percentage of sequencing reads derived from CSF-RNA mapped to Y-RNAs in all 3 experimental groups (with percentages ranging from 44 to 58%); followed by miRNAs (20–36%) with the higher representation in the SRMA group; lncRNAs (4–13%); rRNAs (3–8%) and protein-coding genes (2.6–6%). Of the various RNA biotypes found in circulation, the main focus of this study has been on miRNAs due to their known role in gene-regulation. Nevertheless, a very abundant class of RNAs associated with EV are the Y-RNAs; in fact, in the present study, they were the most abundant class of small non-coding RNAs in CSF. Y-RNAs have been previously found in all biofluids, including CSF (76, 77), which has generated interest in their potential use as disease biomarkers. Intact Y-RNAs have been associated with DNA replication and RNA stability (78). Accumulating evidence suggests that extracellular Y-RNAs present in biofluids may have immune-related functions, and changes in their circulating concentration have been associated with disease (79). Nevertheless, the exact role of fragments derived from Y-RNA genes found in CSF in the present study and their association to MUO and SRMA pathologies, if any, is not clear and further studies are warranted.
In summary, our results should be interpreted as a pilot study setting the bases for larger research projects. We hence recommend for future studies that: (i) they include a larger group of patients, and (ii) they include both serum and CSF for the examination of miRNAs identified as relevant in the SRMA group. If future studies successfully identify miRNAs related to SRMA in serum, there is a potential to develop a less invasive method for diagnosing SRMA. Moreover, miRNA profiling of both the aqueous phase of serum and CSF, as well as their nucleated cellular complement, could provide a more comprehensive view of the origin and biological function of the circulating miRNAs highlighted. In future studies regarding SRMA, we suggest to include two control groups. Control groups ideally include one group of healthy dogs, but also one group of dogs showing symptoms resembling SRMA, i.e., the clinical differential diagnoses, such as dogs suffering from, for instance, immune-mediated polyarthritis or infectious meningitis. This would further qualify miRNAs as potential biomarkers for SRMA in the clinical setting, as this is dependent on their ability to differentiate dogs with SRMA from dogs presenting with similar clinical symptoms, but another underlying pathology. For the MUO group, future studies would potentially benefit from subdivisions based on histopathology of patients (GME, NME and NLE), to identify relevant miRNAs in this disease group. Subdivision could potentially allow a better understanding of the differences between MUO subtypes and help optimize treatment for each.
5. Conclusion
In this pilot study, we have characterized the miRNA composition of the cerebrospinal fluid (CSF) in healthy dogs, dogs with steroid responsive meningitis-arteritis (SRMA) and dogs with meningoencephalitis of unknown origin (MUO). Our analyses based on small RNA-seq and qPCR verified three miRNAs (miR-142-5p, miR-191-5p and miR-92a-3p) present in canine CSF at differential abundance, which would be able to discriminate healthy dogs from dogs with MUO or SRMA. The detected differentially abundant miRNAs were predicted to have biological functions related to the immune response against the ongoing inflammation in dogs suffering from SRMA or MUO. The relatively few miRNAs identified to have discriminatory potential to differentiate healthy from pathological states in CSF might reflect the complex and poorly defined MUO pathology and the diverse subclasses included within this diagnosis, as well as the challenges of accurately profiling circulating RNAs in CSF. Future studies may benefit from a comprehensive differentiation of pathological subtypes within MUO, and the incorporation of additional controls and sequencing information coming from other biofluids, such as serum. Moreover, adding transcriptomes from the cellular fraction present in CSF and serum will help disentangle and differentiate the immune response activated in domestic dogs affected by MUO and SRMA.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: Sequence Read Archive (SRA), accession number PRJNA921866.
Ethics statement
The study was approved by the local ethical committee (permission number 2022–14) at the Department of Veterinary Clinical Sciences, University of Copenhagen. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
EM-S did the bioinformatic analysis on RNA-seq and qPCR data. PH recruited the dogs, performed the clinical work-up, and characterization as well as collection of samples. HG assisted in the clinical work-up and sample collection, secured funding, and helped directing the project. SC did the RNA isolation, qPCR analyses, secured funding, and coordinated the whole project. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by AGRIA Dyreforsikring Fund (Project nr: N2018-0004), Sweden.
Acknowledgments
We thank the Veterinary Technicians at the University Hospital for Companion Animals, University of Copenhagen, Denmark in assisting with sample collection and laboratory analysis, as well as all dog owners for participating in this study.
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1144084/full#supplementary-material
Footnotes
References
1. Lowrie, M, Smith, PM, and Garosi, L. Meningoencephalitis of unknown origin: investigation of prognostic factors and outcome using a standard treatment protocol. Vet Rec. (2013) 172:527. doi: 10.1136/vr.101431
2. Lowrie, M, Penderis, J, Eckersall, PD, McLaughlin, M, Mellor, D, and Anderson, TJ. The role of acute phase proteins in diagnosis and management of steroid-responsive meningitis arteritis in dogs. Vet J. (2009) 182:125–30. doi: 10.1016/j.tvjl.2008.05.001
3. Andersen-Ranberg, E, Berendt, M, and Gredal, H. Biomarkers of non-infectious inflammatory CNS diseases in dogs -where are we now? Part I: meningoencephalitis of unknown origin. Vet J. (2021) 273:105678:273. doi: 10.1016/j.tvjl.2021.105678
4. Coates, JR, and Jeffery, ND. Perspectives on meningoencephalomyelitis of unknown origin. Vet Clin North Am Small Anim Pract. (2014) 44:1157–85. doi: 10.1016/j.cvsm.2014.07.009
5. Granger, N, Smith, PM, and Jeffery, ND. Clinical findings and treatment of non-infectious meningoencephalomyelitis in dogs: a systematic review of 457 published cases from 1962 to 2008. Vet J. (2010) 184:290–7. doi: 10.1016/j.tvjl.2009.03.031
6. Cornelis, I, Volk, HA, van Ham, L, and De Decker, S. Prognostic factors for 1-week survival in dogs diagnosed with meningoencephalitis of unknown aetiology. Vet J. (2016) 214:91–5. doi: 10.1016/j.tvjl.2016.05.008
7. Barnoon, I, Shamir, MH, Aroch, I, Bdolah-Abram, T, Srugo, I, Konstantin, L, et al. Retrospective evaluation of combined mycophenolate mofetil and prednisone treatment for meningoencephalomyelitis of unknown etiology in dogs: 25 cases (2005-2011). J Vet Emerg Crit Care. (2016) 26:116–24. doi: 10.1111/vec.12399
8. Biedermann, E, Tipold, A, and Flegel, T. Relapses in dogs with steroid-responsive meningitis-arteritis. J Small Anim Pract. (2016) 57:91–5. doi: 10.1111/jsap.12418
9. Burgener, I, van Ham, L, Jaggy, A, Vandevelde, M, and Tipold, A. Chemotactic activity and IL-8 levels in the cerebrospinal fluid in canine steroid responsive meningitis-arteriitis. J Neuroimmunol (1998) 89:182–190. doi: 10.1016/s0165-5728(98)00134-9, 182
10. Cizinauskas, S, Jaggy, A, and Tipold, A. Long-term treatment of dogs with steroid-responsive meningitis-arteritis: clinical, laboratory and therapeutic results. J Small Anim Pract. (2000) 41:295–301. doi: 10.1111/j.1748-5827.2000.tb03205.x
11. Günther, C, Steffen, F, Alder, DS, Beatrice, L, Geigy, C, and Beckmann, K. Evaluating the use of cytosine arabinoside for treatment for recurrent canine steroid-responsive meningitis-arteritis. Vet Rec. (2020) 187:e7. doi: 10.1136/vr.105683
12. Burgos, K, Malenica, I, Metpally, R, Courtright, A, Rakela, B, Beach, T, et al. Profiles of extracellular miRNA in cerebrospinal fluid and serum from patients with Alzheimer’s and Parkinson’s diseases correlate with disease status and features of pathology. PLoS One. (2014) 9:e94839. doi: 10.1371/journal.pone.0094839
13. Cortez, MA, Bueso-Ramos, C, Ferdin, J, Lopez-Berestein, G, Sood, AK, and Calin, GA. MicroRNAs in body fluids--the mix of hormones and biomarkers. Nat Rev Clin Oncol. (2011) 8:467–77. doi: 10.1038/nrclinonc.2011.76
14. Singh, RP, Massachi, I, Manickavel, S, Singh, S, Rao, NP, Hasan, S, et al. The role of miRNA in inflammation and autoimmunity. Autoimmun Rev. (2013) 12:1160–5. doi: 10.1016/j.autrev.2013.07.003
15. Alexandrov, PN, Dua, P, Hill, JM, Bhattacharjee, S, Zhao, Y, and Lukiw, WJ. microRNA (miRNA) speciation in Alzheimer’s disease (AD) cerebrospinal fluid (CSF) and extracellular fluid (ECF). Int J Biochem Mol Biol. (2012) 3:365–73.PMC3533883
16. Vansteenkiste, DP, Fenger, JM, Fadda, P, Martin-Vaquero, P, and da Costa, RC. MicroRNA expression in the cerebrospinal fluid of dogs with and without cervical spondylomyelopathy. J Vet Intern Med. (2019) 33:2685–92. doi: 10.1111/jvim.15636
17. Steudemann, C, Bauersachs, S, Weber, K, and Wess, G. Detection and comparison of microRNA expression in the serum of Doberman pinschers with dilated cardiomyopathy and healthy controls. BMC Vet Res. (2013) 9:12. doi: 10.1186/1746-6148-9-12
18. Cirera, S, Andersen-Ranberg, EU, Langkilde, S, Aaquist, M, and Gredal, H. Challenges and standardization of microRNA profiling in serum and cerebrospinal fluid in dogs suffering from non-infectious inflammatory CNS disease. Acta Vet Scand. (2019) 61:57. doi: 10.1186/s13028-019-0492-y
19. Weber, JA, Baxter, DH, Zhang, S, Huang, DY, How Huang, K, Jen Lee, M, et al. The microRNA spectrum in 12 body fluids. Clin Chem. (2010) 56:1733–41. doi: 10.1373/clinchem.2010.147405
20. Correia, CN, Nalpas, NC, McLoughlin, KE, Browne, JA, Gordon, SV, MacHugh, DE, et al. Circulating microRNAs as potential biomarkers of infectious disease. Front Immunol. (2017) 8:118. doi: 10.3389/fimmu.2017.00118
21. Turchinovich, A, Weiz, L, Langheinz, A, and Burwinkel, B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. (2011) 39:7223–33. doi: 10.1093/nar/gkr254
22. Enelund, L, Nielsen, LN, and Cirera, S. Evaluation of microRNA stability in plasma and serum from healthy dogs. MicroRNA. (2017) 6:42–52. doi: 10.2174/2211536606666170113124114
23. Whitney, MS, and Coates, JR. Cerebrospinal fluid analysis in the dog and cat In:. Veterinary Cytology. 1st ed. New Jersey: John Wiley & Sons (2020) 638–54.
24. di Terlizzi, R, and Platt, SR. The function, composition and analysis of cerebrospinal fluid in companion animals: part I –function and composition. Vet J. (2006) 172:422–31. doi: 10.1016/j.tvjl.2005.07.021
25. di Terlizzi, R, and Platt, SR. The function, composition and analysis of cerebrospinal fluid in companion animals: part II -analysis. Vet J. (2009) 180:15–32. doi: 10.1016/j.tvjl.2007.11.024
26. MacNeill, AL, Andre, BG, Zingale, Y, Packer, RA, and McGrath, S. The effects of iatrogenic blood contamination on total nucleated cell counts and protein concentrations in canine cerebrospinal fluid. Vet Clin Pathol. (2018) 47:464–70. doi: 10.1111/vcp.12639
27. Morris, JA, Harrison, LM, and Telford, DR. Postmortem cerebrospinal fluid pleocytosis: a marker of inflammation or postmortem artifact? Int J Pediatr. (2012) 2012:964074. doi: 10.1155/2012/964074
28. Martin, M . Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. (2011) 17:10–2. doi: 10.14806/ej.17.1.200
29. Chen, S, Zhou, Y, Chen, Y, and Gu, J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560
30. Kozomara, A, Birgaoanu, M, and Griffiths-Jones, S. miRBase: from microRNA sequences to function. Nucleic Acids Res. (2019) 47:D155–62. doi: 10.1093/nar/gky1141
31. Langmead, B, Trapnell, C, Pop, M, and Salzberg, SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. (2009) 10:R25. doi: 10.1186/gb-2009-10-3-r25
32. Liao, Y, Smyth, GK, and Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. (2014) 30:923–30. doi: 10.1093/bioinformatics/btt656
33. Cunningham, F, Allen, JE, Allen, J, Alvarez-Jarreta, J, Amode, MR, Armean, IM, et al. Ensembl 2022. Nucleic Acids Res. (2022) 50:D988–95. doi: 10.1093/nar/gkab1049
34. Robinson, MD, McCarthy, DJ, and Smyth, GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. (2010) 26:139–40. doi: 10.1093/bioinformatics/btp616
35. Robinson, MD, and Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. (2010) 11:R25. doi: 10.1186/gb-2010-11-3-r25
36. Benjamini, Y, and Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc: Series B (Methodological). (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
37. Balcells, I, Cirera, S, and Busk, PK. Specific and sensitive quantitative RT-PCR of miRNAs with DNA primers. BMC Biotechnol. (2011) 11:70. doi: 10.1186/1472-6750-11-70
38. Busk, PK . A tool for design of primers for microRNA-specific quantitative RT-qPCR. BMC Bioinformatics. (2014) 15:29. doi: 10.1186/1471-2105-15-29
39. Cirera, S, and Busk, PK. Quantification of miRNAs by a simple and specific qPCR method. Methods Mol Biol. (2014) 1182:73–81. doi: 10.1007/978-1-4939-1062-5_7
40. Andersen, CL, Jensen, JL, and Ørntoft, TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. (2004) 64:5245–50. doi: 10.1158/0008-5472.CAN-04-0496
41. Koenig, EM, Fisher, C, Bernard, H, Wolenski, FS, Gerrein, J, Carsillo, M, et al. The beagle dog MicroRNA tissue atlas: identifying translatable biomarkers of organ toxicity. BMC Genomics. (2016) 17:649. doi: 10.1186/s12864-016-2958-x
42. Penso-Dolfin, L, Swofford, R, Johnson, J, Alföldi, J, Lindblad-Toh, K, Swarbreck, D, et al. An improved microRNA annotation of the canine genome. PLoS One. (2016) 11:e0153453. doi: 10.1371/journal.pone.0153453
43. Fromm, B, Høye, E, Domanska, D, Zhong, X, Aparicio-Puerta, E, Ovchinnikov, V, et al. MirGeneDB 2.1: toward a complete sampling of all major animal phyla. Nucleic Acids Res. (2022) 50:D204–10. doi: 10.1093/nar/gkab1101
44. McInnes, L, Healy, J, Saul, N, and Großberger, L. UMAP: uniform manifold approximation and projection. J Open Source Soft (2018) 3:861. doi: 10.21105/joss.00861
45. Agarwal, V, Bell, GW, Nam, JW, and Bartel, DP. Predicting effective microRNA target sites in mammalian mRNAs. elife. (2015) 4:e05005. doi: 10.7554/eLife.05005
46. Friedman, RC, Farh, KKH, Burge, CB, and Bartel, DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. (2009) 19:92–105. doi: 10.1101/gr.082701.108
47. Bindea, G, Mlecnik, B, Hackl, H, Charoentong, P, Tosolini, M, Kirilovsky, A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101
48. Shannon, P, Markiel, A, Ozier, O, Baliga, NS, Wang, JT, Ramage, D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
49. Gaitero, L, Russell, SJ, Monteith, G, and LaMarre, J. Expression of microRNAs miR-21 and miR-181c in cerebrospinal fluid and serum in canine meningoencephalomyelitis of unknown origin. Vet J. (2016) 216:122–4. doi: 10.1016/j.tvjl.2016.07.014
50. Marioni-Henry, K, Zaho, D, Amengual-Batle, P, Rzechorzek, NM, and Clinton, M. Expression of microRNAs in cerebrospinal fluid of dogs with central nervous system disease. Acta Vet Scand. (2018) 60:80. doi: 10.1186/s13028-018-0434-0
51. Shrestha, A, Mukhametshina, RT, Taghizadeh, S, Vásquez-Pacheco, E, Cabrera-Fuentes, H, Rizvanov, A, et al. MicroRNA-142 is a multifaceted regulator in organogenesis, homeostasis, and disease. Dev Dyn. (2017) 246:285–90. doi: 10.1002/dvdy.24477
52. Cogswell, JP, Ward, J, Taylor, IA, Waters, M, Shi, Y, Cannon, B, et al. Identification of miRNA changes in Alzheimer’s disease brain and CSF yields putative biomarkers and insights into disease pathways. J Alzheimers Dis. (2008) 14:27–41. doi: 10.3233/JAD-2008-14103
53. Narasimhan, M, Patel, D, Vedpathak, D, Rathinam, M, Henderson, G, and Mahimainathan, L. Identification of novel microRNAs in post-transcriptional control of Nrf2 expression and redox homeostasis in neuronal, SH-SY5Y cells. PLoS One. (2012) 7:e51111. doi: 10.1371/journal.pone.0051111
54. Talebi, F, Ghorbani, S, Chan, WF, Boghozian, R, Masoumi, F, Ghasemi, S, et al. MicroRNA-142 regulates inflammation and T cell differentiation in an animal model of multiple sclerosis. J Neuroinflammation. (2017) 14:55. doi: 10.1186/S12974-017-0832-7
55. Prieto-Fernández, E, Lopez-Lopez, E, Martin-Guerrero, I, Bárcena, L, Gonzalez-Lopez, M, Aransay, AM, et al. Variability in cerebrospinal fluid microRNAs through life. Mol Neurobiol. (2020) 57:4134–42. doi: 10.1007/s12035-020-02011-3
56. Nagpal, N, and Kulshreshtha, R. miR-191: an emerging player in disease biology. Front Genet. (2014) 5:99. doi: 10.3389/fgene.2014.00099
57. Petrova, TA, Kondratyev, SA, Kostareva, AA, Rutkovskiy, RV, Savvina, IA, and Kondratyeva, EA. miR-21, miR-93, miR-191, miR-let-7b, and miR-499 expression level in plasma and cerebrospinal fluid in patients with prolonged disorders of consciousness. Neurol Int. (2022) 15:40–54. doi: 10.3390/neurolint15010004
58. Wang, L, Shui, X, Zhang, M, Mei, Y, Xia, Y, Lan, G, et al. MiR-191-5p attenuates tau phosphorylation, Aβ generation, and neuronal cell death by regulating death-associated protein kinase 1. ACS Chem Neurosci. (2022) 13:3554–66. doi: 10.1021/acschemneuro.2C00423
59. Wang, L, Cui, M, Qu, F, Cheng, D, Yu, J, Tang, Z, et al. MiR-92a-3p promotes the malignant progression of hepatocellular carcinoma by mediating the PI3K/AKT/mTOR signaling pathway. Curr Pharm Des. (2021) 27:3244–50. doi: 10.2174/1381612827666210612054156
60. Peña-Bautista, C, Tarazona-Sánchez, A, Braza-Boils, A, Balaguer, A, Ferré-González, L, Cañada-Martínez, AJ, et al. Plasma microRNAs as potential biomarkers in early Alzheimer disease expression. Sci Rep. (2022) 12:15589. doi: 10.1038/S41598-022-19862-6
61. Sanjabi, S, Oh, SA, and Li, MO. Regulation of the immune response by TGF-β: from conception to autoimmunity and infection. Cold Spring Harb Perspect Biol. (2017) 9:a022236. doi: 10.1101/cshperspect.a022236
62. Arthur, JSC, and Ley, SC. Mitogen-activated protein kinases in innate immunity. Nat Rev Immunol. (2013) 13:679–92. doi: 10.1038/nri3495
63. Becker, T, Loch, G, Beyer, M, Zinke, I, Aschenbrenner, AC, Carrera, P, et al. FOXO-dependent regulation of innate immune homeostasis. Nature. (2010) 463:369–73. doi: 10.1038/nature08698
64. Wehbi, VL, and Taskén, K. Molecular mechanisms for cAMP-mediated immunoregulation in T cells -role of anchored protein kinase a signaling units. Front Immunol. (2016) 7:222. doi: 10.3389/fimmu.2016.00222
65. Hong, L, Li, X, Zhou, D, Geng, J, and Chen, L. Role of hippo signaling in regulating immunity. Cell Mol Immunol. (2018) 15:1003–9. doi: 10.1038/s41423-018-0007-1
66. Ye, YP, Wu, P, Gu, CC, Deng, D-l, Jiao, HL, Li, TT, et al. miR-450b-5p induced by oncogenic KRAS is required for colorectal cancer progression. Oncotarget. (2016) 7:61312–24. doi: 10.18632/oncotarget.11016
67. Jiang, SP, and Li, ZR. MiR-503-5p regulates cell epithelial-to-mesenchymal transition, metastasis and prognosis of hepatocellular carcinoma through inhibiting WEE1. Eur Rev Med Pharmacol Sci. (2019) 23:2028–37. doi: 10.26355/eurrev_201903_17242
68. Zhao, YJ, Zhang, J, Wang, YC, Wang, L, and He, XY. MiR-450a-5p inhibits gastric cancer cell proliferation, migration, and invasion and promotes apoptosis via targeting CREB1 and inhibiting AKT/GSK-3β signaling pathway. Front Oncol. (2021) 11:633366. doi: 10.3389/fonc.2021.633366
69. de la Fuente, C, Monreal, L, Cerón, J, Pastor, J, Viu, J, and Añor, S. Fibrinolytic activity in cerebrospinal fluid of dogs with different neurological disorders. J Vet Intern Med. (2012) 26:1365–73. doi: 10.1111/j.1939-1676.2012.00991.x
70. Bathen-Noethen, A, Carlson, R, Menzel, D, Mischke, R, and Tipold, A. Concentrations of acute-phase proteins in dogs with steroid responsive meningitis-arteritis. J Vet Intern Med. (2008) 22:1149–56. doi: 10.1111/j.1939-1676.2008.0164.x
71. Andersen-Ranberg, E, Berendt, M, and Gredal, H. Biomarkers of non-infectious inflammatory CNS diseases in dogs: where are we now? Part 2 -steroid responsive meningitis-arteritis. Vet J. (2021) 273:105692. doi: 10.1016/j.tvjl.2021.105692
72. Lowrie, M, Penderis, J, McLaughlin, M, Eckersall, PD, and Anderson, TJ. Steroid responsive meningitis-arteritis: a prospective study of potential disease markers, prednisolone treatment, and long-term outcome in 20 dogs (2006-2008). J Vet Intern Med. (2009) 23:862–70. doi: 10.1111/j.1939-1676.2009.0337.x
73. Lau, J, Nettifee, JA, Early, PJ, Mariani, CL, Olby, NJ, and Muñana, KR. Clinical characteristics, breed differences, and quality of life in north American dogs with acute steroid-responsive meningitis-arteritis. J Vet Intern Med. (2019) 33:1719–27. doi: 10.1111/jvim.15543
74. Tipold, A, and Schatzberg, SJ. An update on steroid responsive meningitis-arteritis. J Small Anim Pract. (2010) 51:150–4. doi: 10.1111/j.1748-5827.2009.00848.x
75. Cornelis, I, van Ham, L, Gielen, I, de Decker, S, and Bhatti, SFM. Clinical presentation, diagnostic findings, prognostic factors, treatment and outcome in dogs with meningoencephalomyelitis of unknown origin: a review. Vet J. (2019) 244:37–44. doi: 10.1016/j.tvjl.2018.12.007
76. Godoy, PM, Bhakta, NR, Barczak, AJ, Cakmak, H, Fisher, S, MacKenzie, TC, et al. Large differences in small RNA composition between human biofluids. Cell Rep. (2018) 25:1346–58. doi: 10.1016/j.celrep.2018.10.014
77. Lee, KY, Seo, Y, Im, JH, Rhim, J, Baek, W, Kim, S, et al. Molecular signature of extracellular vesicular small non-coding RNAs derived from cerebrospinal fluid of leptomeningeal metastasis patients: functional implication of miR-21 and other small RNAs in cancer malignancy. Cancers. (2021) 13:209. doi: 10.3390/cancers13020209
78. Kowalski, MP, and Krude, T. Functional roles of non-coding Y RNAs. Int J Biochem Cell Biol. (2015) 66:20–9. doi: 10.1016/j.biocel.2015.07.003
Keywords: next generation sequencing, microRNA, dog, cerebrospinal fluid, steroid responsive meningitis-arteritis, meningoencephalitis of unknown origin
Citation: Mármol-Sánchez E, Heidemann PL, Gredal H and Cirera S (2023) MicroRNA profiling of cerebrospinal fluid from dogs with steroid responsive meningitis-arteritis and meningoencephalitis of unknown origin. Front. Vet. Sci. 10:1144084. doi: 10.3389/fvets.2023.1144084
Edited by:
Holger Andreas Volk, University of Veterinary Medicine Hannover, GermanyReviewed by:
Rodrigo Gutierrez-Quintana, University of Glasgow, United KingdomJasmin Nicole Nessler, University of Veterinary Medicine Hannover, Germany
Copyright © 2023 Mármol-Sánchez, Heidemann, Gredal and Cirera. 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: Hanne Gredal, hbg@sund.ku.dk; Susanna Cirera, scs@sund.ku.dk
†These authors have contributed equally to this work and share first authorship