<![CDATA[BOL: Related items]]> https://bioinformaticsonline.com/related/34221?offset=380 https://bioinformaticsonline.com/blog/view/42972/list-of-bioinformatics-workflow-management-tools Sat, 20 Mar 2021 00:15:25 -0500 https://bioinformaticsonline.com/blog/view/42972/list-of-bioinformatics-workflow-management-tools <![CDATA[List of bioinformatics workflow management tools !]]> Here are list of Workflow Managers
  • BigDataScript – A cross-system scripting language for working with big data pipelines in computer systems of different sizes and capabilities. [ paper-2014 | web ]
  • Bpipe – A small language for defining pipeline stages and linking them together to make pipelines. [ web ]
  • Common Workflow Language – a specification for describing analysis workflows and tools that are portable and scalable across a variety of software and hardware environments, from workstations to cluster, cloud, and high performance computing (HPC) environments. [ web ]
  • Cromwell – A Workflow Management System geared towards scientific workflows. [ web ]
  • Galaxy – a popular open-source, web-based platform for data intensive biomedical research. Has several features, from data analysis to workflow management to visualization tools. [ paper-2018 | web ]
  • Nextflow (recommended) – A fluent DSL modelled around the UNIX pipe concept, that simplifies writing parallel and scalable pipelines in a portable manner. [ paper-2018 | web ]
  • Ruffus – Computation Pipeline library for python widely used in science and bioinformatics. [ paper-2010 | web ]
  • SeqWare – Hadoop Oozie-based workflow system focused on genomics data analysis in cloud environments. [ paper-2010 | web ]
  • Snakemake – A workflow management system in Python that aims to reduce the complexity of creating workflows by providing a fast and comfortable execution environment. [ paper-2018 | web ]
  • Workflow Descriptor Language – Workflow standard developed by the Broad. [ web ]
]]> Rahul Nayak https://bioinformaticsonline.com/bookmarks/view/35057/ectools-long-read-correction-and-other-correction-tools Fri, 05 Jan 2018 04:02:22 -0600 https://bioinformaticsonline.com/bookmarks/view/35057/ectools-long-read-correction-and-other-correction-tools <![CDATA[ECTOOLS: Long Read Correction and other Correction tools]]> Long Read Correction and other Correction tools

This package is a loose collection of scripts. To run the correction
routine see the section below. Descriptions of the other scripts
are at the bottom of this file.

Contact: gurtowsk@cshl.edu

In short, the correction algorithm takes as input the unitigs from a short read assembly and uses them to correct long read data. More background information for the algorithm can be found:
http://schatzlab.cshl.edu/presentations/2013-06-18.PBUserMeeting.pdf

Address of the bookmark: https://github.com/jgurtowski/ectools

]]> Jit https://bioinformaticsonline.com/bookmarks/view/36583/eugi-a-novel-resource-for-studying-genomic-islands-to-facilitate-horizontal-gene-transfer-detection-in-eukaryotes Sat, 12 May 2018 07:26:59 -0500 https://bioinformaticsonline.com/bookmarks/view/36583/eugi-a-novel-resource-for-studying-genomic-islands-to-facilitate-horizontal-gene-transfer-detection-in-eukaryotes <![CDATA[EuGI: a novel resource for studying genomic islands to facilitate horizontal gene transfer detection in eukaryotes]]> SWGIS v2.0 along with the EuGI database, which houses GIs identified in 66 different eukaryotic species, and the EuGI web-resource, provide the first comprehensive resource for studying HGT in eukaryotes.

https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-018-4724-8

Address of the bookmark: https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-018-4724-8

]]> Surabhi Chaudhary https://bioinformaticsonline.com/blog/view/38238/list-of-motif-discovery-tools Tue, 20 Nov 2018 03:54:26 -0600 https://bioinformaticsonline.com/blog/view/38238/list-of-motif-discovery-tools <![CDATA[List of motif discovery tools !]]>
In genetics, a sequence motif is a nucleotide or amino-acid sequence pattern that is widespread and has, or is conjectured to have, a biological significance. For proteins, a sequence motif is distinguished from a structural motif, a motif formed by the three-dimensional arrangement of amino acids which may not be adjacent.
 
Following are the list of tools for motif discovery:
 

Perform secondary structure predictions on protein sequences.

Find binding specificity information about DNA-protein complexes.

Find information about the binding specificity of DNA-binding proteins.

Use this web-based sequence motif visualization system to display sequence motif information in its appropriate three-dimensional (3D) context.

