Different soybean genotypes can differ in their tolerance toward aluminum stress depending on their rhizosphere-inhabiting microorganisms. However, there is limited understanding of the response of fungal communities to different aluminum concentrations across different genotypes. Here, we used metabarcoding of fungal ribosomal markers to assess the effects of aluminum stress on the rhizosphere fungal community of aluminum-tolerant and aluminum-sensitive soybean genotypes. Shifts in fungal community structure were related to changes in plant biomass, fungal abundance and soil chemical properties. Aluminum stress increased the difference in fungal community structure between tolerant and sensitive genotypes. Penicillium, Cladosporium and Talaromyces increased with increasing aluminum concentration. These taxa associated with the aluminum-tolerant genotypes were enriched at the highest aluminum concentration. Moreover, complexity of the cooccurrence network associated with the tolerant genotypes increased at the highest aluminum concentration. Collectively, increasing aluminum concentrations magnified the differences in fungal community structure between the two studied tolerant and sensitive soybean genotypes. This study highlights the possibility to focus on rhizosphere fungal communities as potential breeding target to produce crops that are more tolerant toward heavy metal stress or toxicity in general.

Aluminum (Al)-resistant plant cultivars can recruit beneficial microbes to alleviate the stresses. However, the mechanism of how rhizobacterial communities strengthen Al tolerance of wild soybean has not been addressed. The aim of this study was to investigate the bacterial community structure in the rhizosphere of Al-tolerant (BW69) and Al-sensitive (W270) wild soybean germplasm subjected to three Al concentrations. We analyzed the rhizobacterial communities of the two genotypes by high-throughput sequencing of 16S rRNA genes. The results showed that high Al stress recruited different rhizobacterial communities between two genotypes. In total, 49 OTUs, such as OTU15 (Gammaproteobacteria_KF-JG30-C25_norank), OTU23 (Mizugakiibacter), and OTU93 (Alkanibacter), were enriched in the rhizosphere of BW69 at the low and high Al concentrations. Moreover, bacterial community in the rhizosphere of BW69 had a more complex co-occurrence network than did W270 at the high Al concentration. Overall, our findings highlighted that high Al concentration magnified the difference in rhizobacterial community structure between two genotypes. However, the lower modularity of the co-occurrence network in rhizosphere of BW69 than W270 under Al stress may cause the rhizobacterial community to be less resistant and more influenced by disturbance. This study emphasizes the possibility of using rhizobacteria as an improved crop breeding or gene to produce crops that are more resistant to the toxicity of heavy metal.