The grass subfamily Poöideae arose 30–40 million years ago, in conjunction with a monophyletic clade of systemic endophytic fungi in the genus Epichloë. These fungi colonize most aboveground grass tissues, but are usually absent from root tissues. Some Epichloë species cause ‘choke disease’ whereby the fungus suppresses the host’s seed production, using the culms as a site to produce fungal ascospores that disperse to infect new host plants. However, many Epichloë are also able to asymptomatically colonize the grass ovary, ovule and embryo without damaging the seed, and transmit themselves vertically from one host generation to the next. Indeed, some Epichloë exhibit choke disease on some tillers, and seed colonization on other tillers. Thus, these species simultaneously undergo sexual recombination combined with horizontal transmission between hosts, and clonal reproduction coupled with vertical transmission from one host generation to the next. Still other species do not reproduce sexually, and transmission to the next generation is restricted to vertical transmission (Schardl 2010).

The most recent taxonomic revision of this genus (Leuchtmann et al. 2014) shows that it comprises 34 species, 3 subspecies, and 6 varieties (i.e., 43 distinct lineages). Whether these species reproduce: only sexually (2 species plus 1 subspecies), only asexually (23 species plus 5 varieties), or both sexually and asexually (9 species plus 2 subspecies) largely determines their mode(s) of transmission into the next generation. Sexual reproduction in these species provides the ability for horizontal transmission of the fungus from one individual host plant to another. Strictly asexually reproducing species seem to be limited exclusively to vertical transmission into the next generation. See Table 1 for details.

Sexual reproduction often involves the production of a stromatal stage that ‘chokes’ some of the host plant’s reproductive tillers (see Fig. 1). It is often presumed that these interactions are parasitic because of the obvious fitness loss suffered by the host plant. However, most, if not all of the Epichloë lineages produce secondary metabolites, alkaloids, that can defend the host plant against some herbivores and pathogens (see, e.g., Schardl et al. 2013). So it is unclear what the balance of lifetime costs and benefits is in these sexually reproducing lineages. Thus, it is an open question as to whether these lineages are parasites or mutualists, and whether they may be either or both in different environments.

As previously mentioned, the asexual species seem to be restricted to vertical transmission. Endosymbionts that are strictly vertically transmitted ought to evolve into mutualists if they are to maintain themselves in the population (see Section 1.3). Indeed, benefits accruing to the host plants were demonstrated ‘early and often’ (for review of the early literature see, e.g., Saikkonen et al. 1998), and the interaction was widely asserted to be a mutualism. This view started to be challenged in the mid-1990s when Stan Faeth and his students began studying the native grass Festuca arizonica and its asexual endophyte Neotyphodium starii.¹ From early on in their research, Faeth and his colleagues commonly failed to find the widely reported benefits to the host plants of the endophyte infection (see, e.g., Lopez, Faeth, and Miller 1995) and indeed found costs of harboring the endophyte. The lack of agreement needed an explanation, and one possible reason for the difference that Faeth and colleagues hit upon early on was that they were studying a ‘native’ grass–endophyte pair in natural settings, while much of the early work on the epichloids in general was done on introduced agricultural grasses paired with their nonnative Epichloë species in managed agroecosystems. Indeed, one can see hints of this explanation forming in the last lines of Faeth’s earliest paper on the subject:

Our results, however, suggest the effect on herbivores may not be as widespread as currently believed (Clay, Hardy, and Hammond 1985; Carrol 1988; Strong 1988), or the effects are more variable in native grass–herbivore interactions. (Lopez, Faeth, and Miller 1995, 1579, emphasis added)

As Faeth’s body of work began to grow in this field, seemingly so did his conviction that what he and his students were seeing was not the mutualism predicted by evolutionary theory and supported by empirical evidence in the agricultural systems. He began publishing papers with provocative titles such as:

“Neotyphodium in Native Populations of Arizona Fescue: A Nonmutualist?” (Faeth et al. 1997),

“Neotyphodium in Arizona Fescue: A Necessary Symbiont?” (Sullivan and Faeth 1999),

“ What Maintains High Levels of Neotyphodium Endophytes in Native Grasses? A Dissenting View and Alternative Hypotheses” (Faeth, Sullivan, and Hamilton 2000),

