API-CALC 1.0: a computer program for calculating the average probability of identity allowing for substructure, inbreeding and the presence of close relatives

Individual identification via DNA profiling is important in molecular ecology, particularly in the case of noninvasive sampling. A key quantity in determining the number of loci required is the probability of identity (PIave), the probability of observing two copies of any profile in the population. Previously this has been calculated assuming no inbreeding or population structure. Here we introduce formulae that account for these factors, whilst also accounting for relatedness structure in the population. These formulae are implemented in API-CALC 1.0, which calculates PIave for either a specified value, or a range of values, for F-IS and F-ST.

Enhancing and diminishing gene function in human embryonic stem cells

It is widely recognized that gain- and loss-of-function approaches are essential for understanding the functions of specific genes, and such approaches would be particularly valuable in studies involving human embryonic stem (hES) cells. We describe a simple and efficient approach using lipofection to transfect hES cells, which enabled us to generate hES cell lines expressing naturally fluorescent green or red proteins without affecting cell pluripotency. We used these cell lines to establish a means of diminishing gene function using small interfering (si)RNAs, which were effective at knocking down gene expression in hES cells. We then demonstrated that stable expression of siRNA could knock down the expression of endogenous genes. Application of these gain- and loss-of-function approaches should have widespread use, not only in revealing the developmental roles of specific human genes, but also for their utility in modulating differentiation.

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.

Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.

Wormlike micelle formation and flow alignment of a pluronic block copolymer in aqueous solution

The self-assembly into wormlike micelles of a poly(ethylene oxide)-b-poly(propylene oxide)-b-poly(ethylene oxide) triblock copolymer Pluronic P84 in aqueous salt solution (2 M NaCl) has been studied by rheology, small-angle X-ray and neutron scattering (SAXS/SANS), and light scattering. Measurements of the flow curves by controlled stress rheometry indicated phase separation under flow. SAXS on solutions subjected to capillary flow showed alignment of micelles at intermediate shear rates, although loss of alignment was observed for high shear rates. For dilute solutions, SAXS and static light scattering data on unaligned samples could be superposed over three decades in scattering vector, providing unique information on the wormlike micelle structure over several length scales. SANS data provided information on even shorter length scales, in particular, concerning "blob" scattering from the micelle corona. The data could be modeled based on a system of semiflexible self-avoiding cylinders with a circular cross-section, as described by the wormlike chain model with excluded volume interactions. The micelle structure was compared at two temperatures close to the cloud point (47 degrees C). The micellar radius was found not to vary with temperature in this region, although the contour length increased with increasing temperature, whereas the Kuhn length decreased. These variations result in an increase of the low-concentration radius of gyration with increasing temperature. This was consistent with dynamic light scattering results, and, applying theoretical results from the literature, this is in agreement with an increase in endcap energy due to changes in hydration of the poly(ethylene oxide) blocks as the temperature is increased.

Adapting to the effects of climate change on water supply reliability

This paper presents preliminary results from an assessment of the barriers to adaptation to water supply shortage in a case study catchment in south east England with multiple supply companies. The investigation applies a conceptual framework, which distinguishes between generic barriers affecting the ability of supply companies to make adaptation decisions, and specific barriers to the implementation of each option. The preliminary analysis suggests that whilst there is a widespread awareness of the challenge of climate change, and a conceptual understanding of the need for adaptation, some of the generic barriers that will affect detailed evaluations and actual adaptation decisions have yet to be approached. The analysis also shows that different individual adaptation options are assessed differently by different stakeholders, and that there are differences in the barriers to adoption between supply-side and demand-side measures. First, however, the paper develops the general conceptual framework for the characterisation of the barriers to adaptation used in the study.

Effects of dietary fat modification on oxidative stress and inflammatory markers in the LIPGENE study.

