Bölüm 6 / 11

Omics maglumatlary derňewi: RNA-seq, aňlatma, statistika we ýoly baýlaşdyrmak

Girdeji:

  • RNA-seq iş prosesine, diferensial aňlatma derňewine we birnäçe synag düzedişine (FDR) düşünmek we emeli intellekt barlanylýan derňew kody döredýär
  • Wayollary baýlaşdyrmagyň netijelerini statistik we biologiki taýdan tankytlamak ukyby
  • Statistik çaklamalary we birnäçe synag ýalňyşlyklaryny barlamak we biologiki ähmiýetini özbaşdak barlamak düzgünini ýerine ýetirmek ukyby.

“Omics” - öýjügiň ýa-da dokumanyň ähli molekulalaryny bilelikde ölçän çemeleşmeleriň umumy ady: genomika (ähli DNK), transkriptomika (ähli RNK / aňlatma), proteomika (ähli beloklar), metabolomika (ähli metabolitler). Bu maglumatlar ullakan - RNA-seq synagynda on müňlerçe geniň ölçegi bar. Bu bölümde, kod ýazýan, statistika saýlaýan we biologiki düşündiriş taýýarlaýan omics derňewinde kömekçi hökmünde emeli intellekti (AI) ulanmagy öwrenersiňiz; ýöne näme üçin statistiki we biologiki netijäni barlamalydygyňyzy öwrenersiňiz.

Tankytly duýduryş: Omiksde iň howply ýalňyşlyklar statistik we göze görünmeýän. AI birnäçe deňeşdirme düzedişinden (aşakda) geçýän, nädogry synagy saýlaýan ýa-da “ýalan polo positiveitel” gen sanawlaryny döredýän kod ýazyp biler. Mundan başga-da, AI hiç hili çeşmesiz "bu kesel şol gen köpelýär" ýaly biologiki talaplary öňe sürüp biler. Dogry ýol: ýerine ýetirip boljak, eşidilýän kod bilen derňew geçirmek we edebiýatda her biologiki talapy tassyklamak.

Esasy düşünjeler

  • Düşündiriş: RNK-a näçe gen terjime edilýär; "işjeňlik" çäresi.
  • Differensial aňlatma (DE): Iki toparyň arasynda (meselem, hassa / sagdyn) aňlatmasy ep-esli üýtgeýän genler.
  • p-baha: Tapawudyň tötänlik bolmagy ähtimallygy; kiçi bolsa, tapawut "möhüm" hasaplanýar.
  • Birnäçe deňeşdirme düzedişi: Şol bir wagtyň özünde on müňlerçe gen seredilende, tötänleýin "ähmiýetli" käbirleri bolar. Muny düzetmek üçin FDR (ýalňyş açyş tizligi) / Benjamini-Hochberg ýaly usullar ulanylýar. Bu düzediş goýulmasa, ýüzlerçe ýalan netijeler ýüze çykar.
  • log2 esse üýtgeşme (log2FC): Iki toparyň arasyndaky geniň aňlatma gatnaşygy logarifmi; +1 iki esse ýokarlanmagy aňladýar, −1 iki esse azalmagy aňladýar.
  • Wayollary baýlaşdyrmak: Üýtgedilen genleriň haýsy biologiki ýollaryň (meselem, öýjük bölünişi, immunitet) jemlenendigini tapmak; GO we KEGG ýaly maglumat bazalary ulanylýar.
  • Toparyň täsiri: Dürli günlerde / enjamlarda nusgalary gaýtadan işlemekden döreýän biologiki däl ýalňyş tapawut.

Stepdimme-ädim: AI bilen işleýän omikleriň derňewi

1. Synag dizaýnyny we soragyny düşündiriň. Näçe nusga, näçe topar, näçe gaýtalama? Statistik güýç ýeterlikmi? Dizaýny AI-e düşündiriň.

