<![CDATA[BOL: Related items]]> https://bioinformaticsonline.com/related/29638?offset=70 https://bioinformaticsonline.com/bookmarks/view/29586/eforgev12 Fri, 28 Oct 2016 09:06:59 -0500 https://bioinformaticsonline.com/bookmarks/view/29586/eforgev12 <![CDATA[eFORGE.v1.2]]> The eFORGE tool provides a method to view the tissue specific regulatory component of a set of EWAS DMPs. eFORGE analysis takes a set of DMPs, such as those hits above genome-wide significance threshold in an EWAS study, and analyses whether there is enrichment for overlap of putative functional elements compared to matched background DMPs. It assesses enrichment on a per cell type basis, since functional elements are differentially active in different cell types, and hence can expose tissue-specific signals of enrichment for the given test DMP set. This can reveal the sites of action underlying the EWAS signal, and provide confirmation of the validity of the EWAS where a tissue-specific mechanism is known or expected for the phenotype. Conversely unknown tissue involvements can also be revealed.

Address of the bookmark: http://eforge.cs.ucl.ac.uk/eFORGE.v1.2/?documentation

]]> Jit https://bioinformaticsonline.com/bookmarks/view/30076/sga-string-graph-assembler Thu, 08 Dec 2016 05:08:59 -0600 https://bioinformaticsonline.com/bookmarks/view/30076/sga-string-graph-assembler <![CDATA[SGA: String Graph Assembler]]> SGA is a de novo genome assembler based on the concept of string graphs. The major goal of SGA is to be very memory efficient, which is achieved by using a compressed representation of DNA sequence reads.

More at

https://github.com/jts/sga

SGA dependencies:
-google sparse hash library (http://code.google.com/p/google-sparsehash/)
-the bamtools library (https://github.com/pezmaster31/bamtools)
-zlib (http://www.zlib.net/)
-(optional but suggested) the jemalloc memory allocator (http://www.canonware.com/jemalloc/download.html)

Address of the bookmark: https://github.com/jts/sga

]]> Jit https://bioinformaticsonline.com/bookmarks/view/30140/cutadapt Wed, 14 Dec 2016 09:59:52 -0600 https://bioinformaticsonline.com/bookmarks/view/30140/cutadapt <![CDATA[Cutadapt]]> Cutadapt finds and removes adapter sequences, primers, poly-A tails and other types of unwanted sequence from your high-throughput sequencing reads.

Cutadapt helps with these trimming tasks by finding the adapter or primer sequences in an error-tolerant way. It can also modify and filter reads in various ways. Adapter sequences can contain IUPAC wildcard characters. Also, paired-end reads and even colorspace data is supported. If you want, you can also just demultiplex your input data, without removing adapter sequences at all.

Cutadapt comes with an extensive suite of automated tests and is available under the terms of the MIT license.

If you use cutadapt, please cite DOI:10.14806/ej.17.1.200 .

More at https://github.com/marcelm/cutadapt

Address of the bookmark: http://cutadapt.readthedocs.io/en/stable/guide.html

]]> Bulbul https://bioinformaticsonline.com/bookmarks/view/30214/megamerge-a-tool-to-merge-assembled-contigs-long-reads-from-metagenomic-sequencing-runs Mon, 19 Dec 2016 09:42:15 -0600 https://bioinformaticsonline.com/bookmarks/view/30214/megamerge-a-tool-to-merge-assembled-contigs-long-reads-from-metagenomic-sequencing-runs <![CDATA[MeGAMerge: A tool to merge assembled contigs, long reads from metagenomic sequencing runs]]> MeGAMerge

MeGAMerge (A tool to merge assembled contigs, long reads from metagenomic sequencing runs)

Description

MeGAMerge is a perl based wrapper/tool that can accept any number of sequence (FASTA) files containing assembled contigs of any length in Multi-FASTA format to produce an improved contig set based on OLC based assembly. All overlap parameters (Minimum Overlap Length, Identity, etc) are user-declarable at runtime. It is written to run on Linux.

