Aonad 6 / 11

Anailís Sonraí Omics: RNA-seq, Slonn, Staitisticí agus Saibhriú Conair

Gnóthachain:

  • Cumas sreabhadh oibre RNA-seq, anailís ar léiriú difreálach agus ceartúchán ilthástála (FDR) a thuiscint agus cód anailíse infhíoraithe a tháirgeadh le hintleachta saorga
  • Cumas torthaí saibhrithe cosáin a léirmhíniú go criticiúil go staitistiúil agus go bitheolaíoch
  • Cumas an disciplín a fheidhmiú maidir le boinn tuisceana staitistiúla agus iliomad gaistí tástála a sheiceáil agus tábhacht bhitheolaíoch a fhíorú go neamhspleách.

Is é “Omics” an t-ainm comhchoiteann do chur chuige a thomhaiseann gach móilín de chill nó de fhíochán le chéile: géanóm (DNA ar fad), trascríobh (RNA/léiriú ar fad), próitéamaíocht (gach próitéin), meitibilít (gach meitibilít). Tá na sonraí seo ollmhór - baineann turgnamh RNA-seq le tomhais na mílte géinte. San aonad seo, foghlaimeoidh tú conas hintleachta saorga (AI) a úsáid mar chúntóir in anailís omics a scríobhann cód, a roghnaíonn staitisticí, agus a dhréachtaíonn léirmhíniú bitheolaíoch; ach beidh tú ag foghlaim cén fáth go gcaithfidh tú an toradh staidrimh agus bitheolaíoch a fhíorú.

Rabhadh criticiúil: In Omics, is iad na hearráidí staidrimh agus dofheicthe na hearráidí is contúirtí. Is féidir leis an AI cód a scríobh a sheachnaíonn ceartú comparáide iolrach (thíos), a roghnaíonn an tástáil mícheart, nó a tháirgeann liostaí géine “dearfach bréagach”. Ina theannta sin, is féidir le AI éilimh bhitheolaíocha a dhéanamh mar "méadaíonn an ghéin seo sa ghalar sin" gan aon fhoinse. An bealach ceart: an anailís a dhéanamh le cód inrite, in-iniúchta agus gach éileamh bitheolaíoch sa litríocht a dhearbhú.

Coincheapa bunúsacha

  • Slonn: Cé mhéad géine a aistrítear go RNA; tomhas ar "gníomhaíocht".
  • Slonn difreálach (DE): Géin a n-athraíonn a slonn go suntasach idir dhá ghrúpa (m.sh., othar/sláintiúil).
  • p-luach: An dóchúlacht gur comhtharlú í difríocht; má tá sé beag, meastar go bhfuil an difríocht "suntasach".
  • Ceartú comparáide iolracha: Nuair a bhreathnaítear ar na mílte géinte ag an am céanna, de sheans go mbeidh roinnt acu atá “suntasach”. Chun é seo a shocrú, úsáidtear modhanna mar FDR (ráta fionnachtana bréagach) / Benjamini-Hochberg. Má fhágtar an ceartú seo ar lár, tiocfaidh na céadta torthaí bréagacha chun cinn.
  • athrú fillte log2 (log2FC): Logarithm an chóimheas slonn de ghéin idir dhá ghrúpa go dtí an bonn 2; Ciallaíonn +1 méadú faoi dhó, −1 ciallaíonn laghdú faoi dhó.
  • Saibhriú cosáin: A fháil amach ina bhfuil conairí bitheolaíocha (m.sh. ceall-roinnt, díolúine) comhchruinnithe ar na géinte athraithe; Úsáidtear bunachair shonraí mar GO agus KEGG.
  • Éifeacht bhaisc: Difríocht bhréagach neamhbhitheolaíoch a eascraíonn as próiseáil samplaí ar laethanta/feistí éagsúla.

Céim ar chéim: anailís omics atá á gcumhachtú ag AI

1. An dearadh turgnamhach agus an cheist a shoiléiriú. Cé mhéad sampla, cé mhéad grúpa, cé mhéad athrá? An leor an chumhacht staidrimh? Mínigh an dearadh don AI.

