Biopython: freely available Python tools for computational molecular biology and bioinformatics

Summary: The Biopython project is a mature open source international collaboration of volunteer developers, providing Python libraries for a wide range of bioinformatics problems. Biopython includes modules for reading and writing different sequence file formats and multiple sequence alignments, dealing with 3D macro molecular structures, interacting with common tools such as BLAST, ClustalW and EMBOSS, accessing key online databases, as well as providing numerical methods for statistical learning. Availability: Biopython is freely available, with documentation and source code at www.biopython.org under the Biopython license. Contact: All queries should be directed to the Biopython mailing lists, see www.biopython.org/wiki/[email protected]

Buccal swab as a reliable predictor for X inactivation ratio in inaccessible tissues

Background As a result of the epigenetic phenomenon of X chromosome inactivation (XCI) every woman is a mosaic of cells with either an inactive paternal X chromosome or an inactive maternal X chromosome. The ratio between inactive paternal and maternal X chromosomes is different for every female individual, and can influence an X-encoded trait or disease. A multitude of X linked conditions is known, and for many of them it is recognised that the phenotype in affected female carriers of the causative mutation is modulated by the XCI ratio. To predict disease severity an XCI ratio is usually determined in peripheral blood samples. However, the correlation between XCI ratios in peripheral blood and disease affected tissues, that are often inaccessible, is poorly understood. Here, we tested several tissues obtained from autopsies of 12 female individuals for patch size and XCI ratio. Methods XCI ratios were analysed using methylsensitive PCR-based assays for the AR, PCSK1N and SLITRK4 loci. XCI patch size was analysed by testing the XCI ratio of tissue samples with decreasing size. Results XCI patch size was analysed for liver, muscle, ovary and brain samples and was found too small to confound testing for XCI ratio in these tissues. XCI ratios were determined in the easily accessible tissues, blood, buccal epithelium and hair follicle, and compared with ratios in several inaccessible tissues. Conclusions Buccal epithelium is preferable over peripheral blood for predicting XCI ratios of inaccessible tissues. Ovary is the only inaccessible tissue showing a poor correlation to blood and buccal epithelium, but has a good correlation to hair follicle instead

Pointlike structure for super p-branes

