Incomplete lineage sorting (ILS), the interaction between coalescence and speciation, is a 1 process which generates incongruences between gene trees and species trees, as do gene 2 duplication (D), transfer (T) and loss (L). These processes are usually modelled 3 independently, but in reality, ILS can affect gene copy number polymorphism, i.e., interfere 4 with DTL. This has been previously recognised, but not treated in a satisfactory way, 5 mainly because DTL events are naturally modelled forward-in-time, while ILS requires 6 backward-in-time coalescent modelling. Here we consider the joint action of ILS and DTL on 7 the gene tree/species tree problem in all its complexity. In particular, we show that the 8 interaction between ILS and duplications/transfers only can result in patterns usually 9 interpreted as resulting from gene loss, and that the realised rate of D, T and L becomes 10 non-homogeneous in time when ILS is taken into account. We introduce algorithmic 11 solutions to these problems. Our new model, the multilocus multispecies coalescent 12 (MLMSC), which also accounts for any level of linkage between loci, generalises the 13 multispecies coalescent model and offers a versatile, powerful framework for proper 14 simulation and inference of gene family evolution. 15 1 2 Q. LI, N. GALTIER, C. SCORNAVACCA, AND Y. CHAN