“Fungal Endophytes: Common Host Plant Symbionts but Uncommon Mutualists” (Faeth and Fagan 2002),

“Are Endophytic Fungi Defensive Plant Mutualists?” (Faeth 2002),

culminating in the particularly assertive:

“Mutualistic Asexual Endophytes in a Native Grass Are Usually Parasitic” (Faeth and Sullivan 2003).

He continued to develop these ideas throughout these papers, for example:

Endophytic fungi, especially asexual, systemic endophytes in grasses, are generally viewed as plant mutualists, mainly through the action of mycotoxins, such as alkaloids in infected grasses, which protect the host plant from herbivores. Most of the evidence for the defensive mutualism concept is derived from studies of agronomic grass cultivars, which may be atypical of many endophyte–host interactions. I argue that endophytes in native plants, even asexual, seed-borne ones, rarely act as defensive mutualists. In contrast to domesticated grasses where infection frequencies of highly toxic plants often approach 100%, natural grass populations are usually mosaics of uninfected and infected plants. The latter, however, usually vary enormously in alkaloid levels, from none to levels that may affect herbivores. This variation may result from diverse endophyte and host genotypic combinations that are maintained by changing selective pressures, such as competition, herbivory and abiotic factors. Other processes, such as spatial structuring of host populations and endophytes that act as reproductive parasites of their hosts, may maintain infection levels of seed-borne endophytes in natural populations, without the endophyte acting as a mutualist. (Faeth 2002, 25) ²

As Faeth has continued to add to this argument through the years, others too started detecting costs of the endophyte to the host grass (e.g., Cheplick, Perera, and Koulouris 2000; Cheplick 2007, 2008; Härri, Krauss, and Müller 2008a, 2008b; Fuchs et al. 2013). Taken together, Faeth’s arguments and the publication of other papers demonstrating costs have given rise to the (albeit minority) view that the strictly asexual vertically transmitted Epichloë –grass interactions form a ‘mutualism–parasitism continuum’ (Müller and Krauss 2005). In a systematic survey of the literature (323 papers, see Appendix A), we found that for strictly vertically transmitted Epichloë species, 58% of the papers considered the interaction to be a mutualism, but 30% considered it at least possible that the interaction could indeed be a parasitism (see Table 2 for details). And this view persists today (see, e.g., Raman and Suryanarayanan 2017).

Cases of a parasite exhibiting only vertical transmission are interesting from an evolutionary perspective. To see why, let: ρ be the proportional change in host fecundity relative to an uninfected host; ψ be the proportional change in the probability of survival from seed through to reproduction; γ be the probability an infected mother will infect an average seed; and δ be the same probability for the father. Then, from the Fundamental Vertical Transmission Equation, it can be shown that a sufficient condition for the persistence of a parasite exhibiting only vertical transmission is ρψ(γ+δ)>1, and that outside of the possibly rare condition that γ+δ>1.65, it is also a necessary condition (Fine 1975). By definition, a parasite reduces the lifetime fitness of its host relative to an uninfected conspecific, so ρψ<1. The condition for a vertically transmitted parasite to maintain itself in a population thus requires a sufficiently high transmission rate (γ+δ) to counterbalance the fitness detriment it causes to its host. In the case of the Epichloë species, the fungus has never been shown to be transmitted via the pollen, and therefore δ=0. Even assuming maternal vertical transmission were perfect (i.e., γ=1) it is not possible for the fungus to maintain itself in the population while still causing a reduction in the lifetime fitness of the host species.³ That is, the fungus cannot continue to be a parasite, and evolutionary theory thus strongly suggests that such species must evolve into mutualists (Ewald 1987).

In light of such theory, those who would argue that these strictly vertically transmitted fungi are not mutualists must provide a mechanism that explains how the fungus is maintained within a plant population, in the face of lifetime fitness costs to the plant. Faeth and Sullivan (2003) offer four possible explanations for their claim that “mutualistic asexual endophytes in a native grass are usually parasitic” (see also Faeth 2010). We briefly summarize these here. See Faeth and Sullivan for a more detailed presentation.

