Mating disrupts morning anticipation in Drosophila melanogaster females.

After mating, the physiology of Drosophila females undergo several important changes, some of which are reflected in their rest-activity cycles. To explore the hypothesis that mating modifies the temporal organization of locomotor activity patterns, we recorded fly activity by a video tracking method. Monitoring rest-activity patterns under light/dark (LD) cycles indicated that mated females lose their ability to anticipate the night-day transition, in stark contrast to males and virgins. This postmating response is mediated by the activation of the sex peptide receptor (SPR) mainly on pickpocket (ppk) expressing neurons, since reducing expression of this receptor in these neurons restores the ability to anticipate the LD transition in mated females. Furthermore, we provide evidence of connectivity between ppk+ neurons and the pigment-dispersing factor (PDF)-positive ventral lateral neurons (sLNv), which play a central role in the temporal organization of daily activity. Since PDF has been associated to the generation of the morning activity peak, we hypothesized that the mating signal could modulate PDF levels. Indeed, we confirm that mated females have reduced PDF levels at the dorsal protocerebrum; moreover, SPR downregulation in ppk+ neurons mimics PDF levels observed in males. In sum, our results are consistent with a model whereby mating-triggered signals reach clock neurons in the fly central nervous system to modulate the temporal organization of circadian behavior according to the needs of the new status.

Color discrimination properties arising from optimal decoding in the early stages of visual systems.

In order to interpret animal behavior we need to understand how they see the world. As directly testing color discrimination in animals is difficult and time consuming, it is important to develop theoretical models based in the properties of visual systems. One of the most successful for the prediction of color discrimination behavior is the receptor noise-limited (RNL) model, which depends only on the level of noise in photoreceptors and opponent mechanisms. Here a complementary approach to model construction is used, and optimal color discrimination properties are obtained using information theoretical tools, for the early stages of visual systems. It is shown here that, for most biologically relevant conditions the optimal discrimination function of an ideal observer coincides with the one obtained with the RNL model. Furthermore, within this framework the influence of opponency can be studied by considering models with and without that mechanism but with exactly the same parameters at the level of photoreceptors. As an example, it is shown here that opponency is necessary to explain the discrimination of monochromatic stimuli in honeybees, but not in budgerigars. Since this is a consequence of the narrowing of absorption spectra of photoreceptors, produced by the presence of oil droplets, this could also be true for most other species of birds. This suggests that in order to study opponency in birds, stimuli should have a relatively wide spectrum.

Role of the circadian clock in the statistics of locomotor activity in Drosophila.

In many animals the circadian rhythm of locomotor activity is controlled by an endogenous circadian clock. Using custom made housing and video tracking software in order to obtain high spatial and temporal resolution, we studied the statistical properties of the locomotor activity of wild type and two clock mutants of Drosophila melanogaster. We show here that the distributions of activity and quiescence bouts for the clock mutants in light-dark conditions (LD) are very different from the distributions obtained when there are no external cues from the environment (DD). In the wild type these distributions are very similar, showing that the clock controls this aspect of behavior in both regimes (LD and DD). Furthermore, the distributions are very similar to those reported for Wistar rats. For the timing of events we also observe important differences, quantified by how the event rate distributions scale for increasing time windows. We find that for the wild type these distributions can be rescaled by the same function in DD as in LD. Interestingly, the same function has been shown to rescale the rate distributions in Wistar rats. On the other hand, for the clock mutants it is not possible to rescale the rate distributions, which might indicate that the extent of circadian control depends on the statistical properties of activity and quiescence.

Insights About Head Lice Transmission From Field Data and Mathematical Modeling.