2. AI bilen degişli gural / usuly saýlaň. RNA-seq üçin DESeq2 / edgeR ýaly standart usullary saýlaň (maglumatlary sanamak üçin döredilen statistiki paketler); AI “düzen” statistika däl-de, subut edilen usullary ulanyň.

3. Kody işlediň we aralyk netijeleri barlaň. Adatylaşdyrma, partiýa täsirine gözegçilik, hil meýilnamalary (PCA) bolmazdan DE derňewine geçmäň.

4. Birnäçe deňeşdirme düzedişini ýerine ýetiriň. Netijeleri çig p-bahasy däl-de, düzedilen baha (padj / FDR) boýunça süzüň.

5. Edebiýatda biologiki düşündirişi tassyklaň. AI bilen ýoly baýlaşdyrmagyň netijesini düşündiriň, ýöne her talapy çeşmede barlaň.

Maslahat: DE derňewinde ilki hiline we jemleýji täsir grafalaryna serediň. Nusgalar biologiki topar däl-de, gaýtadan işlenen güni jemlenen bolsa, tapýan "möhüm" genleriň köpüsi hakyky biologiýa däl-de, jemleýji täsirlerdir.

üç sany kiçi halat

1-nji hadysa - düzedilmedik p-bahasy. Bir okuwçy AI tarapyndan ýazylan kod bilen 20,000 geni barlady we "540 möhüm gen" tapdy. Kody gözden geçirende, AI-iň birnäçe deňeşdirme düzedişinden geçendigini gördi. FDR düzedişi goşulanda möhüm genleriň sany 32-e çenli azaldy. Düzediş bolmasa 500-den gowrak ýalan gen hekaýa esaslanar.

2-nji hadysa - Biologiki talap. DE sanawyndaky gen üçin bir gözlegçi AI-den "bu gen bu keselde näme edýär?" Diýip sorady. diýip sorady. AI ynandyryjy mehanizmi we makalany düşündirdi. PubMed-de gözläninde ne mehanizmiň, ne-de makalanyň bardygyny gördi. Talap hasabatdan aýryldy.

3-nji waka - Tassyklama gazanyldy. Filosofiýanyň doktory okuwçylaryň nusgalaryň PCA meýilnamasynda iki gün bölünendigini gördi. AI-den koda partiýa effekt üýtgeýjisini goşmagyny haýyş etdi; Düzedişden soň gen sanawy düýbünden üýtgedildi we biologiki taýdan möhüm boldy. Hil gözegçiligi diagrammasy bolmasa galp netije çap edilip bilnerdi.

Mysal: düzedilen p-bahasy bilen süzmek

pandalary pd görnüşinde import ediň # tablisa 'de' DE guralynyň (DESeq2 / edgeR) netijesi bolsun: # sütünler: gen, log2FC, baha, padj (FDR düzedilen) de = pd.read_csv ("de_results.csv") ähmiýetli = de [(de ["padj"] <0.05) & (de) " p <0.05: ", (de [" baha "] <0.05) .sum ()) çap (" FDR düzedilen padj <0.05 & | log2FC |> = 1: ", len (ähmiýetli))

Çig we düzedilen sanyň arasyndaky tapawut, köp deňeşdirmeleriň näme üçin möhümdigini görkezýär.

Göçürip boljak dört şablon

1) Analiz meýilnamasy we usuly saýlamak:

Siziň roluňyz: bioinformatika kömekçisi. Aşakdaky RNA-seq synagy üçin derňew meýilnamasyny düzüň: [toparlaryň sany, nusgalaryň sany / göçürmeler, sorag]. Haýsy standart gural (DESeq2 / edgeR) saýlamaly we näme üçin, haýsy hil gözegçilik ädimleri (PCA, partiýa effekti) talap edilýär, köp deňeşdirme düzedişini nädip ulanmaly? Stepdim-ädim ýazyň.

2) Kod + hökmany barlaglar:

