Motivation: The enormous amount of short reads generated by the new DNA sequencing technologies call for the development of fast and accurate read alignment programs. A first generation of hash table-based methods has been developed, including MAQ, which is accurate, feature rich and fast enough to align short reads from a single individual. However, MAQ does not support gapped alignment for single-end reads, which makes it unsuitable for alignment of longer reads where indels may occur frequently. The speed of MAQ is also a concern when the alignment is scaled up to the resequencing of hundreds of individuals.

Results: We implemented Burrows-Wheeler Alignment tool (BWA), a new read alignment package that is based on backward search with Burrows–Wheeler Transform (BWT), to efficiently align short sequencing reads against a large reference sequence such as the human genome, allowing mismatches and gaps. BWA supports both base space reads, e.g. from Illumina sequencing machines, and color space reads from AB SOLiD machines. Evaluations on both simulated and real data suggest that BWA is ~10–20× faster than MAQ, while achieving similar accuracy. In addition, BWA outputs alignment in the new standard SAM (Sequence Alignment/Map) format. Variant calling and other downstream analyses after the alignment can be achieved with the open source SAMtools software package.

Availability: http://maq.sourceforge.net 

Contact: [email protected]

/pdf/fast-and-accurate-short-read-alignment-with-burrows-wheeler-43uqhgr33c.pdf

Fast and accurate short read alignment with Burrows–Wheeler transform