Head lice infest millions of school-age children every year, both in developed and developing countries. However, little is known about the number of lice transferred among children during school activities, because direct methods to study this are almost impossible to implement. This issue has been addressed following an indirect method, which consist in collecting data of real infestation from several children groups and using a mathematical model of lice colonies to infer how the infestation observed might have evolved. By determining the events that would most likely lead to infestations as those observed, we find that severe infestations are most likely initiated by a relatively large number of lice transferred at the same moment or within relatively short time spans. In turn, analysis of the data obtained from screenings of the same groups of children a few days apart shows evidence of such transmission events. Interestingly, only children with severe infestations could harbor the lice necessary for this type of transmission. Thus, they play the same role as 'superspreaders' in epidemiology. As part of our experimental study it is also shown that a simple procedure of combing can be very effective to remove all mobile lice, and thus could be used as an effective preventive measure against those severe infestations that are responsible for the spread of pediculosis.

Spontaneous circadian rhythms in a cold-adapted natural isolate of Aureobasidium pullulans.

Circadian systems enable organisms to synchronize their physiology to daily and seasonal environmental changes relying on endogenous pacemakers that oscillate with a period close to 24 h even in the absence of external timing cues. The oscillations are achieved by intracellular transcriptional/translational feedback loops thoroughly characterized for many organisms, but still little is known about the presence and characteristics of circadian clocks in fungi other than Neurospora crassa. We sought to characterize the circadian system of a natural isolate of Aureobasidium pullulans, a cold-adapted yeast bearing great biotechnological potential. A. pullulans formed daily concentric rings that were synchronized by light/dark cycles and were also formed in constant darkness with a period of 24.5 h. Moreover, these rhythms were temperature compensated, as evidenced by experiments conducted at temperatures as low as 10 °C. Finally, the expression of clock-essential genes, frequency, white collar-1, white collar-2 and vivid was confirmed. In summary, our results indicate the existence of a functional circadian clock in A. pullulans, capable of sustaining rhythms at very low temperatures and, based on the presence of conserved clock-gene homologues, suggest a molecular and functional relationship to well-described circadian systems.

A mathematical model of communication between groups of circadian neurons in Drosophila melanogaster.

In the fruit fly, circadian behavior is controlled by a small number of specialized neurons, whose molecular clocks are relatively well known. However, much less is known about how these neurons communicate among themselves. In particular, only 1 circadian neuropeptide, pigment-dispersing factor (PDF), has been identified, and most aspects of its interaction with the molecular clock remain to be elucidated. Furthermore, it is speculated that many other peptides should contribute to circadian communication. We have developed a relatively detailed model of the 2 main groups of circadian pacemaker neurons (sLNvs and LNds) to investigate these issues. We have proposed many possible mechanisms for the interaction between the synchronization factors and the molecular clock, and we have compared the outputs with the experimental results reported in the literature both for the wild-type and PDF-null mutant. We have studied how different the properties of each neuron should be to account for the observations reported for the sLNvs in the mutant. We have found that only a few mechanisms, mostly related to the slowing down of nuclear entry of a circadian protein, can synchronize neurons that present these differences. Detailed immunofluorescent recordings have suggested that, whereas in the mutant, LNd neurons are synchronized, in the wild-type, a subset of the LNds oscillate faster than the rest. With our model, we find that a more likely explanation for the same observations is that this subset is being driven outside its synchronization range and displays therefore a complex pattern of oscillation.

Intransitivity and coexistence in four species cyclic games.

Intransitivity is a property of connected, oriented graphs representing species interactions that may drive their coexistence even in the presence of competition, the standard example being the three species Rock-Paper-Scissors game. We consider here a generalization with four species, the minimum number of species allowing other interactions beyond the single loop (one predator, one prey). We show that, contrary to the mean field prediction, on a square lattice the model presents a transition, as the parameter setting the rate at which one species invades another changes, from a coexistence to a state in which one species gets extinct. Such a dependence on the invasion rates shows that the interaction graph structure alone is not enough to predict the outcome of such models. In addition, different invasion rates permit to tune the level of transitiveness, indicating that for the coexistence of all species to persist, there must be a minimum amount of intransitivity.

Modelling the effect of phosphorylation on the circadian clock of Drosophila.