Requirements:

You will need to have the following tools installed and in $PATH, or added to $binpath in the tool:

Newbler (specifically runAssembly)
Minimus2 (part of AMOS, also requires MUMmer)

Address of the bookmark: https://github.com/LANL-Bioinformatics/MeGAMerge

]]> Jit https://bioinformaticsonline.com/bookmarks/view/36880/jvarkit-java-utilities-for-bioinformatics Fri, 08 Jun 2018 09:31:55 -0500 https://bioinformaticsonline.com/bookmarks/view/36880/jvarkit-java-utilities-for-bioinformatics <![CDATA[Jvarkit : Java utilities for Bioinformatics]]> Address of the bookmark: http://lindenb.github.io/jvarkit/

]]> Jit https://bioinformaticsonline.com/pages/view/11399/next-generation-sequencing-in-r-or-bioconductor-environment Mon, 02 Jun 2014 18:03:09 -0500 https://bioinformaticsonline.com/pages/view/11399/next-generation-sequencing-in-r-or-bioconductor-environment <![CDATA[Next generation sequencing in R or bioconductor environment]]> There are many R software and bioconductor packages for NGS data analysis, some of them are as follows

Biostrings

The Biostrings package from Bioconductor provides an advanced environment for efficient sequence management and analysis in R. It contains many speed and memory effective string containers, string matching algorithms, and other utilities, for fast manipulation of large sets of biological sequences. The objects and functions provided by Biostrings form the basis for many other sequence analysis packages. Documentation

IRanges Overview

IRanges provides the low-level infrastructure and containers for handling sets of integer ranges within Bioconductor's BioC-Seq domain. Its classes and methods provide support for many more high-level packages like GenomicRanges, ShortRead, Rsamtools, etc. Documentation

GenomicRanges Overview

The GenomicRanges package serves as the foundation for representing genomic locations within the Bioconductor project. It is built upon the IRanges infrastructure and defines three major data containers - GRanges, GRangesList and GappedAlignments - which are supporting other important BioC-Seq packages including ShortRead, Rsamtools, rtracklayer, GenomicFeatures and BSgenome.  Compared to the IRanges container, the GRanges/GRangesList classes are more flexible and extensible to store additional information about sequence ranges, such as chromosome identifiers (sequence space), strand information and annotation data. Documentation

Motif Discovery

cosmo

The cosmo package allows to search a set of unaligned DNA sequences for a shared motif that may function as transcription factor binding site. The algorithm extends the popular motif discovery tool MEME (Bailey and Elkan, 1995) in that it allows the search to be supervised by specifying a set of constraints that the motif to be discovered must satisfy. Documentation

BCRANK

BCRANK is a method that takes a ranked list of genomic regions as input and outputs short DNA sequences that are overrepresented in some part of the list. The algorithm was developed for detecting transcription factor (TF) binding sites in a large number of enriched regions from high-throughput ChIP-chip or ChIP-seq experiments, but it can be applied to any ranked list of DNA sequences. Documentation

rGADEM: Documentation

MotIV: Documentation

ShortRead

The ShortRead package provides input, quality control, filtering, parsing, and manipulation functionality for short read sequences produced by high throughput sequencing technologies. While support is provided for many sequencing technologies, this package is primairly focused on Solexa/Illumina reads. Documentation

Rsamtools

Rsamtools provides functions for parsing and inspecting samtools BAM formatted binary alignment data. SAM/BAM is quickly becoming a universal standard alignment format, and is now supported by a wide variety of alignment tools. Documentation

Samtools Website
BWA (Burrows-Wheeler Alignment) Website

 

Additional tools for SNP analysis: 

snpMatrix

BSgenome

BSgenome provides an object oriented infrastructure for interacting with a Biostring based genome sequence. BSgenome packages exist for many common genomes, and can be created to represent custom genomes. See the "How to forge a BSgenome data package" Vignette for instructions to create a new BSgenome package if a prebuilt package does not exist for your organism. Documentation

