-
Notifications
You must be signed in to change notification settings - Fork 36
/
02-ngs.Rmd
53 lines (27 loc) · 4.24 KB
/
02-ngs.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
# High throughput sequencing {#ngs}
## Three generations of sequencing technologies
First generation sequencing is Sanger sequencing. It is the technology that was used to obtain the first human genome sequence.
Second generation sequencing is also called next generation sequencing (NGS) and is the start of high throughput sequencing. It is what scientists use most often nowadays, and Illumina is the market leader. Most of the rest of this course will cover data analysis using second generation sequencing.
Third generation sequencing is single-molecule sequencing. There are many new technologies still under active development, although none has reached market penetration.
<iframe width="560" height="315" src="https://www.youtube.com/embed/fVExBR1x4QU" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
## FASTQ and FASTQC
NGS generates FASTQ files. FASTQC is an computational approach to evaluate the quality of your NGS data.
<iframe width="560" height="315" src="https://www.youtube.com/embed/T2QAGtuvq8I" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
## Early sequence alignment (1 with 1)
In the early days (1970s), scientists were not worried about having to align too many sequences. They wanted to find the best alignment between two sequences. Many bioinformatics courses start with learning these, although it is not the main focus of our course. We included two videos in case you are interested.
The Needlemen-Wunsch algorithm is the earliest algorithm to find the alignment between two sequences and score their similarity.
<iframe width="560" height="315" src="https://www.youtube.com/embed/ipp-pNRIp4g" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
When two sequences are long, and only a portion of them can align well with each other, the Smith-Waterman algorithm can find the best local sequence alignment. It is still considered the best alignment approach, although it is slow.
<iframe width="560" height="315" src="https://www.youtube.com/embed/lu9ScxSejSE" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
## Sequence search algorihtms (1 with many)
With more and more sequences available in the public in the 1980s, scientists were interested in finding whether their newly sequenced string has been sequenced before in the public database. Therefore, the fast search algorithm BLAST was developed, using one sequence as the query to find similar sequences from a database.
<iframe width="560" height="315" src="https://www.youtube.com/embed/Ge00nXqIlow" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
## Burrows-Wheeler Aligner (many with many)
With NGS, scientists need much faster search (aka mapping) algorithms in order to align the millions of sequences to the reference genome. The current best algorithm is called Borrow-Wheeler Aligner or BWA.
In order to understand BWA, we first need to introduce Borrows-Wheeler transformation and LF mapping
<iframe width="560" height="315" src="https://www.youtube.com/embed/gqM3j2IRQH4" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
The basic idea of Borrows-Wheeler alignment
<iframe width="560" height="315" src="https://www.youtube.com/embed/RXjc4wDV9BE" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>
## Alignment output
NGS raw data is in FASTQ. Alignment gives you SAM (alignment) or BAM (binary version of SAM) files which contain the sequence information in FASTQ and the mapping locations. BED file is the simpliest, although there is information loss.
<iframe width="560" height="315" src="https://www.youtube.com/embed/nZG6t_NobFQ" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>