Protein function prediction and annotation in an integrated environment powered by web service.

Use to predict function from de novo protein sequences.

Search for short active protein sequences with demonstrated biological activities.

Search for ungapped segments corresponding to the most highly conserved regions of proteins.

Identify and measure surface accessible pockets as well as interior inaccessible cavities, for proteins and other molecules.

To search for catalytic residue annotation for enzymes in the Protein Data Bank.

Predict protein function using Gene Ontology.

Automatically calculate evolutionary conservation scores of key amino acid residues and map them on protein structures.

Mine the protein structure space.

Predict short linear motifs (3-8 residues) in a set of protein sequences.

A web client for visualizing protein sequence feature information using DAS.

Identify the domain architecture within a protein sequence.

Predict functional sites in eukaryotic proteins.

Use a collection of tools for protein analyses.

Predict potential protein post-translational modifications and find potential single amino acid substitutions in peptides.

Search for information related to the catalytic mechanisms of enzymes.

An integrated feature-based function prediction server for vertebrate proteomes.

Identify the closest matching PRINTS sequence motif fingerprints in a protein sequence.

Search for functional annotation of important sites in proteins with known structures.

Produce 3D conformations of small drug compounds.

A database presenting experiment-based results in human proteomics.

Conduct exhaustive intermediate profile searches of a set of homologous protein sequences.

Design protein mutations in site-directed mutagenesis.

Annotate protein using this integrated annotation resource.

Identify protein family (and DNA) domains, patterns, motifs, protein families, and functional sites.

Interactive forecasting of protein interaction hot spots.

Database containing pairs of structural analogs and their alignments.

Find sequence patterns in DNA and protein sequences.

A web server for knowledge-based modeling of protein-peptide complexes, specifically peptides in complex with major histocompatibility complex (MHC) proteins and kinases.

Find structural segments or motifs for protein structures.

Visualize protein sequence motifs on the 3D protein structures.

Find presence of any known protein motif (Prosite and Pfam) in a protein sequence.

Recognize spatial chemical binding patterns common to a set of protein structures.

Analyze proteins for the presence of N-terminal N-myristoylation site.

Find the presence of N-Glycosylation sites in human proteins.

Find the presence of O-GalNAc (mucin type) glycosylation sites in mammalian proteins.

Analyze eukaryotic proteins for the presence of serine, threonine and tyrosine phosphorylation sites.

Find possible kinase specific phosphorylation sites in eukaryotic proteins.

Predict cleavage sites at basic amino acid locations in neuropeptide precursor sequences.

Find information about patented nucleotide and protein sequences.

Search for information about glycoproteins with O-linked and C-linked glycosylation sites.

Find information about protein sequence annotations.

A server to predict protein active site residues.

Search for structural and functional information on the protein functional sites.

Search 3D protein fragments similar in structure to known active, binding and posttranslational modification sites.

Conduct genome wide functional and structural analysis.

Search for phosphorylation data of any protein of interest.

Search for information on prokaryotic proteins that undergo serine, threonine, or tyrosine phosphorylation.

Determine correct names for proteins.

Web application for predicting protein disorder by using physicochemical features and reduced amino acid set of a position-specific scoring matrix.

Find homologous protein-protein interactions across multiple species.

Search your query sequence against PROSITE pattern database for protein motifs.

Find information about protein-RNA complexes from the Protein Data Bank (PDB).

Identify protein families and domains for a given protein sequence.

A comprehensive database of pattern-recognition receptors and their ligands.

Search for short nucleotide or peptide sequences such as cis-elements in nucleotide sequences or small domains and motifs in protein sequences.

Database specialized in documenting human PPBD-containing proteins and PPBD-mediated interactions.

Predicts potential protease cleavage sites and sites cleaved by chemicals in a given protein sequence.

Predict combined transmembrane topology and signal peptides.

Search for eukaryotic phosphorylation sites.

Search for 3D structure and functional annotation of phosphorylation sites in proteins.

Search the database of in vivo phosphorylation sites of human and mouse proteins

Find information about polyglutamine (polyQ) repeats.

Find the presence of protein motifs and patterns in an amino acid sequence.

Predict signal peptide sequences and their cleavage positions in bacterial and eukaryotic amino acid sequences.

Predict protein functions based on known structures.

Predict the location of potential protein-protein binding sites for unbound proteins.

Identify short linear signatures in protein termini.

Analyze and identify newly obtained protein sequences.