Genomic sequencing has made it clear that a large fraction of the genes specifying the core biological functions are shared by all eukaryotes. Knowledge of the biological role of such shared proteins in one organism can often be transferred to other organisms. The goal of the Gene Ontology Consortium is to produce a dynamic, controlled vocabulary that can be applied to all eukaryotes even as knowledge of gene and protein roles in cells is accumulating and changing. To this end, three independent ontologies accessible on the World-Wide Web (http://www.geneontology.org) are being constructed: biological process, molecular function and cellular component.

/pdf/gene-ontology-tool-for-the-unification-of-biology-ppsktrcye3.pdf

Gene Ontology: tool for the unification of biology

There is, I think, something ethereal about i —the square root of minus one. I remember first hearing about it at school. It seemed an odd beast at that time—an intruder hovering on the edge of reality.

Usually familiarity dulls this sense of the bizarre, but in the case of i it was the reverse: over the years the sense of its surreal nature intensified. It seemed that it was impossible to write mathematics that described the real world in …

I and i

Bowtie is an ultrafast, memory-efficient alignment program for aligning short DNA sequence reads to large genomes. For the human genome, Burrows-Wheeler indexing allows Bowtie to align more than 25 million reads per CPU hour with a memory footprint of approximately 1.3 gigabytes. Bowtie extends previous Burrows-Wheeler techniques with a novel quality-aware backtracking algorithm that permits mismatches. Multiple processor cores can be used simultaneously to achieve even greater alignment speeds. Bowtie is open source http://bowtie.cbcb.umd.edu.

/pdf/ultrafast-and-memory-efficient-alignment-of-short-dna-3mrm42jmnc.pdf

Ultrafast and memory-efficient alignment of short DNA sequences to the human genome

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

We study the stationary nonequilibrium states of the van Beijeren/Schulman model of a driven lattice gas in two dimensions. In this model, jumps are much faster in the direction of the driving force than orthogonal to it. Van Kampen's Ω-expansion provides a suitable description of the model in the high-temperature region and specifies the critical temperature and the spinodal curve. We find the rate dependence ofT
c and show that independently of the jump rates the critical exponents of the transition are classical, except for anomalous energy fluctuations. We then study the stationary solution of the deterministic equations (zeroth-orderΩ-expansion). They can be obtained as trajectories of a dissipative dynamical system with a three-dimensional phase space. Within a certain temperature range belowT
c, these equations have a kink solution whose asymptotic densities we identify with those of phase coexistence. They appear to coincide with the results of the “Maxwell construction.” This provides a dynamical justification for the use of this construction in this nonequilibrium model. The relation of the Freidlin-Wentzell theory of small random perturbations of dynamical systems to the steady-state distribution belowT
c is discussed.

The fast rate limit of driven diffusive systems

Motivation:The interaction between transcription factor (TF) and transcription factor binding site (TFBS) is essential for gene regulation. Mutation in either the TF or the TFBS may weaken their interaction and thus result in abnormalities. To maintain such vital interaction, a mutation in one of the interacting partners might be compensated by a corresponding mutation in its binding partner during the course of evolution. Confirming this co-evolutionary relationship will guide us in designing protein sequences to target a specific DNA sequence or in predicting TFBS for poorly studied proteins, or even correcting and rescuing disease mutations in clinical applications.

Results: Based on six, publicly available, experimentally validated TF–TFBS binding datasets for the basic Helix–Loop–Helix (bHLH) family, Homeo family, High-Mobility Group (HMG) family and Transient Receptor Potential channels (TRP) family, we showed that the evolutions of the TFs and their TFBSs are significantly correlated across eukaryotes. We further developed a mutual information-based method to identify co-evolved protein residues and DNA bases. This research sheds light on the dynamic relationship between TF and TFBS during their evolution. The same principle and strategy can be applied to co-evolutionary studies on protein–DNA interactions in other protein families.

Availability: All the datasets, scripts and other related files have been made freely available at: http://jjwanglab.org/co-evo.

Contact: junwen@uw.edu

Supplementary Information:Supplementary data are available at Bioinformatics online.

/pdf/correlated-evolution-of-transcription-factors-and-their-5dv9en4k80.pdf

Correlated evolution of transcription factors and their binding sites

Motivation: Traditional methods to identify potential binding sites of known transcription factors still suffer from large number of false predictions. They mostly use sequence information in a position-specific manner and neglect other types of information hidden in the proximal promoter regions. Recent biological and computational researches, however, suggest that there exist not only locational preferences of binding, but also correlations between transcription factors. Results: In this paper, we propose a novel approach, OSCAR, which utilizes one-class SVM algorithms, and incorporates multiple factors to aid the recognition of transcription factor binding sites. Using both synthetic and real data, we find that our method outperforms existing algorithms, especially in the high sensitivity region. The performance of our method can be further improved by taking into account locational preference of binding events. By testing on experimentally-verified binding sites of GATA and HNF transcription factor families, we show that our algorithm can infer the true co-occurring motif pairs accurately, and by considering the co-occurrences of correlated motifs, we not only filter out false predictions, but also increase the sensitivity. Availability: An online server based on OSCAR is available at http://bioinfo.au.tsinghua.edu.cn/oscar.

/pdf/oscar-one-class-svm-for-accurate-recognition-of-cis-elements-niwpe034wx.pdf

OSCAR: one-class SVM for accurate recognition of cis-elements

Heterochromatin remodeling is critical for various cell processes. In particular, the "loss of heterochromatin" phenotype in cellular senescence is associated with the process of aging and age-related disorders. Although biological processes of senescent cells, including senescence-associated heterochromatin foci (SAHF) formation, chromosome compaction, and redistribution of key proteins, have been closely associated with high-order chromatin structure, the relationship between the high-order chromatin reorganization and the loss of heterochromatin phenotype during senescence has not been fully understood. By using senescent and deep senescent fibroblasts induced by DNA damage harboring the "loss of heterochromatin" phenotype, we observed progressive 3D reorganization of heterochromatin during senescence. Facultative and constitutive heterochromatin marked by H3K27me3 and H3K9me3, respectively, show different alterations. Facultative heterochromatin tends to switch from the repressive B-compartment to the active A-compartment, whereas constitutive heterochromatin shows no significant changes at the compartment level but enhanced interactions between themselves. Both types of heterochromatin show increased chromatin accessibility and gene expression leakage during senescence. Furthermore, increased chromatin accessibility in potential CTCF binding sites accompanies the establishment of novel loops in constitutive heterochromatin. Finally, we also observed aberrant expression of repetitive elements, including LTR (long terminal repeat) and satellite classes. Overall, facultative and constitutive heterochromatin show both similar and distinct multiscale alterations in the 3D map, chromatin accessibility, and gene expression leakage. This study provides an epigenomic map of heterochromatin reorganization during senescence.

The loss of heterochromatin is associated with multiscale three-dimensional genome reorganization and aberrant transcription during cellular senescence.

A core-promoter, approximately from -60 bp upstream to +40 bp downstream of a RNA polymerase (RNAP) II transcription start site (TSS), binds to the preinitiation complex (PIC) and determine the position of TSS. Using position-specific k-tuple feature variables, a quadratic discriminant analysis (QDA) method is shown to be very effective in identifying human core-promoters.

/pdf/a-discrimination-study-of-human-core-promoters-8vjfmydyx3.pdf

Michael Q. Zhang

Papers

The fast rate limit of driven diffusive systems

Correlated evolution of transcription factors and their binding sites

OSCAR: one-class SVM for accurate recognition of cis-elements

The loss of heterochromatin is associated with multiscale three-dimensional genome reorganization and aberrant transcription during cellular senescence.

A discrimination study of human core-promoters.