FASTA Input to Generate SCENIC+ Databases for Rat

[jonhsussman]

Hello,

Thanks for the helpful tutorial. I am trying to generate the relevant database files to use SCENIC+ for rat genome. I've generated the motifs table by using the motifs-v10-nr.mgi-m0.00001-o0.0.tbl file and replacing the gene_name with the homologous gene for rat. But it is a little unclear what to use as the other input files for create_cistarget_motif_databases.py. For the motif IDs I am using all the motifs available under "singletons." But I am not sure what is the most correct to use for the fasta filename. Currently, I am trying this with the whole genome (Rattus_norvegicus.mRatBN7.2.dna.toplevel.fa) but it seems that it would make sense for these regions to be filtered in some way.

Is there an example of what was used to generate the databases for mouse/human or other species?

Thanks,
Jonathan

[ghuls]

You can use a BED file with pseudobulk peaks from your scATAC data and use bedtools to make the FASTA file:

bedtools getfasta \
    -fi Rattus_norvegicus.mRatBN7.2.dna.toplevel.fa \
    -bed pseudobulk_peaks_from_scATAC.bed \
    -fo rattus_norvegicus.pseudobulk_peaks.fa

[jonhsussman]

Thanks for this explanation, this is very helpful!

…

On Wed, Aug 16, 2023 at 6:54 AM Gert Hulselmans ***@***.***> wrote: You can use a BED file with pseudobulk peaks from your scATAC data and use bedtools to make the FASTA file: bedtools getfasta \ -fi Rattus_norvegicus.mRatBN7.2.dna.toplevel.fa \ -bed pseudobulk_peaks_from_scATAC.bed \ -fo rattus_norvegicus.pseudobulk_peaks.fa — Reply to this email directly, view it on GitHub <https://urldefense.com/v3/__https://github.com/aertslab/create_cisTarget_databases/issues/41*issuecomment-1680387133__;Iw!!LIr3w8kk_Xxm!uKOkHRSWEZpQL9tixRkixN6cBEdHQeVN4Mk_Plxu574NBMPhwFPPmdu-BIPpMphcEtrS2QtzaneWWvQMvOeN2cJs$>, or unsubscribe <https://urldefense.com/v3/__https://github.com/notifications/unsubscribe-auth/AKIGYMQ33K66H7HH2NGDBA3XVSRHNANCNFSM6AAAAAA3O5GMFU__;!!LIr3w8kk_Xxm!uKOkHRSWEZpQL9tixRkixN6cBEdHQeVN4Mk_Plxu574NBMPhwFPPmdu-BIPpMphcEtrS2QtzaneWWvQMvKRj8Uwf$> . You are receiving this because you authored the thread.Message ID: ***@***.***>