BackgroundAdolescent depression shows high clinical heterogeneity. Brain functional networks serve as a powerful tool for investigating neural mechanisms underlying depression profiles. A key challenge is to characterize how variation in brain functional organization links to behavioral features and psychosocial environmental influences. MethodsWe recruited 80 adolescents with major depressive disorder (MDD) and 42 healthy controls (HCs). First, we estimated the differences in functional connectivity of resting-state networks (RSN) between the two groups. Then, we used sparse canonical correlation analysis to characterize patterns of associations between RSN connectivity and symptoms, cognition, and psychosocial environmental factors in MDD adolescents. Clustering analysis was applied to stratify patients into homogenous subtypes according to these brain-behavior-environment associations. ResultsMDD adolescents showed significantly hyperconnectivity between the ventral attention and cingulo-opercular networks compared with HCs. We identified one reliable pattern of covariation between RSN connectivity and clinical/environmental features in MDD adolescents. In this pattern, psychosocial factors, especially the interpersonal and family relationships, were major contributors to variation in connectivity of salience, cingulo-opercular, ventral attention, subcortical and somatosensory-motor networks. Based on this association, we categorized patients into two subgroups which showed different environment and symptoms characteristics, and distinct connectivity alterations. These differences were covered up when the patients were taken as a whole group. ConclusionThis study identified the environmental exposures associated with specific functional networks in MDD youths. Our findings emphasize the importance of the psychosocial context in assessing brain function alterations in adolescent depression and have the potential to promote targeted treatment and precise prevention.