Succession of diversity, functions, and interactions of the fungal community in activated sludge under aromatic hydrocarbon stress

Although fungi are regarded as the important degraders of aromatic hydrocarbons (AHs) in various environments, the dynamic succession and interaction of their community under aromatic hydrocarbon stress has been rarely reported. In this study, we systematically investigated the responses of the fungal community and the associations among fungal species when facing the continuous stress of two typical AHs, benzene and naphthalene. Using high-throughput sequencing technology, we demonstrated that fungal diversity displayed a significant downward trend during six weeks of continuous aromatic hydrocarbon treatment. Community succession was observed during the operational period, and the relative abundance of some typical degraders, such as Exophiala sp. and Candida sp., increased during the later period of operation. Meanwhile, by predicting the functions of the fungal community through PICRUSt2, we found that some relevant enzymes, such as peroxidase, dioxygenase, and monooxygenase, may play an important role in the degradation process and maintaining overall community multifunctionality. Furthermore, the measurement of modified normalized stochasticity ratio (MST) indicated that continuous aromatic hydrocarbon stress resulted in a stronger deterministic process in community assembly over time, suggesting environmental selection dominated succession of the fungal community in activated sludge. Finally, molecular ecological network analysis (MENA) demonstrated that, the cooperative behaviors among members, the network keystone genera related to biodegradation, such as Exophiala sp. and Haglerozyma sp., and a well-organized topological structure, together, maintained the structural stability of the fungal community under AH stress. Our study provides new insights for understanding the stability of fungal communities during the degradation of contaminants in activated sludge.