fastSTRUCTURE: Variational Inference of Population Structure in Large SNP Data Sets

Tools for estimating population structure from genetic data are now used in a wide variety of applications in population genetics. However, inferring population structure in large modern data sets imposes severe computational challenges. Here, we develop efficient algorithms for approximate inference of the model underlying the STRUCTURE program using a variational Bayesian framework. Variational methods pose the problem of computing relevant posterior distributions as an optimization problem, allowing us to build on recent advances in optimization theory to develop fast inference tools. In addition, we propose useful heuristic scores to identify the number of populations represented in a data set and a new hierarchical prior to detect weak population structure in the data. We test the variational algorithms on simulated data and illustrate using genotype data from the CEPH–Human Genome Diversity Panel. The variational algorithms are almost two orders of magnitude faster than STRUCTURE and achieve accuracies comparable to those of ADMIXTURE. Furthermore, our results show that the heuristic scores for choosing model complexity provide a reasonable range of values for the number of populations represented in the data, with minimal bias toward detecting structure when it is very weak. Our algorithm, fastSTRUCTURE, is freely available online at http://pritchardlab.stanford.edu/structure.html.     

A social science view on the FRAME symposium: identities and networks

This paper, written by two social scientists, presents a social science perspective on the issues raised at the FRAME symposium on Human Alternatives to Animal Studies. Drawing upon the authors' experience of conducting research with stem cell scientists, issues around access to human tissue for in vitro uses are considered. The paper concludes by raising questions pertinent to both interested social scientists and the Three Rs agenda.     

Evaluation of extraction processes for intracellular metabolite profiling of mammalian cells: matching extraction approaches to cell type and metabolite targets

In this study we report on the optimisation of the technologies for generation of a global metabolomics profile for intracellular metabolites in Chinese hamster ovary (CHO) cells. We evaluated the effectiveness of a range of different extraction methods applied to CHO cells which had been quenched using a previously optimised approach. The extraction methods tested included cold methanol, hot ethanol, acid, alkali and methanol/chloroform plus combinations of these. The extraction of metabolites using two 100% methanol extractions followed by a final water extraction recovered the largest range of metabolites. For the majority of metabolites, extracts generated in this manner exhibited the greatest recovery with high reproducibility. Therefore, this was the best extraction method for attaining a global metabolic profile from a single sample. However, another parallel extraction method (e.g. alkali) may also be required to maximise the range of metabolites recovered (e.g. non-polar metabolites).     

Role of CheB and CheR in the complex chemotactic and aerotactic pathway of Azospirillum brasilense

It has previously been reported that the alpha-proteobacterium Azospirillum brasilense undergoes methylation-independent chemotaxis; however, a recent study revealed cheB and cheR genes in this organism. We have constructed cheB, cheR, and cheBR mutants of A. brasilense and determined that the CheB and CheR proteins under study significantly influence chemotaxis and aerotaxis but are not essential for these behaviors to occur. First, we found that although cells lacking CheB, CheR, or both were no longer capable of responding to the addition of most chemoattractants in a temporal gradient assay, they did show a chemotactic response (albeit reduced) in a spatial gradient assay. Second, in comparison to the wild type, cheB and cheR mutants under steady-state conditions exhibited an altered swimming bias, whereas the cheBR mutant and the che operon mutant did not. Third, cheB and cheR mutants were null for aerotaxis, whereas the cheBR mutant showed reduced aerotaxis. In contrast to the swimming bias for the model organism Escherichia coli, the swimming bias in A. brasilense cells was dependent on the carbon source present and cells released methanol upon addition of some attractants and upon removal of other attractants. In comparison to the wild type, the cheB, cheR, and cheBR mutants showed various altered patterns of methanol release upon exposure to attractants. This study reveals a significant difference between the chemotaxis adaptation system of A. brasilense and that of the model organism E. coli and suggests that multiple chemotaxis systems are present and contribute to chemotaxis and aerotaxis in A. brasilense.     

Turkey erythrocyte membranes as a model for regulation of phospholipase C by guanine nucleotides