Aşakdaky DE derňewini ýerine ýetirýän ýerine ýetirip boljak kod ýazyň: [jikme-jiklik]. KOD kadalaşdyrmak, PCA hil sýu ploteti, partiýa täsirine gözegçilik we FDR (Benjamini-Hochberg) düzedişini öz içine alýar. Her ädimi düşündiriň. Çig p-bahasy däl-de, düzedilen p-bahasy bilen süzüň.

3) wayol baýlaşdyryş düşündirişi:

Möhüm genleriň sanawy üçin ýoly baýlaşdyrmagyň netijelerini düşündirmäge kömek ediň: [list / çykyş]. Haýsy ýollaryň görnüklidigini düşündiriň, ýöne her biologiki talapy tassyklamak üçin haýsy çeşmeden (GO, KEGG, deň-gözden geçirilen makala) aýdyň. Mehanizm / makala FITTING.

4) Statistika barlagy:

Statistik ýalňyşlyklar üçin aşakdaky derňew koduny barlaň: [kod]. Hususan-da: nädogry synag saýlamasy, köp deňeşdirme düzedişiniň ýalňyşlygy, partiýa täsirini äsgermezlik etmek, ýeterlik köpeltmek. Tapan her meseläňizi we düzedişini sanaň.

Gowşak haýyş / Güýçli haýyş

Gowşak: "Bu RNA-seq maglumatlarynda möhüm genleri tapyň."

Mesele: Usul, hil gözegçiligi we birnäçe deňeşdirme görkezilmedi; AI düzedişsiz ýalan pozitiwlerden doly sanawy berip biler.

Güýçli: "DESeq2 logikasy bilen bu RNA-seq sanamak maglumatlary üçin DE derňewini ýerine ýetirýän, PCA bilen hil gözegçiligini we köpçülikleýin täsir gözegçiligini we FDR düzedişli süzgüçleri öz içine alýan kod ýazyň; her ädimiňize düşündiriş beriň we näme üçin bu usuly saýlandygyňyzy düşündiriň."

Näme üçin güýçli: Usul standart, hil we düzediş hökmany, netijesi diňlenýär.

Töwekgelçilik

alamaty

seresaplyk

ýalan polo positiveitel

Örän "manyly" genler

FDR / köp deňeşdirme düzedişi

köpçülikleýin täsir

Nusgalar gün-günden jemlenýär

PCA + partiýa üýtgeýjisi

nädogry synag

Maglumatlary sanamak üçin adaty synag

DESeq2 / edgeR ýaly degişli usul

biologiýany düzdi

Kebşirleýiş mehanizmi

Edebiýaty tassyklamak

güýji ýeterlik däl

1-2 rep

Dizaýnda ýeterlik gaýtalanma

Adaty ýalňyşlyklar

  • Birnäçe deňeşdirme düzedişinden geçmek. Iň ýaýran we iň zyýanly statistiki ýalňyşlyk.
  • Köpçülikleýin täsire ähmiýet bermezlik. Falsealan biologiki tapawut döredýär.
  • Maglumatlary sanamak üçin nädogry synaglary ulanmak. RNA-seq ýörite usullary talap edýär.
  • Biologiki talaplary çeşmesiz kabul etmek. AI mehanizmleri we makalalary düzüp biler.
  • .Eterlik däl gaýtalamak bilen umumylaşdyrmak. Statistik güýç bolmasa, netije ygtybarly däl.
Üns beriň: “Statistik taýdan ähmiýetli” “biologiki taýdan ähmiýetli” däl. Örän kiçi, ýöne tehniki taýdan ähmiýetli ep-esli üýtgeşiklik biologiki taýdan ähmiýetsiz bolup biler; Uly nusgalarda hemme zat "ähmiýetli" bolup biler. Log2FC we p-bahany bilelikde bahalandyryň.

Çuňlugy: ýoly baýlaşdyrmakda fon we iki gezek sanamak

