Heat-shock response in Arabidopsis thaliana explored by multiplexed quantitative proteomics using differential metabolic labeling

We have developed a general method for multiplexed quantitative proteomics using differential metabolic stable isotope labeling and mass spectrometry. The method was successfully used to study the dynamics of heat-shock response in Arabidopsis thaliana. A number of known heat-shock proteins were confirmed, and some proteins not previously associated with heat shock were discovered. The method is applicable in stable isotope labeling and allows for high degrees of multiplexing.

Active-control trials with binary data: a comparison of methods for testing superiority or non-inferiority using the odds ratio

This paper considers methods for testing for superiority or non-inferiority in active-control trials with binary data, when the relative treatment effect is expressed as an odds ratio. Three asymptotic tests for the log-odds ratio based on the unconditional binary likelihood are presented, namely the likelihood ratio, Wald and score tests. All three tests can be implemented straightforwardly in standard statistical software packages, as can the corresponding confidence intervals. Simulations indicate that the three alternatives are similar in terms of the Type I error, with values close to the nominal level. However, when the non-inferiority margin becomes large, the score test slightly exceeds the nominal level. In general, the highest power is obtained from the score test, although all three tests are similar and the observed differences in power are not of practical importance. Copyright (C) 2007 John Wiley & Sons, Ltd.

ONeSAMP: a program to estimate effective population size using approximate Bayesian computation

The estimation of effective population size from one sample of genotypes has been problematic because most estimators have been proven imprecise or biased. We developed a web-based program, ONeSAMP that uses approximate Bayesian computation to estimate effective population size from a sample of microsatellite genotypes. ONeSAMP requires an input file of sampled individuals' microsatellite genotypes along with information about several sampling and biological parameters. ONeSAMP provides an estimate of effective population size, along with 95% credible limits. We illustrate the use of ONeSAMP with an example data set from a re-introduced population of ibex Capra ibex.

Fitting models for the joint action of two drugs using SAS(R)

Combinations of drugs are increasingly being used for a wide variety of diseases and conditions. A pre-clinical study may allow the investigation of the response at a large number of dose combinations. In determining the response to a drug combination, interest may lie in seeking evidence of synergism, in which the joint action is greater than the actions of the individual drugs, or of antagonism, in which it is less. Two well-known response surface models representing no interaction are Loewe additivity and Bliss independence, and Loewe or Bliss synergism or antagonism is defined relative to these. We illustrate an approach to fitting these models for the case in which the marginal single drug dose-response relationships are represented by four-parameter logistic curves with common upper and lower limits, and where the response variable is normally distributed with a common variance about the dose-response curve. When the dose-response curves are not parallel, the relative potency of the two drugs varies according to the magnitude of the desired effect and the models for Loewe additivity and synergism/antagonism cannot be explicitly expressed. We present an iterative approach to fitting these models without the assumption of parallel dose-response curves. A goodness-of-fit test based on residuals is also described. Implementation using the SAS NLIN procedure is illustrated using data from a pre-clinical study. Copyright © 2007 John Wiley & Sons, Ltd.

The best of both worlds: Producing national statistics using participatory methods

This paper presents our experience with combining statistical principles and participatory methods to generate national statistics. The methodology was developed in Malawi during 1999–2002. We demonstrate that if PRA is combined with statistical principles (including probability-based sampling and standardization), it can produce total population statistics and estimates of the proportion of households with certain characteristics (e.g., poverty). It can also provide quantitative data on complex issues of national importance such as poverty targeting. This approach is distinct from previous PRA-based approaches, which generate numbers at community level but only provide qualitative information at national level.

Probing the receptor interactions ofan H5 avian influenza virus using a baculovirus expression system and functionalised poly(acrylic acid) ligands

Influenza viruses attach to host cells by binding to terminal sialic acid (Neu5Ac) on glycoproteins or glycolipids. Both the linkage of Neu5Ac and the identity of other carbohydrates within the oligosaccharide are thought to play roles in restricting the host range of the virus. In this study, the receptor specificity of an H5 avian influenza virus haemagglutinin protein that has recently infected man (influenza strain A/Vietnam/1194/04) has been probed using carbohydrate functionalised poly(acrylic acid) polymers. A baculovirus expression system that allows facile and safe analysis of the Neu5Ac binding specificity of mutants of H5 HA engineered at sites that are predicted to effect a switch in host range has also been developed. (C) 2007 Elsevier Ltd. All rights reserved.

Modelling 2,4-dichlorophenol bioavailability and bioaccumulation by the freshwater fingernail clam Sphaerium corneum using artificial particles and humic acids

