lncRNA组装流程的软件介绍之Stringtie

StringTie 是用于 RNA-seq 的转录本组装和定量软件，StringTie 可以看

做是cufflinks软件的升级版本，其功能和Cufflinks是一样的，包括下面两个主要功能：转录本组装和定量；相比Cuffinks, 其运行速度更快。

conda install stringtie

输入的文件必须是一个根据基因组位置排好序的BAM文件，可以是Tophat或Hisat2的输出文件，在通过samtools排序

注意：一定要使用-dta选项来运行HISAT2，否则结果将会受到影响。

作为选项，可以向StringTie提供GTF / GFF3格式的参考注释基因组文件。在这种情况下，StringTie更喜欢使用注释文件中的这些“已知”基因，对于那些被表达的基因，它将计算coverage，TPM和FPKM值。它还会产生额外的转录本，而注释文件中并没有这些转录本。请注意，如果不使用选项-e，那么参考转录本就需要被reads 完全覆盖，以便包含在StringTie的输出中。在这种情况下，其他通过StringTie从数据中组装的转录本，且不在注释文件中的转录本也会输出。

-o [<path/>]<out.gtf> # 设置StringTie组装转录本的输出GTF文件的路径和文件名。此处可指定完整路径，在这种情况下，将根据需要创建目录。默认情况下，StringTie将GTF写入标准输出。

-p <int>  # 指定组装转录本的线程数（CPU）。默认值是1

-G <ref_ann.gff>  # 使用参考注释基因文件指导组装过程，格式GTF/GFF3。输出文件中既包含已知表达的转录本，也包含新的转录本。选项-B，-b，-e，-C需要此选项（详情如下）

--rf  # 链特异性建库方式：fr-firststrand(最常用的是dUTP测序方式，其他有NSR，NNSR).

--fr  # 链特异性建库方式：fr-secondstrand(如 Ligation,Standard SOLiD).

-l <label>  # 将<label>设置为输出转录本名称的前缀。默认：STRG

-f <0.0-1.0>  #将预测转录本的最低isoform的丰度设定为在给定基因座处组装的丰度最高的转录本的一部分。较低丰度的转录物通常是经加工的转录本的不完全剪接前体的artifacts。默认值为0.1。

-m <int>  # 设置预测的转录本所允许的最小长度.默认值为200

-A <gene_abund.tab> # 输出基因丰度的文件（制表符分隔格式）

-C <cov_refs.gtf> # 输出所有转录本对应的reads覆盖度的文件，此处的转录本是指参考注释基因文件中提供的转录本。(需要参数 -G).

-a <int>  # Junctions that don't have spliced reads that align across them with at least this amount of bases on both sides are filtered out. Default: 10

-j <float>  # 连接点的覆盖度，即设置至少有这么多的spliced reads 比对到连接点(align across a junction)。 这个数字可以是分数, 因为有些reads可以比对到多个地方。 当一个read 比对到 n 个地方是，则此处连接点的覆盖度为1/n 。默认值为1。

-t  # 该参数禁止修剪组装的转录本的末端。默认情况下，StringTie会根据组装的转录本的覆盖率的突然下降来调整预测的转录本的开始和/或停止坐标。

-c <float>  # 设置预测转录本所允许的最小read 覆盖度。 当一个转录本的覆盖度低于阈值，则输出文件中不含该转录本。默认值为 2.5

-g <int>  #设置ga最小值。Reads that are mapped closer than this distance are merged together in the same processing bundle. Default: 50 (bp)

-B  # 应用该选项，则会输出Ballgown输入表文件（* .ctab），其中包含用-G选项给出的参考转录本的覆盖率数据。（有关这些文件的说明，请参阅Ballgown文档。）
    如果选项-o 给出输出转录文件的完整路径，则* .ctab文件与输出GTF文件在相同的目录下。

-b <path>   # 指定 *.ctab 文件的输出路径, 而非由-o选项指定的目录。
    # 注意: 建议在使用-B/-b选项中同时使用-e选项，除非StringTie GTF输出文件中仍需要新的转录本。

-e  # 限制reads比对的处理，仅估计和输出与用-G选项给出的参考转录本匹配的组装转录本。使用该选项，则会跳过处理与参考转录本不匹配的组装转录本，这将大大的提升了处理速度。

