Basic Galaxy Tutorial

• RNA-Seq tutorial based on Trapnell et al. (2012) Nature Protocols

In this tutorial we cover the concepts of RNA-Seq differential gene expression (DGE) analysis using a very small synthetic dataset from a well studied organism.

Advanced Galaxy Tutorial

• RNA-Seq (Advanced) Tutorial

In this tutorial we compare the performance of three statistically-based differential expression tools:

* CuffDiff

* EdgeR

* DESeq2

Advanced Command Line Tutorial

• Graphical Output with CummeRbund introduces some basic commands using the cummeRbund package of the R programming language

You will need to install R, RStudio and cummeRbund on your PC (explained in the Tutorial). You will learn how to produce graphical output from RNA-Seq analysis previously done using a Cuffdiff analysis.

Variant Detection

Basic Galaxy Tutorial

• Variant Detection tutorial

In this tutorial we cover the concepts of detecting small variants (SNVs and indels) in human genomic DNA using a small set of reads from chromosome 22.

Advanced Galaxy Tutorial

• Variant Detection (Advanced) Tutorial

In this tutorial we compare the performance of three statistically-based variant detection tools:

* SAMtools: Mpileup

* GATK: Unified Genotyper

* FreeBayes

Each of these tools takes as its input a BAM file of aligned reads and generates a list of likely variants in VCF format

Pipelines are for those who are comfortable with using the UNIX command line; and often allow more control over branching and iteration logic.

• WGS/exome GATK-based variant calling pipeline

This is a basic variant-calling and annotation pipeline developed at the Victorian Life Sciences Computation Initiative (VLSCI), University of Melbourne. It is based around BWA, GATK and ENSEMBL and was originally designed for human (or similar) data. The master branch is configured for WGS data; there is an exome branch configured for variant calling in exome data.

To run the pipeline you will need Rubra: https://github.com/bjpop/rubra. Rubra uses the python Ruffus library: http://www.ruffus.org.uk/.

Protocols

• Familial Variant Calling

In this protocol we discuss and outline the process of calling familial related mutations.

• Somatic Variant Calling

In this protocol we discuss and outline the process of identifying somatic variants or mutations.

Assembly

Basic Galaxy Tutorial

• Genome assembly tutorial

In this tutorial we carry out de novo assembly of a microbial genome. We have also written a De novo Genome Assembly for Illumina Data Protocol for a more generic description of the method.

Protocol

• De novo Genome Assembly for Illumina Data

In this protocol we discuss and outline the process of de novo assembly for small to medium sized genomes. Use our Genome assembly tutorial to learn a specific case of using Galaxy to carry out de novo assembly of a microbial genome.

Small RNAs

Basic Galaxy Tutorial

• Quality control for small RNA

This tutorial covers initial steps of the workflow for analysis of short RNA expression such as a quality control of the raw reads, processing of the raw reads for the subsequent analysis and initial quality assessment of the library.

ChIP Seq

Protocol

• ChIP-Seq

In this protocol we discuss ChIP-Seq: a method to analyze the interaction between proteins and DNA.

Amplicons

Protocol

• Amplicon Alignment

In this protocol we discuss and outline the process of aligning custom amplicons using primers for high precision.

Learn Galaxy

Introduction to Galaxy, for those who are very new to Galaxy.

Using Histories and Workflows, for those with some Galaxy knowledge.

The Galaxy project website has many tutorials and screencasts about using Galaxy and the tools, and developing new tools.

Address of the bookmark: https://genome.edu.au/wiki/Learn

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

• Lack of a comprehensive graphical user interface, but two do exist: However some do exist: R commander: http://socserv.mcmaster.ca/jfox/Misc/Rcmdr/ and Limma-gui (microarrays) : http://bioinf.wehi.edu.au/limmaGUI/

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

This is an extensive tutorial to take you through the main features and gotchas of configuring GBrowse as a server. This tutorial assumes that you have successfully set up Perl, GD, BioPerl and the other GBrowse dependencies. If you haven't, please see the GBrowse HOWTO During most of the tutorial, we will be using the "in-memory" GBrowse database (no relational database required!) Later we will show how to set up a genome size database using the berkeleydb and MySQL adaptors.

More at http://elp.ucdavis.edu/tutorial/tutorial.html

Address of the bookmark: http://elp.ucdavis.edu/tutorial/tutorial.html

1. Setting up your computing environment
2. Retrieving and storing sequencing files (your own data or from public sources)
3. Checking sequence quality, trimming, general sequence manipulation
4. Mapping reads to a reference genome
5. Manipulating SAM/BAM alignment files
6. Visualizing data in a genome browser

RNA-Seq

1. De novo transcript discovery and differential analysis with Cufflinks
2. Differential expression analysis with R/Bioconductor
3. Clustering of large expression datasets (microarray or RNA-Seq)

Microarray

1. Basic analysis of Affymetrix Gene Expression Arrays using R/Bioconductor

General Tips for Data Analysis

1. Excel workarounds, adding gene annotation, X-Y plots tips, etc.

Address of the bookmark: http://homer.salk.edu/homer/basicTutorial/

More at https://gatb.inria.fr/

Address of the bookmark: https://gatb.inria.fr/

Before you get started be sure to check out:

• devtools Google Group – https://groups.google.com/forum/#!forum/rdevtools
• book on “Advanced R programming” – http://adv-r.had.co.nz/
• GitHub repository – https://github.com/hadley/devtools

Address of the bookmark: https://www.rstudio.com/products/rpackages/devtools/

This page is a HOWTO guide for setting up a SLURM installation, currently focused on a CentOS 7 Linux OS. Please send feedback to Ole.H.Nielsen /at/ fysik.dtu.dk.

See the SLURM homepage (also https://computing.llnl.gov/linux/slurm/).

Address of the bookmark: https://wiki.fysik.dtu.dk/niflheim/SLURM

[uzi@quince-srv2 ~/]$ cp -r /home/opt/MScBioinformatics/linuxTutorial .
[uzi@quince-srv2 ~/]$ cd linuxTutorial
[uzi@quince-srv2 ~/linuxTutorial]$

I have deliberately chosen Awk in the exercises as it is a language in itself and is used more often to manipulate NGS data as compared to the other command line tools such as grep, sed, perl etc. Furthermore, having a command on awk will make it easier to understand advanced tutorials such as Illumina Amplicons Processing Workflow.

In Linux, we use a shell that is a program that takes your commands from the keyboard and gives them to the operating system. Most Linux systems utilize Bourne Again SHell (bash), but there are several additional shell programs on a typical Linux system such as ksh, tcsh, and zsh. To see which shell you are using, type

[uzi@quince-srv2 ~/linuxTutorial]$ echo $SHELL

/bin/bash

Address of the bookmark: http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/linux.html

Address of the bookmark: http://www.homolog.us/Tutorials/index.php

Full documentation on Read the Docs.

A set of tools for viral metagenomics.

virmet is called with a command subcommand syntax: virmet fetch --viral n, for example, downloads the bacterial database. Other available subcommands so far are

• fetch download genomes
• update update viral/bacterial database
• index index genomes
• wolfpack analyze a Miseq run
• covplot plot coverage for a specific organism

A short help is obtained with virmet subcommand -h.

More at https://github.com/ozagordi/VirMet

Address of the bookmark: https://github.com/ozagordi/VirMet