lncRNA-seq数据分析之新lncRNA鉴定和注释视频课程众筹

前面我系统性的总结了：lncRNA的一些基础知识 ，和lncRNA芯片的一般分析流程 ，还有LncRNA-seq的一般分析流程 ,里面提到了一个目前非常小众的分析方向，就是新lncRNA鉴定和注释，因为大部分人研究的物种的human或者mouse，已经被分析的很透彻了，encode计划等资源非常丰富，很少需要鉴定新的lncRNA。

不过对于其它物种，猫狗猪，甚至其它你叫不出来名字的昆虫，鱼类，这个分析策略还是蛮常见的。比如发表在Front. Genet., 18 March 2019 | https://doi.org/10.3389/fgene.2019.00196的文章

• Transcriptome Analysis Suggests the Roles of Long Intergenic Non-coding RNAs in the Growth Performance of Weaned Piglets

就是重新下载一个公共数据，然后进行新lncRNA鉴定和注释，分析部分主要是分成两个大块，首先是hisat2+stringtie流程，然后是组装好的gtf文件的后，细致的进行新lncRNA鉴定和注释。

LncRNA-seq数据分析的两个部分

分析流程如下：

新lncRNA鉴定和注释图解流程

前面的hisat2+stringtie流程流程很简单

就是参考：猪狗的参考基因组构建索引，还有使用ebi数据库直接下载fastq测序数据  ，做好准备工作，然后使用conda安装一些软件，建立好目录

conda create -n lncRNA
conda activate lncRNA
conda install -y -c  bioconda hisat2 stringtie samtools fastp   gffcompare
# conda search gffcompare
mkdir 0.qc 1.raw_fq 2.clean_fq 3.hisat2_bams 4.stringtie_gtfs 5.lncRNA

流程基本上3个软件，衔接一些即可

conda activate lncRNA
index=/home/jmzeng/reference/genome/pig/pig_hisat2
gtf=/home/jmzeng/reference/genome/pig/Sus_scrofa.Sscrofa11.1.99.chr.gtf

fastp -i 1.raw_fq/${id}_1.fastq.gz \
                    -o 2.clean_fq/${id}_1.fastp.fq.gz \
                    -I 1.raw_fq/${id}_2.fastq.gz \
                    -O 2.clean_fq/${id}_2.fastp.fq.gz \
                    -l 36 -q 20 --compression=6 \
                    -R ${id} -h ${id}.html
fq1=2.clean_fq/${id}_1.fastp.fq.gz
fq2=2.clean_fq/${id}_2.fastp.fq.gz
hisat2 -p 4 -x $index -1  $fq1 -2  $fq2 | \
samtools sort -@ 4  -o 3.hisat2_bams/$sample.bam -

stringtie -p 4 -G $gtf \
            -o 4.stringtie_gtfs/$sample.gtf  \
            -l  $sample 3.hisat2_bams/$sample.bam                 

当然，你需要自己去搜索理解软件的参数啦。

后面的新lncRNA鉴定和注释还是蛮耗费时间的

而且不同物种的新lncRNA鉴定和注释细节还不一样，不同的gtf文件版本可以对比印证。

我们研发的步骤是：

新lncRNA鉴定和注释的具体步骤

完整课程思维导图在：https://mubu.com/doc/ISk-Ev1tg

课程录制需要一些反馈和动力，所以采取众筹模式，吸纳部分真正有兴趣的朋友进入微信群参与讨论哈。（毕竟新lncRNA鉴定和注释是一个小众方向，大部分朋友就是看个热闹）

还是老规矩，18.8元进群，一个简单的门槛，隔绝那些营销号！同时，我们也会在群里共享一些lncRNA相关资料，仅此而已，考虑清楚哦！ 

长按识别二维码

添加微信

支付18.8元入学习群

烦请备注姓名学校单位信息

众筹模式

再次强调，你完全无需参与众筹，视频会完全免费共享在B站！！！

• https://space.bilibili.com/338686099/#/

只不过是B站这个途径不方便共享这些学习素材，所以大家统一添加我们的加群小助手二维码然后进入微信群，我们录制视频过程会选择性跟大家互动，整理好资料后腾讯微云发送给群里的朋友！

(注意，不是答疑，不是售后，也不讲解Linux和R基础知识)，你需要自己跟着我们生信技能树的系统性基础入门视频学习背景知识！

首先是LINUX学习

我在《生信分析人员如何系统入门Linux(2019更新版)》把Linux的学习过程分成6个阶段 ，提到过每个阶段都需要至少一天以上的学习：

• 第1阶段：把linux系统玩得跟Windows或者MacOS那样的桌面操作系统一样顺畅，主要目的就是去可视化，熟悉黑白命令行界面，可以仅仅以键盘交互模式完成常规文件夹及文件管理工作。
• 第2阶段：做到文本文件的表格化处理，类似于以键盘交互模式完成Excel表格的排序、计数、筛选、去冗余，查找，切割，替换，合并，补齐，熟练掌握awk,sed,grep这文本处理的三驾马车。
• 第3阶段：元字符，通配符及shell中的各种扩展，从此linux操作不在神秘！
• 第4阶段：高级目录管理：软硬链接，绝对路径和相对路径，环境变量
• 第5阶段：任务提交及批处理，脚本编写解放你的双手
• 第6阶段：软件安装及conda管理，让linux系统实用性放飞自我

然后是R学习

我在在生信分析人员如何系统入门R(2019更新版) 里面给初学者的知识点路线图如下：

• 了解常量和变量概念
• 加减乘除等运算（计算器）
• 多种数据类型（数值，字符，逻辑，因子）
• 多种数据结构（向量，矩阵，数组，数据框，列表）
• 文件读取和写出
• 简单统计可视化
• 无限量函数学习

文末友情宣传

强烈建议你推荐给身边的博士后以及年轻生物学PI，多一点数据认知，让他们的科研上一个台阶：

• 全国巡讲全球听（买一得五） ，你的生物信息学入门课
• 生信技能树的2019年终总结  ，你的生物信息学成长宝藏