CircSplice軟件可變剪切實(shí)現(xiàn)原理

CircSplice軟件可變剪切實(shí)現(xiàn)原理

(1) 通過STAR比對至參考基因組序列，得到chimeric reads；

(2) 保留那些比對上同一染色體和鏈的chimeric reads，通過CIGAR值判斷reads的坐標(biāo)；

(3) 識別BSJ，去除另一端落在供體和受體位點(diǎn)之外的reads，以確保另一端也來自circRNA；

(4) 通過剪切位點(diǎn)GT-AG或者CT-AC對BSJ進(jìn)行過濾和識別；

(5) 對BSJ進(jìn)行基因注釋；

(6) 將BSJ和其另一端reads比對到基因的外顯子或內(nèi)含子，結(jié)合基因注釋進(jìn)行4種剪切類型的判斷。允許2bp容錯。另外要求代表剪切的read的junction一端坐標(biāo)與注釋外顯子的一端坐標(biāo)一致，以減少假陽性預(yù)測；

(7) 通過GT-AG或CT-AC原則對可變剪切事件進(jìn)行過濾；

(8) 計(jì)算可變剪接事件的reads數(shù)，并進(jìn)行標(biāo)準(zhǔn)化（支持可變剪切的reads數(shù)*106/總chimeric reads數(shù)）；

(9) 根據(jù)坐標(biāo)融合不同樣本的剪切事件，比較樣本間的環(huán)狀RNA剪切事件；允許2bp容錯；