2. Roghnaigh an uirlis/modh cuí le AI. Le haghaidh RNA-seq, roghnaigh modhanna caighdeánacha mar DESeq2/edgeR (pacáistí staidrimh deartha le haghaidh sonraí comhairimh); Bain úsáid as modhanna cruthaithe seachas staitistic a bhfuil an AI “déanta suas”.

3. Rith an cód agus seiceáil na haschuir idirmheánacha. Ná téigh ar aghaidh chuig anailís DE gan normalú, rialú éifeacht bhaisc, ceapacha cáilíochta (PCA).

4. Ceartú comparáide iolrach a fhorfheidhmiú. Scag na torthaí de réir luacha ceartaithe (padj/FDR), ní luach-p amh.

5. Deimhnigh an léirmhíniú bitheolaíoch sa litríocht. Léirmhínigh aschur saibhrithe cosáin le AI, ach fíoraigh gach éileamh ag an bhfoinse.

Leid: In anailís RO, féach ar na graif tionchair cháilíochta agus chomhiomlána ar dtús. Má dhéantar samplaí a bhraisliú faoin lá a próiseáladh iad seachas de réir grúpa bitheolaíoch, is éifeachtaí carnacha an chuid is mó de na géinte "suntasacha" a fhaigheann tú, ní fíor-bhitheolaíocht.

trí chás mion

Cás 1 — P-luach neamhcheartaithe. Thástáil mac léinn 20,000 géinte leis an gcód scríofa ag an AI agus fuair sé "540 géinte suntasacha." Nuair a scrúdaigh sé an cód, chonaic sé go raibh an AI tar éis ceartú comparáide iolracha a scipeáil. Nuair a cuireadh ceartú FDR leis, tháinig laghdú ar líon na ngéinte suntasacha go 32. Gan an ceartú, bheadh ​​​​níos mó ná géinte bréagach 500 bunaithe ar an scéal.

Cás 2 — Éileamh bitheolaíoch déanta. Maidir le géine ar an liosta DE, d'iarr taighdeoir ar an AI "cad a dhéanann an géine seo sa ghalar seo?" d'iarr sé; Mhínigh AI meicníocht áititheach agus an t-alt. Nuair a chuardaigh sé ar PubMed, chonaic sé nach raibh an mheicníocht sin ná an t-alt ann. Baineadh an t-éileamh den tuarascáil.

Cás 3 — Daingniú faighte. Thug mac léinn PhD faoi deara go raibh na samplaí scartha le dhá lá sa phlota PCA. D'iarr an AI athróg éifeacht bhaisc a chur leis an gcód; Tar éis an ceartúcháin, athraíodh an liosta géine go hiomlán agus tháinig sé suntasach go bitheolaíoch. Gan an chairt um rialú cáilíochta, d’fhéadfaí toradh bréige a fhoilsiú.

Sampla: scagadh le p-luach ceartaithe

pandaí a allmhairiú mar pd# ligean tábla 'de' a bheith mar thorthaí ó uirlis DE (DESeq2/edgeR): # colúin: géine, log2FC, pvalue, padj (FDR-corrected)de = pd.read_csv("de_results.csv") suntasach = de[(de["padj"] < 0.05s) & (de). 1)]print("Amh p<0.05:"", (de["pvalue"] < 0.05).sum())print("Padj ceartaithe FDR<0.05 & |log2FC|>=1:", len(suntasach))

Léiríonn an difríocht idir an uimhir amh agus uimhir cheartaithe an fáth a bhfuil sé ríthábhachtach comparáidí iolracha.

Ceithre theimpléad inchóipeála

1) Plean anailíse agus roghnú modhanna:

Do ról: cúntóir bithfhaisnéisíochta. Socraigh plean anailíse don turgnamh RNA-seq seo a leanas:[líon grúpaí, líon samplaí/macasamhlaithe, ceist]. Cén uirlis chaighdeánach (DESeq2/edgeR) ba cheart dom a roghnú agus cén fáth, cad iad na céimeanna rialaithe cáilíochta (PCA, éifeacht bhaisc) atá ag teastáil, conas a chuirim ceartúchán ilchomparáide i bhfeidhm? Scríobh céim ar chéim.