Predict protein binding sites in a protein sequence based on geometrical analysis of protein tertiary substructures.

Search for evolutionarily conserved motif-like patterns in protein sequences.

Web-based server for analyzing and predicting RNA binding sites in proteins.

Search for similarities between proteins by simultaneous matching of multiple motifs.

Predict residues in protein sequences that determine the proteins' functional specificity.

Predict specificity-determining residues in protein families.

Find shared motifs in proteins with a common attribute.

Conduct in silico sumoylation sites prediction.

Search for motifs and patterns within protein sequences.

Search for motifs within proteins that are likely to be phosphorylated by specific protein kinases or bind to domains such as SH2 domains, 14-3-3 domains or PDZ domains.

Find information about populations carrying polymorphisms within protein binding pockets that make them susceptible to serious adverse drug reaction (SADR).

Search the presence of a motif in either amino acid sequence or nucleotide sequence.

Predict the presence of tyrosine sulfation sites in protein sequences

Look at protein structure from a ligand and binding site perspective.

Use a collection of bioinformatics tools at this portal site.

Find information about tandem repeats in proteins that carry fundamental biological functions and are related to a number of human diseases.

Find information about functional residues in alpha-helical and beta-barrel membrane proteins.

Database of domains and motifs with conservative location in transmembrane proteins.

Search for highly conserved and specific protein sequence motifs.

Predict functional sites in protein sequence alignments use different methodologies.

Find de novo protein motifs from chromatin immunoprecipitation data.

Scan query structures for functional sites in both proteins and nucleic acids.

Analyze quantitative structure-activity relationship of related protein families.

Search for ungapped alignments of highly conserved regions among a protein family or superfamily.

An expert system for predicting ligand-binding residues in protein structures.

Automatically identify spatially interacting motifs among distantly related proteins sharing similar folds and possessing common ancestral lineage.

]]> Neel https://bioinformaticsonline.com/news/view/42023/encode3-a-collection-of-research-articles-and-related-content-describing-the-encyclopedia-of-dna-elements-its-datasets-and-tools Sat, 08 Aug 2020 08:25:21 -0500 https://bioinformaticsonline.com/news/view/42023/encode3-a-collection-of-research-articles-and-related-content-describing-the-encyclopedia-of-dna-elements-its-datasets-and-tools <![CDATA[ENCODE3: A collection of research articles and related content describing the Encyclopedia of DNA Elements, its datasets and tools.]]> How cells, tissues and organisms interpret the information encoded in the genome has vital implications for our understanding of development, health and disease. Launched in 2003, the ENCyclopedia Of DNA Elements (ENCODE) project has the aim of mapping the functional elements in the human genome (later expanded to include model organisms).

During the first phase of ENCODE, published in 2007, microarray-based technologies were used to detect regions associated with transcription factors, certain histone modifications and open chromatin within a pre-specified 1% of the human genome.

ENCODE’s second phase saw a switch to sequencing-based technologies, the addition of new assay types and the analysis of functional elements genome-wide, described in a collection of research articles in 2012.

The Encyclopedia paper of ENCODE 3, published in Nature, gives an overview of the various assays that were performed in human and mouse cell lines and tissues and describes a Registry of human and mouse candidate cis-regulatory elements (cCREs).

More at https://www.nature.com/immersive/d42859-020-00027-2/index.html

]]> Shruti Paniwala https://bioinformaticsonline.com/blog/view/42936/ancient-whole-genome-duplication-wgd-detection-tools Sun, 07 Mar 2021 00:32:44 -0600 https://bioinformaticsonline.com/blog/view/42936/ancient-whole-genome-duplication-wgd-detection-tools <![CDATA[Ancient whole genome duplication (WGD) detection tools !]]> There are two methods for ancient WGD detection, one is collinearity analysis, and the other is based on the Ks distribution map. Among them, Ks is defined as the average number of synonymous substitutions at each synonymous site, and there is also a Ka corresponding to it, which refers to the average number of non-synonymous substitutions at each non-synonymous site.

At present, some people have posted articles about the analysis process of WGD. I searched for the keyword "wgd pipeline" and found the following:

GenoDup: https:// github.com/MaoYafei/GenoDup-Pipeline
https://peerj.com/articles/6303/
WGDdetector: https:// github.com/yongzhiyang2 012/WGDdetector
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2670-3
wgd: https:// github.com/arzwa/wgd
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2#Sec1
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
GeNoGAP https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
https://github.com/dfguan/purge_dups
https://www.biorxiv.org/content/10.1101/2020.01.24.917997v1

