Inferring population history with DIY ABC: a user-friendly approach to approximate Bayesian computation

Genetic data obtained on population samples convey information about their evolutionary history. Inference methods can extract part of this information but they require sophisticated statistical techniques that have been made available to the biologist community (through computer programs) only for simple and standard situations typically involving a small number of samples. We propose here a computer program (DIY ABC) for inference based on approximate Bayesian computation (ABC), in which scenarios can be customized by the user to fit many complex situations involving any number of populations and samples. Such scenarios involve any combination of population divergences, admixtures and population size changes. DIY ABC can be used to compare competing scenarios, estimate parameters for one or more scenarios and compute bias and precision measures for a given scenario and known values of parameters (the current version applies to unlinked microsatellite data). This article describes key methods used in the program and provides its main features. The analysis of one simulated and one real dataset, both with complex evolutionary scenarios, illustrates the main possibilities of DIY ABC.

An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.

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.

Inhibition of coxsackie B virus infection by soluble forms of its receptors: Binding affinities, altered particle formation, and competition with cellular receptors

We previously reported that soluble decay-accelerating factor (DAF) and coxsackievirus-adenovirus receptor (CAR) blocked coxsackievirus 133 (CVB3) myocarditis in mice, but only soluble CAR blocked CVB3-mediated pancreatitis. Here, we report that the in vitro mechanisms of viral inhibition by these soluble receptors also differ. Soluble DAF inhibited virus infection through the formation of reversible complexes with CVB3, while binding of soluble CAR to CVB induced the formation of altered (A) particles with a resultant irreversible loss of infectivity. A-particle formation was characterized by loss of VP4 from the virions and required incubation of CVB3-CAR complexes at 37 degrees C. Dimeric soluble DAF (DAF-Fc) was found to be 125-fold-more effective at inhibiting CVB3 than monomeric DAF, which corresponded to a 100-fold increase in binding affinity as determined by surface plasmon resonance analysis. Soluble CAR and soluble dimeric CAR (CAR-Fc) bound to CVB3 with 5,000- and 10,000-fold-higher affinities than the equivalent forms of DAF. While DAF-Fc was 125-fold-more effective at inhibiting virus than monomeric DAF, complement regulation by DAF-Fc was decreased 4 fold. Therefore, while the virus binding was a cooperative event, complement regulation was hindered by the molecular orientation of DAF-Fc, indicating that the regions responsible for complement regulation and virus binding do not completely overlap. Relative contributions of CVB binding affinity, receptor binding footprint on the virus capsid, and induction of capsid conformation alterations for the ability of cellular DAF and CAR to act as receptors are discussed.

Comparative evaluation of a new effective population size estimator based on approximate Bayesian computation

We describe and evaluate a new estimator of the effective population size (N-e), a critical parameter in evolutionary and conservation biology. This new "SummStat" N-e. estimator is based upon the use of summary statistics in an approximate Bayesian computation framework to infer N-e. Simulations of a Wright-Fisher population with known N-e show that the SummStat estimator is useful across a realistic range of individuals and loci sampled, generations between samples, and N-e values. We also address the paucity of information about the relative performance of N-e estimators by comparing the SUMMStat estimator to two recently developed likelihood-based estimators and a traditional moment-based estimator. The SummStat estimator is the least biased of the four estimators compared. In 32 of 36 parameter combinations investigated rising initial allele frequencies drawn from a Dirichlet distribution, it has the lowest bias. The relative mean square error (RMSE) of the SummStat estimator was generally intermediate to the others. All of the estimators had RMSE > 1 when small samples (n = 20, five loci) were collected a generation apart. In contrast, when samples were separated by three or more generations and Ne less than or equal to 50, the SummStat and likelihood-based estimators all had greatly reduced RMSE. Under the conditions simulated, SummStat confidence intervals were more conservative than the likelihood-based estimators and more likely to include true N-e. The greatest strength of the SummStat estimator is its flexible structure. This flexibility allows it to incorporate any, potentially informative summary statistic from Population genetic data.

A block in virus-like particle maturation following assembly of murine leukaemia virus in insect cells

Expression of the murine leukaemia virus (MLV) major Gag antigen p65(Gag) using the baculovirus expression system leads to efficient assembly and release of virus-like particles (VLP) representative of immature MLV. Expression of P180(Gag-Pol), facilitated normally in mammalian cells by readthrough of the p65(Gag) termination codon, also occurs efficiently in insect cells to provide a source of the MLV protease and a pattern of p65(Gag) processing similar to that observed in mammalian cells. VLP release from P180(Gag-Pol) expressing cells however remains essentially immature with disproportionate levels of the uncleaved p65(Gag) precursor when compared to the intracellular Gag profile. Changing the p65(Gag) termination codon altered the level of p65(Gag) and p180(Gag-Pol) within expressing cells but did not alter the pattern of released VLP, which remained immature. Coexpression of p65(Gag) with a fixed readthrough p180(Gag-Pol) also led to only immature VLP release despite high intracellular protease levels. Our data suggest a mechanism that preferentially selects uncleaved p65(Gag) for the assembly of MLV in this heterologous expression system and implies that, in addition to their relative levels, active sorting of the correct p65(Gag) and p180(Gag-Pol) ratios may occur in producer cells. (C) 2003 Elsevier Inc. All rights reserved.

