Archaeal communities perform an important role in maintaining microbial stability under long term continuous cropping systems.

Long-term continuous cropping of soybean can generate the development of disease-suppressive soils. However, whether the changes in microbial communities, especially for archaea, contribute to controlling soil sickness and improving crop yields remains poorly understood. Here, real-time PCR and high-throughput sequencing were employed to investigate the changes in soil archaeal communities in both bulk and rhizosphere soils under four cropping systems, including the continuous cropping of soybeans for a short-term of 3 and 5 years (CC3 and CC5, respectively) and for a long-term of 13 years (CC13), as well as a soybean-maize rotation for 5 years (CR5). The results showed that CC13 and CR5 significantly increased archaeal abundance, reduced the alpha-diversity of archaeal communities, and changed soil archaeal community structures compared to CC3 and CC5. Microbial co-occurrence network analysis revealed that CC13 led to the higher resistant microbial community and lower the relative abundance of potential plant pathogens in the network compared to CC3 and CC5. Correlation analysis showed that the microbial resistance index was negatively correlated with the relative abundance of potential plant pathogens and positively correlated with soybean yields in both bulk and rhizosphere soils. Intriguingly, the random forest (RF) analysis showed that archaea contributed the most to soil microbial resistance even though they were not at the core positions of the network. Overall, structural equation models (SEMs) revealed that high resistant microbial community could directly or indirectly improved soybean yields by regulating the relative abundance of plant pathogens and the soil nutrients, suggesting that the regulation of soil microbial taxa may play an important role in maintaining agricultural productivity under continuous cropping of soybean.