Subjects with the metabolic syndrome (MetS) have enhanced oxidative stress and inflammation. Dietary fat quality has been proposed to be implicated in these conditions. We investigated the impact of four diets distinct in fat quantity and quality on 8-iso-PGF2α (a major F2-isoprostane and oxidative stress indicator), 15-keto-13,14-dihydro-PGF2α (15-keto-dihydro-PGF2α, a major PGF2α metabolite and marker of cyclooxygenase-mediated inflammation) and C-reactive protein (CRP). In a 12-week parallel multicentre dietary intervention study (LIPGENE), 417 volunteers with the MetS were randomly assigned to one of the four diets: two high-fat diets (38 % energy (%E)) rich in SFA or MUFA and two low-fat high-complex carbohydrate diets (28 %E) with (LFHCC n-3) or without (LFHCC) 1·24 g/d of very long chain n-3 fatty acid supplementation. Urinary levels of 8-iso-PGF2α and 15-keto-dihydro-PGF2α were determined by RIA and adjusted for urinary creatinine levels. Serum concentration of CRP was measured by ELISA. Neither concentrations of 8-iso-PGF2α and 15-keto-dihydro-PGF2α nor those of CRP differed between diet groups at baseline (P>0·07) or at the end of the study (P>0·44). Also, no differences in changes of the markers were observed between the diet groups (8-iso-PGF2α, P = 0·83; 15-keto-dihydro-PGF2α, P = 0·45; and CRP, P = 0·97). In conclusion, a 12-week dietary fat modification did not affect the investigated markers of oxidative stress and inflammation among subjects with the MetS in the LIPGENE study.

Effects of dietary fat modification on insulin sensitivity and on other risk factors of the metabolic syndrome—LIPGENE: a European randomized dietary intervention study

Background:Excessive energy intake and obesity lead to the metabolic syndrome (MetS). Dietary saturated fatty acids (SFAs) may be particularly detrimental on insulin sensitivity (SI) and on other components of the MetS. Objective:This study determined the relative efficacy of reducing dietary SFA, by isoenergetic alteration of the quality and quantity of dietary fat, on risk factors associated with MetS. Design:A free-living, single-blinded dietary intervention study. Subjects and Methods:MetS subjects (n=417) from eight European countries completed the randomized dietary intervention study with four isoenergetic diets distinct in fat quantity and quality: high-SFA; high-monounsaturated fatty acids and two low-fat, high-complex carbohydrate (LFHCC) diets, supplemented with long chain n-3 polyunsaturated fatty acids (LC n-3 PUFAs) (1.2 g per day) or placebo for 12 weeks. SI estimated from an intravenous glucose tolerance test (IVGTT) was the primary outcome measure. Lipid and inflammatory markers associated with MetS were also determined. Results:In weight-stable subjects, reducing dietary SFA intake had no effect on SI, total and low-density lipoprotein cholesterol concentration, inflammation or blood pressure in the entire cohort. The LFHCC n-3 PUFA diet reduced plasma triacylglycerol (TAG) and non-esterified fatty acid concentrations (P<0.01), particularly in men. Conclusion:There was no effect of reducing SFA on SI in weight-stable obese MetS subjects. LC n-3 PUFA supplementation, in association with a low-fat diet, improved TAG-related MetS risk profiles.

Self-assembly of a modified amyloid peptide fragment: pH-responsiveness and nematic phase formation

The self-assembly of peptide YYKLVFFC based on a fragment of the amyloid beta (A) peptide, A beta 16-20, KLVFF has been studied in aqueous solution. The peptide is designed with multiple functional residues to examine the interplay between aromatic interactions and charge on the self-assembly, as well as specific transformations such as the pH-induced phenol-phenolate transition of the tyrosine residue. Circular dichroism (CD) and Fourier-transform infrared (FTIR) spectroscopies are used to investigate the conditions for beta-sheet self-assembly and the role of aromatic interactions in the CD spectrum as a function of pH and concentration. The formation of well-defined fibrils at pH 4.7 is confirmed by cryo-TEM (transmission electron microscope) and negative stain TEM. The morphology changes at higher pH, and aggregates of short twisted fibrils are observed at pH 11. Polarized optical microscopy shows birefringence at a low concentration (1 wt.-%) of YYKLVFFC in aqueous solution, and small-angle X-ray scattering was used to probe nematic phase formation in more detail. A pH-induced transition from nematic to isotropic phases is observed on increasing pH that appears to be correlated to a reduction in aggregate anisotropy upon increasing pH.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting thatmost effectors represent species-specific adaptations.

Impact of stratospheric variability on tropospheric climate change

