Mapping lncRNAs onto multilevel mRNA co‑expression modules in autism.
TL;DR
The findings establish MSCN as a scalable framework for decomposing ASD mRNA co-expression architecture, provide a hypothesis-generating resource linking coding and noncoding transcriptomic alterations, and help prioritize these relationships for future validation.
Abstract
Autism spectrum disorder (ASD) involves heterogeneous genetic and transcriptomic alterations, but how these changes are organized across hierarchical co-expression scales and linked to long noncoding RNAs (lncRNAs) remains incompletely understood. Here, we present Minimum Span Clustering Network (MSCN), an unsupervised, deterministic framework that constructs traceable multilevel mRNA co-expression hierarchies without requiring soft-thresholding powers, fixed module numbers, cut heights, or stochastic initialization. Applied to two independent ASD brain transcriptome cohorts, MSCN reveals hierarchical gene modules across four resolution levels, enabling the detection of transcriptional patterns ranging from low-level, specific signals to high-level, broader biological pathways. We uncovered modules enriched in neuronal/axonal, developmental, and immune-related pathways, reflecting interconnected neurodevelopmental and immune-dysregulation programs in ASD. Cross-method comparisons showed that MSCN complements WGCNA and MEGENA by preserving biologically concordant modules while providing an explicit parent-child hierarchy; simulation and benchmark analyses supported comparable module recovery, preservation, and enrichment performance. By mapping lncRNAs to MSCN-derived mRNA modules, we identified lncRNAs associated with ASD-related mRNA modules and evaluated them as candidate statistical mediators in downstream transcription factor (TF)-lncRNA-mRNA analyses. This analysis identified over 11,000 candidate TF-lncRNA-mRNA axes, including 644 axes involving 46 SFARI score 1 or 2 genes, four of which were syndromic genes, and 12 named lncRNA mediators. External transcriptomic evaluation further supported 789 axes, including representative HIF1A-STXBP5-AS1-CADPS, E2F1-PART1-SCN2A, and RELA-RFPL1S-GRIN2A relationships. Together, these findings establish MSCN as a scalable framework for decomposing ASD mRNA co-expression architecture, provide a hypothesis-generating resource linking coding and noncoding transcriptomic alterations, and help prioritize these relationships for future validation.