Microbial impacts on plant-herbivore interactions: the indirect effects of a birch pathogen on a birch aphid

The role of indirect interactions in structuring communities is becoming increasingly recognised. Plant fungi can bring about changes in plant chemistry which may affect insect herbivores that share the same plant, and hence the two may interact indirectly. This study investigated the indirect effects of a fungal pathogen (Marssonina betulae) of silver birch (Betula pendula) on an aphid (Euceraphis betulae), and the processes underpinning the interaction. There was a strong positive association between natural populations of the aphid and leaves bearing high fungal infection. In choice tests, significantly more aphids settled on leaves inoculated with the fungus than on asymptomatic leaves. Individual aphids reared on inoculated leaves were heavier, possessed longer hind tibiae and displayed enhanced embryo development compared with aphids reared on asymptomatic leaves; population growth rate was also positively correlated with fungal infection when groups of aphids were reared on inoculated branches. Changes in leaf chemistry were associated with fungal infection with inoculated leaves containing higher concentrations of free-amino acids. This may reflect a plant-initiated response to fungal attack in which free amino acids from the degradation of mesophyll cells are translocated out of infected leaves via the phloem. These changes in plant chemistry are similar to those occurring during leaf senescence, and are proposed as the mechanistic basis for the positive interaction between the fungus and aphid.

Indirect population dynamic benefits of altered life-history trade-offs in response to egg harvesting

Variations in demographic rates due to differential resource allocation between individuals are important considerations in the development of accurate population dynamic models. Systematic harvesting can alter age structure and/or reduce population density, conferring indirect positive benefits on the source population as a result of a consequent redistribution of resources between the remaining individuals. Independently of effects mediated through changes in density and competition, demographic rates can also be influenced by within-individual competition for resources. Harvesting dependent life stages can reduce an individual's current reproductive costs, allowing increased investment in its future fecundity and survival. Although such changes in demographic rates are well known, there has been little exploration of the potential impact on population dynamics. We use empirical data collected from a successfully reintroduced population of the Mauritius kestrel Falco punctatus to explore the population consequences of manipulating reproductive effort through harvesting. Consequent increases in an individual's future fecundity and survival allow source populations to withstand longer and more intensive harvesting regimes without being exposed to an increase in extinction risk, increasing maximum sustainable yields. These effects may also buffer populations against the impacts of stochastic events, but directional shifts in environmental conditions that increase reproductive costs may have detrimental population-level effects.

Direct versus indirect sexual selection: genetic basis of colour, size and recruitment in a wild bird

Indirect and direct models of sexual selection make different predictions regarding the quantitative genetic relationships between sexual ornaments and fitness. Indirect models predict that ornaments should have a high heritability and that strong positive genetic covariance should exist between fitness and the ornament. Direct models, on the other hand, make no such assumptions about the level of genetic variance in fitness and the ornament, and are therefore likely to be more important when environmental sources of variation are large. Here we test these predictions in a wild population of the blue tit (Parus caeruleus), a species in which plumage coloration has been shown to be under sexual selection. Using 3 years of cross-fostering data from over 250 breeding attempts, we partition the covariance between parental coloration and aspects of nestling fitness into a genetic and environmental component. Contrary to indirect models of sexual selection, but in agreement with direct models, we show that variation in coloration is only weakly heritable (h(2) < 0.11), and that two components of offspring fitness-nestling size and fledgling recruitment-are strongly dependent on parental effects, rather than genetic effects. Furthermore, there was no evidence of significant positive genetic covariation between parental colour and offspring traits. Contrary to direct benefit models, however, we find little evidence that variation in colour reliably indicates the level of parental care provided by either males or females. Taken together, these results indicate that the assumptions of indirect models of sexual selection are not supported by the genetic basis of the traits reported on here.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

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.

Bayesian analysis of correlated evolution of discrete characters by reversible-jump Markov chain Monte Carlo

