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).     

Multi-Sample Pooling and Illumina Genome Analyzer Sequencing Methods to Determine Gene Sequence Variation for Database Development

Determination of sequence variation within a genetic locus to develop clinically relevant databases is critical for molecular assay design and clinical test interpretation, so multisample pooling for Illumina genome analyzer (GA) sequencing was investigated using the RET proto-oncogene as a model. Samples were Sanger-sequenced for RET exons 10, 11, and 13–16. Ten samples with 13 known unique variants (“singleton variants” within the pool) and seven common changes were amplified and then equimolar-pooled before sequencing on a single flow cell lane, generating 36 base reads. For comparison, a single “control” sample was run in a different lane. After alignment, a 24-base quality score-screening threshold and 3` read end trimming of three bases yielded low background error rates with a 27% decrease in aligned read coverage. Sequencing data were evaluated using an established variant detection method (percent variant reads), by the presented subtractive correction method, and with SNPSeeker software. In total, 41 variants (of which 23 were singleton variants) were detected in the 10 pool data, which included all Sanger-identified variants. The 23 singleton variants were detected near the expected 5% allele frequency (average 5.17%±0.90% variant reads), well above the highest background error (1.25%). Based on background error rates, read coverage, simulated 30, 40, and 50 sample pool data, expected singleton allele frequencies within pools, and variant detection methods; ≥30 samples (which demonstrated a minimum 1% variant reads for singletons) could be pooled to reliably detect singleton variants by GA sequencing.     

2D-DiGE analysis of the human endometrial secretome reveals differences between receptive and nonreceptive states in fertile and infertile women

Endometrial secretions in the uterine cavity contain mediators important for endometrial receptivity and embryo implantation. Unbiased analysis of uterine fluid from a receptive versus nonreceptive time of the menstrual cycle and in fertile and infertile women will provide new insights into uterine receptivity. We hypothesized that proteomic analysis of human uterine lavages would identify proteins important for the establishment of pregnancy in humans. Lavages collected from fertile (n = 7) and infertile (n = 8) women during the midsecretory (MS) phase, and from fertile women during the midproliferative (MP) (n = 7) phase, were assessed using 2D-differential in gel electrophoresis (2D-DiGE) over a pI 4-7 range. Statistical analysis revealed 7 spots that were significantly decreased in the MP compared to the MS phase, while 18 spots showed differential expression between fertile and infertile women. A number of proteins were identified by mass spectrometry, including antithrombin III and alpha-2-macroglobulin, whose production was confirmed in endometrial epithelium. Their staining pattern suggests roles during embryo implantation. Assessment of the human endometrial secretome has identified differences in the protein content of uterine fluid with respect to receptivity and fertility.     

Glioblastoma multiforme arising in the irradiated spinal cord of a rhesus monkey (Macaca mulatta)

Abstract: An adult female rhesus monkey that had received 44.0 Gy of cobalt 60 radiation to 8 cm of the cervical and upper thoracic spinal cord approximately 2.8 years postirradiation developed a sudden onset of self-mutilation and loss of function of the right arm followed progressively by loss of function of the left arm and terminally bilateral paresis of the legs. Histopathologic examination of the cervical spinal cord revealed a glioblastoma multiforme that extended from the cervical medullary junction to the sixth cervical vertebrae. Because of the infrequent occurrence of spontaneous neoplasia in rhesus monkeys and the location in the radiation field, the glioblastoma is believed to be radiation induced.     

Directed evolution of the thermostable xylanase from Thermomyces lanuginosus

The thermostability of the endo-beta-1,4-xylanase from Thermomyces lanuginosus (xynA) was improved by directed evolution using error-prone PCR. Transformants expressing the variant xylanases were first selected on 0.4% Remazol Brilliant Blue-xylan and then exposed to 80 degrees C. Whereas the wild type XynA lost 90% activity after 10 min at 80 degrees C, five mutants displayed both higher stabilities and activities than XynA. Four mutants were subjected to further mutagenesis to improve the stability and activity of the xylanase. Subsequent screening revealed three mutants with enhanced thermostability. Mutant 2B7-10 retained 71% of its activity after treatment at 80 degrees C for 60 min and had a half-life of 215 min at 70 degrees C, which is higher than that attained by XynA. Sequence analysis of second generation mutants revealed that mutations were not concentrated in any particular region of the protein and exhibited much variation. The best mutant obtained from this study was variant 2B7-10, which had a single substitution (Y58F) in beta-sheet A of the protein, which is the hydrophilic, solvent-accessible outer surface of the enzyme. Most of the mutants obtained in this study displayed a compromise between stability and activity, the only exception being mutant 2B7-10. This variant showed increased activity and thermostability.