From 37f983830d08e49855035986db1fdd3cb2fe0225 Mon Sep 17 00:00:00 2001
From: Viktoria Petrova <vipet103@hhu.de>
Date: Sat, 5 Oct 2024 14:59:02 +0200
Subject: [PATCH] add gene expression comparison assay and protocol

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

diff --git a/assays/GeneExpressionComparison/README.md b/assays/GeneExpressionComparison/README.md
new file mode 100644
index 0000000..e69de29
diff --git a/assays/GeneExpressionComparison/dataset/.gitkeep b/assays/GeneExpressionComparison/dataset/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/GeneExpressionComparison/isa.assay.xlsx b/assays/GeneExpressionComparison/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..b66528c2d53360cc64e0d680fe2abbefccaad75c
GIT binary patch
literal 6903
zcmai3by(AV_XepU-E*T$x}+2k1W_2>FhD{^cOxks0wY9_MnYOZIwh3uP*4$&?ixt>
z8}gz)`n<pQ^T)1T+quvAT=zLU=X~!e%b}rPML|IUpd6Y?>%>g8VC$fwplqR|ppYOx
z)s?cdb%xnG8>oBO!<_WF+-<Cj;}vXMc?g3J#H7D+LU8rfHKbd!UiXma@1cyiOL9<7
z?Ibo&xe>Dzu+bJ6lCUQUZF;7)UsI}zlpP1kh$Kg-2$Et?4zt=&)E=fZ=SI_K0yG~x
zkUN!05y#EiicWiE2PSjQXGIN(%gVwc@st-tp(W~E{Dj{GIkH0>v2o%L+GYYcQX|bD
zf^-Q~p5$r2>2CiD>xn+Na|J})wfo?Gxg+Xlu)gHiIz8R}kXUNKH8)pKL1i8)PN1l)
zIr}1NNRwNd{NUBNYk`NDc!gibJ0o#=xEB*Xz~V?)^);4>pR#M=RhxW8^UnhD3No*b
z7opl<ot`Mlm+RgW!gt);>>~Q+J$LrVJFq_&OUA(rZYzf9e;utSC9zM0M)M>Y>u}wB
zCQjc#cK*!54^OS%CO)U)i|95H-r;xsvy?YgBKy{XbSphC#E3r}t&|yyum`Q6>$-{;
z&4QKxlL?d`nLv{N^I#(1892hMowzS;F0!?gE<3r5{7&u)3JT6o8?JN9;@CHeZ9KRV
zn_kjikDuGh05WsxX|q2u^WSb^v;|j2zGx5F*bIY9;6_?tMLDc3n)R&qiAJtOsv;<H
zISN6f7(M+V{VhFb77he#IG7;}SbX-`xXBBqr}7Mt{I14m9Wt3t%ttYY87@^~3W@Ad
zq}@tA7Lu$3qsyVLEm3;~C7aGiU49FJ`BnWAlXXllPkDz0tZjAemSfh9vS`hqC*3S2
z@>F4uJw`WEg0&Vh*^~zt^))BzPsH-lWCk-cs#45})bbT)d!WrOlvDcHivv&OB&)yI
zQa%szXK>TOet5z+d|hbJ3+0lLzga^=GU9Iild``JurCb$joCoPGR<>lYmv-?em8To
zbF?zHv$Hz?C0r71*~UY7ck}dW*#XoQ3So_LVtuTqdw>@Bx=P+!)?CQu^gshFLkY;l
zN^dWEpY%><weLVUR7rjmt2BcRLb#;xs&SLli+}3ajIe!DR}4M+Dl_BzQMXkgO8!ME
zM{mX$1yNS}-fd|JAO3v>>CmNJ1@Q38n?0jlFZ&+2v1q}*F^+shTftE(tGGJcuSxry
zb=J2*Y6^jA5ODL!(l<yTV_I-8Hz@U%3C4YW6B_u$oR-{Mr_UWuN9m&Yg}S%|5t?6|
zC6D`A#kU&5z8#YH&c*T=(5-jT?OglHe<a=aF8`Xx_S_21e<vFMf}nG1oXlY`=ReD|
zI$WhL5J_Y!l5pT(wm%@e=X^d*P_awqA@rL*V24kGg@qXDt`xYcRLUmv+oU4afFwm%
z?p4KA9q-o*%Y7}$7$b-g2H(a%Zc5orUMz{kmjve@`L}S;#Vpkpm^8C@4ul*Rl7&5i
zd3Qj^BK#WwQnuX$-{OX&bBS8-0hSc*fUH7XD8ZjYqmq(9ccX7q4pUVCSM;qhIko4h
zS0)pgyXZF(*Sf-9jc*#gtV__dCY^T6s=gjd&RLPzWbxIo|8-O|)&%cq4NXbS3g%sF
z82freK<vx8)q(m47F9Oe73oI(he_-HLWhcu#ZL;>aUFa~ypQ{IQaWVVb|JGAlv-lL
zSzMi*MVl<}n%+-+2V7ltw+im|A&O5+(r#ndT5)WsK<F#0i`P0R90)lT$rwGv5Gs|n
zy|B8H2ZF8kdZqr^%kyMlQ=KHMebSk3w2hsXsAGT?(1R;!h%|2Dou}P4CM6$`rkgRE
z07mU}LbuZTf=+a0+c3qAs<j^B?wle-Fw-aYd0&_tZzhM@?smU!;1)SGH`OAil)GGX
zHE^+5(w#lJnsO0F9eo!4%%fo}XmMfVq}KXGAC6AuIjgmJw^k}qyh(ybWt%XYu))zh
zYl*IteoMeVtS*E$rbP&w8A9W)cg27?@ud%77~3n|P!a{echxJ$uoV>$d)2GZP~dT|
zKfFnl$A|#0Jmd&5m~TpCLa@J{`F)bav(d&-Wc%rI>YBo6_~c`IhQ0jgI;%oVw;+wz
z8sh~UlmhRof^Oe?d8asD3~M{=SnB;LAuA{E9k(Nv&WzV%@=!0}J74w`bD@BCQ>hA%
zD(zw1TKCcN*l2kmCFwiy>w0C-NvwT2%_@h~4?)imzP+?fl+E4ZmsNhzU<=`axG>~6
z6hV$d(u>M<GB<LBnW#BCTG*O7{ixbFCAl^U!j4U!aM-4)6#;#0c6@B8VP*FI7YF&=
zzSgra5wJybV={dfbKpW^()3ir2M;iJ6tqmqnXy>*P=S9ASHXha0w3G=)BJ>6N(`3T
z<8n+AAIaH>JvQbY4jDYA21HXWg?~3lK-mP$pIK3@7(kG|LHp7HKsFZofSotJ&o&6t
zLYUpg%BywJsUON94%Is)9Zoa>9Xr>4;+i&|8PRqQ+N38+c@KIF4FN@()vE<OQ?U+F
ztbpmzyGT_mYBs+}&iXVt{<JfkgGO`xRZ*k{!X8w%Pw;%!Qu)@K(R+_K$8X>+sM2|2
ze%i&;Pi{v8^+p-JFnhsX2)2(cB8sKGbz?|nsLaj1`hHbc_(_hoS%>T4Vlup)=aeMS
z9na%XYJOPXzZlj*9l1wL8tJbC$N@14y11jyrDceJ+dL~Ts7nT$Aw2UL5A7RY=oyhk
zr8_g{-YqNzcRvM`il1J&Ov`_5;?CwU8`#fYf04b|eu3*xkh8O$_F|bq4ucD`AFn@J
z&bp40-IO5dcQ*fx5)<8>^bA|09@h~;3;(z4gA5EK<PRbIGHhWgq`=7$E;LN%HZC;R
zDI=h!A%^B!IRQyV$GCT>nCuzDq<V6d#gLH$N^-W(=JvY-Y8r#~!U-iQOUj2_fC#AR
zw{RPf@XiPL-MIu0yFqRNc_|rEJR%s8ZPE>n-qnecEY53R2*I>eVVt9UXsmg<2~1wV
zd$U)QX8~!EK!TV+E*f2_C`Ab?nc{Gpi6X+7PJVh)|KY}=(8KpYMz+439qOP`nQ)Vq
zypt9svtGIE_cPfkC!tAOEJ9I4Hc31K@KoHbe15#}z(6d{5k>J_m-T{5(yVd`!Rr%B
z;!d@-E%vIi^M{A_tV_ZrYF6C)T70L;H4hrKR>u;KR8e-kwJmgBR6N*DxhaU%{T%e^
z)1zvr=GM*AF+8s+g)>qwmmoL6t<$FJawGhBflqi1Qx>ju1wL*;r8ZwyUF%mz$-J&!
zWeH;jn@<&MG<H}%gA=`&g1V2m9PM{gHc9H#E_Hup&>^*=yJg%!;5S+;6Wa-|V5O#`
zv8T@Z3MY9O@&pagz#5+DXH;)|u+HXn-B)oEA{Ji7(=BTz)$|?+IRv&Yny~3Ea1L57
z&~h6z2z#EKtZi?aDW#=|ow?jT$o(i7{??;*Lfm_IfWQFc`yDa2nm4_&O&{&M-=kG}
zcDSpsk|ViVZ&sQn?)mL>!uRw<9nH?kIWdGlCN#(aD;?;j1Hc)T_apdLFN=HSE-5xC
ziYnS#mjvqPB)Wn&msggp$CP|Q`0fn1RcolxMHB9#HBbd!wSUGtX{1-5GS#K`ZN7a-
zTl{A9ZcROaYIOlIPWdq+k8tO1-Z0K9NvCxG>LK8LZBW+@!8h6r2tUG(lKS>V`ByRg
z)PyerLdO}=99i}yrx?+jd!5nWjD_2{Zi1=qxvbP5YD>3PNvP)$t93jPQsz*3340>>
z{Qd^*wIvzq-kMUNp3K%BP|b8NSR0VqOP$W@Y$p{wRPX()Uy`hgDQO8V5?Ur&-oyOB
zuQfm6Mwq`SA~#VYl>TA5EGKsE20WX<M?xjDev^b_+D3Oe05_We{#CZuMfb`L7xDfj
zfKlp>Ai>kM$%VUMs#N!)Km`6<($oyPw{IJhZ;2-Ji7o4*nmn~Bs_5g)erq`#{mPGz
zHQ$Oqib)dt1DdEQhL@|ERhU_D7U4F1pBuAf=iZn2k=$Xi6{UD}h<mxkJ|p68Xn&52
zxCBI-$vZQ8;#%})m*~g6L!O;q+!7wQ``Aw55_TN<U87~;ab70BYuqSPo~l}=%`T)W
zwBC(34s9ffT`KMvfPR?sd0xKwf=racV1+a%e*a1y;DCn1s*WWJ!JESPOrsB`vFlr8
z8IIV2aXRVr=#BHG7bWOC<<>moR4QQN@S49tA5YzRtwV&K(+hp=)QARiM3L>aR~tWG
z#?5raH<33XQ<!EuIP?JlC<;c;9dEU1_i983@xrb^u2(3#K_7L@aDh!<%<s(RFav07
zK_SCdy!b)Sd+%B8k4@F*9K9AI_qbwPuy?}f7?;mj03zOB`6d&;9$`E*`pP_hy1ljA
z0|Wm$8k1d{{`H6GGnjN@@+a|K+zp;HX<U4{b!*|K9tw-CqQ&sepkxR=MKCHQ^=cc%
z+N!DA=1O?Ny~b@`FB_W$Oi;eq96Ca83`5}d^QvB=lM(bj`>RKHWAh2A?sJYjbCPOS
zGEH)%tCvB)#^D`bEY3*m6c9^jEVhXWDz16T?sB}#H1S$oyhmqM@$CyD1nf;`0nPv)
z(_|hEK=Bp_<}lWqt#;vT-sXrFW!uISpuJ3>9@I(3ovTJmQ%_E7xA>)ARZ$g7)(A?U
zQoI-6z!wWo+3VF&@ykaA9I4m@?@l9xU8cESC-CFo!yg%uAG?QmMQLaZbdN36nYQ`t
zvTn~{Nk&1Xl!^=TG@V`b4DY~`OEe9uA(=wrTm0J5P<Ab?+{Q><8&`MqnCT!(fm~G+
z?stlrCkXO`7#(Jr`xGDM)j_=1Id?PHtEft6A2HR>ARvlkWpa4c`ka;e)zFPvSp8JK
z=#>|pv>f#xK6)RTrS~s9Lo6S@?eMN)sVqI@c&!b^(#^Mj<X})jYVZ+^&awMSb6Kdc
zPqBEksOp$0Z^|P6`GaD$mGOp0)zkYQ?^S&7<d0H&HB$0Ov3P=ItUFE1(b842d5B{y
z6_T?)Iu)5V;=^ApP~^2;&tuoUI@x?V<1U65hKSz-SLD)+h8$pc|5-hpJ*;6rc8;8Q
zedMN08oYVH&M~cqhQ|@+XK4|xo}bJn5a+t;CtGxcc5Ng^;doO7lbG0LG8fVPt~5G*
z_vl9La<&9_FTY+PDVEmld-7y;UoD@_mL`8pemJUK0AYJ=%zFD3=Clb->|1Jd&&eC6
zM1<;bvyGMevGCmYH)NC9^>1J>JQa+EE<LBRuxvB`P$R1fA3S(S&#L84%K_UKCzS1T
z()A2V5UBL(!Ke;N3yR@LI1tIf`*hWQ65|FmgDX_Qy=YLbG+gEZ@Y7QF1JwAXf>{jD
zb(cHjf#f9Euk08wh_>5|W^`07U|{tnO`Z4N86qJ)W;xA;AMtKFbaux|;@H50RC~M7
zo!M%IGbq?pDQ{&8wPN0Jd>{Vcgcoo}{Ql6xR!%Z6s~VCn<xD%!7e_k95=l{j&K0n2
zlFCizL7*MJK%$Dh0H%EULG)T`W(e(*e&#*A(NC6Ytr{taITl!$;nN`O;7vqHWgWBa
zAoTqb2kXsG7Scgd)Migm4RM&;<LXPN3vgVu=*GJBCXe|#?6*IyMAHF5ZKyY#DZi%(
zzM77d)_3U=&pei-*ozDQD4Tm2;R!aD%&_SRnqMIZeyrmfH{ml{_s&rV1IScp&c3$5
zpE%2U)(iREy%oRTY4Dn=BKPo3sybmksHlP<a$R>6ls230Ru~<zuTTBR$yeaObb@1t
ze8_TLx?*bVQ7ZSzVDqU9!g6}s^j!!h9L7M~^CW_*Ph0^SK5#4G)aQE(qV_-<c$}Wq
za<F41zGcdr>bZM(YiEZAFtKgzLNPy-bg&fJ9H``P$VbD;`ApuvgKOQcUqhI2=?q5n
z)t%VnsHtN5vOxcR2{W;KsO)HG|L4PB*L)i+5V=A#BOk{Jk#+pj%*fvUe0@%hy<<1U
zgB!GYOcuNmnPSTlETPGp4Su9uh9WyP*V5(8p`L&qBM^If)P+BLYgMscES6v(owZ;}
zYB8Le<ii*CGQIdL2|G!f(<b8RdONLj={vYhm=vO#v!RT(`&gB-bjHRJ>R-{7o|-4k
zK_Gmbt`ZTV2#?{wQPf4gb(lHGYKtpz0$;3WwY(w5-EdD}T+87qX<oa=z;`ed>&Dj5
z$Oh0UP1C^Gc^qXieUxB4TOxC2cJ)n$n&WNto8}P<0y?)>RJ%O#t1u}`oL<?s5V3t5
zS3Ur=^$Qx3Sd8r!+;5{JGx1$4){mH_AjB7bx<se5#i3ZG#%oCt7UzS%90a|~%#^pl
zviP~SDt*t9-^9Tpo<wH+S*S=4{__2b$L-%~AM}Ua$K-wZ8eRM`AZn9vfaaLphpG;6
zdsA6)1yyv-VM}qPexG{}=r)579FS|WF0jGDUOnL&VUWjjNw<gtpRAIfP_I@jo*D{=
z<&At;o$t#!(N36SO7Vm=>DQd0T_*OT^W4i9@iRbXY!2Cdh<?x5Lpw*<kCdtQoGXA0
z2ki9I0VLF9;gaylaJ3k=rv(~zC5?V^#=3}McysA4E$UM>I}wqJId5V8v~bZw)YWWF
z9O7ttV!v6aOr<Wox{oa-u7ys`e;)eAU;RrbU%gzdY)#dcaV2U5vEG-+c)1!KxLB4f
zi;EWlh7&`na@t2rU7o$cjh99zEv^#pRGgJ?B$$DGP;hXr8)ItDm&Tl%lMTFtDU$dY
zEYAGY9nORF?^{1{m?z@E=*sx))gCc>?ss)P&`y|tnZmC+?JF-nK2X~CGE;rOKp4(J
zCOSh7v`lWy4;apYEK8Rxs2Jnx#BE5@)Xq%^*;|i@ZIz+CE9PI?9F5h&_KiYQ&#7wl
zKwKSrW(*LqDNF#tv2HrmAQnGE2nogpR}I+KU8nj{=(j{<C-Ksf4!Y(>Or^^m_|t+N
zNbglBh$@;^nB4ka<k?g}gB!6C<13PeSNE;t42Hn@XEsbZrF`^|*>FZ?gWzHg{H{So
z#wsWwb-v)ulW>-UV$J~6u}N6~-1QEnAN!7~TSNkP$FlZ!G28nwgcMv3YwghnB09xN
z<Dm@qhCgeRCoZaHs5K4QYIRjl<b5b6jc*Y!Y-ocZWE}Mkif((0S~rRV`rXvKydjyr
ze5#(WC7fsmw`Esf1hw0KgptxqN_n=MD7Ypq_VoI@?#Q#u`2w44DqFEiZ$)sW(h=b?
zYAShBnP=(j)1(BsFSW{hfPrOW_wsxKAx0XtA8yhZS=FqnIf%u09Anv3hNpQYFFZLS
z|CTpRKM+3q(Xr98^3B2SC`RBGuIo$VY#xc_vLWnbh@--s__h|U>rB%~aB97D6I|zX
zmzZ0H|GlKyFni(irIXyA!UmS?FR|Eja-;W^#X>?ps!)0B@|Gl<3X|BXXJg;?0~~FW
zFh)@dSWsi*unJ5+NnWr-fsx2QRw**!VOY-X<hFYYj_W@{CIa;eDawBr$C26i?fi>W
z{>#DNE%RK{y*TGxF&vrUOC9o`LI2cvFHtBcB|kAggvkFd*tryXan6IaAWP|oFnSs3
z@7tYAqZj8KFmnO;_Y3|hm0kjoe<>itBI8`ZelObpdEawI@#36^XGN+nNI(D37cZ}V
z(X#$j7%$N%C>noZf9Z|C7`yQDr`q$&3jx-T8~mmH{POU>D9i=?@Ng+p{PX&M3KhQs
z_x%q0C0qQ8`|r!*F9+eso`K}f>F4VB3;MsV%uf8<wEhLVi1_DN{}oa7ci4Z1{GaGQ
xN6arYI#Oo)xguN+oIl}z_QYTCShAm`a;efQ%VGTZtq~&EkGIH*1XG-w{~xp+X?Oqt

literal 0
HcmV?d00001

diff --git a/assays/GeneExpressionComparison/protocols/.gitkeep b/assays/GeneExpressionComparison/protocols/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/GeneExpressionComparison/protocols/ComparisonOfGeneExpressionAcrossGeneSets.md b/assays/GeneExpressionComparison/protocols/ComparisonOfGeneExpressionAcrossGeneSets.md
new file mode 100644
index 0000000..6fb75a5
--- /dev/null
+++ b/assays/GeneExpressionComparison/protocols/ComparisonOfGeneExpressionAcrossGeneSets.md
@@ -0,0 +1,7 @@
+## Comparison of gene expression across gene sets
+
+To compare the mean relative expression between *R*-genes (*R*-gene set size for tomato = 359 and for potato = 581) and non-*R*-genes (the rest of the genome) we generated 100 replicate datasets for each transcriptome by sampling the TPM values of 359 random genes from tomato and 581 random genes from potato. The average TPM of all expressed genes was calculated for each replicate dataset. To compare expression values, four reference genes were used: ubiquitin (*Solyc09g018730.4.1*) and actin4 (*Solyc04g011500.3.1*) for tomato (Müller et al., 2015) an importin subunit (*PGSC0003DMG400007289*) and elongation factor-1 (*PGSC0003DMG400023270*) for potato (Mariot et al., 2015; Tang et al., 2017). TPM values were tested for normality using the Anderson-Darling (>5000 data points; Thode, 2002) or Shapiro test (< 5000 data points; Shapiro and Wilk, 1965) and for equal variances using the test from Kendall (1938). Significant differences in expression were identified using a Mann-Whitney-U test (Mann and Whitney, 1947) for non-normally distributed data or a two-sample t-test for normally distributed data.
+
+We visualized *R*-gene expression using heatmaps created in R (v. 3.6.1). Genes were classified as off (if TPM < 1) or on (if TPM ≥1). In the heatmaps, libraries were clustered by similarity in patterns of expression between libraries and *R*-genes were sorted by the number of libraries expressing the corresponding gene. Correlations between 1) the total number of expressed *R*-genes and the total number of expressed genes, 2) the total number of expressed genes and the number of pseudo-aligned reads, as well as 3) the number of libraries in which an *R*-gene was expressed and the average level of expression of each *R*-gene were performed using a Spearman’s rank correlation test (Hollander et al., 2013).
+
+To investigate the extent to which expression patterns of *R*-genes were similar to wild close relatives of tomatoes, we evaluated additional transcriptomes of four wild tomato species: *S. peruvianum*, *S. chilense*, *S. ochranthum*, and *S. lycopersicoides* (Beddows et al., 2017). A subset of *R*-genes was further analyzed for their patterns of sequence variation within and between these wild species. Standard population genetic parameters including intraspecific variation (Ï€) and interspecific divergence (K) were estimated using DNaSP v. 5.10 (Librado and Rozas, 2009).
\ No newline at end of file
-- 
GitLab