Abstract
Bacteria play a key role in phosphate solubilization, but related genome-centric research on agricultural microbiomes is scarce. Here, we reconstructed 472 metagenome-assembled genomes (MAGs) covering agricultural soils from six long-term field trials across China. A total of 79 MAGs contained gcd encoding quinoprotein glucose dehydrogenase (GCD), which is the key biomarker for phosphate solubilization. Our findings showed that all GCD-MAGs represent potentially novel species, with gcd copy numbers varying from 1 to 10 per genome. Large genome size, a high ratio of glycosyl hydrolase genes, and increased capacity for carbohydrate utilization were specific traits of GCD-MAGs. Notably, the gcd copy number showed a significant and positive correlation with genome size. Generated using a machine learning approach, our findings were validated in a dataset of 692 genotypes covering the 18 bacterial families to which the 79 GCD-MAGs belong. Our results improve the knowledge of both the diversity and the genetic composition of phosphate-solubilizing bacteria. In particular, they reveal a genomic link between phosphate solubilization capacity and increased potential for carbohydrate metabolism, which may accelerate targeted engineering and improve management practices for sustainable agriculture.