The V3 and V4 regions of 16S rRNA gene were amplified using primers (341F 5′-CCTACGGGAGGCAGCAG-3′ and 806R 5′-GGACTACNNGGGTATCTAAT-3′). PCR products were purified with GeneJET Gel Extraction Kit (Thermo Fisher Scientific, Massachusetts, USA). Sequencing libraries were generated using NEB Next Ultra DNA Library Prep Kit for Illumina (NEB, Ipswich, USA) and library quality was assessed on the Qubit@ 2.0 Flourometer (Thermo Fisher Scientific, Massachusetts, USA) and Agilent Bioanalyzer 2100 system. The amplicon library was sequenced on an Illumina HiSeq 2500 platform, generating 250 bp paired-end reads. Paired-end reads were merged using FLASH (V1.2.7, http://ccb.jhu.edu/software/FLASH/)39 (link) and assigned to each sample according to the unique barcodes. The sequence reads were quality filtered using QIIME (V1.7.0)40 (link) and compared with the reference database (Gold database, http://drive5.com/uchime/uchime_download.html) using UCHIME algorithm41 (link) to detect and remove chimeric sequences42 (link). The reads were clustered using Uparse software (Uparse v7.0.1001)43 (link) and OTUs (Operational Taxonomic Units) generated based on 97% sequence homology and OTU representative sequences were then classified taxonomically using the QIIME-based wrapper of the Ribosomal Database Project (Version 2.2)44 (link).
Free full text: Click here