Background: Recent advances in resting-state fMRI allow us to study spatial dynamics, the phenomenon of brain networks spatially evolving over time. However, most dynamic studies still use subject-specific, spatially-static nodes. As recent studies have demonstrated, incorporating time-resolved spatial properties is crucial for precise functional connectivity estimation and gaining unique insights into brain function. Nevertheless, estimating time-resolved networks poses challenges due to the low signal-to-noise ratio, limited information in short time segments, and uncertain identification of corresponding networks within and between subjects. Methods: We adapt a reference-informed network estimation technique to capture time-resolved spatial networks and their dynamic spatial integration and segregation. We focus on time-resolved spatial functional network connectivity (spFNC), an estimate of network spatial coupling, to study sex-specific alterations in schizophrenia and their links to multi-factorial genomic data. Results: Our findings are consistent with the dysconnectivity and neurodevelopment hypotheses and align with the cerebello-thalamo-cortical, triple-network, and frontoparietal dysconnectivity models, helping to unify them. The potential unification offers a new understanding of the underlying mechanisms. Notably, the posterior default mode/salience spFNC exhibits sex-specific schizophrenia alteration during the state with the highest global network integration and correlates with genetic risk for schizophrenia. This dysfunction is also reflected in high-dimensional (voxel-level) space in regions with weak functional connectivity to corresponding networks. Conclusions: Our method can effectively capture spatially dynamic networks, detect nuanced SZ effects, and reveal the intricate relationship of dynamic information to genomic data. The results also underscore the potential of dynamic spatial dependence and weak connectivity in the clinical landscape.