--merge #转录本合并模式。在合并模式下，StringTie将所有样品的GTF/GFF文件列表作为输入，并将这些转录本合并/组装成非冗余的转录本集合。这种模式被用于新的差异分析流程中，用以生成一个跨多个RNA-Seq样品的全局的、统一的转录本。如果提供了-G选项（参考注释基因组文件），则StringTie将从输入的GTF文件中将参考转录本组装到transfrags中。(个人理解：transfrags可能指的是拼接成更大的转录本片段，tanscript fragments)。在此模式下可以使用以下附加选项：
-G <guide_gff> #参考注释基因组文件(GTF/GFF3)
-o <out_gtf>   #指定输出合并的GTF文件的路径和名称 (默认值：标准输出)
-m <min_len>   #合并文件中，指定允许最小输入转录本的长度 (默认值: 50)
-c <min_cov>   #合并文件中，指定允许最低输入转录本的覆盖度(默认值: 0)
-F <min_fpkm>  #合并文件中，指定允许最低输入转录本的FPKM值 (默认值: 0)
-T <min_tpm>   #合并文件中，指定允许最低输入转录本的TPM值 (默认值: 0)
-f <min_iso>   #minimum isoform fraction (默认值: 0.01)
-i             #合并后，保留含retained introns的转录本 (默认值: 除非有强有力的证据，否则不予保留)
-l <label>     #输出转录本的名称前缀 (默认值: MSTRG)

ls /home/data/lihe/lncRNA_project/04.mapping/*.bam > 1
ls /home/data/lihe/lncRNA_project/04.mapping/*.bam | cut -d "/" -f 7 | cut -d "_" -f 1 | cut -d "." -f 1 > 0
paste 0 1 > config

cat > 05.stringtie.sh

config=$1
number1=$2
number2=$3
gtf=/home/data/lihe/reference/human/gtf/gencode.v37.annotation.gtf
cat $1 | while read id
do
    if((i%$number1==$number2))
    then
    arr=(${id})
    input=${arr[1]}
    sample=${arr[0]}
    stringtie -B -p 4 -G $gtf -o ./${sample}.stringtie.gtf -A ./${sample}.tab -l ${sample} $input
    fi    ## end for number1
    i=$((i+1))
done

for i  in {0..4}
do 
(nohup bash 05.stringtie.sh  config 5 $i 1>log.$i.txt 2>&1 & )
done

-G $gtf # 提供的参考gtf文件，指导组装 
-o ./${sample}.stringtie.gtf # 输出文件名
-l ${sample} # 输出的转录本前缀名
$input # 输入bam文件
-A ./${sample}.tab # 对输出的gtf统计基因表达量，并以一个tab分割的文件输出，这里需要提交输出的文件名{sample}

ls ../01.stringtie_gtf/*.gtf > mergelist.txt
gtf=/home/data/lihe/reference/human/gtf/gencode.v36.annotation.gtf
nohup stringtie --merge -p 8 -G $gtf -o ./stringtie_merged.gtf ./mergelist.txt > merge.log 2>&1 &

ls /home/data/lihe/lncRNA_project/04.mapping/*.bam > 1
ls /home/data/lihe/lncRNA_project/04.mapping/*.bam | cut -d "/" -f 7 | cut -d "_" -f 1 | cut -d "." -f 1 > 0
paste 0 1 > config

cat > 05.expression_stringtie.sh

config=$1
number1=$2
number2=$3
cat $1 | while read id
do
    if((i%$number1==$number2))
    then
    arr=(${id})
    input=${arr[1]}
    sample=${arr[0]}
    stringtie -e -B -p 4 -G ~/lncRNA_project/05.stringtie/02.merge_gtf/stringtie_merged.gtf -A ./${sample}_gene_abundances.tsv -o ./${sample}/${sample}.gtf  $input

    fi    ## end for number1
    i=$((i+1))
done

for i  in {0..4}
do 
(nohup bash 05.expression_stringtie  config 5 $i 1>log.$i.txt 2>&1 & )
done

ls -d SRR*|while read id;do stringtie -e -A $id/$id'_gene_abund.tab' -C $id/$id'_cov_refs.gtf' -B -p 10 -G $gtf -o $id/$id'.merged.gtf' $id/$id'.sorted.bam' ;done

1、 GTF文件： 记录组装的转录本信息

2、 Tab文件： 记录基因丰度信息

3、 GTF文件：完全覆盖与参考注释基因组文件所匹配的转录本信息

4、 *.ctab文件：用于下游Ballgown软件做差异表达分析的输入文件

5、 GTF文件：在合并模式下，生成一个合并的GTF文件

如果StringTie与 -C <cov_refs.gtf> 选项一起运行（需要选项-G

如果StringTie与-B选项一起运行，它将返回Ballgown输入文件，包含以下文件：
(1) e2t.ctab 
(2) e_data.ctab
(3) i2t.ctab
(4) i_data.ctab
(5) t_data.ctab。

如果StringTie使用--merge选项运行，它将多个GTF / GFF文件作为输入，并将这些转录本合并和组装成非冗余转录本集合。