Skip to content

Commit 85c7718

Browse files
committed
Added some material on cross-contamination
Added acknowledgements Moved figures to its own folder Modified CSS to remove shadows for images
1 parent aae7272 commit 85c7718

20 files changed

+97
-18
lines changed

Cuernavaca/slides/Cuernavaca.md

+95-15
Original file line numberDiff line numberDiff line change
@@ -1,16 +1,18 @@
11
# Secuenciación de próxima generación para la evolución del virus
22
<!--- Next-generation sequencing analysis for virus evolution -->
3-
## Art Poon [@GitHub](github.com/ArtPoon)
3+
## Art Poon [@GitHub](https://github.com/ArtPoon)
44

5-
#### Departments of Pathology and Laboratory Medicine; Microbiology and Immunology; Applied Mathematics at Western University, Ontario, Canada
5+
#### Department of Pathology and Laboratory Medicine; Microbiology and Immunology; and Applied Mathematics
6+
7+
![](images/WesternLogo.png)
68

79

810

911
# Materiales del seminario
1012

1113
* All slides, code and data are available at
1214
[http://github.com/PoonLab/courses](http://github.com/PoonLab/courses)
13-
![GitHub](github-logo.png)
15+
<center>![GitHub](images/github-logo.png)</center>
1416

1517
* Download a copy from the releases page.
1618
* Both `*.tar.bz` and `*.zip` archives should be available.
@@ -39,7 +41,7 @@
3941

4042
# Cumplir con la l&iacute;nea de comandos
4143

42-
![](es.tar.png)
44+
<center>![](images/es.tar.png)</center>
4345

4446

4547

@@ -106,24 +108,35 @@ CCGGGATCTTGTTGATTGTGAACGCTGCGGTACCTAAGGATGACACGCCTTTCAATCCATGTTTGTCCGTT
106108

107109

108110

109-
# Significaci&oacute;n cl&iacute;nica
110-
111-
* Every subtype B sample in a MiSeq run had a substantial frequency (~3%) of resistance mutation E138A
112-
* Quality scores were fine!
111+
# Otros problemas
113112

114-
![Frequency of E138A declined from 3% to neglible once we removed bad cycles](mutation-barplot-v2.png)
113+
1. Bad tile/cycle combinations
114+
2. Longitudinal (between-run) contamination
115+
3. Cross-sectional (within-run) contamination
115116

116117

117118

118119
# Mala azulejos / ciclos
119-
120-
![Bad tile-cycle combinations in a MiSeq run](bad-cycles.png)
120+
<center>
121+
![Bad tile-cycle combinations in a MiSeq run](images/bad-cycles.png)
122+
</center>
121123

122124
* Not necessarily a big problem for random libraries
123125
* Can be a really big problem for amplicon libraries!
124126

125127

126128

129+
# Significaci&oacute;n cl&iacute;nica
130+
131+
* Every subtype B sample in a MiSeq run had a substantial frequency (~3%) of resistance mutation E138A
132+
* Quality scores were fine!
133+
134+
<center>
135+
![Frequency of E138A declined from 3% to neglible once we removed bad cycles](images/mutation-barplot-v2.png)
136+
</center>
137+
138+
139+
127140
# Archivos de InterOp
128141

129142
* These are binary files (not human readable) that store quality and error metrics associated with the run.
@@ -152,7 +165,49 @@ less ErrorMetricsOut.bin
152165

153166
# Visualizar las tasas de error
154167

155-
![alt text](ErrorMetricsOut.png)
168+
![Plot of error metrics](images/ErrorMetricsOut.png)
169+
170+
171+
172+
# Contaminac&oacute;n cruzada longitudinal
173+
174+
* Reads from previous runs appear in the current run.
175+
176+
* Used modified PCR primer with extra T (TAATG, start 130926), alternated primers between runs.
177+
178+
<center>![Measures of longitudinal run contamination](images/carryover1.png)</center>
179+
<small>LC Swenson et al., Tropism Testing by MiSeq Is Comparable to 454-Based Methods but Exhibits Contamination Issues; Abstract #611, Conference on Retroviruses and Opportunistic Infections, Boston, MA, March 3-6, 2014</small>
180+
181+
182+
183+
# &iexcl;Blanquear al rescate!
184+
185+
* A lab member raised issue [publicly](http://seqanswers.com/forums/archive/index.php/t-29110.html) and with Illumina
186+
187+
* Eventually a revised wash protocol with bleach was [posted](https://drive.google.com/file/d/0B383TG7oJh2CTWp6MzZkMHBDd28/edit) by Illumina
188+
189+
<center>
190+
![Reduction of longitudinal run contamination with bleach wash](images/carryover2.png)
191+
</center>
192+
193+
194+
195+
# Contaminaci&oacute;n cruzada en el interior
196+
197+
* We were running hepatitis C virus (HCV) and human leukocyte antigen (HLA) loci
198+
samples on the same run, different indices
199+
200+
* Some reads in HCV sample were mapping to HLA, and vice versa!
201+
202+
203+
204+
# Asociados con &iacute;ndices compartidos
205+
206+
![](images/cross-contamination.png)
207+
208+
* Can be resolved by using unique indices!
209+
<small>CJ Brumme et al., Within-Run Cross-Contamination in Deep Sequencing Applications on the Illumina MiSeq; Abstract #592, Conference on Retroviruses and Opportunistic Infections, Seattle, WA, February 23-26, 2015</small>
210+
156211

157212

158213

@@ -202,6 +257,7 @@ bowtie2 -x data/zika -1 data/Zika-envelope.n1E4.R1.fastq.gz \
202257
```
203258

204259

260+
205261
# Muy poco lecturas mapeadas a la referencia
206262
```
207263
10000 reads; of these:
@@ -363,8 +419,9 @@ Updated reference with 10 differences
363419

364420
# Echemos un vistazo al resultado de nuevo
365421

366-
![](tablet-updated.png)
422+
![](images/tablet-updated.png)
367423

424+
* Try opening the new SAM file in Tablet
368425

369426

370427
# &iexcl;Esto no es &uacute;til!
@@ -394,7 +451,9 @@ python scripts/slice-sam.py -refname NC_012532.1 \
394451

395452
# &iexcl;Mucho mejor!
396453

397-
![](slice-1700-2000.png)
454+
<center>
455+
![](images/slice-1700-2000.png)
456+
</center>
398457

399458

400459

@@ -412,8 +471,18 @@ sandbox/slice-1700-2000.nwk
412471

413472
#
414473

415-
![A tree](tree-small.png)
474+
<center>
475+
![A tree](images/tree-small.png)
476+
</center>
477+
478+
479+
480+
# &iquest;Qu&eacute; pasa con la metagen&oacute;mica viral?
481+
482+
* We are already doing this with multiple virus targets (HIV, hepatitis C virus)
416483

484+
* Need to be cautious about mapping to similar references because adaptive algorithm
485+
can cause collisions
417486

418487

419488

@@ -428,3 +497,14 @@ or clinical purposes.
428497
* I have skipped many steps, such as read trimming.
429498

430499
* All materials are free to use under the GNU Affero GPL license.
500+
501+
502+
503+
# Expresiones de gratitud
504+
505+
* The laboratory results I have cited were obtained by the *BC Centre for Excellence in HIV/AIDS*
506+
(Chanson J. Brumme, Luke C. Swenson, Hope Lapointe, P. Richard Harrigan)
507+
508+
* My new lab is supported by grants from Genome Canada and the Canadian Institutes of Health Research
509+
510+
* **&iexcl;Posiciones disponibles!**

Cuernavaca/slides/bad-cycles.png

-182 KB
Binary file not shown.

Cuernavaca/slides/custom.css

+2-3
Original file line numberDiff line numberDiff line change
@@ -246,10 +246,9 @@ body {
246246
* IMAGES
247247
*********************************************/
248248
.reveal section img {
249-
margin: 15px 0px;
249+
margin: 0px 0px;
250250
background: rgba(255, 255, 255, 0.12);
251-
border: 4px solid #fff;
252-
box-shadow: 0 0 10px rgba(0, 0, 0, 0.15); }
251+
border: 4px solid #fff; }
253252

254253
.reveal section img.plain {
255254
border: 0;

Cuernavaca/slides/github-logo.png

-30.3 KB
Binary file not shown.
11.7 KB
Loading
27.1 KB
Loading
161 KB
Loading
26.3 KB
Loading
44.4 KB
Loading

Cuernavaca/slides/images/cihr.png

8.43 KB
Loading
49.7 KB
Loading
File renamed without changes.
16.8 KB
Loading
43.7 KB
Loading
360 KB
Loading
File renamed without changes.
-105 KB
Binary file not shown.

Cuernavaca/slides/slice-1700-2000.png

-80.4 KB
Binary file not shown.

0 commit comments

Comments
 (0)