Ewald (1987) probably introduced the term mutualism–parasitism continuum, but Johnson et al. seem to have popularized the term in their consideration of mychorrhizal fungi–host plant interactions. They described it thus (Johnson, Graham, and Smith 1997, 583):

Mutualism and parasitism are extremes of a dynamic continuum of species interactions. Mycorrhizal associations are generally at the mutualistic end of the continuum, but they can be parasitic when the stage of plant development or environmental conditions make costs greater than benefits, or possibly when the genotypes of the symbionts do not form ‘win–win’ associations. In natural systems, plant genotypes exist because they successfully propagate more plants, and fungal genotypes exist because they successfully propagate more fungi. Usually, by living together, plants and mycorrhizal fungi improve each other’s probability for survival and reproductive success. But sometimes, and being anthropomorphic, plant ‘interests’ are in conflict with those of the fungi. In the words of Dawkins (1978) “we must expect lies and deceit, and selfish exploitation of communication to arise whenever the interests of the genes of different individuals diverge.” The ‘interests’ of plants and mycorrhizal fungi are likely to diverge in highly managed agricultural systems, where fertilization eliminates shortages of soil nutrients, and plant genotypes are selected by humans and not by millennia of natural selection.

Although the argument about where on this supposed continuum these epichloid species sit had been going on for quite some time (for a sampling of this literature, see: Clay 1988, 1991; Saikkonen et al. 1999; Faeth et al. 1997; Faeth and Fagan 2002; Faeth, Sullivan, and Hamilton 2000; Faeth, Helander, and Saikkonen 2004; Faeth and Hamilton 2006; Rudgers et al. 2010; Rudgers et al. 2012), Müller and Krauss (2005) seem to have introduced the use of the term to the study of epichloids with their paper titled: “Symbiosis between grasses and asexual fungal endophytes”:

The symbiosis between vertically transmitted asexual endophytic fungi and grasses is common and generally considered to be mutualistic. Recent studies have accumulated evidence of negative effects of endophytes on plant fitness, prompting a debate on the true nature of the symbiosis. Genetic factors in each of the two partners show high variability and have a range of effects (from positive to negative) on plant fitness. In addition, interacting environmental factors might modify the nature of the symbiosis. Finally, competition and multitrophic interactions among grass consumers are influenced by endophytes, and the effects of plant neighbours or consumers could feedback to affect plant fitness. We propose a mutualism–parasitism continuum for the symbiosis between asexual endophytes and grasses, which is similar to the associations between plants and mycorrhizal fungi. (450; emphasis added)

While debate over the usefulness of the concept continues in the mycorrhizal literature (see, e.g., Smith and Smith 2013), in this paper we argue that the concept is useful to the Epichloë situation only when discussing the genus as a whole, or when discussing individual species that possess horizontal transmission. We argue that, contra Müller and Krauss, the term is not appropriate for use with strictly vertically transmitted Epichloë. There are important biological differences between mycorrhizal fungi and Epichloë fungi that make this continuum concept less universal than it is perhaps in the mycorrhizal situation. As Ewald (1987) pointed out, transmission mode is probably the most important difference. In the next section we develop a simple model of each transmission strategy and use this model to draw conclusions about the nature of these Epichloë –grass interactions.

Let 0≤x≤1 be the proportion of a pasture that is infected with the endophyte. Let 0≤y≤1 be the proportion of the pasture that is uninfected. At equilibrium, we denote these as x* and y*. We can then define (1) ϕ*=x*x*+y*, as the proportion of the plants in the pasture that are infected with the endophyte at equilibrium. Note that x+y need not equal 1; the difference simply represents bare ground.

We assume that both birth and death rates are density dependent (e.g., Bullock, Hill, and Silvertown 1994). Let b be the density-independent per capita birth rate, and d be the densityindependent death rate. Note that the model we will develop lacks a stage-structure. Here, ‘birth’ means the rate at which individual plants successfully recruit from seeds (i.e., it includes seed and seedling mortality), while ‘death’ means the death of mature plants. We then model the density-dependent per capita birth and death rates, respectively, as: (2) (1−x−y)b, and (3) (x+y)d.