Wayol baýlaşdyryş netijeleri, köpleriň düşünmeýän iki sany gizlin çaklamasyna daýanýar we AI olary ýuwaşlyk bilen aýlap bilýär. Birinjisi fon / älem saýlamak: baýlaşdyrmak "üýtgän genler" toplumyny "haýsy genlere seredilendigi" bilen deňeşdirýär. Fon tutuş genomdan alnan bolsa, ýöne synagyňyz diňe belli bir dokuma panelini ölçese, netijeler emeli usulda “baýlaşdyrylan” bolup görünýär. Dogry fon, synagda hakykatdanam beýan edilip / ölçelip bilýän genlerdir. Bir topar tutuş genomy ýalňyşlyk bilen fon arkaly “immun ýoluny gaty ähmiýetli baýlaşdyrmagy” tapdy; Derňew dogry fon (ölçenen genler) bilen gaýtalananda, baýlaşdyrmak ýitdi - tapyş usuly artifaktdy.

Ikinjiden, gen toplumynyň ululygy we goşa sanamak: örän uly we umumy ýollar (meselem, “metabolik prosesler”, müňlerçe genler) her sanawda diýen ýaly “möhüm” bolup görünýär; kiçi, aýratyn ýollar has habarlydyr. Mundan başga-da, şol bir gen birnäçe ýolda tapylýandygy sebäpli, biri-birine gabat gelýän ýollara garaşsyz subutnama hökmünde garamak ýalňyşdyr. Üçünji nokat: baýlaşdyrmagyň ugruny görkezmeýär; Bir ýol baýlaşdyrylyp bilner, ýöne içindäki genleriň ýarysy köpelip, galan ýarysy azalmagy mümkin. Muny görmek üçin ugrukdyryjy maglumatlary (GSEA ýaly) aýratyn gözden geçirmeli.

duzak

alamaty

seresaplyk

nädogry fon

Hemme zat baýlaşdyrylan ýaly

Ölçelenen genleriň fonuny alyň

Örän umumy ýol

"Metabolizm" hemişe ýüze çykýar

Kiçijik, aýratyn ýollara üns beriň

goşa sanamak

gabat gelýän ýollar

Garaşsyz subutnama hökmünde kabul etmäň

ugry geçmek

garyşyk ýokarlanmak / aşak düşmek

GSEA bilen ugurlaýyn gözegçilik

Gysgaça

  • Omiki derňewde AI; usullary saýlaýan, kod ýazýan we teswir ýazýan kömekçi; statistika we biologiýa kararlary dogry bolmaly.
  • Adaty, subut edilen usullar (DESeq2 / edgeR) ulanylmaly; Hil gözegçiligi we köpçülikleýin täsir barlagy geçmeli däldir.
  • Birnäçe deňeşdirme düzedişi (FDR) hökmanydyr; netijeler düzedilen baha bilen süzülýär.
  • Her biologiki talap edebiýatda tassyklanmalydyr; "ähmiýetli" "möhüm" -den tapawutlanmalydyr.

Programma meselesi

DE netijeler tablisasyny alyň (ýa-da açyk RNA-seq maglumat toplumy). Çig p-bahasyna we düzedilen padj çäklerine esaslanýan möhüm gen sanlaryny ýokardaky bölek bilen deňeşdiriň. Bir sözlemdäki tapawudy düşündiriň. Soňra aýratyn gen üçin 3-nji şablon bilen biologiki düşündiriş soraň we PubMed-de özüňizi barlaň.

gözegçilik sanawy

  • [] Synag dizaýnyna we statistiki güýjüne baha berdim.
  • [] Adaty, amatly usuly saýladym.
  • [] PCA we partiýa effektiniň hiline gözegçilik etdim.
  • [] Birnäçe deňeşdirme (FDR) düzediş girizdim.
  • [] Netijeleri düzedilen p-baha bilen süzdüm.
  • [] Edebiýatda biologiki talaplary tassykladym.