| 24小時(shí)熱門(mén)版塊排行榜 |
| 2 | 1/1 | 返回列表 |
| 查看: 2593 | 回復(fù): 1 | |||
QQ894064647新蟲(chóng) (初入文壇)
|
[交流]
高通量測(cè)序環(huán)狀(circular RNA)生物信息分析(軟件) 已有1人參與
|
|
高通量測(cè)序環(huán)狀(circular RNA)生物信息分析 環(huán)形RNA是mRNA在剪接的過(guò)程中,上游exon的5’端與下游exon的3’端剪接到一起,從而形成的首尾相接的環(huán)狀RNA分子。早在1993 年,性別決定基因SRY的mRNA就被發(fā)現(xiàn)可以形成環(huán)狀RNA,且在睪丸中具有很高的豐度。近來(lái)的研究表明,高等動(dòng)物中環(huán)狀RNA的種類(lèi)和含量遠(yuǎn)遠(yuǎn)超過(guò)預(yù)期。 通過(guò)高通量測(cè)序分析,我們可以找到circular RNA。在測(cè)序下機(jī)數(shù)據(jù)查找circular RNA的方法主要有5種,即find_circ,circRNA_finder,CIRI,circExplorer,MapSplice。下面我們一一展示這些方法找circular RNA的命令,希望對(duì)大家的研究有所幫助。 #find_circ bowtie2 -p 16 --very-sensitive --mm -M20 --score-min=C,-15,0 -x /path/to/bowtie2_index -q -1 /path/to/file1.fq.gz -2 /path/to/file2.fq.gz | samtools view -hbuS - | samtools sort - output samtools view -hf 4 output.bam | samtools view -Sb - > unmapped.bam python unmapped2anchors.py unmapped.bam | gzip > anchors.qfa.gz bowtie2 -p 16 --reorder --mm -M20 --score-min=C,-15,0 -q -x /path/to/bowtie2_index -U anchors.qfa.gz | python find_circ.py -G /path/to/chomosomes.fa -p prefix -s find_circ.sites.log > find_circ.sites.bed 2> find_circ.sites.reads grep circ find_circ.sites.bed | grep -v chrM | python sum.py -2,3 | python scorethresh.py -16 1 | python scorethresh.py -15 2 | python scorethresh.py -14 2 | python scorethresh.py 7 2 | python scorethresh.py 8,9 35 | python scorethresh.py -17 100000 > finc_circ.candidates.bed #circRNA_finder ./runStar.pl /path/to/file1.fq /path/to/file2.fq /path/to/star_index /path/to/star/output ./postProcessStarAlignment.pl /path/to/star/output /path/to/circRNA_finder/output #CIRI bwa mem -t 16 -T 19 -v 2 /path/to/bwa_index /path/to/file1.fq /path/to/file2.fq >output.sam perl CIRI_v1.2.pl -I output.sam -O /path/to/ciri/output.txt -F /path/to/genome.fa -P #circExplorer tophat2 -a 6 --microexon-search -m 2 -p 16 -o circExplorer_output /path/to/bowtie1_index /path/to/file1.fq.gz /path/to/file2.fq.gz bamToFastq -i circExplorer_output/unmapped.bam -fq circExplorer_output/unmapped.fastq tophat2 -o circExplorer_output -p 16 --fusion-search --keep-fasta-order --bowtie1 --no-coverage-search path/to/bowtie2_index circExplorer_output/unmapped.fastq python CIRCexplorer.py -f circExplorer_output/accepted_hits.bam -g /path/to/genome.fa -r /path/to/gene_annotation_ucsc_hg19.txt -o output.txt #MapSplice python MapSplice-v2.1.8/mapsplice.py -1 /path/to/file1.fq -2 file2.fq -c /path/to/chromosome_sequences -x /path/to/bowtie1_index -p 16 -o Mapsplice_output --min-fusion-distance 200 --gene-gtf Homo_sapiens.GRCh37.66.gtf --fusion 好了,有了這些命令,大家可以隨心所欲處理高通量數(shù)據(jù)了,希望對(duì)家有幫助。如果大家對(duì)我們研究和分析感興趣,可以通過(guò)我的丁香園Id聯(lián)系我們。當(dāng)然,我們也大量生物信息方面的學(xué)習(xí)視頻,感興趣的請(qǐng)進(jìn):https://shop119322454.taobao.com/ 參考文獻(xiàn):1. Comparison of circular RNA prediction tools 2. CIRI: an efficient and unbiased algorithm for denovo circular RNA identification 3. Computational approaches for the analysis of ncRNA through deep sequencing techniques |
新蟲(chóng) (初入文壇)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|
[考研] 材料專(zhuān)業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 材料工程專(zhuān)碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|