1. Introduction
Aging has been defined as a progressive organic functional decline with loss of homeostasis and increasing probability of illness and death [
1]. All age-dependent diseases are connected to aging and these diseases have an influence on one another. Some of the mechanisms of aging have been revealed by modern biology [
1,
2], and the data suggest that passage of time alone is not the best measure of aging. Therefore, development of aging biomarkers is necessary to clearly define the aging and disease processes [
3].
MicroRNAs (miRNAs) are single-stranded non-coding RNAs of approximately 21–23 nucleotides in length that can affect gene expression at the posttranscriptional level by binding to the three prime untranslated region (3′-UTR) of target messenger RNAs [
4]. MiRNAs are present in several body fluids including amniotic fluid, breast milk, bronchial lavage, cerebrospinal fluid, colostrum, peritoneal fluid, plasma, pleural fluid, saliva, seminal fluid, tears and urine [
5]. It has been reported that miRNAs are important regulators of cellular senescence and aging [
6]. It is also known that serum miRNAs can act as diagnostic and prognostic biomarkers related to the biological and pathological processes of aging [
7]. These data suggest that circulating miRNAs may be useful as aging biomarkers. However, there are few reports on the relationships between aging and miRNAs in body fluids other than blood.
Saliva contains a variety of molecular and microbial analytes, and may be an effective indicator of both local and systemic conditions [
8]. It has been reported that blood-derived molecules entering salivary tissues via transcellular or paracellular routes affect the molecular composition of oral fluids [
8]. It is also estimated that miRNA expression profiles in saliva are similar to those in serum [
9]. Thus, miRNAs obtained in human saliva may also be considered as candidate aging biomarkers. In the present study, we hypothesized that salivary exosomal miRNAs are useful aging biomarkers. Therefore, the aim of this pilot study was to examine whether salivary exosomal miRNAs could be identified or not as aging biomarkers.
3. Discussion
To the best of our knowledge, this is the first report to investigate possible aging biomarkers obtained from salivary exosomal miRNAs. The analyses revealed that miR-24-3p was identified as a novel candidate aging biomarker. Moreover, oral conditions could affect the discrimination power of miR-24-3p. Especially in the low salivary flow rate group and high PPD group, using the expression rates of miR-24-3p could correctly discriminate the old group from the young group. This pilot study indicated that salivary exosomal miRNAs could serve as aging biomarkers.
Aging is correlated with increased prevalence and severity of periodontitis [
10]. In this study, we recruited the participants without severe periodontitis. However, there were statistically significant correlations between miR-24-3p and oral conditions in old participants. It is feasible that miR-24-3p could reflect periodontal inflammation at a subclinical level in the elderly.
Inflamm-aging, which is chronic low-grade inflammation during aging, is commonly characterized by pro-inflammatory biomarkers [
11]. Our enrichment analysis indicated that the unique target genes of miR-24-3p could strongly affect MAPK signaling pathway, which is involved in inflammatory cytokine and chemokine gene regulation at both the transcriptional and the post-transcriptional levels [
12]. The increase in miR-24-3p may contribute to increased susceptibility to age-dependent alterations in the immune and inflammatory status.
The present result also showed that expression rate of salivary miR-24-3p was negatively correlated with salivary flow rate in old participants but not in young participants. This suggests that expression rates of salivary miR-24-3p reflect age-related decline in salivary function. In our findings, target genes of miR-24-3p might be associated with apoptosis. It has been reported that an increased number of apoptotic salivary epithelial cells may contribute to the decreased saliva production and excretion in aging mouse [
13]. It is suggested that miR-24-3p could induce apoptosis of salivary glands during aging, and such conditions might decrease salivary functions in the elderly. On the other hand, exosomes are released from many cell types including dendritic cells, lymphocytes, platelets mast cells, epithelial cells, endothelial cells, and neurons, and can be found in several types of body fluids including saliva [
14]. It is still unclear which cell types released miR-24-3p and affect whether systemic aging and/or tissue-specific aging including oral tissue and salivary gland. Further studies are needed to clarify this point.
The unique target genes and modified NOD value could be varied with the set of exosomal miRNAs. Moreover, all of the target genes could be related to aging processes. Therefore, to clarify the more detailed function of miR-24-3p for aging processes, regulatory power of miR-24-3p should be analyzed with the different set of exosomal miRNAs obtained by different body fluids or tissues in the future study.
The reasons why miR-4739 was considered as suitable for internal control miRNA were as follows. In samples of 15 young participants and 13 old participants, miR-4739 represented consistent expression (raw Ct ≤ 35). Mann-Whitney U test indicated that raw Ct value of miR-4739 were stable in young and old groups (p = 0.339). Moreover, the expression of miR-4739 and the amounts of total RNA were significantly correlated with Spearman’s rank correlation test (r = 0.450, p = 0.016).
Studies have investigated the relationships between miRNA profiles and aging. For instance, a clinical study reported that serum expression levels of miR-151a-3p, miR-181a-5p and miR-1248 were significantly lower in old participants (mean age, 64 years) than in young participants (mean age, 30 years) [
15]. Another clinical study using serum samples also showed that PCR analysis identified five down-regulated miRNAs (miR-29b, miR-106b, miR-130b, miR-142-5p and miR-340) and three up-regulated miRNAs (miR-92a, miR-222 and miR-375) with age [
7]. Among our findings, miR-24-3p in saliva was confirmed to show significant age-related increases. The present and previous findings support the notion that circulating miRNAs are useful for aging biomarkers. However, in our study, overlap of the relative expression profile of miR-24-3p were also found in the young and the old groups. This suggests that miR-24-3p might be a non-specific screening biomarker for aging.
Expression of miRNA profiles in saliva may be similar to those in serum. However, the results of the present and previous studies indicate that there seems to be little overlap of age-associated miRNAs between serum and saliva. The miRNAs in saliva could be produced by components of the salivary gland. Therefore, our results may more reflect the aging change of the salivary gland than a whole body.
The limitations of this pilot study are as follows. External validity is limited because all participants were recruited at the Department of Preventive Dentistry, Okayama University Hospital. Correlation between expression level of miR-24-3p and age are still unclear because the participants did not cover wide range of age and this study design was cross-sectional. Moreover, subclinical status might affect these results especially in the old group.
4. Experimental Section
4.1. Study Participants
This study was in accordance with the Declaration of Helsinki [
16] and was approved by the Ethics Committee of Okayama University Hospital (approval number: 1975). After obtaining informed consent, participants who fulfilled the study requirements were enrolled. We considered that the minimum sample size for this pilot study was 12 per group, because there was no prior information to base a sample size on [
17]. Fifteen young participants (median age, 21.0 years) and 13 old participants (median age, 66.0 years) were finally recruited at the Department of Preventive Dentistry, Okayama University Hospital from August 2014 to December 2014. Exclusion criteria were as follows: cancer; diabetes; chronic lung disease other than asthma; cardiovascular event in the past six months; renal dysfunction; hepatic dysfunction; dementia; severe immune alterations; autoimmune diseases; current smoker; and pregnancy at the time of saliva collection. All of these were as reported in previous studies related to aging [
18,
19,
20]. In addition, severe periodontitis [
21] was set as an exclusion criterion.
4.2. Questionnaire
Participants filled out the self-administered questionnaire for assessment of systemic conditions. Weight and height of participants were also recorded and BMI (weight in kilograms per square of height in meters) was calculated for each participant.
4.3. Oral Examination
PPD and CAL were determined at six sites (mesio-buccal, mid-buccal, disto-buccal, mesio-lingual, mid-lingual and disto-lingual) on all teeth using a color-coded probe (Hu-Friedy; Chicago, IL, USA). Sites that bled upon gentle probing (25 g probing force) were recorded, and the proportion of sites with BOP was measured. PCR were measured after staining with erythrosine, and were recorded in terms of presence or absence adjacent to the gingival margin at four sites (mesial, distal, buccal and lingual) around each tooth [
22]. All clinical procedures were performed by four trained dentists (Takaaki Tomofuji, Takayuki Maruyama, Tatsuya Machida, and Yuya Kawabata).
4.4. Saliva Collection
Unstimulated whole saliva was collected as reported previously with minor modification [
23]. Briefly, saliva was collected in the morning (7:00 AM to 12:00 noon). During the collection period, participants were seated straight up and were instructed to refrain from speaking or swallowing. They allowed the saliva to accumulate in the floor of the mouth, and then spit it through a funnel into a tube kept on ice. At least 1 mL of unstimulated whole saliva was collected. Salivary flow rate was calculated by dividing the volume of collected saliva by the duration of collection time. After collection of saliva, it was stored at 4 °C for up to 6 h, after which it was stored at −80 °C until use.
4.5. Exosome Isolation
Exosomes were isolated from saliva samples (0.5–1.0 mL) using total exosome isolation reagent (Invitrogen, Carlsbad, CA, USA), in accordance with the manufacturer’s protocol. Briefly, whole saliva was centrifuged at 2000× g (g indicates gravity) for 10 min at room temperature to remove cells and debris. Supernatant was incubated with total exosome isolation reagent at 2 to 8 °C for 1 h. After incubation, samples were centrifuged at 10,000× g for 1 h at 2 to 8 °C. Residual supernatant was then discarded and the exosome pellet was collected. The pellet was then re-suspended in a 0.1 volume (calculated by salivary supernatant volume) of exosome re-suspension buffer (Invitrogen).
4.6. RNA Extraction
Total exosome RNA and protein isolation kits (Invitrogen) were utilized for extraction of total RNA from exosome samples. The exosome pellet sample was incubated for 5–10 min at room temperature in order to allow the pellet to dissolve and 200 μL of each sample (brought up to volume with exosome re-suspension buffer if necessary) was combined with 200 μL of 2× denaturing solution, and was incubated on ice for 5 min. After incubation, 400 μL of acid-phenol chloroform was added to the mixture, followed by vortexing for 30–60 s. Samples were centrifuged for 5 min at 10,000× g at room temperature to separate the mixture into the aqueous (upper) phase, interphase and lower phase. The 250 μL of aqueous phase, which contains total RNA, was removed and transferred into a new tube. Then, 99.5% ethanol was added to the aqueous phase, up to 700 μL was placed onto a spin column in a collection tube, and it was centrifuged at 10,000× g for 15 s. After discarding the flow-through and saving the spin column, the sample was washed three times with wash solution mixed with 99.5% ethanol (centrifuged for 15 s for each wash). After washing, the filter was dried by a further 1 min centrifuging at 10,000× g and the filter cartridge was transferred into a fresh collection tube. Next, 100 μL of preheated (95 °C) nuclease-free water was applied to the center of the filter and the sample was centrifuged for 30 s at 10,000× g to recover total RNA. We repeated this procedure three times and the elution containing the RNA was collected. It was stored at 4 °C for up to 6 h, after which it was stored at −80 °C until use. After extraction of total RNA from salivary exosome, we confirmed the quality of them using the Agilent 2100 bioanalyzer (Agilent Technologies, Carlsbad, CA, USA) and the Agilent RNA 6000 Pico Kit (Agilent Technologies).
4.7. Microarray Analysis
A microarray analysis using 3D-Gene (Toray, Kanagawa, Japan) was performed. We prepared two pooled RNA samples by mixing equal amounts (20 ng) of total RNA, which were calculated using the Agilent 2100 bioanalyzer (Agilent Technologies), obtained by five young and five old participants, respectively. These samples had at least 100 pg/μL of total RNA and were mostly consisted of small RNA (less than 500 nucleotide). Samples were labeled with the 3D-Gene miRNA labeling kit (Toray) and hybridized onto 3D-Gene human miRNA Oligo chips v20 in accordance with the supplier’s protocols (
http://www.3d-gene.com). Fluorescent signals were scanned with the 3D-Gene Scanner and analyzed using 3D-Gene Extraction software (Toray). Raw data for each spot was normalized by substitution with a mean intensity for background signals, which were determined from blank spots’ signal intensities of 95% confidence intervals. Measurement of spots with signal intensities greater than 2 standard deviations (SD) from the background signal intensity were considered to be valid. Background-subtracted signal intensities were calculated for these spots. We used the global normalization method (ratio median = 1). We calculated
Ct values of old samples for each miRNA using the signals of young samples as reference.
4.8. Selecting Candidate miRNAs for Aging Biomarkers
We used the DIANA miRPath v.2.0. [
24] to predict target genes and pathways of the miRNAs. To evaluate the regulatory power involved in aging processes, we calculated the novel out degree (NOD) values [
25] with modification. NOD values have been proposed to measure the independent regulatory power of individual miRNAs. In other words, NOD values represents the number of uniquely regulated genes by one specific miRNA [
25]. We considered the miRNAs with high regulatory power, which means the miRNAs could uniquely (in our miRNAs set) regulate many target genes involved in aging processes, as suitable for candidate aging biomarkers. We also considered the miRNAs with NOD values (≥4) [
25] as miRNAs with high regulatory power for aging. Several research groups have performed genome-wide analyses of transcriptional profiling of aging in human tissues, including brain, kidney, muscle, skin and blood. Common aging-related signature genes and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [
26] related to these genes were listed in a previous study [
27]. To focus on regulatory power for the aging process, we calculated NOD values using only genes related to these aging pathways. Moreover, the miRNAs commonly detected in microarray and could be analyzed in DIANA miRPath were used to calculate NOD values. These list of miRNAs was uploaded to the online software. The default DIANA-microT-CDS cut-off value of 0.8 was applied for filtering of miRNA-gene interactions. To increase the strictness of the target pathway prediction, a posteriori approach (pathway union) was applied. miRNA and pathway related information was obtained from miRBase 18 and KEGG v58.1, respectively, for homo sapiens.
4.9. RT-qPCR
To validate the microarray data, we performed TaqMan RT-qPCR assays. We used the inventoried TaqMan microRNA Assays (Life Technologies, Carlsbad, CA, USA) for the RT-qPCR analyses performed on Mx3000P Real-time QPCR System (Agilent Technologies) according to the manufacturer’s instructions [
28,
29]. Briefly, reverse transcription (RT) enzymes including specific RT primers for each miRNA and total RNA sample were mixed in a tube. RT run conditions were then set. The plate was incubated at 16 °C for 30 min, 42 °C for 30 min and 85 °C for 5 min, and was then held at 4 °C. Products of the RT reaction were stored at −20 °C until use. RT-qPCR was performed in triplicate using TaqMan Universal Master Mix II, no UNG (Applied Biosystems, Life Technologies, Carlsbad, CA, USA). After holding at 95 °C for 10 min, fifty thermal cycles (95 °C for 15 s and 60 °C for 60 s) were run.
Ct values were determined using background-based threshold (cycle-range 5 to 8) calculated by the MxPro Mx3000P v4.10 software (Stratagene, Edinburgh, UK). Data with a >35 raw
Ct value for each miRNA was excluded from each analysis. We considered the miRNAs with consistent expression in our study as suitable for internal control miRNA because internal control small RNAs for normalization of salivary exosomal miRNAs have not yet clearly decided [
23]. In addition, we considered that the miRNAs with high NOD value, which indicates high regulatory power in aging process, could serve as useful internal control miRNAs. This is based on the concept that identification of multiple miRNA changes could be valuable because specific characteristics of miRNA combinations unique to a normal physiological or pathological state could serve as useful reference [
30]. In our study, miR-4739 was considered as suitable for internal control miRNAs. We calculated the relative expression rates of each miRNA between young and old groups using 2
−∆∆Ct method [
31]. To determine the statistical significance of relative expression levels of each miRNA between young and old groups, we performed the bootstrapping approach (Pair Wise Fixed Reallocation Randomisation Test) using REST 2009 software with 2000 random reallocations of samples [
32]. This statistical test could be suitable for our analysis because it could avoid making any assumptions about distributions [
32].
4.10. Relationships of miRNA and Clinical Parameters
It has been reported that periodontitis and obesity affect the expression of circulating miRNAs [
33]. Therefore, we analyzed the relationships between candidate miRNA and clinical parameters including age, BMI, and oral conditions in young and old groups using Spearman’s rank correlation test. Furthermore, we analyzed the relationship between candidate miRNA and significantly correlated variables using multiple linear regression analysis with stepwise method. The relative expression of candidate miRNA was normalized by log2 transformation for the analysis.
4.11. Discrimination Power of Candidate miRNA Biomarker for Aging
To evaluate the discrimination power of candidate miRNA biomarker for aging, ROC curves [
34] were constructed for detection of old group. The ROC curves were obtained by plots of the sensitivity (true positive rate) and 1-specificity (false positive rate) of the diagnostic test at various cut-off values. The ROC curves graphically show the relationship between both the true positive rate and false positive rate according to the various cut-off values, which are effective to determine the optimal cut-off values. AUC summarizes the entire location of the ROC curve and reflects the discrimination power of the biomarkers. To decide the optimal threshold value, we used the Youden index (sensitivity + specificity − 1) [
35].
4.12. Enrichment Analysis for Biological Process and Pathway of Target Genes
An enrichment analysis for biological process and pathway of target genes was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 (
http://david.abcc.ncifcrf.gov) [
36] using the homosapiens background list. It was pointed out that multiple testing correction techniques are conservative approaches and it could hurt the sensitivity of discovery if overemphasizing them [
36]. Therefore, we considered that the gene ontology (GO) terms and KEGG pathways with
p value (without multiple testing correction) <0.05 also warranted discussion (even though their corrected
p value was not significant).
4.13. Statistical Analyses
Mann-Whitney U test, chi-square test and Fischer’s exact test were used for comparison of the characteristics of young and old groups. Values of two-sided p < 0.05 were considered to represent statistically significant differences. These statistical analyses were performed using SPSS software version 20 (SPSS Inc., Chicago, IL, USA).