Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

Isotopes, wool, and rangeland monitoring: Let the sheep do the sampling

Stable carbon isotope analyses of wool staples provided insight into the vegetation consumed by sheep at a temporal resolution not previously studied. Contemporary Australian and historic South African samples dating back to 1916 were analyzed for their stable carbon isotope ratio, a proxy for the proportion of C-3 and C-4 plant species consumed by animals. Sheep sample vegetation continuously throughout a year, and as their wool grows it integrates and stores information about their diet. In subtropical and tropical rangelands the majority of grass species are C-4. Since sheep prefer to graze, and their wool is an isotopic record of their diet, we now have the potential to develop a high resolution index to the availability of grass from a sheep's perspective. Isotopic analyses of wool suggest a new direction for monitoring grazing and for the reconstruction of past vegetation changes, which will make a significant contribution to traditional rangeland ecology and management. It is recommended that isotopic and other analyses of wool be further developed for use in rangeland monitoring programs to provide valuable feedback for land managers.

Bacteria in shrimp pond sediments: their role in mineralizing nutrients and some suggested sampling strategies

Strategies for sampling sediment bacteria were examined in intensive shrimp, Penaeus monodon (Fabricius), ponds in tropical Australia. Stratified sampling of bacteria at the end of the production season showed that the pond centre, containing flocculated sludge, had significantly higher bacterial counts (15.5 X 10(9) g(-1) dw) than the pond periphery (8.1 X 10(9) g(-1) dw), where the action of aerators had swept the pond floor. The variation in bacterial counts between these two zones within a pond was higher than that between sites within each zone or between ponds. Therefore, sampling effort should be focused within these zones: for example, sampling two ponds at six locations within each of the two zones resulted in a coefficient of variation of approximate to 5%. Bacterial numbers in the sediment were highly correlated with sediment grain size, probably because eroded soil particles and organic waste both accumulated in the centre of the pond. Despite high inputs of organic matter added to the ponds, principally as pelleted feeds, the mean bacterial numbers and nutrient concentrations (i.e. organic carbon, nitrogen and phosphorus) in the sediment were similar to those found in mangrove sediments. This suggests that bacteria are rapidly remineralizing particulates into soluble compounds. Bacterial numbers were highly correlated with organic carbon and total kjeldahl nitrogen in the sediment, suggesting that these were limiting factors to bacterial growth.

Strong genetic structuring in a habitat specialist, the Oxleyan Pygmy Perch Nannoperca oxleyana

This study used allozyme and mitochondrial DNA variation to examine genetic structure in the Oxleyan Pygmy Perch Nannoperca oxleyana. This small-bodied freshwater fish has a very restricted distribution occurring only in some small coastal streams in south-east Queensland and northern New South Wales. It was expected that subpopulations may contain little genetic variation and be highly differentiated from one another. The results, based on allozyme and mitochondrial DNA control region variation were in agreement with these expectations. Allozyme variation was very low overall, with only one locus showing variation at most sites. The high differentiation was because a different locus tended to be polymorphic at each site. Mitochondrial variation within sites was also low, but some sites had unique haplotypes. The patterns of similarity among mitochondrial DNA haplotypes were not as expected from geographical proximity alone. In particular, although some northern sites had unique haplotypes, four sites spread along 200 km of coastline were remarkably similar, sharing the same common haplotype at similar frequencies. We suggest that these four streams may have had a confluence relatively recently, possibly when sea levels were lower, 8000-10 000 BP.

Semiochemicals and social signaling in the wild European rabbit in Australia: II. Variations in chemical composition of chin gland secretion across sampling sites