This article introduces the usage of wgd.

Wgd cannot be installed directly with bioconda at present, so it is a little troublesome to install, because it depends on a lot of software. wgd depends on the following software

BLAST
MCL
MUSCLE/MAFFT/PRANK
PAML
PhyML/FastTree
i-ADHoRe

But the good news is that most of the software it depends on can be installed with bioconda

conda create -n wgd python=3.5 blast mcl muscle mafft prank paml fasttree cmake libpng mpi=1.0=mpich
conda activate wgd

Here mpi=1.0=mpich is selected, because i-adhore depends on mpich. If openmpi is installed, an error will appear while loading shared libraries: libmpi_cxx.so.40: cannot open shared object file: No such file or directory

After that, the installation is much simpler

git clone https://github.com/arzwa/wgd.git
cd wgd
pip install .
pip install git+https://github.com/arzwa/wgd.git
For i-ADHoRe, you need to register at http:// bioinformatics.psb.ugent.be /webtools/i-adhore/licensing/Agree to the license to download i-ADHoRe-3.0

Since my miniconda3 installed ~/opt/, the installation path is so~/opt/miniconda3/envs/wgd/

tar -zxvf i-adhore-3.0.01.tar.gz
cd i-adhore-3.0.01
mkdir -p build && cd build
cmake .. -DCMAKE_INSTALL_PREFIX=~/opt/miniconda3/envs/wgd/
make -j 4
make insatall

Take the sugarcane genome Saccharum spontaneum L as an example. The genome is 8-ploid with 32 chromosomes (2n = 4x8 = 32)

Download the tutorial for CDS and GFF annotation files

mkdir -p wgd_tutorial && cd wgd_tutorial
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.cds.fasta.gz
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.gff3.gz
gunzip *.gz

First conda activate wgdstart our analysis environment, and then start the analysis

Step 1 : Use to wgd mclidentify homologous genes in the genome

wgd mcl -n 20 --cds --mcl -s Sspon.v20190103.cds.fasta -o Sspon_cds.out

Step 2 : Use to wgd ksdbuild Ks distribution

wgd ksd --n_threads 80 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl Sspon.v20190103.cds.fasta

Step 3 : If the quality of the genome is good, then wgd syncollinearity analysis can be used . It can help us find the collinearity block in the genome and the corresponding anchor point

wgd syn --feature gene --gene_attribute ID \
-ks wgd_ksd/Sspon.v20190103.cds.fasta.ks.tsv \
Sspon.v20190103.gff3 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl

 For more reading - There are 9 sub-modules in WGD

  • kde: KDE fitting to the Ks distribution
  • ksd: Ks distribution construction
  • mcl: BLASP comparison of All-vs-ALl + MCL classification analysis.
  • mix: Hybrid modeling of Ks distribution.
  • pre: preprocess the CDS file
  • syn: Call I-ADHoRe 3.0 to use GFF files for collinearity analysis
  • viz: draw histogram and density plot
  • wf1: Ks standard analysis procedure of the whole genome paranome (paranome), call mcl, ksd and syn
  • wf2: Ks standard analysis procedure of one-vs-one homologous gene (ortholog), call wcl and kSD
]]> Rahul Nayak https://bioinformaticsonline.com/pages/view/43728/short-read-assembly-using-spades Mon, 31 Jan 2022 07:18:16 -0600 https://bioinformaticsonline.com/pages/view/43728/short-read-assembly-using-spades <![CDATA[Short-read assembly using Spades !]]> If we only had Illumina reads, we could also assemble these using the tool Spades.

You can try this here, or try it later on your own data.

Get data

We will use the same Illumina data as we used above:

  • illumina_R1.fastq.gz: the Illumina forward reads
  • illumina_R2.fastq.gz: the Illumina reverse reads

Assemble

Run Spades:

spades.py -1 illumina_R1.fastq.gz -2 illumina_R2.fastq.gz --careful --cov-cutoff auto -o spades_assembly_all_illumina
  • -1 is input file of forward reads
  • -2 is input file of reverse reads
  • --careful minimizes mismatches and short indels
  • --cov-cutoff auto computes the coverage threshold (rather than the default setting, “off”)
  • -o is the output directory

Results

Move into the output directory and look at the contigs:

