In this tutorial we will walk through the basics of generating an image. This walk through will be comparing an E. coli genome with five other E. coli genomes and mapping the read coverage from the underlying genome assembly onto the same image. For this walk through, users will need BRIG-Example.zip, which is available from the BRIG website (http://sourceforge.net/projects/brig/files/). This contains all the genomes and files needed to follow along with this walk through. Unzip it somewhere easily accessible, like the home directory or desktop.

About the reference genome

The reference genome used in this walk through is a simulated E. coli genome assembly. We took the published E. coli O157:H7 Sakai genome (Accession number BA000007) sequence and had assembly reads simulated by METASIM and then assembled these using Newbler version 2.3. The resulting contiguous sequences were ordered using Mauve against the published Sakai genome. This simulated E. coli is useful for illustrating some of BRIG’s graphing features for assembly read coverage.

Enterohemorrhagic E. coli are gram-negative, enteric bacterial pathogens. They can cause diarrhea, hemorrhagic colitis, and hemolytic uremic syndrome. This particular genome we are using in this example was based on an E. coli O157:H7 isolated from the Sakai, Japan outbreak.

Step 1: Load in sequences

The walk through will work out of the unzipped BRIG-Example.zip folder. The walk through and related figures will use C:\BRIGEXAMPLE as that location. To keep the final image consistent with the walk through, please open "ExampleProfile.xml" from the unzipped BRIG-Example.zip folder. This file configures BRIG to the same image settings in the walk through.

  1. First, set BRIGExample.fna as the reference sequence.
  2. Set <unzipped BRIG-Example folder>\genomes as the query sequence folder.
  3. Press “add to data pool", this should load several items into the pool list, there should be nine files.
  4. Set the unzipped BRIG-Example folder as the folder.
  5. The BLAST options box should be left blank.
  6. Click next

PRO TIP : Users can add individual files to the data pool too.

Step 2: Configure rings

The next step is to configure what information is shown on each of the concentric rings in BRIG. Create six rings, for each ring:

  1. Set the legend text for each ring
  2. Select the required sequences from the data pool and click on “add data” to add to the ring list.
  3. Choose a colour
  4. Set the upper (90) and lower (70) identity threshold.
  5. Click on “add new ring” and repeat steps for each new ring required.

The values required for each ring are detailed in the table below. Notice that that sequences can be collated into a single ring, like the example of K12 & HS. The ring will show BLAST matches from both HS and K12.

Legend text Required sequences Colour
GC Content GC Content Ignore
GC Skew GC Skew Ignore
Coverage BRIGExample.graph 153,0,0
O157:H7 E_coli_O157H7Sakai.gbk 0,0,153
HS and K12 E_coli_HS.fna + E_coli_K12MG1655.fna 0,153,0
CFT073 and UTI89 E_coli_CFT073.fna + E_coli_UTI89.fna 153,0,153

PRO TIP : Rings can be reordered by dragging them in the Ring List pane.

PRO TIP : You can set default threshold values in “BRIG options”. See section [sec:options] (page ) for more details.

PRO TIP : When using a Genbank/EMBL file as a reference, users can choose whether to use the protein or nucleotide sequence.

Step 3: Review and submit

The last window allows us to change the BLAST options, the location of the image file and set the image title, which will appear in the centre of the ring. For the walkthrough configure the third window as:

  1. Set the image title as “BRIG example image”.
  2. Hit submit.
  3. The image will be created in the specified output directory and should look something like Figure [BRIGfinal].

BRIG will format Genbank files, run BLAST, parse the results and render the image. The final image (Figure [BRIGfinal]) shows GC Content and Skew, the Genome coverage, contig boundaries, and the BLAST results against the other E. coli genomes. The results for HS and K12 have been collated into a single ring, likewise for UTI89 and CFT073.

PRO TIP : * Image settings, like size, fonts, etc can be configured in: Main window >Preferences >Image options..*


Comments

comments powered by Disqus