Throughout the remainder of this paper we refer to the effects of the presence of the endophyte as either an advantage or a disadvantage to the host plant in terms of the effects the endophyte has, directly or indirectly, on the host’s birth rate and death rate. It would also be appropriate to refer to these as benefits and costs of harbouring the endophyte. We use the term (dis)advantage to refer to the effect of the endophyte on the host generally, to denote circumstances in which the these effects might be positive or negative.

Let a be the per capita birth rate (dis)advantage. We then model the endophyte infected proportion of the pasture’s per capita birth rate as: (4) (1−x−y)ab.

If a=1 then there is neither an advantage nor a disadvantage to the host plant of harboring the endophyte, in terms of the plant’s birth rate. If a>1 then the endophyte is an advantage to the plant (at least in terms of its birth rate) because it increases the birth rate. For example, if a=1.2 then the per capita birth rate of infected host plants will be 20% greater than the per capita birth rate of uninfected host plants. And if a<1 then the endophyte imposes a disadvantage on the host plant, at least in terms of its birth rate.

Similarly, let c be the per capita death rate (dis)advantage. We model the endophyte infected proportion of the pasture’s per capita death rate as: (5) (x+y)cd.

If c=1 the endophyte imposes neither an advantage nor a disadvantage to the host plant. If c<1 then the endophyte is an advantage to the host plant, at least in terms of its per capita death rate, and if c>1 then the endophyte confers a disadvantage to the host plant in terms of death rates.⁵ These relationships are summarized in Fig. 2.

Let 0≤ℓ≤1 be the per capita loss rate of the endophyte, and conversely, 1−ℓ is the transmission efficiency rate. Endophyte loss occurs through several non-mutually exclusive mechanisms. The endophyte may fail to colonize a newly produced vegetative tiller, it may fail to colonize the developing seeds prior to dispersal, and it may be lost from successfully colonized seeds prior to germination (Afkhami and Rudgers 2008). Because we are modelling individual plants, we are unconcerned with the failure of the fungus to colonize some tillers on an otherwise infected plant. We are only concerned with the failure of the fungus to give rise to a new, infected plant. This may happen because the seed produced is not successfully colonized prior to formation and dispersal, or the fungus may die in the seed prior to germination, or the fungal infection may fail to establish itself at the seedling stage prior to developing into a mature plant. We subsume all of these mechanisms into the single parameter ℓ. Thus, a positive loss rate means that the proportion ℓ of the per capita birth rate from infected individuals will turn out to be uninfected. And 1−ℓ will turn out to have the infection. The special case of ℓ=0 implies that the endophyte is always perfectly transmitted to the next generation. When ℓ=1 the endophyte is strictly horizontally transmitted.

In the case of strictly horizontal transmission, we assume that all births from infected plants are initially uninfected. Furthermore, we assume that the horizontal transmission rate depends on the interactions between infected and uninfected individuals. For any encounter, the chance of passing on the infection is given by h (conceptually similar to β in standard susceptible– infected (SI) epidemiological models; see, e.g., Anderson and May 1992). Thus the rate of horizontal transmission is given by: hxy. A useful extension to our model would be to relax this assumption by creating a spatially explicit model with a more mechanistically detailed treatment of horizontal transmission.

We can use the above to construct a single model of the Epichloë –grass interaction as follows: (9) dxdt=(1−ℓ)(1−x−y)atbx︸(𝟙)+hxy︸(𝟚)−(x+y)cdx︸(𝟛),dydt=ℓ(1−x−y)atbx︸(𝟜)+(1−x−y)by︸(𝟝)−hxy︸(𝟚)−(x+y)dy︸(𝟞).

Here, (𝟙) denotes the birth rate at which new infected individuals arise from infected plants, (𝟚) denotes the rate of horizontal transmission by which uninfected plants become infected, (𝟛) denotes the density-dependent death rate of infected plants, (𝟜) denotes the birth rate at which new uninfected individuals arise from infected plants, (𝟝) denotes the birth rate at which uninfected individuals arise from uninfected plants, and (𝟞) denotes the density-dependent death rate of uninfected individuals. Note that if ℓ=1 we obtain a model that represents strictly horizontal transmission, and if h=0 we obtain a model that represents strictly vertical transmission. Also note that if we let at=a we obtain a model with temporally constant endophyte (dis)advanatges.

