De novo Aktif Olarak Çevrilmiş Açık Okuma Çerçevelerinin Ribozom Profil Oluşturma Verileri ile Tanımlanması

Özet

Ribozomların çevrilmesi, kodon başına üç nükleotidin peptitlere şifresini çözer. Ribozom profillemesi ile yakalanan mRNA boyunca hareketleri, karakteristik üçlü periyodiklik sergileyen ayak izlerini üretir. Bu protokol, tüm transkriptom düzeyinde aktif olarak çevrilmiş açık okuma çerçevelerini tanımlamak için ribozom profil oluşturma verilerinden bu belirgin özelliği deşifre etmek için RiboCode'un nasıl kullanılacağını açıklar.

Özet

Açık okuma çerçevelerinin (ORF'ler), özellikle küçük peptitleri kodlayan ve belirli fizyolojik bağlamlar altında aktif olarak çevrilenlerin tanımlanması, bağlama bağlı translatomların kapsamlı ek açıklamaları için kritik öneme sahiptir. RNA üzerindeki ribozomları çevirmenin bağlanma yerlerini ve yoğunluklarını tespit etmek için kullanılan bir teknik olan ribozom profilleme, genom çapında çevirinin nerede gerçekleştiğini hızlı bir şekilde keşfetmek için bir yol sunar. Bununla birlikte, biyoinformatikte, ribozom profillemesi için çeviri yapan ORF'leri verimli ve kapsamlı bir şekilde tanımlamak önemsiz bir görev değildir. Burada, ribozom profil oluşturma verilerindeki bozuk ve belirsiz sinyallerden herhangi bir boyuttaki ORF'leri aktif olarak çevirmek için tasarlanmış RiboCode adlı kullanımı kolay bir paket açıklanmaktadır. Daha önce yayımlanmış veri kümemizi örnek alan bu makalede, ham verilerin önceden işlenmesinden nihai çıktı sonucu dosyalarının yorumlanmasına kadar tüm RiboCode işlem hattı için adım adım yönergeler sağlanmaktadır. Ayrıca, açıklamalı ORF'lerin çeviri oranlarını değerlendirmek için, her bir ORF'deki ribozom yoğunluklarının görselleştirilmesi ve ölçülmesi için prosedürler de ayrıntılı olarak açıklanmaktadır. Özetle, bu makale çeviri, küçük ORF'ler ve peptitlerle ilgili araştırma alanları için yararlı ve zamanında bir talimattır.

Giriş

Son zamanlarda, giderek artan bir çalışma grubu, kodlayan genlerin ORF'lerinden ve daha önce açıklanmış genlerden kodlanmamış olarak çevrilmiş peptitlerin yaygın bir şekilde üretildiğini ortaya koymuştur, örneğin uzun kodlamayan RNA'lar (lncRNA'lar)1,2,3,4,5,6,7,8. Bu çevrilmiş ORF'ler, çevresel değişikliklere, strese ve hücre farklılaşmasına yanıt vermek için hücreler tarafından düzenlenir veya uyarılır1,8,9,10,11,12,13. Bazı ORF'lerin çeviri ürünlerinin, gelişim ve fizyolojideki çeşitli biyolojik süreçlerde önemli düzenleyici roller oynadığı gösterilmiştir. Örneğin, Chng ve ^ark.14, kardiyovasküler gelişim için kritik olan Elabela (Ela, Apela / Ende / Toddler olarak da bilinir) adlı bir peptid hormonu keşfetti. Pauli ve ark., Ela'nın ayrıca erken balık embriyosunda hücre göçünü teşvik eden bir mitojen görevi gördüğünü öne ^sürdü15. Magny ve ark., kalsiyum taşınımını düzenleyen ve Drosophila kalbindeki düzenli kas kasılmasını etkileyen 30'dan az amino asitten oluşan iki mikropeptit ^{bildirmiştir10}.

Bu tür peptitlerin genom tarafından kaç tane kodlandığı ve biyolojik olarak alakalı olup olmadıkları belirsizliğini korumaktadır. Bu nedenle, bu potansiyel olarak kodlayan ORF'lerin sistematik olarak tanımlanması oldukça arzu edilir. Bununla birlikte, evrimsel koruma16,17 ve kütle spektrometrisi18,19 gibi geleneksel yaklaşımları kullanarak bu ^ORF'lerin (yani protein veya peptit) ürünlerini doğrudan belirlemek zordur, çünkü her iki yaklaşımın da tespit verimliliği, üretilen proteinlerin veya peptitlerin uzunluğuna, bolluğuna ve amino asit bileşimine bağlıdır. Nükleotid çözünürlüğünde mRNA'lar üzerindeki ribozom doluluğunu tanımlamak için bir teknik olan ribozom profillemesinin ortaya çıkışı, uzunluklarına ve bileşimlerine bakılmaksızın farklı transkriptlerin kodlama potansiyelini değerlendirmek için kesin bir yol ^{sağlamıştır3,20,21}. Ribozom profillemesini kullanarak aktif olarak çevrilen ORF'leri tanımlamak için önemli ve sık kullanılan bir özellik, ribozomun mRNA üzerindeki ayak izlerinin başlangıç kodonundan durdurma kodonuna kadar üç nükleotid (3-nt) periyodikliğidir. Bununla birlikte, ribozom profilleme verileri genellikle ORF'ler boyunca düşük ve seyrek dizileme okumaları, yüksek dizileme gürültüsü ve ribozomal RNA (rRNA) kontaminasyonları dahil olmak üzere çeşitli sorunlara sahiptir. Bu nedenle, bu tür veriler tarafından üretilen çarpık ve belirsiz sinyaller, ribozomların mRNA üzerindeki ayak izlerinin 3-nt periyodiklik kalıplarını zayıflatır ve bu da sonuçta yüksek güvenilirliğe sahip çevrilmiş ORF'lerin tanımlanmasını zorlaştırır.

"RiboCode" adlı bir paket, ORF'nin çerçeve dışı RPF'lerden önemli ölçüde daha fazla kare içi ribozom korumalı parçaya (RPF) sahip olup olmadığını incelemek için değiştirilmiş bir Wilcoxon imzalı sıralama testi ve P değeri entegrasyon stratejisini ^uyarladı22. Simüle edilmiş ve gerçek ribozom profil oluşturma verilerinde translatomun de novo ek açıklaması için son derece verimli, hassas ve doğru olduğu gösterilmiştir. Burada, önceki çalışma tarafından oluşturulan ham ribozom profil oluşturma sıralama veri kümelerinden potansiyel çeviri ORF'lerini tespit etmek için bu aracın nasıl kullanılacağını ^{açıklıyoruz23}. Bu veri kümeleri, kontrol (si-Ctrl) ve EIF3E (si-eIF3e) küçük müdahale eden RNA'lar (siRNA'lar) ile transfekte edilen MCF-10A hücrelerinin ribozom doluluk profillerini karşılaştırarak çeviride EIF3 alt birimi "E" nin (EIF3E) işlevini araştırmak için kullanılmıştır. RiboCode'u bu örnek veri kümelerine uygulayarak, potansiyel olarak küçük peptitleri veya proteinleri kodlayan 5.633 yeni ORF tespit ettik. Bu ORF'ler, yukarı akış ORF'leri (uORF'ler), aşağı akış ORF'leri (dORF'ler), üst üste binen ORF'ler, yeni protein kodlayan genlerden (yeni PCG'ler) ORF'ler ve yeni protein kodlamayan genlerden (yeni NonPCG'ler) ORF'ler dahil olmak üzere kodlama bölgelerine göre konumlarına göre çeşitli tiplerde kategorize edildi. uORF'ler üzerindeki RPF okuma yoğunlukları, kontrol hücrelerine kıyasla EIF3E eksikliği olan hücrelerde önemli ölçüde artmıştır, bu da en azından kısmen aktif olarak çevrilen ribozomların zenginleşmesinden kaynaklanabilir. EIF3E eksikliği olan hücrelerin 25. ila 75. kodonu arasındaki bölgede lokalize ribozom birikimi, erken evrede translasyon uzamasının tıkandığını göstermiştir^. Bu protokol ayrıca, tanımlanan ORF'ler üzerindeki ribozom ayak izlerinin 3-nt periyodiklik modellerini incelemek için istenen bölgenin RPF yoğunluğunun nasıl görselleştirileceğini de göstermektedir. Bu analizler, RiboCode'un ORF'lerin çevirisini tanımlamadaki ve çevirinin düzenlenmesini incelemedeki güçlü rolünü göstermektedir.

Protokol

1. Ortam kurulumu ve RiboCode kurulumu

1. Bir Linux terminal penceresi açın ve bir conda ortamı oluşturun:
   conda create -n RiboCode python=3.8
2. Oluşturulan ortama geçin ve RiboCode ve bağımlılıklarını yükleyin:
   conda RiboCode'u etkinleştir
   conda yüklemek -c bioconda ribocode ribominer sra-tools fastx_toolkit cutadapt papyon yıldız samtools

2. Veri hazırlama

1. Genom referans dosyalarını alın.
   1. Referans sırası için, https://www.ensembl.org/index.html adresindeki Ensemble web sitesine gidin, üst menüyü İndir'i ve sol taraftaki FTP İndir menüsünü tıklayın. Sunulan tabloda, DNA (FASTA) sütununda ve Türlerin İnsan olduğu satırda FASTA'ya tıklayın. Açılan sayfada, Homo_sapiens bağlantısını kopyalayın. GRCh38.dna.primary_assembly.fa.gz dosyasını indirin, ardından terminalde indirin ve sıkıştırmasını açın:
      wget -c \
      http://ftp.ensembl.org/pub/release-104/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz
      gzip -d Homo_sapiens. GRCh38.dna.primary_assembly.fa.gz
   2. Başvuru ek açıklaması için, son açılan web sayfasındaki Gen kümeleri sütununda GTF'yi sağ tıklatın. Homo_sapiens bağlantısını kopyalayın. GRCh38.104.gtf.gz ve aşağıdakileri kullanarak indirin:
      wget -c \
      http://ftp.ensembl.org/pub/release-104/gtf/homo_sapiens/Homo_sapiens.GRCh38.104.gtf.gz
      gzip -d Homo_sapiens. GRCh38.104.gtf.gz
      NOT: GTF dosyasının üç seviyeli bir hiyerarşide düzenlenmiş genom ek açıklamaları içerdiğinden Ensemble web sitesinden alınması önerilir, yani her gen ekzonlar ve isteğe bağlı çeviriler içeren transkriptler içerir (örneğin, kodlama dizileri [CDS], çeviri başlangıç sitesi, çeviri bitiş sitesi). Bir genin veya transkriptin ek açıklamaları eksik olduğunda, örneğin UCSC veya NCBI'den elde edilen bir GTF dosyası, tam üst-alt hiyerarşi ek açıklamalarıyla güncellenmiş bir GTF oluşturmak için GTFupdate'i kullanın: GTFupdate original.gtf > updated.gtf. .gff biçimindeki ek açıklama dosyası için, AGAT ^toolkit24'ü veya .gtf biçimine dönüştürmek üzere başka bir aracı kullanın.
2. rRNA dizileri alın.
   1. UCSC Genom Tarayıcısını https://genome.ucsc.edu'da açın ve Araçlar | Açılır listede Tablo Tarayıcısı .
   2. Açılan sayfada, klade için Memeli , genom için İnsan , grup için Tüm Tablolar , tablo için rmask ve bölge için genom değerlerini belirtin. Filtre için, yeni bir sayfaya gitmek üzere Oluştur'u tıklatın ve repClass'ı rRNA ile eşleştiği gibi ayarlayın.
   3. Gönder'e tıklayın ve ardından çıktı biçimini hg38_rRNA.fa olarak sıraya ve çıktı dosya adına ayarlayın. Son olarak, tıklayın Çıktı al | rRNA dizisini almak için dizi alın.
3. Sıralı Okuma Arşivi'nden (SRA) ribozom profil oluşturma veri kümeleri alın.
   1. si-eIF3e tedavi grubunun çoğaltılmış örneklerini indirin ve yeniden adlandırın:
      fastq-dump SRR9047190 SRR9047191 SRR9047192
      mv SRR9047190.fastq si-eIF3e-1.fastq
      mv SRR9047191.fastq si-eIF3e-2.fastq
      mv SRR9047192.fastq si-eIF3e-3.fastq
   2. Denetim grubunun çoğaltma örneklerini indirin ve bunları yeniden adlandırın:
      fastq-dump SRR9047193 SRR9047194 SRR9047195
      mv SRR9047193.fastq si-Ctrl-1.fastq
      mv SRR9047194.fastq si-Ctrl-2.fastq
      mv SRR9047195.fastq si-Ctrl-3.fastq
      NOT: Bu örnek veri kümeleri için SRA katılım kimlikleri, GSE131074 araması yapılarak Gene Expression Omnibus (GEO) web ^sitesinden25 alınmıştır.

3. Adaptörleri kesin ve rRNA kontaminasyonunu giderin

1. (İsteğe bağlı) Bağdaştırıcıları sıralama verilerinden çıkarın. Bu örnekte olduğu gibi bağdaştırıcı dizileri zaten kırpılmışsa bu adımı atlayın. Aksi takdirde, adaptörleri okumalardan kırpmak için cutadapt kullanın.
   i in si-Ctrl-1 si-Ctrl-2 si-Ctrl-3 si-eIF3e-1 si-eIF3e-2 si-eIF3e-3 için
   yapmak
   cutadapt -m 15 --match-read-wildcards -a CTGTAGGCACCATCAAT \
   -o ${i}_trimmed.fastq ${i}.fastq
   yapılmış
   NOT: -a parametresinden sonraki adaptör dizisi, cDNA kütüphanesinin hazırlanmasına bağlı olarak değişecektir. 15'ten kısa okumalar ( -m ile verilir) atılır, çünkü ribozom korumalı parçalar genellikle bu boyuttan daha uzundur.
2. Aşağıdaki adımları kullanarak rRNA kontaminasyonunu kaldırın:
   1. İndeks rRNA referans dizileri:
      bowtie-build -f hg38_rRNA.fa hg38_rRNA
   2. rRNA'dan kaynaklanan okumaları ekarte etmek için okumaları rRNA referansıyla hizalayın:
      i in si-Ctrl-1 si-Ctrl-2 si-Ctrl-3 si-eIF3e-1 si-eIF3e-2 si-eIF3e-3 için
      yapmak
      papyon -n 0 -y -a --norc --en iyi --tabaka -S -p 4 -l 15 \
      --un=./${i}_noncontam.fastq hg38_rRNA -q ${i}.fastq ${i}.aln
      yapılmış
      -p , görevlerin paralel olarak çalıştırılması için iş parçacığı sayısını belirtir. RPF okumalarının nispeten küçük boyutu göz önüne alındığında, bildirilen hizalamaların en iyi olduğunu garanti etmek için diğer bağımsız değişkenler (örneğin, -n, -y, -a, -norc, --best, --strata ve -l) belirtilmelidir. Daha fazla ayrıntı için Bowtie web ^sitesine26 bakın.

4. Temiz okumaları genomla hizalayın

1. Bir genom indeksi oluşturun.
   mkdir STAR_hg38_genome
   STAR --runThreadN 8 --runMode genomeGenerate --genomeDir ./STAR_hg38_genome --genomeFastaFiles Homo_sapiens. GRCh38.dna.primary_assembly.fa --sjdbGTFfile Homo_sapiens. GRCh38.104.gtf
2. Temiz okumaları (rRNA kontaminasyonu yok) oluşturulan referansla hizalayın.
   i in si-Ctrl-1 si-Ctrl-2 si-Ctrl-3 si-eIF3e-1 si-eIF3e-2 si-eIF3e-3 için
   yapmak
   STAR --runThreadN 8 --outFilterType Normal --outWigType wiggle --outWigStrand Stranded --outWigNorm RPM --outFilterMismatchNmax 1 --outFilterMultimapNmax 1 --genomeDir STAR_hg38_genome --readFilesIn ${i}_noncontam.fastq --outFileNamePrefix ${i}. --outSAMtype BAM SortedByCoordinate --quantMode TranscriptomeSAM GeneCounts --outSAMattributes Tümü
   yapılmış
   NOT: Ters transkriptaz27 tarafından her okunuşun 5' ucuna sıklıkla şablonsuz bir nükleotid eklenir ve varsayılan olarak yumuşak kırpma gerçekleştirdiği için STAR tarafından verimli bir şekilde kesilir. STAR parametreleri STAR ^{kılavuzunda28} açıklanmıştır.
3. Hizalama dosyalarını sıralayın ve dizine ekleyin.
   i in si-Ctrl-1 si-Ctrl-2 si-Ctrl-3 si-eIF3e-1 si-eIF3e-2 si-eIF3e-3 için
   yapmak
   samtools sort -T ${i}. Aligned.toTranscriptome.out.sorted \
   -o ${i}. Aligned.toTranscriptome.out.sorted.bam \
   ${i}. Aligned.toTranscriptome.out.bam
   samtools dizini ${i}. Aligned.toTranscriptome.out.sorted.bam
   samtools dizini ${i}. Aligned.sortedByCoord.out.bam
   yapılmış

5. RPF'lerin boyut seçimi ve P-bölgelerinin tanımlanması

1. Transkript ek açıklamalarını hazırlayın.
   prepare_transcripts -g Homo_sapiens. GRCh38.104.gtf \
   -f Homo_sapiens. GRCh38.dna.primary_assembly.fa -o RiboCode_annot
   NOT: Bu komut, GTF dosyasından mRNA transkriptlerinin gerekli bilgilerini toplar ve FASTA dosyasından tüm mRNA transkriptleri için dizileri çıkarır (her transkript, GTF dosyasında tanımlanan yapılara göre ekzonları birleştirerek birleştirilir).
2. Belirli uzunluklardaki RPF'leri seçin ve P sitesi konumlarını belirleyin.
   i in si-Ctrl-1 si-Ctrl-2 si-Ctrl-3 si-eIF3e-1 si-eIF3e-2 si-eIF3e-3 için
   yapmak
   metaplots -a RiboCode_annot -r ${i}. Aligned.toTranscriptome.out.bam \
   -o ${i} -f0_percent 0,35 -pv1 0,001 -pv2 0,001
   yapılmış
   NOT: Bu komut, her uzunluktaki hizalanmış okumaların 5' ucunun toplam profillerini, açıklamalı çeviri başlangıç (veya durdurma) kodonları etrafında çizer. Okuma uzunluğuna bağımlı P-sitesi, ana okumaların 5' uçları ile başlangıç kodonu arasındaki ofset mesafelerinin dağılım grafikleri (örneğin, Şekil 1B) incelenerek manuel olarak belirlenebilir. RiboCode ayrıca her örnek için, önemli 3-nt periyodiklik desenleri gösteren okumaların P-site konumlarının otomatik olarak belirlendiği bir yapılandırma dosyası oluşturur. -f0_percent, -pv1 ve -pv2 parametreleri, okuma çerçevesinde zenginleştirilmiş RPF okumalarını seçmek için oran eşiğini ve p-değeri kesmelerini tanımlar. Bu örnekte, 29, 30 ve 31 nt okumaların 5' ucundaki +12, +13 ve +13 nükleotidler her yapılandırma dosyasında el ile tanımlanmıştır.
3. Her örnek için yapılandırma dosyalarını düzenleme ve birleştirme
   NOT: Benzersiz ORF'lerden oluşan bir fikir birliği kümesi oluşturmak ve sonraki analizleri gerçekleştirmek üzere yeterli okuma kapsamını sağlamak için, önceki adımdaki tüm örneklerin seçilen okumaları birleştirilir. merged_config.txt dosyada (Ek Dosya 1) tanımlanan belirli uzunluklardaki okumalar ve bunların P-site bilgileri, bir sonraki adımda ORF'lerin çeviri potansiyelini değerlendirmek için kullanılır.

6. De novo ek açıklama çeviri ORF'leri

1. RiboCode'u çalıştırın.
   RiboCode -a RiboCode_annot -c merged_config.txt -l evet -g \
   -o RiboCode_ORFs_result -s ATG -m 5 -A CTG,GTG,TTG
   Bu komutun önemli parametreleri aşağıdaki gibidir:
   -c, giriş dosyalarının yolunu ve seçilen okumaların ve bunların P-sitelerinin bilgilerini içeren yapılandırma dosyası.
   -l, durdurma kodonlarının yukarısında birden fazla başlangıç kodonu bulunan transkriptler için, çeviri potansiyellerini değerlendirmek için en uzun ORF'lerin (en distal başlangıç kodonundan durdurma kodonuna kadar olan bölge) kullanılıp kullanılmadığı. Hayır olarak ayarlanırsa, başlangıç kodonları otomatik olarak belirlenir.
   -s, ORF'lerin tanımlanması için kullanılan kanonik başlangıç kodon(lar)ı.
   -A, (isteğe bağlı olarak) ORF tanımlaması için kullanılan kanonik olmayan başlangıç kodonları (örneğin, insan için CTG, GTG ve TTG), diğer türlerin mitokondrilerinde veya çekirdeğinde farklılık ^{gösterebilir29}.
   -m, ORF'lerin minimum uzunluğu (yani amino asitler).
   -o, tahmin edilen ORF'lerin ayrıntılarını içeren çıktı dosya adının öneki (Ek Dosya 2).
   -g ve -b, tahmin edilen ORF'leri sırasıyla gtf veya yatak formatına çıkarır.

7. (İsteğe bağlı) ORF ölçümü ve istatistikleri

1. Her ORF'de RPF okumalarını sayın.
   i in si-Ctrl-1 si-Ctrl-2 si-Ctrl-3 si-eIF3e-1 si-eIF3e-2 si-eIF3e-3 için
   yapmak
   ORFcount -g RiboCode_ORFs_result_collapsed.gtf \
   -r ${i}. Aligned.sortedByCoord.out.bam -f 15 -l 5 -m 25 -M 35 \
   -o ${i}_ORF.counts -s yes -c intersection-strict
   yapılmış
   NOT: ORF'lerin başlangıç ve bitişlerinin etrafında biriken potansiyel ribozomları hariç tutmak için, ilk 15'te (-f ile belirtilir) ve son 5 kodonda (-l ile spesifik) ayrılan okuma sayısı sayılmaz. İsteğe bağlı olarak, sayılan RPF'lerin uzunlukları 25 ila 35 nt (RPF'lerin ortak boyutları) aralığıyla sınırlıdır.
2. RiboCode kullanarak tespit edilen ORF'lerin temel istatistiklerini hesaplayın:
   Rscript RiboCode_utils. R
   NOT: RiboCode_utils. R (Ek Dosya 3), RiboCode çıktısı için bir dizi istatistik sağlar, örneğin, tanımlanan ORF'lerin sayısını sayma, ORF uzunluklarının dağılımını görüntüleme ve normalleştirilmiş RPF yoğunluklarını hesaplama (yani, RPKM, eşlenen milyon okuma başına kilobaz başına okuma).

8. (İsteğe bağlı) Tahmin edilen ORF'lerin görselleştirilmesi

1. İstediğiniz ORF için başlangıç ve bitiş kodonlarının göreceli konumlarını (örneğin, ENSG00000100902_35292349_35292552_67) transkriptinde RiboCode_ORFs_result_collapsed.txt'dan (Ek dosya 3) alın. Ardından, ORF'de RPF okumalarının yoğunluğunu çizin:
   plot_orf_density -a RiboCode_annot -c merged_config.txt -t ENST00000622405 \
   -s 33 -e 236 --start-kodon ATG -o ENSG00000100902_35292349_35292552_67
   Burada - s ve - e, ORF'yi çizmenin çeviri başlangıç ve bitiş konumunu belirtir. -- start-codon , şekil başlığında görünecek olan ORF'nin başlangıç kodonunu tanımlar. -o , çıktı dosyası adının önekini tanımlar.

9. (İsteğe bağlı) RiboMiner kullanarak metagen analizi

NOT: Aşağıdaki adımları izleyerek, EIF3E knockdown'ın tanımlanan açıklamalı ORF'lerin çevirisi üzerindeki etkisini değerlendirmek için metagen analizini gerçekleştirin:

1. RiboMin için transkript ek açıklamaları oluşturun, bu da RiboCode tarafından oluşturulan ek açıklama dosyasına dayanarak her gen için en uzun transkriptleri ayıklar (adım 5.1).
   OutputTranscriptInfo -c RiboCode_annot/transcripts_cds.txt \
   -g Homo_sapiens. GRCh38.104.gtf -f RiboCode_annot/transcripts_sequence.fa \
   -o longest.transcripts.info.txt -O all.transcripts.info.txt
2. RiboMiner için yapılandırma dosyasını hazırlayın. RiboCode'un metaplots komutu tarafından oluşturulan yapılandırma dosyasını kopyalayın (adım 5.4) ve "RiboMiner_config.txt" olarak yeniden adlandırın. Ardından, Ek dosya 4'te gösterilen biçime göre değiştirin.
3. RiboMiner kullanarak metagen analizleri
   1. Transkriptler arasında RPF'lerin yoğunluklarının toplu ve ortalama profilini oluşturmak için MetageneAnalysis'i kullanın.
      MetageneAnalysis -f RiboMiner_config.txt -c longest.transcripts.info.txt \
      -o MA_normed -U kodon -M RPKM -u 100 -d 400 -l 100 -n 10 -m 1 -e 5 --norm evet \
      -y 100 --tipi UTR
      Önemli parametreler şunlardır: --type, CDS veya UTR bölgelerini analiz etme; --norm, okuma yoğunluğunun normalleştirilip normalleştirilmediği; -y, her transkript için kullanılan kodon sayısı; -U, RPF yoğunluğunu kodon seviyesinde veya nt seviyesinde çizin; -u ve -d, kodonu başlatmak veya kodonu durdurmak için analiz bölgelerinin aralığını tanımlar; -l, CDS'nin minimum uzunluğu (yani kodon sayısı); -M, transkript filtreleme modu, ya sayım ya da RPKM; -n minimum sayım veya analiz için CDS'de RPKM. -m normalleştirilmiş bölgede CDS'nin minimum sayıları veya RPKM'si; -e, normalleştirilmiş bölgeden hariç tutulan kodon sayısı.
   2. Kontrol hücrelerinde ve eIF3 eksikliği olan hücrelerde mRNA üzerindeki ribozom doluluklarını karşılaştırmak için bir dizi pdf dosyası oluşturun.
      PlotMetageneAnalysis -i MA_normed_dataframe.txt -o MA_normed \
      -g si-Ctrl,si-eIF3e -r si-Ctrl-1,si-Ctrl-2,si-Ctrl-3__si-eIF3e-1,si-eIF3e-2,si-eIF3e-3 -u 100 -d 400 --mod ortalama
      NOT: PlotMetageneAnalysis pdf dosyaları kümesini oluşturur. MetageneAnalysis ve PlotMetageneAnalysis kullanımı ile ilgili ayrıntılar RiboMiner web ^sitesinde30 mevcuttur.

Sonuçlar

Örnek ribozom profil oluşturma veri kümeleri, GEO veritabanında GSE131074 katılım numarası altında depolanmıştır. Bu protokolde kullanılan tüm dosyalar ve kodlar Ek dosyalar 1-4'ten edinilebilir. RiboCode'u bir dizi yayınlanmış ribozom profilleme veri ^kümesine23 uygulayarak, kontrol ve EIF3E siRNA'ları ile tedavi edilen MCF-10A hücrelerinde aktif olarak çevrilen yeni ORF'leri tanımladık. Büyük olasılıkla çeviri ...

Tartışmalar

Ribozom profilleme, ribozomların hücrelerdeki etkisini genom ölçeğinde incelemek için benzeri görülmemiş bir fırsat sunar. Ribozom profilleme verileri tarafından taşınan bilgilerin tam olarak deşifre edilmesi, hangi gen bölgelerinin veya transkriptlerin aktif olarak çevrildiğine dair fikir verebilir. Bu adım adım protokol, ribozom profil oluşturma verilerini paket yükleme, veri hazırlama, komut yürütme, sonuç açıklaması ve veri görselleştirme dahil olmak üzere ayrıntılı olarak analiz etm...

Malzemeler

Name	Company	Catalog Number	Comments
A computer/server running Linux	Any	-	-
Anaconda or Miniconda	Anaconda	-	Anaconda: https://www.anaconda.com; Miniconda:https://docs.conda.io/en/latest/miniconda.html
R	R Foundation	-	https://www.r-project.org/
Rstudio	Rstudio	-	https://www.rstudio.com/