From f2897f587966385f998418a0f2f5920c783c700b Mon Sep 17 00:00:00 2001 From: Viktoria Petrova <vipet103@hhu.de> Date: Sat, 5 Oct 2024 14:52:12 +0200 Subject: [PATCH] add calculation of transcript abundance assay and protocol --- .../CalculationOfTranscriptAbundance/README.md | 0 .../dataset/.gitkeep | 0 .../isa.assay.xlsx | Bin 0 -> 6910 bytes .../protocols/.gitkeep | 0 ...culationOfTranscriptAbundanceUsingKallisto.md | 3 +++ 5 files changed, 3 insertions(+) create mode 100644 assays/CalculationOfTranscriptAbundance/README.md create mode 100644 assays/CalculationOfTranscriptAbundance/dataset/.gitkeep create mode 100644 assays/CalculationOfTranscriptAbundance/isa.assay.xlsx create mode 100644 assays/CalculationOfTranscriptAbundance/protocols/.gitkeep create mode 100644 assays/CalculationOfTranscriptAbundance/protocols/CalculationOfTranscriptAbundanceUsingKallisto.md diff --git a/assays/CalculationOfTranscriptAbundance/README.md b/assays/CalculationOfTranscriptAbundance/README.md new file mode 100644 index 0000000..e69de29 diff --git a/assays/CalculationOfTranscriptAbundance/dataset/.gitkeep b/assays/CalculationOfTranscriptAbundance/dataset/.gitkeep new file mode 100644 index 0000000..e69de29 diff --git a/assays/CalculationOfTranscriptAbundance/isa.assay.xlsx b/assays/CalculationOfTranscriptAbundance/isa.assay.xlsx new file mode 100644 index 0000000000000000000000000000000000000000..abd79d7d4d2ce2af98e99be7080b6ce1a10f3bd8 GIT binary patch literal 6910 zcmai3by(Ex(gvwT8dg~8mRM3!x)he~mXM{pk(7{bM38Rjk`9rO?h;TDkOt`m3H=s1 zsIR`~JI^0`?X^4i%yZo{J2Ss~m1L2SZz3Qd01zZirL<$`+R(KT5fFBf5fBLApXy53 z**ZgQoek7H?4eHjtnN0}r3v!3?d-Tg$3jw{Ss<ADYU)z$*{}Ku3l9;--NnJAbNflJ z=iKlaikK*i3<;Q%d3QY1J8vu2M#@Zqqy<vKm3e^Zvtx`lL=7jYuk)j+vj7@T90;8% zCGg`HZ3X8&asyLXmb0TqMPy{4kyuJAf;wettX#NXdBC~Bj_4Q(#~lj+;Iv5dN5r}~ z$`J)xZ~8hvL;ItT#cmMe_Z~caSLKNKkwjnobEBT_QE(hN;I^A9aZybHB1WK~j5+fv zVsNWly4=Xk_}hUes943Hr+OkW`q@?!XQ1%}jQZ;9_|KR%vFc1dBl%|&bMVk_PL&|q zpk16R$W`gy=f!qB-08*p>b>;sv3KB5KAN<HS&Xd^Z0J?Ao`mQT9uh@F3fjrG_d>kB zfz0wZ3qLH?A)AD}>Q91ucvvUj^}nUQsTDY~4y4-X&k@3(Io+r<7GMtAK-P5?DOn^@ z`cEbheq;hp{?CI7e`nwbwRU2=wz<mIeum8KI{Z7?8wdy(KW$hqElcCxD0Hx6itczx zeLj0(D-FoXYog5kK+h%IMr%t_6PeQ)@MR}7MievB0xim6Yt^iOb5JmHBT@xMiU}?z z2BP#21rN3Lf3tAFVZuNSra|Mh&&5nxF};wZffV+(L~GxY?m>MVdy?r=D<q%990lxC z?6(kS9G+MYaczq_EGpY^KJE2e2`sD~5}j?NdwIbz#%*n@YquV|ZIn%ErgPrMU?N8r z`ov>mC-tf3N*0sS$f~}^bkn&|LAvxvR%UIgIi6~v!eYP9YZuZvee~7g2wCyE?+v6c zg8XURw9y}(bB^8N9q~fAX5??ykl>8CTmPi&Zv*rzgMVW-khW6elGz40v&6rfx!E~d z8Qa-eUH%fTiMH%u$Ccc<xLJ9u<EjH;jCEpsqNjU|6!@xE&RWKt*XH6_okW@xkcF1f zS@JHqTY7Wwm@h<8ZUU`5lL>;mCSTFA1N7pWJ2S)WoYfUVj=o7x`)<N*lb4ihmCVtb zHdbDc(SBf03c`u~Kwc_j?LeMn>?MEyMDNSNhi(j-(66-P?~yhz6f3K5jtyy0zF=JR zZI+mWp&A76N34B?1k$EI9bh9)yJv#(K;MKSW_n3e_N~*$E~nEBLF{5(Oq_6yPtM|J zLyRK3&7ogU2nUwp*bS(*d#Uzsf95)s8tg8-?XkDCLGj;-#=auxk{Ty-DAf7SGOY_! zZVZGI*$yWh^q1`q2*)L#&k~jGQrL0*=8u_U=1KT?X{l}$xhmJlq;T1!!L~pIB{%NZ z#@C)5HSx)QF3X(6iRB{^#y)FJJxEzCi^LWuDLnOW15?GWH58e=X6_jdJ}bTz8Ugk0 z(wPkRZw5%%_ThYuAB)b%YrhXzlNTel3U(nS`4|$FoJ=ekeYa+etQxeTZ;i^LwM@P- zn?&DB{UvFuH?(4E$LM8aqMkKy-YvWCP6#1ObyBOvXTzaaQLoXaIWFoc%IY^zC9R>% z+u;FmFPAoln;u%!+U!+l7!92yZ~OC}C^(iz6m4TV_!4-Z4Qi)$Ni!Wl7KunTg~qa3 zdss?#7-H%NJ`5hS_S)Smk{pDUUX-N^qc>QAzbHefYwAk3x`-TbSrl&3dI-UkYitLg zjb#sc+U@nq{d3osZ;_a4CtDo>7y6LC?6*ao0jxkCtjVLm_?2!?yFFA=P7p;OZ8Q#) z-07TZqdkX8aAVIf)s3vZ3Fhvc%8RowDD?4&j}2=fkK9gpNH=hokerQdl|{l`Ho6Y9 zS}N|&99>7c3MG&J7X93#c`|5q<;!`4^|^iw@-5Fr%~i<;i6oI$QFi4$++5sdNAv78 zsvhcHZvW86V9MAwUUYg0g}>en1Nx+wK7cWFuM9(R1OVqvuROzcL_pk4uVO>)Cj<U5 zt%B@EI9OG|r?8PiQ#=!#qwTEkvjm<kHiiOw&(_nn<R`*rpV-qJ7DhK(6{ETZslQU6 zD*8go{jN4h`2I_=(gY#2z0fm>cNe&fEF5BPrwl!ruO{Vmyg=Q)%&F$Q0iC82)gHB4 zW0(!@6IF51av)M*x5yp6N}XA>BUz1FhqRfX=P=&^%2v|XeInOYe$`+%ki?j+;p0#O zJ`RCbmFr|~<OnrUb#}C{HFNq=wFw5Y9qhoa`3oS`4yzLm^_VMATg?7m<LRd*{l|5v z`voicRo<TCRM};~vm{U7JN&&uHdTVegTf>6jy%h})%<4WdBx9!a?fj+%69m25Nc7k z!TINbozG{~#erH0IuH*6aZfhnUR11+fg1BiW;@MKaG!spe0c=m<_>vdUYRs(7la7q zk1wmHS`J4sOH<z=6ZRsi3RFuB1kc+Z)JX<{@w^v`rIGK{-p?G$5hKxZvE1_o7LSl5 z%7#uDyby9T4fY9h-Mq&Z`>fblOFhfFWc-Ugape)(i$zPNdv7N0ORRO@#UzB#c%nKj zkeu>_qK#CA3rId^JL?OeTURI9*C&!Nrn@P{E>K;alfTqObs;!os!qtZ{8-&KSl-&0 zdWl(3EMg|+xq({Nxo<^H(M*oI*pqx^`&%7_4?^V$M4mb*+It*iTr`(?lOFL^&oh?} zu?h7(Q&>7tFf)d9vbsoo#ouF#@xF1Lo&Va)oz0;((4Q^;Dv5DJ+;<}2Beaw9YNvtE zge$WjuRl7^wzi_(91rk2lmC}86W#rcOk1OV*Kr;T|F_#CG&JLck04yqOrgp^(Cj!X z5~_0tE0XJ+5lF)jMPsW9haj_S$~#0z<{QndT1u_OsF4FgO0Lh&UiUFMg+XWWw4#J1 z=_6J^xQ^-9FdJgN{h1iar9=<A5jJi)327h}9u&_u`7U^1bGj^><@P6B5=ydAmI+QI z#sb|$Ixo=u#hWRMfb>WZPHZ46g|0-Df~b{rX_(D)32tl;7d6m-tYtLh<Q<5XX)te} zJg8hc%%rX0yiL(;KsNW?LT>7LNb)WNZxo(QGW&2$8s=^x7gktcAR5cKf=Is0c2Ny5 zyGoSj&a|S4Q$s_Wy^758$%#GV8ef^J72A;}=S52W!xqiW$;49?gne%<3+<fhhkL2~ zJZOC{h(COIT&JV4%YQM6<uxb&4d~?(<i@jm(OOq!gq^_s0jqh=!nLu;$1SMb=F_HY z)8@o2ubVd+Lg`7&=StOEx-6f^;N{HexR1M>9`%v7ifcEl^?jz%23k?wGj7K5n`n@Z z>xrpmB&VXVC(r&IL+~g#0twK}7?$K`)MR|T&E$2*S78<+6jsaLCu1hj`VIs+0kyB1 zFzK$aj99KvvKch<d7htd?d_N;rl$&ha}hqyf6o*4)}vus#QR_v$AH-PJ8WsQV18qd zI@<TBU$gw%$$|Vvp7>^yS$VpM=husA--~lKBs(YPq+lHBkRS)N43L{P0AoVV59e!> z4CW0<AUY611!=2S6!Bv+RZ)jag=PC0DQ6J2JB_ePJvp*qq9jr?S>R3k=Nz*}dQGWw zy?S4lJ4dxd_@fW%n*d~+E3hfj_u&P&`;rA?7!~498UA&npa)vSy?1%uXwks@aJ$Ny zI#=Z?V!6n1a{@x9XptNlj>P9^kzWruBfpspvvJ)aA-nIg(R89E)m|&AR)DYC6~U_n zR(=VM5P$LD3+3%KY4U;ka*&?%?jcCk^zf+`AZ>s=gVEVe;^}CU_wym~TfKD2YcT>L zm4a3M^bh^o3lr~#`U}GHlSD(PA7#j}pyz*y$;I&zRnBVKApp<Y=*|aV=HkSBmKku- zy>Zt?WM~awl(x@<^K5H&MUsRp&AlWLhW!?pmPz&YZA;2M!4yuRbzMZ0XEr6(gDkmk zEytoO{5Tm4t+=A-#L;Jv1Wi%AT+OUP&7NlC?okiA(OdQ$eo7e69~0V8Oi+WkS6Lj< z!X!h6@?1njAzF0aS<%zCqd&StKN%SH?D^!D_@vXvb`BG_>(uWyB?G(jI-#U-i*!|* zN~IPvuL|#WAJ&vk3qjmkY1gpM%#zQGs@0rZf<y)zz`TT`8wG%43b0ipLllf7mG-&% zAXNRpx5P3Gwhv`-((c!r;>;*X)PBaM@r^~Xhz{(ve3v?byyte805yvj^45hB1?sp0 z(<`qIF04%c423t5{E#_RvwaNefB*z}qZf|%I<yAVqk}l0Hz0SamE3e5cP+4znC2|+ zFXqt$C>w}_$E-N8gI)~Ww>p}fYsx!)#Y^aM!?x(~oYpbEkhX{z|7hc@bi#JH@#sW_ zdBS{Wd!Gji_8lZTyAJ(3kB}Eosf6Ut6MEU2Jr~khIrAI0!c0BnSK9?kV|s#8Ak;)p z5lP86JBYS6O;vX`!V>Sd>~VP6*sP!u7YZ#Q!}KOmxP@QT4&a@SBM;i&Je7<q#3g&c zGXC62;<ch_vLjWKH1ch*cS5NMExuDg9Imm@4k~eJ{WE5lvje*6S0W<)+M5b*bMRo$ zH$6od!<=-p1rz{<dtlTtv^Tq*e7PL2!`qZ>TT(&x(t&z9PSWnI^_m)bvYH2_FZF6m zY8kS}5e5|#yf}wHS$N9asf$WjKP>{Mq2qL4g!8$~v%X5?!oZGsY(#kG9_$sRu0Gs1 zxzcFb;d8*aw}2)dr6ZwOT2!Fn?5bxd7L!t@VOR&r;uYED(u&q$*3`^ziPW`mbw`e! z53=OWS21DhR>(Ss5gx~C(@Q@fnpsvO=D5RhkhNV)R=)U{u4w@VQJAch#j4Y1snM^~ z`O*Mwnky8%k<&v7Zkl=TePWg|wDKIbe)6`<yPlz@`~>_;O9xH2(EhQ5K^f5CJqa@S zph9Dvw|G#YbfTp8jILnLBH_iuQq_&A=ErsONAK@ffA8UnQmq&-d#q48%`n-QuIXs$ zD*k#DJedZ`+n$(<Odt2*s^c#4+G}FB>)V`teLdr@h8K!}-$PgU?u-N<U|9dzK%70S zp+8oSyaavt!VG-6bIc5$S4F}Chx%DsgsBy#FmcDbZu-fToFd&GkCi{$5kSSqcbUzH z^>vp=XB?c~ZCKA0WgFnqD+Z!z3g4Hz)%e-+`C@s>`;<o$T1611SH_IO6{zzj6mf6K zkv(Vcn&RQA#V@wh9L2@tzq>1w!mNK6h2|MgoX*+{G7HNN^O<@Xm6(y^N7Rg(?v!BY zo(QhYpp&j=P$GAYS3gQ!aC%TIIPq8@59`BC`&pE`I+?5?^6n)gvgKjY4?!Q+`W_-C ztQ9Sycy7Ch5e5<xpjX(@py2Iw7%gb4SU^eC)-<%=d1nd)_nYOtE`E>2@6gj1Cyrqg z6QnZGi|owQz?VtHq(XWxi?<zB%<)~=!*dRRn8<_CN9`=PysYX8dX=*5h(9?}DU^u| z0<>?C*d{CUJCA_uu(^|z?YU9qGLECS(z1dnBZlY?u_iuPs<x}A!j>4Ip@z?b(4X$W z%4!<vZAWz8t$`W&KUhcwNsyaGAR1!Ocg8oB&lh31YEn)1>CK*TcG>TJ*odYA5qBWo zbte6u%2P2PDW&hyE0T34PIMR__Fg9cB;1q4SUl6FKWKRa=jjt|*Z66liN<b6Z4?k) zu{rbB3Rlu1<F^6G$G+W!qaK4-WYzg6Z_?Con}|!QaU!>MCy3J*Gu(=!!;kdIA3OPS zADd2t_X$TWx23A*CLgD<osYb}aDiFQ@0oT7qsBmKDElMA$p%H_b;5@41zh-iZ-X@) zOM%WZvfGaLtweTBInq22PVVjRGXSRdtX+tfN0X1&B3}n8`Wtdmu&_Lrv+rWvwi{CC zqh0$3#ry1zZ*tmNJ%3%G|GtEo*gaBmw6p*78L)S`0~!e5r0L;L<+$)V{%K}pZ-2Q( zr^Sic&9P$!?VR0u`Xw^emf@+W21hQ*W35UAnYpF5UT3gcB62Ku+{I}x_Ts%wg(jgm zoZ$?{qB)7xFmi&KPt28i3E84{;x-qp_|Z*vni*1Jn60Qpf*Ok<w6;fRH8NDj#^Gw8 zkrkhrCoe%DoGh-Q;es%avA_w$RnBdwIkD9)Ytl5fQ2%CCbF91JA@`J~!%bj8r~2@B z5*@U=yQAY@Ku+lz2FA`)2&)+rI8(WzSqqDsZ!%RKh1K}Y!&kVqg*Q}sJql}4Nz0rn zY}@ddzD_9}13HFy3<)eI4~iaiP~9@|T`koQUnIiC=6kkArM(MQs8!{#Bnpl9!Cnv2 zk))?9SYcTG*if5s=*VT_V39x|J@q_9pdWkvLG_c)@01URLLZ=VJbHyJ@&pjIgF8%d z#_U7Z2oT;;QrJKg+;Z4e*l0RpI|R8cV1tHbTdj+1Fwi&8S;uMQuv{`M;z`ao3D3zl zt5+`!`9cfEXEv7yv(L2>m*`SGA+7rL-;k~od)0Xor915n;2B$jcOSgpGxo^N5&9!# zD*cy2VDm9E^?U#UgE~nNWzDSuYIUZU>UL!<ezL~8u(6oe<p-3w&s6OM1ge+3`O+xD z1y2w+b2aYbM^oebE$T?uq{SXt;(d&7qf+%>)_D`5_Nhm(Nv1)jzINBL1~D98?^9%g zY<+SpS+)#=ix(P-6HU2N`g?O-_QMg|LL;x?TG4K0=LZ;~^dL?hlBfA@v}t*tTJrec zHFMymiesZNICGJA9hOK{1U2c5uoLx(8$&-89ZJGBa^hPm!2Qbn7P+n8lYGDiA73a- z6mF<HP--)|x3jQM@aPt|jO9VrrOH;6jj^@kzlc*bbWHQwTaOFvRw8tla;@!5#03P3 zbdwjU=q?Bfh#uvE+BN2i@`ic>Qu*~xo$P2sx;Nf$8Zd3U&J7_@?+VDwVx=b?_l_S} zid8vqwPGzXMC2E{r^@BfK0Yb&B%jE@3jY%8E0B)a_!ad#=NZaBvti02;iCu7hBG`H zI9GGvcPT0|PF@i%_dVS?4`VnkWeGr>oRtB@xQdbbG4HFmg(tFgt!sT3vV9PXOT_B1 z)fsIdpk1ms6+&}=?4x>B(yB_PYU`-2W^dhe!AvPIp^e+HxdQ@|cGNd05%w0eZV>?t zxvBMfL$U@qRXksbI#G=5$!z8Xb=tm%0;$C%JUdO~U6WV)2mD?4<rtQHL9I45?P%rq z!dcU(@UUn#6g|nzv$c;XQiI&r+T}b*K$VmC3w#10M(PbS{1ir3^_!{=La`oaXf`!r z>0T)-5vPP-3+AbZ!xrB=wpiA@IX;*`3EahWeQBJ_F1lVhik<>-lwT6r)1-7=XdQo= z)+E&$qkVCJ&!)`vPTXvax%kD}d47L!GehpDIP@jii3dtT!NKp9$-H$r%2G`E2yE4I z(S`j0r+WmnQKZ}!<mecT-18ABE0zdQ0+}buB_`}N>-jxwc5g}IhfXn+WD#!w5&pYI z4$sDK=U)`_Uk?6mnU~V;)w%47VekxJ%aQ*K`lr-;jY2>u`-%CXNB)1oE)~(Ma~Z4+ zUP?do(d$Tm-|kW&y*ihGg)6|nU+_=G^co2NX8|4-9_I@7d(r;S`(6r*SLZT3BV2@m z`}v2=czyM&mi4F5c#TFtQ2z`2OLqLl*p-(*MW0_@aL|6-;4ca2mxupFVJYB;hijSQ zpV$9WsQ49l;&<3D+2U8+e_s}VIf#Y#3^;dAKUc?J(EoL1PW<1d^)J{}#6QRSuZY^e z!~Qem|3v>eVt%19;A-2?72$f|{0aZFC;oyb-}+f9*J8brEXt4H8ZLbOcnhycYNAW? F{{uAhR3QKW literal 0 HcmV?d00001 diff --git a/assays/CalculationOfTranscriptAbundance/protocols/.gitkeep b/assays/CalculationOfTranscriptAbundance/protocols/.gitkeep new file mode 100644 index 0000000..e69de29 diff --git a/assays/CalculationOfTranscriptAbundance/protocols/CalculationOfTranscriptAbundanceUsingKallisto.md b/assays/CalculationOfTranscriptAbundance/protocols/CalculationOfTranscriptAbundanceUsingKallisto.md new file mode 100644 index 0000000..1b55321 --- /dev/null +++ b/assays/CalculationOfTranscriptAbundance/protocols/CalculationOfTranscriptAbundanceUsingKallisto.md @@ -0,0 +1,3 @@ +## Calculation of transcript abundance using Kallisto + +The program Kallisto (v.0.46.0) was used to estimate the relative expression of genes in tomato and potato (Bray et al., 2016). As a first step, the raw sequence reads were compared to the transcript sequences. This step in Kallisto is designated as the pseudoalignment step. To improve the quality of the pseudoalignment, low-quality reads and adapters were removed from the transcriptomes using Trimmomatic under the following settings: seed mismatch = 2; palindrome clip threshold = 30; simple clip threshold = 10; LEADING = 3; TRAILING = 3; SLIDINGWINDOW= 4:15; MINLEN =36 (Bolger et al., 2014; Figure S2). Subsequent quality controls were performed using FastQC (Andrews, 2010). As Kallisto requires information on fragment length for single-end sequenced transcriptomes, the fragment length denoted by the authors was used. If this information was not available, the recommended fragment length of the reported RNA isolation kit was used. The standard deviation was set to ±17.5 bp. Kallisto indices (used for generating the pseudoalignments) were based on the tomato ITAG4.0 and the potato PGSC_DM_v4.03 genome releases. *R*-genes missing from the current genome releases were manually added to the list of transcripts (indices in Kallisto). Transcript abundance was calculated as transcripts per million (TPM; Wagner et al., 2012). We chose to use TPM since it normalizes the transcript abundance for gene length and library size, making TPM values comparable across experiments. Genes for which the TPM values were less than 1 were treated as “off†and for these genes, TPM was set to zero. All scripts and settings used for these analyses are available at the following site: https://github.com/LauraERose/LargeScaleTranscriptomeAnalysis. \ No newline at end of file -- GitLab