Shruti Singh Kakan1
Srikanth R. Janga2
Benjamin Cooperman1
David W. Craig3Maria C. Edman2
Curtis T. Okamoto1
Sarah F. Hamm-Alvarez1,2*- 1Department of Pharmacology and Pharmaceutical Sciences, School of Pharmacy, University of Southern California, Los Angeles, CA, United States
- 2Department of Ophthalmology, Keck School of Medicine, Roski Eye Institute, University of Southern California, Los Angeles, CA, United States
- 3Department of Translational Genomics, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
Sjögren's Syndrome (SS) is an autoimmune disease characterized by lymphocytic infiltration and loss of function of moisture-producing exocrine glands as well as systemic inflammation. SS diagnosis is cumbersome, subjective and complicated by manifestation of symptoms that overlap with those of other rheumatic and ocular diseases. Definitive diagnosis averages 4–5 years and this delay may lead to irreversible tissue damage. Thus, there is an urgent need for diagnostic biomarkers for earlier detection of SS. Extracellular vesicles called exosomes carry functional small non-coding RNAs which play a critical role in maintaining cellular homeostasis via transcriptional and translational regulation of mRNA. Alterations in levels of specific exosomal miRNAs may be predictive of disease status. Here, we have assessed serum exosomal RNA using next generation sequencing in a discovery cohort of the NOD mouse, a model of early-intermediate SS, to identify dysregulated miRNAs that may be indicative of SS. We found five miRNAs upregulated in serum exosomes of NOD mice with an adjusted p < 0.05—miRNA-127-3p, miRNA-409-3p, miRNA-410-3p, miRNA-541-5p, and miRNA-540-5p. miRNAs 127-3p and 541-5p were also statistically significantly upregulated in a validation cohort of NOD mice. Pathway analysis and existing literature indicates that differential expression of these miRNAs may dysregulate pathways involved in inflammation. Future studies will apply these findings in a human cohort to understand how they are correlated with manifestations of SS as well as understanding their functional role in systemic autoimmunity specific to SS.
Introduction
Sjögren's Syndrome (SS) is a chronic and systemic autoimmune disease marked by lymphocytic infiltration and loss of function of the body's moisture producing exocrine glands (e.g., lacrimal and salivary glands) as its defining manifestation. It is the second most common rheumatic autoimmune disease, affecting about 0.5–1% of the general population (1, 2). The progressive inflammation of lacrimal and salivary exocrine glands is associated with their loss of function, leading to debilitating dry eye and dry mouth, respectively (3). SS is associated with increased inflammation of internal organs including brain, lung and liver (4) as well as a 44-fold increased risk of developing B-cell lymphoma (5, 6). SS can occur in the absence of another autoimmune disease (primary SS) or concurrently with another autoimmune disease such as rheumatoid arthritis or systemic lupus erythematosus (secondary SS).
Diagnosis of SS relies on the weighted score obtained from a series of criteria established in 2016 by the American College of Rheumatology (ACR) in collaboration with the European League Against Rheumatism (EULAR) (7). These criteria include: (1) labial salivary gland biopsy showing focal lymphocytic sialadenitis with a focus score ≥1; (2) anti-SSA (Ro) positivity (serum autoantibody); (3) ocular surface staining score ≥ 5 (or van Bijsterveld score ≥ 4) on at least one eye; (4) a Schirmer's value (tear flow) ≤ 5 mm/5 min on at least one eye; and (5) an unstimulated whole saliva flow rate ≤ 0.1 ml/min. Although, commonly used for inclusion in clinical trials, these criteria are not always practically applicable for clinical diagnosis. In particular, the labial salivary gland biopsy is painful, impractical and error-prone (8). Furthermore, these criteria have been developed primarily for patients with primary SS, and are not extensively validated in the far greater numbers of patients suffering from secondary SS. Therefore, many patients experience a substantial delay in diagnosis while some are never formally diagnosed. This delay in diagnosis may also delay treatment with anti-inflammatory agents to the point when irreversible damage to exocrine glands and other internal organs may already have occurred (9). Hence, there is an urgent need for an early, sensitive and non-invasive diagnostic test for SS.
In this study we utilized the Non Obese Diabetic mouse model of SS. The NOD/Shi strain originated from inbreeding of the Cataract Shionogi (CTS) strain, based on elevated fasting blood glucose level in cataract-free mice for the development of a model for insulin dependent diabetes mellitus. It was later also shown to develop features of exocrinopathy consistent with SS (10). Despite limitations to any animal model, this model exhibits several features of SS in humans including reduced tear and salivary secretion (11, 12), alterations in tear and salivary composition (13–15), lymphocytic infiltration of lacrimal and salivary glands (dacryoadenitis and sialadenitis, respectively) (11, 12), and the presence of many of the serum autoantibodies that are often present in human patients, such as autoantibodies to Ro/SSA and La/SSB (16), the M3 muscarinic acetylcholine receptor (17, 18), salivary gland protein 1 (19), carbonic anhydrase 6 (19) and parotid secretory protein (PSP) (19). Finally, the lacrimal glands (LG) of these mice show characteristic changes in specific proteins involved in the secretory process of exocrine glands typical of SS patients (20, 21). We have previously demonstrated that tear biomarkers identified in this murine model are also identified in SS patients (13–15).
Although, SS is more prevalent in women than in men at a 9:1 ratio (3), only male NOD mice were used in this study, because the male mice have been extensively characterized to exhibit the features of autoimmune dacryoadenitis and systemic disease prior to the development of diabetes. Females of this strain instead develop autoimmune sialadenitis concurrent with diabetes (22), complicating interpretation of any results. Thus, use of male NOD mice allows us to avoid confounding effects associated with the concurrent development of diabetes. Our mice were chosen at the age just after lymphocytic infiltration of the LG is typically established (14-weeks), representing an early-intermediate stage disease model of autoimmune dacryoadenitis in SS. As the disease development in this strain is polygenetic, and many of the diabetes resistant sub-strains that have been developed as controls for studying diabetes development still develop autoimmune exocrinopathy (11, 23, 24) there has been a lack of closely related healthy control strains for studies of SS disease development and treatment. With this said, the Balb/c strain has been the most commonly used for studies of SS exocrinopathy by multiple groups beyond our own (25–28). Therefore, we considered it the most prudent choice for use as a control strain.
MicroRNAs (miRNA) are evolutionarily conserved short non-coding RNA that function in gene silencing and post-transcriptional gene regulation (29), regulating nearly 60% of messenger RNA (mRNA) (30). A single miRNA can target several mRNA and any given mRNA may be targeted by more than one miRNA. Cooperativity between a group of dysregulated miRNAs targeting one or more mRNAs of a given signaling pathway or cellular process may substantially upregulate or downregulate that pathway and lead to development and progression of disease. Indeed, miRNA dysregulation is associated with cancer (31), obesity (32), heart disease (33), kidney disease (34), and diseases of the nervous system (35). Their diagnostic potential has been explored with high fidelity in various cancers (36), as well as neurodegenerative (37), autoimmune (38), and metabolic diseases (39). Thus, assessing the level of expression of a panel of functional mature miRNAs can be diagnostic of a given disease. Compared to proteins, miRNAs have lower inter-individual variation and less sequence heterogeneity, allowing for high specificity as biomarkers (40). Further, as miRNAs are highly evolutionarily conserved in mammals, results from an animal model are typically readily applicable to human subjects (41).
miRNAs circulate in a stable, cell-free form in all biofluids and are particularly enriched in serum (42). Relative to other biofluids such as saliva, urine, and cerebrospinal fluid (CSF), serum has a higher concentration of miRNA (42, 43). Most extracellular miRNAs as well as other non-coding RNAs (ncRNA) can be found in exosomes which are nano-sized extracellular vesicles, generated as intraluminal vesicles in multivesicular bodies (MVBs). Exosomes are actively shed from nearly every cell type and engage in intercellular signaling. Literature suggests that aberrant exosome-based intercellular communication plays a role in infectious and inflammatory diseases and various cancers. Exosomes are an extremely reliable source of miRNA as their contents are resistant to degradation by nucleases such as RNase (36, 44–47). Moreover, exosomal miRNAs are relatively stable following storage at −80°C, have a very low inter-individual variance, and also exhibit a low intra-individual variance over time (48).
The goal of this study was to understand how changes in regulatory small RNAs in serum, particularly miRNAs, might have utility as a source of diagnostic biomarkers for SS. Here we identify a subset of dysregulated miRNAs that may be specific to SS. We identified 5 exosomal miRNAs that showed significant changes in concert with establishment of SS-like symptoms in NOD mice, representing putative biomarkers for early disease diagnosis. Additionally, we have also assessed differential expression of other small non-coding RNA such as piRNA, that protect the genome by silencing transposons.
Materials and Methods
Mice
Age-matched male NOD/ShiltJ (Stock No. 001976) and BALB/cJ (Stock No. 000651) mice were purchased from Jackson Laboratories (Bar Harbor, ME) and housed with a 12 h light, 12 h dark cycle with ad libitum access to food and water until 14-weeks of age, when SS-like ocular symptoms are established in the NOD strain. All procedures performed on the mice were in accordance with protocols approved by the University of Southern California's Institutional Animal Care and Use Committee (IACUC) and the Guide for Care and Use of Laboratory Animals 8th edition (49).
LG Histology and Quantitative Analysis of Lymphocytic Infiltration
Lymphocytic infiltration in mouse lacrimal glands, indicative of autoimmune dacryoadenitis, was confirmed and quantified with hematoxylin & eosin staining as described (50). Briefly, lacrimal glands from NOD and BALB/c mice were fixed in 10% NBF (Richard-Allan Scientific, Kalamazoo, MI), fixed in paraffin, then cut into 5 μ horizontal sections and stained with hematoxylin and eosin. Sections were imaged using an Aperio Digital ScanScope (Leica Biosystems Inc., Buffalo Grove, IL) using the 40x objective lens. The percentage of lymphocytic infiltration in the tissue was determined by calculating the area of infiltrates manually using ImageJ (National Institutes of Health, http://imagej.nih.gov/ij). Data were analyzed by GraphPad Prism using one-way non-parametric ANOVA (Kruskal-Wallis).
Isolation of Serum Exosomes
Mice were anesthetized by intraperitoneal injection with ketamine/xylazine (60–70 mg + 5–10 mg/kg, respectively), and blood was collected by cardiac puncture using a 1 mL syringe (BD Biosciences, San Jose, CA) into MiniCollect 0.8 mL gold cap Z Serum Separator tubes (Greiner Bio-One, Kremsmünster, Austria). Thereafter the mice were euthanized by cervical dislocation. Blood was allowed to clot for 20 min at room temperature followed by centrifugation at 4°C, 2,000 × g for 15 min. Serum was collected and spun at 2,000 × g for 20 min at 4°C to pellet cellular debris. The supernatant was collected and spun at 12,000 × g for 45 min to remove microvesicles. Approximately 2000 μL of pooled supernatant from 5 mice was concentrated using 10 kDa Millipore Amicon Ultra concentrators (Burlington MA) to 500 μL and then loaded on an equilibrated iZON qEV original size exclusion column (Christchurch, New Zealand). Fractions 7–9 containing 1.5 mL of exosomes were collected and concentrated by 10-fold using 100 kDa concentrators (MilliporeSigma, Burlington, MA). Alternatively, exosomes were enriched from the supernatant obtained by centrifugation of serum as above and resolved by differential ultracentrifugation as previously described (51) with some modifications using a Beckman Coulter Optima LE-80k with a Beckman Coulter Type 50.2 Fixed Angle Rotor. Briefly, after the 12,000 × g spin of the clarified serum, the supernatant was centrifuged at 110,000 × g for 120 min. The pellet was resuspended in 2 mL of PBS containing 0.25 mM Trehalose (PBST) and centrifuged at 110, 000 × g for 70 min. The exosome pellet was resuspended in 200 μL of PBST. Purified exosomes from each protocol were used directly for RNA isolation, Western blotting or flash frozen and stored at −80°C for later analysis. A total of 5 groups per strain and 5 mice per group were utilized as biological replicates for the discovery as well as validation cohort.
Total RNA Isolation
RNA was isolated using the miRNeasy Serum/Plasma Mini Kit (Qiagen, Hilden, Germany). The manufacturer's protocols were followed as written, except for the final collection step which was performed sequentially in two steps. First, 25 μL of nuclease-free water was added to the spin column for 10 min before elution of sample. Then, this step was repeated using 15 μL of nuclease-free water to increase the recovery yield. Combination of both eluates yielded around 30 μL of total RNA collected per pooled exosome sample from five mice. The amount and quality of RNA was analyzed using a Nanodrop to assess initial concentration, and TapeStation (Agilent) to assess sample quality utilizing RNA integrity number (RIN).
Transmission Electron Microscopy
Exosomes stored at −80°C were thawed and fixed on 150 mesh copper carbon formvar grids (Electron Microscopy Sciences, Hatfield, PA). With high precision negative forceps (Electron Microscopy Sciences, Hatfield, PA), 10 μL of exosome samples were incubated with grids for 5 min. Excess liquid was absorbed using filter paper. The grid was incubated with 1% aqueous uranyl acetate (Electron Microscopy Sciences, Hatfield, PA) for 5 min. After rinsing with 10 μL ultrapure water, the grid was air dried for 30 min before storage or immediate viewing in a JEM1400 transmission electron microscope operating at 100 keV.
Western Blotting
Equal volumes of exosome samples were heated for 5 min at 95°C under reducing conditions and resolved over 8–16% Novex WedgeWell Tris-Glycine Polyacrylamide Gels (ThermoFisher, Waltham, MA) for 90 min at 125 volts, under constant voltage. Proteins in gels were transferred to nitrocellulose membrane using an iBLOT 2 device and Invitrogen iBLOT 2 NC stacks (ThermoFisher, Waltham, MA). Membranes were rinsed in Phosphate Buffered Saline (PBS) and blocked in Rockland Blocking Buffer for Fluorescent Western Blotting (Pottstown, PA) for 1 h at room temperature. Membranes were incubated with rabbit primary polyclonal antibodies to TSG101 [Abcam—EPR7130(B), 1:250 dilution], and primary monoclonal antibody to Cathepsin L (Abcam—EPR8011, 1:500 dilution) overnight at 4°C. After six 5 min washes in 1x PBS, membranes were incubated in goat-anti rabbit IR800 secondary antibody for 1 h at RT and rinsed again with 1x PBS, 6 times for 5 min each before imaging on a LI-COR Odyssey Fluorescent Imager. Images were analyzed using ImageStudio v5.2.5.
Particle Size Analysis
Size and concentration of exosomes was measured by Nanoparticle Tracking Analysis (NTA) using a ZetaView (Particle Metrix, Meerbusch, Germany). Some samples were also shipped to Alpha Nano Tech LLC (Chapel Hill, NC) for analysis by a ZetaView S/N 17-332 running the software ZetaView 8.04.02. After calibration with 100 nm standards (Applied Microspheres, The Netherlands), samples were diluted in varying amounts of PBS to reach the optimal concentration for analysis, then injected into the ZetaView cell for measurement. Eleven cell positions were sampled for two cycles each, with outliers automatically removed by the software. Measurements were taken at 22°C, using a sensitivity of 75, a frame rate of 30, and a shutter speed of 100. These measurements were analyzed using a minimum brightness of 20, a maximum size of 500 pixels, and a minimum size of 10 pixels. As it is the best determinant of particle size, the mode was selected as the main sizing parameter (52). Total particle count was calculated to account for varying resuspension volumes.
Particle size was also analyzed by Dynamic Light Scattering (DLS) using a Wyatt Dyna-Pro Plate reader II (Wyatt Technologies, Santa Barbara, CA). Briefly, 60 μL of exosome samples were run in triplicates at 25°C in a 384-well clear bottom plate (Greiner Bio One, Monroe, NC). The hydrodynamic radius of isolated exosomes was measured and presented as a normalized diameter. Data was analyzed using Dynamics V7 software (Wyatt, Santa Barbara, CA).
Small RNA Deep Sequencing
Library preparation and sequencing on exosome fractions were performed by GeneWiz (South Plainfield, NJ). Total RNA containing the small RNA fraction was converted into cDNA using the Illumina TruSeq Small RNA library prep kit according to the manufacturer's instructions. Briefly, 3' adapter “RA3” and 5' adapter “RA5,” were ligated to total RNA which was then reverse transcribed. Adapter ligated cDNA library was enriched by PCR using primers that selectively anneal to the adapter sequence and then purified by gel electrophoresis. Quality of the cDNA library was assessed using a DNA chip on bioanalyzer. Barcodes were added to each sample and all 10 samples were sequenced on a single lane of an Illumina HiSeq system set to a 2 × 150 bp configuration. The output generated a total of ~414 million reads which were then demultiplexed with the added barcode separating the files according to the samples into FASTQ files.
Bioinformatics
The raw FASTQ files obtained from Genewiz were assessed for their quality using FastQC v0.11.9. Adapter trimming was performed using Cutadapt v2.8 (https://github.com/marcelm/cutadapt). High quality reads of minimum length 15 nucleotides (nt) were mapped to whole genome (mm10 assembly GRCm38) using Bowtie v1.2.3 (53) and annotated using featureCounts v2.0.0 (54), using GENCODE (Release M24, GRCm38.p6) comprehensive gene annotation GTF file (PRI) to obtain distribution of reads over genome and raw counts for various non-coding RNA such as pre-miRNA, scRNA, scaRNA, snRNA, tRNA, rRNA, snoRNA, and lncRNA. Reads were then mapped to small RNA transcriptomes (miRbase v22, piRdbv2.0). The output files in the SAM file format were sorted to mapped reads that had an alignment CIGAR string of 18M or higher. Sam2counts, a python program (https://github.com/vsbuffalo/sam2counts), was used to acquire counts of reads aligned to transcriptomes (55).
Raw reads were also aligned to the piRNA and miRNA transcriptomes using an in-house aligner “miRGrep” (https://github.com/singhkakan/miRGrep) that applies brute force to count the number of reads containing a given miRNA or piRNA sequence. As miRNA are 19–26 bp long and piRNA are 24–30 bp long, the entirety of their sequence is read during sequencing. As a result, the reads (75 to 150 bp) are longer than the miRNA or piRNA of interest and contain their complete sequence. miRGrep yielded a final count table which was assembled in RStudio using the dplyr package for further processing. After the miRNA or piRNA counts table was generated, differential gene expression analysis was conducted using three statistical R packages DESEq2 (56), EdgeR (57), and LimmaVoom (58) in RStudio. Statistical significance was determined by adjusted p < 0.05 by DESEq2 or Limma or False Discovery Rate (FDR) <0.1 by EdgeR. We have included miRNAs considered significant by at least 1 statistical package for downstream analyses. The experimental procedures and analysis pipeline are detailed in Supplementary Figure 1.
miRNA Validation Assays
In a separate validation cohort of 5 groups with five mice per group, we isolated serum exosomal RNA. The differential expression of miRNAs of interest was validated by qRT-PCR using individual Taqman Advanced miRNA Assays (Applied Biosystems). Briefly, poly-A tailing and adapter ligation was performed on 2 μL of total RNA isolated from serum exosomes of the validation cohort using the Taqman Advanced cDNA synthesis kit. Following this cDNA synthesis, miRNA amplification was conducted using the same kit. The amplified cDNA was diluted 1:10 and set up in triplicate qRT-PCRs with 1 μL of specific Taqman Advanced miRNA primer and run on a Quant-Studio Flex 6 (Applied Biosystems, Foster City, CA), using the assay's recommended cycling conditions. Results were analyzed by the ΔΔCt method with BALB/c serum exosomal small RNA as reference and miR-16-5p as housekeeping miRNA. 16-5p is identified in the literature as a suitable housekeeping miRNA (59–62) and was unchanged between serum exosomes in BALB/c and NOD mice in the sequencing data obtained in this study (p = 0.995, DESeq2).
Pathway and Functional Enrichment Analysis
Pathway analyses were conducted using miTALOS v2.0 (63) using StarBase2, a database of experimentally validated miRNA targets. Pathway data were extracted from KEGG, Reactome, and WikiPathways by miTALOS. Pathways with a corrected p < 0.05 and Enrichment score > 1 are expected to contain over-represented miRNA targets. Functional enrichment analysis was done using the custom heatmap calculator of miRPathDB v2.0 (64) using GeneOntology (Biological Enrichment) with at least two miRNA per pathway and two pathways per miRNA. The settings were chosen such that signaling pathways in which at least 2 of our miRNAs of interest have an mRNA target, would be identified. Additionally, the program was directed to identify at least two pathways targeted by each miRNA. With these constraints we may identify pathways that have a high probability of being dysregulated with the aberrant expression of the miRNA of interest. Pathway analysis was also conducted using Ingenuity Pathway Analysis (IPA) to visualize the interaction of miRNA “hits” with their targets in the signaling pathways relevant to autoimmunity, which were also identified by miTALOS and miRPathDB.
Results
Characterization of Exosomes From Mouse Serum
We chose exosomes as the principal source of serum miRNAs based on findings that miRNAs are concentrated in these organelles in extracellular biofluids. Exosomes were isolated by differential ultracentrifugation (UC) for discovery experiments and by size exclusion chromatography (SEC) for validation experiments. Transmission electron microscopy (TEM) images showed ~100 nm sized vesicles with the characteristic cup shaped morphology typical of exosomes (65) for both UC (Figure 1A) and SEC exosomes (Figure 1C). UC exosomes had a median diameter of 136 nm from male NOD mouse serum and 131 nm for BALB/c mouse serum, while SEC exosomes had a median diameter of 122 nm for NOD and BALB/c mouse serum (Figures 1B,D). There was no significant difference in the sizes of the exosomes isolated by the two methods, although exosomes isolated by UC were slightly larger in accord with previous reports (66). SEC has been reported to give a better yield than UC (67), consistent with our findings. There were no significant differences in the size or concentration of exosomes between the two strains of mice. Western blotting also showed the presence of the universal exosome membrane protein marker, TSG101, in exosomes isolated from both strains (Figure 1E). Various cathepsins (S, D, K and L) have been found in exosomes derived from plasma (68), macrophages (69), and microglia (70). Interestingly, we found enrichment of cathepsin L in exosomes isolated from both strains (Figure 1E). DLS found that SEC exosomes ranged from 60 to 130 nm with a mean diameter of 122 nm (Figure 1F). There was no strain specific difference in size, concentration or marker expression between exosomes isolated from the two strains. For study design, UC exosomes were used for miRNA “hit” identification and SEC exosomes were used for “hit” validation.
Figure 1. Characterization of mouse serum-derived exosomes by differential Ultracentrifugation (UC) and Sice Exclusion Chromatography (SEC). (A) Transmission Electron Microscopy (TEM) images of exosomes isolated from NOD (left) and BALB/c (right) mouse serum using differential UC. The cup and saucer shaped morphology typical of exosomes is visible. (B) NTA showed median diameters of 136 and 131 nm for serum exosomes from NOD and BALB/c mice, respectively. The graph is representative of four separate experiments per mouse strain. No significant differences in exosome median diameter were observed between the two strains. (C) TEM images of SEC exosomes isolated from NOD (left) and BALB/c (right) mouse serum. (D) NTA of SEC exosomes. (E) Western blotting of NOD and BALB/c serum exosomes isolated by SEC show enrichment of TSG101 and Cathepsin L. (F) Exosome particle size of SEC exosomes analyzed by DLS showed a median diameter of 122 nm from both strains. The graph is representative of three separate experiments.
Mouse Serum Exosomes Contain Several Small RNA Biotypes
Since most extracellular ncRNA are associated with extracellular vesicles (71), UC serum exosomes were used as a source for small RNAs from NOD and BALB/c mice. Establishment of autoimmune dacryoadenitis, the most notable characteristic of SS in these mice, was confirmed in each NOD cohort relative to BALB/c by H & E staining of lacrimal gland sections and quantitation of lymphocytic infiltrates (Supplementary Figure 2). While we did not observe any lymphocytes infiltrating LG in the BALB/c, NOD mice had infiltration between 8 and 15% (p < 0.0001, one-way ANOVA). Using next-generation sequencing (NGS), we profiled small RNA in serum exosomes from each mouse strain to identify both novel and differentially expressed miRNA. A total of 417 million raw reads were generated which were mapped to the mouse genome and various small ncRNA transcriptomes. Of the reads that mapped to the mouse genome, roughly two-thirds of the reads mapped to intergenic or intronic regions while a third of the reads mapped to exons. Less than 10% of the reads mapped to transcription start or end sites (Figure 2A). Read distribution was fairly uniform across samples with no strain specific difference. More than 50% of mapped reads were comprised of miRNA (both mature and precursor) and piRNA. Nearly 27% of reads mapped to lncRNA, and 20% to rRNA, whereas <1% mapped to small nucleolar RNA (snoRNA) (Figures 2B,C) in serum exosomes of both NODs and BALB/c mice. Reads that mapped to miRNA, piRNA, lncRNA, snRNA, snoRNA, scaRNA and rRNA are provided in Supplementary Table 1. We did not observe any strain-specific differences in proportion of RNA sub-types within the small ncRNA libraries of NOD and BALB/c serum exosomes (p = 0.925, ordinary two-way ANOVA).
Figure 2. Graphical presentations of the distribution of reads from NGS sequencing of RNAs. (A) Shows the regions of the mouse genome that reads of all samples were aligned to. The mapping percentage is >100 as multi-mapping reads were also counted. No significant sample to sample variation was observed. (B,C) Show the distribution of reads aligning to various ncRNAs in NOD & BALB/c serum UC exosomes. Results were obtained from n = 5 UC serum exosome samples per mouse strain, with each sample comprising pooled serum exosomes from five mice. TSS, Transcription Start Site; TES, Transcription End Site; kb, kilobases; pre miRNA, precursor microRNA; rRNA, ribosomal RNA; lnc RNA, long non-coding RNA; snRNA, small nuclear RNA; snoRNA, small nucleolar RNA; scaRNA, small cajal body specific RNA; mature miRNA, mature microRNA; piRNA, piwi interacting RNA.
piRNA are involved in gene silencing of transposons by forming complexes with argonaute proteins (72), and appear to provide RNA-mediated adaptive immunity against transposons (73). In our analysis we found that 208 piRNA were expressed in NOD mouse UC exosome samples and 238 in BALB/c mouse samples. Of these, 171 were found in both but 37 were unique to NOD mouse samples, and 67 unique to BALB/c mouse samples (Supplementary Figure 3B). The 10 most highly expressed piRNA did not appear to be dysregulated in the NOD strain (Supplementary Figure 3A). There were no significant differences in the total number of distinct piRNA or the total number of reads aligning to piRdb between the two strains (Supplementary Figure 3C). We found 13 piRNA to be upregulated in NOD mouse samples with a log2 fold change >3, while 15 were downregulated with a log2 fold change < −3. Of these, only mmu-piR-58696 was significantly upregulated in the NODs as determined by LimmaVoom (padj = 0.047) (Supplementary Figure 3D).
NOD Serum Exosomes Contain a Subset of Dysregulated miRNA
FASTQ reads were preprocessed using FASTqc which identified the presence of adapter in ~90% of the reads. Cutadapt was used to remove adapter sequences and exclude reads of quality <20 and length <15. Pre-processed trimmed reads were then mapped to the miRNA transcriptome from miRbase v22.0 using Bowtie v1.2.3, as well as our in-house aligner miRGrep, which utilizes brute force and was written specifically for the alignment of RNA < 30 nt in length such as miRNA and piRNA. With Bowtie, we identified 550 distinct miRNAs in NODs and 255 in BALB/c. Using miRGrep, we identified 251 miRNAs in the NODs and 242 miRNAs in BALB/c. This is to be expected because miRGrep uses brute-force to align reads to miRNA, with no mismatch allowed in alignment and therefore, the miRNA identified by miRGrep are a subset of those identified by Bowtie (74). Of these, read counts for 38 miRNAs were found only in NOD serum exosomes while 29 miRNAs were found only in BALB/c serum exosomes (Supplementary Figure 4A) in at least 3 out of 5 sample groups per strain. The top 20 expressed miRNA in NOD serum exosomes were overrepresented to the same extent in BALB/c (Supplementary Figure 4B) with no discernible strain specific differences. Of these miRNAs, miR-191-5p, miR-92a-3p, miR-22-3p, miR-16-5p, let-7f-5p, let-7i-5p, miR-26a-5p, miR-30e-5p, miR-186-5p, miR-30d-5p, miR-451a, miR-181a-5p, miR-148a-3p, miR-423-5p, let-7a-5p, and miR-25-3p have been previously reported to be abundant in serum exosomes (75). miR-486-5p (not shown) was the top over-represented miRNA as reported previously (75). It is possible that these miRNAs serve important regulatory functions that are evolutionarily conserved.
The volcano plot in Figure 3A shows the level of differential expression for all expressed miRNA identified in NOD mouse serum. miRNA that had an adjusted p > 0.05 and log2 fold change <3 or >-3 were not considered significant. DESeq2 and Limma determine significance when an adjusted p < 0.05 is reached whereas EdgeR considers a hit significant only when the FDR is <0.1 in RStudio (Table 1). We have reported and assessed miRNA that met our significance criteria by at least one of the statistical packages. Unsupervised hierarchical clustering analysis of top hits using the Euclidean method clustered NOD mouse miRNA samples in the same group, separate from BALB/c samples, as shown by the top tree (Figure 3B), suggesting that the differential expression observed may be attributed to dysregulation in the NOD strain. We performed an additional unsupervised technique to visualize the variability between the two groups. Principal component analysis (PCA) of the 10 samples with ~500 expressed miRNAs revealed that 22% of the variance could be explained by the differences in strain (Figure 3C). All three packages determined that miR-127-3p, miR-409-3p, and miR-540-3p were significantly overexpressed in NOD mouse (Figure 4A). Limma identified miR-410-3p and miR-541-5p to be overexpressed in NOD serum exosomes (Figure 4B).
Figure 3. Differential miRNA expression analysis. (A) The volcano plot highlights miRNAs that were identified in exosomes isolated by differential ultracentrifugation as upregulated or downregulated by more than 3-fold, shown in orange. Upregulated miRNAs that are statistically significantly different are shown in red. (B) Unsupervised hierarchically clustered heatmaps of the most differentially expressed miRNA (with highly expressed miRNAs in red and minimally expressed miRNAs in blue) in the 10 samples. Samples are labeled either N for NOD or B for BALB/c. Individual samples for the top differentially expressed miRNA cluster by strain. (C) PCA plot of all miRNAs shows that 22% of the variance was accounted for by strain.
Table 1. Summary of the 7 most differentially-expressed miRNAs detected in serum UC exosomes from 14-week male BALB/c and NOD mice using three statistical packages in RStudio.
Figure 4. Boxplots of differentially expressed miRNAs in exosomes isolated from 14-week male BALB/c and NOD mouse serum. Reads were aligned with Bowtie and an in-house brute force-based aligner, and differential gene expression analysis was conducted using DESEq2, EdgeR, and LimmaVoom in RStudio. (A) Both DESeq2 and LimmaVoom identified miR-127-3p, miR-409-3p, and miR-540-3p as “hits” with padj < 0.05. (B) Additionally, LimmaVoom identified miR-410-3p and miR-541-5p to be dysregulated with padj < 0.05. miR-329-5p failed to reach statistical significance due to one outlier but exhibited the highest fold change with Log2FC > 5 in EdgeR and DESEq2. Results were obtained from data sets generated from small RNA sequencing of 5 sets of RNA from exosomes from each mouse cohort, each isolated from purified pooled serum exosomes from five groups of age-, gender-, and strain-matched mice, with 5 mice per group.
Furthermore, we found an additional 19 miRNAs that displayed an at least 3-fold higher expression in NOD serum exosomes than in Balb/c and 11 miRNAs that were under-expressed by at least 3-fold in NOD mice (Supplementary Figure 4C). Despite a meaningful fold change these did not reach statistical significance due to an outlier. Of these, miR-329-5p was found to be over 5-fold overexpressed in NOD serum exosomes (DESeq2, EdgeR). As NOD mice, even when age matched, show a variation in disease progression, it is possible that the outlier group may have progressed further in disease than other groups and vice versa. Thus, we included miR-329-5p in our downstream analyses as its differential expression may have biologically significance.
Validation of miRNA Differential Expression
To validate our deep sequencing findings of known miRNAs that were differentially expressed in NOD vs. BALB/c mouse serum exosomes, we purified serum exosomes from independent cohorts of mice, 5 samples per strain with each sample comprised of serum exosomes from each of five mice. For this round of exosome isolation, we used the SEC method (Figures 1C–F) which yielded exosomes of greater particle homogeneity. qRT-PCR was performed using Advanced Taqman miRNA Assays. In agreement with the sequencing data, miRNA miR-127-3p, miR-409-3p, miR-540-3p, miR-410-3p, miR-541-5p, and miR-329-5p were expressed at a higher level in NOD mouse serum exosomes (Figure 5). Five out of the six miRNAs were more than 25-fold over-expressed in the NODs on average. Of these, over-expression of mmu-miR-127-3p and mmu-miR-541-5p were found to be statistically significant (p < 0.01, Mann-Whitney U-test). miR-410-3p and miR-329-5p were upregulated in 4/5 groups whereas miR-409-3p was upregulated in 3/5 groups.
Figure 5. Validation of miRNA ‘hits’ in Size Exclusion Chromatography (SEC)-purified serum exosomes from 14-week male BALB/c and NOD mice. Fold-change data represent mean ± SD normalized to miR-16-5p using the ΔΔCt method. Statistical significance was determined by the Mann-Whitney U-test; **p < 0.01. Results were obtained from data sets generated from qRT-PCR from RNA isolated from SEC-purified pooled serum exosomes from five groups of mice per strain, with 5 mice per group. All mice were age matched.
Pathway Analysis Identifies Several Pathways Involved in Lymphocyte Activation
To understand the signaling pathways the “hit” miRNA may be involved in, we used three methods (miTALOS, miRPathDB, and IPA) which utilize a range of databases of predicted and experimentally validated mRNA-miRNA interactions (such as Tarbase, TargetScan, miRanda, Starbase2) and databases of known pathways (KEGG, Reactome, WikiPathways) in their pathway analysis algorithm. These databases do not have information on all known miRNA for every species, and so we used multiple tools in our pathway analysis to be as comprehensive as possible. Functional over-representation analysis using Gene Ontology (Biological Functions) in miTALOS allows miRNA target prediction using the latest versions of TargetScan (6.2) and miRanda. Additionally, a database of experimentally validated targets—StarBase2- was also implemented in the pathway analysis and is available for both human and murine microRNAs. StarBase2 catalogs 3 miRNA hits identified in this study –mmu-miR-127-3p, mmu-miR-409-3p, mmu-miR-410-3p. As these have identical sequences for human and mice, pathway analysis for both species (Supplementary Figures 5A,B) was conducted using these three hits. Databases in miRPathDB catalog mmu-miR-541-5p in addition to the three miRNAs above. Using miRPathDB, both KEGG (not shown) and WikiPathways identified B cell receptor signaling (Supplementary Figure 5C) while Gene ontology—Biological functions identified lymphocyte activation (not shown). Other pathways identified are associated with cellular proliferation, pluripotency, apoptosis, and p53 signaling, all pathways that may be relevant to cancer. IPA identified several pathways involved in immune regulation (Supplementary Figure 5D). Pathways common to all three analyses include B-Cell Receptor (miR-targeted genes in lymphocyte in miTALOS), TGF-beta and IL-6 Signaling.
Discussion
Here we report that an unbiased screen of exosomal serum miRNA in a murine model of early-intermediate stage autoimmune dacryoadenitis and SS, identified the significant upregulation of multiple miRNAs. qPCR analysis of serum exosomes from a separate set of disease-model vs. healthy control mice validated two significant hits, miR-127-3p and miR-541-3p, while confirming marked elevation with some variance across animal groups for three additional miRNAs, miR-409-3p, miR-410-3p, and miR-329-5p. SS is a complicated autoimmune disease, with treatment hindered by a poorly understood pathogenesis and a lack of reliable diagnostics. Given that miRNAs are master regulators of gene expression with their dysregulation implicated in many diseases (76), identification of this group of dysregulated miRNAs in the male NOD mice at the initial stages of autoimmune dacryoadenitis and development of other indicators of established systemic disease in SS may be useful in establishment of future diagnostic biomarkers.
mmu-miR-127-3p, a significant and validated hit, has been shown to be necessary for the self-renewal and differentiation of hematopoietic stem cells (HSC) in a mouse model of HSC self-renewal defect (77). It has also been proposed to be a regulator of senescence as a tumor suppressor by directly targeting BCL6 (78, 79), a known protooncogene in human cell lines. As BCL6 is a transcriptional repressor and inhibits the production of IL10, its downregulation by miR-127-3p can lead to an increase in IL-10 (80). Increased levels of IL-10 are well-documented in SS patients (81, 82) and are also reported in the NOD mouse lacrimal gland in association with development of autoimmune dacryoadenitis (15). Our pathway analysis shows that mmu-miR-127-3p, is involved in the regulation of TGF-beta and B-Cell receptor Signaling (Supplementary Figure 5C) through its targeting of several MAP kinases and BCL6 (Supplementary Figure 5D). Appropriate regulation of BCL6 is also necessary for the development of germinal center B cell and follicular helper T cells (83). Upregulated levels of hsa-miR-127-3p are reported in testicular and nodal diffused large B-cell lymphoma, with an inverse correlation to BCL6 levels (84). Thus, regulated levels of miR-127-3p are necessary for appropriate control of lymphoproliferation and B cell homeostasis and elevated levels may be indicative of immune dysfunction/autoimmunity. This finding in the NOD mice is of particular relevance because a subset of SS patients develop B cell lymphoma (5) and there is great interest in biomarkers that may distinguish these patients from others with SS so that earlier interventions may be applied to suppress development of B cell lymphoma. It will be of great interest to study the potential dysregulation of miR-127-3p longitudinally in SS patients to explore its relationship to this debilitating and most destructive manifestation of SS.
mmu-miR-541-5p, a second significant and validated hit, may work in concert with mmu-miR-127-3p. In a mouse model of multiple sclerosis, miR-541-5p and miR-127-3p were upregulated in lymph nodes indicating that they may be involved in pathogenic neuro-inflammation (85). Interestingly, knockout of TNF-α in a mouse model led to downregulation of both miR-541-5p and miR-127-3p in epidermal skin, hinting at a close involvement of these miRNAs with pro-inflammatory cytokines (86).
Another hit, miR-409-3p, is broadly implicated in autoimmune disease in animal models and patients with chronic fatigue syndrome/myalgic encephalomyelitis (87), multiple sclerosis (88) and systemic lupus erythematosus (89). Our pathway analysis using miRPathDB and IPA indicates miR-409-3p's involvement in the B Cell Receptor, STAT3 and IL-6 signaling pathways (Supplementary Figures 5C,D). Studies have found that in mice with experimental autoimmune encephalomyelitis (EAE, a murine model of multiple sclerosis), mmu-miR-409-3p targets suppressor of cytokine signaling protein 3 (SOCS3). Upregulation of mmu-miR-409-3p in astrocytes of EAE mice silences SOCS3, leading to an increase in phosphorylation of STAT3 and increased production of inflammatory cytokines such as IL-1β, CXCL10, IL-6, MC-P1 (90). Another study found that in a co-culture of NOD mice salivary gland acinar cells (SGAC) and B-lymphocyte (an in-vitro model of salivary gland disease in SS), there was a significant increase in production of cytokines IL-6 and IL-1β by B-lymphocytes and increased phosphorylation of STAT3 in SGAC (91). IL-1β is upregulated in diseased NOD mouse LG, while its injection into murine LG further impairs tear production (92). Secretion of these cytokines by circulating lymphocytes is also increased in SS patients (93). Thus, mmu-miR-409-3p may be pro-inflammatory in nature and its upregulation may increase cytokine production via the SOCS3/STAT3 signaling pathway.
According to TargetScan, mammals including mice and humans have an 8-mer conserved site on the Stat3 gene for miR-410 (30). hsa-miR-410-3p appears to directly target STAT3, leading to a reduction of IL-10 in T cells of patients with systemic lupus erythematosus (94), and was also shown to be elevated in the plasma of these patients (95). Increased expression of miR-410-3p was also observed in males with relapsing remitting multiple sclerosis, which is characterized by cycling of autoimmune inflammatory status (96). Expression of hsa-miR-410-3p was decreased in the synovial fluid and synoviocytes of rheumatoid arthritis (RA) patients. On the other hand, in an in-vitro model of RA, overexpression of hsa-miR-410-3p decreased the pro-inflammatory cytokines, TNF-α, IL-6, IL-1β, and MMP-9 (97) and was anti-proliferative and apoptotic in nature through targeting of transcription factor YY1 (98). Based on these results, upregulated miR-410-3p may have an immune-protective effect. If validated, use of this miRNA may have value as a potential therapeutic.
Of interest, five of the identified miRNAs (miR-127-3p, miR-329-5p, miR-409-3p, miR-410-3p and miR-541-5p) are encoded on the highly evolutionarily conserved Dlk1-Gtl2 locus on the maternally inherited allele of mouse chromosome 12, which is analogous to the locus Dlk1-Dio3 on the maternally inherited allele of human chromosome 14. The miRNAs from this locus seem to regulate ground state pluripotency in embryonic stem cells (99). Genes from this imprinted locus play a critical role in embryonic and fetal development and appear to be dysregulated in several diseases, including blood cancers such as lymphoma, acute myeloid and acute promyelocytic leukemia, as well as autoimmune diseases such as lupus nephritis (100) and multiple sclerosis (96). It is also interesting that the sequences of mature miRNA 127-3p, 409-3p, and 410-3p are identical in human and mouse. Sequences of miRNA-541-5p and miRNA-329-5p vary only by 3 nucleotides between human and mouse. This further highlights the applicability of miRNA-based biomarkers arising from murine model in this study to humans.
Studies of extracellular vesicles rely on particle size analysis which cannot discriminate between functional vesicles and lipid droplets of similar size, and the presence of these can be a confounding factor in our studies. Future studies will aim to validate these results utilizing different methods of exosome isolation. Although our choice of serum over plasma and the use of serum-separator tubes was aimed at reducing hemolysis, it may also be a confounding factor in our study. While miRNA stability studies in plasma have shown that levels of miR-127-3p are not altered by hemolysis (87), strengthening its potential as a biomarker candidate, similar investigations are needed for the other miRNAs identified in this study.
While we have identified the potential role of some miRNA in inflammatory pathways and SS, further study is required to better understand the relationships of these miRNA as well as the temporal changes that occur with disease progression. The panel of miRNA identified in this study reflect early changes in SS progression, while different patterns of dysregulation may be observed longitudinally as disease advances and/or as it impacts different organs. Since the secretory contents from both lacrimal and salivary glands are also influenced by the autoimmune inflammation characteristic of SS and have been used as sources of biomarkers to reflect both systemic and local inflammation, evaluating these biofluids in later studies as additional sources of potential biomarkers will be of importance. Future studies will also focus on assessing the utility of this panel of miRNA in identifying SS patients as well as exploring the utility of these miRNA in longitudinal disease progression in combination with existing blood-based biomarkers such as anti-La and anti-Ro antibodies and others.
Data Availability Statement
The original contributions presented in the study are publicly available. The datasets can be found in the Sequence Research Archive under ID PRJNA622527 and are available for download here https://www.ncbi.nlm.nih.gov/bioproject/622527.
Ethics Statement
The animal study was reviewed and approved by University of Southern California, Institutional Animal Care and Use Committee (IACUC).
Author Contributions
SK performed the bioinformatics, statistical analysis of next generation sequencing data, isolated and characterized mouse serum exosomes from the validation mouse cohort, isolated exosomal RNA, performed qPCR validation, and prepared the manuscript and supplementary files. SJ and BC isolated and characterized mice serum exosomes and exosomal RNA from the test cohort. SJ assisted with blood collection and performed qPCR validation of miRNA in the validation cohort. BC and SK optimized serum exosome isolation protocols. ME contributed to writing the manuscript and characterization of serum exosomes. DC contributed to data analysis and revised final manuscript. CO and SH-A secured research funding and contributed to experimental design and writing the manuscript. SH-A revised the final manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Institute of Health under Grant number NIH EY011386. Research reported in this publication was also supported by the National Eye Institute of the National Institutes of Health under Award Number P30EY029220 and Unrestricted Grant to the Department of Ophthalmology from Research to Prevent Blindness, New York, NY.
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.
Acknowledgments
The authors thank Anthony Rodriguez for his assistance with TEM grid preparation and imaging, Alexander Yao for his assistance with light microscopy and Hugo Avila in the preparation of the supplementary figure draft.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2020.01475/full#supplementary-material
References
1. Bowman SJ, Ibrahim GH, Holmes G, Hamburger J, Ainsworth JR. Estimating the prevalence among caucasian women of primary Sjogren's syndrome in two general practices in Birmingham UK. Scand J Rheumatol. (2004) 33:39–43. doi: 10.1080/03009740310004676
2. Thomas E, Hay EM, Hajeer A, Silman AJ. Sjogren's syndrome, a community-based study of prevalence and impact. Br J Rheumatol. (1998) 37:1069–76. doi: 10.1093/rheumatology/37.10.1069
3. Brito-Zeron P, Baldini C, Bootsma H, Bowman SJ, Jonsson R, Mariette X, et al. Sjogren syndrome. Nat Rev Dis Primers. (2016) 2:16047. doi: 10.1038/nrdp.2016.47
5. Zufferey P, Meyer OC, Grossin M, Kahn MF. Primary Sjogren's syndrome (SS) and malignant lymphoma. A retrospective cohort study of 55 patients with SS. Scand J Rheumatol. (1995) 24:342–5. doi: 10.3109/03009749509095178
6. Kassan SS, Thomas TL, Moutsopoulos HM, Hoover R, Kimberly RP, Budman DR, et al. Increased risk of lymphoma in sicca syndrome. Ann Intern Med. (1978) 89:888–92. doi: 10.7326/0003-4819-89-6-888
7. Shiboski CH, Shiboski SC, Seror R, Criswell LA, Labetoulle M, Lietman TM, et al. International Sjogren's syndrome criteria working, 2016 American College of Rheumatology/European League against rheumatism classification criteria for primary Sjogren's syndrome, a consensus and data-driven methodology involving three international patient cohorts. Arthritis Rheumatol. (2017) 69:35–45. doi: 10.1136/annrheumdis-2016-210571
8. Martin-Martin LS, Latini A, Pagano A, Ragno A, Stasi R, Coppe A, et al. A new mathematical model based on clinical and laboratory variables for the diagnosis of Sjogren's syndrome. Clin Rheumatol. (2003) 22:123–6. doi: 10.1007/s10067-002-0696-8
9. Pavlidis NA, Karsh J, Moutsopoulos HM. The clinical picture of primary Sjogren's syndrome, a retrospective study. J Rheumatol. (1982) 9:685–90.
10. Humphreys-Beher MG. Animal models for autoimmune disease-associated xerostomia and xerophthalmia. Adv Dent Res. (1996) 10:73–5. doi: 10.1177/08959374960100011501
11. Doyle ME, Boggs L, Attia R, Cooper LR, Saban DR, Nguyen CQ, et al. Autoimmune dacryoadenitis of NOD/LtJ mice and its subsequent effects on tear protein composition. Am J Pathol. (2007) 171:1224–36. doi: 10.2353/ajpath.2007.070388
12. Delaleu N, Nguyen CQ, Peck AB, Jonsson R. Sjogren's syndrome, studying the disease in mice. Arthritis Res Ther. (2011) 13:217. doi: 10.1186/ar3313
13. Edman MC, Janga SR, Meng Z, Bechtold M, Chen AF, Kim C, et al. Increased cathepsin S activity associated with decreased protease inhibitory capacity contributes to altered tear proteins in Sjögren's syndrome patients. Sci Rep. (2018) 8:11044. doi: 10.1038/s41598-018-29411-9
14. Hamm-Alvarez SF, Janga SR, Edman MC, Madrigal S, Shah M, Frousiakis SE, et al. Tear cathepsin S as a candidate biomarker for Sjogren's syndrome. Arthritis Rheumatol. (2014) 66:1872–81. doi: 10.1002/art.38633
15. Li X, Wu K, Edman M, Schenke-Layland K, MacVeigh-Aloni M, Janga SR, et al. Increased expression of cathepsins and obesity-induced proinflammatory cytokines in lacrimal glands of male NOD mouse. Invest Ophthalmol Vis Sci. (2010) 51:5019–29. doi: 10.1167/iovs.09-4523
16. Lavoie TN, Lee BH, Nguyen CQ. Current concepts, mouse models of Sjogren's syndrome. J Biomed Biotechnol. (2011) 2011:549107. doi: 10.1155/2011/549107
17. Sumida T, Tsuboi H, Iizuka M, Asashima H, Matsumoto I. Anti-M3 muscarinic acetylcholine receptor antibodies in patients with Sjogren's syndrome. Mod Rheumatol. (2013) 23:841–5. doi: 10.3109/s10165-012-0788-5
18. Nguyen KH, Brayer J, Cha S, Diggs S, Yasunari U, Hilal G, et al. Evidence for antimuscarinic acetylcholine receptor antibody-mediated secretory dysfunction in nod mice. Arthritis Rheum. (2000) 43:2297–3063. doi: 10.1002/1529-0131(200010)43:10<2297::AID-ANR18>3.0.CO;2-X
19. Shen L, Suresh L, Lindemann M, Xuan J, Kowal P, Malyavantham K Jr, et al. Novel autoantibodies in Sjogren's syndrome. Clin Immunol. (2012) 145:251–5. doi: 10.1016/j.clim.2012.09.013
20. Bahamondes V, Albornoz A, Aguilera S, Alliende C, Molina C, Castro I, et al. Changes in Rab3D expression and distribution in the acini of Sjogren's syndrome patients are associated with loss of cell polarity and secretory dysfunction. Arthritis Rheum. (2011) 63:3126–35. doi: 10.1002/art.30500
21. Meng Z, Edman MC, Hsueh PY, Chen CY, Klinngam W, Tolmachova T, et al. Imbalanced Rab3D versus Rab27 increases cathepsin S secretion from lacrimal acini in a mouse model of Sjogren's syndrome. Am J Physiol Cell Physiol. (2016) 310:C942–54. doi: 10.1152/ajpcell.00275.2015
22. Yasunami R, Bach JF. Anti-suppressor effect of cyclophosphamide on the development of spontaneous diabetes in NOD mice. Eur J Immunol. (1988) 18:481–4. doi: 10.1002/eji.1830180325
23. Ju Y, Janga SR, Klinngam W, MacKay JA, Hawley D, Zoukhri D, et al. NOD and NOR mice exhibit comparable development of lacrimal gland secretory dysfunction but NOD mice have more severe autoimmune dacryoadenitis. Exp Eye Res. (2018) 176:243–51. doi: 10.1016/j.exer.2018.09.002
24. Nguyen C, Singson E, Kim JY, Cornelius JG, Attia R, Doyle ME, et al. Sjogren's syndrome-like disease of C57BL/6.NOD-Aec1 Aec2 mice, gender differences in keratoconjunctivitis sicca defined by a cross-over in the chromosome 3 Aec1 locus. Scand J Immunol. (2006) 64:295–307. doi: 10.1111/j.1365-3083.2006.01828.x
25. Lee HJ, Shin S, Yoon SG, Cheon EJ, Chung SH. The effect of chloroquine on the development of dry eye in sjogren syndrome animal model. Invest Ophthalmol Vis Sci. (2019) 60:3708–16. doi: 10.1167/iovs.19-27469
26. Morthen MK, Tellefsen S, Richards SM, Lieberman SM, Rahimi Darabad R, Kam WR, et al. Testosterone influence on gene expression in lacrimal glands of mouse models of sjogren syndrome. Invest Ophthalmol Vis Sci. (2019) 60:2181–97. doi: 10.1167/iovs.19-26815
27. Su YC, Xiang RL, Zhang Y, Ding C, Cong X, Guo XH, et al. Decreased submandibular adiponectin is involved in the progression of autoimmune sialoadenitis in non-obese diabetic mice. Oral Dis. (2014) 20:744–55. doi: 10.1111/odi.12197
28. Humphreys-Beher MG, Brinkley L, Purushotham KR, Wang PL, Nakagawa Y, Dusek D, et al. Characterization of antinuclear autoantibodies present in the serum from nonobese diabetic (NOD) mice. Clin Immunol Immunopathol. (1993) 68:350–6. doi: 10.1006/clin.1993.1137
29. Bartel DP. MicroRNAs, genomics, biogenesis, mechanism, and function. Cell. (2004) 116:281–97. doi: 10.1016/S0092-8674(04)00045-5
30. Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. (2009) 19:92–105. doi: 10.1101/gr.082701.108
31. Di Leva G, Croce CM. miRNA profiling of cancer. Curr Opin Genet Dev. (2013) 23:3–11. doi: 10.1016/j.gde.2013.01.004
32. Jordan SD, Kruger M, Willmes DM, Redemann N, Wunderlich FT, Bronneke HS, et al. Obesity-induced overexpression of miRNA-143 inhibits insulin-stimulated AKT activation and impairs glucose metabolism. Nat Cell Biol. (2011) 13:434–46. doi: 10.1038/ncb2211
33. Sucharov C, Bristow MR, Port JD. miRNA expression in the failing human heart, functional correlates. J Mol Cell Cardiol. (2008) 45:185–92. doi: 10.1016/j.yjmcc.2008.04.014
34. Li JY, Yong TY, Michael MZ, Gleadle JM. Review: The role of microRNAs in kidney disease. Nephrology. (2010) 15:599–608. doi: 10.1111/j.1440-1797.2010.01363.x
35. Hebert SS, De Strooper B. Alterations of the microRNA network cause neurodegenerative disease. Trends Neurosci. (2009) 32:199–206. doi: 10.1016/j.tins.2008.12.003
36. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. (2008) 105:10513–8. doi: 10.1073/pnas.0804549105
37. Gui Y, Liu H, Zhang L, Lv W, Hu X. Altered microRNA profiles in cerebrospinal fluid exosome in Parkinson disease and Alzheimer disease. Oncotarget. (2015) 6:37043–53. doi: 10.18632/oncotarget.6158
38. Churov AV, Oleinik EK, Knip M. MicroRNAs in rheumatoid arthritis, altered expression and diagnostic potential. Autoimmun Rev. (2015) 14:1029–37. doi: 10.1016/j.autrev.2015.07.005
39. Zhao C, Dong J, Jiang T, Shi Z, Yu B, Zhu Y, et al. Early second-trimester serum miRNA profiling predicts gestational diabetes mellitus. PLoS ONE. (2011) 6:e23925. doi: 10.1371/journal.pone.0023925
40. Rotival M, Siddle KJ, Silvert M, Pothlichet J, Quach H, Quintana-Murci L. Population variation of miRNAs and isomiRs and their impact on human immunity to infection. bioRxiv [Preprint]. (2020) doi: 10.1101/2020.01.31.928580
41. Warnefors M, Liechti A, Halbert J, Valloton D, Kaessmann H. Conserved microRNA editing in mammalian evolution, development and disease. Genome Biol. (2014) 15:R83. doi: 10.1186/gb-2014-15-6-r83
42. Weber JA, Baxter DH, Zhang S, Huang DY, Huang KH, Lee MJ, et al. The microRNA spectrum in 12 body fluids. Clin Chem. (2010) 56:1733–41. doi: 10.1373/clinchem.2010.147405
43. Rozowsky J, Kitchen RR, Park JJ, Galeev TR, Diao J, Warrell J, et al. exceRpt a comprehensive analytic platform for extracellular RNA profiling. Cell Syst. (2019) 8:352–7.e3. doi: 10.1016/j.cels.2019.03.004
44. Koga Y, Yasunaga M, Moriya Y, Akasu T, Fujita S, Yamamoto S, et al. Exosome can prevent RNase from degrading microRNA in feces. J Gastrointest Oncol. (2011) 2:215–22. doi: 10.3978/j.issn.2078-6891.2011.015
45. Lv LL, Cao Y, Liu D, Xu M, Liu H, Tang RN, et al. Isolation and quantification of microRNAs from urinary exosomes/microvesicles for biomarker discovery. Int J Biol Sci. (2013) 9:1021–31. doi: 10.7150/ijbs.6100
46. Miranda KC, Bond DT, McKee M, Skog J, Paunescu TG, Da Silva N, et al. Nucleic acids within urinary exosomes/microvesicles are potential biomarkers for renal disease. Kidney Int. (2010) 78:191–9. doi: 10.1038/ki.2010.106
47. Skog J, Wurdinger T, van Rijn S, Meijer DH, Gainche L, Sena-Esteves M, et al. Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers. Nat Cell Biol. (2008) 10:1470–6. doi: 10.1038/ncb1800
48. Sanz-Rubio D, Martin-Burriel I, Gil A, Cubero P, Forner M, Khalyfa A, et al. Stability of circulating exosomal miRNAs in healthy subjects. Sci Rep. (2018) 8:10306. doi: 10.1038/s41598-018-28748-5
49. National Research Council (US) Committee for the Update of the Guide for the Care and Use of Laboratory Animals. Guide for the Care and Use of Laboratory Animals, 8th ed. Washington, DC. National Academies Press (2011).
50. Shah M, Edman MC, Janga SR, Shi P, Dhandhukia J, Liu S, et al. A rapamycin-binding protein polymer nanoparticle shows potent therapeutic activity in suppressing autoimmune dacryoadenitis in a mouse model of Sjogren's syndrome. J Control Release. (2013) 171:269–79. doi: 10.1016/j.jconrel.2013.07.016
51. Thery C, Amigorena S, Raposo G, Clayton A. Isolation and characterization of exosomes from cell culture supernatants and biological fluids. Curr Protoc Cell Biol. (2006) 30:3.221–29. doi: 10.1002/0471143030.cb0322s30
52. Helwa I, Cai J, Drewry MD, Zimmerman A, Dinkins MB, Khaled ML, et al. A comparative study of serum exosome isolation using differential ultracentrifugation and three commercial reagents. PLoS One. (2017) 12:e0170628. doi: 10.1371/journal.pone.0170628
53. Langmead B, Trapnell C, Pop M, 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
54. Liao Y, Smyth GK, Shi W. Feature counts, an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. (2014) 30:923–30. doi: 10.1093/bioinformatics/btt656
55. Buffalo VS. Sam2Counts. In: Github (2013). Available online at: https://github.com/vsbuffalo/sam2counts
56. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
57. Robinson MD, McCarthy DJ, 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
58. Phipson B, Lee S, Majewski IJ, Alexander WS, Smyth GK. Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression. Ann Appl Stat. (2016) 10:946–63. doi: 10.1214/16-AOAS920
59. Mi QS, Weiland M, Qi RQ, Gao XH, Poisson LM, Zhou L. Identification of mouse serum miRNA endogenous references by global gene expression profiles. PLoS ONE. (2012) 7:e31278. doi: 10.1371/journal.pone.0031278
60. Song J, Bai Z, Han W, Zhang J, Meng H, Bi J, et al. Identification of suitable reference genes for qPCR analysis of serum microRNA in gastric cancer patients. Dig Dis Sci. (2012) 57:897–904. doi: 10.1007/s10620-011-1981-7
61. Matouskova P, Bartikova H, Bousova I, Hanusova V, Szotakova B, Skalova L. Reference genes for real-time PCR quantification of messenger RNAs and microRNAs in mouse model of obesity. PLoS ONE. (2014) 9:e86033. doi: 10.1371/journal.pone.0086033
62. Rinnerthaler G, Hackl H, Gampenrieder SP, Hamacher F, Hufnagl C, Hauser-Kronberger C, et al. miR-16-5p Is a stably-expressed housekeeping microRNA in breast cancer tissues from primary tumors and from metastatic sites. Int J Mol Sci. (2016) 17:156. doi: 10.3390/ijms17020156
63. Kowarsch A, Preusse M, Marr C, Theis FJ. miTALOS, analyzing the tissue-specific regulation of signaling pathways by human and mouse microRNAs. RNA. (2011) 17:809–19. doi: 10.1261/rna.2474511
64. Kehl T, Kern F, Backes C, Fehlmann T, Stockel D, Meese E, et al. miRPathDB 2.0, a novel release of the miRNA pathway dictionary database. Nucleic Acids Res. (2020) 48:D142–7. doi: 10.1093/nar/gkz1022
65. Thery C, Boussac M, Veron P, Ricciardi-Castagnoli P, Raposo G, Garin J, et al. Proteomic analysis of dendritic cell-derived exosomes, a secreted subcellular compartment distinct from apoptotic vesicles. J Immunol. (2001) 166:7309–18. doi: 10.4049/jimmunol.166.12.7309
66. Mol EA, Goumans MJ, Doevendans PA, Sluijter JPG, Vader P. Higher functionality of extracellular vesicles isolated using size-exclusion chromatography compared to ultracentrifugation. Nanomedicine. (2017) 13:2061–65. doi: 10.1016/j.nano.2017.03.011
67. Lobb RJ, Becker M, Wen SW, Wong CS, Wiegmans AP, Leimgruber A, et al. Optimized exosome isolation protocol for cell culture supernatant and human plasma. J Extracell Vesicles. (2015) 4:27031. doi: 10.3402/jev.v4.27031
68. Goetzl EJ, Boxer A, Schwartz JB, Abner EL, Petersen RC, Miller BL, et al. Altered lysosomal proteins in neural-derived plasma exosomes in preclinical Alzheimer disease. Neurology. (2015) 85:40–7. doi: 10.1212/WNL.0000000000001702
69. Zhu Y, Chen X, Pan Q, Wang Y, Su S, Jiang C, et al. A comprehensive proteomics analysis reveals a secretory path- and status-dependent signature of exosomes released from tumor-associated macrophages. J Proteome Res. (2015) 14:4319–31. doi: 10.1021/acs.jproteome.5b00770
70. Potolicchio I, Carven GJ, Xu X, Stipp C, Riese RJ, Stern LJ, et al. Proteomic analysis of microglia-derived exosomes, metabolic role of the aminopeptidase CD13 in neuropeptide catabolism. J Immunol. (2005) 175:2237–43. doi: 10.4049/jimmunol.175.4.2237
71. Wei Z, Batagov AO, Schinelli S, Wang J, Wang Y, El Fatimy R, et al. Coding and noncoding landscape of extracellular RNA released by human glioma stem cells. Nat Commun. (2017) 8:1145. doi: 10.1038/s41467-017-01196-x
72. Seto AG, Kingston RE, Lau NC. The coming of age for Piwi proteins. Mol Cell. (2007) 26:603–9. doi: 10.1016/j.molcel.2007.05.021
73. Duc C, Yoth M, Jensen S, Mouniee N, Bergman CM, Vaury C, et al. Trapping a somatic endogenous retrovirus into a germline piRNA cluster immunizes the germline against further invasion. Genome Biol. (2019) 20:127. doi: 10.1186/s13059-019-1736-x
74. Kakan SS. miRGrep. Github (2019). Available online at: https://github.com/singhkakan/TRGN_510_Final_Project
75. Juzenas S, Venkatesh G, Hubenthal M, Hoeppner MP, Du ZG, Paulsen M, et al. A comprehensive, cell specific microRNA catalogue of human peripheral blood. Nucleic Acids Res. (2017) 45:9290–301. doi: 10.1093/nar/gkx706
76. Ardekani AM, Naeini MM. The role of microRNAs in human diseases. Avicenna J Med Biotechnol. (2010) 2:161–79.
77. Crisafulli L, Muggeo S, Uva P, Wang Y, Iwasaki M, Locatelli S, et al. MicroRNA-127-3p controls murine hematopoietic stem cell maintenance by limiting differentiation. Haematologica. (2019) 104:1744–55. doi: 10.3324/haematol.2018.198499
78. Chen J, Wang M, Guo M, Xie Y, Cong YS. miR-127 regulates cell proliferation and senescence by targeting BCL6. PLoS ONE. (2013) 8:e80266. doi: 10.1371/journal.pone.0080266
79. Saito Y, Liang G, Egger G, Friedman JM, Chuang JC, Coetzee GA, et al. Specific activation of microRNA-127 with downregulation of the proto-oncogene BCL6 by chromatin-modifying drugs in human cancer cells. Cancer Cell. (2006) 9:435–43. doi: 10.1016/j.ccr.2006.04.020
80. Sawant DV, Sehra S, Nguyen ET, Jadhav R, Englert K, Shinnakasu R, et al. Bcl6 controls the Th2 inflammatory activity of regulatory T cells by repressing Gata3 function. J Immunol. (2012) 189:4759–69. doi: 10.4049/jimmunol.1201794
81. Perrier S, Serre AF, Dubost JJ, Beaujon G, Plazonnet MP, Albuisson E, et al. Increased serum levels of interleukin 10 in Sjogren's syndrome; correlation with increased IgG1. J Rheumatol. (2000) 27:935–9.
82. Bertorello R, Cordone MP, Contini P, Rossi P, Indiveri F, Puppo F, et al. Increased levels of interleukin-10 in saliva of Sjogren's syndrome patients. Correlation with disease activity. Clin Exp Med. (2004) 4:148–51. doi: 10.1007/s10238-004-0049-9
83. Nurieva RI, Chung Y, Martinez GJ, Yang XO, Tanaka S, Matskevitch TD, et al. Bcl6 mediates the development of T follicular helper cells. Science. (2009) 325:1001–5. doi: 10.1126/science.1176676
84. Robertus JL, Harms G, Blokzijl T, Booman M, de Jong D, van Imhoff G, et al. Specific expression of miR-17-5p and miR-127 in testicular and central nervous system diffuse large B-cell lymphoma. Mod Pathol. (2009) 22:547–55. doi: 10.1038/modpathol.2009.10
85. Bergman P, James T, Kular L, Ruhrmann S, Kramarova T, Kvist A, et al. Next-generation sequencing identifies microRNAs that associate with pathogenic autoimmune neuroinflammation in rats. J Immunol. (2013) 190:4066–75. doi: 10.4049/jimmunol.1200728
86. Singh A, Willems E, Hafeez BB, Ong IM, Mehta SL, Verma AK. Ultraviolet radiation-induced tumor necrosis factor alpha, which is linked to the development of cutaneous SCC, modulates differential epidermal microRNAs expression. Oncotarget. (2016) 7:17945–56. doi: 10.18632/oncotarget.7595
87. Brenu EW, Ashton KJ, Batovska J, Staines DR, Marshall-Gradisnik SM. High-throughput sequencing of plasma microRNA in chronic fatigue syndrome/myalgic encephalomyelitis. PLoS ONE. (2014) 9:e102783. doi: 10.1371/journal.pone.0102783
88. Ebrahimkhani S, Vafaee F, Young PE, Hur SSJ, Hawke S, Devenney E, et al. Exosomal microRNA signatures in multiple sclerosis reflect disease status. Sci Rep. (2017) 7:14293. doi: 10.1038/s41598-017-14301-3
89. Chen JQ, Papp G, Poliska S, Szabo K, Tarr T, Balint BL, et al. MicroRNA expression profiles identify disease-specific alterations in systemic lupus erythematosus and primary Sjogren's syndrome. PLoS ONE. (2017) 12:e0174585. doi: 10.1371/journal.pone.0174585
90. Liu X, Zhou F, Yang Y, Wang W, Niu L, Zuo D, et al. MiR-409-3p and MiR-1896 co-operatively participate in IL-17-induced inflammatory cytokine production in astrocytes and pathogenesis of EAE mice via targeting SOCS3/STAT3 signaling. Glia. (2019) 67:101–12. doi: 10.1002/glia.23530
91. Xu T, Xie W, Ma Y, Zhou S, Zhang L, Chen J, et al. Leptin/OB-R pathway promotes IL-4 secretion from B lymphocytes and induces salivary gland epithelial cell apoptosis in Sjögren's syndrome. Oncotarget. (2017) 8:63417–29. doi: 10.18632/oncotarget.18823
92. Zoukhri D, Hodges RR, Byon D, Kublin CL. Role of proinflammatory cytokines in the impaired lacrimation associated with autoimmune xerophthalmia. Invest Ophthalmol Vis Sci. (2002) 43:1429–36.
93. Roescher N, Tak PP, Illei GG. Cytokines in Sjogren's syndrome. Oral Dis. (2009) 15:519–26. doi: 10.1111/j.1601-0825.2009.01582.x
94. Liu D, Zhang N, Zhang X, Qin M, Dong Y, Jin L. MiR-410 down-regulates the expression of interleukin-10 by targeting stat3 in the pathogenesis of systemic lupus erythematosus. Cell Physiol Biochem. (2016) 39:303–15. doi: 10.1159/000445625
95. Navarro-Quiroz E, Pacheco-Lugo L, Navarro-Quiroz R, Lorenzi H, Espana-Puccini P, Diaz-Olmos Y, et al. Profiling analysis of circulating microRNA in peripheral blood of patients with class IV lupus nephritis. PLoS ONE. (2017) 12:e0187973. doi: 10.1371/journal.pone.0187973
96. Baulina N, Osmak G, Kiselev I, Popova E, Boyko A, Kulakova O, et al. MiRNAs from DLK1-DIO3 imprinted locus at 14q32 are associated with multiple sclerosis, gender-specific expression and regulation of receptor tyrosine kinases signaling. Cells. (2019) 8:133. doi: 10.3390/cells8020133
97. Wang Y, Xu N, Zhao S, Jiao T, Fu W, Yang L, et al. miR-410-3p suppresses cytokine release from fibroblast-like synoviocytes by regulating NF-kappaB signaling in rheumatoid arthritis. Inflammation. (2019) 42:331–341. doi: 10.1007/s10753-018-0896-2
98. Wang Y, Jiao T, Fu W, Zhao S, Yang L, Xu N, et al. miR-410-3p regulates proliferation and apoptosis of fibroblast-like synoviocytes by targeting YY1 in rheumatoid arthritis. Biomed Pharmacother. (2019) 119:109426. doi: 10.1016/j.biopha.2019.109426
99. Moradi S, Sharifi-Zarchi A, Ahmadi A, Mollamohammadi S, Stubenvoll A, Gunther S, et al. Small RNA sequencing reveals Dlk1-Dio3 locus-embedded microRNAs as major drivers of ground-state pluripotency. Stem Cell Rep. (2017) 9:2081–96. doi: 10.1016/j.stemcr.2017.10.009
Keywords: Sjögren's Syndrome, diagnostic miRNA biomarkers, extracellular vesicles, small-RNA sequencing, microRNA, piwi-RNA
Citation: Kakan SS, Janga SR, Cooperman B, Craig DW, Edman MC, Okamoto CT and Hamm-Alvarez SF (2020) Small RNA Deep Sequencing Identifies a Unique miRNA Signature Released in Serum Exosomes in a Mouse Model of Sjögren's Syndrome. Front. Immunol. 11:1475. doi: 10.3389/fimmu.2020.01475
Received: 01 April 2020; Accepted: 05 June 2020;
Published: 17 July 2020.
Edited by:
Ger Rijkers, University College Roosevelt, NetherlandsReviewed by:
Janine Adele Lamb, University of Manchester, United KingdomTao Li, National Center of Biomedical Analysis (NCBA), China
Copyright © 2020 Kakan, Janga, Cooperman, Craig, Edman, Okamoto and Hamm-Alvarez. 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: Sarah F. Hamm-Alvarez, shalvar@usc.edu