Selection-mutation balance models with epistatic selection

We present an application of birth-and-death processes on configuration spaces to a generalized mutation4 selection balance model. The model describes the aging of population as a process of accumulation of mu5 tations in a genotype. A rigorous treatment demands that mutations correspond to points in abstract spaces. 6 Our model describes an infinite-population, infinite-sites model in continuum. The dynamical equation which 7 describes the system, is of Kimura-Maruyama type. The problem can be posed in terms of evolution of states 8 (differential equation) or, equivalently, represented in terms of Feynman-Kac formula. The questions of interest 9 are the existence of a solution, its asymptotic behavior, and properties of the limiting state. In the non-epistatic 10 case the problem was posed and solved in [Steinsaltz D., Evans S.N., Wachter K.W., Adv. Appl. Math., 2005, 11 35(1)]. In our model we consider a topological space X as the space of positions of mutations and the influence of epistatic potentials

Relative proportions of polycyclic aromatic hydrocarbons differ between accumulation bioassays and chemical methods to predict bioavailability

Chemical methods to predict the bioavailable fraction of organic contaminants are usually validated in the literature by comparison with established bioassays. A soil spiked with polycyclic aromatic hydrocarbons (PAHs) was aged over six months and subjected to butanol, cyclodextrin and tenax extractions as well as an exhaustive extraction to determine total PAH concentrations at several time points. Earthworm (Eisenia fetida) and rye grass root (Lolium multiflorum) accumulation bioassays were conducted in parallel. Butanol extractions gave the best relationship with earthworm accumulation (r2 ≤ 0.54, p ≤ 0.01); cyclodextrin, butanol and acetone–hexane extractions all gave good predictions of accumulation in rye grass roots (r2 ≤ 0.86, p ≤ 0.01). However, the profile of the PAHs extracted by the different chemical methods was significantly different (p < 0.01) to that accumulated in the organisms. Biota accumulated a higher proportion of the heavier 4-ringed PAHs. It is concluded that bioaccumulation is a complex process that cannot be predicted by measuring the bioavailable fraction alone. The ability of chemical methods to predict PAH accumulation in Eisenia fetida and Lolium multiflorum was hindered by the varied metabolic fate of the different PAHs within the organisms.

Genotypic variations in the accumulation of Cd, Cu, Pb and Zn exhibited by six commonly grown vegetables

Metal contaminants in garden and allotment soils could possibly affect human health through a variety of pathways. This study focused on the potential pathway of consumption of vegetables grown on contaminated soil. Five cultivars each of six common vegetables were grown in a control and in a soil spiked with Cd, Cu, Pb and Zn. Highly significant differences in metal content were evident between cultivars of a number of vegetables for several of the contaminants. Carrot and pea cultivars exhibited significant differences in accumulated concentrations of Cd and Cu with carrot cultivars also exhibiting significant differences in Zn. Distinctive differences were also identified when comparing one vegetable to another, legumes (Leguminosae) tending to be low accumulators, root vegetables (Umbelliferae and Liliaceae) tending to be moderate accumulators and leafy vegetables (Compositae and Chenopodiaceae) being high accumulators. (c) 2006 Elsevier Ltd. All rights reserved.

A new Snow-SVAT to simulate the accumulation and ablation of seasonal snow cover beneath a forest canopy

A new snow-soil-vegetation-atmosphere transfer (Snow-SVAT) scheme, which simulates the accumulation and ablation of the snow cover beneath a forest canopy, is presented. The model was formulated by coupling a canopy optical and thermal radiation model to a physically-based multi-layer snow model. This canopy radiation model is physically-based yet requires few parameters, so can be used when extensive in-situ field measurements are not available. Other forest effects such as the reduction of wind speed, interception of snow on the canopy and the deposition of litter were incorporated within this combined model, SNOWCAN, which was tested with data taken as part of the Boreal Ecosystem-Atmosphere Study (BOREAS) international collaborative experiment. Snow depths beneath four different canopy types and at an open site were simulated. Agreement between observed and simulated snow depths was generally good, with correlation coefficients ranging between r^2=0.94 and r^2=0.98 for all sites where automatic measurements were available. However, the simulated date of total snowpack ablation generally occurred later than the observed date. A comparison between simulated solar radiation and limited measurements of sub-canopy radiation at one site indicates that the model simulates the sub-canopy downwelling solar radiation early in the season to within measurement uncertainty.

The effects of triazole and strobilurin fungicide programmes on nitrogen uptake, partitioning, remobilization and grain N accumulation in winter wheat cultivars

