search

GoekeLab / xpore

Identification of differential RNA modifications from nanopore direct RNA sequencing
https://xpore.readthedocs.io/
MIT License
137 stars 22 forks source link

RNA004 Inverted DRACH Results #225

Open BrendanBeahan opened 3 months ago

BrendanBeahan commented 3 months ago

My team recently ran xPore on both an RNA002 and RNA004 dataset. The RNA004 dataset, understandably, resulted in some pretty strange outcomes. In particular, the RNA004 output seemed to have an inverse relationship between the statistical significance of the predicted modification and the amount of DRACH enrichment said site was associated with. This is the opposite of what we observed in RNA002.

I suppose I'm just curious if your team had any insights into how you underlying method may lead to such an outcome or if this seems like a plausible outcome?

yuukiiwa commented 2 months ago

Hi @BrendanBeahan,

The main branch of xpore hasn't been updated to support RNA004. You can provide the RNA004 kmer model (https://github.com/GoekeLab/xpore/blob/RNA004_kmer_model/xpore/diffmod/RNA004_5mer_model.txt) during in the xpore diffmod config.yml when running xpore diffmod:

data:
    <CONDITION_NAME_1>:
        <REP1>: <DIR_PATH_TO_DATA_JSON>
        ...

    <CONDITION_NAME_2>:
        <REP1>: <DIR_PATH_TO_DATA_JSON>
        ...

    ...

model_kmer: <PATH_TO_KMER_MODEL>

out: <DIR_PATH_FOR_OUTPUTS>

Thanks!

Best wishes, Yuk Kei