rtracklayer

rtracklayer provides an interface for exporting annotation feature data to various genome browsers and file formats (such as GFF). See the Small RNA Profiling exercise for an example of using rtracklayer to visualize alignment coverage. Documentation

biomaRt

The biomaRt package, provides an interface to a growing collection of databases implementing the BioMart software suite (http:// www.biomart.org). The package enables online retrieval of large amounts of data in a uniform way without the need to know the underlying database schemas. This data is retrieved automatically via the Internet, so it's recommended that you cache the data locally, or check versions if your code will be adversely affected by updates to these data. Documentation

ChIP-Seq Analysis Packages

Bioconductor provides various packages for analyzing and visualizing ChIP-Seq data. Only a small selection of these packages is introduced here. Additional useful introductions to this topic are: BioC ChIP-seq Case Study and BioC ChIP-Seq.

chipseq

The chipseq package combines a variety of HT-Seq packages to a pipeline for ChIP-Seq data analysis. Documentation

BayesPeak

BayesPeak is a peak calling package for identifying DNA binding sites of proteins in ChIP-Seq experiments. Its algorithm uses hidden Markov models (HMM) and Bayesian statistical methods. The following sample code introduces the identification of peaks with the BayesPeak package as well as the incorporation of read coverage information obtained by the chipseq package. Documentation [ Publication ]

PICS

The PICS package applies probabilistic inference to aligned-read ChIP-Seq data in order to identify regions bound by transcription factors. PICS identifies enriched regions by modeling local concentrations of directional reads, and uses DNA fragment length prior information to discriminate closely adjacent binding events via a Bayesian hierarchical t-mixture model. The following sample code uses the test data set from the above BayesPeak package in order to compare the results from both methods by identifying their consensus peak set. Documentation [ Publication ]

ChIPpeakAnno

The ChIPpeakAnno package provides. batch annotation of the peaks identified from either ChIP-seq or ChIP-chip experiments. It includes functions to retrieve the sequences around peaks, obtain enriched Gene Ontology (GO) terms, find the nearest gene, exon, miRNA or custom features such as most conserved elements and other transcription factor binding sites supplied by users. The package leverages the biomaRt, IRanges, Biostrings, BSgenome, GO.db, multtest and stat packages. Documentation

Additional ChIP-Seq Packages

DiffBind: Documentation

MOSAICS: Documentation

iSeq: Documentation

ChIPseqR: Documentation

ChiPsim: Documentation

CSAR: Documentation

ChIP-Seq Pipeline: PICS, rGADEM and MotIV (developer web site)

SPP: ChIP-seq processing pipeline

SPP Tutorial

MACS

SIPeS

RNA-Seq Analysis

Counting Reads that Overlap with Annotation Ranges 

The GenomicRanges package provides support for importing into R short read alignment data in BAM format (via Rsamtools) and associating them with genomic feature ranges, such as exons or genes. This way one can quantify the number of reads aligning to annotated genomic regions. The package defines general purpose containers for storing genomic intervals as well as more specialized containers for storing alignments against a reference genome. The two main functions for read counting provided by this infrastructure are countOverlaps and summarizeOverlaps. For their proper usage, it is important to read the corresponding PDF manual. Documentation

Differential Gene Expression Analysis with DESeq

The DESeq package contains functions to call differentially expressed genes (DEGs) in count tables based on a model using the negative binomial distribution. It expects as input a data frame with the raw read counts per region/gene of interest (rows) for each test sample (columns).  Such a count table can be imported into R or generated from BAM alignment files using the countOverlaps function as introduced above. Documentation

Differential Gene Expression Analysis with edgeR

The edgeR package uses empirical Bayes estimation and exact tests based on the negative binomial distribution to call differentially expressed genes (DEGs) in count data. 

Documentation

A variety of additional R packages are available for normalizing RNA-Seq read count data and identifying differentially expressed genes (DEG):

easyRNASeq (simplifies read counting per genome feature)

DEXSeq (Inference of differential exon usage); parathyroidSE explains how to generate exon read counts in R

DEGseq

baySeq (also see: segmentSeq)

Genominator (Bullard et al. 2010)

Detection of Alternative Splice Junctions

Another utility of RNA-Seq experiments is the analysis of splice junctions. The following software suggestions provide this utility:

ERANGE
TopHat

SpliceMap

SplitSeek

DNA-Methylation Data Analysis

HT-Seq Data Visualization

ggbio: ggplot2 extension for genomics data (online manual) Gviz: Plotting data and annotation information along genomic coordinates HilbertVis: Hilbert genome plots

GenomeGraphs: Plotting genomic information from Ensembl

TileQC: Flow Cell Quality Visualization

rtracklayer: R interface to genome browsers

genoPlotR: Plotting maps of genes and genomes

Genominator: Tools for storing, accessing, analyzing and visualizing genomic data.

 

To install all packages

source("http://bioconductor.org/biocLite.R")
biocLite()
biocLite(c("ShortRead", "Biostrings", "IRanges", "BSgenome", "rtracklayer", "biomaRt", "chipseq", "ChIPpeakAnno", "Rsamtools", "BayesPeak", "PICS", "GenomicRanges", "DESeq", "edgeR", "leeBamViews", "GenomicFeatures", "BSgenome.Celegans.UCSC.ce2"))

]]> John Parker https://bioinformaticsonline.com/pages/view/21365/a-guide-for-complete-r-beginners Fri, 20 Feb 2015 23:36:46 -0600 https://bioinformaticsonline.com/pages/view/21365/a-guide-for-complete-r-beginners <![CDATA[A guide for complete R beginners !]]> This tutorial is intended to introduce users quickly to the basics of R, focusing on a few common tasks that  biologists need to perform  some basic analysis:  load a table, plot some graphs, and perform some basic statistics. More extensive tutorials can be found on the project website and via bioconductor (not covered here).

