$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Virüs stokları oluşturmak 1.
Not: Bu protokol ıslak tezgah yönünün bir akış şemasıdır, Şekil 1 'de gösterilen viral hazır üretimi ve doku kültürü hücrelerinin sonraki enfeksiyon detayları genel retrovirüslerin farklı için geçerli olacaktır.. Bazı deneyler için, hedef hücre, endojen virüs reseptörü (ler) i eksprese olabilir, ve bu gibi durumlarda heterolog viral zarf glikoproteinini barındıran pseudotyped retroviral parçacıklar yapılar örn veziküler stomatit virüsü (VSV-G) G glikoprotein olacak enfeksiyon 44,45 için gerekli.
Not: önlemlerin alınması gerekir, HIV-1 ile çalışırken. Belirli kurallar kurumdan kuruma değişir rağmen, tüm virüs tabanlı çalışma (tipik olarak bir doku kültürü kaputu olarak anılacaktır) özel bir operatör kısıtlı biyolojik güvenlik kabini yapılmalıdır. Uygun kişisel koruyucu donanımlarO yüz koruması, galoş, bir çift eldiven tabakası ve bir tam vücut tulum elbise her zaman giyilmelidir içerir. Virüs ile ilgili deneyler kaynaklanan tüm sıvı atık çamaşır suyu (% 10 nihai konsantrasyon) ile inaktive edilmelidir ve katılar dahil olmak üzere tüm atık önce otoklavlanmalıdır.
- Bir gün, transfeksiyondan önce,% 10 (h / h) cenin sığır serumu ve% 1 ile takviye edilmiş Dulbecco Modified Eagle Medium (DMEM), 10 ml plakası 3.3 x 10 6 HEK293T hücreleri (h / h), penisilin / streptomisin (10,000 U / beş 100 mm yemekleri her ml stok).
Not: desteklenmiş DMEM bu noktadan sonra DMEM-FPS olarak ifade edilir. - Sonraki gün, ticari olarak temin edilebilir transfeksiyon reaktifleri veya kalsiyum fosfat kullanılarak VSV-G sentezleme konstruktunun 1 ug ile zarf-silinmiş tek tur vektörler 10 tam uzunlukta retroviral moleküler klonlar taşıyan plazmidin ug veya 9 ug hücreleri transfekte.
- c inkübe% 5 CO2 ile nemlendirilmiş bir hücre kültürü inkübatöründe 37 ° C'de arşın (bu durum buradan sonra "doku kültür inkübatörü" olarak anılacaktır). Yaklaşık olarak 48 saat, volümetrik pipet kullanarak virüs içeren hücre ortamı hasat ve yerçekimi akışı ile 0.45 um filtreden geçirilir.
- 4 ° C 'de 1 saat süre ile 200,000 x g'de ultrasantrifugasyon sureti ile virüsü konsantre edilir. 20 u ihtiva eden DNaz 500 ul DMEM-FPS virüs pelletini, ve 37 ° C'de 1 saat boyunca inkübe edilir.
Not: DNaz adım transfeksiyon prosedüründen devam plazmid DNA'nın yükünü ortadan kaldırarak istenmeyen plazmid dizilerinin kurtarma azaltmaya yardımcı olur.
- Üreticinin talimatlarına göre HIV-1 p24 antijen yakalama kiti kullanılarak p24 konsantrasyonu 46 belirleyin.
Not: Virüs konsantrasyonu da ters transkriptaz aktivite deneyi 47,48 belirlenebilir. Alternatif olarak, fonksiyonel virüs düzeyi canMOI ölçülmesiyle belirlenebilir. Bu en kolay gibi gelişmiş yeşil floresan proteini gibi flüoresan raportör genleri ifade virüslerle ayırma floresans ile aktive edilmiş hücre kullanılarak yapılır. Optimize edilmiş hücre çizgileri gibi enfeksiyon aynı düzeyde desteklemez primer hücreler çalışırken MOI belirlenmesi özellikle faydalı olabilir.
Virüs ile 2. Infect Hücreleri
- Plaka 3.0 X 10 5 HEK293T 2.5 ml DMEM-FPS 6 oyuklu bir plaka içerisinde oyuk başına hücre ve doku kültürü kuluçka makinesi içinde bir gece boyunca inkübe edilir.
Not: Bu protokol ile geri entegre olarak sitelerinin sayısı kötü kullanılan hücreler ve aktif virüs miktarı sayısı ile doğru orantılıdır. - bir doku kültürü kuluçka makinesi içinde 2 saat süre ile 500 ul taze DMEM-FPS bir nihai hacim içinde 500 ng / ml değerinde nihai viral p24 konsantrasyonuna sahip hücreleri enfekte, daha sonra 2 ml DMEM-FPS oyuk başına 37 ° C önceden ısıtılmış ekleme ve inkübasyon devam edin.
- atEnfeksiyondan 48 saat, ortamını çıkarın ve 2 mi fosfat-tamponlu salin (PBS) hücreleri yıkamak. 0.5 ml tripsin-EDTA, 37 ° C'ye kadar önceden ısıtılmış ve birkaç saniye sonra görsel olarak hücre yer değiştirme için kuyu kontrol ekleyin.
- 2 ml önceden ısıtılmış DMEM-FPS ekleyin ve ~ bir hacimsel pipet yardımıyla 10 kez / aşağı pipet yukarı yumuşak hücreleri tekrar süspansiyon. DMEM-FPS önceden ısıtılmış 18 ml içeren 75 cm2'lik doku kültür şişesi çözüm aktarın ve bir doku kültürü kuluçka makinesi içinde hücreleri inkübe edin.
- enfeksiyonun baştan az beş gün sonra, 5 ml PBS ile yıkayın, ortam kaldırarak hücreleri toplamak 2 mi önceden ısıtılmış tripsin-EDTA ekleyin ve pipetleme 5 ml önceden ısıtılmış DMEM-FPS ile yeniden süspanse edin. 2.500 x g'de oda sıcaklığında 5 dakika boyunca çözelti santrifüj ve süpernatant atılır.
Not: yaklaşık 48 saat sonrası enfeksiyon 49,50 olarak bu koşullar platolarda altında bütünleşme, kültür ek 3 gün sufficientl için gerekli olmasına rağmenY, hücre bazlı bir DNA rekombinasyon ya da viral aracılı autointegration sonucu olmamış DNA moleküllerinin seyreltilmesi. - Bir ticari olarak temin edilebilen bir kit (örneğin, 51 bakınız) kullanılarak hücre topağından genomik DNA ekstrakte edin. 10 mM Tris-HCl, 200 ul, pH 8.5 ile birlikte, iyon değişim kolonundan DNA Zehir.
Not: hücrelerin bir alikotu önce NGS uygun virüs enfeksiyonu sağlamak için enfekte analizi için 48 saat sonra enfeksiyondan (adım 2.3) de paylaştırılır olmalıdır.
Sonikleştirme ile veya Kısıtlama Enzim Digest tarafından 3. Fragman Genomik DNA
Not: neredeyse sekansa bağımsız bir şekilde Sonikasyon fragmanları genomik DNA ve dolayısıyla parçalanma tercih edilen biçimi eklendiğinde, düşük beklenen geri kazanım oranı ile dizilenerek numuneler (örneğin, nispeten düşük MOl'de başlatılan bulaşmış bir hasta hücreleri ya da enfeksiyonları). Ayrıca, sonikasyon biri Parti PCR çiftleri ayırt edilmelerine olanakkritik aynı yerde özel entegrasyonları arasında vasküler entegre edilen sekansı, enfekte olmuş hastalar (aşağıdaki Adım 11) 39,52-54 bölgesindeki provirüs içeren hücrelerinin klonal genişlemesini ayırt.
Not: DNA, LM-PCR sırasında iç viral dizilerinin amplifikasyonu azaltmak için yukarı LTR hemen alt baş ayrılabilir olmalıdır. Yukan U5 dizisine 43 bp'lik alt yatmaktadır ve bu Msel oluşturulan DNA daha sonra ligasyon için uyumlu olan kısıtlama enzimi BgIII birçok HIV-1 suşlarının (Şekil 1B) ile çalışır sona erer. Sonikasyon ile DNA hazırlarken, iç-yarma restriksiyon enzimi bağlayıcı bağlanmasından sonra uygulanmalıdır (- E ve aşağıda 4.3 Adım Şekil 1C bakınız).
- sonikasyon için, 120 ul'lik nihai bir hacme kadar nükleaz içermeyen su içinde genomik DNA, 10 ug karıştırın. (500 bp ortalama mola boyutu parametreleri kullanarak aşağıdaki paragraf iki tur ses dalgalarına maruzmetre: görev döngüsü:% 5; yoğunluk: 3; patlama başına devir: 200; Zaman: 80 sn).
- PCR saflaştırma kiti kullanılarak sonikasyona tabi tutulmuş DNA saflaştırılır. DNA DNA sonu tamir kiti ile uçlarını tamir ve PCR saflaştırma kiti kullanılarak DNA arındırmak. Bir kuyruk Klenow ekzo kullanılarak DNA - enzimi ve bir PCR saflaştırma kiti kullanılarak bir kuyruklu DNA arındırmak. Kit kullanımı ek ayrıntılar için 51,52 bakın.
- kısıtlama endonükleaz sindirimi, bir yanı sıra, üreticiler ve 5'-ta çıkıntılar oluşturur (her biri 100 U) enzimler kokteyline tarafından sağlanan tampon maddesi ile 100 ul bir hacim içinde gece boyunca 37 ° C 'de, genomik DNA, 10 ug kesilmiş bu BglII uyumsuz enzim yukan virüs LTR aşağısına bölen. PCR saflaştırma kiti kullanılarak ertesi gün DNA saflaştınlır.
Not: restriksiyon enzimleri hiçbiri terminali içinde kesmek gerekir ~ LM-PCR protokolü tarafından yükseltilir viral DNA sonu 30 bp. Bu protokol, özellikle U5 güçlendirirHIV-1 DNA sonu.
4. Tavlama bağlayıcı oligonükleotidlerden ve Parçalanmış genomik DNA ligate
Not: Yukarıdaki bir DNA fragmanları (Bu protokolde kullanılan oligonükleotidlerin sekansları için Tablo 1 'e bakınız) ile uyumlu olan bir çıkıntı ihtiva eden bir asimetrik bağlayıcıyı hazırlayın. Msel-sindirilmiş DNA bağlayıcı uyumlu 5'-TA fazlalığı (Şekil 1) içermelidir ise bağlayıcı, uyumlu T-3 'çıkıntı içermelidir sonikasyona uğratılmış DNA ile birlikte kullanılmak üzere. kısa bir bağlayıcı şerit ek olarak ilgili DNA daha sonra doğru amplifikasyon reaksiyonlarını sınırlamak için, örneğin 3'-amin olmayan bir uzatılabilir kimyasal modifikasyon içermelidir.
Not: paralel olarak birden çok farklı entegrasyon Alanı kitaplıkları hazırlanmasında ve / veya aynı sıralama kaçak çoğullama özgü örnekler, örnek çapraz Kirlenmiş potansiyelinin sınırlandırılması için her bir örnek için özel bağlayıcıları kullanmak tavsiye edilirPCR sırasında tirme. Bu ek olarak yarı iç içe PCR (aşağıda tarif edilmiştir) içinde, her bir örnek için özel bağlayıcı primerlerin kullanımı anlamına gelir. Benzersiz bağlayıcı iplikler ve bağlayıcı primerler benzer bir genel% GC içeriğine ve ilgili çıkıntı pozisyonları muhafaza ederken, Tablo 1 'de listelenen bağlayıcı oligonükleotid dizileri karıştırıcı ile tasarlanabilir.
- 90 ° C'de ısıtılarak 10 mM Tris-HCI, pH 8,0-0,1 mM EDTA (her oligonükleotidin 10 uM son konsantrasyon) 35 ul, kısa ve uzun bağlayıcı lifler tavlanması yavaşça 1 ° 'lik adımlarla Oda sıcaklığına soğutulduktan ve dakika başına ° C.
- 1.5 uM bağlanmış bağlayıcı, 1 ug parçalanmış DNA ve 50 ul 800 u T4 DNA ligaz içeren genomik DNA, örnek başına en az dört paralel ligasyon reaksiyonları, hazırlayın. 12 ° C de bir gece Arter. Bir PCR saflaştırma kiti ile ertesi gün arındırın.
- sonikasyon ile hazırlanan örnekler için, bir sınırlamalara 100 U saflaştınldı ligasyon reaksiyonu sindirmeküreticinin göre alt (HIV-1 için, örneğin, BglII) yukan LTR'den böler gece boyunca önerdiği koşullar altında bu yon enzimi. PCR saflaştırma kiti kullanılarak DNA saflaştırılır.
Yarı iç içe PCR ile 5. yükseltin Viral LTR-Host Genomik DNA Junctions
Not: geri kazanılan ligasyon reaksiyonunun DNA konsantrasyonuna bağlı olarak, uygun bir kütüphane çeşitliliği için en az 4-8 paralel PCR'ler sağlamak için, her iki PCR tur için her bir örnek için hazırlanmalıdır. DNA şablonu konsantrasyonunun spektrofotometri ile nicelenen gerekmektedir. Bu protokol, PCR ilk ve ikinci tur iç içe LTR-spesifik primerler kullanılmıştır, ancak, aynı bağlayıcı-spesifik primer her iki mermi (Tablo 1) kullanılır. İkinci tur LTR-spesifik primer ve DNA kümelenme ve dizileme primeri bağlayıcı siteler için bağlayıcı-spesifik primer kodlamak adaptör dizileri. İç içe LTR-spesifik primer de, WH 6 nt endeksi dizisini kodlayanich aynı sıralama çalışmasında çoklayıcı kütüphaneleri için farklı primerler arasında değiştirilebilir.
- Tablo 2'de listelenen tüp başına bileşenler içeren birinci tur PCRlar hazırlayın.
Not: bağlayıcı-spesifik primer bağlayıcı, terminaline 3 '15-16 bp yukanya doğru yer alan bir son' 53 ° C,% 45 bir GC içeriğine ve onun 3 arasında bir erime sıcaklığına tamamlayıcılık 22 mt barındıran farklı bağlayıcı uzun şeritler (Tablo 1). İlk tur 27 nt LTR astar 59 ° C,% 48 bir GC muhtevasının bir erime sıcaklığına sahiptir ve 3 'ucu, HIV-1 U5 terminalinden 34 bp üst baş bulunmaktadır. HIV-1 LTR tamamlayıcı olan ikincil 26 nt LTR primer bölgesi, 60 ° C,% 50 bir GC içeriğine ve 3 'ucuna bir erime sıcaklığına sahip, viral U5 18 bp üst baş bulunduğu terminus. Oligonükleotid erime sıcaklığı ve GC içeriği bu parametreleri kullanıcıları eğer taklit gerektiğini tavsiye edilirdeğiştirilmiş diziler 21 (diğer retrovirüsler ile kullanım için olan) tasarım, PCR primerleri. - Aşağıdaki PCR parametreleri altında ilk PCR yuvarlak çalıştırın: Bir döngü: 2 dakika süreyle 94 ° C; 30 döngü: 15 saniye boyunca 94 ° C, 30 saniye, 45 saniye boyunca 68 ° C, 55 ° C; Bir çevrim: 10 dakika boyunca 68 ° C.
- reaksiyonları bir havuz ve PCR saflaştırma kiti kullanılarak saflaştırılması. . Tablo 3'de belirtilen tüp başına bileşenler içeren ikinci tur PCRlar hazırlayın Adım 5.2'de açıklanan PCR parametreleri kullanılarak PCR ikinci turda çalıştırın. reaksiyonları havuzu ve üreticinin talimatlarına göre bir ticari PCR saflaştırma kiti kullanılarak DNA arındırmak.
Not: DNA, kümelenme NGS uyumlu önerilen göstergesi çeşitli diziler 71 mevcuttur.
6. gerçekleştirin QC ve NGS (Tipik Ardışık Tesisi tarafından Tamamlandı)
- (QC tahlil # 1) Onayla Adım 5.3 kütüphane DNA konsantrasyonu floro kullanarakmetre 55. Kısaca, 200 ul nükleaz içermeyen su bir son hacim içinde standartları ve deneysel numunelerin hazırlanması. Vorteks tüpler 2-3 saniye, 2 dakika boyunca oda sıcaklığında inkübe ve florometre örnekleri okuyun.
Not: Örnekler 15 ul'lik bir minimum hacimde 2 nM kütüphanesi DNA minimum konsantrasyonu içermelidir. - (QC tahlil # 2) bir teyp tabanlı tahlil 56 kullanılarak DNA parçası boyutu dağılımını onaylayın.
Not: Bir ideal dağıtım uzunluğu yaklaşık 500 bp merkezleme nispeten geniş bir DNA zirvedir. malzemenin önemli bir miktarı, 1 kb daha büyükse, o zaman kümelenme sırasında köprü amplifikasyonu engelleyecektir uzun DNA türleri, ortadan kaldırmak için bir boyut seçim prosedürü bulunması tavsiye edilir. önemli bir pik yaklaşık 100 bp ila 200 belirgindir Bunun tersine, bir primer dimer PCR sırasında meydana gelmiş olabilirler. Bu durumda, prosedür, primer dimerlerin oluşumunu en aza indirmek için optimize edilmelidir. - (QC tahlil # 3) Cokantitatif PCR 57 DNA kütüphaneye adaptörlerin uygun birleşme teyid edildi.
- Üreticinin uygulama literatürde aşağıdaki NGS gerçekleştirin. % 10 başak-in yararlanın (a / a) sıralama vadede dengeli bir baz kompozisyonu sağlayarak gerçek zamanlı kalite ölçütleri optimize ΦX174 DNA.
Not: Entegrasyon sitesi sıralama deneyler genellikle tek uç 150 bp (SE150) ya da eşleştirilmiş-uç 150 bp (PE150) sıralama tabi tutulmaktadır. PE150 Her bir DNA molekülü üzerinde bağlayıcı bağlanma noktasını yakalamak için özellikle yararlıdır (örneğin, konakçı hücre klonal genişleme için kanıt entegrasyon sitesi incelenmesi sırasında).
7. LTR içeren Diziler, Bitki uzak LTR ve bağlayıcı Diziler için Sıralama Verilerini ayrıştırma özelleştirilmiş Python veya PERL komut dosyası kullanın ve Blat Referans Genom için Harita
- Sekans okur LTR içeren uzaklıkta konukçu genomik DNA dizisinden bitki LTR ve bağlayıcı dizileri ve tara FAŞTA dosyalarıYeni bir FAŞTA dosyasına bu dizileri ihracat. Harita referans genom hem okur kırpılmış (örneğin insan genomu versiyonları hg19 veya GRCh38) ve çıkış entegrasyon sitesi ile Blat 58 kullanılarak viral genomun, aşağıdaki ayarları kullanarak, ayrı bir .txt dosyası ihraç koordinatları:
stepSize = 6, minIdentity = 97 ve maxIntron = 0 - HIV-1 genomunun ve diğer diziler haritalama (LTR ucu dahili viral DNA genom bölgesi içine entegre olduğunu, yani bir kanıt) autointegrations kaldırmak dosya txt Blat çıkış ayrıştırmak ve dosya txt ayrı çıktı oluşturmak hangi tüm yinelenen entegrasyon siteleri tek ve benzersiz koordinat hit yoğunlaşmış edilmiştir.
8. FAŞTA Dosyalara Bu dönüştürme, Entegrasyonlar Etraftaki 15-Nt aralıkları İçeren .Yatakta Dosyaları oluşturun ve Taban Tercihleri Entegrasyon Siteleri Etraftaki Display Dizi Logolar Construct
- üsleri bir aralık için liste .Yatakta dosyaları oluşturunHer entegrasyon sitesi. En az 15 baz (5 yukarı ve aşağı 10) dizi logosu üretimi için önerilmektedir. BEDTools 59 fastaFromBed fonksiyonu ve bu komutu kullanarak bu .Yatakta dosyalarından bir FASTA dosyası oluşturun:
fastaFromBed -fi / dizin / için / referans / genom / -name -s -bed 15_base_pair_file.bed -fo output_file.fasta
Not: değişmeyen viral 5'-CA-3 'dinükleotid entegrasyonu sırasında DNA barındırmak için katıldı ve hücresel DNA LTR terminalinin kavşağı doğrulayarak iyi niyetli entegrasyon siteleri belirlemek için önemli bir başlangıç filtre edilir. Biz ayrıca deney sonuçlarını doğrulamak için bu konak DNA dizisi nüfus dizi logoları derlemek. Retrovirüsler kendi entegrasyon sitelerini 14,15 çevreleyen imza baz tercihlerini görüntülemek üzere, dizi logolar gibi homolog olmayan DNA gibi diğer rekombinasyon mekanizmaları ile karşılaştırıldığında IN-aracılı entegrasyonu ile eşleştirilmiş genomik siteleri ortaya çıktığını doğrulamak için hizmet60,61 katılmadan sona erer. - Kullanım WebLogo 3 (http://weblogo.threeplusone.com/create.cgi) FASTA dosyaları dizi logolar oluşturmak için. Click FAŞTA dosya yüklemek için 'Dosya Seç', ve aşağıdaki ayarları kullanın: Çıktı biçimi PDF (vektör); Logo boy, büyük; İlk pozisyon numarası, -5; Logo aralığı, -5 ila 5; Y ekseni ölçeği, 0.1, Y ekseni tic aralık, 0.5, Renk düzeni, klasik (NA).
9. Örnek Çapraz Kontaminasyon için kontrol edin, Merkez Bankası Çifti Dosyaları .bed oluşturun ve ihtiyaç duyulan Genomik Özellikleri Bağıl Benzersiz Entegrasyon Siteleri Dağılımı Harita
- retroviral entegrasyonu T-DNA şerit boyunca kademeli bir şekilde gerçekleştiğinden, genomik özelliklere genomik dağıtım nisbetle doğru eşleme için hedef site çoğaltılması merkez bp'lik yansıtmak için entegrasyon sitelerinden tam koordinatlarını ayarlayın.
- Bu nedenle, 5 bp HIV-1 gibi virüs çoğaltma için, i ofset merkezi bp ile .Yatakta dosya oluşturmakartı sarmalına entegrasyonlar haritalama için aşağı iki üsleri ve eksi strand entegrasyonlar haritalama için yukarı iki üsleri tarafından ntegration sitesi.
- Merkezi bp iki farklı örnek için ve bu komutu takip ederek dosyaları .bed kesişir BEDTools işlevini kesişir kullanarak farklı kütüphaneleri arasında ortak entegrasyon sitelerin sayısını hesaplamak, örnek çapraz kontaminasyon kontrol etmek için:
bedtools kesiştiği -a central_basepair_1.bed -b central_basepair_2.bed -f 1.00 -r -s> overlap1v2.txt - Aşağıdaki komutu kullanarak iki kütüphaneler arasında ortak sitelerin tam sayısını ölçmek için çıktı overlap1v2.txt dosyası içinde satır sayısını saymak:
wc -l overlap1v2.txt - Örneğin http://hgdownload.cse.ucsc.edu/goldenPath/hg38/da (UCSC Genom Açıklaması Veritabanı entegrasyonu sitesi haritalama için kullanılan referans genomunun sürümü için RefSeq açıklama .Yatakta dosyasını indirinkurumların veri tabanı) 62.
- RefSeq Bu komutun aşağıdaki dosya .bed ile numune için oluşturulan merkezi baz çifti .bed dosyayı kesiştiği için BEDTools işlevini kesiştiği kullanarak RefSeq genler giren entegrasyon sitelerin sayısını hesaplayın:
bedtools kesiştiği -a central_basepair_1.bed -b RefSeq_hg38.bed -u> RefSeq_sample1.bed
- aşağıdaki komutu kullanarak RefSeq genlerde düşen sitelerin tam sayısını ölçmek için çıktı RefSeq_sample1.bed dosyası içinde satır sayısını saymak:
wc -l RefSeq_sample1.bed - Yineleyin bir aralık dosyası mevcut .bed olan ilgi başka bir açıklama eşleme entegrasyon siteleri için 9.3 ve 9.4 adımları tekrarlayın. Adım 9.4 belirtildiği gibi UCSC Genom Açıklaması Database ilgi referans genomunun en güncel CpG ada açıklama .Yatakta dosyasını indirin.
- Belirli bir di giren entegrasyon sitelerin sayısını hesaplayınduruş BEDTools pencere işlevini kullanarak ve bu komutu takip ederek CpG adalarının (bu örnekte gösterildiği bir 5 kb penceresidir):
central_basepair_1.bed -b CpG_hg38.bed -u> CpG_sample1.bed 2500 -w bedtools penceresi
- aşağıdaki komutu kullanarak yukarı veya aşağı CpG adalarının 2.5 kb giren sitelerin tam sayısını ölçmek için çıktı CpG_sample1.bed dosyası içinde satır sayısını saymak:
wc -l CpG_sample1.bed - Tekrarlayın 9.6 ve haritalama entegrasyon siteleri için 9.7 yakındaki TSSler adımlar. genomik birden fazla gen haritalama koordinatları RefSeq.bed dosyasının, alternatif bir sürümünü oluşturmak bu konumda sadece tek bir gen hediye yansıtacak şekilde ayarlandı. Bu entegrasyon sitesi etrafındaki gen yoğunluğu abartılmış önler. BEDTools pencere işlevini kullanarak ve bu komutu takip ederek, her entegrasyon alanını çevreleyen 1 Mb bölgesinde gen yoğunluğunu hesaplamak:
Bedtools penceresi -w 500000 central_basepair_1.bed -b RefSeq_hg38_NonRedundant.bed -u> GeneDensity_sample1.bed - Bu komutu takip ederek veri kümesindeki tüm entegrasyonları için ortalama gen yoğunluğu hesaplayın:
awk '(toplamı + = $ 7) END (print "Ortalama =" toplamı / NR)' GeneDensity_sample1.bed
10. istatistiksel Ar Fisher Exact Testi ve İki kuyruklu Wilcoxon Rank Sum Testi İki kuyruklu kullanma Örnekleri arasında entegrasyon Sitesi Dağılımlar karşılaştır
RefSeq genler içinde veya CpG adaları ya da TSSlerin bir pencere içinde entegrasyon sitelerin oranını karşılaştırmak için kullanın Fisher exact testi, ancak entegrasyon siteleri çevreleyen gen yoğunluğu dağılımını karşılaştırmak için Wilcoxon rank sum testini kullanın: Not. R programı http://www.r-project.org/ mevcuttur.
Fisher testi iki uçlu:
- cr, Adımlar 9.4 ve 9.7 belirtildiği gibi hesaplanır numaralarını kullanarakgözlenen olaylar bu komutu izleyerek siteleri kalan karşı (bir açıklama içinde veya bir ek açıklama çevreleyen bir pencere içinde entegrasyonları) R her karşılaştırma için matrisler eate:
(Annotation_of_interest <- Matris (c (SampleA # içinde, SampleA # kalan, SampleB # kalan SampleB # içinde), sat_say = 2, dimnames = liste (c ( 'Merkez', 'Kalan'), c ( 'SampleA' 'SampleB')))) - Aşağıdaki komutu Fisher kesin testi iki kuyruklu tarafından karşılaştırma için p değerini hesaplayın:
fisher.test (annotation_of_interest alternatif = 'two.sided') $ p.value
Wilcoxon rank sum testi iki uçlu: - , Her sütun üst hücredeki örnek adını içeren bir sekme ile sınırlandırılmış .txt dosyası oluşturun (Adım 9.9 oluşturulan .Yatakta dosyasından alınan) o kütüphanedeki tüm entegrasyon siteleri için gen yoğunluk değerleri ile aşağıda izledi. Aşağıdaki komutu kullanarak Ar içine bu sekme ile sınırlandırılmış .txt dosyasını almak ve nadoğru dosya dizinine vigating:
FILENAME <- as.data.frame (read.delim (file.choose (), başlık = T, check.names = FALSE) = DOĞRU Eylül = ' t' doldurun) - Aşağıdaki komutu iki kuyruklu Wilcoxon rank sum testi ile karşılaştırma için p değerini hesaplayın:
wilcox.test (FILENAME $ eşleştirilmiş SampleA, DOSYAADI $ SampleB, alternatif = 'two.sided' = F, tam = T) $ p.value
Not: P değerleri program tarafından döndürülür sıfır sonra, sadece R belli bir (çok düşük) sınırın aşağı hesaplanabilir. R bir P = 0 verim kitlesel farklı numuneler için, <2.2 x 10 -308 olarak P değerini tahmin.
11. Entegre Viral DNA'yı içeren hücrelerin klonal genişleme için kanıt Ham Dizileme verileri inceleyin
Not: Küçük bir potansiyel referansı genomunda aynı nt birden fazla entegre ihtiyaç vardır. Alternatif olarak, tek birbütünleşme etkinliği nedeniyle kütüphane hazırlanması ve / veya DNA hazırlama öncesinde hücre çoğaltma PCR kullanımı dizileme veri artık olarak mevcut olabilir. HIV ile enfekte hastalardan genomik DNA'nın Son analizler aynı entegrasyon siteleri 52-54 ihtiva eden DNA dizilerine içinde (sadece öncesinde PCR ortaya çıkabilir) eşsiz sonication kesme noktaları / bağlayıcı bağlantı noktaları belirleyerek bu olanakları ayırt var. klonal genişletilmiş hücrelerin içinde barındırdığı proviruses latent viral rezervuar katkıda ve böylece insan hastalarda entegrasyon siteleri okuyan zaman genişleme düzeylerini karakterize özel ilgi olup olmadığı konusunda bir tartışma şu anda yok.
- Adım 8.1 listelenen prosedüre benzer şekilde, bu durumda uzanan üsleri bir aralık listeleyen .Yatakta dosyaları oluşturmak, 25 nt aşağısında her benzersiz entegrasyon sitesinden (yukarı bazlar burada gereksizdir). belirtildiği gibi (bu .Yatakta dosyalarından bir FASTA dosyası oluşturunBEDTools gelen fastaFromBed işlevini kullanarak ve bu komutu takip ederek adım 8.1):
fastaFromBed -fi / dizin / için / referans / genom / -name -s -bed 25_base_pair_file.bed -fo output_file.fasta
Not: klonal büyüme analizi için, her entegrasyon sitesinden nt aşağısında, en az 25 elde etmek için önerilmektedir arama, her özgüllüğünün geliştirilmesi için. - Tercihen özelleştirilmiş komut dosyası kullanarak, her benzersiz entegrasyon sitesinden nt aşağı 25 tam eşleme içeren tüm dizeleri için ham dizi veri FASTA dosyasını aramak, ve yeni bir dosya içine bu dizileri mevduat. ham dizeleri LTR ve bağlayıcı dizileri Trim. dizeleri nt örtüşen en az 20 paylaşıyorsa, ters tamamlayıcı okur kendi read1 çifti read2 dizeleri atayarak sonra LTR ve bağlayıcı dizileri kırparak ve dönüştürerek okur PE dizisi Birleştirme.
- Her entegrasyon sitesi bloğunun bağlayıcı bağlantı noktaları tarayın. klonal genişletilmiş "Her entegrasyonu sınıflandırır &# 34; Bağlayıcı bağlantı noktaları ise ≥3 ayrı bp.
Not: dizisi birleştirme olmadan klonal genişleme analizi için bir protokol 52 tarif edilmiştir okur.
Not: sonikasyon ile tam aynı yerde genomun Parçalanma klonal genişleme ölçüde küçümsenmesi yol açar ve yöntemler 63,64 tarif edilmiştir Elde edilen deneysel önyargı düzeltmek için.