Scalable heuristics for design of DNA probe arrays.

Design of DNA arrays for very large-scale immobilized polymer synthesis (VLSIPS) (Fodor et al., 1991) seeks to minimize effects of unintended illumination during mask exposure steps. Hannenhalli et al. (2002) formulate this requirement as the Border Minimization Problem and give an algorithm for placement of probes at array sites under the assumption that the array synthesis is synchronous; i.e., nucleotides are synthesized in a periodic sequence (ACGT)(k) and every probe grows by exactly one nucleotide with every group of four masks. Drawing on the analogy with VLSI placement, in this paper we describe and experimentally validate the engineering of several scalable, high-quality placement heuristics for both synchronous and asynchronous DNA array design. We give empirical results on both randomly generated and industry test cases confirming the scalability and improved solution quality enjoyed by our methods. In general, our techniques improve on state-of-the-art industrial results by over 4% and surpass academically published results by up to 35%. Finally, we give lower bounds that offer insights into the amount of available further improvements.     

Correcting base-assignment errors in repeat regions of shotgun assembly

Accurate base-assignment in repeat regions of a whole genome shotgun assembly is an unsolved problem. Since reads in repeat regions cannot be easily attributed to a unique location in the genome, current assemblers may place these reads arbitrarily. As a result, the base-assignment error rate in repeats is likely to be much higher than that in the rest of the genome. We developed an iterative algorithm, EULER-AIR, that is able to correct base-assignment errors in finished genome sequences in public databases. The Wolbachia genome is among the best finished genomes. Using this genome project as an example, we demonstrated that EULER-AIR can 1) discover and correct base-assignment errors, 2) provide accurate read assignments, 3) utilize finishing reads for accurate base-assignment, and 4) provide guidance for designing finishing experiments. In the genome of Wolbachia, EULER-AIR found 16 positions with ambiguous base-assignment and two positions with erroneous bases. Besides Wolbachia, many other genome sequencing projects have significantly fewer finishing reads and, hence, are likely to contain more base-assignment errors in repeats. We demonstrate that EULER-AIR is a software tool that can be used to find and correct base-assignment errors in a genome assembly project     

Are There Rearrangement Hotspots in the Human Genome?

In a landmark paper, Nadeau and Taylor [18] formulated the random breakage model (RBM) of chromosome evolution that postulates that there are no rearrangement hotspots in the human genome. In the next two decades, numerous studies with progressively increasing levels of resolution made RBM the de facto theory of chromosome evolution. Despite the fact that RBM had prophetic prediction power, it was recently refuted by Pevzner and Tesler [4], who introduced the fragile breakage model (FBM), postulating that the human genome is a mosaic of solid regions (with low propensity for rearrangements) and fragile regions (rearrangement hotspots). However, the rebuttal of RBM caused a controversy and led to a split among researchers studying genome evolution. In particular, it remains unclear whether some complex rearrangements (e.g., transpositions) can create an appearance of rearrangement hotspots. We contribute to the ongoing debate by analyzing multi-break rearrangements that break a genome into multiple fragments and further glue them together in a new order. In particular, we demonstrate that (1) even if transpositions were a dominant force in mammalian evolution, the arguments in favor of FBM still stand, and (2) the "gene deletion" argument against FBM is flawed.     

Systematics of the Platyrrhinus helleri species complex (Chiroptera: Phyllostomidae), with descriptions of two new species

Platyrrhinus is a diverse genus of small to large phyllostomid bats characterized by a comparatively narrow uropatagium thickly fringed with hair, a white dorsal stripe, comparatively large inner upper incisors that are convergent at the tips, and three upper and three lower molars. Eighteen species are currently recognized, the majority occurring in the Andes. Molecular, morphological, and morphometric analyses of specimens formerly identified as Platyrrhinus helleri support recognition of Platyrrhinus incarum as a separate species and reveal the presence of two species from western and northern South America that we describe herein as new ( Platyrrhinus angustirostris sp. nov. from eastern Colombia and Ecuador, north‐eastern Peru, and Venezuela and Platyrrhinus fusciventris sp. nov. from Guyana, Suriname, French Guiana, Trinidad and Tobago, northern Brazil, eastern Ecuador, and southern Venezuela). These two new species are sister taxa and, in turn, sister to Platyrrhinus incarum. © 2010 The Linnean Society of London, Zoological Journal of the Linnean Society, 2010, 159 , 785–812.     

Linguistics of Nucleotide Sequences II: Stationary Words in Genetic Texts and the Zonal Structure of DNA

Abstract

