MotifCut: regulatory motifs finding with maximum density subgraphs

MOTIVATION: DNA motif finding is one of the core problems in computational biology, for which several probabilistic and discrete approaches have been developed. Most existing methods formulate motif finding as an intractable optimization problem and rely either on expectation maximization (EM) or on local heuristic searches. Another challenge is the choice of motif model: simpler models such as the position-specific scoring matrix (PSSM) impose biologically unrealistic assumptions such as independence of the motif positions, while more involved models are harder to parametrize and learn. RESULTS: We present MotifCut, a graph-theoretic approach to motif finding leading to a convex optimization problem with a polynomial time solution. We build a graph where the vertices represent all k-mers in the input sequences, and edges represent pairwise k-mer similarity. In this graph, we search for a motif as the maximum density subgraph, which is a set of k-mers that exhibit a large number of pairwise similarities. Our formulation does not make strong assumptions regarding the structure of the motif and in practice both motifs that fit well the PSSM model, and those that exhibit strong dependencies between position pairs are found as dense subgraphs. We benchmark MotifCut on both synthetic and real yeast motifs, and find that it compares favorably to existing popular methods. The ability of MotifCut to detect motifs appears to scale well with increasing input size. Moreover, the motifs we discover are different from those discovered by the other methods. AVAILABILITY: MotifCut server and other materials can be found at motifcut.stanford.edu.     

R type anion channel: A multifunctional channel seeking its molecular identity

Plant genomes code for channels involved in the transport of cations, anions and uncharged molecules through membranes. Although the molecular identity of channels for cations and uncharged molecules has progressed rapidly in the recent years, the molecular identity of anion channels has lagged behind. Electrophysiological studies have identified S-type (slow) and R-type (rapid) anion channels. In this brief review, we summarize the proposed functions of the R-type anion channels which, like the S-type, were first characterized by electrophysiology over 20 years ago, but unlike the S-type, have still yet to be cloned. We show that the R-type channel can play multiple roles.Key words: R-type anion channel, nitrate, sulphate, guard cell, action potentialAnion channels play a central role in signal transduction, nutrient transport and cell turgor regulation.¹ By far, their function was particularly well investigated in the guard cells of stomata using a combination of electrophysiological, pharmacological and genetic tools. In this system, anion channel activation was shown to be one of the limiting steps in the loss of cell turgor leading to stomatal closure.² In algal cells, anion channels were shown to contribute to membrane excitability through the generation of action potential.^1,3With the burst of molecular biology in the nineties, the genes coding for plant ion channels started to be unveiled. The first channel gene to be cloned in plant was the shaker-like potassium channel identified in a yeast functional expression screen.^4,5 More than ten years later, TaALMT1 and AtCLCa were characterized as the first members of two important anion channel families.^6,7 This growing group of newly identified channels, accounting for electrophysiological activity described long ago, includes the MSLs anion selective mechanosensitive channels.⁸ Recently, the well known S-type channel has been finally recognized to be encoded by members of the SLAC1 (and other SLAH) family (Slow Anion Channel-Associated 1).⁹ In agreement with electrophysiological data,^10–13 it requires phosphorylation by a Protein Kinase in order to be functional.^14,15 In contrast, the molecular identity of the R-type anion channel remains unknown. Therefore, this candidate, which has been functionally known since twenty years, remains the next challenge for plant channel physiologists.     

A new computational efficient approach for trabecular bone analysis using beam models generated with skeletonized graph technique

Micro-finite element (FE) analysis is a well established technique for the evaluation of the elastic properties of trabecular bone, but is limited in its application due to the large number of elements that it requires to represent the complex internal structure of the bone. In this paper, we present an alternative FE approach that makes use of a recently developed 3D-Line Skeleton Graph Analysis (LSGA) technique to represent the complex internal structure of trabecular bone as a network of simple straight beam elements in which the beams are assigned geometrical properties of the trabeculae that they represent. Since an enormous reduction of cputime can be obtained with this beam modeling approach, ranging from approximately 1,200 to 3,600 for the problems investigated here, we think that the FE modeling technique that we introduced could potentially constitute an interesting alternative for the evaluation of the elastic mechanical properties of trabecular bone.     

