From 2382e65eb4c10bc1258885b02e3d9fa34de59270 Mon Sep 17 00:00:00 2001
From: Viktoria Petrova <vipet103@hhu.de>
Date: Sun, 10 Nov 2024 12:28:04 +0100
Subject: [PATCH] add Quantify metabolic activities using PAS assay and
protocol
---
.../README.md | 0
.../dataset/.gitkeep | 0
.../isa.assay.xlsx | Bin 0 -> 6912 bytes
.../protocols/.gitkeep | 0
...uantifyMetabolicActivitiesUsingPASProtocol.md | 13 +++++++++++++
5 files changed, 13 insertions(+)
create mode 100644 assays/QuantifyMetabolicActivitiesUsingPAS/README.md
create mode 100644 assays/QuantifyMetabolicActivitiesUsingPAS/dataset/.gitkeep
create mode 100644 assays/QuantifyMetabolicActivitiesUsingPAS/isa.assay.xlsx
create mode 100644 assays/QuantifyMetabolicActivitiesUsingPAS/protocols/.gitkeep
create mode 100644 assays/QuantifyMetabolicActivitiesUsingPAS/protocols/QuantifyMetabolicActivitiesUsingPASProtocol.md
diff --git a/assays/QuantifyMetabolicActivitiesUsingPAS/README.md b/assays/QuantifyMetabolicActivitiesUsingPAS/README.md
new file mode 100644
index 0000000..e69de29
diff --git a/assays/QuantifyMetabolicActivitiesUsingPAS/dataset/.gitkeep b/assays/QuantifyMetabolicActivitiesUsingPAS/dataset/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/QuantifyMetabolicActivitiesUsingPAS/isa.assay.xlsx b/assays/QuantifyMetabolicActivitiesUsingPAS/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..fc2e300bc97ad14f89a423aef6262fb02939ec5c
GIT binary patch
literal 6912
zcmai3by(Ex(gvxe7pawQX{13qg{8X%l%=~{q<d*W8W98mK{`c3>26e{Te?BMMGoq#
z@A=O2$6kBw&OP&7_sq`B?_Om&1Vl7AI5+@Yims(j>~u4#4m=#(HX<AxF6>i12|HV7
zh^@1sx`#c)$$-_}#=1B`!M5cdFz84~YJmleVW6%d)sj=uO;E58H|#FXOftQj{BGLq
zHhm!@MWGQcQwr~vS4P_%rK%{|F_4TvYJ>_89_rK(gAHNraoW2)C{;E<^Qi-YQ<=o=
z_&Hm_8IRn+RF=h@=s^)#Sx6M7@{*u#i8?D6@PvmsH^dPYJ>jTzHh?)j%KR~r9=1wk
zzV`dBwgpHx^hoR$(e2K?hZAod;lB_Yh%eOX>m7u|kpb?wxe^st=EI{03d))@EyIU2
zxn;-?pvB(_JVwTRy)f1hh2G7!oHPlE$7L|kSiPOdq=i{+vVh>9L&VNQw?0+`Z-a7i
zt|<RjPns9Yaeu25=fr3JyQfcJUml8#gISEN&|zN%R9`~$00)6QG8N@`(`Po`z)*Ja
zyM-U7TAxkA%kr;+J2;rfmj>U{-d716SO-$Bb-xh0J$brT_DFy!Xbn-%RitQ+SouGh
zfcudN82LXBGVGn9BgEQ??Z)OhTf3REQ>(D=<Zi*iq5rgDy|OHhd#~7f4?}dzTWaC#
zxvdN!`(-^v?q@nK;bt0J;>xHOZ2{l5!eT@*qAXCN9X6KDy4QOJqt>ES4@oeXUlZXW
zb@zqzHFtlvaKL6nM-HJz;jqueNL@0$kf#P0bT&eD@MSuXJ!6lvT&jc=l9{6Mx|F&t
z#2NZWRzqEzqxTC-wwzBp{gwg?s`^By>S$kGun%!t+v?e^#%>ztP?+hSchQ^3lZHL@
z7}-h-)>_JDR32D1&>XKn7s}6&8OYA6N;Ah%D^Q&4)_vzfGHrmm+#e|?UVT|h@;u0&
z+D!-b@j1uPUETq2xEn_PW(@(xh`aSq%KkP$y*BtaW&>%;G_RPgg)vL?yP2Dvqtzoj
zJFBZ-!VS@ut@nVETNh|$N4l=MV1`&HhNt>^M+kuxRr1!d=Dao+M;gR3B!Fy`%(kM5
zl#ep&y+?eZO7bHpZ?YJ{z!ink#w|Q=uIV!~VB3_Q5F!+fj%H%SZJn2dYnjy1hbC4*
zkiouZM+(e=^*})?bY)M0c<2>>_ekff-iL1VT96Z(;ZF!_=t^bfXhVIP6wevvo-{~I
zA0itD@JFtkfCFhVf_vDA((jufJuomKj~SoWlKbHFrQPW?QxNO59tL)V=2vI&vpxoq
z?S`<EV}hRfxO;|_o1K)qcNVx#rFuUW-0|3%UnBqTL}OhObVZGmIRxVTXPH)qtJDR;
zh-`ro4*JXX2Za5K&%`7ZyVQF?znLSZm>FU|UK+|<g{~@<vZ-7)>4zI2+@f33Rq<74
z2lae%3nf{j*s*-X!dPcbX?v;5B~e)7#097R&CHasE476t@0dFJL(X2~hebks+I2@G
z{2Krgwq4jK@k7u&oEB-oih>xCRfr1-@t4r(loTRKC}-smX*p=kz#5rFdy#ByDw(d6
z>Ra+gXISajmhr2)Bz<eV8MmD3yP*Uu<;hJJ3r2kv(eF^k*)M9yOKR4TC9NS$n-KwV
zujbeL>mOQF+3b{O8uuNiZ2I#aD>@cO7H(oVJi+xj>(xnXmtouk&k>Sn2@U12cCZv}
z(Z|&EeC|DB?X<gJDA{{hd{L4ij9P2O{7nT+RasrU(N5?9WKqPY@en#xskH5Z)RjEs
zX|dOT<Da{_h)-;)lVWv%H`|5qZMQl43}6NFU`-jsi(mTaWw(P&!T}=hqJd&V$ehk8
z*IHgs3a;%KrMZ!|)E~M#r}1LX_6mJD;A6v_eMx2~+@}|~O+dy*y38WsE(fg!Ef<Ts
zGeN6KmLX)&@6cx+4WmKJOW)3Gt<Mc&5b?d{w3a1nC6Yy&MDMBW0CRy2j^;TllpR#t
z-2P#8Ar!IAyr^_wa)13>hIGlVd;vqK-kC<?Z~zW8@0Ug`@PIfp@7G4$Pka1hngs6|
zV`IJzIXxUGFvT&!KG@8@oWk{Lv@sIcNnB0eP#6iHdTLL-UjVJMdX4NBq*0+UR``vC
zd!i~xSo)P%ae@%aPS}~m#08Lng<Z_;l)fXYVpLw&8}#uBQ<^z%K%1#Vxkr`u5Js*0
z$lExmJctDEqsU$TGTkYZ13AqqhxEyyXNOOED4Iy#b&1?m`E`S}L_L<YgN;LJ*f_+y
zu3RT`V@HUInzN&Yt(nu0s!iyHby{HHmSE7LfXga8WdfqqtJoNGR9SxePxq%k#61Zg
z7wes#rwH_u#u4&7N_F2Ly<00l)ca*H-jQcftDN7={H0aG=lpYsn*JAiX?SadDIIwa
zXxrA?h(P*M`pSe|p1{^(AdMcED5sk!v4FrbW?;>}mlz6kKsf>wNf2>zt_xC5;#|E?
ztlFe)Kc;gqHhbYlDJKJw?bxmBGsMMtD~S>Gh6-ojI!HAAK>$yw!jtAp73+X6rTXC9
z*#Ya%DJ-W>6_x7^lA0_L4x`8rt8~{h-<C1*6%wr%L!9!7L@oC+97>gmCW(5&>_AgA
z83LtZo2U556ZeZvSOYw@#3&n;BNbA{#c_IscP3SacT6&kbBpvf+tn&{H4kc_9Kjq9
z+o;=eqb5(J@VhQi^a_H7CRXsd6#ToBm4svyY`n3DGrsOGh@5#0AM<Tmm{zOliiO4|
zFEHjioXTak@iPdFMoG!Hk&7L|-z4b2wsU85hz;at*T2qVTp#z{NZ1r@qqtscU}NFh
z?8obm=Ci4zWH-%&cgg7gt;9rcH#5uDxZ8D@$HM=^<^VPIFu`LmmkeW=3La=`m=yup
zxs?^cb=nxDX@sP?@fI65t9{HTR7mza^^|&QmBpa316*pZ@7B)8BQkQsw%6lI5|$*7
zSpgBcrYGSxM0~rGF_QC19(Dt4-0~7Kc$hd49NQF5=AQNOk{p&hUxCCFq+u*090&~g
zdP%h2An7@@)HybL&ukd<6dB3e<@N~Sp6X1oX(+rdSJ=Redq7<xPbqG9ZPxl0!G
zMkd^(Isd#_$*e~%cVad-?L0JPo1Qls$0p@oe@r^Yb^#Y=cwit3%dnzIp37!oC0@>3
zQJ%ZwN+M3RwaxabvWv&Z_6#e0C2Ce|2U;8#sWlH9wbn<IPF3M{eY7ofUX(xFN#o~1
z>3UA|`Lk!WuI4uX#VDrtw8D2hZ<io9p6!dK>bJ&N3EZDC8>TH>>k56{g5KDCU3aZt
zAHnxVL!%F)BQ~Ed)@W?Ed=`WAVp`XI*yZ$~i=;_hr*@@lfm#R8it_%W25i5PTA8?x
zm~sX(N^*O$oP`+N$03mjfCh%}WIyBjM@O5C-gloUPJxBOtL}Blnn^TGfWXI~mSq!0
zy(N|b%Owgn!v;RD^Ye|JEi<KzG@<V<!bf?Zc)~w;)Q*ez?Db<C5<R&*oL|qMS=*t4
zJ~`;tdh`8wPhst)_<FtBn+y@JlZ){u7w75-c23U8A=omZK@KRHAU7QV`iQ(A_DQ`g
z#w|%aR6ICUgpE#7_%A7xg{>~7mMv!_96?y_)WWJYWQc-Ek_Zi?foS&6*r$y3>(iz?
z^-mVt2DL@_p?fv;0MhlP!!eRi5&6Jf$^0SoQgNqD|LQ@|18t&CPM-JL)Q5h+_LBOx
zW%<%rE;8VYfY31-1V{P<@o5^wcRkLC??=OJT(^ixrCrwQkF}**szlZEZ>zOO@+vc{
zyn;lEKY#Fz;?9Z;Sx?OykiN|JK1j`UKUf=(-b0qj;A|%mJXr7ZtWO-jlQv}~Mj*6I
z@NGBUL%)`SB+f8@!Na^{(NL<#nX)XXdEa7kv3*5Vvg@~SnP+VDW&$vBv11lwdtCHx
zak_}~tpJSEcX_ZAH>Q>(iAmGlivka^KH#NiQGWQ)n0j9@l|yJ%58fovrl`D^CHI5n
z5VX{fgQ38RE1FgubrM0)6v^Aw%qq+*I0v{x)$2xQ*|GmMVK{F{XiX_W9qj(r;(+E*
zGPLidi-;&#o7N{AI(`TG#RdAbXV9zTt6S33HecIm3}E}I-yI72d(Nu_l8+i?-lnUT
zX*2Pv@@{rvj_Eey#;p{$_v=p1`#yiW`~qK)&~Od!Wx~O&e83SovsE2^^dWm1%`=T&
zh{oQNBFpf@T?mVlPPhITM`lryP9mG;cNV2WT4wJ>PO1d5jyvrFR4m?z8yCjp$is?^
z72d5}m|6UpitnTN!PCfQyXaH_0dNY&&mHf#YWHYBgV-Urz<0})-E=+MXIY6&Uo7s<
zy`%$B)DndZS+QdUJ@1jWIvAa<e|cKLOW<+Kws8NP#xcHtrjY3N!P<#T!e+#y!I4t)
zgqgOME)OKEy9l&)tp;}=BhDgI3dx@*bh0&g&1SH2<kfA2n|dfLw+I%;bOfb>sR)DN
zNyye)2{+bF)wb5clcXDW*u8CRmXL`Ggys<s^+%Dog`ZdT;G7R5_S&PJO2!ocNguEb
zKXa0Jr(~MqNLep~c!$|1p;(0Gwo^bH@R86KGEs3&B9qJ69_@IAh)B22y5ffyIERq;
z9fj!q9JEvU<N(F{%*aD1@3-6da@pTSG%MRSrh)8b0`+yBWZYS6v^4eQwDyW$=~oq1
z(dP`q^(rNJbM$|;@RGe-9i6axTF9J^iv95-g3o1!wIYcN9V^DunBdGk#5-C;qrYo(
zsm`?3caLFb7DYT-S3;?{FkjQzRo_S~CbdM<s2ZHjE3(a{4b^4R(#mU$(z9`OM~s~b
zvgFQFHDUXxn0<aoa1^UUC-Z=Ca#5X#{VvO1_GT67n>kO~`q@LU;%J#1X0-uJr9rjs
zw^~U3bb;Wl7abJL^^>1`j?FUrmYyB19)D=}siCiYbIe?!tqU7r_MQ%gC3uFPh!L6h
zN;Ow`U-v2&j}%p%(dJKEBs_muthP4R;8{I$@JYJ-vV$vHt#r7=Q?YoQezYq?%hA$R
z{M{h)Xgc`i=E!tZ#;`9}HFuHsPW?T*uJx&RH#6>fcp(Y+J#>Yw&j_#qhWVc*#M#3d
z@?-mWnP33hnel?Rj+mHd)DSS4!~85Q!qp2>8M)(K*ZpLRP7&@5$10p{2_WCT?J|{j
z*!A%ZG;{BivvxICl&yzL|1}<pmaw!ue%*rQv$;2^pHd%>XcvMRD;_Zjmm<%YkjH%>
zL-d;BG{pg`$ImrZ9>m4uO>oMlG8u3pQ77`m>8?B{wXke8pRAEpjTtz4OvRw(PQeV>
z5dq5fI_Y@@C2?1JcOzAYWCX=BCmjjA#QcnAKZV4po5dQc;9fK!_a<ECA?Wi;*F*S(
zmBKkBuT2**f<OXX)KWWYB%GaA<5?Y53kb3Lil)wlPnJMPx7o{guRme(J9Ko#iKE-Z
z1gZ9PB04kH@?{Y+s*>E#=50Y1bDRi&c+L(G6L~QBxP=AZ+o}e)Q#sp?=&K{8Vu`pQ
zK<5^*ZHfxN^8m;Wi#u7xo*P*{^9Z_;o*hCF*+;jJIr7<3twkg4aGo9#W|SC&8oYH_
zQdvi5JD@wU!py+`*+MEvg3K%u-UywpEx!KEOd-0f7UgJ{{?r*qyZz4RH7F&Bs1=^m
zndCB!r*tMt%D|;lB>PO9a6dl$lWg8`gctE6@hqF}pv5)p;HNsS@#DTDbsrsdkU+Gr
z&6zfqxRU1>zW0E?bZsXbbQo5UmggP6Pge)l6BU(XM{Vkj5M|6|y1j-*92k&!Iz8b&
zG972$B^b2alq#Pd^-O0wA9#1+a%efTWBM@!IR-*a(H$8<+AE@<8{U6E;KKK^`LOm#
z3Urp4(|ojRC9-YGp6<1Ge1CVB9x%RR?LxRXm~yld^)67!--v^ph2@#NeLL%>U7rRY
z&B}KO&Vu`Glhda1nVSOr_a)54?y<6?o&BGWft`!3kU-cnO$U202g2(3r<t+6{na9!
z9w%lueGemO>kL2mTU44YeXyt|doHo3b{U-P^n7!tGqZXUVk~#u#c3zj-2HXMdZ9S%
z{!E6#X^G`<GTh0pOlA5BIihyrHWy8|q4joJnNngHO~{0TnscEvwg)JcvXqY=MW`<z
zDkYky%!9!kEUuyvf`=YMfg|wC9GehxBCBoI<Z&#a?)A3~vF=9u++$h}Xn6T;8vU2V
zx+t96gTvoIP8phpkDSNgmNQ4N$8trpXXn=6XQ?>~tMi*jEOF}yuc>x=6jULTlsJ{z
zHsdg!j42-hTKjm6a4kmn3Lmsm;+s5KE;fjmBLrgcC9Y8FY%?oXsj*uUhQ<40tp@2z
z($VHG(Jz0gt;*bY<T7!vNWhgDdlo9tjkWrq{At@I#lya^2gvM?D-cDV0;0Eo{p4p%
zzNB>k;VoswHF&`dhi%2R`UAFoklQR4s9&zhy3ht4b^V-mm|7mwCDS6F_<Wt<oNT>(
z`ND`VEPr@%eX%#^Tsvu=Hq8s%WKi=R;U=-you?G3w)G(_WAm`?gY$dF9@{xWexyvb
z`|8=JVGloct6-i07(S$aYXmnj#vEzRbi>B#`9xA9ZSTyy;Pa}SyS4%X<?}v#sia+0
zr|`>(gK90(s9~b$)QMl@J1zO4<T4AfR-NcCjKP-MCP_0UEqN_`?k*qvC6I*Cq^I&*
z6gvq9^AbGEg5u6plvhMcqL}!RS9nGo%jSu*=4op4X(^tO@Y`TK7T}Bfq8k2w=$pYH
zk&IZj-)%0>F6qdO3P8!e7P{_|ybfa#i+&<mr<l*SB~Vj?yO6G!YrlWx<J@8~DYs^w
zP9j{KnQCqe{Db-3#4i?VtC1OtL}d5`G1aTW7_l=?d;*$Th5UeW+bs+cvx=l1v4G_X
zF@rmIeT>+Yofx7^*-O5?!bQvaN=7*Q$qH$=O!(mRS=Z@jS5Cjk7R4dkfCo!HyZa4k
z@!n~Xc8@<iN@Bhv_fgjKqi|{l+nQ1^<g()(DHQp_oA8q366v4GFlCYO)rTd+8I}y}
z>oM@V7!?(#pac{Ag164Y>5qz80^modWC1a*VkCY{yQ*#xNo?(_+LuDM4`P9YtPUG(
zP(uNoVx_TAYU!ab8gG-ARkPHZ25q%EtH<*vi}4bgxs4iH!G|)A28Kn#K7!VbB7i<O
z^-dpfb`OWDSB0n(`M{3s`ir18+fNWYDsc&~HWLNcl;!Rof7e}k`uQiICY#C@lsESy
zSkozSFljWEyhzP+bPmYVg4|bH<UNQ%WuwyhzJXw4joL|ma$~ESbu|Z}SdTLlo67JE
z@6@HpQ-YKH8LIy9xlfLbmX+_1_C}Bbw=rB_J<7c&x>`1fnhJJQm>1d6qHvvU8V*jc
zmuia9x!AkSrouHLZZ^d9`uWOvUia$;`rNN^sPl3o50r&MLO!XG`slHjq?+>K+N$TG
z3i|<0cW`N<Nw_V@P|+E<XChOVEa4!yvQJftOzu&y=5?^yeISnSJ4IKPgTI9b_uo}=
zST=q;|Du`ya`1P@yb^b>&s9?lhh_Liko;%RKgHe~6dYX1Ps|TR^8X8VrHNjjt6<Hr
zQu?8g-bDKQc2^qd^|=DfUIYI9f`4kJH$d2b3$U=TIM=Y>i}rus_exs4K3CxxU@{EM
z&p(95o2y@UtUsm38#El8#$VW9!s9Q-uD$#z`~32Pjq>9Le~Cc9Jp3;TO94MT+^7`)
zy#Aj$#jn8Qzr%j17Qf>D`>y!QK?1C2z_@ezxjFuV{;w-@-~Me_|AJjd{Bx@Rim3lP
z>_0R9PxPM?<`)_Zrnmjv5N_tppYT6>;xBkA{?Agmk?WP^kbeBufUxc31FRxx2(Qfl
E59?WH{{R30
literal 0
HcmV?d00001
diff --git a/assays/QuantifyMetabolicActivitiesUsingPAS/protocols/.gitkeep b/assays/QuantifyMetabolicActivitiesUsingPAS/protocols/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/QuantifyMetabolicActivitiesUsingPAS/protocols/QuantifyMetabolicActivitiesUsingPASProtocol.md b/assays/QuantifyMetabolicActivitiesUsingPAS/protocols/QuantifyMetabolicActivitiesUsingPASProtocol.md
new file mode 100644
index 0000000..856ea3e
--- /dev/null
+++ b/assays/QuantifyMetabolicActivitiesUsingPAS/protocols/QuantifyMetabolicActivitiesUsingPASProtocol.md
@@ -0,0 +1,13 @@
+## Quantify metabolic activities using PAS
+
+Pathway activity score (PAS) was introduced to quantify the activity of different metabolic pathways in single-cell transcriptomes (Xiao, Z. et al., 2019) (Kim, J.-Y. et al., 2021). It is designed with a permutation test along with a P-value to examine whether the gene expression of a pathway at a particular cell cluster is significantly higher or lower than the sample average. Since we are working with bulk RNA-seq and not scRNA-seq in the current study, this algorithm is modified, and the permutation test is no longer suitable and discarded.
+
+There are 3 sources for pathways used in our analysis. (A) we obtained the pathway data of Arabidopsis from PlantCyc (Hawkins, C. et al., 2021) (link to the tables https://pmn.plantcyc.org/organism-summary?object=ARA, downloaded on 18-NOV-2022). Pathways with patterns ‘glucos’, ‘galactos’, ‘fructos’, ‘xylos’, ‘sucros’, or ‘maltos’ in their name are considered relevant to sugar metabolism. (B) we are also interested in the potential to convert sucrose coming from the phloem into different sugars for secretion at different root segments, and these pathways are not clearly defined in PlantCyc. Thus, we applied flux balanced analysis to define the pathways of relevant metabolic reactions. See section ‘Flux balance analysis (FBA) to define pathways’ for the details of the algorithm. Table S6 shows the genes involved in each pathway. (C) we also assign the SWEET and SUC genes to their own single-gene-pathway to facilitate our analysis. This means the gene SWEET1 is assigned to a new pathway ‘SWEET1’, and so on.
+
+Brady et al. performed bulk transcriptomics on different segments of Arabidopsis root (Brady, S.M. et al., 2007). The experiment involved two replicate roots. In the published dataset, the 13 segments of root 1 are labeled as: *LCOLUMELLASB, L1SB, L2SB, …, L11SB, L12SB*; the 12 segments of root 2 are labeled as: *Slice1JW, Slice2JW, …, Slice11JW, Slice12JW*. We found that the two roots are in different developmental stages, therefore we dropped root 2 and used only root 1 in our analysis.
+
+Given the matrix of gene expression across different samples, we normalized the data using trimmed mean of M values (TMM) normalization (Robinson, M.D. and Oshlack, A., 2010). In practice, this is implemented by the function *calcNormFactors* within the R package *edgeR* (Robinson, M.D., 2010). We used the function argument *method="TMM"* to call TMM and set *logratioTrim=0.3*.
+
+Let us denote gi,j to be the normalized read count of gene i in sample j. The read count of a gene is normalized to give the relative transcript level, which is 1 when averaged over different samples. Mathematically, the relative transcript level of gene i at sample j is denoted as ri,j, and is defined as ð‘Ÿð‘–,ð‘—=ð‘”ð‘–,ð‘—/(1/ð‘*â¢âˆ‘ð‘˜ð‘”ð‘–,ð‘˜), where N is the total number of samples, and the label k goes over all samples. The PAS of pathway t at sample j is denoted as pt,j, which is a weighted average of the relative transcript levels across the genes of the pathway: ð‘ð‘¡,ð‘—=∑ð‘šð‘¡
+ð‘–=1ð‘¤ð‘–â¢ð‘Ÿð‘–,ð‘—/∑ð‘šð‘¡
+ð‘–=1ð‘¤ð‘–. Here mt is the number of genes in pathway t, and wi is the weight of gene i, defined as the reciprocal of the number of pathways that gene i is involved in. Because ri,j is centered around 1, so do pt,j. Thus, if pt,j>1, the expression of genes associated with pathway t in sample j is higher than the average over all samples, and vice versa.
\ No newline at end of file
--
GitLab