Status survey of the critically endangered Iberian lynx <Emphasis Type="Italic">Lynx pardinus</Emphasis> in Portugal

A survey of the Iberian lynx was conducted in Portugal between January 2002 and November 2003 in order to define lynx status and distribution baselines. Intensive search for lynx scats, scat DNA analysis, and camera trapping were used in areas of potential lynx presence. Over 4,200 km were investigated during a global searching effort of 1,975 man-hours. DNA obtained from 168 scats was analyzed, producing no positive lynx amplifications. Camera trapping represented a total effort of 5,647 camera-days in three potential lynx areas, producing no positive detections. Although results cannot confirm the species’ extinction, the scenario is rather pessimistic and the Iberian lynx is probably no longer present in Portugal. Presently, considering the Portuguese lynx’s historical range, the only significantly suitable areas for the species are located in the southeastern part of the country bordering Andalusia where no recent evidence of lynx presence was recorded.     

Stable long-period cycling and complex dynamics in a single-locus fertility model with genomic imprinting

Although long-period population size cycles and chaotic fluctuations in abundance are common in ecological models, such dynamics are uncommon in simple population-genetic models where convergence to a fixed equilibrium is most typical. When genotype-frequency cycling does occur, it is most often due to frequency-dependent selection that results from individual or species interactions. In this paper, we demonstrate that fertility selection and genomic imprinting are sufficient to generate a Hopf bifurcation and complex genotype-frequency cycling in a single-locus population-genetic model. Previous studies have shown that on its own, fertility selection can yield stable two-cycles but not long-period cycling characteristic of a Hopf bifurcation. Genomic imprinting, a molecular mechanism by which the expression of an allele depends on the sex of the donating parent, allows fitness matrices to be nonsymmetric, and this additional flexibility is crucial to the complex dynamics we observe in this fertility selection model. Additionally, we find under certain conditions that stable oscillations and a stable equilibrium point can coexist. These dynamics are characteristic of a Chenciner (generalized Hopf) bifurcation. We believe this model to be the simplest population-genetic model with such dynamics.     

Seven decades of southern range dynamics of Canada lynx

The range of the Canada lynx (Lynx canadensis) has contracted substantially from its historical range. Using harvest records, we found that the southern range of the lynx in Ontario in the late 1940s collapsed and then, in a short period of time, increased to its largest extent in the mid‐1960s when the lynx range spread south of the boreal forest for a decade. After this expansion, the southern range contracted northwards beginning in the 1970s. Most recently, there has been a slight expansion between 2010 and 2017. We have attributed these dynamics on the southern range periphery to the fluctuation of the boreal lynx population in the core of the species'' range. In addition, connectivity to boreal lynx populations and snow depth seemed to condition whether the lynx expanded into an area. However, we did not find any evidence to suggest that these changes were due to anthropogenic landscape disturbances or competition. The boreal lynx population does not reach the peak abundance it once did, without which we would not expect to see large expansions of the southern lynx range as in the mid‐1960s. Our results suggest that the southern lynx range in Ontario has been driven by the magnitude of the boreal lynx population cycle, connectivity to the boreal forest, and snow conditions. Future persistence of lynx in the southern range periphery will likely depend on dynamics in the range core.     

Quantifying noise-induced stability of a cortical fast-spiking cell model with Kv3-channel-like current

Population oscillations in neural activity in the gamma (>30 Hz) and higher frequency ranges are found over wide areas of the mammalian cortex. Recently, in the somatosensory cortex, the details of neural connections formed by several types of GABAergic interneurons have become apparent, and they are believed to play a significant role in generating these oscillations through synaptic and gap-junctional interactions. However, little is known about the mechanism of how such oscillations are maintained stably by particular interneurons and by their local networks, in a noisy environment with abundant synaptic inputs. To obtain more insight into this, we studied a fast-spiking (FS)-cell model including Kv3-channel-like current, which is a distinctive feature of these cells, from the viewpoint of nonlinear dynamical systems. To examine the specific role of the Kv3-channel in determining oscillation properties, we analyzed basic properties of the FS-cell model, such as the bifurcation structure and phase resetting curves (PRCs). Furthermore, to quantitatively characterize the oscillation stability under noisy fluctuations mimicking small fast synaptic inputs, we applied a recently developed method from random dynamical system theory to estimate Lyapunov exponents, both for the original four-dimensional dynamics and for a reduced one-dimensional phase-equation on the circle. The results indicated that the presence of the Kv3-channel-like current helps to regulate the stability of noisy neural oscillations and a transient-period length to stochastic attractors.     