Influence of particle characteristics and organic matter content on the bioavailability and bioaccumulation of pyrene by clams

Hydrophobic chemicals are known to associate with sediment particles including those from both suspended particulate matter and bottom deposits. The complex and variable composition of natural particles makes it very difficult therefore, to predict the bioavailability of sediment-bound contaminants. To overcome these problems we have previously devised a test system using artificial particles, with or without humic acids, for use as an experimental model of natural sediments. In the present work we have applied this experimental technique to investigate the bioavailability and bioaccumulation of pyrene by the freshwater fingernail clam Sphaerium corneum. The uptake and accumulation of pyrene in clams exposed to the chemical in the presence of a sample of natural sediment was also investigated. According to the results obtained, particle surface properties and organic matter content are the key factors for assessing the bioavailability and bioaccumulation of pyrene by clams. (C) 2002 Elsevier Science Ltd. All rights reserved.

Analysis of ochres from Clearwell Caves: the role of particle size in determining colour

Three ochre samples (A (orange-red in colour), B (red) and C (purple)) from Clearwell Caves, (Gloucestershire, UK) have been examined using an integrated analytical methodology based on the techniques of IR and diffuse reflectance UV-visible-NIR spectroscopy, X-ray diffraction, elemental analysis by ICP-AES and particle size analysis. It is shown that the chromophore in each case is haematite. The differences in colour may be accounted for by (i) different mineralogical and chemical composition in the case of the orange ochre, where hi,,her levels of dolomite and copper are seen and (ii) an unusual particle size distribution in the case of the purple ochre. When the purple ochre was ground to give the same particle size distribution as the red ochre then the colours of the two samples became indistinguishable. An analysis has now been completed of a range of ochre samples with colours from yellow to purple from the important site of Clearwell Caves. (C) 2004 Elsevier B.V. All rights reserved.

Silica-encapsulated nanomagnetic particle as a new recoverable biocatalyst carrier

Enzymes are versatile biocatalysts with major advantages of ultrahigh reaction selectivity and specificity under mild conditions, which currently find increasing applications. However, their applications are often hampered by difficulties in recovery and recycling. As a result, we carried out detailed investigations on the synthesis and characterization of silica-encapsulated iron oxide magnetic nanoparticles of controlled dimension as an enzyme carrier. It is shown that the relatively smaller sized silica-coated magnetic nanoparticle prepared by the microemlusion technique can a carry bulky enzyme, beta-lactamase, via chemical linkages on the silica overlayer without severely blocking the enzymatic active center ( which is commonly encountered in conventional solid supports). An activity study by Michalis-Menten kinetics reflects that this new type of immobilization allows enzyme isolation with accessibility as good as free enzyme. The recovery and reusability of the nanoparticle-supported enzyme upon application of magnetic separation are also demonstrated.

Ultra-thin porous silica coated silver-platinum alloy nano-particle as a new catalyst precursor

Templated sol-gel encapsulation of surfactant-stabilised micelles containing metal precursor(s) with ultra-thin porous silica coating allows solvent extraction of organic based stabiliser from the composites in colloidal state hence a new method of preparing supported alloy catalysts using the inorganic silica-stabilised nano-sized, homogenously mixed, silver - platinum (Ag-Pt) colloidal particles is reported.

Chylomicron particle size and number, factor VII activation and dietary monounsaturated fatty acids