The volatile components of the chin gland secretion of the wild European rabbit, Oryctolagus cuniculus (L.), were investigated with the use of gas chromatography. Studies of the chemical nature of this secretion by previous workers demonstrated that it was important in the maintenance of social structure in this species. This study identified 34 different volatile components that consist primarily of aromatic and aliphatic hydrocarbons. Especially common are a series of alkyl-substituted benzene derivatives that provide most of the compound diversity in the secretion. Samples of chin gland secretion collected from animals at three different geographical locations, separated by more than 100 km, showed significant differences in composition. This work suggests that variation among populations needs to be considered when undertaking semiochemical research. Alternate nonparametric methods are also used for the analysis of chromatographic data.

The fMRI response of the LGN and V1 to cardinal red-green, blue-yellow, and achromatic visual stimuli

We compared the responsiveness of the LGN and the early retinotopic cortical areas to stimulation of the two cone-opponent systems (red - green and blue - yellow) and the achromatic system. This was done at two contrast levels to control for any effect of contrast. MR images were acquired on seven subjects with a 4T Bruker MedSpec scanner. The early visual cortical areas were localised by phase encoded retinotopic mapping with a volumetric analysis (Dumoulin et al, 2003 NeuroImage 18 576 - 587). We initially located the LGN in four subjects by using flickering stimuli in a separate scanning session, but subsequently identified it using the experimental stimuli. Experimental stimuli were sine-wave counterphasing rings (2 Hz, 0.5 cycle deg-1), cardinal for the selective activation of the L/M cone-opponent (RG), S cone-opponent (BY), and achromatic (Ach) systems. A region of interest analysis was performed. When presented at equivalent absolute contrasts (cone contrast = 5% - 6%), the BOLD response of the LGN is strongest to isoluminant red - green stimuli and weakest to blue - yellow stimuli, with the achromatic response falling in between. Area V1, on the other hand, responds best to both chromatic stimuli, with the achromatic response falling below. The key change from the LGN to V1 is a dramatic boost in the relative blue - yellow response, which occurred at both contrast levels used. This greatly enhanced cortical response to blue - yellow relative to the red - green and achromatic responses may be due to an increase in cell number and/or cell response between the LGN and V1. We speculate that the effect might reflect the operation of contrast constancy across colour mechanisms at the cortical level.

Local dispersion of the Eucalyptus leaf beetle Chrysophtharta bimaculata (Coleoptera: Chrysomelidae), and implications for forest protection

1. Chrysophtharta bimaculata is a native chrysomelid species that can cause chronic defoliation of plantation and regrowth Eucalyptus forests in Tasmania, Australia. Knowledge of the dispersion pattern of C. bimaculata was needed in order to assess the efficiency of an integrated pest management (IPM) programme currently used for its control. 2. Using data from yellow flight traps, local populations of C. bimaculata adults were monitored over a season at spatial scales relevant to commercial forestry: within a 50-ha operational management unit (a forestry 'coupe') and between coupes. In addition, oviposition was monitored over a season at a subset of the between-coupe sites. 3. Dispersion indices (Taylor's Power Law and Iwao's Mean Crowding regression method) demonstrated that C. bimaculata adults were spatially aggregated within and between coupes, although the number of egg-batches laid at the between-coupe scale was uniform. Spatial autocorrelation analysis showed that trap-catches at the within-coupe level were similar (positively autocorrelated) to a radius distance of approximately 110 m, and then dissimilar (negatively autocorrelated) at approximately 250 m. At the between-coupe scale, no repeatable spatial autocorrelation patterns were observed. 4. For any individual site, rapid changes in beetle density were observed to be associated with loosely aggregated flights of beetles into and out of that site. Peak adult catches (> the weekly mean plus standard deviation trap-catch) for a site occurred for a period of 2.0 +/- 0.22 weeks at a time (n = 37), with normally only one or two peaks per site per season. Peak oviposition events for a site occurred on average 1.4 +/- 0.11 times per season and lasted 1.5 +/- 0.12 weeks. 5. Analysis of an extensive data set (n = 417) demonstrated that adult abundance at a site was positively correlated with egg density, but negatively correlated with tree damage (caused by conspecifics) and the presence of conspecific larvae. There was no relationship between adult abundance and a visual estimate of the amount of young foliage on trees. 6. Adults of C. bimaculata are show n to occur in relatively small, mobile aggregations. This means that pest surveys must be both regular (less than 2 weeks apart) and intensive (with sampling points no more than 150 m apart) if beetle populations are to be monitored with confidence. Further refinement of the current IPM strategy must recognize the problems posed by this temporal and spatial patchiness, particularly with regard to the use of biological insecticides, such as Bacillus thuringiensis, for which only a very short operational window exists.