Molecular phylogeny of the Haplosporidia based on two independent gene sequences

The phylogenetic position of the Haplosporidia has confounded taxonomists for more than a century because of the unique morphology of these parasites. We collected DNA sequence data for small subunit (SSU) ribosomal RNA and actin genes from haplosporidians and other protists for conducting molecular phylogenetic analyses to help elucidate relationships of taxa within the group, as well as placement of this group among Eukaryota. Analyses were conducted using DNA sequence data from more than 100 eukaryotic taxa with various combinations of data sets including nucleotide sequence data for each gene separately and combined, as well as SSU ribosomal DNA data combined with translated actin amino acids. In almost all analyses, the Haplosporidia was sister to the Cercozoa with moderate bootstrap and jackknife support. Analysis with actin amino acid sequences alone grouped haplosporidians with the foraminiferans and cercozoans. The haplosporidians Minchinia and Urosporidium were found to be monophyletic, whereas Haplosporidium was paraphyletic. "Microcell" parasites, Bonamia spp. and Mikrocytos roughleyi, were sister to Minchinia, the most derived genus, with Haplosporidium falling between the "microcells" and the more basal Urosporidium. Two recently discovered parasites, one from abalone in New Zealand and another from spot prawns in British Columbia, fell at the base of the Haplosporidia with very strong support, indicating a taxonomic affinity to this group.     

Randomized sequence databases for tandem mass spectrometry peptide and protein identification

Tandem mass spectrometry (MS/MS) combined with database searching is currently the most widely used method for high-throughput peptide and protein identification. Many different algorithms, scoring criteria, and statistical models have been used to identify peptides and proteins in complex biological samples, and many studies, including our own, describe the accuracy of these identifications, using at best generic terms such as "high confidence." False positive identification rates for these criteria can vary substantially with changing organisms under study, growth conditions, sequence databases, experimental protocols, and instrumentation; therefore, study-specific methods are needed to estimate the accuracy (false positive rates) of these peptide and protein identifications. We present and evaluate methods for estimating false positive identification rates based on searches of randomized databases (reversed and reshuffled). We examine the use of separate searches of a forward then a randomized database and combined searches of a randomized database appended to a forward sequence database. Estimated error rates from randomized database searches are first compared against actual error rates from MS/MS runs of known protein standards. These methods are then applied to biological samples of the model microorganism Shewanella oneidensis strain MR-1. Based on the results obtained in this study, we recommend the use of use of combined searches of a reshuffled database appended to a forward sequence database as a means providing quantitative estimates of false positive identification rates of peptides and proteins. This will allow researchers to set criteria and thresholds to achieve a desired error rate and provide the scientific community with direct and quantifiable measures of peptide and protein identification accuracy as opposed to vague assessments such as "high confidence."     

An algorithm for automatic evaluation of the spot quality in two-color DNA microarray experiments

Background  

Although DNA microarray technologies are very powerful for the simultaneous quantitative characterization of thousands of genes, the quality of the obtained experimental data is often far from ideal. The measured microarrays images represent a regular collection of spots, and the intensity of light at each spot is proportional to the DNA copy number or to the expression level of the gene whose DNA clone is spotted. Spot quality control is an essential part of microarray image analysis, which must be carried out at the level of individual spot identification. The problem is difficult to formalize due to the diversity of instrumental and biological factors that can influence the result.     

Towards understanding the first genome sequence of a crenarchaeon by genome annotation using clusters of orthologous groups of proteins (COGs)

Background  

Standard archival sequence databases have not been designed as tools for genome annotation and are far from being optimal for this purpose. We used the database of Clusters of Orthologous Groups of proteins (COGs) to reannotate the genomes of two archaea, Aeropyrum pernix, the first member of the Crenarchaea to be sequenced, and Pyrococcus abyssi.     

Encapsulation enhances protoplast fusant stability

