Hot water recirculation systems (HWRS) in hotels and nursing homes, which are common in countries such as Spain, have been related to outbreaks of legionellosis. To establish the relationships of microbial and physicochemical parameters, especially protozoa, with the occurrence of Legionella in HWRS, 231 samples from hotels and nursing homes were analysed for Legionella, protozoa, heterotrophic plate counts (HPC) at 22 and 37 °C, Pseudomonas, metals, temperature and others. Legionella pneumophila was the dominant species isolated, and 22 % were sg. 1. The sampling method became particularly important in order to define which factors were involved on the occurrence of Legionella. Results showed that the bacteria and the accompanying microbiota were more abundant in the first flush water whose temperature was lower. The bacteria occurred in those samples with high HPC and were inversely correlated with high temperatures. Multivariate regression showed that a concentration above 1 × 10(5) CFU/100 mL of HPC at 37 °C, Fe above 0.095 ppm and the presence of protozoa increased significantly the risk of Legionella colonization, while univariant regression showed that the presence of Cu above 0.76 ppm and temperature above 55 °C diminished it. Therefore, to reduce the risk associated with Legionella occurrence in HWRS these parameters should be taken into consideration.