Worker demography in a large-colony, swarm-founding wasp

Neotropical swarm-founding wasps build nests enclosed in a covering envelope, which makes it difficult to count individual births and deaths. Thus, knowledge of worker demography is very limited for swarm-founding species compared with that for independent-founding species. In this study, we explored the worker demography of the swarm-founding wasp Polybia paulista, the colony size of which usually exceeds several thousand adults. We considered each wasp colony as an open-population and estimated the survival probability, recruitment rate, and population size of workers using the developments of the Cormack-Jolly-Seber model. We found that capture probability varied considerably among the workers, probably due to age polyethism and/or task specialization. The daily survival rate of workers was high (around 0.97) throughout the season and was not related to the phase of colony development. On the other hand, the recruitment rate ranged from 0 to 0.37, suggesting that worker production was substantially less important than worker survival in determining worker population fluctuations. When we compared survival rates among worker groups of one colony, the mean daily survival rate was lower for founding workers than for progeny workers and tended to be higher in progeny workers that emerged in winter. These differences in survivorship patterns among worker cohorts would be related to worker foraging activity and/or level of parasitism.

Maned wolf survival rate in central Brazil

Although many carnivores are of conservation concern, most are poorly studied. The maned wolf Chrysocyon brachyurus Illiger, 1811 is the largest South American canid with a broad distribution; however, the largest portion of its range is in the Brazilian Cerrado savannah, where due to intensive agricultural expansion, it is threatened by habitat loss. Maned wolf population trends are virtually unknown. We analyzed radio telemetry data from a 13-year study in Emas National Park, central Brazil, with Burnham`s live recapture/dead recovery models in the program MARK to obtain the first analytically sound estimate of the apparent maned wolf survival rate. We constructed 16 candidate models including variation in survival rate and resighting probability associated with an individual`s sex or age and year of study. Apparent adult survival rate throughout the study ranged from 0.28 (se=0.08) to 0.97 (se=0.06). There was no evidence for sex specificity but strong support for time variation. Model weights supported an age effect and the subadult survival rate was 0.63 (se=0.15). Results indicate similar life patterns for male and female maned wolves and similar mortality risks for adults and subadults in the study area. The observed temporal fluctuations of adult survival rate are important for population dynamics as they decrease average population growth rates. Population dynamics are central for conservation planning and our results are an important step towards a better understanding of the maned wolf`s ecology.

Density of the Near Threatened jaguar Panthera onca in the caatinga of north-eastern Brazil

We report the first estimate of jaguar density in the semi-arid caatinga biome of north-eastern Brazil. During August-October 2007, in the Serra da Capivara National Park, we used camera traps to identify and count jaguars. Jaguar abundance and density were calculated using mark-recapture models. In a sampling effort of 1,249 camera-trap-nights we identified 12 adult jaguars and estimated an 2 abundance of 14 +/- 3.6 jaguars in an area of 524 km(2), i.e. a density of 2.67 +/- 1.00 jaguars per 100 km(2). This estimate is higher than in most other Brazilian biomes and indicates Serra da Capivara National Park as an important reserve for protecting jaguars in north-eastern Brazil.

Maned Wolf Density in a Central Brazilian Grassland Reserve

The maned wolf (Chrysocyon brachyurus) inhabits the savannah-like Cerrado of Brazil. Although 80% of this biome has already been affected by human activity, little is known about maned wolf abundance. Using mark-recapture models, we obtained the first density estimate from central Brazil, the core of maned wolf distribution. With 5.19 individuals/100 km(2), even large reserves support only small maned wolf populations. Therefore, long-term conservation of the maned wolf depends on land management outside of reserves. ( JOURNAL OF WILDLIFE MANAGEMENT 73( 1): 68-71; 2009)

Conditional male dimorphism and alternative reproductive tactics in a Neotropical arachnid (Opiliones)

In arthropods, most cases of morphological dimorphism within males are the result of a conditional evolutionarily stable strategy (ESS) with status-dependent tactics. In conditionally male-dimorphic species, the status` distributions of male morphs often overlap, and the environmentally cued threshold model (ET) states that the degree of overlap depends on the genetic variation in the distribution of the switchpoints that determine which morph is expressed in each value of status. Here we describe male dimorphism and alternative mating behaviors in the harvestman Serracutisoma proximum. Majors express elongated second legs and use them in territorial fights; minors possess short second legs and do not fight, but rather sneak into majors` territories and copulate with egg-guarding females. The static allometry of second legs reveals that major phenotype expression depends on body size (status), and that the switchpoint underlying the dimorphism presents a large amount of genetic variation in the population, which probably results from weak selective pressure on this trait. With a mark-recapture study, we show that major phenotype expression does not result in survival costs, which is consistent with our hypothesis that there is weak selection on the switchpoint. Finally, we demonstrate that switchpoint is independent of status distribution. In conclusion, our data support the ET model prediction that the genetic correlation between status and switchpoint is low, allowing the status distribution to evolve or to fluctuate seasonally, without any effect on the position of the mean switchpoint.

Dispersal and survival of Nyssomyia intermedia and Nyssomyia neivai (Diptera: Psychodidae: Phlebotominae) in a cutaneous leishmaniasis endemic area of the speleological province of the Ribeira Valley, state of São Paulo, Brazil

The dispersal and survival of the phlebotomines Nyssomyia intermedia and Nyssomyia neivai (both implicated as vectors of the cutaneous leishmaniasis agent) in an endemic area was investigated using a capture-mark-release technique in five experiments from August-December 2003 in municipality of Iporanga, state of São Paulo, Brazil. A total of 1,749 males and 1,262 females of Ny. intermedia and 915 males and 411 females of Ny. neivai were marked and released during the five experiments. Recapture attempts were made using automatic light traps, aspiration in natural resting places and domestic animal shelters and Shannon traps. A total of 153 specimens (3.48%) were recaptured: 2.59% (78/3,011) for Ny. intermedia and 5.35% (71/1,326) for Ny. neivai. Both species were recaptured up to 144 h post-release, with the larger part of them recaptured within 48 h. The median dispersion distances for Ny. intermedia and Ny. neivai, respectively, were 109 m and 100 m. The greatest dispersal range of Ny. intermedia was 180 m, while for Ny. neivai one female was recaptured in a pasture at 250 m and another in a pigsty at 520 m, showing a tendency to disperse to more open areas. The daily survival rates calculated based on regressions of the numbers of marked insects recaptured on the six successive days after release were 0.746 for males and 0.575 for females of Ny. intermedia and 0.649 for both sexes of Ny. neivai. The size of the populations in the five months ranged from 8,332-725,085 for Ny. intermedia males, 2,193-104,490 for Ny. intermedia females, 1,687-350,122 for Ny. neivai males and 254-49,705 for Ny. neivai females.

Survival, population pize, and gonotrophic cycle duration of Nyssomyia neivai (Diptera: Psychodidae) at an endemic area of American Cutaneous Leishmaniasis in Southeastern Brazil

The survival, absolute population size, gonotrophic cycle duration, and temporal and spatial abundance of Nyssomyia neivai (Pinto) were studied in a rural area endemic for American cutaneous leishmaniasis (ACL) in Conchal, Sõo Paulo State, southeastern Brazil, using mark-release-recapture techniques and by monitoring population fluctuation. The monthly abundance exhibited a unimodal pattern, with forest and domicile habitats having the highest relative abundances. A total of 1,873 males and 3,557 females were marked and released during the six experiments, of which 4.1-13.0 per cent of males and 4.1-11.8 per cent of females were recaptured. Daily survivorship estimated from the decline in recaptures per day was 0.681 for males and 0.667 for females. Gonotrophic cycle duration was estimated to be 4.0 d. Absolute population size was calculated using the Lincoln Index and ranged from 861 to 4,612 males and from 2,187 to 19,739 females. The low proportion of females that reach the age when they are potentially infective suggests that N. neivai has a low biological capacity to serve as a vector and that factors such as high biting rates and opportunistic feeding behavior would be needed to enable Leishmania (Viannia) braziliensis Vianna transmission. This agreed with the epidemiological pattern of ACL in southeastern Brazil that is characterized by low incidence, with isolated cases acquired principally within domiciliary habitats

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Merging Resource Availability with Isotope Mixing Models: The Role of Neutral Interaction Assumptions

Background: Bayesian mixing models have allowed for the inclusion of uncertainty and prior information in the analysis of trophic interactions using stable isotopes. Formulating prior distributions is relatively straightforward when incorporating dietary data. However, the use of data that are related, but not directly proportional, to diet (such as prey availability data) is often problematic because such information is not necessarily predictive of diet, and the information required to build a reliable prior distribution for all prey species is often unavailable. Omitting prey availability data impacts the estimation of a predator's diet and introduces the strong assumption of consumer ultrageneralism (where all prey are consumed in equal proportions), particularly when multiple prey have similar isotope values. Methodology: We develop a procedure to incorporate prey availability data into Bayesian mixing models conditional on the similarity of isotope values between two prey. If a pair of prey have similar isotope values (resulting in highly uncertain mixing model results), our model increases the weight of availability data in estimating the contribution of prey to a predator's diet. We test the utility of this method in an intertidal community against independently measured feeding rates. Conclusions: Our results indicate that our weighting procedure increases the accuracy by which consumer diets can be inferred in situations where multiple prey have similar isotope values. This suggests that the exchange of formalism for predictive power is merited, particularly when the relationship between prey availability and a predator's diet cannot be assumed for all species in a system.

Bayesian online algorithms for learning in discrete Hidden Markov Models

We propose and analyze two different Bayesian online algorithms for learning in discrete Hidden Markov Models and compare their performance with the already known Baldi-Chauvin Algorithm. Using the Kullback-Leibler divergence as a measure of generalization we draw learning curves in simplified situations for these algorithms and compare their performances.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

Discharge-calcium concentration relationships in streams of the Amazon and Cerrado of Brazil: soil or land use controlled

Stream discharge-concentration relationships are indicators of terrestrial ecosystem function. Throughout the Amazon and Cerrado regions of Brazil rapid changes in land use and land cover may be altering these hydrochemical relationships. The current analysis focuses on factors controlling the discharge-calcium (Ca) concentration relationship since previous research in these regions has demonstrated both positive and negative slopes in linear log(10)discharge-log(10)Ca concentration regressions. The objective of the current study was to evaluate factors controlling stream discharge-Ca concentration relationships including year, season, stream order, vegetation cover, land use, and soil classification. It was hypothesized that land use and soil class are the most critical attributes controlling discharge-Ca concentration relationships. A multilevel, linear regression approach was utilized with data from 28 streams throughout Brazil. These streams come from three distinct regions and varied broadly in watershed size (< 1 to > 10(6) ha) and discharge (10(-5.7)-10(3.2) m(3) s(-1)). Linear regressions of log(10)Ca versus log(10)discharge in 13 streams have a preponderance of negative slopes with only two streams having significant positive slopes. An ANOVA decomposition suggests the effect of discharge on Ca concentration is large but variable. Vegetation cover, which incorporates aspects of land use, explains the largest proportion of the variance in the effect of discharge on Ca followed by season and year. In contrast, stream order, land use, and soil class explain most of the variation in stream Ca concentration. In the current data set, soil class, which is related to lithology, has an important effect on Ca concentration but land use, likely through its effect on runoff concentration and hydrology, has a greater effect on discharge-concentration relationships.

Bayesian updating rules in continuous opinion dynamics models

Here, I investigate the use of Bayesian updating rules applied to modeling how social agents change their minds in the case of continuous opinion models. Given another agent statement about the continuous value of a variable, we will see that interesting dynamics emerge when an agent assigns a likelihood to that value that is a mixture of a Gaussian and a uniform distribution. This represents the idea that the other agent might have no idea about what is being talked about. The effect of updating only the first moments of the distribution will be studied, and we will see that this generates results similar to those of the bounded confidence models. On also updating the second moment, several different opinions always survive in the long run, as agents become more stubborn with time. However, depending on the probability of error and initial uncertainty, those opinions might be clustered around a central value.

Modeling associations between genetic markers using Bayesian networks

Motivation: Understanding the patterns of association between polymorphisms at different loci in a population ( linkage disequilibrium, LD) is of fundamental importance in various genetic studies. Many coefficients were proposed for measuring the degree of LD, but they provide only a static view of the current LD structure. Generative models (GMs) were proposed to go beyond these measures, giving not only a description of the actual LD structure but also a tool to help understanding the process that generated such structure. GMs based in coalescent theory have been the most appealing because they link LD to evolutionary factors. Nevertheless, the inference and parameter estimation of such models is still computationally challenging. Results: We present a more practical method to build GM that describe LD. The method is based on learning weighted Bayesian network structures from haplotype data, extracting equivalence structure classes and using them to model LD. The results obtained in public data from the HapMap database showed that the method is a promising tool for modeling LD. The associations represented by the learned models are correlated with the traditional measure of LD D`. The method was able to represent LD blocks found by standard tools. The granularity of the association blocks and the readability of the models can be controlled in the method. The results suggest that the causality information gained by our method can be useful to tell about the conservability of the genetic markers and to guide the selection of subset of representative markers.