We describe a Bayesian method for investigating correlated evolution of discrete binary traits on phylogenetic trees. The method fits a continuous-time Markov model to a pair of traits, seeking the best fitting models that describe their joint evolution on a phylogeny. We employ the methodology of reversible-jump ( RJ) Markov chain Monte Carlo to search among the large number of possible models, some of which conform to independent evolution of the two traits, others to correlated evolution. The RJ Markov chain visits these models in proportion to their posterior probabilities, thereby directly estimating the support for the hypothesis of correlated evolution. In addition, the RJ Markov chain simultaneously estimates the posterior distributions of the rate parameters of the model of trait evolution. These posterior distributions can be used to test among alternative evolutionary scenarios to explain the observed data. All results are integrated over a sample of phylogenetic trees to account for phylogenetic uncertainty. We implement the method in a program called RJ Discrete and illustrate it by analyzing the question of whether mating system and advertisement of estrus by females have coevolved in the Old World monkeys and great apes.

A pressure-jump time-resolved x-ray diffraction study of cubic-cubic transition kinetics in monoolein

In the past two decades, the geometric pathways involved in the transformations between inverse bicontinuous cubic phases in amphiphilic systems have been extensively theoretically modeled. However, little experimental data exists on the cubic-cubic transformation in pure lipid systems. We have used pressure-jump time-resolved X-ray diffraction to investigate the transition between the gyroid Q(II)(G) and double-diamond Q(II)(D) phases in mixtures of 1-monoolein in 30 wt% water. We find for this system that the cubic-cubic transition occurs without any detectable intermediate structures. In addition, we have determined the kinetics of the transition, in both the forward and reverse directions, as a function of pressure-jump amplitude, temperature, and water content. A recently developed model allows (at least in principle) the calculation of the activation energy for lipid phase transitions from such data. The analysis is applicable only if kinetic reproducibility is achieved, at least within one sample, and achievement of such kinetic reproducibility is shown here, by carrying out prolonged pressure-cycling. The rate of transformation shows clear and consistent trends with pressure-jump amplitude, temperature, and water content, all of which are shown to be in agreement with the effect of the shift in the position of the cubic-cubic phase boundary following a change in the thermodynamic parameters.

Pressure-jump X-ray studies of liquid crystal transitions in lipids

In this paper, we give an overview of our studies by static and time-resolved X-ray diffraction of inverse cubic phases and phase transitions in lipids. In 1, we briefly discuss the lyotropic phase behaviour of lipids, focusing attention on non-lamellar structures, and their geometric/topological relationship to fusion processes in lipid membranes. Possible pathways for transitions between different cubic phases are also outlined. In 2, we discuss the effects of hydrostatic pressure on lipid membranes and lipid phase transitions, and describe how the parameters required to predict the pressure dependence of lipid phase transition temperatures can be conveniently measured. We review some earlier results of inverse bicontinuous cubic phases from our laboratory, showing effects such as pressure-induced formation and swelling. In 3, we describe the technique of pressure-jump synchrotron X-ray diffraction. We present results that have been obtained from the lipid system 1:2 dilauroylphosphatidylcholine/lauric acid for cubic-inverse hexagonal, cubic-cubic and lamellar-cubic transitions. The rate of transition was found to increase with the amplitude of the pressure-jump and with increasing temperature. Evidence for intermediate structures occurring transiently during the transitions was also obtained. In 4, we describe an IDL-based 'AXCESS' software package being developed in our laboratory to permit batch processing and analysis of the large X-ray datasets produced by pressure-jump synchrotron experiments. In 5, we present some recent results on the fluid lamellar-Pn3m cubic phase transition of the single-chain lipid 1-monoelaidin, which we have studied both by pressure-jump and temperature-jump X-ray diffraction. Finally, in 6, we give a few indicators of future directions of this research. We anticipate that the most useful technical advance will be the development of pressure-jump apparatus on the microsecond time-scale, which will involve the use of a stack of piezoelectric pressure actuators. The pressure-jump technique is not restricted to lipid phase transitions, but can be used to study a wide range of soft matter transitions, ranging from protein unfolding and DNA unwinding and transitions, to phase transitions in thermotropic liquid crystals, surfactants and block copolymers.

Analyzing diffusion tensor images with ghosting artifacts: the effects of direct and indirect normalization