This study evaluated the effects of substituting dietary saturated fatty acids (SFAs) with monounsaturated fatty acids (MUFAs) on postprandial chylomicron (triacylglycerol (TAG), apolipoprotein B-48 (apo B-48) and retinyl ester (RE)), chylomicron particle size and factor VII (FVII) response when subjects were given a standard meal. In a controlled sequential design, 51 healthy young subjects followed an SFA-rich diet (Reference diet) for 8 weeks after which half of the subjects followed a moderate MUFA diet (n = 25) and half followed a high MUFA diet (n = 26) for 16 weeks. Fasting lipoprotein and lipid measurements were evaluated at baseline and at 8-week intervals during the Reference and MUFA diets. In 25 of the subjects (n = 12 moderate MUFA, n = 13 high MUFA), postprandial responses to a standard test meal containing RE and 13 C-tripalmitin were investigated at the end of the Reference and the MUFA diet periods. Although there were no differences in the postprandial lipid markers (TAG, RE, C-13-TAG) on the two diets, the postprandial apo B-48 response (incremental area under the curve (IAUC) was reduced by 21% on the moderate MUFA diet (NS) and by 54% on the high MUFA diet (P < 0.01). The postprandial peak concentrations of apo B-48 were reduced by 33% on the moderate MUFA diet (P < 0.01) and 48% on the high MUFA diet (P < 0.001). Fasting values for factor VII activity (FVIIc), activated factor VII (FVIIa) or factor VII antigen (FVIIag) did not differ significantly when subjects were transferred from Reference to MUFA diets. However, the postprandial increases in coagulation FVII activity (FVIIc) were 18% lower and of activated FVII (FVIIa) were 17% lower on the moderate MUFA diet (NS). Postprandial increases in FVIIc and FVIIa were 50% (P < 0.05) and 29% (P < 0.07) lower on the high MUFA diet and the area under the postprandial FVIIc response curve (AUC) was also lower on the high MUFA diet (P < 0.05). Significantly higher ratios of RE:apo B-48 (P < 0.001) and 13 C-palmitic acid:apo B-48 (P < 0.01) during both MUFA diets suggest that the CMs formed carry larger amounts of dietary lipids per particle, reflecting an adaptation to form larger lipid droplets in the enterocyte when increased amounts of dietary MUFAs are fed. Smaller numbers of larger chylomicrons may explain attenuated activation of factor VII during the postprandial state when the background diet is rich in MUFA. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Saturated fat-induced changes in S-f 60-400 particle composition reduces uptake of LDL by HepG2 cells

The ability of human postprandial triacylglycerol-rich lipoproteins (TRLs), isolated after meals enriched in saturated fatty acids (SFAs), n-6 PUFAs, and MUFAs, to inhibit the uptake of I-125-labeled LDL by the LDL receptor was investigated in HepG2 cells. Addition of TRLs resulted in a dose-dependent inhibition of heparin-releasable binding, cell-associated radioactivity, and degradation products of I-125-labeled LDL (P < 0.001). SFA-rich Svedberg flotation rate (S-f) 60-400 resulted in significantly greater inhibition of cell-associated radioactivity than PUFA-rich particles (P = 0.016) and total uptake of I-125-labeled LDL compared with PUFA- and MUFA-rich particles (P = 0.02). Normalization of the apolipoprotein (apo)E but not apoC-III content of the TRLs removed the effect of meal fatty acid composition, and addition of an anti-apoE antibody reversed the inhibitory effect of TRLs on the total uptake of I-125-labeled LDL. Real time RT-PCR showed that the SFA-rich Sf 60-400 increased the expression of genes involved in hepatic lipid synthesis (P < 0.05) and decreased the expression of the LDL receptor-related protein 1 compared with MUFAs (P = 0.008). In conclusion, these findings suggest an alternative or additional mechanism whereby acute fat ingestion can influence LDL clearance via competitive apoE-dependent effects of TRL on the LDL receptor.-Jackson, K. G., V. Maitin, D. S. Leake, P. Yaqoob, and C. M. Williams. Saturated fat-induced changes in Sf 60 400 particle composition reduces uptake of LDL by HepG2 cells.

Time course of oculomotor inhibition revealed by saccade trajectory modulation

Selecting a stimulus as the target for a goal-directed movement involves inhibiting other competing possible responses. Both target and distractor stimuli activate populations of neurons in topographic oculomotor maps such as the superior colliculus. Local inhibitory interconnections between these populations ensure only one saccade target is selected. Suppressing saccades to distractors may additionally involve inhibiting corresponding map regions to bias the local competition. Behavioral evidence of these inhibitory processes comes from the effects of distractors on oculomotor and manual trajectories. Individual saccades may initially deviate either toward or away from a distractor, but the source of this variability has not been investigated. Here we investigate the relation between distractor-related deviation of trajectory and saccade latency. Targets were presented with, or without, distractors, and the deviation of saccade trajectories arising from the presence of distractors was measured. A fixation gap paradigm was used to manipulate latency independently of the influence of competing distractors. Shorter- latency saccades deviated toward distractors and longer-latency saccades deviated away from distractors. The transition between deviation toward or away from distractors occurred at a saccade latency of around 200 ms. This shows that the time course of the inhibitory process involved in distractor related suppression is relatively slow.

Regression, developmental trajectory and associated problems in disorders in the autism spectrum: the SNAP study

We report rates of regression and associated findings in a population derived group of 255 children aged 9-14 years, participating in a prevalence study of autism spectrum disorders (ASD); 53 with narrowly defined autism, 105 with broader ASD and 97 with non-ASD neurodevelopmental problems, drawn from those with special educational needs within a population of 56,946 children. Language regression was reported in 30% with narrowly defined autism, 8% with broader ASD and less than 3% with developmental problems without ASD. A smaller group of children were identified who underwent a less clear setback. Regression was associated with higher rates of autistic symptoms and a deviation in developmental trajectory. Regression was not associated with epilepsy or gastrointestinal problems.