Comparative evaluation of 16S rRNA primer pairs in identifying nitrifying guilds in soils under long-term organic fertilization and water management

Compared with 454 sequencing technology, short-read sequencing (e.g., Illumina) technology generates sequences of high accuracy, but limited length (<500 bp). Such a limitation can prove that studying a target gene using a large amplicon (>500 bp) is challenging. The ammonia monooxygenase subunit A (amoA) gene of ammonia-oxidizing archaea (AOA), which plays a crucial part in the nitrification process, is such a gene. By providing a full overview of the community of a functional microbial guild, 16S ribosomal ribonucleic acid (rRNA) gene sequencing could overcome this problem. However, it remains unclear how 16S rRNA primer selection influences the quantification of relative abundance and the identification of community composition of nitrifiers, especially AOA. In the present study, a comparison was made between the performance of primer pairs 338F-806R, 515F-806R, and 515F-907R to a shotgun metagenome approach. The structure of nitrifier communities subjected to different long-term organic matter amendment and water management protocols was assessed. Overall, we observed higher Chao1 richness diversity of soil total bacteria by using 515F-806R compared to 338F-806R and 515F-907R, while higher Pielou’s evenness diversity was observed by using 515F-806R and 515F-907R compared to 338F-806R. The studied primer pairs revealed different performances on the relative abundance of Thaumarchaeota, AOB, and NOB. The Thaumarchaeota 16S rRNA sequence was rarely detected using 338F-806R, while the relative abundances of Thaumarchaeota detected using 515F-806R were higher than those detected by using 515F-907R. AOB showed higher proportions in the 338F-806R and 515F-907R data, than in 515F-806R data. Different primers pairs showed significant change in relative proportion of NOB. Nonetheless, we found consistent patterns of the phylotype distribution of nitrifiers in different treatments. Nitrosopumilales (NP) and Nitrososphaerales (NS) clades were the dominant members of the AOA community in soils subject to controlled irrigation, whereas Ca. Nitrosotaleales (NT) and NS clades dominated the AOA community in soils subject to flooding irrigation. Nitrospira lineage II was the dominant NOB phylotype in all samples. Overall, ideal 16S rRNA primer pairs were identified for the analysis of nitrifier communities. Moreover, NP and NT clades of AOA might have distinct environmental adaptation strategies under different irrigation treatments.