An improved stratospheric representation has been included in simulations with the Hadley Centre HadGEM1 coupled ocean atmosphere model with natural and anthropogenic forcings for the period 1979–2003. An improved stratospheric ozone dataset is employed that includes natural variations in ozone as well as the usual anthropogenic trends. In addition, in a second set of simulations the quasi biennial oscillation (QBO) of stratospheric equatorial zonal wind is also imposed using a relaxation towards ERA-40 zonal wind values. The resulting impact on tropospheric variability and trends is described. We show that the modelled cooling rate at the tropopause is enhanced by the improved ozone dataset and this improvement is even more marked when the QBO is also included. The same applies to warming trends in the upper tropical troposphere which are slightly reduced. Our stratospheric improvements produce a significant increase of internal variability but no change in the positive trend of annual mean global mean near-surface temperature. Warming rates are increased significantly over a large portion of the Arctic Ocean. The improved stratospheric representation, especially the QBO relaxation, causes a substantial reduction in near-surface temperature and precipitation response to the El Chichón eruption, especially in the tropical region. The winter increase in the phase of the northern annular mode observed in the aftermath of the two major recent volcanic eruptions is partly captured, especially after the El Chichón eruption. The positive trend in the southern annular mode (SAM) is increased and becomes statistically significant which demonstrates that the observed increase in the SAM is largely subject to internal variability in the stratosphere. The possible inclusion in simulations for future assessments of full ozone chemistry and a gravity wave scheme to internally generate a QBO is discussed.

Dietary prebiotics: current status and new definition

In November 2008, a group of scientists met at the 6th Meeting of the International Scientific Association of Probiotics and Prebiotics (ISAPP) in London, Ontario, Canada, to discuss the functionality of prebiotics. As a result of this, it was concluded that the prebiotic field is currently dominated by gastrointestinal events. However, in the future, it may be the case that other mixed microbial ecosystems may be modulated by a prebiotic approach, such as the oral cavity, skin and the urogenital tract. Therefore, a decision was taken to build upon the current prebiotic status and define a niche for ‘dietary prebiotics’. This review is co-authored by the working group of ISAPP scientists and sets the background for defining a dietary prebiotic as ‘‘a selectively fermented ingredient that results in specific changes in the composition and/or activity of the gastrointestinal microbiota, thus conferring benefit(s) upon host health’’.

Humoral immunity following chickenpox is influenced by geography and ethnicity

Objectives To investigate the contribution of ethnicity and geographical location to varicella-zoster virus (VZV) serostatus and antibody concentrations. Methods The presence and concentrations of antibodies to VZV were measured in 639 Bangladeshi women born in Bangladesh (BBB), 94 Bangladeshi women born in the UK (BUK) and 262 White women born in the UK (WUK). The results were anaylsed in relation to demographic and social data. Results BBB women were significantly less likely to be VZV seropositive at all ages than both BUK and WUK women. However, the odds of a Bangladeshi-born woman being seropositive increased by 1.04 for each year under the age of 15 spent in the UK. In contrast, antibody concentrations were significantly lower in ethnic Bangladeshi women, irrespective of country of birth. White, but not Bangladeshi women, showed evidence of antibody boosting over time despite the latter having more exposure to children. Conclusion Geographical location during childhood is the major influence on age of primary infection with VZV while the level of antibody is related to ethnicity. Since the risk of re-infection with VZV following both natural infection and vaccination is increased as antibody concentrations fall, these results have implications for VZV vaccination programmes particularly in non-White populations.

Measuring departures from Hardy-Weinberg: a Markov chain Monte Carlo method for estimating the inbreeding coefficient

Many well-established statistical methods in genetics were developed in a climate of severe constraints on computational power. Recent advances in simulation methodology now bring modern, flexible statistical methods within the reach of scientists having access to a desktop workstation. We illustrate the potential advantages now available by considering the problem of assessing departures from Hardy-Weinberg (HW) equilibrium. Several hypothesis tests of HW have been established, as well as a variety of point estimation methods for the parameter which measures departures from HW under the inbreeding model. We propose a computational, Bayesian method for assessing departures from HW, which has a number of important advantages over existing approaches. The method incorporates the effects-of uncertainty about the nuisance parameters--the allele frequencies--as well as the boundary constraints on f (which are functions of the nuisance parameters). Results are naturally presented visually, exploiting the graphics capabilities of modern computer environments to allow straightforward interpretation. Perhaps most importantly, the method is founded on a flexible, likelihood-based modelling framework, which can incorporate the inbreeding model if appropriate, but also allows the assumptions of the model to he investigated and, if necessary, relaxed. Under appropriate conditions, information can be shared across loci and, possibly, across populations, leading to more precise estimation. The advantages of the method are illustrated by application both to simulated data and to data analysed by alternative methods in the recent literature.