2) Cód + seiceálacha éigeantacha:

Scríobh cód inrite a dhéanann an anailís DE seo a leanas: [mion]. NÍ MÓR a áireamh sa chód normalú, plota cáilíochta PCA, rialú éifeacht bhaisc, agus ceartú FDR (Benjamini-Hochberg). Déan trácht ar gach céim. Scag le p-luach ceartaithe, ní le p-luach amh.

3) Trácht saibhrithe cosán:

Cabhraigh le torthaí saibhrithe cosáin a léirmhíniú don liosta seo de ghéinte suntasacha: [liosta/aschur]. Mínigh cé na conairí atá feiceálach, ach inis dom cén fhoinse (GO, KEGG, alt piarmheasúnaithe) chun gach éileamh bitheolaíoch a dhearbhú. Meicníocht/Feistiú alt.

4) Iniúchadh staitisticí:

Seiceáil an cód anailíse seo a leanas le haghaidh earráidí staidrimh: [cód]. Go sonrach: roghnú mícheart tástála, ceartú comparáide iolracha a fhágáil ar lár, neamhaird a dhéanamh ar éifeacht bhaisc, macasamhlú neamhleor. Déan liosta de gach fadhb a fuair tú agus a réiteach.

Leid lag / Pras láidir

Lag: "Faigh géinte suntasacha sna sonraí RNA-seq seo."

Fadhb: Ní shonraítear modh, rialú cáilíochta, agus comparáidí iolracha; Is féidir le AI liosta iomlán de rudaí dearfacha bréagacha a thabhairt gan aon cheartú.

Láidir: "Scríobh cód a dhéanann anailís DE ar na sonraí comhairimh RNA-seq seo le loighic DESeq2, lena n-áirítear rialú cáilíochta agus rialú bulcéifeacht le PCA, agus scagairí le ceartú FDR; déan trácht ar gach céim agus mínigh cén fáth ar roghnaigh tú an modh seo."

Cén fáth go bhfuil sé cumhachtach: Tá an modh caighdeánach, tá cáilíocht agus ceartú éigeantach, is féidir an toradh a iniúchadh.

Riosca

symptom

réamhchúram

bréagach dearfach

An iomarca géinte "brí".

FDR/ceartúchán comparáide iolrach

tionchar comhchoiteann

Cnuasaítear samplaí de réir an lae

Athróg bhaisc PCA +

tástáil mícheart

Gnáth-thástáil chun sonraí a chomhaireamh

Modh cuí amhail DESeq2/edgeR

déanta suas bitheolaíocht

Meicníocht gan táthú

Daingniú litríochta

cumhacht neamhleor

1-2 ionadaí

Go leor athrá sa dearadh

Botúin choitianta

  • Gan bacadh le ceartú comparáide iolrach. An earráid staidrimh is coitianta agus is díobhálaí.
  • Neamhaird a thabhairt ar an tionchar comhchoiteann. Táirgeann sé difríocht bitheolaíoch bréagach.
  • Ag cur tástáil mhícheart i bhfeidhm chun sonraí a chomhaireamh. Éilíonn RNA-seq modhanna speisialta.
  • Glacadh le héilimh bhitheolaíocha gan foinse. Is féidir le AI meicníochtaí agus earraí a dhéanamh.
  • Ginearálú gan dóthain athrá. Gan cumhacht staitistiúil, tá an toradh neamhiontaofa.
Aire: Ní hionann “suntasach ó thaobh staitistice” agus “suntasach go bitheolaíoch”. D'fhéadfadh athrú fillte an-bheag ach suntasach go teicniúil a bheith neamhshuntasach ó thaobh na bitheolaíochta de; I samplaí móra is féidir le gach rud a bheith "suntasach". Déan log2FC agus p-value a mheas le chéile.

Doimhneacht: cúlra agus gaistí comhaireamh dúbailte i saibhriú cosáin

