Salivary miRNA profiles identify children with autism spectrum disorder, correlate with adaptive behavior, and implicate ASD candidate genes involved in neurodevelopment

Steven D. Hicks, Cherry Ignacio, Karen Gentile, Frank A. Middleton

Autism spectrum disorder (ASD) is a common neurodevelopmental disorder that lacks adequate screening tools, often delaying diagnosis and therapeutic interventions. Despite a substantial genetic component, no single gene variant accounts for >1 % of ASD incidence. Epigenetic mechanisms that include microRNAs (miRNAs) may contribute to the ASD phenotype by altering networks of neurodevelopmental genes. The extracellular availability of miRNAs allows for painless, noninvasive collection from biofluids. In this study, we investigated the potential for saliva-based miRNAs to serve as diagnostic screening tools and evaluated their potential functional importance. Salivary miRNA was purified from 24 ASD subjects and 21 age- and gender-matched control subjects. The ASD group included individuals with mild ASD (DSM-5 criteria and Autism Diagnostic Observation Schedule) and no history of neurologic disorder, pre-term birth, or known chromosomal abnormality. All subjects completed a thorough neurodevelopmental assessment with the Vineland Adaptive Behavior Scales at the time of saliva collection. A total of 246 miRNAs were detected and quantified in at least half the samples by RNA-Seq and used to perform between-group comparisons with non-parametric testing, multivariate logistic regression and classification analyses, as well as Monte-Carlo Cross-Validation (MCCV). The top miRNAs were examined for correlations with measures of adaptive behavior. Functional enrichment analysis of the highest confidence mRNA targets of the top differentially expressed miRNAs was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID), as well as the Simons Foundation Autism Database (AutDB) of ASD candidate genes. Fourteen miRNAs were differentially expressed in ASD subjects compared to controls (p <0.05; FDR <0.15) and showed more than 95 % accuracy at distinguishing subject groups in the best-fit logistic regression model. MCCV revealed an average ROC-AUC value of 0.92 across 100 simulations, further supporting the robustness of the findings. Most of the 14 miRNAs showed significant correlations with Vineland neurodevelopmental scores. Functional enrichment analysis detected significant over-representation of target gene clusters related to transcriptional activation, neuronal development, and AutDB genes. Measurement of salivary miRNA in this pilot study of subjects with mild ASD demonstrated differential expression of 14 miRNAs that are expressed in the developing brain, impact mRNAs related to brain development, and correlate with neurodevelopmental measures of adaptive behavior. These miRNAs have high specificity and cross-validated utility as a potential screening tool for ASD.