High throughput sequencing-based analysis of the soil bacterial community structure and functions of Tamarix shrubs in the lower reaches of the Tarim River.

Abstract

Tamarix is a dominant species in the Tarim River Basin, the longest inland river in China. Tamarix plays an important role in the ecological restoration of this region. In this study, to investigate the soil bacterial community diversity in Tamarix shrubs, we collected soil samples from the inside and edge of the canopy and the edge of nebkhas and non-nebkhas Tamarix shrubs located near the Yingsu section in the lower reaches of Tarim River. High throughput sequencing technology was employed to discern the composition and function of soil bacterial communities in nebkhas and non-nebkhas Tamarix shrubs. Besides, the physicochemical properties of soil and the spatial distribution characteristics of soil bacteria and their correlation were analyzed. The outcomes of this analysis demonstrated that different parts of Tamarix shrubs had significantly different effects on soil pH, total K (TK), available K (AK), ammonium N (NH4+), and available P (AP) values (P < 0.05), but not on soil moisture (SWC), total salt (TDS), electrical conductivity (EC), organic matter (OM), total N (TN), total P (TP), and nitrate N (NO3−) values. The soil bacterial communities identified in Tamarix shrubs were categorized into two kingdoms, 71 phyla, 161 classes, 345 orders, 473 families, and 702 genera. Halobacterota, unidentified bacteria, and Proteobacteria were found to be dominant phyla. The correlation between the soil physicochemical factors and soil bacterial community was analyzed, and as per the outcomes OM, AK, AP, EC, and NH4+ were found to primarily affect the structure of the soil bacterial community. SWC, TK and pH were positively correlated with each other, but negatively correlated with other soil factors. At the phyla level, a significantly positive correlation was observed between the Halobacterota and AP, OM as well as Bacteroidota and AK (P < 0.01), but a significantly negative correlation was observed between the Chloroflexi and AK, EC (P < 0.01). The PICRUSt software was employed to predict the functional genes. A total of 6,195 KEGG ortholog genes were obtained. The function of soil bacteria was annotated, and six metabolic pathways in level 1, 41 metabolic pathways in level 2, and 307 metabolic pathways in level 3 were enriched, among which the functional gene related to metabolism, genetic information processing, and environmental information processing was found to have the dominant advantage. The results showed that the nebkhas and canopy of Tamarix shrubs had a certain enrichment effect on soil nutrients content, and bacterial abundance and significant effects on the structure and function of the soil bacterial community.