The complex and variable composition of natural sediments makes it very difficult to predict the bioavailability and bioaccumulation of sediment-bound contaminants. Several approaches have been proposed to overcome this problem, including an experimental model using artificial particles with or without humic acids as a source of organic matter. For this work, we have applied this experimental model, and also a sample of a natural sediment, to investigate the uptake and bioaccumulation of 2,4-dichlorophenol (2,4-DCP) by Sphaerium corneum. Additionally, the particle-water partition coefficients (K-d) were calculated. The results showed that the bioaccumulation of 2,4-DCP by clams did not depend solely on the levels of chemical dissolved, but also on the amount sorbed onto the particles and the characteristics and the strength of that binding. This study confirms the value of using artificial particles as a suitable experimental model for assessing the fate of sediment-bound contaminants. (c) 2006 Elsevier Ltd. All rights reserved.

Soil disturbed using a strip tillage implement on a range of soil types and the effects on sugar beet establishment

The area of soil disturbed using a single tine is well documented. However, modern strip tillage implements using a tine and disc design have not been assessed in the UK or in mainland Europe. Using a strip tillage implement has potential benefits for European agriculture where economic returns and sustainability are key issues. Using a strip tillage system a narrow zone is cultivated leaving most of the straw residue on the soil surface. Small field plot experiments were undertaken on three soil types and the operating parameters of forward speed, tine depth and tine design were investigated together with measurements of seedbed tilth and crop emergence. The type of tine used was found to be the primary factor in achieving the required volume of disturbance within a narrow zone whilst maintaining an area of undisturbed soil with straw residue on the surface. The winged tine produced greater disturbance at a given depth compared with the knife tine. Increasing forward speed did not consistently increase the volume of disturbance. In a sandy clay loam the tilth created and emergence of sugar beet by strip tillage and ploughing were similar but on a sandy loam the strip tillage treatments generally gave a finer tilth but poorer emergence particularly at greater working depth.

Quantitative proteomics using uniform 15N-labeling, mascot and the trans-proteomic pipeline

Stable isotope labeling combined with MS is a powerful method for measuring relative protein abundances, for instance, by differential metabolic labeling of some or all amino acids with 14N and 15N in cell culture or hydroponic media. These and most other types of quantitative proteomics experiments using high-throughput technologies, such as LC-MS/MS, generate large amounts of raw MS data. This data needs to be processed efficiently and automatically, from the mass spectrometer to statistically evaluated protein identifications and abundance ratios. This paper describes in detail an approach to the automated analysis of uniformly 14N/15N-labeled proteins using MASCOT peptide identification in conjunction with the trans-proteomic pipeline (TPP) and a few scripts to integrate the analysis workflow. Two large proteomic datasets from uniformly labeled Arabidopsis thaliana were used to illustrate the analysis pipeline. The pipeline can be fully automated and uses only common or freely available software.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Evidence for complex multigenic inheritance of radiation AML susceptibility in mice revealed using a surrogate phenotypic assay

The mapping of genes which affect individual cancer risk is an important but complex challenge. A surrogate assay of susceptibility to radiation-induced acute myeloid leukaemia (AML) in the mouse based on chromosomal radiosensitivity has been developed and validated. This assay was applied to the mapping of radiation-induced AML risk modifier loci by association with microsatellite markers. A region on chromosome (chr) 18 with strong association is identified and confirmed by backcross analysis. Additional loci on chrs 8 and 13 show significant association. A key candidate gene Rbbp8 on chr18 is identified. Rbbp8 is shown to be upregulated in response to X-irradiation in the AML sensitive CBA strain but not AML resistant C57BL/6 strain. This study demonstrates the strength of utilizing surrogate endpoints of cancer susceptibility in the mapping of mouse loci and identifies additional loci that may affect radiation cancer risk.

An extensible automated protein annotation tool: standardizing input and output using validated XML

Motivation: There is a frequent need to apply a large range of local or remote prediction and annotation tools to one or more sequences. We have created a tool able to dispatch one or more sequences to assorted services by defining a consistent XML format for data and annotations. Results: By analyzing annotation tools, we have determined that annotations can be described using one or more of the six forms of data: numeric or textual annotation of residues, domains (residue ranges) or whole sequences. With this in mind, XML DTDs have been designed to store the input and output of any server. Plug-in wrappers to a number of services have been written which are called from a master script. The resulting APATML is then formatted for display in HTML. Alternatively further tools may be written to perform post-analysis.

Using phylogeny to investigate the origins of the Cape flora: the importance of taxonomic, gene and genome sampling strategies

Phylogenetic methods hold great promise for the reconstruction of the transition from precursor to modern flora and the identification of underlying factors which drive the process. The phylogenetic methods presently used to address the question of the origin of the Cape flora of South Africa are considered here. The sampling requirements of each of these methods, which include dating of diversifications using calibrated molecular trees, sister pair comparisons, lineage through time plots and biogeographical optimizations are reviewed. Sampling of genes, genomes and species are considered. Although increased higher-level studies and increased sampling are required for robust interpretation, it is clear that much progress is already made. It is argued that despite the remarkable richness of the flora, the Cape flora is a valuable model system to demonstrate the utility of phylogenetic methods in determining the history of a modern flora.

Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.