RNA-seq简介
2021-10-12 本文已影响0人
Hello育种
![](https://img.haomeiwen.com/i24959989/9ceee0899f41da3e.png)
研究对象:mRNA,其来自特定的细胞和组织(在特定条件)
RNA使用试剂盒提取后,需要-80℃保存。质量查看:28S与18S比值,接近10为好
研究方法:
![](https://img.haomeiwen.com/i24959989/960efc70daed7e51.png)
得到fastQ文件,格式:
![](https://img.haomeiwen.com/i24959989/1ea69e889aed670b.png)
解释:
![](https://img.haomeiwen.com/i24959989/6d32b0a668abe252.png)
而后需要将RNA匹配到genome上(找基因)
查看基因的表达:
![](https://img.haomeiwen.com/i24959989/d32ddb8bffd9c837.png)
2(基因大小不同)和3项(测序深度不同)表明,需要将reads数标准化
4则是为也校正p值,一般采用:
![](https://img.haomeiwen.com/i24959989/9a0c49800ff1e02c.png)
查看基因表达不同
使用DESeq包,先做PCA,查看分别,接着clustering和heatmap查看样本分组情况。需要重点研究差异表达基因,对其做GO分析。
WGCNA做全部分析
文档及例子:
https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/index.html
soft threshould:
![](https://img.haomeiwen.com/i24959989/3349ee05389e37a5.png)
需要为network选择β,
![](https://img.haomeiwen.com/i24959989/9b6ad9fcb75ec99b.png)
对不同基因进行module分组:
![](https://img.haomeiwen.com/i24959989/629344d9dc58ad8e.png)
查看各module与各性状的关系:
![](https://img.haomeiwen.com/i24959989/3219a96bd77ed875.png)
最后查看显著差异基因与module内部基因的关系:
模块成员 (MM) = 模块特征基因(模块的平均值)与基因表达谱的相关性
![](https://img.haomeiwen.com/i24959989/f7fc9920b38423cc.png)