A barrier to cost-efficient biomanufacturing is the instability of engineered genetic elements, such as plasmids. Instability can also manifest at the whole-genome level, when fungal dikaryons revert to parental species due to nuclear segregation during cell division. Here, we show that by encapsulating Saccharomyces cerevisiae-Pichia stipitis dikaryons in an alginate matrix, we can limit cell division and preserve their expanded metabolic capabilities. As a proxy to cellulosic ethanol production, we tested the capacity of such cells to carry out ethanologenic fermentation of glucose and xylose, examining substrate use, ploidy, and cell viability in relation to planktonic fusants, as well as in relation to planktonic and encapsulated cell cultures consisting of mixtures of these species. Glucose and xylose consumption and ethanol production by encapsulated dikaryons were greater than planktonic controls. Simultaneous co-fermentation did not occur; rather the order and kinetics of glucose and xylose catabolism by encapsulated dikaryons were similar to cultures where the two species were encapsulated together. Over repeated cycles of fed-batch culture, encapsulated S. cerevisiae-P. stipitis fusants exhibited a dramatic increase in genomic stability, relative to planktonic fusants. Encapsulation also increased the stability of antibiotic-resistance plasmids used to mark each species and preserved a fixed ratio of S. cerevisiae to P. stipitis cells in mixed cultures. Our data demonstrate how encapsulating cells in an extracellular matrix restricts cell division and, thereby, preserves the stability and biological activity of entities ranging from genomes to plasmids to mixed populations, each of which can be essential to cost-efficient biomanufacturing.     

ProServer: a simple, extensible Perl DAS server

SUMMARY: The increasing size and complexity of biological databases has led to a growing trend to federate rather than duplicate them. In order to share data between federated databases, protocols for the exchange mechanism must be developed. One such data exchange protocol that is widely used is the Distributed Annotation System (DAS). For example, DAS has enabled small experimental groups to integrate their data into the Ensembl genome browser. We have developed ProServer, a simple, lightweight, Perl-based DAS server that does not depend on a separate HTTP server. The ProServer package is easily extensible, allowing data to be served from almost any underlying data model. Recent additions to the DAS protocol have enabled both structure and alignment (sequence and structural) data to be exchanged. ProServer allows both of these data types to be served. AVAILABILITY: ProServer can be downloaded from http://www.sanger.ac.uk/proserver/ or CPAN http://search.cpan.org/~rpettett/. Details on the system requirements and installation of ProServer can be found at http://www.sanger.ac.uk/proserver/.     

Persistence and decontamination of Bacillus atrophaeus subsp. globigii spores on corroded iron in a model drinking water system

Persistence of Bacillus atrophaeus subsp. globigii spores on corroded iron coupons in drinking water was studied using a biofilm annular reactor. Spores were inoculated at 10(6) CFU/ml in the dechlorinated reactor bulk water. The dechlorination allowed for observation of the effects of hydraulic shear and biofilm sloughing on persistence. Approximately 50% of the spores initially adhered to the corroded iron surface were not detected after 1 month. Addition of a stable 10 mg/liter free chlorine residual after 1 month led to a 2-log(10) reduction of adhered B. atrophaeus subsp. globigii, but levels on the coupons quickly stabilized thereafter. Increasing the free chlorine concentration to 25 or 70 mg/liter had no additional effect on inactivation. B. atrophaeus subsp. globigii spores injected in the presence of a typical distribution system chlorine residual (approximately 0.75 mg/liter) resulted in a steady reduction of adhered B. atrophaeus subsp. globigii over 1 month, but levels on the coupons eventually stabilized. Adding elevated chlorine levels (10, 25, and 70 mg/liter) after 1 month had no effect on the rate of inactivation. Decontamination with elevated free chlorine levels immediately after spore injection resulted in a 3-log(10) reduction within 2 weeks, but the rate of inactivation leveled off afterward. This indicates that free chlorine did not reach portions of the corroded iron surface where B. atrophaeus subsp. globigii spores had adhered. B. atrophaeus subsp. globigii spores are capable of persisting for an extended time in the presence of high levels of free chlorine.