We introduce spatially explicit stochastic processes to model multispecies host-symbiont interactions. The host environment is static, modeled by the infinite percolation cluster of site percolation. Symbionts evolve on the infinite cluster through contact or voter type interactions, where each host may be infected by a colony of symbionts. In the presence of a single symbiontspecies, the condition for invasion as a function of the density of the habitat of hosts and the maximal size of the colonies is investigated in details. In the presence of multiple symbiont species, it is proved that the community of symbionts clusters in two dimensions whereas symbiont species may coexist in higher dimensions. From a joint work with N.Lanchier and F.Zucca.