JustOrthologs: a fast, accurate and user-friendly ortholog identification algorithm.
Ontology highlight
ABSTRACT: Motivation
Orthologous gene identification is fundamental to all aspects of biology. For example, ortholog identification between species can provide functional insights for genes of unknown function and is a necessary step in phylogenetic inference. Currently, most ortholog identification algorithms require all-versus-all BLAST comparisons, which are time-consuming and memory intensive.Results
In contrast to existing approaches, JustOrthologs exploits the conservation of gene structure by using the lengths of coding sequence regions and dinucleotide percentages to identify orthologs. In comparison to OrthoMCL, OMA and OrthoFinder, JustOrthologs decreases ortholog identification runtime by more than 96% and achieves comparable precision and recall scores. The computational speedup allowed us to conduct pairwise comparisons of 1197 complete genomes (780 eukaryotes and 417 archaea). We confirmed gene annotations for 384 120 genes, grouped 1 675 415 genes in previously unreported ortholog groups, and identified 51 429 potentially mislabeled genes across 622 843 ortholog groups.Availability and implementation
JustOrthologs is an open source collaborative software package available in the GitHub repository: https://github.com/ridgelab/JustOrthologs/. All test FASTA files used for comparisons are freely available at https://github.com/ridgelab/JustOrthologs/comparisonFastaFiles/. Reference genomes used in this work are available for download from the NCBI repository: ftp://ftp.ncbi.nih.gov/genomes/.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Miller JB
PROVIDER: S-EPMC6378933 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:SummaryReduced representation bisulfite sequencing (RRBS) is a powerful yet cost-efficient method for studying DNA methylation on a genomic scale. RRBS involves restriction-enzyme digestion, bisulfite conversion and size selection, resulting in DNA sequencing data that require special bioinformatic handling. Here, we describe RRBSMAP, a short-read alignment tool that is designed for handling RRBS data in a user-friendly and scalable way. RRBSMAP uses wildcard alignment, and avoids the need for any preprocessing or post-processing steps. We benchmarked RRBSMAP against a well-validated MAQ-based pipeline for RRBS read alignment and observed similar accuracy but much improved runtime performance, easier handling and better scaling to large sample sets. In summary, RRBSMAP removes bioinformatic hurdles and reduces the computational burden of large-scale epigenome association studies performed with RRBS.Availabilityhttp://rrbsmap.computational-epigenetics.org/ http://code.google.com/p/bsmap/Contactwl1@bcm.tmc.eduSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC3268241 | biostudies-literature
Project description:TILLING (Targeting Induced Local Lesions IN Genomes) is a powerful reverse genetics method in plant functional genomics and breeding to identify mutagenized individuals with improved behavior for a trait of interest. Pooled high throughput sequencing (HTS) of the targeted genes allows efficient identification and sample assignment of variants within genes of interest in hundreds of individuals. Although TILLING has been used successfully in different crops and even applied to natural populations, one of the main issues for a successful TILLING experiment is that most currently available bioinformatics tools for variant detection are not designed to identify mutations with low frequencies in pooled samples or to perform sample identification from variants identified in overlapping pools. Our research group maintains the Next Generation Sequencing Experience Platform (NGSEP), an open source solution for analysis of HTS data. In this manuscript, we present three novel components within NGSEP to facilitate the design and analysis of TILLING experiments: a pooled variants detector, a sample identifier from variants detected in overlapping pools and a simulator of TILLING experiments. A new implementation of the NGSEP calling model for variant detection allows accurate detection of low frequency mutations within pools. The samples identifier implements the process to triangulate the mutations called within overlapping pools in order to assign mutations to single individuals whenever possible. Finally, we developed a complete simulator of TILLING experiments to enable benchmarking of different tools and to facilitate the design of experimental alternatives varying the number of pools and individuals per pool. Simulation experiments based on genes from the common bean genome indicate that NGSEP provides similar accuracy and better efficiency than other tools to perform pooled variants detection. To the best of our knowledge, NGSEP is currently the only tool that generates individual assignments of the mutations discovered from the pooled data. We expect that this development will be of great use for different groups implementing TILLING as an alternative for plant breeding and even to research groups performing pooled sequencing for other applications.
| S-EPMC7886796 | biostudies-literature
Project description:The discovery of the small regulatory RNA populations has changed our vision of cellular regulations. Indeed, loaded on Argonaute proteins they formed ribonucleoprotein complexes that target complementary sequences and achieved widespread silencing mechanisms conserved in most eukaryotes. The recent development of deep sequencing approaches highly contributed to their detection. Small RNA isolation form cells and/or tissues remains a crucial stage to generate robust and relevant sequencing data. In 2006, a novel strategy based on anion-exchange chromatography has been purposed as an alternative to the standard size-isolation purification procedure. However, the eventual biases of such a method have been poorly investigated. Moreover, this strategy not only relies on advanced technical skills and expensive material but is time consuming and requires an elevated starting biological material amount. Using bioinformatic comparative analysis of six independent small RNA-sequencing libraries of Drosophila ovaries, we here demonstrate that anion-exchange chromatography purification prior to small RNA extraction unbiasedly enriches datasets in bona fide reads (small regulatory RNA reads) and depletes endogenous contaminants (ribosomal RNAs and degradation products). The resulting increase of sequencing depth provides a major benefit to study rare populations. We then developed a fast and basic manual procedure to purify loaded small non coding RNAs using anion-exchange chromatography at the bench. We validated the efficiency of this new method and used this strategy to purify small RNAs from various tissues and organisms. We moreover determined that our manual purification increases the output of the previously described anion-exchange chromatography procedure.
2013-06-13 | GSE40748 | GEO
Project description:BackgroundHigh-throughput, image-based screens of cellular responses to genetic or chemical perturbations generate huge numbers of cell images. Automated analysis is required to quantify and compare the effects of these perturbations. However, few of the current freely-available bioimage analysis software tools are optimized for efficient handling of these images. Even fewer of them are designed to transform the phenotypic features measured from these images into discriminative profiles that can reveal biologically meaningful associations among the tested perturbations.ResultsWe present a fast and user-friendly software platform called "cellXpress" to segment cells, measure quantitative features of cellular phenotypes, construct discriminative profiles, and visualize the resulting cell masks and feature values. We have also developed a suite of library functions to load the extracted features for further customizable analysis and visualization under the R computing environment. We systematically compared the processing speed, cell segmentation accuracy, and phenotypic-profile clustering performance of cellXpress to other existing bioimage analysis software packages or algorithms. We found that cellXpress outperforms these existing tools on three different bioimage datasets. We estimate that cellXpress could finish processing a genome-wide gene knockdown image dataset in less than a day on a modern personal desktop computer.ConclusionsThe cellXpress platform is designed to make fast and efficient high-throughput phenotypic profiling more accessible to the wider biological research community. The cellXpress installation packages for 64-bit Windows and Linux, user manual, installation guide, and datasets used in this analysis can be downloaded freely from http://www.cellXpress.org.
| S-EPMC3853218 | biostudies-literature
Project description:Sweet potato (Ipomoea batatas) is one of the ten most important staple crops and provides a livelihood for many people around the globe. To adapt to ever-changing circumstances farmers and breeders need to have access to a broad diversity of germplasm. This study focuses on the development of a cryopreservation protocol that allows the long term storage of different sweet potato cultivars. For this, a droplet vitrification protocol was optimized, comparing several parameters; preculture method (0.3 M sucrose vs no preculture); meristem position (axillary vs apical); plant age (3 to 9 weeks); regeneration medium (MS?+?2.22 µM BA, Hirai and MS); and length of loading solution treatment (20 to 360 min). Two months after cryopreservation, the regeneration rates of the meristems were compared, which resulted in significant differences for the preculture method, meristem position and loading solution. With these new insights an optimized droplet vitrification protocol was developed with the following parameters: use of 3-9 week old axillary meristems, no preculture phase, 20 min LS treatment, 30 min PVS2 treatment, exposure to liquid nitrogen by droplet vitrification, warming treatment in RS for 15 min, 1 day 0.3 M sucrose recuperation culture, 1 month MS?+?2.22 µM BA followed by 1 month of MS cultures. This protocol was subsequently tested on 10 representative accessions resulting in a post cryopreservation regeneration rate of more than 40% for 70% of the tested cultivars, showing that this protocol could be implemented for a large portion of existing sweet potato collections.
| S-EPMC7477159 | biostudies-literature
Project description:MotivationPhylogenies are important for fundamental biological research, but also have numerous applications in biotechnology, agriculture and medicine. Finding the optimal tree under the popular maximum likelihood (ML) criterion is known to be NP-hard. Thus, highly optimized and scalable codes are needed to analyze constantly growing empirical datasets.ResultsWe present RAxML-NG, a from-scratch re-implementation of the established greedy tree search algorithm of RAxML/ExaML. RAxML-NG offers improved accuracy, flexibility, speed, scalability, and usability compared with RAxML/ExaML. On taxon-rich datasets, RAxML-NG typically finds higher-scoring trees than IQTree, an increasingly popular recent tool for ML-based phylogenetic inference (although IQ-Tree shows better stability). Finally, RAxML-NG introduces several new features, such as the detection of terraces in tree space and the recently introduced transfer bootstrap support metric.Availability and implementationThe code is available under GNU GPL at https://github.com/amkozlov/raxml-ng. RAxML-NG web service (maintained by Vital-IT) is available at https://raxml-ng.vital-it.ch/.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC6821337 | biostudies-literature
Project description:The discovery of the small regulatory RNA populations has changed our vision of cellular regulations. Indeed, loaded on Argonaute proteins they formed ribonucleoprotein complexes that target complementary sequences and achieved widespread silencing mechanisms conserved in most eukaryotes. The recent development of deep sequencing approaches highly contributed to their detection. Small RNA isolation form cells and/or tissues remains a crucial stage to generate robust and relevant sequencing data. In 2006, a novel strategy based on anion-exchange chromatography has been purposed as an alternative to the standard size-isolation purification procedure. However, the eventual biases of such a method have been poorly investigated. Moreover, this strategy not only relies on advanced technical skills and expensive material but is time consuming and requires an elevated starting biological material amount. Using bioinformatic comparative analysis of six independent small RNA-sequencing libraries of Drosophila ovaries, we here demonstrate that anion-exchange chromatography purification prior to small RNA extraction unbiasedly enriches datasets in bona fide reads (small regulatory RNA reads) and depletes endogenous contaminants (ribosomal RNAs and degradation products). The resulting increase of sequencing depth provides a major benefit to study rare populations. We then developed a fast and basic manual procedure to purify loaded small non coding RNAs using anion-exchange chromatography at the bench. We validated the efficiency of this new method and used this strategy to purify small RNAs from various tissues and organisms. We moreover determined that our manual purification increases the output of the previously described anion-exchange chromatography procedure. Comparison of small regulatory RNA populations obtained after three different small RNA purification procedures
2013-06-13 | E-GEOD-40748 | biostudies-arrayexpress
Project description:Supertree methods merge a set of overlapping phylogenetic trees into a supertree containing all taxa of the input trees. The challenge in supertree reconstruction is the way of dealing with conflicting information in the input trees. Many different algorithms for different objective functions have been suggested to resolve these conflicts. In particular, there exist methods based on encoding the source trees in a matrix, where the supertree is constructed applying a local search heuristic to optimize the respective objective function. We present a novel heuristic supertree algorithm called Bad Clade Deletion (BCD) supertrees. It uses minimum cuts to delete a locally minimal number of columns from such a matrix representation so that it is compatible. This is the complement problem to Matrix Representation with Compatibility (Maximum Split Fit). Our algorithm has guaranteed polynomial worst-case running time and performs swiftly in practice. Different from local search heuristics, it guarantees to return the directed perfect phylogeny for the input matrix, corresponding to the parent tree of the input trees, if one exists. Comparing supertrees to model trees for simulated data, BCD shows a better accuracy (F1 score) than the state-of-the-art algorithms SuperFine (up to 3%) and Matrix Representation with Parsimony (up to 7%); at the same time, BCD is up to 7 times faster than SuperFine, and up to 600 times faster than Matrix Representation with Parsimony. Finally, using the BCD supertree as a starting tree for a combined Maximum Likelihood analysis using RAxML, we reach significantly improved accuracy (1% higher F1 score) and running time (1.7-fold speedup).
| S-EPMC5850620 | biostudies-literature
Project description:BackgroundWhile RNA is often created from linear splicing during transcription, recent studies have found that non-canonical splicing sometimes occurs. Non-canonical splicing joins 3' and 5' and forms the so-called circular RNA. It is now believed that circular RNA plays important biological roles such as affecting susceptibility of some diseases. During the past several years, multiple experimental methods have been developed to enrich circular RNA while degrade linear RNA. Although several useful software tools for circular RNA detection have been developed as well, these tools are based on reads mapping may miss many circular RNA. Also, existing tools are slow for large data due to their dependence on reads mapping.MethodIn this paper, we present a new computational approach, named CircMarker, based on k-mers rather than reads mapping for circular RNA detection. CircMarker takes advantage of transcriptome annotation files to create the k-mer table for circular RNA detection.ResultsEmpirical results show that CircMarker outperforms existing tools in circular RNA detection on accuracy and efficiency in many simulated and real datasets.ConclusionsWe develop a new circular RNA detection method called CircMarker based on k-mer analysis. Our results on both simulation data and real data demonstrate that CircMarker runs much faster and can find more circular RNA with higher consensus-based sensitivity and high accuracy ratio compared with existing tools.
| S-EPMC6101086 | biostudies-literature
Project description:In quantitative magnetic resonance mapping, the variable flip angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution weighted images in a clinically feasible time. Fast, linear methods that estimate maps from these weighted images have been proposed, such as DESPOT1 and iterative re-weighted linear least squares. More accurate, non-linear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this paper, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared with conventional gradient-based NLLS estimators while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than the efficient implementations of the variable projection (VARPRO) method. Furthermore, NOVIFAST is shown to be robust against initialization.
| S-EPMC6277233 | biostudies-literature