Python package to determine best clustering solutions with stability-based relative clustering validation Highlights d Implementation of a stability-based relative validation technique for clustering d Best number of clusters selection for multiple tasks d Automation of stable clustering labels on new data via supervised learning d Complement or alternative to internal validation measures for clustering

Social-communication (SC) and restricted repetitive behaviors (RRB) are autism diagnostic symptom domains. SC and RRB severity can markedly differ within and between individuals and may be underpinned by different neural circuitry and genetic mechanisms. Modeling SC-RRB balance could help identify how neural circuitry and genetic mechanisms map onto such phenotypic heterogeneity. Here, we developed a phenotypic stratification model that makes highly accurate (97–99%) out-of-sample SC = RRB, SC > RRB, and RRB > SC subtype predictions. Applying this model to resting state fMRI data from the EU-AIMS LEAP dataset (n = 509), we find that while the phenotypic subtypes share many commonalities in terms of intrinsic functional connectivity, they also show replicable differences within some networks compared to a typically-developing group (TD). Specifically, the somatomotor network is hypoconnected with perisylvian circuitry in SC > RRB and visual association circuitry in SC = RRB. The SC = RRB subtype show hyperconnectivity between medial motor and anterior salience circuitry. Genes that are highly expressed within these networks show a differential enrichment pattern with known autism-associated genes, indicating that such circuits are affected by differing autism-associated genomic mechanisms. These results suggest that SC-RRB imbalance subtypes share many commonalities, but also express subtle differences in functional neural circuitry and the genomic underpinnings behind such circuitry.