Estimation of population growth or decline in genetically monitored populations

This article introduces a new general method for genealogical inference that samples independent genealogical histories using importance sampling (IS) and then samples other parameters with Markov chain Monte Carlo (MCMC). It is then possible to more easily utilize the advantages of importance sampling in a fully Bayesian framework. The method is applied to the problem of estimating recent changes in effective population size from temporally spaced gene frequency data. The method gives the posterior distribution of effective population size at the time of the oldest sample and at the time of the most recent sample, assuming a model of exponential growth or decline during the interval. The effect of changes in number of alleles, number of loci, and sample size on the accuracy of the method is described using test simulations, and it is concluded that these have an approximately equivalent effect. The method is used on three example data sets and problems in interpreting the posterior densities are highlighted and discussed.

Joint effects of population density and toxicant exposure on population dynamics of Capitella sp I

Very few studies have analyzed the dependence of population growth rate on population density, and even fewer have considered interaction effects of density and other stresses, such as exposure to toxic chemicals. Yet without such studies we cannot know whether chemicals harmful at low density have effects on carrying capacity or, conversely, whether chemicals reducing carrying capacity are also harmful at low density, impeding a population's capacity to recover from disturbance. This study examines the combined effects of population density and a toxicant (fluoranthene) on population growth rate (pgr) and carrying capacity using the deposit-feeding polychaete Capitella sp. I as a test organism. Populations were initiated with a stable age distribution, and population density and age/size distribution were followed during a period of 28 wk. Fluoranthene (FLU), population density, and their interaction influenced population growth rate. Population growth rate declined linearly with the logarithm of population biomass, but the slope of the relationship was steeper for the control populations than for populations exposed to 50 mug FLU/(g sediment dry mass). Populations exposed to 150 mug FLU/(g sediment dry mass) went extinct after 8 wk of exposure. Despite concerns that toxicant effects would be exacerbated at high density, we found the reverse to be the case, and effects of fluoranthene on population growth rate were much reduced in the region of carrying capacity. Fluoranthene did. reduce carrying capacity by 46%, and this could haven important implications for interacting species and/or sediment biogeochemical processes.

Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity

1. We studied a reintroduced population of the formerly critically endangered Mauritius kestrel Falco punctatus Temmink from its inception in 1987 until 2002, by which time the population had attained carrying capacity for the study area. Post-1994 the population received minimal management other than the provision of nestboxes. 2. We analysed data collected on survival (1987-2002) using program MARK to explore the influence of density-dependent and independent processes on survival over the course of the population's development. 3.We found evidence for non-linear, threshold density dependence in juvenile survival rates. Juvenile survival was also strongly influenced by climate, with the temporal distribution of rainfall during the cyclone season being the most influential climatic variable. Adult survival remained constant throughout. 4. Our most parsimonious capture-mark-recapture statistical model, which was constrained by density and climate, explained 75.4% of the temporal variation exhibited in juvenile survival rates over the course of the population's development. 5. This study is an example of how data collected as part of a threatened species recovery programme can be used to explore the role and functional form of natural population regulatory processes. With the improvements in conservation management techniques and the resulting success stories, formerly threatened species offer unique opportunities to further our understanding of the fundamental principles of population ecology.

Population subdivision and molecular sequence variation: Theory and analysis of Drosophila ananassae data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to Mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter Theta to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of Theta, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.

Attention processes in children with movement difficulties, reading difficulties or both

Reading difficulties (RD) and movement difficulties (MD) co-occur more often in clinical populations than expected for independent disorders. In this study, we investigated the pattern of association between attentional processes, RD and MD in a population of 9 year old school children. Children were screened to identify index groups with RD, MD or both, plus a control group. These groups were then tested on a battery of cognitive attention assessments (TEA-Ch). Results confirmed that the occurrence of RD and MD was greater than would be predicted for independent disorders. Additionally, children with MD, whether or not combined with RD, had poor performance on all attention measures when compared with typically developing children. Children with RD only, were no poorer on measures of attention than typical children. The results are discussed with respect to approaches proposed to account for the co-occurrence of disorders.

Climate change and the risks associated with delayed breeding in a tropical wild bird population

There is growing evidence of changes in the timing of important ecological events, such as flowering in plants and reproduction in animals, in response to climate change, with implications for population decline and biodiversity loss. Recent work has shown that the timing of breeding in wild birds is changing in response to climate change partly because individuals are remarkably flexible in their timing of breeding. Despite this work, our understanding of these processes in wild populations remains very limited and biased towards species from temperate regions. Here, we report the response to changing climate in a tropical wild bird population using a long-term dataset on a formerly critically endangered island endemic, the Mauritius kestrel. We show that the frequency of spring rainfall affects the timing of breeding, with birds breeding later in wetter springs. Delays in breeding have consequences in terms of reduced reproductive success as birds get exposed to risks associated with adverse climatic conditions later on in the breeding season, which reduce nesting success. These results, combined with the fact that frequency of spring rainfall has increased by about 60 per cent in our study area since 1962, imply that climate change is exposing birds to the stochastic risks of late reproduction by causing them to start breeding relatively late in the season.

