Effects of temperature heterogeneity on freshwater turtle habitat selection at their northern range limit.

Environmental temperature is a crucial resource for ectotherms, affecting their physiology, behaviour and fitness. To maintain body temperatures within a suitable performance range, ectotherms select thermally-favourable locations, but this selection may be challenging in environments with high spatio-temporal heterogeneity. We assessed thermal habitat selection in two freshwater turtles (Emydoidea blandingii; Chrysemys picta) within a thermally heterogeneous environment at two spatial scales (selection of home ranges within the landscape, selection of locations within home ranges) and across seasons, by comparing temperatures at turtle locations vs. those available in the environment. Turtles selected warmer locations compared to those available in aquatic and terrestrial habitats only within home ranges, but did not show any temperature preferences when selecting home ranges at the larger scale. Turtles selected locations that were less thermally-variable than their surroundings, both at the home range scale and within home ranges. Thermal habitat selection was strongest during colder and more thermally-variable pre-nesting season compared to later periods. Despite differences in thermal mass between species, both species responded similarly to temperature variation. We conclude that freshwater turtles at their northern range margin select suitable microclimates within the suite of conditions that are naturally available.

Diversity and composition of mixed-ploidy unisexual salamander assemblages reflect the key influence of host species.

Understanding processes that govern and sustain biological diversity is a central goal of community ecology. Unisexual complexes, where reproduction depends on sperm from males of one or more bisexual host species, are rare and the processes driving their diversity and structure remain poorly understood. Unisexual Ambystoma salamanders produce distinct biotypes ('genomotypes') depending on which bisexual species they 'steal' sperm from. This reproductive mode should generate distinct assemblages depending on the locally available bisexual host species. Yet, how availability and relative abundance of multiple bisexual hosts influences composition and diversity of natural unisexual assemblages at local or regional scales remains unknown. We hypothesize that host identity most directly drives local assemblage composition, with host variation associated with increased beta and gamma diversity within unisexuals. We collected genetic samples from Ambystoma salamanders across Pelee Island, Ontario, Canada (2015-2022). Two host species were identified (A. texanum and A. laterale) with nine sites having a single host and one site having both. Unisexual assemblages were grouped into four clusters by similarity, with host identity being a key determinant. Gamma diversity increased as a result of distinct host-specific assemblages forming at different sites on the island (i.e., high beta diversity). Assemblage composition, but not diversity, was correlated with relative host abundance, which may reflect matching niche requirements between host and unisexual forms they produce. Our results demonstrate that diversity and structure of unisexual assemblages are clearly shaped by their host(s) and such systems may serve as models for studying how biotic interactions shape ecological communities.

When death comes: linking predator-prey activity patterns to timing of mortality to understand predation risk.

The assumption that activity and foraging are risky for prey underlies many predator-prey theories and has led to the use of predator-prey activity overlap as a proxy of predation risk. However, the simultaneous measures of prey and predator activity along with timing of predation required to test this assumption have not been available. Here, we used accelerometry data on snowshoe hares (Lepus americanus) and Canada lynx (Lynx canadensis) to determine activity patterns of prey and predators and match these to precise timing of predation. Surprisingly we found that lynx kills of hares were as likely to occur during the day when hares were inactive as at night when hares were active. We also found that activity rates of hares were not related to the chance of predation at daily and weekly scales, whereas lynx activity rates positively affected the diel pattern of lynx predation on hares and their weekly kill rates of hares. Our findings suggest that predator-prey diel activity overlap may not always be a good proxy of predation risk, and highlight a need for examining the link between predation and spatio-temporal behaviour of predator and prey to improve our understanding of how predator-prey behavioural interactions drive predation risk.

Does coat colour influence survival? A test in a cyclic population of snowshoe hares.