infoseq contigs.fasta
]]> Abhimanyu Singh https://bioinformaticsonline.com/pages/view/44377/mitochondrial-genome-assembly-tools Wed, 06 Sep 2023 00:37:18 -0500 https://bioinformaticsonline.com/pages/view/44377/mitochondrial-genome-assembly-tools <![CDATA[Mitochondrial genome assembly tools !]]> Mitochondrial genome assembly tools are specialized software and algorithms designed to accurately reconstruct the mitochondrial genome (mitogenome) from sequencing data, typically obtained through techniques like next-generation sequencing (NGS). The mitochondrial genome is relatively small compared to the nuclear genome, making it an ideal target for assembly. Here are some commonly used mitochondrial genome assembly tools:

MitoFinder: Mitofinder is a pipeline to assemble mitochondrial genomes and annotate mitochondrial genes from trimmed read sequencing data.

MitoHiFi: a python pipeline for mitochondrial genome assembly from PacBio high fidelity reads

MITObim: MITObim is a tool specifically developed for the iterative assembly of mitochondrial genomes. It starts with a reference mitogenome and iteratively refines the assembly using the read data.

MITOS: MITOS is a web-based platform that provides a pipeline for annotating mitochondrial genomes. It integrates multiple software tools for assembly, annotation, and visualization of mitogenomes.

MIRA: MIRA (Mimicking Intelligent Read Assembly) is a versatile genome assembly tool that can be used for mitochondrial genome assembly. It supports various sequencing technologies and allows for reference-based or de novo assembly.

NOVOPlasty: NOVOPlasty is a user-friendly tool designed for de novo assembly of organelle genomes, including mitochondria. It utilizes a seed-and-extend algorithm and is suitable for both short-read and long-read data.

MITOS2: MITOS2 is an updated version of the MITOS pipeline, which automates the annotation of mitochondrial genomes. It provides improved accuracy and additional features for mitochondrial genome analysis.

GetOrganelle: While primarily designed for chloroplast genome assembly, GetOrganelle can also be used for mitochondrial genome assembly. It is particularly useful for dealing with high-throughput sequencing data.

SPAdes: SPAdes (St. Petersburg genome assembler) is a versatile genome assembly tool that can be employed for mitochondrial genome assembly, especially when dealing with complex datasets that may contain nuclear mitochondrial DNA sequences (numts).

IDBA-UD: IDBA-UD (Iterative De Bruijn Graph De Novo Assembler) is another de novo assembly tool that can be used for mitochondrial genome assembly, especially in cases with relatively low coverage.

Velvet: Velvet is a de novo assembly tool that can be applied to mitochondrial genome assembly, especially when working with short-read data.

When selecting a mitochondrial genome assembly tool, it's important to consider the specific characteristics of your sequencing data, such as read length and coverage, as well as the complexity of the mitochondrial genome. Additionally, some tools are better suited for specific organisms or research objectives, so choosing the right tool will depend on your particular project requirements.

]]> Abhi https://bioinformaticsonline.com/bookmarks/view/31278/metapred2cs Fri, 03 Mar 2017 05:15:07 -0600 https://bioinformaticsonline.com/bookmarks/view/31278/metapred2cs <![CDATA[MetaPred2CS]]> MetaPred2CS Web server is a meta-predictor based on Support Vector Machine (SVM) that combines 6 individual sequence based protein-protein interaction prediction methods to predict prokaryotic two-component system protein-protein interactions (PPIs). The methods implemented in MetaPred2CS are 2 co-evolutionary methods: in-silico two hybrid (i2h) and mirror tree (MT) methods and 4 genomics context based methods: phylogenetic profiling (PP), gene fusion (GF), gene neighbourhood (GN) and and gene operon methods (GO).

 http://metapred2cs.ibers.aber.ac.uk/

Address of the bookmark: https://github.com/martinjvickers/MetaPred2CS

]]> Manisha Mishra https://bioinformaticsonline.com/bookmarks/view/27847/anvio Thu, 16 Jun 2016 18:15:41 -0500 https://bioinformaticsonline.com/bookmarks/view/27847/anvio <![CDATA[Anvio]]> In a nutshell

Anvi’o is an analysis and visualization platform for ‘omics data.

Please find the methods paper here: https://peerj.com/articles/1319/

Anvi’o would not have been possible without the help of many people who directly or indirectly contributed to its development. Here is the acknowledgements section of our methods paper

An analysis and visualization platform for 'omics data http://merenlab.org/projects/anvio

Paper https://peerj.com/articles/1839/

Address of the bookmark: https://github.com/meren/anvio

]]> Shruti Paniwala