【哈佛大学:计算生物学 & 生物信息学】学习记录(二)

笔记主要内容:1~3代测序技术,fastq文件 & FASTQC

(1)测序技术

1、一代测序技术:Sanger Sequencing

测序条件:需要有足够的量的单链DNA,即相同序列需要达到多少数量才能进行测序。测序过程:以一条链作为模板,DNA聚合酶将环境中的材料(dNTP & ddNTP),进行结合,即合成另一条链(ddNTP结合之后,DNA合成反应终止)。下图所示,SEQUENCE (END)代表对应DNA序列的最后一个碱基是什么。

image

参考阅读资料:https://zhuanlan.zhihu.com/p/94183808

2、二代测序技术:Illumina Sequencing Cluster Generation

产业巨头:Illumina(e.g. NovaSeq 6000)
测序过程:对DNA序列加上接头序列 & 引物,经桥式PCR,不同的DNA序列形成不同的cluster。在反应环境中加入ATCG(荧光),摄像机对不同的光信号进行捕捉,记录对应的碱基信息,随后洗去目前碱基的荧光组分,进行下一步的DNA链合成反应。
【需要注意的技术】桥式PCR

参考视频:https://www.youtube.com/watch?v=fCd6B5HRaZ8

3、三代测序技术:Single Molecule Sequencing

技术代表:PacBio & Nanopore(两者技术不同)
与NGS对比,其不需要对序列进行扩增。同时,三代测序技术的测序片段长度也远远长于二代测序技术。
参考视频:

(2)Fastq & FASTQC

老生常谈了,生物信息学入门必备知识:fastq & 用于检查测序文件质量的fastqc

1、FastQ文件格式

  1. Sequence ID
  2. Sequence
  3. Quality ID
  4. Quality score

图示:

image

质量值表示方式:Phred quality(一般为Phred 33),其数学含义代表某一个碱基测错的概率,计算公式为-10log_{10}Pr

质量值从低到高:

image

关于fastq文件,需要注意的是:

  • 开头几个bp,测定结果是不稳定的,因为该测序阶段机器正在预热。
  • 最后的几个bp,可能也会出现问题(测序环境问题)

2、为什么要进行质量控制?

  • 检查公司返回的测序结果,如果不行当然要求重测了~
  • 为了后续分析的准确性

3、FASTQC结果文件解读

这边当然得推荐一下,生信菜鸟团的FASTQC笔记:http://www.bio-info-trainee.com/95.html

Per Base Sequence Quality

虽然测序片段是250bp,但是经FASTQC查看过后,最终想保留下来用于后续分析的reads长度可能就缩减到了150~200bp。
下图类似一个箱线图,用于表示测序片段某一个位置的测序质量范围。

image
Per Sequence Quality Distribution

碱基对应质量值的分布(主要集中在什么质量区间,比如Q20~Q30)

image
Nucleotide Content Per Position

对应的就是100bp的序列在每一个位置上的碱基占比(A,T,C,G)。如果是在全基因组水平进行测序,每一种类型的碱基其占比应在25%左右,若非上述情况(比如链特异性测序,则碱基占比会改变)

下图“Good quality”的情况,可以选择将每一条序列的前几个bp给trim掉。


image
Per Sequence GC Content

此部分就是FASTQC将实际的GC content与期望的GC content进行对比,如果像下图右边一样,就代表测序可能出现了一些问题(e.g. 样品被污染)

参考阅读资料:https://www.bioinformatics.babraham.ac.uk/projects/fastqc/Help/3%20Analysis%20Modules/5%20Per%20Sequence%20GC%20Content.html

image
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,588评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,456评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,146评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,387评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,481评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,510评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,522评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,296评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,745评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,039评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,202评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,901评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,538评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,165评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,415评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,081评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,085评论 2 352