The current study aims to assess the applicability of direct or indirect normalization for the analysis of fractional anisotropy (FA) maps in the context of diffusion-weighted images (DWIs) contaminated by ghosting artifacts. We found that FA maps acquired by direct normalization showed generally higher anisotropy than indirect normalization, and the disparities were aggravated by the presence of ghosting artifacts in DWIs. The voxel-wise statistical comparisons demonstrated that indirect normalization reduced the influence of artifacts and enhanced the sensitivity of detecting anisotropy differences between groups. This suggested that images contaminated with ghosting artifacts can be sensibly analyzed using indirect normalization.

Evidence that indirect inhibition of saccade initiation improves saccade accuracy

Saccadic eye-movements to a visual target are less accurate if there are distracters close to its location (local distracters). The addition of more distracters, remote from the target location (remote distracters), invokes an involuntary increase in the response latency of the saccade and attenuates the effect of local distracters on accuracy. This may be due to the target and distracters directly competing (direct route) or to the remote distracters acting to impair the ability to disengage from fixation (indirect route). To distinguish between these we examined the development of saccade competition by recording saccade latency and accuracy responses made to a target and local distracter compared with those made with an addition of a remote distracter. The direct route would predict that the remote distracter impacts on the developing competition between target and local distracter, while the indirect route would predict no change as the accuracy benefit here derives from accessing the same competitive process but at a later stage. We found that the presence of the remote distracter did not change the pattern of accuracy improvement. This suggests that the remote distracter was acting along an indirect route that inhibits disengagement from fixation, slows saccade initiation, and enables more accurate saccades to be made.

Kinetic studies of reactions of the nitrate radical (NO3) with peroxy radicals (RO2): an indirect source of OH at night?

A discharge-flow system, coupled to cavity-enhanced absorption spectroscopy (CEAS) detection systems for NO3 at lambda = 662 nm and NO2 at lambda = 404 nm, was used to investigate the kinetics of the reactions of NO3 with eight peroxy radicals at P similar to 5 Torr and T similar to 295 K. Values of the rate constants obtained were (k/10(-12) cm(3) molecule(-1) s(-1)): CH3O2 (1.1 +/- 0.5), C2H5O2 (2.3 +/- 0.7), CH2FO2 (1.4 +/- 0.9), CH2ClO2 (3.8(-2.6)(+1.4)), c-C5H9O2 (1.2(-0.5)(+1.1)), c-C6H11O2 (1.9 +/- 0.7), CF3O2 (0.62 +/- 0.17) and CF3CFO2CF3 (0.24 +/- 0.13). We explore possible relationships between k and the orbital energies of the reactants. We also provide a brief discussion of the potential impact of the reactions of NO3 with RO2 on the chemistry of the night-time atmosphere.

The indirect global warming potential and global temperature change potential due to methane oxidation

Methane is the second most important anthropogenic greenhouse gas in the atmosphere next to carbon dioxide. Its global warming potential (GWP) for a time horizon of 100 years is 25, which makes it an attractive target for climate mitigation policies. Although the methane GWP traditionally includes the methane indirect effects on the concentrations of ozone and stratospheric water vapour, it does not take into account the production of carbon dioxide from methane oxidation. We argue here that this CO2-induced effect should be included for fossil sources of methane, which results in slightly larger GWP values for all time horizons. If the global temperature change potential is used as an alternative climate metric, then the impact of the CO2-induced effect is proportionally much larger. We also discuss what the correction term should be for methane from anthropogenic biogenic sources.

An analysis of the impact of R&D on productivity using Bayesian model averaging with a reversible jump algorithm.

A Bayesian Model Averaging approach to the estimation of lag structures is introduced, and applied to assess the impact of R&D on agricultural productivity in the US from 1889 to 1990. Lag and structural break coefficients are estimated using a reversible jump algorithm that traverses the model space. In addition to producing estimates and standard deviations for the coe¢ cients, the probability that a given lag (or break) enters the model is estimated. The approach is extended to select models populated with Gamma distributed lags of di¤erent frequencies. Results are consistent with the hypothesis that R&D positively drives productivity. Gamma lags are found to retain their usefulness in imposing a plausible structure on lag coe¢ cients, and their role is enhanced through the use of model averaging.