How parasite virulence evolves is arguably the most important question in both the applied and fundamental study of host-parasite interactions. Typically, this research area has been progressing through the formalization of the problem via mathematical modelling. This is because the question is a complex one, as virulence is both affected and affects several aspects of the host-parasite interaction. Moreover, the evolution of virulence is a problem in which ecology (epidemiology) and evolution (changes in trait values through time) are tightly intertwined, generating what is now known as eco-evolutionary dynamics. Therefore, intuition is not sufficient to address how virulence may evolve. In their classical model, Anderson and May (1982)[1] predict that the optimal virulence level results from a trade-off between increasing parasite load within hosts and promoting transmission between hosts. Although very useful and foundational, this model incurs into several simplifying assumptions. One of the most obvious is that it considers that hosts are infected by a single parasite strain/species. Some subsequent models have thus accounted for multiple infections, generally predicting that this will select for higher virulence, because it increases the strength of selection in the within-host compartment. Usually, when attacked, hosts deploy defences to combat their parasites. In many systems, however, parasites can suppress the immune response of their hosts. This leads to prolonged infection, which is beneficial for the parasite. However,