A Wiener-Hopf type factorization for the exponential functional of Levy processes

For a Lévy process ξ=(ξt)t≥0 drifting to −∞, we define the so-called exponential functional as follows: Formula Under mild conditions on ξ, we show that the following factorization of exponential functionals: Formula holds, where × stands for the product of independent random variables, H− is the descending ladder height process of ξ and Y is a spectrally positive Lévy process with a negative mean constructed from its ascending ladder height process. As a by-product, we generate an integral or power series representation for the law of Iξ for a large class of Lévy processes with two-sided jumps and also derive some new distributional properties. The proof of our main result relies on a fine Markovian study of a class of generalized Ornstein–Uhlenbeck processes, which is itself of independent interest. We use and refine an alternative approach of studying the stationary measure of a Markov process which avoids some technicalities and difficulties that appear in the classical method of employing the generator of the dual Markov process.

The structure of paranoia in the general population

Background: Psychotic phenomena appear to form a continuum with normal experience and beliefs, and may build on common emotional interpersonal concerns. Aims: We tested predictions that paranoid ideation is exponentially distributed and hierarchically arranged in the general population, and that persecutory ideas build on more common cognitions of mistrust, interpersonal sensitivity and ideas of reference. Method: Items were chosen from the Structured Clinical Interview for DSM-IV Axis II Disorders (SCID-II) questionnaire and the Psychosis Screening Questionnaire in the second British National Survey of Psychiatric Morbidity (n = 8580), to test a putative hierarchy of paranoid development using confirmatory factor analysis, latent class analysis and factor mixture modelling analysis. Results: Different types of paranoid ideation ranged in frequency from less than 2% to nearly 30%. Total scores on these items followed an almost perfect exponential distribution (r = 0.99). Our four a priori first-order factors were corroborated (interpersonal sensitivity; mistrust; ideas of reference; ideas of persecution). These mapped onto four classes of individual respondents: a rare, severe, persecutory class with high endorsement of all item factors, including persecutory ideation; a quasi-normal class with infrequent endorsement of interpersonal sensitivity, mistrust and ideas of reference, and no ideas of persecution; and two intermediate classes, characterised respectively by relatively high endorsement of items relating to mistrust and to ideas of reference. Conclusions: The paranoia continuum has implications for the aetiology, mechanisms and treatment of psychotic disorders, while confirming the lack of a clear distinction from normal experiences and processes.

Anthropogenic natal environmental effects on life histories in a wild bird population

Recent work suggests that the environment experienced in early life can alter life histories in wild populations [1, 2, 3, 4 and 5], but our understanding of the processes involved remains limited [6 and 7]. Since anthropogenic environmental change is currently having a major impact on wild populations [8], this raises the possibility that life histories may be influenced by human activities that alter environmental conditions in early life. Whether this is the case and the processes involved remain unexplored in wild populations. Using 23 years of longitudinal data on the Mauritius kestrel (Falco punctatus), a tropical forest specialist, we found that females born in territories affected by anthropogenic habitat change shifted investment in reproduction to earlier in life at the expense of late life performance. They also had lower survival rates as young adults. This shift in life history strategy appears to be adaptive, because fitness was comparable to that of other females experiencing less anthropogenic modification in their natal environment. Our results suggest that human activities can leave a legacy on wild birds through natal environmental effects. Whether these legacies have a detrimental effect on populations will depend on life history responses and the extent to which these reduce individual fitness.

Agriculture modifies the seasonal decline of breeding success in a tropical wild bird population

1.Habitat conversion for agriculture is a major driver of biodiversity loss, but our understanding of the demographic processes involved remains poor. We typically investigate the impacts of agriculture in isolation even though populations are likely to experience multiple, concurrent changes in the environment (e.g. land and climate change). Drivers of environmental change may interact to affect demography but the mechanisms have yet to be explored fully in wild populations. 2.Here, we investigate the mechanisms linking agricultural land-use with breeding success using long-term data for the formerly Critically Endangered Mauritius kestrel Falco punctatus; a tropical forest specialist that also occupies agricultural habitats. We specifically focused on the relationship between breeding success, agriculture and the timing of breeding because the latter is sensitive to changes in climatic conditions (spring rainfall), and enables us to explore the interactive effects of different (land and climate) drivers of environmental change. 3.Breeding success, measured as egg survival to fledging, declines seasonally in this population, but we found that the rate of this decline became increasingly rapid as the area of agriculture around a nest site increased. If the relationship between breeding success and agriculture was used in isolation to estimate the demographic impact of agriculture it would significantly under-estimate breeding success in dry (early) springs, and over-estimate breeding success in wet (late) springs. 4.Analysis of prey delivered to nests suggests that the relationship between breeding success and agriculture might be due, in part, to spatial variation in the availability of native, arboreal geckos. 5.Synthesis and applications. Agriculture modifies the seasonal decline in breeding success in this population. As springs are becoming wetter in our study area and since the kestrels breed later in wetter springs, the impact of agriculture on breeding success will become worse over time. Our results suggest that forest restoration designed to reduce the detrimental impacts of agriculture on breeding may also help reduce the detrimental effects of breeding late due to wetter springs. Our results therefore highlight the importance of considering the interactive effects of environmental change when managing wild populations.