It is by now well known that, at the molecular level, the core of the circadian clock of most living species is a negative feedback loop where some proteins inhibit their own transcription. However, it has recently been shown that post-translational processes, such as phosphorylations, are essential for a correct timing of the clock. Depending on which sites of a circadian protein are phosphorylated, different properties such as degradation, nuclear localization and repressing power can be altered. Furthermore, phosphorylation domains can be related in a positive way, giving rise to consecutive phosphorylations, or in a negative way, hindering phosphorylation at other domains. Here we present a simple mathematical model of a circadian protein having two mutually exclusive domains of phosphorylation. We show that the system has limit cycles that arise from a unique fixed point through a Hopf bifurcation. We find a set of parameters, with realistic values, for which the limit cycle has the same period as the wild type circadian oscillations of the fruit fly. The domains act as a switch, in the sense that alterations in their phosphorylation can alter the period of circadian oscillation in opposite ways, increasing or decreasing the period of the wild type oscillations. In particular, we show that our model is able to reproduce some of the experimental results found for switch-like phosphorylations of the PER protein of the circadian clock of the fly Drosophila melanogaster.

Influence of network dynamics on the spread of sexually transmitted diseases.

Network epidemiology often assumes that the relationships defining the social network of a population are static. The dynamics of relationships is only taken indirectly into account by assuming that the relevant information to study epidemic spread is encoded in the network obtained, by considering numbers of partners accumulated over periods of time roughly proportional to the infectious period of the disease. On the other hand, models explicitly including social dynamics are often too schematic to provide a reasonable representation of a real population, or so detailed that no general conclusions can be drawn from them. Here, we present a model of social dynamics that is general enough so its parameters can be obtained by fitting data from surveys about sexual behaviour, but that can still be studied analytically, using mean-field techniques. This allows us to obtain some general results about epidemic spreading. We show that using accumulated network data to estimate the static epidemic threshold lead to a significant underestimation of that threshold. We also show that, for a dynamic network, the relative epidemic threshold is an increasing function of the infectious period of the disease, implying that the static value is a lower bound to the real threshold. A practical example is given of how to apply the model to the study of a real population.

Of lice and math: using models to understand and control populations of head lice.

In this paper we use detailed data about the biology of the head louse (pediculus humanus capitis) to build a model of the evolution of head lice colonies. Using theory and computer simulations, we show that the model can be used to assess the impact of the various strategies usually applied to eradicate head lice, both conscious (treatments) and unconscious (grooming). In the case of treatments, we study the difference in performance that arises when they are applied in systematic and non-systematic ways. Using some reasonable simplifying assumptions (as random mixing of human groups and the same mobility for all life stages of head lice other than eggs) we model the contagion of pediculosis using only one additional parameter. It is shown that this parameter can be tuned to obtain collective infestations whose characteristics are compatible with what is given in the literature on real infestations. We analyze two scenarios: One where group members begin treatment when a similar number of lice are present in each head, and another where there is one individual who starts treatment with a much larger threshold ("superspreader"). For both cases we assess the impact of several collective strategies of treatment.

Contact switching as a control strategy for epidemic outbreaks.

We study the effects of switching social contacts as a strategy to control epidemic outbreaks. Connections between susceptible and infective individuals can be broken by either individual, and then reconnected to a randomly chosen member of the population. It is assumed that the reconnecting individual has no previous information on the epidemiological condition of the new contact. We show that reconnection can completely suppress the disease, both by continuous and discontinuous transitions between the endemic and the infection-free states. For diseases with an asymptomatic phase, we analyze the conditions for the suppression of the disease, and show that-even when these conditions are not met-the increase of the endemic infection level is usually rather small. We conclude that, within some simple epidemiological models, contact switching is a quite robust and effective control strategy. This suggests that it may also be an efficient method in more complex situations.

Infection spreading in a population with evolving contacts.

We study the spreading of an infection within an SIS epidemiological model on a network. Susceptible agents are given the opportunity of breaking their links with infected agents. Broken links are either permanently removed or reconnected with the rest of the population. Thus, the network coevolves with the population as the infection progresses. We show that a moderate reconnection frequency is enough to completely suppress the infection. A partial, rather weak isolation of infected agents suffices to eliminate the endemic state.

Retrieval-time properties of the Little-Hopfield model and their physiological relevance.

