The bacterial and fungal community in the pilot-reactor were tested using Illumina Miseq (Novogene Co., Ltd., Beijing, China) sequencing of V4-V5 16S hypervariable regions and ITS1 region respectively. To analyze fungal and bacterial taxa distribution, sampling was performed at the start-up and at the end of treatment; in the last case, samples were collected for each PUFs size. Obtained sequences were analyzed using a combination of Qiime 1.9.1 [22] and R packages. Raw data were firstly assembled in to paired-end reads using fastq-join ( (accessed on 4 March 2019). The paired reads were quality filtered and reads with a phred score Q < 20 were discarded. The files containing paired reads were converted into Fasta files using Fastx toolkit ( (accessed on 4 March 2019) and the obtained fasta files were imported into Qiime. Chimera check was performed using vsearch [23] and the reference chimera checked database ‘Gold’ in the Broad Microbiome Utilities ( (accessed on 4 March 2019) operational taxonomic unit (OTU) table was obtained using the Qiime script, method sortmerna and sumaclust. Silva non-reduntant SSU database (128 release) was used with a confidence threshold of 0.97. OTU table, taxonomy, metadata and tree files were imported in R for alpha and beta diversity analysis using Phyloseq package [24]. Data preprocessing was carried out filtering sequences less than 0.1% in abundance and present in less than 2/3 of replicates. For statistical analysis, the OTU table was log transformed as previously described [25] via the implemented script in Qiime [26]. Alpha diversity indices were calculated using diversity function in Vegan 2.5.4 package and the estimated richness function in phyloseq. One-way ANOVA and a Tukey test were used to estimate the statistical differences of alpha diversity indices among metadata categories. Weighted Unifrac distant matrix was used to perform principal coordinates analysis (PCoA) applying permutational analyses of variance (PERMANOVA statistical test). A similarity percentage analysis (SIMPER) was performed using simper function in Vegan package, to explore the OTU that contribute to differences among metadata categories. Summarized taxa table at phylum and genus level (where possible) were reported to highlight the phylogenetic groups that contribute to dissimilarities. To estimate the fold-change at genus level among metadata categories, DESeq2 package was used [27]. The shared OTUs in the different reactor phases were estimated by using Venn diagram, a freely online tool available at (accessed on 2 April 2019) [28]. Detailed steps for the analysis of core microbiota are given along with raw script (accessed on 6 April 2019).

Note: The content above has been extracted from a research article, so it may not display correctly.

Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.

We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.