BOL: Related items

MyPro: A seamless pipeline for automated prokaryotic genome assembly and annotation

Neel — Thu, 15 Dec 2016 05:47:35 -0600

MyPro is an improved genomics software pipeline for prokaryotic genomes. MyPro is user-friendly and requires minimal programming skills. High-quality prokaryotic genome assembly and annotation can be obtained with ease. It performed better than de novo assemblers and contig integration software. Produces more contiguous assemblies, higher N50 values and lower number of contigs.

More at https://sourceforge.net/projects/sb2nhri/files/MyPro/

Address of the bookmark: http://www.sciencedirect.com/science/article/pii/S0167701215001207

Gene Synteny Database

Jit — Fri, 16 Dec 2016 11:09:39 -0600

Comparative genomics remains a pivotal strategy to study the evolution of gene organization, and this primacy is reinforced by the growing number of full genome sequences available in public repositories. Despite this growth, bioinformatic tools available to visualize and compare genomes and to infer evolutionary events remain restricted to two or three genomes at a time, thus limiting the breadth and the nature of the question that can be investigated. Here we present Genomicus, a new synteny browser that can represent and compare unlimited numbers of genomes in a broad phylogenetic view. In addition, Genomicus includes reconstructed ancestral gene organization, thus greatly facilitating the interpretation of the data.

Availability: Genomicus is freely available for online use at http://www.dyogen.ens.fr/genomicus while data can be downloaded at ftp://ftp.biologie.ens.fr/pub/dyogen/genomicus

Contact: rf.sne.eigoloib@crh

Address of the bookmark: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2853686/

YAHA

Jit — Fri, 20 Jan 2017 05:38:05 -0600

YAHA, a fast and flexible hash-based aligner. YAHA is as fast and accurate as BWA-SW at finding the single best alignment per query and is dramatically faster and more sensitive than both SSAHA2 and MegaBLAST at finding all possible alignments. Unlike other aligners that report all, or one, alignment per query, or that use simple heuristics to select alignments, YAHA uses a directed acyclic graph to find the optimal set of alignments that cover a query using a biologically relevant breakpoint penalty. YAHA can also report multiple mappings per defined segment of the query. We show that YAHA detects more breakpoints in less time than BWA-SW across all SV classes, and especially excels at complex SVs comprising multiple breakpoints.

Availability: YAHA is currently supported on 64-bit Linux systems. Binaries and sample data are freely available for download from http://faculty.virginia.edu/irahall/YAHA.

Contact:

http://genome.wustl.edu/people/groups/detail/hall-lab/

Address of the bookmark: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3463118/

MEME suite

Bulbul — Fri, 23 Dec 2016 08:49:55 -0600

Motif based sequence analysis suits

The MEME Suite allows the biologist to discover novel motifs in collections of unaligned nucleotide or protein sequences, and to perform a wide variety of other motif-based analyses.

The MEME Suite supports motif-based analysis of DNA, RNA and protein sequences. It provides motif discovery algorithms using both probabilistic (MEME) and discrete models (MEME), which have complementary strengths. It also allows discovery of motifs with arbitrary insertions and deletions (GLAM2). In addition to motif discovery, the MEME Suite provides tools for scanning sequences for matches to motifs (FIMO, MAST and GLAM2Scan), scanning for clusters of motifs (MCAST), comparing motifs to known motifs (Tomtom), finding preferred spacings between motifs (SpaMo), predicting the biological roles of motifs (GOMo), measuring the positional enrichment of sequences for known motifs (CentriMo), and analyzing ChIP-seq and other large datasets (MEME-ChIP).

The MEME Suite is comprised of a collection of tools that work together, as shown below. Not all the tools are available as webservices, so to get the full power of the MEME Suite you will need to download and install a local copy of the software. To see what has changed recently you can peruse the release notes.

http://meme-suite.org/

Address of the bookmark: http://meme-suite.org/

DnaSP v5: a software for comprehensive analysis of DNA polymorphism data

Jit — Mon, 06 Feb 2017 04:45:37 -0600

DnaSP is a software package for a comprehensive analysis of DNA polymorphism data. Version 5 implements a number of new features and analytical methods allowing extensive DNA polymorphism analyses on large datasets. Among other features, the newly implemented methods allow for: (i) analyses on multiple data files; (ii) haplotype phasing; (iii) analyses on insertion/deletion polymorphism data; (iv) visualizing sliding window results integrated with available genome annotations in the UCSC browser.

Address of the bookmark: http://www.ub.edu/dnasp/

Enrichr: a comprehensive gene set enrichment analysis

Jit — Thu, 27 Apr 2017 05:42:09 -0500

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

