# Resource competition shapes biological rhythms and promotes temporal niche differentiation in a community simulation

Vance Gao, Sara Morley-fletcher, Stefania Maccari, Martha Hotz Vitaterna,
Fred W. Turek

► To cite this version:

Vance Gao, Sara Morley-fletcher, Stefania Maccari, Martha Hotz Vitaterna, Fred W. Turek. Resource competition shapes biological rhythms and promotes temporal niche differentiation in a community simulation. Ecology and Evolution, 2020, Ecology and Evolution, 10 (20), pp.11322-11334. 10.1002/ece3.6770. hal-03179288

HAL Id: hal-03179288

https://hal.univ-lille.fr/hal-03179288

Submitted on 24 Mar 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Resource competition shapes biological rhythms and promotes temporal niche differentiation in a community simulation

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

Vance Difan Gao<fnref n="1" />,<fnref n="2" /> ID

Sara Morley-Fletcher<fnref n="1" />,<fnref n="3" /> ID

Stefania Maccari<fnref n="1" />,<fnref n="3" />,<fnref n="4" /> ID

Martha Hotz Vitaterna<fnref n="2" /> ID

Fred W. Turek<fnref n="2" />

## Correspondence

Vance Difan Gao, UMR 8576 Unité de Glycobiologie Structurale et Fonctionnelle, Campus Cité Scientifique, CNRS, University of Lille, Lille, France.

Email: gaodifan@gmail.com

## Funding information

National Heart, Lung, and Blood Institute, Grant/Award Number: HL 007909; I-SITE

## Abstract

1. Competition for resources often contributes strongly to defining an organism's ecological niche. Endogenous biological rhythms are important adaptations to the temporal dimension of niches, but how other organisms influence such temporal niches has not been much studied, and the role of competition in particular has been even less examined. We investigated how interspecific competition and intraspecific competition for resources shape an organism's activity rhythms.

2. To do this, we simulated communities of one or two species in an agent-based model. Individuals in the simulation move according to a circadian activity rhythm and compete for limited resources. Probability of reproduction is proportional to an individual's success in obtaining resources. Offspring may have variance in rhythm parameters, which allow for the population to evolve over time.

3. We demonstrate that when organisms are arrhythmic, one species will always be competitively excluded from the environment, but the existence of activity rhythms allows niche differentiation and indefinite coexistence of the two species. Two species which are initially active at the same phase will differentiate their phase angle of entrainment over time to avoid each other. When only one species is present in an environment, competition within the species strongly selects for niche expansion through arrhythmicity, but the addition of an interspecific competitor facilitates evolution of increased rhythmic amplitude when combined with additional adaptations for temporal specialization. Finally, if individuals preferentially mate with others who are active at similar times of day, then disruptive selection by intraspecific competition can split one population into two reproductively isolated groups separated in activity time.

4. These simulations suggest that biological rhythms are an effective method to temporally differentiate ecological niches and that competition is an important ecological pressure promoting the evolution of rhythms and sleep. This is the first study to use ecological modeling to examine biological rhythms.

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2020 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd

<footnote n="1">¹UMR 8576 Unité de Glycobiologie Structurale et Fonctionnelle, CNRS, University of Lille, Lille, France</footnote>

<footnote n="2">²Center for Sleep and Circadian Biology, Northwestern University, Evanston, IL, USA</footnote>

<footnote n="3">³International Associated Laboratory (LIA) "Perinatal Stress and Neurodegenerative Diseases", University of Lille, Lille, France</footnote>

<footnote n="4">⁴Department of Medico-Surgical Sciences and Biotechnologies, University Sapienza of Rome, Rome, Italy</footnote>
KEYWORDS

biological rhythms, circadian rhythms, competition, ecological modeling, sleep, speciation, temporal niche

# 1 | INTRODUCTION

An ecological niche is the particular "space" a species occupies in an ecosystem, comprising a diversity of concepts including the physical environment where it lives, the resources it uses, its interactions with other species, and the adaptations it has evolved to exploit the niche. Both theory and evidence indicate that competition is a major factor which determines species' niches (Caro & Stoner, 2003; Tilman, 2004). Competitors for a resource generally will adjust to avoid each other, differentiating their niches by adapting to utilize different types of resources, thus avoiding the depletion of shared resources (exploitative competition) (Maurer, 1984) as well as costly acts of aggression and territorial exclusion (interference competition) (Amarasekare, 2002).

Time is one dimension which defines a niche. Most species live in environments dominated by a day-night cycle and a seasonal cycle, though other cycles such as tides can also strongly define an ecosystem (Gibson, 1967). A temporal niche is defined by both abiotic elements such as light intensity and temperature (Ouyang, Andersson, Kondo, Golden, & Johnson, 1998; Yerushalmi & Green, 2009), as well as by predictable temporal patterns of other organisms (Castillo-Ruiz, Paul, & Schwartz, 2012). Adaptions for a specific temporal niche include specializations in vision, coat color, temperature regulation, and importantly, biological rhythms, that is, physiological and behavioral traits which cycle regularly according to an internal biological clock. In this article, we consider activity rhythms through the lens of ecological niche theory, viewing them as a way partition niches among competitors and examining how competition shapes the expression of these biological clocks. Time can thus be viewed as a sort of resource or space which is partitioned by organisms in response to competition (Kronfeld-Schor & Dayan, 2003), both interspecific (between different species) and intraspecific (between members of the same species) (Begon, Townsend, & Harper, 2006). While some work has examined predation's effect on animals' activity rhythms (Bakker, Reiffers, Olff, & Gleichman, 2005; McCauley, Hoffmann, Young, & Micheli, 2012; Reebs, 2002), the role of competition has been less examined.

A variety of field studies have observed evidence of temporal partitioning between putative competitors. For example, a community of three sympatric species of insectivorous bats were observed to have activity profiles which are displaced from each other, suggesting that they avoid each other (Kunz, 1973). While Argentinian pampas foxes are normally nocturnal, they show a modified half-nocturnal half-diurnal profile in a region where they overlap with crab-eating foxes (Di Bitetti, Di Blanco, Pereira, Paviolo, & Pírez, 2009). Other examples of apparent temporal partitioning have been seen at watering holes (Adams & Thibault, 2006; Valeix, Chamaillé-Jammes, & Fritz, 2007), in desert gerbils (Ziv, Abramsky,

Kotler, & Subach, 1993), carabid beetles (Kamenova, Tougeron, Cateine, Marie, & Plantegenest, 2015), and large African carnivores (Hayward & Slotow, 2009). Stronger evidence comes from seminatural manipulation experiments. Common spiny mice are nocturnal and golden spiny mice are diurnal in areas where the two species coexist, but when the common spiny mice are experimentally removed, the originally diurnal golden spiny mice shift into nocturnality, suggesting that this species was pushed into diurnality by the competing species (Kronfeld-Schor & Dayan, 2003). It should be recognized that these examples all occur at the behavioral and ecological timescale, that is, they are "temporary" displacements from the species' genetically determined preferred activity time. Although certain pairs of taxa are suggestive of evolutionary-level temporal partitioning—nocturnal owls and diurnal hawks, for instance, or the hypothesized nocturnal origin of mammals (Heesy & Hall, 2010)—the role of competition in causing these evolutionary divergences is near impossible to establish conclusively. Nonetheless, it is probable that temporal shifts at the behavioral/ecological level cause selection for alleles suited for the new time period, eventually translating into long-term changes at the genetic/evolutionary level.

