diff --git a/assays/AmpliconSequencingDataAnalysis/README.md b/assays/AmpliconSequencingDataAnalysis/README.md new file mode 100644 index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 diff --git a/assays/AmpliconSequencingDataAnalysis/dataset/.gitkeep b/assays/AmpliconSequencingDataAnalysis/dataset/.gitkeep new file mode 100644 index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 diff --git a/assays/AmpliconSequencingDataAnalysis/isa.assay.xlsx b/assays/AmpliconSequencingDataAnalysis/isa.assay.xlsx new file mode 100644 index 0000000000000000000000000000000000000000..419398f43d7d6f9ede46238daf8a729f41599c39 Binary files /dev/null and b/assays/AmpliconSequencingDataAnalysis/isa.assay.xlsx differ diff --git a/assays/AmpliconSequencingDataAnalysis/protocols/.gitkeep b/assays/AmpliconSequencingDataAnalysis/protocols/.gitkeep new file mode 100644 index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 diff --git a/assays/AmpliconSequencingDataAnalysis/protocols/AmpliconSequencingDataAnalysisProtocol.md b/assays/AmpliconSequencingDataAnalysis/protocols/AmpliconSequencingDataAnalysisProtocol.md new file mode 100644 index 0000000000000000000000000000000000000000..057ea25347d00bffb5aee95e6b3ce34222f10d4e --- /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