We next use this model to develop each of the three transmission cases separately. For each case we first develop the results for the case of at=a (i.e., temporally constant birth and death (dis)advantages) and then we numerically solve for the temporally varying conditions.

This solution occurs in the case of h/d<c. In this case the infection cannot sustain itself because infected plants die faster than the infection can be passed on to susceptible plants. The equilibrium pasture densities are given by Eq. (10). Again, note that x+y need not equal unity.

When c≤h/d, the equilibrium proportions of infected and uninfected plants in the pasture are given by Eq. (11).

Substituting Eq. (11) into Eq. (1) and simplifying yields: (12) ϕ*=1−cdh.

Note that this solution is independent of any birth rate (dis)advantage. Since there is no vertical transmission, any effect the endophyte has on the host plant’s birth rate only influences the birth rate of uninfected plants. Equation 12 is shown in Fig. 4. We see that the larger the value of horizontal transmission rate (h) or the smaller the value of the death rate (dis)advantage (c), the larger the equilibrium fraction of the pasture that is infected (ϕ*). This matches our intuition.

This solution requires that the denominator in Eq. (11) must not be zero. Given the biological constraints on the parameters listed above, this condition is always satisfied. It can also be shown (see Appendix B) that these solutions are all stable fixed points.

We simulated Eq. (9) for the three patterns of temporal variation in at depicted in Fig. 3i. We simulated 100 years of population dynamics, which was more than sufficient to flush the initial conditions and for the system to settle into its equilibrium dynamics. These results are shown in Fig. 4 and examples of the dynamics are shown in Fig. 5a, 5d and 5g. From these figures we see that, although the equilibrium dynamics are cyclic (due to the temporally varying at) the cycles oscillate around the equilibrium from the stable point solutions found in the case where at=a is a constant. We can see this more generally in Fig. 4, where the point equilibria for the constant (dis)advantage case overlay the stable cycles for the temporally varying at.

If ℓ=0 it can be shown that the following three solutions exist.

If, and only if, the per capita birth rate (dis)advantage equals the per capita death rate (dis)advantage, the system yields a coexistence solution—one where both infected and uninfected plants stably coexist in the population. That is, when a=c we obtain: (15) y*=b−bx*−dx*b+d.

In the case of imperfect endophyte transmission (i.e., 0<ℓ) the system has two just solutions.

We again simulated Eq. (9) for three patterns of variation in at depicted in Fig. 3i. Again, we simulated the dynamics for 100 years which was sufficient to flush the initial conditions and for the system to settle into its equilibrium dynamics. These results are shown in Fig. 6a, 6d and 6g and examples of the dynamics are shown in Fig. 5b, 5e and 5h. Note particularly that in Fig. 6a and 5b (i.e., the endophyte is a parasite because both 1<c and a¯t<1) the endophyte goes locally extinct, even when at is not constant, demonstrating the conclusion that a strictly vertically transmitted endophyte cannot persist as a parasite. As was the conclusion for Case A, we can see from the figures that, although the equilibrium dynamics are cyclic in the case of coexistence (due to the temporally varying at) the cycles oscillate around the equilibrium from the stable point solutions found in the case where at=a is a constant. We see this more generally in Fig. 6a, 6d and 6g, where the equilibria point solutions for the constant (dis)advantage case (solid and dash-dotted line) overlay the cycles for the temporally varying at (yellow and purple shading).

Again, we evaluated Eq. (9) numerically for the values of at shown in Fig. 3i, for a period of 100 years, which was more than adequate to flush the initial conditions and for the system to settle into its long-term dynamics. Examples of these dynamics are shown in Fig. 5c, 5f and 5i. There, we see that the solution to Eq. (9) results in a stable cycle while the solution to Eq. (9), where a=a¯t, results in a stable point around which the cyclic solutions oscillate. More generally, the results are shown in Fig. 6. Note that the solutions to Eq. (9) are cyclic, and are thus illustrated in Fig. 6 as shaded areas that denote the minimum and maximum value of the oscillations. The solutions to Eq. (9) using at are nearly identical to those found using a¯t and the same general conclusions hold.

