sing-group / bicycle

bicycle (bisulfite-based methylcytosine caller) is a next-generation sequencing bioinformatics pipeline able to perform a full DNA methylation level analysis
http://www.sing-group.org/bicycle
GNU Lesser General Public License v3.0
8 stars 1 forks source link

Differential methylation for each annotated region (Bed file) header interpretation #8

Open freedomq8 opened 4 years ago

freedomq8 commented 4 years ago

Hi there,

thanks for your efforts to make my analysis easy along with lengthy explanation page. I am looking at the final differential methylation results using BED file and I have a question regarding headers and what do they mean:

#region sequence start stop SRR2052492 (treatment) SRR2052493 (treatment) SRR2052494 (treatment) SRR2052495 (treatment) SRR2052496 (treatment) SRR2052487 (control) SRR2052488 (control) SRR2052489 (control) SRR2052490 (control) SRR2052491 (control) treatment average control average log2FC(treament/control) p-value q-value MEG3 chr14 101291953 101292257 1455768/2043201 1636631/3258147 1497615/2643988 1278819/2677687 1699837/1958100 960649/2260743 1655354/3623531 1234101/2867390 1729232/3490690 2112253/4771083 0.6015893811705044 0.45208907524094044 0.4121720465419593 0.032485239405786 0.054142065676310004 CDKN1A chr6 36645463 36645696 106164/169351 67574/108667 39848/111440 70332/114278 56000/82614 184466/257701 135047/184961 153820/212228 90107/139840 97084/153940 0.5797185981069327 0.6962631895179566 -0.264279982267956 0.06970890125229327 0.0871361265653666 IRS1 chr2 227659612 227659781 471/1525 493/1510 556/2015 240/861 534/2047 326/818 325/1045 455/1336 680/1600 607/2070 0.2882633827594873 0.34837676517688165 -0.27326082007437236 0.031291095111098274 0.054142065676310004 INS chr11 2182552 2182775 133094/209138 291969/567402 119026/425194 152310/651761 177366/278973 80460/481158 165062/765934 174273/633564 169227/724203 288122/1036615 0.40974354597583645 0.2408760847942344 0.7664300613080263 0.022056895484407948 0.054142065676310004 PDE7B chr6 136172766 136172917 358560/2258785 151970/1578967 368147/2086541 104049/1062921 506169/2399992 119657/960632 197487/1205221 168450/1002387 197266/1376956 112174/1444852 0.15860896202767896 0.1327258145510687 0.2570252958195374 0.6929437289294518 0.6929437289294518

For example (i am using the results output from the manual section) i am not sure if i understand what is written beneath each sample at a certain target region

Thanks

osvaldogc commented 4 years ago

Hello !

I am answering your questions between your lines below.

osvaldo