Phosphoinositides of human, rabbit, rat, and turkey erythrocytes were radiolabeled by incubation of intact cells with [32P]Pi. Guanosine 5'-O-(thiotriphosphate) (GTP gamma S) and NaF, which are known activators of guanine nucleotide regulatory proteins, caused a large increase in [32P]inositol phosphate release from plasma membranes derived from turkey erythrocytes, but had no effect on inositol phosphate formation by plasma membranes prepared from the mammalian erythrocytes. High performance liquid chromatography analysis indicated that inositol bisphosphate, inositol 1,3,4-trisphosphate, inositol 1,4,5-trisphosphate, and inositol 1,3,4,5-tetrakisphosphate all increased by 20-30-fold during a 10-min incubation of turkey erythrocyte membranes with GTP gamma S. The increase in inositol phosphate formation was accompanied by a similar decrease in radioactivity in phosphatidylinositol 4-phosphate (PIP) and phosphatidylinositol 4,5-bisphosphate (PIP2). GTP gamma S increased inositol phosphate formation with a K0.5 of 600 nM; guanosine 5'-(beta, gamma-imido)trisphosphate was 50-75% as efficacious as GTP gamma S and expressed a K0.5 of 36 microM. Although GTP alone had little effect on inositol phosphate formation, it blocked GTP gamma S-stimulated inositol phosphate formation, as did guanosine 5'-O-(2-thiodiphosphate). Turkey erythrocytes were also shown to express phosphatidylinositol synthetase activity in that incubation of cells with [3H] inositol resulted in incorporation of radiolabel into phosphatidylinositol, PIP, and PIP2. Incubation of membranes derived from [3H]inositol-labeled erythrocytes with GTP gamma S resulted in large increases in [3H] inositol phosphate formation and corresponding decreases in radiolabel in PIP and PIP2. The data suggest that, in contrast to mammalian erythrocytes, the turkey erythrocyte expresses a guanine nucleotide-binding protein that regulates phospholipase C, and as such, should provide a useful model system for furthering our understanding of hormonal regulation of this enzyme.     

Accounting for Carbon Flux to Mycorrhizal Fungi May Resolve Discrepancies in Forest Carbon Budgets

Ecosystems - Carbon (C) fluxes among different components of plant growth are important to forest ecosystem C cycling and are strongly influenced by species composition and resource availability....     

Multiple metabolic pools of phosphoinositides and phosphatidate in human erythrocytes incubated in a medium that permits rapid transmembrane exchange of phosphate.

1. A Hepes-based medium has been devised which allows rapid Pi exchange across the plasma membrane of the human erythrocyte. This allows the metabolically labile phosphate pools of human erythrocytes to come to equilibrium with [32P]Pi in the medium after only 5 h in vitro. 2. After 5-7 h incubation with [32P]Pi in this medium, only three phospholipids, phosphatidic acid (PtdOH), phosphatidylinositol 4-phosphate (PtdIns4P) and phosphatidylinositol 4,5-bisphosphate (PtdIns4,5P2) are radioactively labelled. The concentrations of PtdIns4P and PtdIns4,5P2 remain constant throughout the incubation, so this labelling process is a reflection of the steady-state turnover of their monoester phosphate groups. 3. During such incubations, the specific radioactivities of the monoesterified phosphates of PtdIns4, PtdIns4,5P2 and PtdOH come to a steady value after 5 h that is only 25-30% of the specific radioactivity of the gamma-phosphate of ATP at that time. We suggest that this is a consequence of metabolic heterogeneity. This heterogeneity is not a result of the heterogeneous age distribution of the erythrocytes in human blood. Thus it appears that there is metabolic compartmentation of these lipids within cells, such that within a time-scale of a few hours only 25-30% of these three lipids are actively metabolized. 4. The phosphoinositidase C of intact human erythrocytes, when activated by Ca2+-ionophore treatment, only hydrolyses 50% of the total PtdIns4,5P2 and 50% of 32P-labelled PtdIns4,5P2 present in the cells: this enzyme does not discriminate between the metabolically active and inactive compartments of lipids in the erythrocyte membrane. Hence at least four metabolic pools of PtdIns4P and PtdIns4,5P2 are distinguishable in the human erythrocyte plasma membrane. 5. The mechanisms by which multiple non-mixing metabolic pools of PtdOH, PtdIns4P and PtdIns4,5P2 are sustained over many hours in the plasma membranes of intact erythrocytes are unknown, although some possible explanations are considered.     

Gene organisation, sequence variation and isochore structure at the centromeric boundary of the human MHC.

We have mapped and sequenced the region immediately centromeric of the human major histocompatibility complex (MHC). A cluster of 13 genes/pseudogenes was identified in a 175 kb PAC linking the TAPASIN locus with the class II region. It includes two novel human genes (BING4 and SACM2L) and a thus far unnoticed human leucocyte antigen (HLA) class II pseudogene, termed HLA-DPA3. Analysis of the G+C content revealed an isochore boundary which, together with the previously reported telomeric boundary, defines the MHC class II region as one of the first completely sequenced isochores in the human genome. Comparison of the sequence with limited sequence from other cell lines shows that the high sequence variation found within the classical class II region extends beyond the identified isochore boundary leading us to propose the concept of an "extended MHC". By comparative analysis, we have precisely identified the mouse/human synteny breakpoint at the centromeric end of the extended MHC class II region between the genes HSET and PHF1.