The cost of migratory prey: seasonal changes in semi‐domestic reindeer distribution influences breeding success of Eurasian lynx in northern Norway

Migratory prey is a widespread phenomenon that has implications for predator–prey interactions. By creating large temporal variation in resource availability between seasons it becomes challenging for carnivores to secure a regular year‐round supply of food. Some predators may respond by following their migratory prey, however, most predators are sedentary and experience strong seasonal variation in resource availability. Increased predation on alternative prey may dampen such seasonal resource fluctuations, but reduced reproduction rates in predators is a predicted consequence of migratory primary prey behavior that has received little empirical attention. We used data from 23 GPS collared Eurasian lynx Lynx lynx monitored during 2007–2013 in northern Norway, to examine how spatio‐temporal variation in the migratory behavior of semi‐domestic reindeer Rangifer tarandus influences lynx spatial organization and reproductive success using estimates of seasonal home range overlap and breeding success. We found that lynx of both sexes maintained seasonally stable home ranges and exhibited site fidelity across years, independent of whether they had access to reindeer throughout the year or experienced a scarcity of reindeer in winter due to migration. However, lynx without access to reindeer in winter showed a decreased probability of reproducing and a tendency for lowered kitten survival into their first winter, when compared to female lynx with reindeer available year around. This supports the hypothesis that sedentary predators experience demographic costs in systems with migratory primary prey. Changes in the migratory behavior of ungulates, including disrupted migrations, is therefore likely to have bottom–up effects on the population dynamics of sedentary predators as well as the previously documented consequences for ungulate population dynamics.     

Stochastic amplification of fluctuations in cortical up-States

Cortical neurons are bistable; as a consequence their local field potentials can fluctuate between quiescent and active states, generating slow [Formula: see text] Hz oscillations which are widely known as transitions between Up and Down States. Despite a large number of studies on Up-Down transitions, deciphering its nature, mechanisms and function are still today challenging tasks. In this paper we focus on recent experimental evidence, showing that a class of spontaneous oscillations can emerge within the Up states. In particular, a non-trivial peak around [Formula: see text] Hz appears in their associated power-spectra, what produces an enhancement of the activity power for higher frequencies (in the [Formula: see text] Hz band). Moreover, this rhythm within Ups seems to be an emergent or collective phenomenon given that individual neurons do not lock to it as they remain mostly unsynchronized. Remarkably, similar oscillations (and the concomitant peak in the spectrum) do not appear in the Down states. Here we shed light on these findings by using different computational models for the dynamics of cortical networks in presence of different levels of physiological complexity. Our conclusion, supported by both theory and simulations, is that the collective phenomenon of "stochastic amplification of fluctuations" - previously described in other contexts such as Ecology and Epidemiology - explains in an elegant and parsimonious manner, beyond model-dependent details, this extra-rhythm emerging only in the Up states but not in the Downs.     

Cryptic population structure in a large,mobile mammalian predator: the Scandinavian lynx

The Eurasian lynx (Lynx lynx) is an example of a species that has gone through a severe bottleneck, leading to near extinction in Scandinavia around 1930-- a pattern shared with several other large carnivorous mammals. Here we extend previous genetic analyses of northern European lynx, confirming that lynx from the Scandinavian Peninsula represent a distinct clade differing clearly from European conspecifics. Furthermore, and despite a recent bottleneck and subsequent range expansion, we detect marked genetic differentiation within Scandinavia. This differentiation is largely manifested as a north-south gradient, with a linear increase in the quantity FST/(1 - FST). Aided by computer simulations we find that this pattern is unlikely to have arisen by random genetic drift in the short time since lynx started to expand in the 1950s, suggesting that the spatial structure may predate the bottleneck. Individual-based analyses indicate that, instead of a continuous gradient, Scandinavian lynx may be structured into three more or less distinct groups, possibly corresponding to northern, central and southern subpopulations. The presence of such structuring was unknown previously and was unexpected from general considerations on the mobility of the species, historical data and the absence of geographical barriers. Our study demonstrates how molecular markers may be used to detect cryptic population structure, invisible using traditional methods.     