R-language: http://www.r-project.org

BioConductor: http://www.bioconductor.org

Advantages of R

  • Free!
  • Powerful, many libraries have been created to perform application specific tasks. e.g. analysis of microarray experiments and Next-Gen sequencing (bioconductor: including Bioseq group).
  • Presentation quality graphics
    • Save as a png, pdf or svg
  • History
    • What you do can be saved for the next time you use R.
    • Ability to turn it into an automated script to perform again and again on different data

Disadvantages

Preparation

  • (Optional) Download and save the tutorial data set from
    • http://bioinformatics.knowledgeblog.org/wp-content/uploads/bioinf/kerr/data.tsv
    • Start R (type R on a Linux or Mac terminal, or find the starting link from PC)

Getting More Help

  • Project Home page
    • http://www.r-project.org/
    • Check out the ‘introduction to R’, which is a much more in depth guide .
    • Also R has a built-in help system (see later)

Working directory

This is the directory used to store your data and results. It is useful if it is also the directory where your input data is stored.

  • Mac/Linux: this is the directory where you typed in R
  • PC: Change using the change working directory option
]]> Archana Malhotra https://bioinformaticsonline.com/bookmarks/view/22961/bioscripts Sun, 28 Jun 2015 07:46:14 -0500 https://bioinformaticsonline.com/bookmarks/view/22961/bioscripts <![CDATA[BioScripts]]> You are requested to please bookmark collection of bioinformatics tools, scripts, codes that can be pieced together in a very easy and flexible manner to perform both simple and complex bioinformatics tasks.

The next-generation sequencing included whole genome sequencing(WGS), transcriptome sequencing (whole cDNA sequencing, RNA-seq), digital gene expression sequencing (Tag-Seq), ChIP-Seq, and so on. And there are many sequencing platform to generate sequece, as well know Sanger/ABi(the frist generation), Solexa/illumina, SOLiD/ABi, 454/Roche. But thier sequence format is different, also they have different error type. High quality data is very important for further analysis or data mining. There are many pipeline for raw sequence quality analysis and control with few of process for reporting reads quality statistical details, trimming, filtering, and error correction. Please bookmarks them for the benefits of bioinformatics community.