We think that the confusion surrounding the ‘mutualism–parasitism continuum’ is illustrated in the questions we asked earlier about Fig. 3. Here we return to those questions and provide responses. For the purposes of this discussion, imagine that the (dis)advantages shown in Fig. 3 are observations from actual plants and endophytes. Here we explore the inferences that can be drawn from such ‘observations.’ It is important to note that Fig. 3 does not specify the mode of transmission of the endophyte, but knowing the mode is critical to the inferences we can draw from such ‘results.’ The ‘data’ must be interpreted in light of the theoretical considerations. The two jointly determine the appropriate inferences.

All models are simplifications or abstractions. We think our abstraction captures the most important general mechanisms and therefore offers a reasonable approximation to the biological system it is meant to represent. Others might quibble with this assertion. Certainly there are important aspects of the biology that are not incorporated in our model. For example, there is no representation of the seed-bank dynamics. The model lacks a stage or age structure. The model also has not been explicitly parameterized for any specific host–endophyte pair, nor for any specific geographical location. By keeping the model non-specific in these ways, we think it demonstrates some useful generalizations. Whether adding any of these biological realisms affects the general conclusions remains to be seen, but our intuition is that the general requirements for coexistence will hold, at least in a qualitative way.

One well-known mechanism by which these endophytes may confer an advantage to the host plant is through increased drought tolerance (White et al. 1992). A drought gradient has long been hypothesized to account for the distribution of endophyte-infected perennial ryegrass (Lolium perenne) in Europe (Lewis et al. 1997). We will explore this mechanism in future work.

Another useful extension, not previously considered in any rigorous way, would be to explore the role of environmental stochasticity on coexistence. There are plenty of examples in the literature in which stochastic models lead to different, sometimes deeper, understandings than the corresponding deterministic model (see, e.g., Newman 1991; Black and Mckane 2012).

Ecologist Richard Levins (1966) famously stated:

However, even the most flexible models have artificial assumptions. There is always room for doubt as to whether a result depends on the essentials of a model or on the details of the simplifying assumptions ….

Therefore, we attempt to treat the same problem with several alternative models each with different simplifications but with a common biological assumption. Then, if these models, despite their different assumptions, lead to similar results we have what we can call a robust theorem which is relatively free of independent lies. (423)

…

Unlike the theory, models are restricted by technical considerations to a few components at a time, even in systems which are complex. Thus a satisfactory theory is usually a cluster of models. These models are related to each other in several ways: as coordinate alternative models for the same set of phenomena, they jointly produce robust theorems; as complementary models they can cope with different aspects of the same problem and give complementary as well as overlapping results; as hierarchically arranged “nested” models, each provides an interpretation of the sufficient parameters of the next higher level where they are taken as given. (431)

In this subsection, we review previous modelling attempts for the Epichloë –grass interaction, and compare and contrast our model to the previous work. We conclude that, taken together, this set of models produces a robust theorem about this system: that strictly vertically transmitted Epichloë endophytes must be mutualists.

Clay (1993) developed a very simple fitness-based model, without endophyte loss, to show that for strictly vertically transmitted fungal endophytes with perfect transmission efficiency, any net fitness reduction caused by the endophyte results in the loss of endophyte-infected individuals from the population, echoing earlier work (Ewald 1987; Fine 1975). Clay used this model to show that the equilibrium proportion of infected individuals in the pasture is either 0 or 1, and the only means of obtaining a mixed solution was for the benefit of the endophyte to exactly balance the costs, a condition Clay thinks is unlikely to occur. These results are confirmed again by our work, sensu Eqs. (13) to (15). The addition that our work contributes here is that by placing our model in a population dynamics framework, we are able to demonstrate how, for this very restricted case (a=c) it is possible to get other proportions in the mixture apart from the 12 that Clay finds. We show that the particular mixed result achieved depends on the host plants’ birth and death rates (Eq. (15)).

