Citation: Moeller, H. V., M. G. Neubert, and M. D. Johnson. 2019. Intraguild predation enables coexistence of competing phytoplankton in a well-mixed water column. Ecology 100(12):Abstract. Resource competition theory predicts that when two species compete for a single, finite resource, the better competitor should exclude the other. However, in some cases, weaker competitors can persist through intraguild predation, that is, by eating their stronger competitor. Mixotrophs, species that meet their carbon demand by combining photosynthesis and phagotrophic heterotrophy, may function as intraguild predators when they consume the phototrophs with which they compete for light. Thus, theory predicts that mixotrophy may allow for coexistence of two species on a single limiting resource. We tested this prediction by developing a new mathematical model for a unicellular mixotroph and phytoplankter that compete for light, and comparing the model's predictions with a laboratory experimental system. We find that, like other intraguild predators, mixotrophs can persist when an ecosystem is sufficiently productive (i.e., the supply of the limiting resource, light, is relatively high), or when species interactions are strong (i.e., attack rates and conversion efficiencies are high). Both our mathematical and laboratory models show that, depending upon the environment and species traits, a variety of equilibrium outcomes, ranging from competitive exclusion to coexistence, are possible. W Varied December 2019 MIXOTROPHY IN PLANKTONIC COMMUNITIES Article e02874; page 3 FIG. 4. Dependence of mathematical model predicted dynamics on surface light I in and attack rate a. Legends and notation as in Fig. 1. As available light increases, the system transitions from competitive exclusion by prey, to alternate competitive exclusion states or coexistence, to bistability of coexistence and competitive exclusion by the mixotroph, and finally to competitive exclusion by the mixotroph. Nullclines are shown for increasing values of I in at a fixed value of a = 1 9 10 À7 (right column); Roman numerals are used to indicate specific parameter values and are chosen for consistency with Roman numerals in Fig. 1. Note that the basins of attraction for the bistable equilibria shift. For example, as I in increases, the basin of attraction for the coexistence equilibrium (purple square) shrinks (nullclines corresponding to locations IVa and IVb). Parameters are as listed in Table 1, with p m = 0.3, h m = 250, ' m = 0.1, k m = 5 9 10 À7 , and b = 0.005.