Motivation: Incorporating the temporal dimension into multimorbidity studies has shown to be crucial for achieving a better understanding of the disease associations. Furthermore, due to the multifactorial nature of human disease, exploring disease associations from different perspectives can provide a holistic view to support the study of their aetiology. Results: In this work, a temporal systems-medicine approach is proposed for identifying time-dependent multimorbidity patterns from patient disease trajectories, by integrating data from electronic health records with genetic and phenotypic information. Specifically, the disease trajectories are clustered using an unsupervised algorithm based on dynamic time warping and three disease similarity metrics: clinical, genetic and phenotypic. An evaluation method is also presented for quantitatively assessing, in the different disease spaces, both the cluster homogeneity and the respective similarities between the associated diseases within individual trajectories. The latter can facilitate exploring the origin(s) in the identified disease patterns. The proposed integrative methodology can be applied to any longitudinal cohort and disease of interest. In this article, prostate cancer is selected as a use case of medical interest to demonstrate, for the first time, the identification of temporal disease multimorbidities in different disease spaces.