Ravel, Michalakis, and Charmet (1997) develop a discrete time, stage structured population dynamics model. They divide the endophyte effects on the host plant into effects on birth rate, death rate, and competitive ability. They do not examine cases where the endophyte might be a disadvantage to the host plant in terms of birth and death rates, as we do, although they do briefly admit the possibility that the endophyte might leave the host plant at a competitive disadvantage. The result of this possibility, however, is never explored in their results. In our model, we do not represent competition as distinct from the birth and death rate (dis)advantages. In our model, the competitive advantage of the endophyte, if there is one, arises from how it affects the host plants’ birth and death rates. Finally, Ravel et al. consider a stochastic version of their model where the birth rates, death rates, and competitive advantages are sampled each year from a normal distribution but they provide no details on the mean or variance of these normal distributions, nor do they explore how the mean and variance interact to affect the results. They merely conclude that stochastic effects are only important in small populations, something that has long been known in population ecology. Mostly, Ravel et al. use their model to explore the trade-off between the endophyte loss rate (their μ, conceptually equivalent to our ℓ) and the birth and death rate advantages of the endophyte. In this paper we are primarily interested in how birth and death rate advantages and disadvantages can trade-off against each other to influence the equilibrium infection rate. For us, 0<ℓ is simply a necessary condition to obtain stable mixtures of infected and uninfected plants.

Gundel et al. (2008) modeled a closed annual grass population using a stage-structured, density-independent, periodic, non-stochastic matrix model. We used a non-stage-structured, density-dependent model of a perennial grass population. Like Gundel et al., we provided an example of an environment of deterministic seasonality (see Fig. 3). Also like Gundel et al., we assumed that the grass’ vital rates (b and d) and transmission efficiency (1−ℓ) were constant through time. Gundel et al. used their model to examine how the endophyte loss rate interacts with the ratio of reproductive rates of the infected and uninfected plants to determine the equilibrium infection frequency. Echoing the earlier work, they conclude that:

This condition implies that, in a closed population, the endophyte can only persist in the long term if the infection results in some increase in the host plant fitness …. However, such an increase could be very small …. (900)

Making different assumptions, and structuring our model differently, we find the same result, the combined benefits have to outweigh the costs (i.e., there has to be an increase in the host plant’s net fitness). However, we show that just because one measures a cost in terms of birth or death rates, this does not imply that the endophyte can persist in a population as a parasite. Furthermore, Gundel et al. (2008, 902) make the important point that this net fitness benefit of the endophyte can be so small that it would be difficult to measure under field conditions, yet this small difference could result in high infection frequencies in pastures. We find that this conclusion is true, but only for perfect transmission (i.e., ℓ=0). Once the endophyte loss rate is positive, small net positive effects of the endophyte on the host lead to coexistence, but at a relatively low frequency of infection in the pasture.

Gundel et al. (2008, 897–8) observed that “relatively high infection frequencies have been observed in natural grass populations exhibiting little or no evidence that the fungus confers a reproductive advantage to its host (Saikkonen et al. 1998; Faeth and Hamilton 2006).” Our model provides two possible explanations for this phenomenon. First, such a solution can occur when the endophyte provides a survivorship advantage instead of a reproductive advantage. Second, if costs and benefits of endophyte infection vary within and between years, it would be possible for the pasture to still maintain a high infection rate even if in some seasons, or across some years, the endophyte provides no benefit, or even a short-term net cost. Notice in Fig. 6 there are examples where a<1, i.e., the endophyte is a disadvantage to the host in terms of reproduction, and yet the endophyte still persists in the population, even possibly at relatively high proportions of the pasture depending on the death rate (dis)advantage (c).

Finally, Saikkonen, Ion, and Gyllenberg (2002) use meta-population dynamics to understand how the endophyte can be a parasite in some places and still persist on the landscape. They find that this is indeed possible so long as the seed dispersal rates are high enough, and there exist sufficient ‘empty’ patches into which the parasite-infected host plant can disperse. As we stated in the introduction, we do not find the meta-population dynamics explanation convincing. While long-distance dispersal of grass seeds is possible, aided by animals for instance, the vast majority of seed dispersal in these grasses is quite local. To us it seems unlikely that such restricted dispersal would be sufficient to drive the meta-population explanation. Nevertheless, our model cannot be directly compared to this model as we consider closed populations and ignore immigration and emigration.