Abstract
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 co-occurrence 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.