Some mammal species inhabiting high-latitude biomes have evolved a seasonal moulting pattern that improves camouflage via white coats in winter and brown coats in summer. In many high-latitude and high-altitude areas, the duration and depth of snow cover has been substantially reduced in the last five decades. This reduction in depth and duration of snow cover may create a mismatch between coat colour and colour of the background environment, and potentially reduce the survival rate of species that depend on crypsis. We used long-term (1977-2020) field data and capture-mark-recapture models to test the hypothesis that whiteness of the coat influences winter apparent survival in a cyclic population of snowshoe hares (Lepus americanus) at Kluane, Yukon, Canada. Whiteness of the snowshoe hare coat in autumn declined during this study, and snowshoe hares with a greater proportion of whiteness in their coats in autumn survived better during winter. However, whiteness of the coat in spring did not affect subsequent summer survival. These results are consistent with the hypothesis that the timing of coat colour change in autumn can reduce overwinter survival. Because declines in cyclic snowshoe hare populations are strongly affected by low winter survival, the timing of coat colour change may adversely affect snowshoe hare population dynamics as climate change continues.

Developing a classification system to assign activity states to two species of freshwater turtles.

Research in ecology often requires robust assessment of animal behaviour, but classifying behavioural patterns in free-ranging animals and in natural environments can be especially challenging. New miniaturised bio-logging devices such as accelerometers are increasingly available to record animal behaviour remotely, and thereby address the gap in knowledge related to behaviour of free-ranging animals. However, validation of these data is rarely conducted and classification model transferability across closely-related species is often not tested. Here, we validated accelerometer and water sensor data to classify activity states in two free-ranging freshwater turtle species (Blanding's turtle, Emydoidea blandingii, and Painted turtle, Chrysemys picta). First, using only accelerometer data, we developed a decision tree to separate motion from motionless states, and second, we included water sensor data to classify the animal as being motionless or in-motion on land or in water. We found that accelerometers separated in-motion from motionless behaviour with > 83% accuracy, whereas models also including water sensor data predicted states in terrestrial and aquatic locations with > 77% accuracy. Despite differences in values separating activity states between the two species, we found high model transferability allowing cross-species application of classification models. Note that reducing sampling frequency did not affect predictive accuracy of our models up to a sampling frequency of 0.0625 Hz. We conclude that the use of accelerometers in animal research is promising, but requires prior data validation and development of robust classification models, and whenever possible cross-species assessment should be conducted to establish model generalisability.

Balancing food acquisition and predation risk drives demographic changes in snowshoe hare population cycles.

Snowshoe hare cycles are one of the most prominent phenomena in ecology. Experimental studies point to predation as the dominant driving factor, but previous experiments combining food supplementation and predator removal produced unexplained multiplicative effects on density. We examined the potential interactive effects of food limitation and predation in causing hare cycles using an individual-based food-supplementation experiment over-winter across three cycle phases that naturally varied in predation risk. Supplementation doubled over-winter survival with the largest effects occurring in the late increase phase. Although the proximate cause of mortality was predation, supplemented hares significantly decreased foraging time and selected for conifer habitat, potentially reducing their predation risk. Supplemented hares also lost less body mass which resulted in the production of larger leverets. Our results establish a mechanistic link between how foraging time, mass loss and predation risk affect survival and reproduction, potentially driving demographic changes associated with hare cycles.

Food availability and long-term predation risk interactively affect antipredator response.

Food availability and temporal variation in predation risk are both important determinants of the magnitude of antipredator responses, but their effects have rarely been examined simultaneously, particularly in wild prey. Here, we determine how food availability and long-term predation risk affect antipredator responses to acute predation risk by monitoring the foraging response of free-ranging snowshoe hares (Lepus americanus) to an encounter with a Canada lynx (Lynx canadensis) in Yukon, Canada, over four winters (2015-2016 to 2018-2019). We examined how this response was influenced by natural variation in long-term predation risk (2-month mortality rate of hares) while providing some individuals with supplemental food. On average, snowshoe hares reduced foraging time up to 10 h after coming into close proximity (≤75 m) with lynx, and reduced foraging time an average of 15.28 ± 7.08 min per lynx encounter. Hares tended to respond more strongly when the distance to lynx was shorter. More importantly, the magnitude of hares' antipredator response to a lynx encounter was affected by the interaction between food-supplementation and long-term predation risk. Food-supplemented hares reduced foraging time more than control hares after a lynx encounter under low long-term risk, but decreased the magnitude of the response as long-term risk increased. In contrast, control hares increased the magnitude of their response as long-term risk increased. Our findings show that food availability and long-term predation risk interactively drive the magnitude of reactive antipredator response to acute predation risk. Determining the factors driving the magnitude of antipredator responses would contribute to a better understanding of the indirect effects of predators on prey populations.

Increasing availability of palatable prey induces predator-dependence and increases predation on unpalatable prey.

Understanding the factors governing predation remains a top priority in ecology. Using a dragonfly nymph-tadpole system, we experimentally varied predator density, prey density, and prey species ratio to investigate: (i) whether predator interference varies between prey types that differ in palatability, (ii) whether adding alternate prey influences the magnitude of predator interference, and (iii) whether patterns of prey selection vary according to the predictions of optimal diet theory. In single-prey foraging trials, predation of palatable leopard frog tadpoles was limited by prey availability and predator interference, whereas predation of unpalatable toad tadpoles was limited by handling time. Adding unpalatable prey did not affect the predator's kill rate of palatable prey, but the presence of palatable prey increased the influence of predator density on the kill rate of unpalatable prey and reduced unpalatable prey handling time. Prey selection did not change with shifts in the relative abundance of prey types. Instead, predators selected easy-to-capture unpalatable prey at low total densities and harder-to-capture palatable prey at high densities. These results improve our understanding of generalist predation in communities with mobile prey, and illustrate that characteristics of the prey types involved govern the extent to which alternate prey influence the predator's kill rate.

Spatial and environmental influences on selection in a clock gene coding trinucleotide repeat in Canada lynx (Lynx canadensis).

Clock genes exhibit substantial control over gene expression and ultimately life-histories using external cues such as photoperiod, and are thus likely to be critical for adaptation to shifting seasonal conditions and novel environments as species redistribute their ranges under climate change. Coding trinucleotide repeats (cTNRs) are found within several clock genes, and may be interesting targets of selection due to their containment within exonic regions and elevated mutation rates. Here, we conduct inter-specific characterization of the NR1D1 cTNR between Canada lynx and bobcat, and intra-specific spatial and environmental association analyses of neutral microsatellites and our functional cTNR marker, to investigate the role of selection on this locus in Canada lynx. We report signatures of divergent selection between lynx and bobcat, with the potential for hybrid-mediated gene flow in the area of range overlap. We also provide evidence that this locus is under selection across Canada lynx in eastern Canada, with both spatial and environmental variables significantly contributing to the explained variation, after controlling for neutral population structure. These results suggest that cTNRs may play an important role in the generation of functional diversity within some mammal species, and allow for contemporary rates of adaptation in wild populations in response to environmental change. We encourage continued investment into the study of cTNR markers to better understand their broader relevance to the evolution and adaptation of mammals.

Assessing connectivity and the contribution of private lands to protected area networks in the United States.

Current targets for protected area network coverage call for increased protection but lack specificity in terms of criteria for parcel type, placement, and landscape connectivity. We assessed land conservation achieved by protected area networks in the contiguous United States, and assessed whether private lands contributed substantially to network coverage and connectivity given species dispersal abilities. On average, states have 4.1% (range: 0.2% to 15.8%, n = 48) protected area coverage with connectivity ≤10 km. Terrain ruggedness, percent farmland, and population density are the primary determinants of protected area placement, leading to biased representation of land features currently under protection. On average, private protected areas contribute <1.1% (range 0.0% to 9.2%) to current protected area coverage across the United States. We conclude that current levels of protection and connectivity are inadequate to support a functional protected area network, and that increased and strategic selection of private protected areas is necessary to meet conservation planning goals.

Experimental increase in predation risk causes a cascading stress response in free-ranging snowshoe hares.

Extensive research confirms that environmental stressors like predation risk can profoundly affect animal condition and physiology. However, there is a lack of experimental research assessing the suite of physiological responses to risk that may arise under realistic field conditions, leaving a fragmented picture of risk-related physiological change and potential downstream consequences on individuals. We increased predation risk in free-ranging snowshoe hares (Lepus americanus) during two consecutive summers by simulating natural chases using a model predator and monitored hares intensively via radio-telemetry and physiological assays, including measures designed to assess changes in stress physiology and overall condition. Compared to controls, risk-augmented hares had 25.8% higher free plasma cortisol, 15.9% lower cortisol-binding capacity, a greater neutrophil:lymphocyte skew, and a 10.4% increase in glucose. Despite these changes, intra-annual changes in two distinct condition indices, were unaffected by risk exposure. We infer risk-augmented hares compensated for changes in their stress physiology through either compensatory foraging and/or metabolic changes, which allowed them to have comparable condition to controls. Although differences between controls and risk-augmented hares were consistent each year, both groups had heightened stress measures during the second summer, likely reflecting an increase in natural stressors (i.e., predators) in the environment. We show that increased predation risk in free-ranging animals can profoundly alter stress physiology and that compensatory responses may contribute to limiting effects of such changes on condition. Ultimately, our results also highlight the importance of biologically relevant experimental risk manipulations in the wild as a means of assessing physiological responses to natural stressors.

Frog Virus 3 Genomes Reveal Prevalent Recombination between Ranavirus Lineages and Their Origins in Canada.

Ranaviruses are pathogens associated with the decline of amphibian populations across much of their distribution. In North America, frog virus 3 (FV3) is a widely distributed pathogen with wild populations of amphibians harboring different lineages and putative recombinants between FV3 and common midwife toad virus (CMTV). These recombinants have higher pathogenicity, and CMTV-derived genes associated with virulence are reported in wild strains in Canada. However, while FV3 is linked to amphibian die-offs in North America, CMTVs have been reported only in commercial frog farms in North America. We sequenced complete genomes of 18 FV3 isolates from three amphibian species to characterize genetic diversity of the lineages in Canada and infer possible recombinant regions. The 18 FV3 isolates displayed different signals of recombination, varying from none to interspersed recombination with previously isolated CMTV-like viruses. In general, most recombination breakpoints were located within open reading frames (ORFs), generating new ORFs and proteins that were a mixture between FV3 and CMTV. A combined spatial and temporal phylogeny suggests the presence of the FV3 lineage in Canada is relatively contemporary (<100 years), corroborating the hypothesis that both CMTV- and FV3-like viruses spread to North America when the international commercial amphibian trade started. Our results highlight the importance of pathogen surveillance and viral dynamics using full genomes to more clearly understand the mechanisms of disease origin and spread.IMPORTANCE Amphibian populations are declining worldwide, and these declines have been linked to a number of anthropogenic factors, including disease. Among the pathogens associated with amphibian mortality, ranaviruses have caused massive die-offs across continents. In North America, frog virus 3 (FV3) is a widespread ranavirus that can infect wild and captive amphibians. In this study, we sequenced full FV3 genomes isolated from frogs in Canada. We report widespread recombination between FV3 and common midwife toad virus (CMTV). Phylogenies indicate a recent origin for FV3 in Canada, possibly as a result of international amphibian trade.

Genome-scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx.

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high-throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20-33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non-model systems.

Low intraspecific variation of Frog virus 3 with evidence for novel FV3-like isolates in central and northwestern Canada.

Frog virus 3 (FV3) and FV3-like ranaviruses can infect a variety of cold-blooded aquatic species and present a primary threat to amphibians across the globe. Previous studies of FV3-like viruses have largely investigated higher-level phylogenetic distinctions of these pathogens via portions of the conserved major capsid protein (MCP), and the putative virulence gene vIF-2α. Few studies, however, have investigated the spatial distribution of FV3 variants at the population level3-data that can be used to further understand the spatial epidemiology of this disease. In this study, we sequenced the MCP and vIF-2α of 127 FV3-positive amphibians sampled from Canadian water bodies in Ontario, northeastern Alberta, and southern Northwest Territories to explore whether intraspecific genetic variation exists within FV3. There was a lack of variation at the 2 markers across these regions, suggesting that there is a lack of FV3 sequence diversity in Canada, which may hint at a single source of infection that has spread. However, an undocumented variant termed Wood Buffalo ranavirus (WBRV) was detected in samples from 3 sites in Alberta and Northwest Territories that clustered within the FV3-like lineage with 99.3% sequence homology for MCP. For vIF-2α, all sequences were the expected truncated variant except for 6 samples in Ontario. These latter sequences were suggestive of recombination with common midwife toad virus (CMTV). The lack of variation suggests that higher-resolution genome analyses will be required to further explore the spatial spread and intraspecific variation of the disease.

Threat of predation alters aggressive interactions among spotted salamander (Ambystoma maculatum) larvae.

Intraspecific aggression represents a major source of mortality for many animals and is often experienced alongside the threat of predation. The presence of predators can strongly influence ecological systems both directly by consuming prey and indirectly by altering prey behavior or habitat use. As such, the threat of attack by higher level predators may strongly influence agonistic interactions among conspecifics via nonconsumptive (e.g., behaviorally mediated) predator effects. We sought to investigate these interactions experimentally using larval salamanders (Ambystoma maculatum) as prey and dragonfly nymphs (Anax junius) as predators. Specifically, we quantified salamander behavioral responses to perceived predation risk (PPR) from dragonfly nymphs and determined the degree to which PPR influenced intraspecific aggression (i.e., intraspecific biting and cannibalism) among prey. This included examining the effects of predator exposure on the magnitude of intraspecific biting (i.e., extent of tail damage) and the resulting change in performance (i.e., burst swim speed). Salamander larvae responded to PPR by reducing activity and feeding, but did not increase refuge use. Predator exposure did not significantly influence overall survival; however, the pattern of survival differed among treatments. Larvae exposed to PPR experienced less tail damage from conspecifics, and maximum burst swim speed declined as tail damage became more extensive. Thus, escape ability was more strongly compromised by intraspecific aggression occurring in the absence of predation risk. We conclude that multitrophic indirect effects may importantly modulate intraspecific aggression and should be considered when evaluating the effects of intraspecific competition.

Spatiotemporal heterogeneity in prey abundance and vulnerability shapes the foraging tactics of an omnivore.

Prey abundance and prey vulnerability vary across space and time, but we know little about how they mediate predator-prey interactions and predator foraging tactics. To evaluate the interplay between prey abundance, prey vulnerability and predator space use, we examined patterns of black bear (Ursus americanus) predation of caribou (Rangifer tarandus) neonates in Newfoundland, Canada using data from 317 collared individuals (9 bears, 34 adult female caribou, 274 caribou calves). During the caribou calving season, we predicted that landscape features would influence calf vulnerability to bear predation, and that bears would actively hunt calves by selecting areas associated with increased calf vulnerability. Further, we hypothesized that bears would dynamically adjust their foraging tactics in response to spatiotemporal changes in calf abundance and vulnerability (collectively, calf availability). Accordingly, we expected bears to actively hunt calves when they were most abundant and vulnerable, but switch to foraging on other resources as calf availability declined. As predicted, landscape heterogeneity influenced risk of mortality, and bears displayed the strongest selection for areas where they were most likely to kill calves, which suggested they were actively hunting caribou. Initially, the per-capita rate at which bears killed calves followed a type-I functional response, but as the calving season progressed and calf vulnerability declined, kill rates dissociated from calf abundance. In support of our hypothesis, bears adjusted their foraging tactics when they were less efficient at catching calves, highlighting the influence that predation phenology may have on predator space use. Contrary to our expectations, however, bears appeared to continue to hunt caribou as calf availability declined, but switched from a tactic of selecting areas of increased calf vulnerability to a tactic that maximized encounter rates with calves. Our results reveal that generalist predators can dynamically adjust their foraging tactics over short time-scales in response to changing prey abundance and vulnerability. Further, they demonstrate the utility of integrating temporal dynamics of prey availability into investigations of predator-prey interactions, and move towards a mechanistic understanding of the dynamic foraging tactics of a large omnivore.

The Origin and Ecological Function of an Ion Inducing Anti-Predator Behavior in Lithobates Tadpoles.

In aquatic environments, chemical cues are believed to be associated with prey response to predation risk, yet few basic cue compositions are known despite the pronounced ecological and evolutionary significance of such cues. Previous work indicated that negatively-charged ions of m/z 501 are possibly a kairomone that induces anti-predator responses in amphibian tadpoles. However, work described here confirms that this specific ion species m/z 501.2886 is produced by injured tadpoles, exhibits increased spectral intensity with higher tadpole biomass, and is not produced by starved predators. These results indicate the anion is an alarm cue released from tadpoles. High resolution mass spectrometry (HR-MS) revealed a unique elemental composition for [M-H]-, m/z 501.2886, of C26H45O7S- which could not be determined in previous studies using low resolution instruments. Collision induced dissociation of m/z 501 ions formed product ions of m/z 97 and m/z 80, HSO4- and SO3-, respectively, showing the presence of sulfate. Green frog tadpoles, Lithobates clamitans, exposed to the m/z 501 anion or sodium dodecyl sulfate exhibited similar anti-predator responses, suggesting organic sulfate is a tadpole behavior modifier.

Climate change can alter predator-prey dynamics and population viability of prey.

For many organisms, climate change can directly drive population declines, but it is less clear how such variation may influence populations indirectly through modified biotic interactions. For instance, how will climate change alter complex, multi-species relationships that are modulated by climatic variation and that underlie ecosystem-level processes? Caribou (Rangifer tarandus), a keystone species in Newfoundland, Canada, provides a useful model for unravelling potential and complex long-term implications of climate change on biotic interactions and population change. We measured cause-specific caribou calf predation (1990-2013) in Newfoundland relative to seasonal weather patterns. We show that black bear (Ursus americanus) predation is facilitated by time-lagged higher summer growing degree days, whereas coyote (Canis latrans) predation increases with current precipitation and winter temperature. Based on future climate forecasts for the region, we illustrate that, through time, coyote predation on caribou calves could become increasingly important, whereas the influence of black bear would remain unchanged. From these predictions, demographic projections for caribou suggest long-term population limitation specifically through indirect effects of climate change on calf predation rates by coyotes. While our work assumes limited impact of climate change on other processes, it illustrates the range of impact that climate change can have on predator-prey interactions. We conclude that future efforts to predict potential effects of climate change on populations and ecosystems should include assessment of both direct and indirect effects, including climate-predator interactions.

Signatures of selection in mammalian clock genes with coding trinucleotide repeats: Implications for studying the genomics of high-pace adaptation.

Climate change is predicted to affect the reproductive ecology of wildlife; however, we have yet to understand if and how species can adapt to the rapid pace of change. Clock genes are functional genes likely critical for adaptation to shifting seasonal conditions through shifts in timing cues. Many of these genes contain coding trinucleotide repeats, which offer the potential for higher rates of change than single nucleotide polymorphisms (SNPs) at coding sites, and, thus, may translate to faster rates of adaptation in changing environments. We characterized repeats in 22 clock genes across all annotated mammal species and evaluated the potential for selection on repeat motifs in three clock genes (NR1D1,CLOCK, and PER1) in three congeneric species pairs with different latitudinal range limits: Canada lynx and bobcat (Lynx canadensis and L. rufus), northern and southern flying squirrels (Glaucomys sabrinus and G. volans), and white-footed and deer mouse (Peromyscus leucopus and P. maniculatus). Signatures of positive selection were found in both the interspecific comparison of Canada lynx and bobcat, and intraspecific analyses in Canada lynx. Northern and southern flying squirrels showed differing frequencies at common CLOCK alleles and a signature of balancing selection. Regional excess homozygosity was found in the deer mouse at PER1 suggesting disruptive selection, and further analyses suggested balancing selection in the white-footed mouse. These preliminary signatures of selection and the presence of trinucleotide repeats within many clock genes warrant further consideration of the importance of candidate gene motifs for adaptation to climate change.

Moving forward in circles: challenges and opportunities in modelling population cycles.

Population cycling is a widespread phenomenon, observed across a multitude of taxa in both laboratory and natural conditions. Historically, the theory associated with population cycles was tightly linked to pairwise consumer-resource interactions and studied via deterministic models, but current empirical and theoretical research reveals a much richer basis for ecological cycles. Stochasticity and seasonality can modulate or create cyclic behaviour in non-intuitive ways, the high-dimensionality in ecological systems can profoundly influence cycling, and so can demographic structure and eco-evolutionary dynamics. An inclusive theory for population cycles, ranging from ecosystem-level to demographic modelling, grounded in observational or experimental data, is therefore necessary to better understand observed cyclical patterns. In turn, by gaining better insight into the drivers of population cycles, we can begin to understand the causes of cycle gain and loss, how biodiversity interacts with population cycling, and how to effectively manage wildly fluctuating populations, all of which are growing domains of ecological research.