RNA-seq入门实战(零):RNA-seq流程前的准备——Linux与R的环境创建

RNA-seq入门实战整体分析流程

前言:
进行RNA-seq入门实战首先需要有一定的linux与R基础,推荐跟着B站生信技能树-jimmy老师学习打牢基础:
【生信技能树】生信人应该这样学linux(更新至第14集)_哔哩哔哩_bilibili
【生信技能树】生信人应该这样学R语言_哔哩哔哩_bilibili

本节概览:

  1. Linux下RNA-seq环境创建:
    Ubuntu子系统下载安装、Mniconda3与上游分析软件下载
  2. R下RNA-seq环境创建
    R与Rstudio下载安装、Bioconductor与R包下载

1. Linux环境设置

1.1 Linux系统的创建——Ubuntu

运行Linux系统一般使用服务器或者个人电脑的虚拟机(Virtualbox、VMware)和子系统,下面简单介绍Windows子系统的安装配置,
详细说明请参阅Windows子系统WSL的体验与配置——Ubuntu-22.04

  • Ubuntu子系统的下载安装
    首先在win10中搜索“ 启用或关闭Windows功能 ”,进入该程序,勾选“适用于Linux的Windows子系统”;之后去微软商店搜索Ubuntu下载安装,一般安装默认版本或者最新的22.04LTS



  • 用户权限设置
    设置好用户名和密码进入Ubuntu后,需要设置一下用户权限
    启用root需要设置密码:sudo passwd root
    添加用户至root组: usermod -aG sudo username
    切换root用户:su
    退出root:exit
  • 软件镜像源设置
    一般选择国内的清华镜像,ubuntu | 镜像站使用帮助 | 清华大学开源软件镜像站 | Tsinghua Open Source Mirror,Ubuntu 的软件源配置文件是 /etc/apt/sources.list。将系统自带的该文件做个备份,再编辑(* 注意要先切换为root用户*才能编辑该文件)
su
cd /etc/apt/
cp sources.list sources.list.bak
vim sources.list

Ubuntu版本 22.04 LTS的设置(注意要找到对应的Ubuntu版本设置)如下:

# 默认注释了源码镜像以提高 apt update 速度,如有需要可自行取消注释
deb https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy main restricted universe multiverse
# deb-src https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy main restricted universe multiverse
deb https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-updates main restricted universe multiverse
# deb-src https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-updates main restricted universe multiverse
deb https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-backports main restricted universe multiverse
# deb-src https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-backports main restricted universe multiverse
deb https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-security main restricted universe multiverse
# deb-src https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-security main restricted universe multiverse
# 预发布软件源,不建议启用
# deb https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-proposed main restricted universe multiverse
# deb-src https://mirrors.tuna.tsinghua.edu.cn/ubuntu/ jammy-proposed main restricted universe multiverse
  • 更新一下所有软件
sudo  apt update       #更新可用软件包列表
sudo  apt upgrade      #更新已安装的包

1.2 Mniconda3下载安装

一般使用Mniconda3软件进行创建分析环境和管理软件,下面简单介绍Mniconda3的安装,其详细使用说明请参阅Miniconda3的安装、配置和使用

  • 下载与安装
    切换到安装位置(一般为主目录~), 下载最新的miniconda3,bash启动安装,一直enter、yes就可以了。
wget https://mirrors.ustc.edu.cn/anaconda/miniconda/Miniconda3-latest-Linux-x86_64.sh
#wget -c https://repo.continuum.io/miniconda/Miniconda3-latest-Linux-x86_64.sh
bash  Miniconda-latest-Linux-x86_64.sh
  • 设置conda软件镜像源
    依次输入以下命令设置软件镜像源,并展示镜像源地址,这里设置了中科大镜像:
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/pkgs/pro
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/pkgs/r
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/pkgs/free
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/pkgs/main
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/cloud/bioconda
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/cloud/conda-forge
conda config  --set show_channel_urls yes
  • 创建分析环境
    创建rna流程环境rna_p3,指定python版本为3,同时下载sra-tools 软件,之后每次要运行rna流程都要进入该环境
conda create -n rna_p3 python=3 sra-tools               
conda env list                                         #查看环境
conda activate rna_p3                                  #进入conda 环境
conda deactivate                                    #退出当前conda环境
  • 上游分析软件下载
    本次RNA-seq流程涉及软件如下所示,使用conda install -y 软件名=版本号依次进行下载即可,可同时下载多个软件,但不建议同时下载太多,容易报错

质控清洗:fastqc multiqc trim-galore
比对计数: hisat2 subread samtools=1.6 salmon


2. R环境设置

注意,安装使用R首先需要你的电脑用户名不能含有中文,否则后期运行程序会很容易报错,如果用户名含有中文请进行更改(很难。。。)或新建一个用户;
Rstudio是R的编译器,能提升用户交互体验,安装Rstudio前一定要先安装R