https://academic.oup.com/nar/article-lookup/doi/10.1093/nar/gkw377

Address of the bookmark: http://amp.pharm.mssm.edu/Enrichr/

RNA-seq Analysis Workshop Course Materials

Jit — Tue, 03 Jul 2018 08:14:14 -0500

RNAseq can be roughly divided into two "types": Reference genome-based - an assembled genome exists for a species for which an RNAseq experiment is performed. It allows reads to be aligned against the reference genome and significantly improves our ability to reconstruct transcripts. This category would obviously include humans and most model organisms but excludes the majority of truly biologically intereting species (e.g., Hyacinth macaw); Reference genome-free - no genome assembly for the species of interest is available. In this case one would need to assemble the reads into transcripts using de novo approaches. This type of RNAseq is as much of an art as well as science because assembly is heavily parameter-dependent and difficult to do well. In this lesson we will focus on the Reference genome-based type of RNA seq. http://chagall.med.cornell.edu/RNASEQcourse/

Address of the bookmark: http://chagall.med.cornell.edu/RNASEQcourse/

RNA-Seq Analysis: A Guide for Bioinformaticians

LEGE — Sat, 07 Dec 2024 22:22:24 -0600

RNA sequencing (RNA-Seq) has revolutionized transcriptomics, offering unprecedented insights into gene expression, splicing, and transcript diversity. For bioinformaticians, RNA-Seq analysis is a gateway to exploring the complexity of RNA biology and its implications in health and disease. This blog post provides an overview of RNA-Seq analysis, key computational steps, and tools for bioinformaticians eager to delve into this powerful technique.

What is RNA-Seq?

RNA-Seq is a next-generation sequencing (NGS) technology used to study the transcriptome—the complete set of RNA molecules in a cell. It quantifies gene expression, detects novel transcripts, and captures alternative splicing events with high sensitivity and resolution.

Workflow for RNA-Seq Analysis

RNA-Seq analysis involves several stages, each requiring computational tools and expertise.

1. Experimental Design and Data Acquisition

Before diving into analysis, bioinformaticians should consider:

Biological Replicates: Ensure statistical power to detect meaningful differences.
Sequencing Depth: Align sequencing depth to study objectives (e.g., higher depth for low-abundance transcripts).
Paired-End vs. Single-End: Paired-end sequencing provides more detailed information on transcript structure.

Once sequencing is complete, raw data is provided in FASTQ format, containing sequence reads and quality scores.

2. Quality Control and Preprocessing

Quality control (QC) ensures data integrity. Tools such as FastQC evaluate metrics like base quality, GC content, and adapter contamination.

Preprocessing Steps:

Trimming: Tools like Trimmomatic or Cutadapt remove low-quality bases and adapter sequences.
Filtering: Discard reads below a certain quality threshold or length.

3. Read Alignment

Reads are mapped to a reference genome or transcriptome to determine their origin. Alignment tools include:

HISAT2: Handles large genomes efficiently and supports spliced alignments.
STAR: High-speed aligner optimized for RNA-Seq.
Bowtie2: Suitable for short-read alignment.

Output: A SAM/BAM file containing aligned reads.

4. Transcript Assembly and Quantification

This step involves identifying transcripts and quantifying their expression levels. Tools used include:

StringTie: Assembles and quantifies transcripts from aligned reads.
Salmon/Kallisto: Perform pseudo-alignment for rapid and accurate quantification.

Expression levels are typically measured as TPM (transcripts per million) or FPKM (fragments per kilobase of transcript per million mapped reads).

5. Differential Expression Analysis

To identify genes with altered expression between conditions, bioinformaticians use tools such as:

DESeq2: Accounts for data normalization and variability.
edgeR: Handles overdispersed count data efficiently.
Limma-voom: Combines linear modeling with RNA-Seq count data.

The output includes a list of differentially expressed genes (DEGs) with statistical significance and fold-change values.

6. Functional Annotation and Pathway Analysis

Understanding the biological significance of DEGs involves:

Gene Ontology (GO) Analysis: Tools like DAVID or clusterProfiler categorize genes based on their biological functions.
Pathway Enrichment Analysis: Identifies pathways enriched in DEGs using tools like KEGG, Reactome, or GSEA.

7. Visualization

Visualizing results enhances interpretability. Common visualizations include:

Heatmaps: Show expression patterns across samples (e.g., pheatmap).
Volcano Plots: Highlight significant DEGs (e.g., ggplot2).
PCA/UMAP: Assess sample clustering and variability (e.g., Seurat).

Challenges in RNA-Seq Analysis