Response to temporal parameter fluctuations in biochemical networks

Metabolic response coefficients describe how variables in metabolic systems, like steady state concentrations, respond to small changes of kinetic parameters. To extend this concept to temporal parameter fluctuations, we define spectral response coefficients that relate Fourier components of concentrations and fluxes to Fourier components of the underlying parameters. It is also straightforward to generalize other concepts from metabolic control theory, such as control coefficients with their summation and connectivity theorems. The first-order response coefficients describe forced oscillations caused by small harmonic oscillations of single parameters: they depend on the driving frequency and comprise the phases and amplitudes of the concentrations and fluxes. Close to a Hopf bifurcation, resonance can occur: as an example, we study the spectral densities of concentration fluctuations arising from the stochastic nature of chemical reactions. Second-order response coefficients describe how perturbations of different frequencies interact by mode coupling, yielding higher harmonics in the metabolic response. The temporal response to small parameter fluctuations can be computed by Fourier synthesis. For a model of glycolysis, this approximation remains fairly accurate even for large relative fluctuations of the parameters.     

猞猁在云南省野外分布的首例影像记录

<正>猞猁(Lynx lynx Linnaeus, 1758)为食肉目(CARNIVORA)猫科(Felidae)猞猁属(Lynx)的中等体型哺乳动物,主要分布于欧洲、北亚、西亚及中亚地区(蒋志刚等,2015)。猞猁分布较广,数量多,为IUCN红色名录的无危物种(Least Concern, LC)(Breitenmoser et al., 2015)。亚洲地区的猞猁主要栖息在山地森林、无林裸岩地带、密集灌木丛和草原,分布海拔可达4 500 m (汪松,1998)。在中国,猞猁分布于东北、华北、西北和西南地区(高耀亭等,1987;王应祥,2003),野外种群数量较少,属于国家二级重点保护野生动物,     

Hybridization Between Canada Lynx and Bobcats: Genetic Results and Management Implications

Hybridization between taxonomically similar species is an often-overlooked mechanism limiting the recovery of threatened and endangered species. We present molecular genetic data for the first time demonstrating that Canada lynx and bobcats hybridize in the wild. We verify that two microsatellite loci Lc106 and Lc110 have non-overlapping allele ranges between Canada lynx and bobcats, and that three putative lynx from Minnesota contain DNA from both bobcats and lynx. Additionally, we use a published test for the 16S rRNA region of mitochondrial DNA (mtDNA) to determine the maternal species; all hybrids had lynx mothers. Fifteen per cent (3/20) of our putative lynx samples were hybrids, although these data are not from a representative sampling effort. Hybridization may be an under-appreciated factor limiting the distribution and recovery of lynx. The presence of hybrids is thus a new factor in the population management of both species with potential implications for hunting and trapping of bobcats.     

Hysteresis Dynamics,Bursting Oscillations and Evolution to Chaotic Regimes

This article describes new aspects of hysteresis dynamics which have been uncovered through computer experiments. There are several motivations to be interested in fast-slow dynamics. For instance, many physiological or biological systems display different time scales. The bursting oscillations which can be observed in neurons, beta-cells of the pancreas and population dynamics are essentially studied via bifurcation theory and analysis of fast-slow systems (Keener and Sneyd, 1998; Rinzel, 1987). Hysteresis is a possible mechanism to generate bursting oscillations. A first part of this article presents the computer techniques (the dotted-phase portrait, the bifurcation of the fast dynamics and the wave form) we have used to represent several patterns specific to hysteresis dynamics. This framework yields a natural generalization to the notion of bursting oscillations where, for instance, the active phase is chaotic and alternates with a quiescent phase. In a second part of the article, we emphasize the evolution to chaos which is often associated with bursting oscillations on the specific example of the Hindmarsh-Rose system. This evolution to chaos has already been studied with classical tools of dynamical systems but we give here numerical evidence on hysteresis dynamics and on some aspects of the wave form. The analytical proofs will be given elsewhere.     

An empirically based model for latitudinal gradient in vole population dynamics