Braitheann torthaí saibhrithe cosáin ar dhá thoimhde folaithe nach dtuigeann an chuid is mó daoine, agus is féidir le AI iad a sheachbhóthar go ciúin. Is é an chéad rogha cúlra/cruinne: déanann saibhriú an tacar “géinte a d’athraigh” a chur i gcomparáid leis an tacar “ar breathnaíodh ar ghéinte”. Má bhaintear an cúlra ón ngéanóm ar fad ach nach ndéanann do thurgnamh ach painéal fíocháin ar leith a thomhas, feictear go bhfuil na torthaí “saibhrithe” go saorga. Is éard atá sa chúlra ceart ná géinte ar féidir iad a chur in iúl/a thomhas sa turgnamh. Fuair ​​foireann amháin "saibhriú an-suntasach ar an gcosán imdhíonachta" trí chúlra an ghéanóm ar fad a chur ina chúlra; Nuair a rinneadh an anailís arís agus arís eile leis an gcúlra ceart (géinte tomhaiste), d'imigh an saibhriú - déantán modha a bhí sa toradh.

Ar an dara dul síos, méid tacair ghéinte agus comhaireamh dúbailte: feictear go bhfuil cosáin an-mhóra agus ginearálta (m.sh. “próisis meitibileach”, na mílte géinte) “suntasach” i mbeagnach gach liosta; tá bealaí beaga sonracha níos faisnéiseach. Ina theannta sin, toisc go bhfuil an ghéin chéanna le fáil ar ilbhealaí, tá sé míthreorach bealaí forluiteacha a chóireáil mar fhianaise neamhspleách. An tríú pointe: ní léiríonn sé treo an saibhrithe; Is féidir cosán a shaibhriú, ach féadfar leath de na géinte laistigh di a mhéadú agus an leath eile a laghdú. Chun é seo a fheiceáil, is gá faisnéis treorach a scrúdú go leithleach (amhail GSEA).

gaiste

symptom

réamhchúram

cúlra mícheart

Is cosúil go bhfuil gach rud saibhrithe

Faigh cúlra géinte tomhaiste

Cosán an-ghinearálta

Tagann "meitibileacht" suas i gcónaí

Fócas ar bhealaí beaga sonracha

comhaireamh dúbailte

cosáin forluiteacha

Ná glac leis mar fhianaise neamhspleách

skip treo

measctha ag dul suas/ag dul síos

Rialú treorach le GSEA

Go hachomair

  • AI in anailís omics; is cúntóir é a roghnaíonn modhanna, a scríobhann cód, agus a dhréachtaíonn tuairimí; ní mór údar a bheith le cinntí maidir le staidreamh agus bitheolaíocht.
  • Ba cheart modhanna caighdeánacha cruthaithe (DESeq2/edgeR) a úsáid; Níor cheart gan bacadh le rialú cáilíochta agus iniúchadh comhthionchair.
  • Tá ceartú comparáide iolrach (FDR) éigeantach; déantar na torthaí a scagadh leis an luach ceartaithe.
  • Ní mór gach éileamh bitheolaíoch a dheimhniú sa litríocht; ba cheart "suntasach" a idirdhealú ó "tábhachtach".

Tasc feidhmchláir

Tóg tábla torthaí samplach DE (nó tacar sonraí RNA-seq oscailte). Déan comparáid idir comhairimh ghéinte suntasacha bunaithe ar amhluach-p agus tairseacha ceartaithe padj leis an mblúire thuas. Déan trácht ar an difríocht in abairt amháin. Ansin iarr léirmhíniú bitheolaíoch le teimpléad 3 le haghaidh géine faoi thrácht agus seiceáil an t-éileamh tú féin i PubMed.

seicliosta

  • [ ] Mheas mé an dearadh turgnamhach agus an chumhacht staidrimh.
  • [ ] Roghnaigh mé modh caighdeánach áisiúil.
  • [ ] Rinne mé PCA agus rialú cáilíochta éifeacht bhaisc.
  • [ ] Chuir mé ceartúchán iolrachchomparáide (FDR) i bhfeidhm.
  • [ ] Scag mé na torthaí le p-luach ceartaithe.
  • [ ] Dheimhnigh mé na héilimh bhitheolaíocha sa litríocht.