From: freedomq8 [notifications@github.com] Sent: Monday, August 10, 2020 5:05 PM To: sing-group/bicycle Cc: Subscribed Subject: [sing-group/bicycle] Differential methylation for each annotated regionBed file) header interpretation (#8)

CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe

Hi there,

thanks for your efforts to make my analysis easy along with lengthy explanation page. I am looking at the final differential methylation results using BED file and I have a question regarding headers and what do they mean:

region sequence start stop SRR2052492 (treatment) SRR2052493 (treatment) SRR2052494 (treatment) SRR2052495 (treatment) SRR2052496 (treatment) SRR2052487 (control) SRR2052488 (control) SRR2052489 (control) SRR2052490 (control) SRR2052491 (control) treatment average control average log2FC(treament/control) p-value q-value

MEG3 chr14 101291953 101292257 1455768/2043201 1636631/3258147 1497615/2643988 1278819/2677687 1699837/1958100 960649/2260743 1655354/3623531 1234101/2867390 1729232/3490690 2112253/4771083 0.6015893811705044 0.45208907524094044 0.4121720465419593 0.032485239405786 0.054142065676310004 CDKN1A chr6 36645463 36645696 106164/169351 67574/108667 39848/111440 70332/114278 56000/82614 184466/257701 135047/184961 153820/212228 90107/139840 97084/153940 0.5797185981069327 0.6962631895179566 -0.264279982267956 0.06970890125229327 0.0871361265653666 IRS1 chr2 227659612 227659781 471/1525 493/1510 556/2015 240/861 534/2047 326/818 325/1045 455/1336 680/1600 607/2070 0.2882633827594873 0.34837676517688165 -0.27326082007437236 0.031291095111098274 0.054142065676310004 INS chr11 2182552 2182775 133094/209138 291969/567402 119026/425194 152310/651761 177366/278973 80460/481158 165062/765934 174273/633564 169227/724203 288122/1036615 0.40974354597583645 0.2408760847942344 0.7664300613080263 0.022056895484407948 0.054142065676310004 PDE7B chr6 136172766 136172917 358560/2258785 151970/1578967 368147/2086541 104049/1062921 506169/2399992 119657/960632 197487/1205221 168450/1002387 197266/1376956 112174/1444852 0.15860896202767896 0.1327258145510687 0.2570252958195374 0.6929437289294518 0.6929437289294518>

region sequence start stop SRR2052492 (treatment) SRR2052493 (treatment) SRR2052494 (treatment) SRR2052495 (treatment) SRR2052496 (treatment) SRR2052487 (control) SRR2052488 (control) SRR2052489 (control) SRR2052490 (control) SRR2052491 (control) treatment average control average log2FC(treament/control) p-value q-value MEG3 chr14 101291953 101292257 1455768/2043201 1636631/3258147 1497615/2643988 1278819/2677687 1699837/1958100 960649/2260743 1655354/3623531 1234101/2867390 1729232/3490690 2112253/4771083 0.6015893811705044 0.45208907524094044 0.4121720465419593 0.032485239405786 0.054142065676310004 CDKN1A chr6 36645463 36645696 106164/169351 67574/108667 39848/111440 70332/114278 56000/82614 184466/257701 135047/184961 153820/212228 90107/139840 97084/153940 0.5797185981069327 0.6962631895179566 -0.264279982267956 0.06970890125229327 0.0871361265653666 IRS1 chr2 227659612 227659781 471/1525 493/1510 556/2015 240/861 534/2047 326/818 325/1045 455/1336 680/1600 607/2070 0.2882633827594873 0.34837676517688165 -0.27326082007437236 0.031291095111098274 0.054142065676310004 INS chr11 2182552 2182775 133094/209138 291969/567402 119026/425194 152310/651761 177366/278973 80460/481158 165062/765934 174273/633564 169227/724203 288122/1036615 0.40974354597583645 0.2408760847942344 0.7664300613080263 0.022056895484407948 0.054142065676310004 PDE7B chr6 136172766 136172917 358560/2258785 151970/1578967 368147/2086541 104049/1062921 506169/2399992 119657/960632 197487/1205221 168450/1002387 197266/1376956 112174/1444852 0.15860896202767896 0.1327258145510687 0.2570252958195374 0.6929437289294518 0.6929437289294518

For example (i am using the results output from the manual section) i am not sure if i understand what is written beneath each sample at a certain target region

If you open the corresponding *DMC.tsv file, you are going to find several cytosine methylation lines for chr14 (i.e. MEG3). If you add the numerator and denominator values of the chr14 lines for that sample (column), you are going to get the numbers that you showed here for MEG3 in that sample 1455768/2043201.

Those represent the "C / CTdepth" ratio.

The log2FC simply tells you whether the level of methylation is higher in the treatment (+) samples or in the control (-) samples. The numerical value shows to what extent the expression is higher in treatment or control samples. It does not necessarily ranges from -1 to 1, it can be higher, as it depends on the samples expression.

Thanks

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/sing-group/bicycle/issues/8, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ADDOQK33K7RJ6IOZPY2RZOTSAAECHANCNFSM4P2BUVIA.

ADVERTENCIA LEGAL: Este correo electrónico, y en su caso los ficheros adjuntos, pueden contener información protegida para el uso exclusivo de su destinatario. Se prohíbe la distribución, reproducción o cualquier otro tipo de transmisión por parte de otra persona que no sea el destinatario. Si usted recibe por error este correo, se ruega comunicarlo al remitente y borrar el mensaje recibido. De conformidad con lo dispuesto en el Reglamento (UE) 2016/679 relativo a la protección de los datos personales de las personas físicas, la información personal que nos pueda facilitar a través de este correo electrónico quedará registrada por la Fundación CNIO con la finalidad de tramitar el objeto del presente correo electrónico. El tratamiento de sus datos personales se encuentra legitimado por ser necesario para gestionar el objeto del presente mensaje. Estos datos personales no serán comunicados a ningún destinatario salvo a aquellos que usted nos autorice o así venga exigido por una ley. Ud. podrá ejercer los derechos de acceso, rectificación, supresión, limitación de tratamiento, portabilidad y oposición en la siguiente dirección: c/Melchor Fernandez Almagro 3, 28029 (Madrid). Podrá ponerse en contacto con el Delegado de Protección de Datos en: delegado_lopd@cnio.es. Para el caso de que Ud. precise conocer información adicional sobre el tratamiento de sus datos personales, puede consultar dicha información adicional en el siguiente enlace dentro de nuestra página web: https://www.cnio.es/es/privacidad/index.asp

LEGAL NOTICE: This email and any attached files may contain protected information for the sole use of its intended recipient or addressee. Anyone other than the intended recipient or addressee is strictly prohibited from distributing, reproducing or transmitting the email and its contents in any way. If you receive this email in error, please notify the sender and delete the message. Pursuant to the provisions of EU Regulation 2016/679 regarding the protection of personal data, any personal information you provide through this email will be registered by the CNIO Foundation in order to deal with content of this email. Your personal data must be processed in order to be able to deal with the content and purpose of this message. Your personal details will not be passed on to anyone else unless you authorise us to do so or we are required to do so by law. You may exercise your rights regarding access, rectification, suppression, limitation of processing, portability and opposition by writing to the following address: c/Melchor Fernandez Almagro 3, 28029 (Madrid). You may contact the Data Protection Delegate (Delegado de Protección de Datos) at: delegado_lopd@cnio.es. If you require further information about the processing of your personal data, go to the following link on our webpage: https://www.cnio.es/es/privacidad/index.asp