In this work, we explore ways in which competition in ecological niche theory can be applied to biological rhythms, using agent-based simulation of a community of two species who compete for limited resources in an environment. We first demonstrate that circadian rhythms are a form of niche differentiation, which allows two species to stably coexist without competitive exclusion, and that two rhythmic species will quickly differentiate their circadian phases to avoid competition. We show that arrhythmicity is preferred when there is no interspecific competitor present, but the presence of interspecific competition greatly facilitates the development of strong activity rhythms. Finally, we show that if organisms have assortative mating based on circadian timing, then intraspecific competition can drive the division of one population into two mating groups, contributing to speciation.

# 2 | MATERIALS AND METHODS

## 2.1 | Environment and organisms

A community of two species in an environment was simulated using Python (code available at https://datadryad.org/stash/dataset/doi:10.5061/dryad.gf1vhhmmj). The environment consisted of 60 spaces, where each space may contain a resource, contain an individual organism, contain multiple individuals, or be empty (Figure 1a). An individual moves by jumping to another random space, and if it moves to a space which contains a resource, the resource is consumed. Six hundred resources per day are added randomly to the
<figure><img src="image_4.png" /><figcaption>FIGURE 1 Simulation overview. (a) The environment consists of 60 spaces, which are occupied by resources and moving organisms of one or two species. (b) Activity profiles are modeled as cosine waves with two properties, amplitude and phase angle. Two example profiles are shown. Regardless of changes in circadian parameters, the number of moves per day is fixed at 20 (equivalent to the area under each curve). (c) Individuals gain energy when consuming resources and expend energy on movement, metabolism, and reproduction. Reproductive fitness is proportional to an individual's stored energy. (d) Offspring's amplitude and phase-angle parameters are the mean of the parents and receive 1/3 of the stored energy from each parent. If testing trait evolution, a small normally distributed variance ε is added to the offspring's traits. (e) For some simulations, we added a specialization parameter, representing evolutionary adaptations to a specific time of day. Resource acquisition is more efficient at circadian peak hours and less efficient at the circadian trough. (f) For some simulations, we added assortative mating to test speciation. Individuals only can mate with others who have similar phase angles, and the probability of two individuals mating is proportional to how much their mating windows overlap</figcaption></figure>

spaces of the environment at a constant rate. These resources can represent food but can also represent resources in general including water, territory, and the access to those resources.

An organism's amount of movement varies across the day following a circadian activity rhythm, which is modeled as a cosine wave with two parameters: amplitude, ranging from 0 to 1, and phase angle, ranging from 0° to 360° (Figure 1b). The number of moves performed between times t₀ and t₁ is:

$$m = \int_{t_0}^{t_1} \frac{M_{\text{total}}}{360^{\circ}} (1 + A * \cos(t - \varphi)) dt$$

where m is the number of moves, Mtotal is the total number of moves per day = 20, A is the amplitude, and φ is the phase angle. The number of moves were calculated in 10 time steps per day (t₁-t₀ = 36). We believe a sinusoidal activity rhythm is a reasonable approximation, and a more flexible activity pattern would in fact make temporal segregation easier to achieve. Though an organism's activity profile varies according to its circadian rhythm parameters, the total number of moves is fixed at 20 moves per day; we hold total movement and energy expenditure constant so that the primary variable of interest is the timing of movement. Each animal has the capacity to store energy; 2 units of energy is expended per day for movement and metabolism, and 0.5 energy is gained when the animal consumes a resource (Figure 1c).

Animals die if their energy reaches zero, or if they reach 12 days of age. Values for arena size, energy expenditure, resource value, and death age are somewhat arbitrary and were chosen at reasonable numbers which limited computational load (mostly determined by population size).

## 2.2 | Reproduction

Depending on the scenario, simulations are run either with one or two species in the environment. Throughout the day at a constant rate, pairs of animals belonging to the same species are randomly chosen to mate. The probability of an individual being chosen for mating is proportional to its energy; animals which have acquired more energy thus tend to have better fitness. During each mating, one offspring is produced, to whom each parent donates 1/3 of its energy. In some scenarios, the circadian amplitude and phase of all animals remain fixed at the same values, while in other scenarios, there is variation in the circadian parameters, which are acted on by natural selection to cause evolution. In the case of variable parameters, the child's amplitude is the mean of the parents' ±0.04 SD, and the child's phase angle is the circular mean of the parents' ±0.04*360° SD (Figure 1d). Matings occur at a frequency such that the population of animals would increase by 10% each day if no animals died. In our two-species simulation scenarios, when not testing
competitive exclusion, we stabilized the population levels so that random population fluctuations would not cause one species to disappear by chance. For this, a target population ratio was set (usually 1:1), and the species with a population below target reproduces slightly more, and the species with a higher population reproduces slightly less, proportional to the difference in population. When testing the effect of uneven population sizes on amplitude strength, the population ratio was set to 2:1 and 9:1.

$$r_{\text{adj.}} = r * k * \frac{P_{\text{total}}}{p}$$

For each species, $r_{adj.}$ is the adjusted population growth rate, which depends on the default growth rate $r = 0.1$, $k$ is the target fraction of the population which is composed of this species ($k = 0.5$ for a 1:1 ratio), $p$ is the population of the this species, and $P_{total}$ is the population of all species. The number of progeny in a time step $t_0$ to $t_1$ is thus:

$$N_{\text{progeny}} = p * r_{\text{adj.}} * \frac{t_1 - t_0}{360^{\circ}}$$

These were calculated in 10 time steps per day ($t_1-t_0 = 360^{\circ}$).

## 2.3 | Specialization

Besides circadian activity rhythms, organisms evolve physiological and behavioral traits to specialize for a certain time of day. These adaptations make resource acquisition more energy efficient during their specialized hours, but often at the cost of less efficient activity during nonspecialized hours. We introduced a specialization term to model this, such that the efficiency of gathering resources varies over the day according to a cosine wave and that the peak time of efficiency coincides with peak time of activity (Figure 1e).

$$E_{\text{adj.}} = E_{\text{orig.}} * (1 + S * \cos(t - \varphi))$$

where $E_{adj.}$ is the energy acquired from a resource adjusted by specialization, $E_{orig.} = 0.5$ is the normal energy value of a resource, $S$ is the specialization value, $t$ is the time, and $\varphi$ is the organism's phase angle. Thus, for example when $S = 0.2$, the organism acquires 20% more energy from a resource at its circadian peak and 20% less at its circadian trough.

## 2.4 | Assortative mating

To study speciation, the standard simulation parameters were modified to model assortative mating. We assume that mating must occur in a 9-hr (135°) time window which is fixed relative to the organism's overall phase angle, and the probability of mating between two individuals is proportional to the amount of overlap of their mating

windows, ranging from 0 to 1 (Figure 1f). In other words, individuals which have a phase difference of nine hours or more do not mate, and the probability of mating increases linearly as the phase difference approaches zero. We define a “cross-group mating” metric (CGM) to measure reproductive segregation in the population. To calculate this, we first define an overlap matrix of mating windows, that is, the matrix where each cell $M_{ij}$ is the mating-window overlap between individuals i and j. K-means clustering is used on the overlap matrix to define two clusters of individuals. The CGM is the mean of the overlap values between individuals of different clusters; this is equivalent to the probability that a given mating is between individuals from different clusters. Cross-group mating of 0 thus indicates complete reproductive isolation of two groups.

# 3 | RESULTS

## 3.1 | Methods summary

Detailed explanations of the methods are provided in Section 2 above. Briefly, we simulated communities of one or two species in an environment, into which resources are inputted at a constant rate (Figure 1a). Each individual organism moves around the environment randomly. The amount of movement follows a circadian activity rhythm, which is modeled as a cosine wave and with parameters amplitude and phase angle. Amplitude represents how strong the circadian rhythm is, where high values indicate that the organism has clear high- and low-activity phases, while low values indicate constant moderate activity throughout the day. Phase angle represents the time of day of maximum activity (0° could represent midnight and 180° could represent noon, for example) (Figure 1b). If an organism encounters a resource, it consumes that resource and gains energy, which is in turn used for movement, basal metabolism, and reproduction (Figure 1c). Death occurs if an individual's energy drops to zero or if the individual reaches an age limit. Organisms reproduce sexually, and the probability of an individual reproducing is proportional to its energy. During reproduction, the resulting offspring's amplitude and phase are the means of the parents' traits. When testing the evolution of a trait, the offspring's trait is the parents' mean plus or minus a small random variance (Figure 1d), which naturally results in natural selection and evolution.

## 3.2 | Rhythmic activity prevents competitive exclusion

The competitive exclusion principle asserts that two species cannot occupy the same ecological niche indefinitely; if two species exist in exactly the same niche, then one species will eventually outcompete the other and become the sole occupier of the niche. We modeled this scenario by simulating two species with no circadian rhythms, that is, with amplitude = 0. The two species are constantly active at a moderate level and thus have no differentiation of temporal niches.
In this situation, one species invariably dies out (Figure 2a,b). If both species are rhythmic (amplitude = 1) and are active at the same time of day (same phase angle), there is again no niche differentiation, and one species will eventually die out (Figure 2c,d). However, if both species are rhythmic but are active during opposite times of day, then both species are likely to continue existing (Figure 2e,f). In 2 of the 40 trials of this scenario, one species died out after some time, demonstrating that random fluctuations can still cause a species to disappear. In the remaining 38 trials, both species survived at the end of 2000 days.

## 3.3 | Character displacement of phase angle

Having shown that being active at opposite phases is favorable for the coexistence of two species, we ask whether species can evolve over time to reach such a state of temporal separation. When two similar species occupying similar niches are found in the same geographical area, the differences between them tend to become accentuated in order to distinguish their niches, in a phenomenon known as character displacement (Grant & Grant, 2006). We tested whether two circadian-rhythmic species would display character displacement on a temporal axis. In this scenario, all individuals started with approximately the same phase angle at the beginning of the simulation, and the offsprings' phase angles were allowed to vary and thus evolve, while amplitudes were fixed at 1.

The two species' phase angles of activity did indeed diverge over time in the simulation (Figure 3a). After starting at 0° separation, the mean phase-angle difference between the two species

grew and eventually reached the maximum of 180° and remained so thereafter (Figure 3b). The source of selective pressure is the difference in resource availability between peak and peak-adjacent times. During times of peak activity, available resources are quickly consumed and resource density is thus very low, while in times adjacent to the peak, less resources are consumed, and resource density is higher (Figure 3c). The result of this time-resource gradient is that individuals whose phase angles are at the edges of the population distribution have access to a more resource-dense environment and thus gain more energy (Figure 3d). Since stored energy is proportional to reproduction probability, the individuals at distribution's edges reproduce more, and the species' mean phase angle is shifted.

## 3.4 | Intraspecific competition causes temporal niche expansion by arrhythmicity

Previous work has shown that when interspecific competition is absent, the population's niche width tends to expand, as competition within the species pushes individuals to exploit a greater range of resources (Sjödin, Ripa, & Lundberg, 2018). We tested whether this pattern would be true for the temporal niche. Indeed, when rhythm amplitude was allowed to vary, the amplitude steadily decreased to near-zero levels, representing the transition from a temporal specialist to a generalist strategy. To clarify, a broad circadian niche does not necessarily imply the absence of sleep or sleep-like rest; these organisms could follow a cathemeral pattern in which sleep is scattered throughout the day.

<figure><img src="image_5.png" /></figure>

**FIGURE 2** Circadian rhythms prevent competitive exclusion. (a) When two species with no activity rhythms (amplitude = 0) are placed in the environment, one species will always die out. An example of an individual trial is displayed. (b) Summary of 40 trials. Since the two species' population levels tend to be symmetric, only the species which is excluded is plotted for clarity. The trials were sorted by the day that competitive exclusion is reached, and the median (20th) trial is represented by the dark line, while the 4th, 12th, 28th, and 36th trials are represented by faint lines. (c) When two species have circadian rhythms (amplitude = 1) and are active at the same time of day, one species will always die out, as before. Example individual trial. (d) Summary of 40 trials. (e) When two species have circadian rhythms are active 180° out of phase, the two species can coexist indefinitely. Example individual trial. (f) Summary of multiple trials. Since competitive exclusion was rare in this scenario, five population profiles were selected at random to be plotted
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** Character displacement of circadian phase. (a) Circular histogram displaying the separation of phase angles of activity in two species over time in one trial. The difference in the two species' phase angles quickly diverges to 180°. (b) Summary of 40 trials. Trials were sorted by days until a phase difference of 150° was reached; the median 20th trial is plotted in a dark line, while the 4th-12th, 28th, and 36th trials are plotted in faint lines. (c) When the two species are at similar phase angles of entrainment, there is a large difference in total resources in the environment during different times of day. When the two species are separated in phase, there is no longer a predictable cycle in resource availability. (d) Scatterplot of organisms' stored energy versus their phase angle of entrainment relative to the population mean. The parabolic best fit line is shown in red. Individuals at the edges of the population phase distribution acquire more energy and consequently are more likely to reproduce</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**FIGURE 4** Intraspecific competition selects for arrhythmicity. (a) Even with a positive specialization value, amplitude is driven toward low levels when there is only one species. (b) Amplitude settles on an intermediate value when specialization = 0.4, demonstrating a positive correlation between specialization and resulting amplitude. (c) In the single-species scenario, specialization must be set quite high at 0.75 for amplitude to reach near 1. (d) Standard deviation of phase angle as a function of amplitude. This illustrates an inverse relationship between expanded interindividual niche width (high phase angle variance) and expanded intraindividual niche width (low amplitude)</figcaption></figure>

However, other factors may prevent dampening of rhythms. To adapt to a temporal niche, organisms do not only develop activity rhythms but additional morphological and physiological

specializations as well, such as in vision, temperature regulation, and nonlocomotor circadian rhythms (Hut, Kronfeld-Schor, van der Vinne, & De la Iglesia, 2012; Kronfeld-Schor & Dayan, 2003). These
adaptations allow an animal to acquire more resources per unit energy expenditure during their specialized "on" phase, but often at the cost of being less efficient during their "off" phase. A temporally specialized organism would thus incur costs when adopting a generalist activity rhythm. To model these temporal adaptations, we added a specialization term to the simulation. The degree of specialization over the day follows a cosine wave such that at specialization = 0.1, an animal acquires an additional 10% energy from consuming a resource during its circadian peak and 10% less energy during the trough. For simplicity, we assume that the phase angle of the specialization wave matches exactly that of the cosine wave of circadian activity, as opposed to modeling them as separate phenotypes, which would naturally converge to match phase angles anyway so that organisms are most specialized when they are most active.

We found that specialization must be quite high to overcome
competition's negative selection for amplitude (Figure 4). At low
levels of specialization, the equilibrium amplitude is positively cor-
related with the degree of specialization. Only when specialization
is at least 0.75 does the resulting amplitude approach 1 (Figure 4c).
This is rather intense specialization, indicating that resource acqui-
sition is seven times more efficient during the peak than during the
trough. This suggests that, in an environment where resource avail-
ability is roughly constant throughout the day, intraspecific competi-
tion is a strong factor selecting for arrhythmicity.

When a population is freed from competition, the subsequent
niche expansion could occur through increasing the niche width
within individuals, or through increasing the variance between in-
dividuals who are still relatively specialized. Whereas low rhythm
amplitude represents within-individual niche expansion, we also
observed between-individual niche expansion, represented the
greater range of phase angles when amplitude is high (Figure 4d).
When amplitude is low, there is no selective pressure to promote
interindividual variance, so sexual reproduction pulls the population
toward the mean.

3.5 | Intraspecific competition allows rhythmicity to develop

While intraspecific competition acts to prevent rhythmicity, we asked whether a competitor species would cause rhythmicity to develop. When specialization = 0, organisms did not show a trend toward either high or low rhythmicity (Figure 5a). Since the two species' activity rhythms were in approximately opposite phase and similar amplitude, the summed activity of all individuals in the environment was approximately constant throughout the day, which in turn caused resource availability to be approximately constant. This situation is much like the final condition of Figure 3a,c. With no time-resource gradient, there is no evolutionary pressure for rhythmicity to evolve.

While interspecific competition alone does not positively select for amplitude, it also effectively neutralizes intraspecific competition's negative selection for amplitude. In this state, if the organisms have even a small degree of specialization, it can drive amplitude to increase steadily over time in both species (Figure 5b), with higher specialization values causing amplitude to increase more quickly. The amplitudes of the two species tended to be correlated (Figure 5a,b). The correlation of the two species' amplitudes can be explained by the fact that a chance increase in amplitude in one species causes a resource density difference in time, which encourages the other species to increase their amplitude due to interspecific competition and the first species to decrease their amplitude due to intraspecific competition.

The relative strength of interspecific versus intraspecific com-
petition influences the equilibrium amplitude of each species. A
high-population species is mainly competing against other conspe-
cifics, while a low-population species is mainly competing with allo-
specifics. To simulate asymmetric populations, we defined a target
population ratio, and the populations were stabilized around this
ratio by automatically adjusting their reproduction rate. We saw that
when species A is set to be two times as numerous as species B, spe-
cies B develops high rhythmicity while species A's amplitude hovers

<figure><img src="image_8.png" /></figure>

**FIGURE 5** Interspecific competition with specialization allows development of biological rhythms. (a) With two species, amplitude drifts randomly without pressure toward either 0 or 1. (b) When a specialization value is added to the model, amplitudes rise steadily over time. (c) The final equilibrium amplitudes depend on relative population numbers. When the blue species is twice as numerous as the orange, its final amplitude remains lower than the orange's. (d) When the blue species is nine times as numerous as the orange, blue's amplitude remains low while orange's is high
around 0.7 (Figure 5c). If species A is set to be nine times as numerous as species B, species A's amplitude remains only at 0.2, demonstrating that intraspecific rather than interspecific competition is the much stronger force in determining its equilibrium amplitude. This also demonstrates that a dominant species which has even mild rhythms can strongly determine the rhythms of minor competitors which have smaller populations. In total, these results suggest that the presence of competitors is a strong factor promoting the expression of biological rhythms.

## 3.6 | Population splitting by differentiation of phase

The single-species scenario illustrated that intraspecific competition causes expansions in both within-individual and between-individual niche widths. Could an increase in between-individual variation cause the population to split into two distinct groups? This would in effect cause speciation, or at least contribute strongly to the ecological divergence and reproductive isolation which are precursors to full speciation. To test this, assortative mating was introduced to the model. Individuals only mate with a partner whose phase differs by less than 9 hr (135°), and the probability of choosing a given partner increases linearly as the phase difference decreases to zero. Amplitudes were fixed at 1 for this test.

We found that, with the addition of assortative mating, the population did indeed segregate into two groups (Figure 6a). We defined a cross-group mating value (CGM) to quantify reproductive isolation, where low CGM indicates high reproductive isolation. The CGM tended to stay constant for some time and then drop suddenly (Figure 6b). This suggests that there is a threshold effect, such that once the conditions occur for a secondary population to escape “gene flow” with the main population, full reproductive isolation quickly follows. To separate from the main population, there must be enough individuals at the far edges of the phase distribution to sustain a “breeding population,” such that a tipping point is reached where individuals in the incipient secondary population become more likely to breed within its own population rather than the main population. Once this tipping point is reached, the secondary

population quickly splits off from the main group. Given the random nature of mate selection and progeny's phase variance in our model, the conditions where this happens occurs stochastically, explaining the large variance in speciation time in the simulations.

Population splitting was not guaranteed for all parameter settings. Conditions necessary to achieve speciation in this model are somewhat constrained, which is unsurprising given the high barriers to sympatric speciation in nature (Via, 2001). The time-resource gradient must be quite strong to reach the necessary levels of density-dependent disruptive selection, and thus, population splitting will not occur for low resource regeneration rates or low-population levels. A smaller mating window of course also promotes reproductive isolation. Higher variance in offsprings' phases also promotes reproductive isolation by increasing the chance of a breakaway group to escape “gene flow” with the rest of the population. A long lifespan also favors speciation, since individuals at the edges of the distribution then have more time to build up an energy surplus.

### 3.6.1 | Population splitting versus individual expansion

Intraindividual niche expansion will tend to oppose speciation, since a decrease of rhythm amplitude would weaken the time-resource gradient and broaden mating windows. Conversely, speciation into two populations occupying opposite niches will remove selective pressure for broad individual niches since the time-resource gradient would be erased (analogous to the last case in Figure 3c). Faced with intraspecific competition, which of these two mutually exclusive paths will a population take? Results suggest that it depends heavily on the relative flexibility of the amplitude or phase traits. We set standard deviations of phase angle and amplitude to various test values, with a higher standard deviation resulting in quicker evolution of that trait. When phase standard deviation was high and amplitude standard deviation was low, populations were more likely to speciate; conversely, when phase standard deviation was low and amplitude standard deviation was high, populations tended to become arrhythmic (Table 1).

<figure><img src="image_9.png" /><figcaption>**FIGURE 6** Population splitting by phase. (A) Circular histograms showing distributions of phase angles of entrainment in one trial. On day 340, the population was a single mating group (cross-group mating = 0.67). On day 580, the population is in the process of splitting into two mating groups (CGM = 0.38). By day 630, there are two populations which are almost completely reproductively isolated (CGM = 0.05). (b) Summary of 40 trials. Trials were sorted by days to reach CGM < 0.1; the 20th trial is plotted in dark blue, while the 4th, 12th, 28th, and 36th trials are plotted in light blue</figcaption></figure>
<table><caption>TABLE 1 Number of trials achieving speciation depend on the variances of amplitude versus phase angle</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="4">Amplitude standard deviation</th></tr><tr><th>0.02</th><th>0.03</th><th>0.04</th><th>0.05</th></tr></thead><tbody><tr><td rowspan="4">Phase standard deviation</td><td>0.04</td><td>10/20</td><td>0/20</td><td>0/20</td><td>0/20</td></tr><tr><td>0.045</td><td>20/20</td><td>14/20</td><td>9/20</td><td>5/20</td></tr><tr><td>0.05</td><td>20/20</td><td>20/20</td><td>17/20</td><td>12/20</td></tr><tr><td>0.055</td><td>20/20</td><td>20/20</td><td>20/20</td><td>17/20</td></tr></tbody></table>

## 4 | DISCUSSION

In summary, we show that circadian rhythms can be a way by which organisms differentiate ecological niches and that competition for resources can shape the evolution of rhythms. While two arrhythmic species in an environment will always result in the competitive exclusion of one species, circadian rhythmicity allows them to differentiate their temporal niches and thus coexist indefinitely. Two rhythmic species which are active at the same time of day will evolve to be active at opposite times, demonstrating character displacement of circadian phase. When only one species occupies an environment, intraspecific competition strongly selects for arrhythmicity; however, the presence of another competing species allows circadian rhythms to develop much more easily, when combined with other temporal specializations. Finally, if individuals preferentially mate with others who have similar times of activity, then intraindividual competition can cause a population to speciate into two populations separated in time.

This is the first modeling study to how selective pressure from other organisms favors the evolution of biological rhythms. While other modeling studies have examined the “internal” forces which make rhythms and sleep physiologically efficient (Hellwegera, 2010; Schmidt, 2014; Schmidt, Swang, Hamilton, & Best, 2017), none have looked at interindividual interactions.

### 4.1 | Timescales of analysis

Though we framed our model in terms of organisms evolving genotypically over time, niche partitioning can equally occur at the timescale of individual behavior or ecological community (Mrosovsky & Hattar, 2005), and the conclusions of this study can easily be generalized to subevolutionary timescales and mechanisms. Many animals are known to switch their temporal niches depending on various environmental factors (Hut et al., 2012); for example, the aforementioned studies on minks, spiny mice (Gutman & Dayan, 2005), and foxes (Di Bitetti et al., 2009) demonstrate temporal partitioning at the subevolutionary level. Indeed, every clear observation of competition's effect on temporal niches has occurred at the behavioral/ecological level, and inferring the role of competition in a species' evolution is difficult if not impossible, since those historical interspecies interactions are unobservable. However, if temporal

displacement continues for an extended time, it seems likely that organisms would eventually evolve to match their internal clock with their expressed behavior. Research on the health consequences of misalignment between the circadian clock and activity rhythms (Baron & Reid, 2014; Ouyang et al., 1998; Rüger & Scheer, 2009) suggests clear selective mechanisms by which behavioral displacements can translate into genetic/evolutionary changes.

The simulation was framed in terms of circadian rhythms, but the results could equally apply to other endogenous biological rhythms like seasonal (circannual) rhythms. A major difference between the daily and seasonal timescales is that individuals can adjust their daily behavior to a new time through nongenetic mechanisms, whereas this is probably less common at longer timescales, though it certainly depends heavily on the particular mechanisms governing these rhythms.

### 4.2 | Competitive exclusion and character displacement

In our model, we demonstrated that competitive exclusion is unavoidable if two species share an ecological niche. Two species which are temporally separated are much more likely to coexist for extended periods of time. Additionally, we show that two circadian-rhythmic species which initially are active at the same phase evolve to be active at opposite times of day in order to avoid competition for resources. The most likely direct natural analog of this scenario is the introduction of a new species to a region already occupied by another similar species. This could happen, for example, through the breakdown of geographical barriers, chance migration to a new area (Lomolino, 2000), introduction by humans, or climate change-induced range shifts (Thomas, 2010). The invasive American mink in the United Kingdom offers a good real-life example: Originally nocturnal when first introduced, the mink largely switched to diurnality when the native otters and polecats recovered their population levels (Harrington et al., 2009). The previously mentioned studies in Argentinian foxes (Di Bitetti et al., 2009) and desert spiny mice (Kronfeld-Schor & Dayan, 2003) also demonstrate niche switching.

### 4.3 | Intraspecific competition and niche expansion

Intraspecific competition is hypothesized to increase the niche width in a population, especially in situations where a population is released from competition from other species, which could occur after migration or introduction to a new region, for instance (Lister, 1976; Bolnick, 2001). Our single-species simulations demonstrate this phenomenon for temporal niches.

Niche expansion can happen either by increasing the variance of resource type usage between individuals (Svanbäck & Bolnick, 2006; Van Valen, 1965) or by increasing the range of resource usage of each individual (Bolnick et al., 2010). Interindividual expansion appears to be the more common pattern in nature in general (Bolnick, Svanbäck,
Araújo, & Persson, 2007), especially if intraindividual generalization has associated costs (Sjödin et al., 2018). Our single-species simulation displayed both interindividual variation, visible as a broadened distribution of phase-angle, and intraindividual generalization, visible as flattened amplitude. The model's design does not allow for very extreme interindividual expansion, since activity timing was determined strictly by genetics and not by individuals' behavioral adjustments, and sexual reproduction tended to pull genotypes toward the population mean. Programming a nongenetic behavioral component may be interesting for future studies.

A few studies in fish have observed offset patterns of diel feeding activity between individuals, with larger dominant individuals feed at preferred times while subordinates being pushed to other times (Alanärä, Burns, & Metcalfe, 2001; Reebs, 2002). This suggests that preexisting heterogeneities in the population such as social dominance could predispose niche expansion to occur between rather than within individuals. We are aware of only one study to observe temporal niche expansion in the wild; American minks which were introduced to a remote island with no competitors were observed to have expanded temporal niches compared to their mainland counterparts (Crego, Jiménez, & Rozzi, 2018); it is unknown whether this was due more to inter- or intraindividual expansion. Intraspecific competition is also likely to play a role in instances of temporal niche expansion after release from predation (McCaulley et al., 2012). Predator-free space may be seen as a type of resource, which becomes more abundant after removal of the predator. Though never specifically examined, it is reasonable to suspect that intraspecific competition could then drive the subsequent niche expansion observed after predator release.

Arrhythmic/cathemeral activity patterns appear to be somewhat common in nature (Roll, Dayan, & Kronfeld-Schor, 2006; van Schaik & Griffiths, 1996). However, besides the aforementioned observation of invasive minks (Crego et al., 2018), there have not been any other experiments or observations to directly test competition's role in causing temporal niche expansion. We suggest that this phenomenon is not uncommon and may explain the behavior and evolution of a number of these arrhythmic animals.

## 4.4 | Competition facilitates emergence of rhythms

The presence of an interspecific competitor erases the advantage that arrhythmicity had in a single-species scenario. A small amount of time-of-day specialization (theoretically, any amount of specialization) can then cause organisms to evolve high amplitudes. We only modeled exploitative competition here, where fitness differences were only driven by differential access to resources. If interference competition were also present, for example if a cost was incurred every time an individual encountered an allospecific, then development of rhythm amplitude would most likely occur even more quickly.

As previously mentioned, the comparison between invasive minks on an island and mainland minks is an instance of

apparent competition-mediated niche expansion/contraction (Crego et al., 2018), but we are not aware of any studies which has observed competition inducing increased amplitude. Again, we believe this phenomenon does exist broadly, but has simply not been observed yet. Another potential real-life analog is if two arrhythmic species existing in a constant environment—cave-dwelling fishes for example (Kavanau, 1998)—are suddenly exposed to a cycling environment and are thus pushed to temporally differentiate.

## 4.5 | Speciation

When individuals were limited to mating with other individuals who were active at similar times of day, intraspecific competition could cause a single population to separate into two groups which were largely reproductively isolated, in effect causing speciation. This is specifically a case of *sympatric* speciation, in which new species arise without geographic separation. Speciation required high levels of density-dependent selection (i.e., a strong time-resource gradient) to occur. This is not surprising given that sympatric speciation in general is believed to be quite difficult to achieve in nature, such that its existence at all was questioned before recent decades (Via, 2001). A theoretical problem against sympatric speciation is that the phenotype which improves fitness usually has no effect on mate choice, and thus even when disruptive selection exists, there is no mechanism for it to translate to mating preferences. However, speciation by time (allochronic speciation) has an advantage here because it is a “magic” trait (Servedio, Doorn, Kopp, Frame, & Nosil, 2011), that is, it has fitness effects which are acted on by natural selection, and it is also directly responsible for assortative mating.

A recent review (Taylor & Friesen, 2017) finds good evidence for dozens of allochronic speciation events, with most cases occurring at the seasonal rather than the circadian time span, but although intraspecific competition is known to play a role in promoting speciation in general (Bolnick, 2004; Laland & Sterelny, 2006), it is unclear how important intraspecific competition was in these specific cases of allochronic speciation. In a study of Madeiran storm petrels, the authors hypothesize that competition for food and nest space may have led to allochronic separation into hot- and cool-season breeding populations (Monteiro & Furness, 1998). Besides competition, allochronic speciation can occur by other reasons too, such as temporal isolation (Hendry & Day, 2005), environmental changes (Yamamoto & Sota, 2012), or a founder effect (Santos et al., 2007); temporal separation can also be secondary to another adaptation (Ragland, Sim, Goudarzi, Feder, & Hahn, 2012) or serve to reinforce reproductive isolation. Intraspecific competition can act alongside these other mechanisms to produce allochronic speciation.

## 4.6 | Sleep

The question of why sleep exists is one of the big questions of biology. Traditionally framed something of a paradox that an animal
would abandon motor and sensory activity for several hours per day, it is perhaps not so mysterious after all when considering that temporal niches exist. External pressures, whether they be competition or other environmental forces, push animals to specialize for a time of day. Internal pressures, such as energy conservation (Siegel, 2009), the need to perform restorative functions, and the separation or coordination of clashing or synergistic physiological processes (Baron & Reid, 2014; Schmidt, 2014; Schmidt et al., 2017), also encourage the evolution of sleep, or at least a rest period. Especially with external and internal pressures combined, sleep/rest appears to be an unsurprising evolutionary consequence. Stated colloquially, "if I can't be active right now, then I might as well rest. If rest is efficient, then I might as well do it when I can't be active."

As to why a defined sleep state exists rather than simply circadian-mediated inactivity, we propose that true sleep can be viewed as an adaptation to delay rest. Urgent circumstances sometimes interrupt an animal's regularly scheduled rest, and thus, the defining feature of true sleep, sleep homeostasis ("sleepiness"), encourages animals to later pay back the lost rest. Another feature of true sleep is that it is a relatively distinct mode with a well-defined sleep-wake boundary. In unpredictable environments, such a clear distinction allows animals to maximize rest when the opportunity is available, rather than a slow transition between rest and sleep.

## 4.7 | Methodology and future directions

We intentionally created a model without distinct day and night phases, so that it would be more generalizable to a variety of hypothetical conditions. Nevertheless, for most organisms, the changes between day and night represent a major division in the environment (van Schaik & Griffiths, 1996), and thus, alternative models may explicitly include this. Activity patterns other than cosine waves, or testing the variable of daily total activity, could also be interesting.

We decided not to use a genetic model of circadian rhythms, opting instead to model only phenotypes with a constant distribution. We do not believe this is a weakness because it is more or less equivalent to a situation of polygenic inheritance with many genes. However, a genetic model may have some advantages such as allowing a clearer look at population genetics. Additional models have been developed to explain the theoretical benefits of sleep and circadian rhythms from the point of view of internal physiological organization (Schmidt, 2014; Schmidt et al., 2017), and it could be interesting to integrate models of these "internal" benefits of rhythms with the "external" ecological benefits modeled here.

Many questions regarding this topic remain open for investigation. Among the many possible axes by which competing species can differentiate, what promotes partitioning of time over partitioning of other resources? How does competition interact with other forces, like predation, to influence rhythms? To what extent does competition shape activity rhythms in nature? In addition to additional field studies, even well-designed experiments with lab rodents could help

us understand the circumstances in which activity rhythms are modified by competition.

## ACKNOWLEDGMENTS

Vance Gao was supported by a NIH T32 grant (HL 007909) and the I-SITE ULNE ERC Generator grant.

## CONFLICT OF INTEREST

The authors report no conflicts of interest.

## AUTHOR CONTRIBUTION

**Vance Difan Gao:** Conceptualization (lead); Data curation (lead); Investigation (lead); Methodology (lead); Project administration (equal); Software (lead); Visualization (lead); Writing-original draft (lead). **Sara Morley-Fletcher:** Formal analysis (equal); Funding acquisition (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing-review & editing (equal). **Stefania Maccari:** Formal analysis (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing-review & editing (equal). **Martha Hotz Vitaterna:** Conceptualization (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing-review & editing (equal). **Fred W. Turek:** Conceptualization (equal); Funding acquisition (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing-review & editing (equal).

## OPEN RESEARCH BADGES

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

This article has been awarded Open Materials, Open Data Badges. All materials and data are publicly accessible via the Open Science Framework at https://datadryad.org/stash/dataset/doi:10.5061/dryad.gf1vhhmmj.

## DATA AVAILABILITY STATEMENT

The Python code and all relevant data for this project may be found at https://datadryad.org/stash/dataset/doi:10.5061/dryad.gf1vhmmj (https://doi.org/10.5061/dryad.gf1vhhmmj) or by request to the author.

## ORCID

Vance Difan Gao ID https://orcid.org/0000-0003-4492-1154

Sara Morley-Fletcher ID https://orcid.org/0000-0003-0423-5563

Stefania Maccari ID https://orcid.org/0000-0002-8336-5990

Martha Hotz Vitaterna ID https://orcid.org/0000-0001-6243-2093

## REFERENCES

Adams, R. A., & Thibault, K. M. (2006). Temporal resource partitioning by bats at water holes. *Journal of Zoology*, 270(3), 466-472. https://doi.org/10.1111/j.1469-7998.2006.00152.x

Alanärä, A., Burns, M. D., & Metcalfe, N. B. (2001). Intraspecific resource partitioning in brown trout: The temporal distribution of foraging is determined by social rank. *Journal of Animal Ecology*, 70(6), 980-986. https://doi.org/10.1046/j.0021-8790.2001.00550.x
Amarasekare, P. (2002). Interference competition and species coexistence. *Proceedings of the Royal Society of London. Series B: Biological Sciences*, 269(1509), 2541-2550. https://doi.org/10.1098/rspb.2002.2181

Bakker, E. S., Reiffers, R. C., Olff, H., & Gleichman, J. M. (2005). Experimental manipulation of predation risk and food quality: Effect on grazing behaviour in a central-place foraging herbi-vore. *Oecologia*, 146(1), 157-167. https://doi.org/10.1007/s0044-2-005-0180-7

Baron, K. G., & Reid, K. J. (2014). Circadian misalignment and health. *International Review of Psychiatry*, 26(2), 139-154. https://doi.org/10.3109/09540261.2014.911149

Begon, M., Townsend, C. R., & Harper, J. L. (2006). *Ecology: from individuals to ecosystems*. Malden, MA: Blackwell Publishing.

Bolnick, D. I. (2001). Intraspecific competition favours niche width expansion in *Drosophila melanogaster*. *Nature*, 410(6827), 463-466. https://doi.org/10.1038/35068555

Bolnick, D. I. (2004). Can intraspecific competition drive disruptive selection? An experimental test in natural populations of sticklebacks. *Evolution*, 58(3), 608-618. https://doi.org/10.1111/j.0014-3820.2004.tb01683.x

Bolnick, D. I., Ingram, T., Stutz, W. E., Snowberg, L. K., Lau, O. L., & Pauli, J. S. (2010). Ecological release from interspecific competition leads to decoupled changes in population and individual niche width. *Proceedings of the Royal Society B: Biological Sciences*, 277(1689), 1789-1797. https://doi.org/10.1098/rspb.2010.0018

Bolnick, D. I., Svanbäck, R., Araújo, M. S., & Persson, L. (2007). Comparative support for the niche variation hypothesis that more generalized populations also are more heterogeneous. *Proceedings of the National Academy of Sciences of the United States of America*, 104(24), 10075-10079. https://doi.org/10.1073/pnas.0703743104

Caro, T. M., & Stoner, C. J. (2003). The potential for interspecific competition among African carnivores. *Biological Conservation*, 110(1), 67-75. https://doi.org/10.1016/S0006-3207(02)00177-5

Castillo-Ruiz, A., Paul, M. J., & Schwartz, W. J. (2012). In search of a temporal niche: Social interactions. In *Progress in brain research* (vol. 199, pp. 267-280). https://doi.org/10.1016/B978-0-444-59427-3.00016-2

Crego, R. D., Jiménez, J. E., & Rozzi, R. (2018). Potential niche expansion of the American mink invading a remote island free of native-predatory mammals. *PLoS One*, 13(4), 1-18. https://doi.org/10.1371/journal.pone.0194745

Di Bitetti, M. S., Di Blanco, Y. E., Pereira, J. A., Paviolo, A., & Pírez, I. J. (2009). Time partitioning favors the coexistence of sympatric crab-eating foxes (*Cerdocyon thous*) and pampas foxes (*Lycalopex gymnocercus*). *Journal of Mammalogy*, 90(2), 479-490. https://doi.org/10.1644/08-mamm-a-113.1

Gibson, R. N. (1967). Experiments on the tidal rhythm of *Blennius pholis*. *Journal of the Marine Biological Association of the United Kingdom*, 47(1), 97-111.

Grant, P. R., & Grant, B. R. (2006). Evolution of character displacement in Darwin's finches. *Science*, 313(5784), 224-226. https://doi.org/10.1126/science.1128374

Gutman, R., & Dayan, T. (2005). Temporal partitioning: An experiment with two species of spiny mice. *Ecology*, 86(1), 164-173. https://doi.org/10.1890/03-0369

Harrington, L. A., Harrington, A. L., Yamaguchi, N., Thom, M. D., Ferreras, P., Windham, T. R., & Macdonald, D. W. (2009). The impact of native competitors on an alien invasive: Temporal niche shifts to avoid interspecific aggression? *Ecology*, 90(5), 1207-1216. https://doi.org/10.1890/08-0302.1

Hayward, M. W., & Slotow, R. (2009). Temporal partitioning of activity in large African carnivores: Tests of multiple hypotheses. *African Journal of Wildlife Research*, 39(2), 109-125. https://doi.org/10.3957/056.039.0207

Heesy, C. P., & Hall, M. I. (2010). The nocturnal bottleneck and the evolution of mammalian vision. *Brain, Behavior and Evolution*, 75(3), 195-203. https://doi.org/10.1159/000314278

Hellwegera, F. L. (2010). Resonating circadian clocks enhance fitness in cyanobacteria in silico. *Ecological Modelling*, 221(12), 1620-1629. https://doi.org/10.1016/j.ecolmodel.2010.03.015

Hendry, A. P., & Day, T. (2005). Population structure attributable to reproductive time: Isolation by time and adaptation by time. *Molecular Ecology*, 14(4), 901-916. https://doi.org/10.1111/j.1365-294X.2005.02480.x

Hut, R. A., Kronfeld-Schor, N., van der Vinne, V., & De la Iglesia, H. (2012). *In search of a temporal niche: Environmental factors. Progress in Brain Research* (1st ed., vol. 199). Amsterdam, The Netherlands: Elsevier B.V.

Kamenova, S., Tougeron, K., Cateine, M., Marie, A., & Plantegenest, M. (2015). Behaviour-driven micro-scale niche differentiation in carbid beetles. *Entomologia Experimentalis et Applicata*, 155(1), 39-46. https://doi.org/10.1111/eea.12280

Kavanau, J. L. (1998). Vertebrates that never sleep: Implications for sleep's basic function. *Brain Research Bulletin*, 46(4), 269-279. https://doi.org/10.1016/S0361-9230(98)00018-5

Kronfeld-Schor, N., & Dayan, T. (2003). Partitioning of time as an ecological resource. *Annual Review of Ecology, Evolution, and Systematics*, 34(1), 153-181. https://doi.org/10.1146/annurev.ecolsys.34.011180.132435

Kunz, T. H. (1973). Resource utilization: Temporal and spatial components of bat activity in central Iowa. *Journal of Mammalogy*, 54(1), 14-32. https://doi.org/10.2307/1378869

Laland, K. N., & Sterelnky, K. (2006). The conditions for speciation through intraspecific competition. *Evolution. International Journal of Organic Evolution*, 60(9), 1751-1762. https://doi.org/10.1554/05-570.1

Lister, B. C. (1976). The nature of niche expansion in West Indian Anolis lizards I: Ecological consequences of reduced competition. *Evolution*, 659-676.

Lomolino, M. V. (2000). A call for a new paradigm of island biogeography. *Global Ecology and Biogeography*, 9(1), 1-6. https://doi.org/10.1046/j.1365-2699.2000.00185.x

Maurer, B. A. (1984). Interference and exploitation in bird communities. *The Wilson Bulletin*, 380-395.

McCauley, D. J., Hoffmann, E., Young, H. S., & Micheli, F. (2012). Night shift: Expansion of temporal niche use following reductions in predator density. *PLoS One*, 7(6), e38871. https://doi.org/10.1371/journal.pone.0038871

Monteiro, L. R., & Furness, R. W. (1998). Speciation through temporal segregation of Madeiran storm petrel (Oceanodroma castro) populations in the Azores? *Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences*, 353(1371), 945-953.

Mrosovsky, N., & Hattar, S. (2005). Diurnal mice (Mus musculus) and other examples of temporal niche switching. *Journal of Comparative Physiology A: Neuroethology, Sensory, Neural, and Behavioral Physiology*, 191(11), 1011-1024. https://doi.org/10.1007/s0035-9-005-0017-1

Ouyang, Y., Andersson, C. R., Kondo, T., Golden, S. S., & Johnson, C. H. (1998). Resonating circadian clocks enhance fitness in cyanobacteria. *Proceedings of the National Academy of Sciences of the United States of America*, 95(15), 8660-8664. https://doi.org/10.1016/j.ecolmodel.2010.03.015

Ragland, G. J., Sim, S. B., Goudarzi, S., Feder, J. L., & Hahn, D. A. (2012). Environmental interactions during host race formation: Host fruit environment moderates a seasonal shift in phenology in host races of *Rhagoletis pomonella*. *Functional Ecology*, 26(4), 921-931. https://doi.org/10.1111/j.1365-2435.2012.01992.x

Reebs, S. G. (2002). Plasticity of diel and circadian activity rhythms in fishes. *Reviews in Fish Biology and Fisheries*, 12(4), 349-371. https://doi.org/10.1023/A:1025371804611
Roll, U., Dayan, T., & Kronfeld-Schor, N. (2006). On the role of phylogeny in determining activity patterns of rodents. *Evolutionary Ecology*, 20(5), 479-490. https://doi.org/10.1007/s10682-006-0015-y

Rüger, M., & Scheer, F. A. J. L. (2009). Effects of Circadian disruption on the cardiometabolic system. *Reviews in Endocrine and Metabolic Disorders*, 10(4), 245-260. https://doi.org/10.1007/s11154-009-9122-8

Santos, H., Rousselet, J., Magnoux, E., Paiva, M.-R., Branco, M., & Kerdelhué, C. (2007). Genetic isolation through time: Allochronic differentiation of a phenologically atypical population of the pine processionary moth. *Proceedings of the Royal Society B: Biological Sciences*, 274(1612), 935-941.

Schmidt, M. H. (2014). The energy allocation function of sleep: A unifying theory of sleep, torpor, and continuous wakefulness. *Neuroscience and Biobehavioral Reviews*, 47, 122-153. https://doi.org/10.1016/j.neubiorev.2014.08.001

Schmidt, M. H., Swang, T. W., Hamilton, I. M., & Best, J. A. (2017). State-dependent metabolic partitioning and energy conservation: A theoretical framework for understanding the function of sleep. *PLoS One*, 12(10), 1-15. https://doi.org/10.1371/journal.pone.0185746

Servedio, M. R., Van Doorn, G. S., Kopp, M., Frame, A. M., & Nosil, P. (2011). Magic traits in speciation: 'magic' but not rare? *Trends in Ecology and Evolution*, 26(8), 389-397. https://doi.org/10.1016/j.tree.2011.04.005

Siegel, J. M. (2009). Sleep viewed as a state of adaptive inactivity. *Nature Reviews Neuroscience*, 10(10), 747-753. https://doi.org/10.1038/nrn2697

Sjödin, H., Ripa, J., & Lundberg, P. (2018). Principles of niche expansion. *Proceedings of the Royal Society B: Biological Sciences*, 285(1893), 20182603. https://doi.org/10.1098/rspb.2018.2603

Svanbäck, R., & Bolnick, D. I. (2006). Intraspecific competition drives increased resource use diversity within a natural population. *Proceedings of the Royal Society B: Biological Sciences*, 274(1611), 839-844.

Taylor, R. S., & Friesen, V. L. (2017). The role of allochrony in speciation. *Molecular Ecology*, 26(13), 3330-3342. https://doi.org/10.1111/mec.14126

Thomas, C. D. (2010). Climate, climate change and range boundaries. *Diversity and Distributions*, 16(3), 488-495. https://doi.org/10.1111/j.1472-4642.2010.00642.x

Tilman, D. (2004). Niche tradeoffs, neutrality, and community structure: A stochastic theory of resource competition, invasion, and community assembly. *Proceedings of the National Academy of Sciences of the United States of America*, 101(30), 10854-10861. https://doi.org/10.1073/pnas.0403458101

Valeix, M., Chamaillé-Jammes, S., & Fritz, H. (2007). Interference competition and temporal niche shifts: Elephants and herbivore communities at waterholes. *Oecologia*, 153(3), 739-748. https://doi.org/10.1007/s00442-007-0764-5

van Schaik, C. P., & Griffiths, M. (1996). Activity periods of Indonesian rain forest mammals. *Biotropica*, 28(1), 105-112. https://doi.org/10.2307/2388775

Van Valen, L. (1965). Morphological variation and width of ecological niche. *The American Naturalist*, 99(908), 377-390. https://doi.org/10.1086/282379

Via, S. (2001). Sympatric speciation in animals: The ugly duckling grows up. *Trends in Ecology and Evolution*, 16(7), 381-390. https://doi.org/10.1002/eqe.4290130210

Yamamoto, S., & Sota, T. (2012). Parallel allochronic divergence in a winter moth due to disruption of reproductive period by winter harshness. *Molecular Ecology*, 21(1), 174-183. https://doi.org/10.1111/j.1365-294X.2011.05371.x

Yerushalmi, S., & Green, R. M. (2009). Evidence for the adaptive significance of circadian rhythms. *Ecology Letters*, 12(9), 970-981. https://doi.org/10.1111/j.1461-0248.2009.01343.x

Ziv, Y., Abramsky, Z., Kotler, B. P., & Subach, A. (1993). Interference competition and temporal and habitat partitioning in two gerbil species. *Oikos*, 237-246. https://doi.org/10.2307/3544810

**How to cite this article:** Gao VD, Morley-Fletcher S, Maccari S, Vitaterna MH, Turek FW. Resource competition shapes biological rhythms and promotes temporal niche differentiation in a community simulation. *Ecol Evol.* 2020;10:11322-11334. https://doi.org/10.1002/ece3.6770