Late-life depression (LLD) is a common and serious mental disorder, whose neural mechanisms are not yet fully understood. In this study, we aimed to characterize LLD-related changes in intrinsic functional brain networks using a large, multi-site sample.Using resting-state functional magnetic resonance imaging, the edge-based functional connectivity (FC) as well as multiple topological brain network metrics at both global and nodal levels were compared between 206 LLD patients and 210 normal controls (NCs).Compared with NCs, the LLD patients had extensive alterations in the intrinsic brain FCs, especially significant decreases in FCs within the default mode network (DMN) and within the somatomotor network (SMN). The LLD patients also showed alterations in several global brain network metrics compared with NCs, including significant decreases in global efficiency, local efficiency, clustering coefficient, and small-worldness, as well as a significantly increased characteristic path length. Moreover, significant alterations in nodal network metrics (increased nodal betweenness and decreased nodal efficiency) were found in patients with LLD, which mainly involved the DMN and SMN. Post-hoc subgroup analyses indicated that the above changes in FC strengths were present in both first-episode, drug-naïve (FEDN) and non-FEDN patients, and were correlated with depression severity in the FEDN patients. Moreover, changes in FC strengths were found in both the early/late-onset (depression starts before/after the age of 50) patients, while altered topological metrics were found in only the late-onset patients.These results may help to strengthen our understanding of the underlying neural mechanisms and biological heterogeneity in LLD.