Viral Genome Sequencing
This tutorial is about determining the genome sequence of a virus, using comparisons to a reference sequence.
Murray Valley encephalitis virus is classified as part of the Flavivirus genus of viruses, all of which are positive strand RNA viruses. The type species for this genus is yellow fever virus, and the genus includes dengue virus, west nile virus and zika virus.
Viruses from this genus have a single-segment genome of about 11 kb. Within this, there is a 10 kb open reading frame that encodes a polyprotein. Following translation, this protein is modified to yield the various, mature structural and non-structural viral proteins as illustrated in Fig. 1.
Fig. 1 - An overview of a typical flavivirus genome and the viral proteins that are part of the viral polyprotein (from Viralzone).
There are 14 complete Murray Valley encephalitis virus genome sequences available at NCBI Viral Genomes.
- Go to NCBI Viral Genomes.
Browse viral genomes by familyand click on the family Flaviviridae: Complete Genomes
- Note: you may have to widen your screen to see all the columns of viral family names.
- Look at the row for Murray Valley encephelitis virus. In the Neigbours column, we can see there are 14.
In this tutorial we will use the prototype strain 1-151 as the reference genome sequence.
- This strain was isolated in the early 1950s - see AF161266.
The isolate we are looking at has been sequenced using the Illumina platform.
- The genomic cDNA was prepared using random hexamers to prime the reverse transcription of the viral genomic RNA.
- After second strand synthesis, the cDNA was used to prepare an Illumina sequencing library and run on an Illumina MiSeq instrument.
In this activity we will use a read mapping approach to determine the sequence of a new Murray Valley encephalitis virus isolate.
- Log in to your Galaxy server
- Import files required for the activity
- View imported files
Go to the Galaxy Page
Import files to Galaxy
- Click on the
Analyze Datamenu at the top of the page.
- Click on the
Historymenu button (the on the top right of the history pane)
Import from File(at the bottom of the list)
- A new page will appear with a text box for the URL of the history to import.
- Copy the following URL into the text box:
- Galaxy will download the data files from the internet and will be available as an additional history (takes about one minute).
To make the newly imported history appear as the current history:
- Click on the View all Histories button (the on the top right of the history pane.)
- If the history has finished downloading it will appear with the title:
“imported from archive: MVEVmapping“
- Click on the
Switch tobutton above this history and then the Donebutton.
You should now have 4 files in the history pane as follows:
Reference sequence files :
MVEV.gbk- genbank format MVEV.fna- fasta format
Illumina sequence reads (R1 and R2) from the new isolate:
MVE_R1.fq- forward reads MVE_R2.fq- reverse reads
- Find variants in the isolate using the tool Snippy.
Snippy is a fast variant caller for haploid genomes. The software is available on GitHub at https://github.com/tseemann/snippy. For this activity, we are using Snippy as installed on Galaxy.
Run FastQC: How many reads in in each of the fastq files
MVEV.fnaand MVEV.gbkeach contain the genome sequence of Murray Valley encephalitis virus strain 1-151 - how many bases in the genome? (Hint: use Fasta Statistics)
Snippy maps reads from the new Murray Valley encephalitis virus isolate (the
- Find Snippy in the tool menu (in NGS: Variant Analysis)
- Select appropriate files (see screenshot below) and Execute (use default settings).
Files cataloging SNP differences:
snippy on data 4, data 3, and data 2 snps vcf file
snippy on data 4, data 3, and data 2 snps gff file
snippy on data 4, data 3, and data 2 snps table
snippy on data 4, data 3, and data 2 snps summary
A log of the progress of the run:
snippy on data 4, data 3, and data 2 log file
Regions where reads aligned (NNNN —, indicate regions where there was low or no read data):
snippy on data 4, data 3, and data 2 aligned fasta
A consensus genome sequence for the new isolate:
snippy on data 4, data 3, and data 2 consensus fasta
Summary of the read depth:
snippy on data 4, data 3, and data 2 mapping depth
A compressed version of the above files (and more that can be downloaded):
snippy on data 4, data 3, and data 2 out dir
Download these files to your local computer (click on the file name and then the disk icon in the lower left hand corner):
snippy on data 4, data 3, and data 2 out dir(and unzip)
Also download these bam files from these URLs (open each URL in a new tab and the file should download automatically):
Note: if you have previously downloaded these files, the new downloads may be renamed. Remove any spaces in the names.
- View the reads from the new isolate mapped against the reference sequence, using the tool Artemis.
Artemis is a tool to view genome sequences and mapped reads, including variants (SNPs).
If Artemis is not installed, go to http://www.sanger.ac.uk/science/tools/artemis.
View the reference sequence
- Open Artemis.
- Go to
File: Openand select MVEV.gbk. The file will probably have a “Galaxy” prefix, e.g. Galaxy2-[MVEV.gbk].genbank.
The Artemis window:
- panes 1 and 2 are the same, but can be scaled differently
- each pane has the double-stranded sequence in the centre, with amino acid translations above and below
- there is a third lower pane with feature information
- coding sequences are highlighted in blue
- other features are highlighted in green
- clicking on one of these will select it in all panels
- black vertical lines are stop codons (when zoomed out)
- move left and right with horizontal scroll bar
- zoom in and out with right-hand scroll bar
Add a plot
- Go to
Graph: Add User Plot, select snps.depth.gz
- A graph should display at the top of the screen
Producing a draft genome sequence
- Produce draft genome sequence for the new viral isolate.
Examine the mapped reads
What is the minimum read depth used by Snippy to call a SNP?
- Hint: Go back to Snippy on Galaxy - look at the information below the ‘Execute’ button
Would Snippy call a SNP at positions 1 → 5 ?
What is the maximum read depth?
What do we know about the sequence of the new isolate in the regions where there is low read coverage?
17: snippy on data 4, data 3, and data 2 out dir
This makes an “out” folder containing some files including the consensus file.
In Artemis, open
This is a file that is based on the reference sequence and includes any confirmed SNPs called by Snippy.
If we were going to use this sequence to produce the draft sequence of the new isolate, what bases would you have at positions 1→ 5?
Getting an Overview of the Difference between strain 1-151 and our new isolate
Simplest overview: view the bam file with Artemis
Once loaded, differences between reads and the reference sequence can be highlighted by right/command click in the bam view window. Select ‘Show > SNP marks’.
Getting more detail: looking at the table of SNPs
Located in the out folder there is a html file that contains a table with information about each of the SNPs called by Snippy.
Included in the table is a column that provides a prediction about the impact each SNP will have on annotated protein coding regions. The genbank file provides the annotation information used by Snippy to make the predictions.
A total of 790 SNP differences were call by Snippy
Open snps.html in your web browser
Summary: 663/790 SNPs do not result a difference in the encoded polyprotein
Is this pattern of variable genome sequence and more conserved protein sequence normal in viruses? What might be the cause?