Vole dynamics in northern Europe exhibit a well-defined geographical gradient, with oscillatory populations being confined to high latitudes. It has been proposed that oscillations in northern vole populations are driven by their interaction with specialist predators (weasels), while the more southern rodent populations are relatively stable because of regulation by generalist predators. We tested this generalist/specialist predation hypothesis by constructing an empirically based model for vole population dynamics, estimating its parameters, and making predictions about the quantitative pattern of the latitudinal shift in vole dynamics. Our results indicated that the model accurately predicted the latitudinal shift in the amplitude and periodicity of population fluctuations. Moreover, the model predicted that vole dynamics should shift from stable to chaotic as latitude is increased, a result in agreement with nonlinear time-series analysis of the data. The striking success of the model at predicting the shifts in amplitude and stability along the geographical gradient in northern Europe provides strong support for the key role of specialist and generalist predators in vole population dynamics.     

Oscillatory dynamics in rock-paper-scissors games with mutations

We study the oscillatory dynamics in the generic three-species rock-paper-scissors games with mutations. In the mean-field limit, different behaviors are found: (a) for high mutation rate, there is a stable interior fixed point with coexistence of all species; (b) for low mutation rates, there is a region of the parameter space characterized by a limit cycle resulting from a Hopf bifurcation; (c) in the absence of mutations, there is a region where heteroclinic cycles yield oscillations of large amplitude (not robust against noise). After a discussion on the main properties of the mean-field dynamics, we investigate the stochastic version of the model within an individual-based formulation. Demographic fluctuations are therefore naturally accounted and their effects are studied using a diffusion theory complemented by numerical simulations. It is thus shown that persistent erratic oscillations (quasi-cycles) of large amplitude emerge from a noise-induced resonance phenomenon. We also analytically and numerically compute the average escape time necessary to reach a (quasi-)cycle on which the system oscillates at a given amplitude.     

Transient dynamics of genetic regulatory networks

We present an approximation scheme for deriving reaction rate equations of genetic regulatory networks. This scheme predicts the timescales of transient dynamics of such networks more accurately than does standard quasi-steady state analysis by introducing prefactors to the ODEs that govern the dynamics of the protein concentrations. These prefactors render the ODE systems slower than their quasi-steady state approximation counterparts. We introduce the method by examining a positive feedback gene regulatory network, and show how the transient dynamics of this network are more accurately modeled when the prefactor is included. Next, we examine the repressilator, a genetic oscillator, and show that the period, amplitude, and bifurcation diagram defining the onset of the oscillations are better estimated by the prefactor method. Finally, we examine the consequences of the method to the dynamics of reduced models of the phage lambda switch, and show that the switching times between the two states is slowed by the presence of the prefactor that arises from protein multimerization and DNA binding.     

Mathematical investigation of IP<Subscript>3</Subscript>-dependent calcium dynamics in astrocytes

We study evoked calcium dynamics in astrocytes, a major cell type in the mammalian brain. Experimental evidence has shown that such dynamics are highly variable between different trials, cells, and cell subcompartments. Here we present a qualitative analysis of a recent mathematical model of astrocyte calcium responses. We show how the major response types are generated in the model as a result of the underlying bifurcation structure. By varying key channel parameters, mimicking blockers used by experimentalists, we manipulate this underlying bifurcation structure and predict how the distributions of responses can change. We find that store-operated calcium channels, plasma membrane bound channels with little activity during calcium transients, have a surprisingly strong effect, underscoring the importance of considering these channels in both experiments and mathematical settings. Variation in the maximum flow in different calcium channels is also shown to determine the range of stable oscillations, as well as set the range of frequencies of the oscillations. Further, by conducting a randomized search through the parameter space and recording the resulting calcium responses, we create a database that can be used by experimentalists to help estimate the underlying channel distribution of their cells.     

Global dynamics and bifurcation analysis of a host–parasitoid model with strong Allee effect

In this paper, we study the global dynamics and bifurcations of a two-dimensional discrete time host–parasitoid model with strong Allee effect. The existence of fixed points and their stability are analysed in all allowed parametric region. The bifurcation analysis shows that the model can undergo fold bifurcation and Neimark–Sacker bifurcation. As the parameters vary in a small neighbourhood of the Neimark–Sacker bifurcation condition, the unique positive fixed point changes its stability and an invariant closed circle bifurcates from the positive fixed point. From the viewpoint of biology, the invariant closed curve corresponds to the periodic or quasi-periodic oscillations between host and parasitoid populations. Furthermore, it is proved that all solutions of this model are bounded, and there exist some values of the parameters such that the model has a global attractor. These theoretical results reveal the complex dynamics of the present model.