https://code.google.com/p/biowiki/

https://code.google.com/p/ngs-pipeline/source/browse/#svn%2Ftrunk

NGSand Perl scripts https://code.google.com/hosting/search?q=NGS+perl&projectsearch=Search+projects

NGS and Python scripts https://code.google.com/hosting/search?q=NGS+Python&projectsearch=Search+projects

Address of the bookmark: https://code.google.com/hosting/search?q=bioinformatics&sa=Search

]]> Rahul Nayak https://bioinformaticsonline.com/researchlabs/view/23209/bisr-jaipur Tue, 07 Jul 2015 23:12:26 -0500 <![CDATA[BISR Jaipur]]> The Bioinformatics Centre at BISR has created an infrastructure for providing facilities to the users working in the field of Biological Sciences. The users of Rajasthan, Jaipur in particular, are using facilities available at the Bioinformatics Centre extensively. The centre has leased line Internet connection as well latest Bioinformatics software for sequence and structure analysis. The centre provides the following services:

Bioinformatics supports to researchers
Customized training in Bioinformatics for researchers and faculty members
Support in Installing, implementing and maintaining software on computer.
Create awareness for taking preventive measure against data security
Organize workshops on thrust ares of Bioinformatics
Research Training to students of Biotechnology and Bioinformatics

More at http://bioinfo.bisr.res.in/index.php

]]> https://bioinformaticsonline.com/opportunity/view/25284/rajiv-gandhi-centre-for-biotechnology-rgcb-invites-applications-for-the-following-three-faculty-scientist Tue, 24 Nov 2015 22:13:16 -0600 <![CDATA[Rajiv Gandhi Centre for Biotechnology (RGCB) invites applications for the following three faculty scientist]]> Scientist Positions
Advt. No.RGCB Advt./SCI 2015/1

November 11, 2015

Rajiv Gandhi Centre for Biotechnology (RGCB) invites applications for the following three faculty scientist positions:

Scientist E-II or F in Bioinformatics & Computational Biology

SCIENTIST E-II OR F IN COMPUTATIONAL BIOLOGY & BIOINFORMATICS

Highly motivated and innovative individual who will pursue basic research, solve biological problems with emphasis on computational and quantitative experimental methods and build active bridges to translational research. The scientist will also provide computational biology support to ongoing research programs in disease biology, provide assistance to analyze complex data sets generated by RGCB scientists and collaborators inclusive of including high dimensional “omics” data and next generation sequencing data, such as whole genome, exome, RNA-seq and ChIP-seq as well as provide leadership for high quality training for junior scientists and regular teaching programs of the institute. Areas of research of interest to RGCB include but are not limited to computational, systems, or quantitative biology with applications to cell biology, developmental biology, metabolism, genomics, proteomics, biophysics, biological information systems, network pharmacology, drug design and cancer research. The scientist’s responsibilities include efforts for the integration of DNA variant annotation with statistical genetic analysis methods including linkage, imputation and association methods, adopting novel and innovative methodologies to analyze, integrate and interpret high dimensional data sets, provision of annotation to robust genetics and genomics findings using several data sources and methods, data management of exploratory clinical and R&D studies in partnership with other lines of genetic data generated from internal and external studies, delivery and documentation of genomic information to support genetic studies, ensuring high-quality genetic and genomic data is incorporated into exploratory- clinical research programs, developing tools that make maximum use of multiple data sources to support annotation of DNA variation and contributes to systems biology initiatives within RGCB

More at http://rgcb.res.in/scientist-positions/

Application Form http://rgcb.res.in/wp-content/uploads/2015/11/APPLICATION-FORMAT-FOR-SCIENTISTS.docx

]]>