From aa6cf59bf0ea2740e9e7fb5e7dc990daeebb4611 Mon Sep 17 00:00:00 2001
From: Viktoria Petrova <vipet103@hhu.de>
Date: Mon, 20 Jan 2025 14:11:15 +0100
Subject: [PATCH] add peak calling assay and protocol

---
 assays/PeakCalling/README.md                     |   0
 assays/PeakCalling/dataset/.gitkeep              |   0
 assays/PeakCalling/isa.assay.xlsx                | Bin 0 -> 6895 bytes
 assays/PeakCalling/protocols/.gitkeep            |   0
 .../PeakCalling/protocols/PeakCallingProtocol.md |   7 +++++++
 5 files changed, 7 insertions(+)
 create mode 100644 assays/PeakCalling/README.md
 create mode 100644 assays/PeakCalling/dataset/.gitkeep
 create mode 100644 assays/PeakCalling/isa.assay.xlsx
 create mode 100644 assays/PeakCalling/protocols/.gitkeep
 create mode 100644 assays/PeakCalling/protocols/PeakCallingProtocol.md

diff --git a/assays/PeakCalling/README.md b/assays/PeakCalling/README.md
new file mode 100644
index 0000000..e69de29
diff --git a/assays/PeakCalling/dataset/.gitkeep b/assays/PeakCalling/dataset/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/PeakCalling/isa.assay.xlsx b/assays/PeakCalling/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..5bcbf6bf34768c7f6cc12854425e08d68a8c3049
GIT binary patch
literal 6895
zcmai3by(D0)23N^NtaIP5TsN}%4O+p0cGiKBt*Jf6lo-+MVh4~l<pE#Qo50)^;_gY
zef4?2_x!QfUb}OjIoCb&o0)U}ROC_7u#u3E07xX3{CaWIEm*q9NJu+qNJzwpU-hNz
z?OdRCE=C%j4p3)9ZVy|Vl0-$j)*HZ}Ls6MUE(ornhNetwc4ZG)!B?aa4@pj{>AmO8
z)9ys9h3s^N#>5=SLfc;H?c~bUQF7xT@ST(hRY4N0sbMzTYjsDd&3Q2lSpcmkj%3c|
zQbY;!b|N#LIe{r$U$Ucz#O36mQFtm#B6_77-2A{}LC%~ICoG)A!?w8q&a^0thZOn*
zs*(9RZ@b$Up*=B&cTp&ay7uo+zHvhSOlc^&SZ|<z5E4%VAa{48D6Gmy#t9UWv*1`p
z4ry{vR~W=jAP+pkz$;oD?~KCf;aN`l08JofGt^uqddi`VS7W+}>Yq))E6BV)UW{yu
zd48s(@J3%&2;b@Jb{FBX&%(DyK7svtm|#cqSUXX8e`Sn;l*9ocDs5y6=Fz6lT!NvI
z+?Q{bet7Erwu!kFvm(2Mct_t2zoovdzH?v`NWa$eQk3Y!$y)h?I~+l4X!>sA#q*RZ
z|3m`mM<fvB|6CY|Z$?f~8)u#?o6BhJWynpfBL0&{K|;d$X~TVCSrY$NsqF@?#J0E0
z;^_-JFd!=zMwj!6nP04h$&Rur>ScSt)^=E|1a6ciX0+qRvU$&XpGeeNlp350m$Qh1
z1ihy}q`#%-o24TGI}Sz&BPO3i4sObl*|`EEq@b%YMwb-aiSa1zDAToCRPi}SG)cE|
zkEJBr!02kITTArU!qRP*lP<rdz=G<2iK+VQug-ag1#Il}?N{SAO|t3C_0GCkO%<ra
zo_LOKrv_^;WwEOaE*ol1z|KVT)4_vTnboNlgz5!K^F4aau2j>8Sjz*E@{%>*>!@A?
z`7^rfVm&<L8)gw2^hUa3<Zsqc5sY}){G{w}1FTDfe`7X~sa)%V**XNX6u+Cf+dEl5
zu(!9qcoMFNwraZpl-@qaE<e<B(}S?ZIkP=6&_6^CtgKeBk+Tr8JwMc>1XBUBFf-bV
zCzIcU*ZU5ILzNXqG0QU9A;1;I*NxjG-u%<2=D_wTeNnU+Y-Xm(QTKHrD*k0^Cm*Ib
zMG-cK-d!08AO1Z>nb4JeMatn<w|YjqUiIB~XVr!tGmU&iUBgi>ufQJe*P?sDHvhOm
zY8sAV6mTnY<rosklpfs6Ly>md6#br|DQ)b;g0}oS=g%F^CmAC6Mf$h|5n8h@lBfM_
z;yVpt$46wn3-LFM=r_CQ_sAFdPh|Su7m$1IF09f1ccSqx3A&)h*#Zi6`7=#x!d2@7
z5k$5k2nYRT`vbyz!RON?Rr{11K);zoj@TJWVId}ZltMSvD!CMX+cfwFh`1O<wmPBu
z^Z+I-zgU_%Mi3`VDTaUAl)9g?TpEQhNm+2>-@-{Bw^CPV+RV{85OP{X8Wsul>ChXC
z@NWP}*>w{fCk)5r5w^+#Ruu12SckY$QGO1MPEMwfj=5PiOkDw5Gql0r()mKOHuapj
zi(%{eMpxMD@oke=^+^UcBs1>WH7ucITouooEEkRYE2En+CwR|mX-jL@Fr;mu9GejV
z@vj!v2VnOtt8I5HGEDl9k~jT@j+C5AA`3Th9Ul|>oc8Iac7WOUA@kR$v_*%rxjVUv
zw^?Irdq4FZa(CI^E|l(rmz<ZTi(%DSb8e|Z7^-SYHaf030=bk(nLI_|s#SKq(E8H*
zf~^h)W&Sy<Uq~s<bd#+QNangxxAt11PXX2-PwwO)l7yxAUiP~fRD2-XZl)LlD2?+O
z{aWiwdXcqV<5YL*Rv6sFB~^%Eu21yyfiMr=TrQ2hSigSY4jBy(^)i=~hkQ&8Xt_kv
zgCnMfY8grs^DX9?XTw;~^3v8>oz0nHEE=iTy!Ntmoz!#jCW#xWyTBY^gOf$}3VkQT
zj(~qyeF$A#ix3twgx22x#fbU&D__7cmUo7+BocrR+dJ2|6&VnZ?OkLn@TAv2wn^lM
z2?5@lkQ4Y|ff=DG!NF$M_bFnpMqA@MyH8isHWWw0r=B=4el3Wpw=TkP57MmE9539W
z5}2$G5|e#(w<J*%b2sc%YVsV&#>IQr{e-nMvvN#9&l~jqF-NL}P(ZtxRE1}?&M<DB
z$LO2*7zGd&$$N1YgL1to%maC?YR9w>LC@fid+C~}n!ClXvi!2ZlCWST-AA-TZbUmI
zxy)Q=3lk@(sk)1krJcF+kE~5lmT!{)c5M5GL$}?m2^ivY5{Lbat-)(=7oIDmbeyZZ
z3l4jG3CTz+0+woTHrJd6ds1=@=#>+=FqMQI_u2H5zZljUrh53TLAh*On2{h41&a5s
zL1w;Zk)2_W6O7l<0RQH~gZ5OwDl=I$3o9+amBL%N{$*x@)T>vrU-iN|ynQ+N!EL&h
zlnW1=bR5P($De0UV6oCn37srO**cBVtt@dX*hVr879Yr-P=ruKnV-S~zrUQMNi!>Y
zBfS=E29s^hPRZu7S%8=xGgE7AmKR59!Y89APf!&d1z6LJ+0|!Z2e<(R8l7&LUQJon
z37x>t<Cz-Z&{(tzc~3Vjm#t3|CDG6n5ufgwlm?J556&D8s<<S0bg-iNl}7z{V}0Fn
zOvXkunU}N7iz-n>y3|hMTUaLv-#5w%yMs!iZ-xYD5P9pIvp&-gS`ab}WwTsc56)(4
z4N(?DzKX|xt==vcP+REFB7YgS_<jMFNJP_Yr@NeD5MAKX?8oPiTC=IEY(FhX@}1p(
ztJG9~FC){=q{nSU(9-|i<{%^E2-!mjKbSpCl>{_3!i|dI(#DPIHf;jZGDg?hctb#(
z*)i@DDk}GlaY`el+H%Om5h*3dcYF8!Aq}ljd(nilloi!OZa{>d*>Sinh49{oSm}i%
zPy0b00R<^A2_7Mo&@TBVXYcw%X*L)6ERd3pI*e<S50x!nKk2$RNOm4OWgd_o1tN$G
z<fhe^idK@a2A71}P80*<I{6t${D&KdLXRdvOzeHRdo)31;BeEH{IeEi^IrL!$+?`=
zv(V%nR-tG@+vFPqv1zzF1^js7fq|G@BTC|VuA7BbB-w8y1X(7O#hvTwS{&5mz8oDn
zu&oG}s$26MX!D(?)ZTB@ULQ+3QA67E(XrHhS#f_i^_C!J_X~<opB~ldY3<xPAH(yW
zR{Tcd?Hc4RxO3iA^Tq@}QQ#9^!?dMaeW9;=P?_!Qx*KeLl++s=n>CD?(qg(qv$4bK
zSuEkpX+4h-*OP;8swPR@x|QxlMqLtX`r8j02>eFt!10~26>K#0v<@`ci?PHHLn2WD
z4Q%1h{Y+pF4ma7oSsp7*K}5r=Z*<F<OEpb`AV;9qWm9(jC9XlMB|09X24SzWvyI(t
zbLI3@(QmF|hj||b!{2$<O^EyK4-gnpJpK+}SkIqX+hvG(e9)s^_U&k2aV=MJ9cEsZ
zF79=FKJobcOas;4+2wf%0XQ_s5i<kit_#2!Rq!J?hRNZgNRwcZAgQ5lbV(q8PNpwx
zbA4^qdP>C?gzv#9rdCUXCXytL+CUwM?eL6u%ESPcI^AV({H1+JNBmaIek}|@y}kq=
zr}`L?58RW^AI5nt>73zTGX%P)L(z3p@U0Fb+z;4M3Tt0hcpb-21AG|}I?jaZ#Cjk(
z&4kw6>w@-nEZo*@o03}Abq#i;Bhy+fp^;Cd-VrIJ!m0WS8Y%hW-WDDC3Yez1whUwd
z-uVhrH~Si_14!$o$zXG_mkJ(&`8?~FB<;GMyb^mSv|Qv(5A%J$)`Fy)Vg4fUyyp_3
z3=cEpxUllJVsi+5B~-Iu+r*qRw)!&xxH$x|i*miL`Y1PD#rszPCTV+u1Wz}nmZT}E
z(>#g;;rQ=J(lY7ay=zRlEt0|~x~h+C`qZ|#qK_-*oz-y6Yd=1=0&D*0>ylU>P({qp
zz1_^M!_0%Tfx8TS?#x!5UuP3X@`gp%loK@|9&aoUnBdZ({kg8<5)hs1K3Opn<T0OJ
zW1jR5d3DaZCp~HRwVTETcAWT;)3M%gStXNx&<K8$rdF=QA*3d>*^M`@*GL?{QqnP?
z_hG^J#hc}qq$1ah)<|*_4^Z*}hqRp5^{mlw-c+V%ntf2s{l~>t;qW~um$Pn<!8l(=
zagy#+9<6U&%7xcCy}#UKNTlf`@3_Oj<&Cy+ZbFMOqQqY5-NuiXc`HNdZPYEuG=}*e
z4nsfylA_5Ar`v5hy_zvWyigPfONENN-lL8=Zc4M4U-ss6nE`Zl6d}Xby!b&cdS$H-
z#-?GpCzV2Eo+x&OU(cAF5(=0KDTof%j=_nW5f6q&Ut1*3w6}J9qT{omUbk;EWO;}-
zhe0o@aF*D`)8I9i&drxszY%Wcskq!KQWD!4lmcP67K}_qv)*=XW8F-Bdo4UkwsDu&
z+tzjogQ7rm0S#_2hAtrXqPmyxYy_>(0sBNcz5qylk89+avsAOPS+Wy742(w3>62I@
z&P3!K5D$DHx{X0mQu~y{^>qLGM5VZRkM6qCyO)G;=-bXhoB_V;Q~9(2rQ4hs!<cV(
z+J$p?n<H9O>>5)+4&XooJ!h~7cdfRTfxPy9$t#2E;%e6H5u`rlL~p)<SxYZDmYV3q
z)ssTbG%SMm=MloLGu)L){5bfrk4(r;Jwm*rH8lsi$Cm2N+I;ufcIPl9qxGbeOA7O~
zT-*$d@5ZK-Y8lr+vV_EU_;q6RIJC9%8l&`W-8|6ZW`e8)^3+Uu-YaFD!O0HebeX~T
zu6_8TLBY$ywV$<FO<gws=sIi;4pAB_m&dCy<f<~P(c7wn!lnyEP+oS@al$@)^f@xm
z=wEsUUp;!);Zw_6Rd&Q#siTLfU*Pb_(WsQf=p!W>=l*N0RiUCjrIOL&>eK7_)0T-Z
z?w6>qjW;~1nK}3<Tk*Y<KU)3uNa-V`k_pza?sRP@D>upJA<nThNbct7bX59?FMo|d
zvG*?QhJE+?RP$AjyKG+Qcl_?VA*N+iL<7V7XXbG6w1NH@HgXdU5hF55@b)1G=Zrck
z9%q=Jm1Vd_K?=J-g4?>ET=5Ai`AD4N>GmBAA|ltRJb3r}vY3qhlbdy`ITAd*{02oN
znA&2p3Z(UmR?p_kQa+|U9Mvg=uvb1{6MKy@V@ezUjt0$Z>ZTbXP$OZ!vFac`HgEE#
zTndNbO?1Ymg7JDQFQ_f8+AKcQ%BjT;9zJAX)Apd_gzkz1<@%iUy@HYis=Ry9YeLe4
z;y9BI@8sfr!giQKzp0nW9jfS2JSbll4!#fiw9<VaIdP?M9^Grx^)6W;88Oytdq#A^
z-8PdsT{TN6rN)Yu?xauVosb^$+~%T>c()uoyW=HsY-5Aedb`kE*z1Hdud%C9-Odtf
z#klJ<8Giqa7jRen-q6EVE>dsnTH-F1EPIMsCwiq)NfCf93Z-4L>MfT+kUhS@b5#cc
z426utn2ofo5W2{I=C638pRCkdHB;dWtk5vyr$Jc3+wjt=dS<&py~!0$wp*VpWrCz=
z%p;MFahTf^U}ZCfIBweXW8DT*r+gg_yPwu#=s^^1$TwZ6zNZSlo{5q%bnOz)I+eWk
zH6i??T;5TH7v%%VOxvEIFKYzBPjuZ9CVWTh-#h7|gRU1@aBM8`Kc8p&)(iREy_0y*
zX;evFk$3bqO#=v{D6Sxg+SDJVNT1JeFN%pcFr<0p{8-@7Y=U!-Y{+U;reb>RQ5w(L
zVDq^v+-hdm?0pDEER>P1Co+P%Ph3$ieBgG#x$pNDc-^54=rkj{<#5khe8-G8&1?VY
z_TC;VU}D$C_1c%A<inMy=0IhCV?J6gu4f7k9o(Dt{hGo|E8n1miylO#CruSIR|)#>
zLzt=kLlq}`hd*z9U0>Rufr#0e8F3E>BJ%jBnTdnL#oU}0f7gEc25!*yDQWOlRH_|o
zu!I(G4&@`AawNIwg_bTCPK_kAIDz={lP>)E+v`d&(RhM^47S2)spW7Q;t#VN<pzn_
z684g|=S@U0FnjF`nY*}67}rF!=0ln64lt|a=pQ_Y&{#xMerl1t0D<svxk*Hbz&(cp
zN0FELHlY?2);rwKC-6mk*55S5c^H2c7}s{hCdqHt9QaPDhk0{nXk-iIoUUc`z-1h1
zIb)PyJVzpHZhrl3rn-}u#x0A8B>`QrHMK6!f@%z^Qs>upErjgH<0^-MwthilV#~4p
z!h3D>q^6ITOAI6CuL1FepRUmB?r<tqtMgi23rq0DUk%ceX1<=k#Jc>st~%qZ6ThjW
zWg;<n{8{Lo9{kmN6;Il~)7|e6yNAL1uo6xD2_SkKI6!;K;Y(c)5ZhK!T0<7uaNJQ^
zgB|dE1-Z}Rg9hZAYzl30u-4DGM;I0GTr(^aD9_f(&S=&vmd}la!}3QytbggtKGR8B
zxSr|-X)>(+hI$p)%gU2y$q{*R+gU(VAHv@w_R!u5`Xgj&JqS?=*sxC;w;h&1rb3!q
zfcA=oon94f-fY9x`o&~Y<MqCo1(6qm*(`Q<?o=%J2<ubzFdQMT=V&>{OY>q$ICAAo
z^><A_BJT&!#<IP|E-_ZJaxjgXG5wfhpzl#W<tKy1y(8P5`?0qRe<&yMc2hunT_;B(
zi_laOu4tvc6=<e{Yfz4AW`rmfD>zaxIE$Q0&Ytba*b=|m0uL=;MFI<(S}{8?y$Vq6
zqK%`sB&iE2#xr=D-g(!@#R;gTNvb)bwRWn!@_y~OidrDCP?r=l-rPU`BZ`#mzEG8R
zy&8FAr+^5Alg@s`XNA%I=maC$@DY+2*=+4ojPtbHPKi6qlTwD39I58_Ue%FEO<L)@
zkCJ<;+z8=#ycLDEtCHb$E)PAogON+1`LbBnyd~Zj^YPZ>krNAJaW_c~Oz=ok81YE@
z{4PUr;N?k$$p`U0`}YZV#b;KGmi}l1W?WLf28dv|Ac8?~*#>^so}%Ivl@a=0@b+0a
z>tP930P@(B93a;1E|nk0o|=0^5>Lmf&UaC}dvU;P+>RUVF-CWEOO(e$8D)n*Yrc8D
ztd^<XG-RjURWp(Qp@bx{MZmbB4FU%{85$Lf`H0vwiUa!HHM)EtS-pH}UX>Egw1c~H
z>o0@a?LI<D7$l{<+D#SRl9zjW{oVEySQj3Hnry3DG0Sd8aHr7|;xTC{dr@0t>mJai
z26?QsDtJ<Y%Ex5$eFGsTnspy;(VAG-uB$tW#(AD%+E#_9d#5Z#o{$~q&oB&x&wq4k
zw5oc0xIc;>xP$BV>Osy8iPiEUtQ3fo;)3|DHl5pC(@1a{Or|MT_k5p-N0on4(tMbs
z=*7xeUQba2YtC#u)`I-#Jr&WAkdLa=KKi_+DQ3dNb{aWYVt#;=U1FwaDgjFxEF3n0
znaGqSD<mkf+!NJe(;JMdd7V7=?<f=cPjFP^kx@vH{<|oSh{kX0UnKHh8~j}{FSOjt
zbx{<<5fQ%9AOF+opW5yf3JIz7C+3F$`Tse)kU%fjg|ik!D*X^duRQ&Ix(iA4a$NxC
zE&>03z&|C@D<I-k0>UlA&n4{lr2U`gy-*S_*M)mFgzAFW=Rb7CtD|35tUr~-D>M?4
z=3m%fI^!?KE_eA;<@vP>0p^br{H6K)+Tnjum;v~)!<9<$pU3}Gr}*Xg-S4nps>LtA
z|2`}J+8_W?GZ5T4|C}6uLI2m0nTURC*1uqv9{=p>zdXwS4*RFa|B3#y!~8;@Bi=Us
poDi=1&7bf;YvM0>IO)$+xl-p<<k5dTYe2;G@eYxZ#Mds&{|}UWM2Y|a

literal 0
HcmV?d00001

diff --git a/assays/PeakCalling/protocols/.gitkeep b/assays/PeakCalling/protocols/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/PeakCalling/protocols/PeakCallingProtocol.md b/assays/PeakCalling/protocols/PeakCallingProtocol.md
new file mode 100644
index 0000000..bc6e831
--- /dev/null
+++ b/assays/PeakCalling/protocols/PeakCallingProtocol.md
@@ -0,0 +1,7 @@
+## Peak calling
+
+Peak calling on predictions and the experimental data was performed with MACS3 (Zhang et al. 2008). The sample bam files of the experimental data per species and dataset were merged. Then peaks were called on the merged bam files with MACS3’s “callpeak” command. The parameters for calling ATAC-seq peaks were the BAMPE format, a q-value of 0.01, keeping all duplicates, using the background lambda as local lambda (“no-lambda”) and the ungapped genome size of the species’ genome assembly (see Supplementary Table S2) as mappable genome size. For ChIP-seq peak calling two parameters, broad and a broad cutoff of 0.1, were added. The chosen q-value was the default 0.05. The ChIP-seq peaks of the species *S.polyrhiza* and *Chlamydomonas reinhardtii* were called using the format BAM instead of BAMPE. MACS3’s “bdgpeakcall” was used to call peaks on the test species predictions in bedGraph file format. The parameters for peak calling were the same MACS3’s “callpeak” determined for the experimental data, i.e. for paired end reads the minimum length and maximum gap are set to the predicted fragment size (Table 5). The cutoff value, threshold of the minimum read coverage to call a peak, was estimated by plotting the average read coverage of predictions around the TSS (see Fig. 5b).
+
+Different cutoff values were also examined. For the ATAC-seq predictions of *A. thaliana*, cutoffs in the range of 1 to 25 with a step of 1 and for *O. sativa* cutoffs in the range of 5 to 200 with a step of 5 and including a cutoff of 1 at the start were chosen. For the ChIP-seq predictions of both species, cutoffs in the range of 5 to 100 with a step of 5 and including a cutoff of 1 at the start were chosen.
+
+The selected parameters of MACS3’s “bdgpeakcall” for each test species and dataset are listed.
\ No newline at end of file
-- 
GitLab