Mitochondria have a central role in cellular (patho)physiology, and they display a highly variable morphology that is probably coupled to their functional state. Here we present a protocol that allows unbiased and automated quantification of mitochondrial 'morphofunction' (i.e., morphology and membrane potential), cellular parameters (size, confluence) and nuclear parameters (number, morphology) in intact living primary human skin fibroblasts (PHSFs). Cells are cultured in 96-well plates and stained with tetramethyl rhodamine methyl ester (TMRM), calcein-AM (acetoxy-methyl ester) and Hoechst 33258. Next, multispectral fluorescence images are acquired using automated microscopy and processed to extract 44 descriptors. Subsequently, the descriptor data are subjected to a quality control (QC) algorithm based upon principal component analysis (PCA) and interpreted using univariate, bivariate and multivariate analysis. The protocol requires a time investment of ∼4 h distributed over 2 d. Although it is specifically developed for PHSFs, which are widely used in preclinical research, the protocol is portable to other cell types and can be scaled up for implementation in high-content screening.