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.

Exploring sequential data: Commentary on Bowers, Jentsch, Salas, and Braun (1998)

In an investigation intended to determine training needs of night crews, Bowers et al. (1998, this issue) report two studies showing that the patterning of communication is a better discriminator of good and poor crews than is the content of communication. Bowers et al. characterize their studies as intended to generate hypotheses for training needs and draw connections with Exploratory Sequential Data Analysis (ESDA). Although applauding the intentions of Bowers ct al., we point out some concerns with their characterization and implementation of ESDA. Our principal concern is that the Bowers et al. exploration of the data does not convincingly lead them back to a better fundamental understanding of the original phenomena they are investigating.

Characterization of the sequential non-covalent and covalent interactions of the antitumour antibiotic hedamycin with double stranded DNA by NMR spectroscopy

Hedamycin, a member of the pluramycin class of antitumour antibiotics, consists of a planar anthrapyrantrione chromophore to which is attached two aminosugar rings at one end and a bisepoxide-containing sidechain at the other end, Binding to double-stranded DNA is known to involve both reversible and non-reversible modes of interaction. As a part of studies directed towards elucidating the structural basis for the observed 5'-pyGT-3' sequence selectivity of hedamycin, we conducted one-dimensional NMR titration experiments at low temperature using the hexadeoxyribonucleotide duplexes d(CACGTG)(2) and d(CGTACG)(2). Spectral changes which occurred during these titrations are consistent with hedamycin initially forming a reversible complex in slow exchange on the NMR timescale and binding through intercalation of the chromophore. Monitoring of this reversible complex over a period of hours revealed a second type of spectral change which corresponds with formation of a non-reversible complex. Copyright (C) 1999 John Wiley & Sons, Ltd.

Sequential increase in Egr-1 and AP-1 DNA binding activity in the dentate gyrus following the induction of long-term potentiation

Establishment of long-term potentiation (LTP) at perforant path synapses is highly correlated with increased expression of Egr and AP-1 transcription factors in rat dentate gyrus granule cells. We have investigated whether increased transcription factor levels are reflected in increased transcription factor activity by assessing Egr and AP-I DNA binding activity using gel shift assays. LTP produced an increase in binding to the Egr element, which was NMDA receptor-dependent and correlated closely with our previously reported increase in Egr-1 (zif/268) protein levels. Supershift analysis confirmed involvement of Egr-1, but not Egr-2 in the DNA binding activity. AP-1 DNA binding was also rapidly elevated in parallel with protein levels, however, the peak increase in activity was delayed until 4 h, a time point when we have previously shown that only jun-D protein was elevated. These data indicate that binding of Egr-1 and AP-1 to their response elements is increased in two phases. This may result in activation of distinct banks of target genes which contribute to the establishment of persistent LTP. (C) 2000 Elsevier Science B.V. All rights reserved.

A simplified, sequential, phosphorus fractionation method

Hedley er al. (1982) developed what has become the most widely used land modified), phosphorus (P) fractionation technique. It consists of sequential extraction of increasingly less phytoavailable P pools. Extracts are centrifuged at up to 25000 g (RCF) and filtered to 0.45 mu m to ensure that soil is not lost between extractions. In attempting to transfer this method to laboratories with limited facilities, it was considered that access to high-speed centrifuges, and the cost of frequent filtration may prevent adoption of this P fractionation technique. The modified method presented here was developed to simplify methodology, reduce cost, and therefore increase accessibility of P fractionation technology. It provides quantitative recovery of soil between extractions, using low speed centrifugation without filtration. This is achieved by increasing the ionic strength of dilute extracts, through the addition of NaCl, to flocculate clay particles. Addition of NaCl does not change the amount of P extracted. Flocculation with low speed centrifugation produced extracts comparable with those having undergone filtration (0.025 mu m). A malachite green colorimetric method was adopted for inorganic P determination, as this simple manual method provides high sensitivity with negligible interference from other anions. This approach can also be used for total P following digestion, alternatively non-discriminatory methods, such as inductively coupled plasma atomic emission spectroscopy, may be employed.

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 preparation and readiness for voluntary movement: a high-field event-related fMRI study of the Bereitschafts-BOLD response

Activity within motor areas of the cortex begins to increase 1 to 2 s prior to voluntary self-initiated movement (termed the Bereitschaftspotential or readiness potential). There has been much speculation and debate over the precise source of this early premovement activity as it is important for understanding the roles of higher order motor areas in the preparation and readiness for voluntary movement. In this study, we use high-field (3-T) event-related fMRI with high temporal sampling (partial brain volumes every 250 ms) to specifically examine hemodynamic response time courses during the preparation, readiness, and execution of purely self-initiated voluntary movement. Five right-handed healthy volunteers performed a rapid sequential finger-to-thumb movement performed at self-determined times (12-15 trials). Functional images for each trial were temporally aligned and the averaged time series for each subject was iteratively correlated with a canonical hemodynamic response function progressively shifted in time. This analysis method identified areas of activation without constraining hemodynamic response timing. All subjects showed activation within frontal mesial areas, including supplementary motor area (SMA) and cingulate motor areas, as well as activation in left primary sensorimotor areas. The time courses of hemodynamic responses showed a great deal of variability in shape and timing between subjects; however, four subjects clearly showed earlier relative hemodynamic responses within SMA/cingulate motor areas compared with left primary motor areas. These results provide further evidence that the SMA and cingulate motor areas are major contributors to early stage premovement activity and play an important role in the preparation and readiness for voluntary movement. (C) 2003 Elsevier Inc. All rights reserved.

Sequential asymmetric auctions with endogenous participation

In this paper we suggest a model of sequential auctions with endogenous participation where each bidder conjectures about the number of participants at each round. Then, after learning his value, each bidder decides whether or not to participate in the auction. In the calculation of his expected value, each bidder uses his conjectures about the number of participants for each possible subgroup. In equilibrium, the conjectured probability is compatible with the probability of staying in the auction. In our model, players face participation costs, bidders may buy as many objects as they wish and they are allowed to drop out at any round. Bidders can drop out at any time, but they cannot come back to the auction. In particular we can determine the number of participants and expected prices in equilibrium. We show that for any bidding strategy, there exists such a probability of staying in the auction. For the case of stochastically independent objects, we show that in equilibrium every bidder who decides to continue submits a bid that is equal to his value at each round. When objects are stochastically identical, we are able to show that expected prices are decreasing.

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.