Mount Kisco, NY, United States
Mount Kisco, NY, United States

Time filter

Source Type

Engstrom P.G.,European Bioinformatics Institute | Engstrom P.G.,University of Stockholm | Steijger T.,European Bioinformatics Institute | Sipos B.,European Bioinformatics Institute | And 39 more authors.
Nature Methods | Year: 2013

High-throughput RNA sequencing is an increasingly accessible method for studying gene structure and activity on a genome-wide scale. A critical step in RNA-seq data analysis is the alignment of partial transcript reads to a reference genome sequence. To assess the performance of current mapping software, we invited developers of RNA-seq aligners to process four large human and mouse RNA-seq data sets. In total, we compared 26 mapping protocols based on 11 programs and pipelines and found major performance differences between methods on numerous benchmarks, including alignment yield, basewise accuracy, mismatch and gap placement, exon junction discovery and suitability of alignments for transcript reconstruction. We observed concordant results on real and simulated RNA-seq data, confirming the relevance of the metrics employed. Future developments in RNA-seq alignment methods would benefit from improved placement of multimapped reads, balanced utilization of existing gene annotation and a reduced false discovery rate for splice junctions. © 2013 Nature America, Inc.


Earl D.,University of California at Santa Cruz | Nguyen N.,University of California at Santa Cruz | Hickey G.,McGill University | Harris R.S.,Pennsylvania State University | And 30 more authors.
Genome Research | Year: 2014

Multiple sequence alignments (MSAs) are a prerequisite for a wide variety of evolutionary analyses. Published assessments and benchmark data sets for protein and, to a lesser extent, global nucleotide MSAs are available, but less effort has been made to establish benchmarks in the more general problem of whole-genome alignment (WGA). Using the same model as the successful Assemblathon competitions, we organized a competitive evaluation in which teams submitted their alignments and then assessments were performed collectively after all the submissions were received. Three data sets were used: Two were simulated and based on primate and mammalian phylogenies, and one was comprised of 20 real fly genomes. In total, 35 submissions were assessed, submitted by 10 teams using 12 different alignment pipelines. We found agreement between independent simulation-based and statistical assessments, indicating that there are substantial accuracy differences between contemporary alignment tools. We saw considerable differences in the alignment quality of differently annotated regions and found that few tools aligned the duplications analyzed. We found that many tools worked well at shorter evolutionary distances, but fewer performed competitively at longer distances. We provide all data sets, submissions, and assessment programs for further study and provide, as a resource for future benchmarking, a convenient repository of code and data for reproducing the simulation assessments. © 2014 Earl et al.


Earl D.,University of California at Santa Cruz | Bradnam K.,University of California at Davis | St. John J.,University of California at Santa Cruz | Darling A.,University of California at Davis | And 83 more authors.
Genome Research | Year: 2011

Low-cost short read sequencing technology has revolutionized genomics, though it is only just becoming practical for the high-quality de novo assembly of a novel large genome. We describe the Assemblathon 1 competition, which aimed to comprehensively assess the state of the art in de novo assembly methods when applied to current sequencing technologies. In a collaborative effort, teams were asked to assemble a simulated Illumina HiSeq data set of an unknown, simulated diploid genome. A total of 41 assemblies from 17 different groups were received. Novel haplotype aware assessments of coverage, contiguity, structure, base calling, and copy number were made. We establish that within this benchmark: (1) It is possible to assemble the genome to a high level of coverage and accuracy, and that (2) large differences exist between the assemblies, suggesting room for further improvements in current methods. The simulated benchmark, including the correct answer, the assemblies, and the code that was used to evaluate the assemblies is now public and freely available from http://www.assemblathon.org/. © 2011 by Cold Spring Harbor Laboratory Press.


PubMed | RAS Institute of Cytology and Genetics, Novosibirsk State University, Softberry Inc. and Russian Academy of Medical Sciences
Type: | Journal: Oncotarget | Year: 2016