Batch Effects: Technical variability can confound biological signals. Combat this with normalization techniques or batch-correction tools like ComBat.
Low-Quality Samples: Poor-quality RNA impacts downstream analyses.
Computational Complexity: RNA-Seq generates massive datasets, requiring robust computing resources and optimized pipelines.

Key Tools and Resources

Bioconductor: A treasure trove of R packages for RNA-Seq analysis.
Galaxy: A web-based platform for running RNA-Seq workflows.
Nextflow/Snakemake: Workflow management tools to streamline analyses.

Applications of RNA-Seq

RNA-Seq is used in diverse research areas, including:

Cancer Transcriptomics: Identifying tumor-specific expression profiles.
Developmental Biology: Studying dynamic transcriptome changes.
Drug Discovery: Screening genes modulated by therapeutic compounds.

Conclusion

RNA-Seq analysis is a cornerstone of modern transcriptomics, offering bioinformaticians a versatile toolkit for unraveling gene expression and regulation. Mastering RNA-Seq workflows and tools empowers researchers to transform raw sequencing data into biological discoveries.

Whether you’re investigating disease mechanisms, exploring cellular pathways, or developing new therapeutics, RNA-Seq is a powerful ally in your bioinformatics arsenal.

Bioinformatics tools developed for Oxford Nanopore data analysis !

biogeek — Wed, 27 Dec 2017 20:47:30 -0600

MinION is the only portable real-time device for DNA and RNA sequencing. Each consumable flow cell can now generate 10–20 Gb of DNA sequence data. Ultra-long read lengths are possible (hundreds of kb) as you can choose your fragment length. One of the technical advantages of ONT data is the read length, which offers great prospects for genome assembly. Generally, assemblers are based on several different types of algorithms, such as greedy, overlap-layout-consensus (OLC), de Bruijn graph (DBG), and string graph.

List of analysis tools developed for Oxford Nanopore data

BWA
Fast nanopore data tuned alignment tool
https://github.com/lh3/bwa

GraphMap
Mapper for long and error-prone reads
https://github.com/isovic/graphmap

LAST
Nanopore tuned alignment tool
http://last.cbrc.jp/

LINKS
Software tool for long read scaffolding
https://github.com/warrenlr/LINKS/

marginAlign
Tools to align nanopore reads to a reference
https://github.com/benedictpaten/marginAlign

minoTour
Real time analysis tools
http://minotour.nottingham.ac.uk/

nanoCORR
Error-correction tool for nanopore sequence data
https://github.com/jgurtowski/nanocorr

NanoOK
Software for nanopore data, quality and error profiles
https://documentation.tgac.ac.uk/display/NANOOK/NanoOK

Nanopolish
Nanopore analysis and genome assembly software
https://github.com/jts/nanopolish

nanopore
Variant-detection tool for nanopore sequence data
https://github.com/mitenjain/nanopore

Nanocorrect
Error-correction tool for nanopore sequence data
https://github.com/jts/nanocorrect/

npReader
Real-time conversion and analysis of nanopore reads
https://github.com/mdcao/npReader

poRe
Tool for analyzing and visualizing nanopore data
https://sourceforge.net/p/rpore/wiki/Home/

PoreSeq
Error-correction and variant-calling software
https://github.com/tszalay/poreseq

Poretools
Nanopore sequence analysis and visualization software
https://github.com/arq5x/poretools

SSPACE-LongRead
Genome scaffolding tool
http://www.baseclear.com/genomics/bioinformatics/basetools/SSPACE-longread

SMIS
Genome scaffolding tool
https://sourceforge.net/projects/phusion2/files/smis/

List of assemblers for Oxford Nanopore MinION long reads

LQS
DALIGNER, Celera OLC Nanocorrect,
Nanopolish corrector
https://github.com/jts/nanopolish

PBcR
HGAP or BLASR, Celera OLC
PBcR corrector
http://wgs-assembler.sourceforge.net/wiki/index.php/PBcR
–
Canu
MHAP, Celera OLC
Canu corrector
https://github.com/marbl/canu

Falcon
String graph, Celera OLC
Falcon corrector
https://github.com/PacificBiosciences/falcon

Miniasm
OLC
https://github.com/lh3/miniasm

ra-integrate
OLC
https://github.com/mariokostelac/ra-integrate/

ALLPATHS-LG
de Bruijn graph
ALLPATHS-L corrector
https://www.broadinstitute.org/software/allpaths-lg/blog/?page_id=12

SPAdes
de Bruijn graph
SPAdes corrector
http://bioinf.spbau.ru/spades

SciLifeLab tutorial for bioinformatics analysis !

Jit — Tue, 17 Apr 2018 04:33:00 -0500

SciLifeLab is a national center for molecular biosciences with focus on health and environmental research.

Courses

Old courses (2012-2014)