We present an efficient method to understand the p-brane dynamics in a unified framework. For this purpose, we reformulate the action for super p-branes in the form appropriate to incorporate the pointlike (parton) structure of higher dimensional p-branes and intend to interpret the p-brane dynamics as the collective dynamics of superparticles. In order to examine such a parton picture of super p-branes, we consider various superparticle configurations that can be reduced from super p-branes, especially, a supermembrane, and study the partonic structure of classical p-brane solutions.Comment: 22 pages, corrected typos, to appear in Phys. Rev. D58, 085018 (1998

The probability to initiate X chromosome inactivation is determined by the X to autosomal ratio and x chromosome specific allelic properties.

Background: In female mammalian cells, random X chromosome inactivation (XCI) equalizes the dosage of X-encoded gene products to that in male cells. XCI is a stochastic process, in which each X chromosome has a probability to be inactivated. To obtain more insight in the factors setting up this probability, we studied the role of the X to autosome (X:A) ratio in initiation of XCI, and have used the experimental data in a computer simulation model to study the cellular population dynamics of XCI. Methodology/Principal Findings: To obtain more insight in the role of the X:A ratio in initiation of XCI, we generated triploid mouse ES cells by fusion of haploid round spermatids with diploid female and male ES cells. These fusion experiments resulted in only XXY triploid ES cells. XYY and XXX ES lines were absent, suggesting cell death related either to insufficient X-chromosomal gene dosage (XYY) or to inheritance of an epigenetically modified X chromosome (XXX). Analysis of active (Xa) and inactive (Xi) X chromosomes in the obtained triploid XXY lines indicated that the initiation frequency of XCI is low, resulting in a mixed population of XaXiY and XaXaY cells, in which the XaXiY cells have a small proliferative advantage. This result, and findings on XCI in diploid and tetraploid ES cell lines with different X:A ratios, provides evidence that the X:A ratio determines the probability for a given X chromosome to be inactivated. Furthermore, we found that the kinetics of the XCI process can be simulated using a probability for an X chromosome to be inactivated that is proportional to the X:A ratio. These simulation studies re-emphasize our hypothesis that the probability is a function of the concentration of an X-encoded activator of XCI, and of X chromosome specific allelic properties determining the threshold for this activator. Conclusions: The present findings reveal that the probability for an X chromosome to be inactivated is proportional to the X:A ratio. This finding supports the presence of an X-encoded activator of the XCI process. © 2009 Monkhorst et al

De-Trending Time Series for Astronomical Variability Surveys

We present a de-trending algorithm for the removal of trends in time series. Trends in time series could be caused by various systematic and random noise sources such as cloud passages, changes of airmass, telescope vibration or CCD noise. Those trends undermine the intrinsic signals of stars and should be removed. We determine the trends from subsets of stars that are highly correlated among themselves. These subsets are selected based on a hierarchical tree clustering algorithm. A bottom-up merging algorithm based on the departure from normal distribution in the correlation is developed to identify subsets, which we call clusters. After identification of clusters, we determine a trend per cluster by weighted sum of normalized light-curves. We then use quadratic programming to de-trend all individual light-curves based on these determined trends. Experimental results with synthetic light-curves containing artificial trends and events are presented. Results from other de-trending methods are also compared. The developed algorithm can be applied to time series for trend removal in both narrow and wide field astronomy.Comment: Revised version according to the referee's second revie

One-way multigrid method in electronic structure calculations

We propose a simple and efficient one-way multigrid method for self-consistent electronic structure calculations based on iterative diagonalization. Total energy calculations are performed on several different levels of grids starting from the coarsest grid, with wave functions transferred to each finer level. The only changes compared to a single grid calculation are interpolation and orthonormalization steps outside the original total energy calculation and required only for transferring between grids. This feature results in a minimal amount of code change, and enables us to employ a sophisticated interpolation method and noninteger ratio of grid spacings. Calculations employing a preconditioned conjugate gradient method are presented for two examples, a quantum dot and a charged molecular system. Use of three grid levels with grid spacings 2h, 1.5h, and h decreases the computer time by about a factor of 5 compared to single level calculations.Comment: 10 pages, 2 figures, to appear in Phys. Rev. B, Rapid Communication

Simcluster: clustering enumeration gene expression data on the simplex space

Transcript enumeration methods such as SAGE, MPSS, and sequencing-by-synthesis EST "digital northern", are important high-throughput techniques for digital gene expression measurement. As other counting or voting processes, these measurements constitute compositional data exhibiting properties particular to the simplex space where the summation of the components is constrained. These properties are not present on regular Euclidean spaces, on which hybridization-based microarray data is often modeled. Therefore, pattern recognition methods commonly used for microarray data analysis may be non-informative for the data generated by transcript enumeration techniques since they ignore certain fundamental properties of this space.

Here we present a software tool, Simcluster, designed to perform clustering analysis for data on the simplex space. We present Simcluster as a stand-alone command-line C package and as a user-friendly on-line tool. Both versions are available at: http://xerad.systemsbiology.net/simcluster.

Simcluster is designed in accordance with a well-established mathematical framework for compositional data analysis, which provides principled procedures for dealing with the simplex space, and is thus applicable in a number of contexts, including enumeration-based gene expression data

Maxwell Chern-Simons Solitons from Type IIB String Theory

We study various three-dimensional supersymmetric Maxwell Chern-Simons solitons by using type IIB brane configurations. We give a systematic classification of soliton spectra such as topological BPS vortices and nontopological vortices in $\mathcal{N}=2,3$ supersymmetric Maxwell Chern-Simons system via the branes of type IIB string theory. We identify the brane configurations with the soliton spectra of the field theory and obtain a nice agreement with field theory aspects. We also discuss possible brane constructions for BPS domain wall solutions.Comment: 23 pages, Latex, 4 figures; (q_1,q_2)-string convention changed, minor correction

SO/Sp Monopoles and Branes with Orientifold 3 Plane

We study BPS monopoles in 4 dimensional N=4 SO(N) and $Sp(N)$ super Yang-Mills theories realized as the low energy effective theory of $N$ (physical and its mirror) parallel D3 branes and an {\it Orientifold 3 plane} with D1 branes stretched between them in type IIB string theory. Monopoles on D3 branes give the natural understanding by embedding in SU(N) through the constraints on both the asymptotic Higgs field (corresponding to the horizontal positions of D3 branes) and the magnetic charges (corresponding to the number of D1 branes) imposed by the O3 plane. The compatibility conditions of Nahm data for monopoles for these groups can be interpreted very naturally through the D1 branes in the presence of O3 plane.Comment: 18 pages, Latex with RevTex, 1 table, 4 figures, v2: Clarified the introduction and improved on the supersymmetric theory on D1 branes in page 7 and 8 and this final version to appear in Phys.Rev.