A Comparison of Bayesian and Sampling Theory Inferences in a Probit Model

HE PROBIT MODEL IS A POPULAR DEVICE for explaining binary choice decisions in econometrics. It has been used to describe choices such as labor force participation, travel mode, home ownership, and type of education. These and many more examples can be found in papers by Amemiya (1981) and Maddala (1983). Given the contribution of economics towards explaining such choices, and given the nature of data that are collected, prior information on the relationship between a choice probability and several explanatory variables frequently exists. Bayesian inference is a convenient vehicle for including such prior information. Given the increasing popularity of Bayesian inference it is useful to ask whether inferences from a probit model are sensitive to a choice between Bayesian and sampling theory techniques. Of interest is the sensitivity of inference on coefficients, probabilities, and elasticities. We consider these issues in a model designed to explain choice between fixed and variable interest rate mortgages. Two Bayesian priors are employed: a uniform prior on the coefficients, designed to be noninformative for the coefficients, and an inequality restricted prior on the signs of the coefficients. We often know, a priori, whether increasing the value of a particular explanatory variable will have a positive or negative effect on a choice probability. This knowledge can be captured by using a prior probability density function (pdf) that is truncated to be positive or negative. Thus, three sets of results are compared:those from maximum likelihood (ML) estimation, those from Bayesian estimation with an unrestricted uniform prior on the coefficients, and those from Bayesian estimation with a uniform prior truncated to accommodate inequality restrictions on the coefficients.

Effects of sampling opportunities on preference for adults with severe disabilities

A variable that appears to affect preference development is the exposure to a variety of options. Providing opportunities for systematically sampling different options is one procedure that can facilitate the development of preference, which is indicated by the consistency of selections. The purpose of this study was to evaluate the effects of providing sampling opportunities on the preference development for two adults with severe disabilities. Opportunities for sampling a variety of drink items were presented, followed by choice opportunities for selections at the site where sampling occurred and at a non-sampling site (a grocery store). Results show that the participants developed a definite response consistency in selections at both sites. Implications for sampling practices are discussed.

Low-density koala (Phascolarctos cinereus) populations in the mulgalands of south-west Queensland. I. Faecal pellet sampling protocol

Koala (Phascolarctos cinereus) populations in eastern Australia are threatened by land clearing for agricultural and urban development. At the same time, conservation efforts are hindered by a dearth of information about inland populations. Faecal deposits offer a source of information that is readily available and easily collected non-invasively. We detail a faecal pellet sampling protocol that was developed for use in a large rangeland biogeographic region. The method samples trees in belt transects, uses a thorough search at the tree base to quickly identify trees with koala pellets under them, then estimates the abundance of faecal pellets under those trees using 1-m(2) quadrats. There was a strong linear relationship between these estimates and a complete enumeration of pellet abundance under the same trees. We evaluated the accuracy of our method in detecting trees where pellets were present by means of a misclassification index that was weighed more heavily for missed trees that had high numbers of pellets under them. This showed acceptable accuracy in all landforms except riverine, where some trees with large numbers of pellets were missed. Here, accuracy in detecting pellet presence was improved by sampling with quadrats, rather than basal searches. Finally, we developed a method to reliably age pellets and demonstrate how this protocol could be used with the faecal-standing-crop method to derive a regional estimate of absolute koala abundance.