Field experiments were conducted over 3 years to assess the effect of a triazole fungicide programme, and additions of strobilurin fungicides to it, on nitrogen uptake, accumulation and partitioning in a range of winter wheat cultivars. Commensurate with delayed senescence, fungicide programmes, particularly when including strobilurins, improved grain yield through improvements in both crop biomass and harvest index, although the relationship with green area duration of the flag leaf (GFLAD) depended on year and in some cases, cultivar. In all years fungicide treatments significantly increased the amount of nitrogen in the above-ground biomass, the amount of nitrogen in the grain and the nitrogen harvest index. All these effects could be linearly related to the fungicide effect on GFLAD. These relationships occasionally interacted with cultivar but there was no evidence that fungicide mode of action affected the relationship between GFLAD and yield of nitrogen in the grain. Fungicide treatments significantly reduced the amount of soil mineral N at harvest and when severe disease had been controlled, the net remobilization of N from the vegetation to the grain after anthesis. Fungicide maintained the filling of grain with both dry matter and nitrogen. The proportionate accumulation of nitrogen in the grain was later than that of dry matter and this difference was greater when fungicide had been applied. Effects of fungicide on grain protein concentration and its relationship with GFLAD were inconsistent over year and cultivar. There were several instances where grain protein concentration was unaffected despite large (1(.)5 t/ha) increases in grain yield following fungicide use. Dilution of grain protein concentration following fungicide use, when it did occur, was small compared with what would be predicted by adoption of other yield increasing techniques such as the selection of high yielding cultivars (based on currently available cultivars) or by growing wheat in favourable climates.

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.

Detection of unusual mutation within the VP1 region of different re-isolates of poliovirus Sabin vaccine

In the present study, a genomic analysis of full VP1 sequence region of 15 clinical re-isolates (14 healthy vaccinees and one bone marrow tumor patient) was conducted, aiming to the identification of mutations and to the assessment of their impact on virus fitness, providing also insights relevant with the natural evolution of Sabin strains. Clinical re-isolates were analyzed by RT-PCR, sequencing and computational analysis. Some re-isolates were characterized by an unusual mutational pattern in which non-synonymous mutations outnumbered the synonymous ones. Furthermore, the majority of amino-acid substitutions were located in the capsid exterior, specifically in N-Ags, near N-Ags and in the north rim of the canyon. Also mutations, which are well-known determinants of attenuation, were identified. The results of this study propose that some re-isolates are characterized by an evolutionary pattern in which non-synonymous mutations with a direct phenotypic impact on viral fitness are fixed in viral genomes, in spite of synonymous ones with no phenotypic impact on viral fitness. Results of the present retrospective characterization of Sabin clinical re-isolates, based on the full VP1 sequence, suggest that vaccine-derived viruses may make their way through narrow breaches and may evolve into transmissible pathogens even in adequately immunized populations. For this reason increased poliovirus laboratory surveillance should be permanent and full VP1 sequence analysis should be conducted even in isolates originating from healthy vaccinees.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

Seasonal changes in olive fruit characteristics and oil accumulation during ripening process

Olive fruits of three of the most important Spanish and Italian cultivars, 'Picual', `Hojiblanca' and 'Frantoio', were harvested at bi-weekly periods during three crop seasons to study their development and ripening process. Fresh and dry weights and ripening index were determined for fruit, while dry matter, oil and moisture contents were determined in both fruit and pulp (flesh). Fruit growth rate and oil accumulation were calculated. Each olive cultivar showed a different ripening pattern, 'Hojiblanca' being the last one to maturate. Fruit weight increased, decreasing its growth rate from the middle of November. Dry matter and moisture contents decreased during ripening in pulp and fruit, 'Hojiblanca' showing the highest values for both. Oil content, when expressed on a fresh weight basis, increased in all cultivars, although for the last time period showed variations due to climatic conditions. During ripening, oil content on a dry weight basis increased in fruit, but oil biosynthesis in flesh ceased from November. Olive fruits presented lower oil and higher dry matter contents in the year of lowest rainfall. Therefore fruit harvesting should be carried out from the middle of November in order to obtain the highest oil yield and avoid natural fruit drop. (C) 2004 Society of Chemical Industry.

An intronic growth hormone receptor mutation causing activation of a pseudoexon is associated with a broad spectrum of growth hormone insensitivity phenotypes

Context: Inherited GH insensitivity (GHI) is usually caused by mutations in the GH receptor (GHR). Patients present with short stature associated with high GH and low IGF-I levels and may have midfacial hypoplasia ( typical Laron syndrome facial features). We previously described four mildly affected GHI patients with an intronic mutation in the GHR gene (A.(1) -> G.(1) substitution in intron 6), resulting in the activation of a pseudoexon (6 Psi) and inclusion of 36 amino acids. Objective: The study aimed to analyze the clinical and genetic characteristics of additional GHI patients with the pseudoexon (6 Psi) mutation. Design/Patients: Auxological, biochemical, genetic, and haplotype data from seven patients with severe short stature and biochemical evidence of GHI were assessed. Main Outcome Measures: We assessed genotype-phenotype relationship. Results: One patient belongs to the same extended family, previously reported. She has normal facial features, and her IGF-I levels are in the low-normal range for age. The six unrelated patients, four of whom have typical Laron syndrome facial features, have heights ranging from -3.3 to -6.0 SD and IGF-I levels that vary from normal to undetectable. We hypothesize that the marked difference in biochemical and clinical phenotypes might be caused by variations in the splicing efficiency of the pseudoexon. Conclusions: Activation of the pseudoexon in the GHR gene can lead to a variety of GHI phenotypes. Therefore, screening for the presence of this mutation should be performed in all GHI patients without mutations in the coding exons.