Words are irregularly distributed in genetic texts. The analysis of this irregularity leads to the notion of stationary and non-stationary words. The polyW and polyS tracts are shown to be the most non-stationary words in genetic texts (here W-(AT), S-{G,C}, a polyW tract is a sequence of AT nucleotides and a polyS tract is a sequence of G,C nucleotides. The distribution of stationary words suggests a method for partitioning DNA into zones. The zones obtained in the case of the phage are interpreted in the light of the Dowe hypothesis of the modular structure of bacteriophage genomes.     

The phenoptosis problem: What is causing the death of an organism? Lessons from acute kidney injury

Programmed execution of various cells and intracellular structures is hypothesized to be not the only example of elimination of biological systems — the general mechanism can also involve programmed execution of organs and organisms. Modern rating of programmed cell death mechanisms includes 13 mechanistic types. As for some types, the mechanism of actuation and manifestation of cell execution has been basically elucidated, while the causes and intermediate steps of the process of fatal failure of organs and organisms remain unknown. The analysis of deaths resulting from a sudden heart arrest or multiple organ failure and other acute and chronic pathologies leads to the conclusion of a special role of mitochondria and oxidative stress activating the immune system. Possible mechanisms of mitochondria-mediated induction of the signaling cascades involved in organ failure and death of the organism are discussed. These mechanisms include generation of reactive oxygen species and damage-associated molecular patterns in mitochondria. Some examples of renal failure-induced deaths are presented with mechanisms and settings determined by some hypothetical super system rather than by the kidneys themselves. This system plays the key role in the process of physiological senescence and termination of an organism. The facts presented suggest that it is the immune system involved in mitochondrial signaling that can act as the system responsible for the organism’s death.     

How many signal peptides are there in bacteria?

Over the last 5 years proteogenomics (using mass spectroscopy to identify proteins predicted from genomic sequences) has emerged as a promising approach to the high‐throughput identification of protein N‐termini, which remains a problem in genome annotation. Comparison of the experimentally determined N‐termini with those predicted by sequence analysis tools allows identification of the signal peptides and therefore conclusions on the cytoplasmic or extracytoplasmic (periplasmic or extracellular) localization of the respective proteins. We present here the results of a proteogenomic study of the signal peptides in Escherichia coli K‐12 and compare its results with the available experimental data and predictions by such software tools as SignalP and Phobius. A single proteogenomics experiment recovered more than a third of all signal peptides that had been experimentally determined during the past three decades and confirmed at least 31 additional signal peptides, mostly in the known exported proteins, which had been previously predicted but not validated. The filtering of putative signal peptides for the peptide length and the presence of an eight‐residue hydrophobic patch and a typical signal peptidase cleavage site proved sufficient to eliminate the false‐positive hits. Surprisingly, the results of this proteogenomics study, as well as a re‐analysis of the E. coli genome with the latest version of SignalP program, show that the fraction of proteins containing signal peptides is only about 10%, or half of previous estimates.     

De novo peptide sequencing via tandem mass spectrometry.

Peptide sequencing via tandem mass spectrometry (MS/MS) is one of the most powerful tools in proteomics for identifying proteins. Because complete genome sequences are accumulating rapidly, the recent trend in interpretation of MS/MS spectra has been database search. However, de novo MS/MS spectral interpretation remains an open problem typically involving manual interpretation by expert mass spectrometrists. We have developed a new algorithm, SHERENGA, for de novo interpretation that automatically learns fragment ion types and intensity thresholds from a collection of test spectra generated from any type of mass spectrometer. The test data are used to construct optimal path scoring in the graph representations of MS/MS spectra. A ranked list of high scoring paths corresponds to potential peptide sequences. SHERENGA is most useful for interpreting sequences of peptides resulting from unknown proteins and for validating the results of database search algorithms in fully automated, high-throughput peptide sequencing.     

Immune Response and the B Blood Group Locus in Chickens

Chickens of three blood group genotypes, B1B1, B1B2 and B1B19, were compared in their ability to produce antibodies. The B locus is the major histocompatibility locus in this species. Homozygous B1 pullets had significantly higher adult mortality than did the B1 heterozygotes. In two separate experiments the B1B1 females produced significantly fewer antibodies to Salmonella pullorum than did the B1 heterozygotes. Also the B1B1 pullets responded with lower antibody titers following immunization with ferritin, bovine serum albumin and parainfluenza-3 virus, although the differences were not significant. The results of this study suggest that an immune response gene is associated with the major histocompatibility locus in chickens, paralleling the H-2 locus and Ir genes in mice.