The source population for the cusp and cleft/LLBL for southward IMF

The distinction between plasma properties in different dayside regions in the Earth's magnetosphere is of strong interest as it is often indicative of specific physical processes. This is certainly true for the distinction between low latitude boundary layer (LLBL) and cusp plasma, which has been attributed to the effects of plasma diffusion across the magnetopause (LLBL) versus more direct entry of magnetosheath plasma(cusp). It is also the case, however, that quite different plasma regions can result more simply from a common source plasma, and from different stages of temporal evolution of the plasma associated with magnetospheric convection. In this paper, we show that, for southward interplanetary magnetic field (IMF) conditions, the distinction between the cusp and cleft/LLBL at low altitudes may result from;the single process of magnetosheath plasma entry into the magnetosphere on reconnected field lines. The different plasma characteristics of the two regions result from the properties of the source magnetosheath ion distribution and the effects of magnetic reconnection. Using well known properties of the magnetosheath, several predictions concerning the cusp and cleft/ LLBL precipitation are readily derived.

The uniqueness of signature problem in the non-Markov setting

We establish a general framework for a class of multidimensional stochastic processes over [0,1] under which with probability one, the signature (the collection of iterated path integrals in the sense of rough paths) is well-defined and determines the sample paths of the process up to reparametrization. In particular, by using the Malliavin calculus we show that our method applies to a class of Gaussian processes including fractional Brownian motion with Hurst parameter H>1/4, the Ornstein–Uhlenbeck process and the Brownian bridge.

Population dynamics of genetically diverse Plasmodium falciparum lineages: community-based prospective study in rural Amazonia

Temporal changes in the prevalence of antigenic variants in Plasmodium falciparum populations have been interpreted as evidence of immune-mediated frequency-dependent selection, but evolutively neutral processes may generate similar patterns of serotype replacement. Over 4 years, we investigated the population dynamics of P. falciparum polymorphisms the community level by using 11 putatively neutral microsatellite markers. Plasmodium falciparum Populations were less diverse than sympatric P. vivax isolates, with less multiple-clone infections, lower number of alleles per locus and lower Virtual heterozygosity, but both species showed significant multilocus linkage disequilibrium. Evolutively neutral P. falciparum polymorphisms showed a high turnover rate, with few lineages persisting for several months in the population. Similar results had previously been obtained, in the same community, for sympatric P. vivax isolates. In contrast, the prevalence of the 2 dimorphic types of a major antigen, MSP-2, remained remarkably stable throughout the Study period. We Suggest that the relatively fast turnover of parasite lineages represents the typical population dynamics of neutral polymorphisms in small populations, with clear implications for the detection of frequency-dependent selection of polymorphisms.

Recording from Two Neurons: Second-Order Stimulus Reconstruction from Spike Trains and Population Coding

We study the reconstruction of visual stimuli from spike trains, representing the reconstructed stimulus by a Volterra series up to second order. We illustrate this procedure in a prominent example of spiking neurons, recording simultaneously from the two H1 neurons located in the lobula plate of the fly Chrysomya megacephala. The fly views two types of stimuli, corresponding to rotational and translational displacements. Second-order reconstructions require the manipulation of potentially very large matrices, which obstructs the use of this approach when there are many neurons. We avoid the computation and inversion of these matrices using a convenient set of basis functions to expand our variables in. This requires approximating the spike train four-point functions by combinations of two-point functions similar to relations, which would be true for gaussian stochastic processes. In our test case, this approximation does not reduce the quality of the reconstruction. The overall contribution to stimulus reconstruction of the second-order kernels, measured by the mean squared error, is only about 5% of the first-order contribution. Yet at specific stimulus-dependent instants, the addition of second-order kernels represents up to 100% improvement, but only for rotational stimuli. We present a perturbative scheme to facilitate the application of our method to weakly correlated neurons.

RUMOR PROCESSES ON N

We study four discrete-time stochastic systems on N, modeling processes of rumor spreading. The involved individuals can either have an active or a passive role, speaking up or asking for the rumor. The appetite for spreading or hearing the rumor is represented by a set of random variables whose distributions may depend on the individuals. Our goal is to understand-based on the distribution of the random variables-whether the probability of having an infinite set of individuals knowing the rumor is positive or not.