From 11238617acfcd341685557373e42325cc80d947d Mon Sep 17 00:00:00 2001
From: Viktoria Petrova <vipet103@hhu.de>
Date: Sun, 10 Nov 2024 12:14:19 +0100
Subject: [PATCH] add amplicon sequencing data analysis assay and protocol
---
assays/AmpliconSequencingDataAnalysis/README.md | 0
.../dataset/.gitkeep | 0
.../isa.assay.xlsx | Bin 0 -> 6907 bytes
.../protocols/.gitkeep | 0
.../AmpliconSequencingDataAnalysisProtocol.md | 3 +++
5 files changed, 3 insertions(+)
create mode 100644 assays/AmpliconSequencingDataAnalysis/README.md
create mode 100644 assays/AmpliconSequencingDataAnalysis/dataset/.gitkeep
create mode 100644 assays/AmpliconSequencingDataAnalysis/isa.assay.xlsx
create mode 100644 assays/AmpliconSequencingDataAnalysis/protocols/.gitkeep
create mode 100644 assays/AmpliconSequencingDataAnalysis/protocols/AmpliconSequencingDataAnalysisProtocol.md
diff --git a/assays/AmpliconSequencingDataAnalysis/README.md b/assays/AmpliconSequencingDataAnalysis/README.md
new file mode 100644
index 0000000..e69de29
diff --git a/assays/AmpliconSequencingDataAnalysis/dataset/.gitkeep b/assays/AmpliconSequencingDataAnalysis/dataset/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/AmpliconSequencingDataAnalysis/isa.assay.xlsx b/assays/AmpliconSequencingDataAnalysis/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..419398f43d7d6f9ede46238daf8a729f41599c39
GIT binary patch
literal 6907
zcmai3by(DE(*~)fmqxloV(FF?5ZR@>o29!EP^7y>kVZmUS{ezZyHSyDDOr%vZ;=;0
z$8+BAeg4>Muid$4p6j03nfcwTEQf@Qfq;MjK+w{$(TSgLN7q3_K-fn{Kp=pBswd^>
z;0ASYGg9|<g1Q>8d)eESB`G*`Z~}u*L?5iML9h(eH6C>2*7Oq>eMK1ek_1uAAEva;
zd*U+{v(OY96R@TV?D=GN-BzlL22YX73a3S?@DrlXjWOGkG@Pcl6~xlz05qRE6T4PQ
z;U_LRh%9*L1*NgA=EjVQgTc^fY~?i(-EwtyUf?%=P+q7DI%d*I=VBlzBihQ0Ob=Hj
zs!;oFU)Kt>KlbGQEi(MxBg2_07sOBG29hg{`g+Hq2~>dFo*rbywS|b7K_Xx))^)_t
zR?kfN5sbv!L8qwLB`Z@s(U|=l>&dgwL;_|5jZOS#tXkOh<|{}6xnx}YjN4PCi1uig
z7mD&#dNKkyE?@V0@xJ*je}Ci`G*o~l>uee4APO6*iPe{qIL1Swj!Hv2-St~cG%y0M
zez*3=RvWTU%CG(`a)5_@`orLR`rA6;W4j>Qt^Svy__JqQm8Qb1!CT0B9^$1-<jVia
z1j5fuz{&r0P~q>4T%dNY95*)C**eSu&uzlLle>k0fceXY{mQZ|;jLmPCziyX?}L@|
z7Y?$3oct!5ypN2$V(knL<h9W+y8^%Lg~v%?MO&lAIPa`m_HPf0L~li_!YHsnC1ivs
z{X?Nc?fu`aopD((QA6p`xSjH_($*|4<>?_sy)CgiM6x}okK#|W-RndZQdnaM`;_{v
zC7FjOHp4vHW4;!b@421z`mY5Q)eT9^HQsr3$u-7j=b-1f8NX|iOJk{f(Z^&iPZ|Ey
zdtxs=L~AXFMR{c1Ky$k3LbNbbb|fdeF5L=Gtw?dHU$@PjV%`9KeK<-^vi?T{#f#to
zdQTm6vkUIAy8<J=2se!UlQkqbBVKmDDEp@Y`nAEoF&o5Csd>e01Dsj1znXbEy4adJ
zI@(_S5^jjL>Er}T?_FY4p6Gh$LYU)SnV;(Gogf9()XCd{tpw~ZPc+D7DF8WWSzV<w
zsqbaC2TugUl;kJSDzaH1zzv1hEqjE%yz}Rlz^*wxQRG+*MuwRQ&usw;-gQbBKZbY(
z5oV`>g9i|9oQDbz!ZwZ+$j4p@^-uJ^8Z`7|(t>_t82^B@g{f3ojWIT)N%Mkv=}EKH
zJPg$+P$+8S8zhJ!Gh~2+EaRRz%0mNl>bU7;ExC8DpSoSovP5u7^ssOvH9xyao)0mL
z?>C2kJ0%`iPT(}6-R-44yuHGE_F(XR(QWU8<t^&}PBhLnL08ncT0x<1f0k)|gi2!&
zoX8G1;iUht{R!c^;`3Ruienll(0}2CHEw}iP=JB<R<ValEjW$WJ_EKxN>F-BrY^DW
z{J2R_ZlyeX5;tCuTny*DHT@`Uy*wI6lDz0FpdCaTztK=^-p1N99C}_t6dncj>(-r&
z3}^;OIrQOvOB{<W!0V6!Y$)6(vki5pApaB=lbT8<9m`WYMp;d|WnhQOroBqFHJ8HJ
zOZO#Zr#Jld)Sk(!#$<gv!UfOV`nzGoY}F~P)+@$CH8E{y(_EKtsLS7Mp-S69S$87?
z6J9NE4>uWF*V!LbXPFG0rtSs^oGQAMMHTO2IX@xrJ0H|Z@0Mjbf-I3xXo-&HviGo+
z?lHx^8TdGO!rtq6uUL8zR(4sQDTdx)3;Lo0p{uPg+vz582C^v<F?fr@RB9atppE5*
z{2fmE6#;pht3>1$I;pnDgo}MhUk=-2&H=Wh-t4KPgo$hKeH;%^DY!|g`xs(zp;WFH
zv|Al7X+^dUjMF_SJDOl#Zs`KJi-V${js-ce7xSqc#fJ2P_KB%DDA(DfyyRl*N!QCH
zy;x)GDb}G>vEO5#dpA!8udjW%Xt295h(jjwS<+gUZjee5Z<XLwIRNGXn_aAOH)wn4
z_W1(B8$)U0+Xc`WA=Cl-w~QE5UOfhkq5EbTOCkWcF?{olI}ia07``ROd`|}g;#x&G
zO>nWRLeF3$MHYDGxW~IWKjsL0TI`L551wsi>?lk`%sq9Y|5_B=Xj_8n8LUyGF;)D9
zf^ViSSWM>C{jwxcw1e<-shLY4GaJ`^&oic;?3zh=U0>4oPgv8f1OmG(q^iB^w8yX-
zye6s=V&zFG2;Ym})vwf@Lpzq!taHwo4So)LGC<Qx(bgw^Q{~qUw!iU3n+<#%-iMDv
z!t2U)wK8#mnyb0FSUXs{{;b-hL3pPH2JMLiuL`?wBGM)yKNyjftx*RHIq9d{oArvn
z-OgUyb!G|jm%$T#bdu(^LwUDBm~>EVB+-R`RjXRa(kj2?*~h{QsGj~OCmBRL>^U8I
z@7S)rs>mRw*G#oZhm1imO+W%RS$3-J=)0jIZzNjwOGD*tJR2Ke;HW8AmJudWp9`ab
zrBojY=Fq*8Pr}VV*v-+4K<GM#8sBW`B$LB(_B9v+GRuIkj;uyd>vr`nsJrVnu>13+
z$#4dzGQ6dp#MHTDJn!j)a+S~J+LewwEGW&KA!RJvs6PrRiVoq7vZMed$<(2`E|G`u
zhkvNi1zWO49aV9ot|H!AA$C)W2%<}&?47+ozGxn4^1`kf6uItF3o3EOOiOw1g{B0Q
z;$K@K(O~*PZI^EYv(iw`w#L_gXSD@39q%gS+(lPZ`d*ziT^k#iY2hK9PZX3A0xowW
zD{OOP7($;ZasV07eMhQNy2;Id?c{D&P<!aFhJT&Hgdx7WQSbrUMRUE;z$e1B+0WOX
zeP>rk$#I^a@CQr4mvVEx!>nuvlYWnJe(QjDyCd}U<HTkVURjoK6++UvadsqBw@!8>
zk9iYPO=A?zohn>{?CvSQFj4S#`Z@KqI_psrXN0u8$9o6wPpGJkx=N;%q--e6*a4Bc
z7T+T5$pjB)<D{39y&Xq5_~fNz39<2@cn+yNpn>h_@?5stpMm5wl;La>+(^uYddYWu
zNoAHW(v|?3(WJQXLG09eQZb4Wwz6dr_S2=n_#R$5!ho@s(Xi7QQU;d6{6nhX3fTzr
z_QH#HCCdT1yqU$k^oy|6eI|hzJo{A6;kXQ}{UToMh@c=ewsA%A0{7kGTEg5a3I4m&
zO5&~!4ed^<;MLPpC*}>oay46yV=eB>v^Rz=THBM!XQ~K?e%jVLFRKj?(uMfZ`d*NI
z{P?I|S94$JauVBjUg0~TuY0g3|NdobeU%AL65mJc=6P$6#^T4G!4>wOw>_G+Cy0D8
zFqpy_$*tzgG+MfCp2y+6oY(alcRxGsqiB`XY1rsnq1Pd_rM+j`jO#zqAe+z=SItaC
zOYKCJyAnrW78->FXl9N`@i%EQJ=taPz57IQ4k8*+$JqzAlxm$Jg`ARhtedmwt+9>R
ztkG~7H4FM&T<jd|St@0wi+*<(J1O|UAMwt+VOrepXc*Us?8y(<@^<0E)&X7YljDA^
zitndK3S0S-+f9}gnc_a*E~lSdUZ^8Ey1J!=;>v~vJELWhdg=f$C*=KczcqoeZb=iO
z6C$W0?et0@eoCb+?sR``({WD09gO2eFQ)p23Rxss8mXBw2*c?)*PMxdQ~G?b{<qbx
zQEhRd*rPX10Ltw(*c8Qw$U@+ubm180Yf0CvfcjC=huUPlJp6C9>0$oB?((Ltb@|uv
zyi~xKfnie&NG?prlJgA6Z3Aw|Zzm({J@&{cW!$%#PPHF&)Jdop;;VH>2`Gb9UO}TI
zUp)Lmb9+OUYT!);slM#~S5h^LuOZrii~*`FW;aKvkkKZ;=R=Z2y?0VK;)KI0MXLH4
z4gEWcl6k@dL|_Fe5@B>^SztEwf-iA-xQ``Na+>xCKnwPI3xQa9xN$4s0e8JyJnrH{
z8vv7xLw?+6J9BH&<dhj+r9m*9cZ3<)wC~=vq}>xq;}+f2Lo|P8Us^rLmiNwPEcUfO
zH*=9KZ_FJ@^jRbk3lv`uOWSbEkX+ya-JmC<P0!cQN#g}$qFYKy>JYCg>thC(bl6b7
zySM~I`;K2u?DXx}Pwuf#2S$B*K6@rV?RxAmj|J>L^S@2Q#Obz4EN$8%Ta}?&sm&^&
zDzMv!J*C@1kg!qKJ*+#s{P;!H`b#1a5~D4`{G{Vsg@6-kkZmJV42&zC;km{jRO9GL
zsZ9jz5X$DN)2~0pomHBw^Nd6DJDXDR9gy!T4_y*f&+TquIyPVAol6sH)Nw_Y8sAP{
z>};Vd#kbKyka<+gLrl8BKm-Mo7cTcYwFfj}gSnu$Aa|>kJ#`;-FS3(cyj(q8%4Y=7
zG?0ak*>d3ozZj6QJ)WFz%0H_SAojlHQ2h0R!6mVXp_mN+c<Y;N(r%>b=)`NQq=l}I
zK5rD9yGVB&I}PrdAupoRippOk^>Q@(EM~HE7c}lfSa>U}cZih5^#rFu=tx2kDX6wP
zNp`j^)b_R_l4V*BxP0yH*HFodM3<3a`jaSpVlV0j@Gi!Y2c0m^q!WsOln>d)pSw!6
zDOset&^E~;-v;?5m5DRpy9OozO-1)m$;#e5V|712x-(rPF5a)Rt@!RG9t`@nrx<gX
z`_5b;H9+wm2z3nY?S7YF9#>mryRt(|I;oRvkiM>~tQY$mElqtnt)sG6`gNssOu6F-
zgGx!h+{2%(eZY6?W0E${ia{CZxbH6`1>G0eYm#{}apE4C5TAR6`o?Hz4EIg0HCl8&
zK4LytM3aosl~O7zF4T1M&^NvxmsYN6To1_+5Z~w3j@4z=(kf_)*0cBULXKYuw&5#K
zHRpJ*m~#OmK8e?1lzm7tyQ)scb(ifZXSa^BV(HPHrbQS;ak5elyWW7U)}UVZO9QlN
zzDVTO%N`m~)9eSoQ_HNOwdb(S(|6r|Z<uN;PC+%=x@dYuPLG_8$_b4=kRyYRUTba&
zlng4CO_bK1-zl88PI_TjrnWWJ{HT86_=8OKj~?C_wb$e2j}*(MnI`)(wOnjGB-=(o
zlNpfw-HG|=%<;#(^?aqi2Th!gecN+wH#6>fc%caU8+yPuXC(Ll!~W0y;pS}z{keGL
zCmF!kWx|lX6IRfI8WJ`r+~3AJLcJ)Bg)h-#+aFwdhID&8Ug3OC7!@DieXao3_r4-F
z>*$Q9VKYyHV}Mt`gb+<jOh%rlamD8OQbpQ_G_wiqVhBr(DYMvX)CF_ugm+ZPK65-4
zctG{UrIy;`gt&qk9&j3~0S^lOGyVkKjTe;GHl0?pZ@{W?BPV8b%vxSFAn1WO5IpFr
z=M$XFSL@r4QXiTb91lu95zfc{h~YGc!lRqb9;V<`IwDsQA!|tbaih-=F=?ZC3B_mE
z{XTIJF#-B(M|u>zgHDr09aU>6x%!5t&WvBSaA?0}ep|^0Y$4~Kz6423`?z4$fnH=c
zmIlFW5*Ag8dpQCfsP|oFA`CCM0Qbcoj+%9_5&7D_A?Q`kaU}cfLaSIVDFV>BMedNQ
zBIGtg>WIUaqT<AdDxY-{yOWU<N)t82_!WENqm5dJMmlVn2^wzvEEqjx4_02=$mlSl
zJF@{|7W!!YAXtjZG78Zcld&tYsbZlR(?g4PvQK~RoV(lU;Nw;-Eh$+iB99xzk97Xm
z3(*e@+<V1y&Lv5{CPsV!7o0}=kef<o+xG{rZsCSJ)$vH2emv3m-bDw6^iGKt>&_Z)
z$`bSU0m!Gm{iNd_qZ-QUg44Gd>cA$l(rVo3UA+mi%%v>PlGw;&1FA=^PxwwOra_0q
zqc*z_s^=#kWpG@Kv|YNxY!(hI-iM;bLFsAwqarB>#T9fThVKPlKK{`TYdCp8dY+Zr
zesX9lzHh;m;d69)@9>ZbFnwU>PO>_hda@DS7Nitl%uUV4_FUenn|;@DNJEfe<2w{@
z#S7p3thIXKra=FF2{U&zQ+9E5`tu2}ceN851mB|>;m_nicpd*TGjVdd+MzQN?mNzN
zVg>J=6NP+<PIq7mk<jGIBY&h_i2$BoZtryisV5`H^CeuK_2Mku+g5B6O~4(_VlJMS
zT92S2nElLJsh^Z9;V5Z;*@_?A<fxVP;67F>Dv5~ZQW%57F<LE{*3>joeFa(RnN{jC
z1j5bcArUD8^BxPDKwRhEg<6r>?z5*%<B0ZeS2f3b8Gq%Q(sITiEbP)4{z0yb#<M><
z{)N;vQ`5-QZ3<yMYXWyFPa<b=Y5Q%qnv0makX7UwpN`m;YOi-u9V$h+>uZO0JeF@$
z$|r!%A%0^5>&c_yhn=)U=1<nk3?i3EfH;EBHfVMBL5g*1Ts9=(iH~tMgLS1D?-Z^v
zt$%8$%lhiVYwm2FL?Ao$JWRMBXY*n8)2<&hhC|^GQMt@&kj0+@V)lT;)aR^^DH{P|
zd&-Jih$1`A`-)pl#~fctJr{9Ehvi!BitRDcw=dYo>E*HAv#b-zFSdy<sJ5%uFO3Dm
z3&&@-R|j)1w3C<br29Zx4c>f5x=HMH=RwQut2BgXY#H8t@cx=HGe;Nb&y=b5U&(;Y
zC#-Y}fo2k~A|&JHBGlqpo)v32mbduJnd-sD;@T>XX#D$%>N7HC55&^hVxV}J{zrWq
zLMAVRq_czLA&%=S_AOC+;FWmWUd)#!iZ=V^@e9z7f(}71=XU{8C}J2gbA=s4y*Q%<
zNslT66V6qXHzXXQScOnG>Ly$&*Y&cO8B#~_<K$w}`9mrxDZq}*-Nx3q<yP2v*~${=
zWTH@dnPFXIt+O-~HrE^w(G;-*L}?{k>=7TJq?k-2QIqV<@~heL#as*->&H?ELN5=q
z>Z1Qk=Wl7BY}GcS3|Gjgh=}9rH^s2x#eMw(+u23^fpYsD%#n+Vl-|z(>oZaYx9|EH
zS17wO$Gqk%|MH3lBl|Oz*5U_pl*1gcle6c2XCFOy{G&P)#~cFhuK69F))*%ErbiWW
znQ_TYMdgFQh1J6a>m_H^CbQy)j_(uii!V&J4Lm{odonE8q#o<Tli>zW2JZD3_^S>T
zouHru*ZM;CE+Uvt%Gd%CC+EO`IFI`j{;Y?po{`BM-J9A!L>(T+14-DOce-MYgmuc4
zro!lD#y)9OrL3!FtF?|gX!X`l7tWRuCbjb!H+Mo{vMvTjrDA>}b}iz7Ay4&QKS<61
zx2jK#ge&#P0eJgmaF@deC?TDsluwtrf=B9l|3HApp*+*_6Vg`u+77gedy(uJw0PJI
zno2&DR=GOI)ak)q8y)iA<fN68GKG(WASN0OvqID+wr{r8oJHfk&(Z8_BQkx{)}qdc
zzZEXf4M!|}aA~oreS2~=ffBTj<?+fik5gi^auhuc;-auDexOC;vDi8ulF{^_HBRU9
z2%kfRcSh23jJ4#&#zjGYNi$R4=LGa+xrv9$qM@N5R4D!QxXRNk1PL6}^U%fo0cQsU
z3^5dZ)>P=2%zO(`X=^qJPy+B%l~Qv~`ptqK4##)oi9=_Y%5sRe2oe6fKn~BwU(UbC
z<-Z;LvtwTAy4UBbDMr9Ee4|DFJ?Ni0?+pq8q5K!-rx^MF1-p_%ug_Jmc6cfM6i06&
z{quHL^62%s0xVtw{{4b~%B43z_}>cfu<$t7u)h}V|Ge*&s(5{_!ZX7a7`UImX^b~l
zzwTIns*E>i1O$zLV1H|lzZtvs@~7hS+Y2t*&l~(h|M~6Ve^FQh`03$BrTF{x|I{gd
z2Oj$??6+$1JMO>lioYGi!g~gsJJ(;E<8SEyx-uvJUxxK>*mcA|r~2=RT7QN8XU6}D
y{&T|oM#CShe{Tpk^X5<ZpFQz6JelZMsoW^_%5o?_e``Sa_VEs0k<=tt=KlwIgkFFE
literal 0
HcmV?d00001
diff --git a/assays/AmpliconSequencingDataAnalysis/protocols/.gitkeep b/assays/AmpliconSequencingDataAnalysis/protocols/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/AmpliconSequencingDataAnalysis/protocols/AmpliconSequencingDataAnalysisProtocol.md b/assays/AmpliconSequencingDataAnalysis/protocols/AmpliconSequencingDataAnalysisProtocol.md
new file mode 100644
index 0000000..057ea25
--- /dev/null
+++ b/assays/AmpliconSequencingDataAnalysis/protocols/AmpliconSequencingDataAnalysisProtocol.md
@@ -0,0 +1,3 @@
+## Amplicon sequencing data analysis
+
+Sequencing reads were demultiplexed using Qiime2 (qiime demux emp_paired) (Bolyen, E. et al., 2019), and merged using Flash2 (MagoÄ, T. and Salzberg, S.L., 2011). Reads were denoised and dereplicated using Dada2 (Callahan, B.J. et al., 2016) and remaining individual reads were denoted as ASVs. Chimeras were removed using Qiime2 (vsearch uchime-denovo). Taxonomic classification was done via the Qiime feature classifier using the silva_138 and sequences classified as mitochondrial or chloroplast were removed from the dataset. Remaining ASVs were included in count tables. Bray-Curtis dissimilarities between samples were calculated normalized ASV tables and performed PCAs using the cmdscale function (*vegan* R package). To quantify the effects of genotype and root region we used a constrained PCoA using the *capscale* function (*vegan* R package). To quantify the contribution of different variables and their interactions to the variance in pairwise Bray-Curtis dissimilarities, we analysed the Bray-Curtis distance matrix between pairs of samples with 999 iterations of a permutation-based test (permutational multivariate analysis of variance (PERMANOVA), adonis function, vegan R package), and removed technical and batch effects, using the formula as follows: *Bray-curtis ∼ VariableX + Condition(Technical_replicates + Biological_replicates)*. We further inspected the effects of the variables using a constrained PCoA using the capscale function (vegan R package). Sequencing data from the SynCom experiment was pre-processed similarly as natural community 16S rRNA data. Quality-filtered merged paired-end reads were aligned to a reference set of sequences extracted from the whole-genome assemblies of the 60 strains included in the SynCom using Rbec79 (v1.0.0). A count table that was employed for downstream analyses of diversity was generated in R (v4.0.3) with the R package vegan (v2.5–6). Scripts for microbiota analysis and data visualization can be found at https://github.com/duranpa/sweet_collaboration
\ No newline at end of file
--
GitLab