Little is known about the changes in soil microbial phosphorus (P) cycling potential during terrestrial ecosystem management and restoration, although much research aims to enhance soil P cycling. Here, we used metagenomic sequencing to analyse 18 soil microbial communities at a P-deficient degraded mine site in southern China where ecological restoration was implemented using two soil ameliorants and eight plant species. Our results show that the relative abundances of key genes governing soil microbial P-cycling potential were higher at the restored site than at the unrestored site, indicating enhancement of soil P cycling following restoration. The gcd gene, encoding an enzyme that mediates inorganic P solubilization, was predominant across soil samples and was a major determinant of bioavailable soil P. We reconstructed 39 near-complete bacterial genomes harboring gcd, which represented diverse novel phosphate-solubilizing microbial taxa. Strong correlations were found between the relative abundance of these genomes and bioavailable soil P, suggesting their contributions to the enhancement of soil P cycling. Moreover, 84 mobile genetic elements were detected in the scaffolds containing gcd in the 39 genomes, providing evidence for the role of phage-related horizontal gene transfer in assisting soil microbes to acquire new metabolic potential related to P cycling.

Phosphate‐solubilizing microbes (PSMs) drive the biogeochemical cycling of phosphorus (P) and hold promise for sustainable agriculture. However, their global distribution, overall diversity and application potential remain unknown. Here, we present the first synthesis of their biogeography, diversity and utility, employing data from 399 papers published between 1981 and 2017, the results of a nationwide field survey in China consisting of 367 soil samples, and a genetic analysis of 12986 genome‐sequenced prokaryotic strains. We show that at continental to global scales, the population density of PSMs in environmental samples is correlated with total P rather than pH. Remarkably, positive relationships exist between the population density of soil PSMs and available P, nitrate‐nitrogen and dissolved organic carbon in soil, reflecting functional couplings between PSMs and microbes driving biogeochemical cycles of nitrogen and carbon. More than 2704 strains affiliated with at least nine archaeal, 88 fungal and 336 bacterial species were reported as PSMs. Only 2.59% of these strains have been tested for their efficiencies in improving crop growth or yield under field conditions, providing evidence that PSMs are more likely to exert positive effects on wheat growing in alkaline P‐deficient soils. Our systematic genetic analysis reveals five promising PSM genera deserving much more attention.

Phosphorus (P) limitation of ecosystem processes is widespread in terrestrial habitats. While a few auxiliary metabolic genes (AMGs) in bacteriophages from aquatic habitats are reported to have the potential to enhance P-acquisition ability of their hosts, little is known about the diversity and potential ecological function of P-acquisition genes encoded by terrestrial bacteriophages. Here, we analyze 333 soil metagenomes from five terrestrial habitat types across China and identify 75 viral operational taxonomic units (vOTUs) that encode 105 P-acquisition AMGs. These AMGs span 17 distinct functional genes involved in four primary processes of microbial P-acquisition. Among them, over 60% (11/17) have not been reported previously. We experimentally verify in-vitro enzymatic activities of two pyrophosphatases and one alkaline phosphatase encoded by P-acquisition vOTUs. Thirty-six percent of the 75 P-acquisition vOTUs are detectable in a published global topsoil metagenome dataset. Further analyses reveal that, under certain circumstances, the identified P-acquisition AMGs have a greater influence on soil P availability and are more dominant in soil metatranscriptomes than their corresponding bacterial genes. Overall, our results reinforce the necessity of incorporating viral contributions into biogeochemical P cycling.

Nucleocytoplasmic large DNA viruses (NCLDVs; also called giant viruses), constituting the phylum Nucleocytoviricota, can infect a wide range of eukaryotes and exchange genetic material with not only their hosts but also prokaryotes and phages. A few NCLDVs were reported to encode genes conferring resistance to beta‑lactam, trimethoprim, or pyrimethamine, suggesting that they are potential vehicles for the transmission of antibiotic resistance genes (ARGs) in the biome. However, the incidence of ARGs across the phylum Nucleocytoviricota, their evolutionary characteristics, their dissemination potential, and their association with virulence factors remain unexplored. Here, we systematically investigated ARGs of 1416 NCLDV genomes including those of almost all currently available cultured isolates and high-quality metagenome-assembled genomes from diverse habitats across the globe. We reveal that 39.5% of them carry ARGs, which is approximately 37 times higher than that for phage genomes. A total of 12 ARG types are encoded by NCLDVs. Phylogenies of the three most abundant NCLDV-encoded ARGs hint that NCLDVs acquire ARGs from not only eukaryotes but also prokaryotes and phages. Two NCLDV-encoded trimethoprim resistance genes are demonstrated to confer trimethoprim resistance in Escherichia coli. The presence of ARGs in NCLDV genomes is significantly correlated with mobile genetic elements and virulence factors.