The novel multi-million read generating sequencing technologies are very promising for resolving the immense soil 16S rRNA gene bacterial diversity. Yet they have a limited maximum sequence length screening ability, restricting studies in screening DNA stretches of single 16S rRNA gene hypervariable (V) regions. The aim of the present study was to assess the effects of properties of four consecutive V regions (V3-6) on commonly applied analytical methodologies in bacterial ecology studies. Using an in silico approach, the performance of each V region was compared with the complete 16S rRNA gene stretch. We assessed related properties of the soil derived bacterial sequence collection of the Ribosomal Database Project (RDP) database and concomitantly performed simulations based on published datasets. Results indicate that overall the most prominent V region for soil bacterial diversity studies was V3, even though it was outperformed in some of the tests. Despite its high performance during most tests, V4 was less conserved along flanking sites, thus reducing its ability for bacterial diversity coverage. V5 performed well in the non-redundant RDP database based analysis. However V5 did not resemble the full-length 16S rRNA gene sequence results as well as V3 and V4 did when the natural sequence frequency and occurrence approximation was considered in the virtual experiment. Although, the highly conserved flanking sequence regions of V6 provide the ability to amplify partial 16S rRNA gene sequences from very diverse owners, it was demonstrated that V6 was the least informative compared to the rest examined V regions. Our results indicate that environment specific database exploration and theoretical assessment of the experimental approach are strongly suggested in 16S rRNA gene based bacterial diversity studies.

The present research is the first comprehensive study regarding the thermophilic anaerobic degradation of cheese wastewater, which combines the evaluation of different reactor configurations (i.e. single and two-stage continuous stirred tank reactors) on the process efficiency and the in-depth characterization of the microbial community structure using genome-centric metagenomics. Both reactor configurations showed acidification problems under the tested organic loading rates (OLRs) of 3.6 and 2.4 g COD/L-reactor day and the hydraulic retention time (HRT) of 15 days. However, the two-stage design reached a methane yield equal to 95% of the theoretical value, in contrast with the single stage configuration, which reached a maximum of 33% of the theoretical methane yield. The metagenomic analysis identified 22 new population genomes and revealed that the microbial compositions between the two configurations were remarkably different, demonstrating a higher methanogenic biodiversity in the two-stage configuration. In fact, the acidogenic reactor of the serial configuration was almost solely composed by the lactose degrader Bifidobacterium crudilactis UC0001. The predictive functional analyses of the main population genomes highlighted specific metabolic pathways responsible for the AD process and the mechanisms of main intermediates production. Particularly, the acetate accumulation experienced by the single stage configuration was mainly correlated to the low abundant syntrophic acetate oxidizer Tepidanaerobacter acetatoxydans UC0018 and to the absence of aceticlastic methanogens.