Using the ability of poorly differentiated cells to natively internalize fragments of extracellular double-stranded DNA as a marker, we isolated a tumorigenic subpopulation present in Krebs-2 ascites that demonstrated the features of tumor-inducing cancer stem cells. Having combined TAMRA-labeled DNA probe and the power of RNA-seq technology, we identified a set of 168 genes specifically expressed in TAMRA-positive cells (tumor-initiating stem cells), these genes remaining silent in TAMRA-negative cancer cells. TAMRA+ cells displayed gene expression signatures characteristic of both stem cells and cancer cells. The observed expression differences between TAMRA+ and TAMRA- cells were validated by Real Time PCR. The results obtained corroborated the biological data that TAMRA+ murine Krebs-2 tumor cells are tumor-initiating stem cells. The approach developed can be applied to profile any poorly differentiated cell types that are capable of immanent internalization of double-stranded DNA.


PubMed | Howard Hughes Medical Institute, McGill University, Lawrence Berkeley National Laboratory, University of Technology, Sydney and 9 more.
Type: Journal Article | Journal: Genome research | Year: 2014

Multiple sequence alignments (MSAs) are a prerequisite for a wide variety of evolutionary analyses. Published assessments and benchmark data sets for protein and, to a lesser extent, global nucleotide MSAs are available, but less effort has been made to establish benchmarks in the more general problem of whole-genome alignment (WGA). Using the same model as the successful Assemblathon competitions, we organized a competitive evaluation in which teams submitted their alignments and then assessments were performed collectively after all the submissions were received. Three data sets were used: Two were simulated and based on primate and mammalian phylogenies, and one was comprised of 20 real fly genomes. In total, 35 submissions were assessed, submitted by 10 teams using 12 different alignment pipelines. We found agreement between independent simulation-based and statistical assessments, indicating that there are substantial accuracy differences between contemporary alignment tools. We saw considerable differences in the alignment quality of differently annotated regions and found that few tools aligned the duplications analyzed. We found that many tools worked well at shorter evolutionary distances, but fewer performed competitively at longer distances. We provide all data sets, submissions, and assessment programs for further study and provide, as a resource for future benchmarking, a convenient repository of code and data for reproducing the simulation assessments.


Shahmuradov I.A.,Computer | Shahmuradov I.A.,Azerbaijan Academy of Sciences | Solovyev V.V.,Computer | Solovyev V.V.,Softberry Inc.
Bioinformatics | Year: 2014

Summary: Gene transcription is mostly conducted through interactions of various transcription factors and their binding sites on DNA (regulatory elements, REs). Today, we are still far from understanding the real regulatory content of promoter regions. Computer methods for identification of REs remain a widely used tool for studying and understanding transcriptional regulation mechanisms. The Nsite, NsiteH and NsiteM programs perform searches for statistically significant (non-random) motifs of known human, animal and plant one-box and composite REs in a single genomic sequence, in a pair of aligned homologous sequences and in a set of functionally related sequences, respectively. Availability and implementation: Pre-compiled executables built under commonly used operating systems are available for download by visiting http://www.molquest.kaust.edu.sa and http://www.softberry.com. © The Author 2015. Published by Oxford University Press.


Kufareva I.,University of California at San Diego | Rueda M.,University of California at San Diego | Katritch V.,University of California at San Diego | Stevens R.C.,Scripps Research Institute | And 105 more authors.
Structure | Year: 2011

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe. © 2011 Elsevier Ltd. All rights reserved.


PubMed | International Rice Research Institute, Softberry Inc., Lawrence Berkeley National Laboratory and University of Arizona
Type: Journal Article | Journal: Nucleic acids research | Year: 2016

We describe updates to the Rice SNP-Seek Database since its first release. We ran a new SNP-calling pipeline followed by filtering that resulted in complete, base, filtered and core SNP datasets. Besides the Nipponbare reference genome, the pipeline was run on genome assemblies of IR 64, 93-11, DJ 123 and Kasalath. New genotype query and display features are added for reference assemblies, SNP datasets and indels. JBrowse now displays BAM, VCF and other annotation tracks, the additional genome assemblies and an embedded VISTA genome comparison viewer. Middleware is redesigned for improved performance by using a hybrid of HDF5 and RDMS for genotype storage. Query modules for genotypes, varieties and genes are improved to handle various constraints. An integrated list manager allows the user to pass query parameters for further analysis. The SNP Annotator adds traits, ontology terms, effects and interactions to markers in a list. Web-service calls were implemented to access most data. These features enable seamless querying of SNP-Seek across various biological entities, a step toward semi-automated gene-trait association discovery. URL: http://snp-seek.irri.org.

Loading Softberry Inc. collaborators
Loading Softberry Inc. collaborators