2.1 R与Rstudio下载安装

  • 先在清华镜像源下载 R,地址:The Comprehensive R Archive Network (tsinghua.edu.cn)
    依次进入:Download R for Windows——base——Previous releases;
    一般选择下载R 4.1.3,使用最新版可能会有不适配的问题,使用默认选项一直确定安装即可

  • 安装好 R后之后,再去Rstudio官网下载Rstudio,地址:Download the RStudio IDE - RStudio

  • Rstudio中设置清华镜像源
    依次进入:Tools——Global Options——Pakages——Change;选择镜像源


2.2 Bioconductor与R包下载

新建一个Rscript,运行以下内容,下载Bioconductor和本次实战所需所有R包等,代码修改自jimmy老师

###Bioconductor 下载
install.packages("BiocManager") 

###设置好清华镜像
rm(list = ls())   
options()$repos 
options()$BioC_mirror
#options(BioC_mirror="https://mirrors.ustc.edu.cn/bioc/") 
options(BioC_mirror="http://mirrors.tuna.tsinghua.edu.cn/bioconductor/")
options("repos" = c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/"))
options()$repos 
options()$BioC_mirror

###安装需要的包
BiocManager::install(c("GSEABase","GSVA","msigdbr","clusterProfiler" ),ask = F,update = F)
BiocManager::install(c("GEOquery","limma","impute" ),ask = F,update = F)
BiocManager::install(c("org.Hs.eg.db","org.Mm.eg.db"),ask = F,update = F)
BiocManager::install(c("DESeq2","edgeR" ),ask = F,update = F)
BiocManager::install("enrichplot",ask = F,update = F)
BiocManager::install("devtools",ask = F,update = F)
BiocManager::install("WGCNA",ask = F,update = F) 
BiocManager::install("data.table",ask = F,update = F)
BiocManager::install("tximport",ask = F,update = F)
BiocManager::install("tidyverse",ask = F,update = F)
BiocManager::install("DOSE",ask = F,update = F)
BiocManager::install("patchwork",ask = F,update = F)
BiocManager::install("RBGL",ask = F,update = F)  #Vennerable依赖包
BiocManager::install("pathview",ask = F,update = F)
BiocManager::install(c("STRINGdb","ggraph","igraph"),ask = F,update = F)
install.packages("Vennerable", repos="http://R-Forge.R-project.org") #安装Vennerable包
install.packages("statmod")

#其他一些基础包安装
options()$repos 
install.packages(c("FactoMineR", "factoextra")) 
install.packages(c("ggplot2", "pheatmap","ggpubr","ggthemes","ggstatsplot","ggsci","ggsignif")) 
install.packages("rvcheck")
(.packages())  #查看当前加载运行包

#更新所有包
rvcheck::update_all(check_R = F,
                    which = c("CRAN", "BioC", "github"))
  • 永久设置bioconductor镜像
#bioconductor | 镜像站使用帮助 | 清华大学开源软件镜像站 | Tsinghua Open Source Mirror
file.edit('~/.Rprofile')

输入以下内容:

# 清华源的镜像
options(BioC_mirror="https://mirrors.tuna.tsinghua.edu.cn/bioconductor")

以上就是进行RNA-seq流程前的全部准备了,下一步就可以开始进行上游数据下载、格式转化和质控清洗等步骤啦

参考资料
GitHub - jmzeng1314/GEO (强烈推荐学习)
【生信技能树】生信人应该这样安装软件_哔哩哔哩_bilibili
【生信技能树】生信分析入门环境搭建_哔哩哔哩_bilibili


RNA-seq实战系列文章:
RNA-seq入门实战(零):RNA-seq流程前的准备——Linux与R的环境创建
RNA-seq入门实战(一):上游数据下载、格式转化和质控清洗
RNA-seq入门实战(二):上游数据的比对计数——Hisat2+ featureCounts 与 Salmon
RNA-seq入门实战(三):从featureCounts与Salmon输出文件获取counts矩阵
RNA-seq入门实战(四):差异分析前的准备——数据检查
RNA-seq入门实战(五):差异分析——DESeq2 edgeR limma的使用与比较
RNA-seq入门实战(六):GO、KEGG富集分析与enrichplot超全可视化攻略
RNA-seq入门实战(七):GSEA——基因集富集分析
RNA-seq入门实战(八):GSVA——基因集变异分析
RNA-seq入门实战(九):PPI蛋白互作网络构建(上)——STRING数据库的使用
RNA-seq入门实战(十):PPI蛋白互作网络构建(下)——Cytoscape软件的使用
RNA-seq入门实战(十一):WGCNA加权基因共表达网络分析——关联基因模块与表型

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

推荐阅读更多精彩内容