We perform an extensive numerical investigation on the retrieval dynamics of the synchronous Hopfield model, also known as Little-Hopfield model, up to sizes of 2(18) neurons. Our results correct and extend much of the early simulations on the model. We find that the average convergence time has a power law behavior for a wide range of system sizes, whose exponent depends both on the network loading and the initial overlap with the memory to be retrieved. Surprisingly, we also find that the variance of the convergence time grows as fast as its average, making it a non-self-averaging quantity. Based on the simulation data we differentiate between two definitions for memory retrieval time, one that is mathematically strict, tau(c), the number of updates needed to reach the attractor whose properties we just described, and a second definition correspondent to the time tau(eta) when the network stabilizes within a tolerance threshold eta such that the difference of two consecutive overlaps with a stored memory is smaller that eta. We show that the scaling relationships between tau(c) and tau(eta) and the typical network parameters as the memory load alpha or the size of the network N vary greatly, being tau(eta) relatively insensitive to system sizes and loading. We propose tau(eta) as the physiological realistic measure for the typical attractor network response.

Topology, phase transitions, and the spherical model.

The topological hypothesis states that phase transitions should be related to changes in the topology of configuration space. The necessity of such changes has already been demonstrated. We characterize exactly the topology of the configuration space of the short range Berlin-Kac spherical model, for spins lying in hypercubic lattices of dimension d. We find a continuum of changes in the topology and also a finite number of discontinuities in some topological functions. We show, however, that these discontinuities do not coincide with the phase transitions which happen for d > or = 3, and conversely, that no topological discontinuity can be associated with them. This is the first short range, confining potential for which the existence of special topological changes are shown not to be sufficient to infer the occurrence of a phase transition.

Properties of dense partially random graphs.

We study the properties of random graphs where for each vertex a neighborhood has been previously defined. The probability of an edge joining two vertices depends on whether the vertices are neighbors or not, as happens in small-world graphs (SWG's). But we consider the case where the average degree of each node is of order of the size of the graph (unlike SWG's, which are sparse). This allows us to calculate the mean distance and clustering, which are qualitatively similar (although not in such a dramatic scale range) to the case of SWG's. We also obtain analytically the distribution of eigenvalues of the corresponding adjacency matrices. This distribution is discrete for large eigenvalues and continuous for small eigenvalues. The continuous part of the distribution follows a semicircle law, whose width is proportional to the "disorder" of the graph, whereas the discrete part is simply a rescaling of the spectrum of the substrate. We apply our results to the calculation of the mixing rate and the synchronizability threshold.

Exploratory behavior, trap models, and glass transitions.

A random walk is performed on a disordered landscape composed of N sites randomly and uniformly distributed inside a d-dimensional hypercube. The walker hops from one site to another with probability proportional to exp[-betaE(D)], where beta=1/T is the inverse of a formal temperature and E(D) is an arbitrary cost function which depends on the hop distance D. Analytic results indicate that, if E(D)=D(d) and N--> infinity, there exists a glass transition at beta(d)=pi(d/2)/[(d/2)Gamma(d/2)]. Below T(d), the average trapping time diverges and the system falls into an out-of-equilibrium regime with aging phenomena. A Lévy flight scenario and applications of exploratory behavior are considered.

Escaping from cycles through a glass transition.

A random walk is performed over a disordered media composed of N sites random and uniformly distributed inside a d-dimensional hypercube. The walker cannot remain in the same site and hops to one of its n neighboring sites with a transition probability that depends on the distance D between sites according to a cost function E(D). The stochasticity level is parametrized by a formal temperature T. In the case T=0, the walk is deterministic and ergodicity is broken: the phase space is divided in a O(N) number of attractor basins of two-cycles that trap the walker. For d=1, analytic results indicate the existence of a glass transition at T(1)=1/2 as N--> infinity. Below T1, the average trapping time in two-cycles diverges and an out-of-equilibrium behavior appears. Similar glass transitions occur in higher dimensions when the right cost function is chosen. We also present some results for the statistics of distances for Poisson spatial point processes.