Unraveling the functional instability of bacterial consortia in crude oil degradation via integrated co-occurrence networks

Introduction Soil ecosystems are threatened by crude oil contamination, requiring effective microbial remediation. However, our understanding of the key microbial taxa within the community, their interactions impacting crude oil degradation, and the stability of microbial functionality in oil degradation remain limited. Methods To better understand these key points, we enriched a crude oil-degrading bacterial consortium generation 1 (G1) from contaminated soil and conducted three successive transfer passages (G2, G3, and G4). Integrated Co-occurrence Networks method was used to analyze microbial species correlation with crude oil components across G1-G4. Results and discussion In this study, G1 achieved a total petroleum hydrocarbon (TPH) degradation rate of 32.29% within 10 days. Through three successive transfer passages, G2-G4 consortia were established, resulting in a gradual decrease in TPH degradation to 23.14% at the same time. Specifically, saturated hydrocarbon degradation rates ranged from 18.32% to 14.17% among G1-G4, and only G1 exhibited significant aromatic hydrocarbon degradation (15.59%). Functional annotation based on PICRUSt2 and FAPROTAX showed that functional potential of hydrocarbons degradation diminished across generations. These results demonstrated the functional instability of the bacterial consortium in crude oil degradation. The relative abundance of the Dietzia genus showed the highest positive correlation with the degradation efficiency of TPH and saturated hydrocarbons (19.48, 18.38, p < 0.05, respectively), Bacillus genus demonstrated the highest positive correlation (21.94, p < 0.05) with the efficiency of aromatic hydrocarbon degradation. The key scores of Dietzia genus decreased in successive generations. A significant positive correlation (16.56, p < 0.05) was observed between the Bacillus and Mycetocola genera exclusively in the G1 generation. The decline in crude oil degradation function during transfers was closely related to changes in the relative abundance of key genera such as Dietzia and Bacillus as well as their interactions with other genera including Mycetocola genus. Our study identified key bacterial genera involved in crude oil remediation microbiome construction, providing a theoretical basis for the next step in the construction of the oil pollution remediation microbiome.


Introduction
Crude oil is a mixture of liquid hydrocarbons encompassing aliphatic, aromatics, resins, and asphaltenes (Varjani, 2017).When crude oil spills into the environment, it presents a potential hazard to both the ecosystem and biotic entities (Tang et al., 2011;Wang et al., 2021).Biological methods are more environment-friendly and costeffective compared to physical and chemical methods (Vidonish et al., 2016;Johnson and Affam, 2018;Zhou et al., 2019).Microbial consortia are preferred over single microorganisms, because they possess distinct strengths in terms of degradation efficiency and adaptability to intricate environmental conditions (Abbasian et al., 2015;Gurav et al., 2017;Cui et al., 2020).
The constructing approaches of bacterial consortia include bottom-up and top-down (Lawson et al., 2019).Bottom-up is a random combination of single bacteria, which is also a common method to build crude oil degradation consortium.For example, several researchers constructed hydrocarbon-degradation consortia by combining biosurfactant-producing strains and crude oildegrading strains (Lee et al., 2018;Chen et al., 2020;Cui et al., 2020;Dai et al., 2020), non-alkane-consuming and alkane degrader (Hu et al., 2020).The approach is full of randomness based on the complementary functions.By contrast, the top-down approach enriched specific function consortium by applying specific environmental pressure on the natural consortium (Lee et al., 2013).However, due to varying complexities of pollutants, functionality stability from top-down approach varies across successive generations.For microbial communities involved in degrading single-component contaminants, their degradation functionalities remain unaffected by the process of successive generations (Feng et al., 2019;Guo et al., 2021;Yin et al., 2021).However, for microbial communities engaged in degrading mixed pollutants, the stability of their degradation potential diminishes across generations (Chen et al., 2019;Gao et al., 2022;Lewin et al., 2022).Petroleum constitutes a complex mixture, adopting the top-down approach for selecting degradative microbial communities may lead to functional instability.To acquire microbial communities capable of crude oil degradation, it is vital to uncover the factors contributing to instability within top-down populations and identify crucial stabilizing species.
Keystone taxa were not numerically dominant in the communities (Banerjee et al., 2018), they had an effect on the structure and function of the consortium (Berry and Widder, 2014).By constructing species interaction networks, researchers have investigated microbial communities in groundwater and soil contaminated with hydrocarbons.They have identified key species exhibiting distinct traits: an elevated affinity for carbon sources, an abundance of genes related to hydrocarbon degradation (Ma et al., 2021), and the ability to produce metabolites that nourish other microorganisms (Geng et al., 2022;Jia et al., 2023).Recently, cooccurrence network analyses using high-throughput metagenomic data have been developed to identify key species and syntrophic relationships in microbial consortia.Involved measures include Pearson correlation, Spearman correlation, Bray-Curtis dissimilarity, Mutual information, and GBLM.Compared with Spearman and Pearson correlations, GBLM can mitigate spurious correlations among non-independent measurements as the increase in one relative abundance must be accompanied by a compositional decrease in another (Faust et al., 2012;Mac Aogain et al., 2021).Mac Aogain et al. ( 2021) integrated these five similarity measures to study diseaserelated microorganisms.To find key species and better understand the interrelationships within the microbial consortium, it would be beneficial to generate co-occurrence networks by merging the five similarity measures.
The objectives of this study were as follows: (1) enrich the crude oil-degrading bacterial consortia; (2) identify the key genera present and syntrophic relationships in the consortia; and (3) analyze the functional stability of the consortia through successive acclimation.

Enrichment and successive transfer culture of crude oil-degrading bacterial consortia
Soil samples polluted with crude oil were collected near the Luliang oil field in Xinjiang, China.The soil samples showed that approximately 70% of the crude oil degraded after 2-3 months in a natural environment.To enrich the crude oil-degrading bacterial consortia, 5 g of soil was cultivated for 0.5 h at 200 rpm in a 100 mL flask containing 45 mL 0.9% NaCl under aerobic conditions, at 37°C, then allowed to stand for 1 h to obtain the leachate.The leachate was transferred to the shake flasks containing 100 mL AHDM and the culture was generation 1 (G1) culture.AHDM without any bacteria was used as the control.The AHDM was autoclaved at 121°C and 30 min before the cultivation experiment.The G1 cultures was incubated on an electrical shaker at 150 rpm and 37°C for 10d.At the end of cultivation, G1 culture was continuously transferred (10% v: v) into fresh 100 mL AHDM every 10 days.The subsequent culture generations were referred to by their consecutive batch culture numbers (G2, G3, and G4).

Separation of saturates, aromatics, resins, and asphaltenes
TPH content was assessed using a gravimetric method (Capelli et al., 2001).Briefly, residual oil from the culture (100 mL) was recovered by adding 40 mL carbon tetrachloride.The lower organic phase was filtered and dehydrated with anhydrous sodium sulfate.The TPH content was quantified gravimetrically, after solvent evaporation.
The residual oil was extracted and divided into saturates, aromatics, resins, and asphaltenes, based on the Chinese National Standard SY/T 5119-2016.Briefly, n-hexane was gradually added to the residual oil.The oil sample was left to stand for more than 12 h to fully precipitate the asphaltenes and measured.The filtrate was then subjected to chromatographic separation using a silica gel-neutral alumina column.Successive elution was performed using distinct organic solvents to sequentially recover the saturated hydrocarbon, aromatic hydrocarbon, and resin fractions.Saturated hydrocarbons were eluted with n-hexane, while a mixture of dichloromethane and n-hexane was employed for eluting aromatic hydrocarbons.Resin components were separated using anhydrous ethanol and chloroform as eluents.Following complete solvent evaporation, the weights of these three fractions were individually recorded.

DNA extraction and 16S rRNA gene-based community analysis
The medium was filtered through a sterile 0.45 μm organic membrane filter (Tianjin Jinteng experimental equipment Co. Ltd., Tianjin, China).Filtered membrane was utilized for DNA extraction.DNA extraction, PCR amplification, and sequencing were performed by Guangdong Magigene Biotechnology Co., Ltd.(Guangzhou, China).The V4 and V5 hypervariable regions of microbial 16S rRNA were amplified using the primers 515F (5′-GTGCCAGCMGCCG CGGTAA-3′) and 909R (5′-CCCCGYCAATTCMTTTRAGT-3′).The constructed libraries were pooled and sequenced using 250PE (paired-end) sequencing on an Illumina novaseq platform (Illumina, USA).Analysis of bacterial consortium composition and diversity was performed in triplicate.Demultiplexing, quality filtering, clustering into Amplicon-sequence variant (ASV), and construction of the ASV table were performed in QIIME2, v.2020.2(Bolyen et al., 2019).The analysis yielded 1,503,255 high-quality reads (94.7% of reads averaged ≥ Q30 scores) distributed across 12 samples, with the minimum and maximum number of reads per sample being 88,513 and 171,446, respectively.The deblurring algorithm was used to construct the ASV table (Amir et al., 2017).The final ASV abundance table was then rarefied at 88513 sequences per sample through the alpha-rarefaction subcommand in QIIME2.Rarefaction curves based on observed ASV and the Good's coverage index were used to evaluate the adequacy of the sampling depth, which was established at 88513 quality-filtered reads per sample.The feature classifier script implemented in QIIME2 was employed for the taxonomic assessment using the SILVA reference database, v. 138 (Quast et al., 2013).Alpha diversity was assessed using diversity (Shannon) and evenness indices.Beta diversity was assessed by computing the Bray-Curtis distance between the samples and by Principal coordinate analysis (PCoA).Core OTUs within G1-G4 were discerned based on specific criteria: OTUs with high frequency, present in over 80% of samples, and abundant OTUs with relative abundances exceeding 0.2% across the entirety of samples (Wu et al., 2019).The neighbor-joining phylogenetic tree of the 16S rRNA gene was constructed using MEGA11 software, employing representative sequences of highly abundant bacteria (relative abundances exceeding 1%).The resulting tree was then visualized using the Interactive Tree of Life (iTOL) platform. 1In this study, the potential functions of microbial consortia within G1-G4 were predicted by using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2, v2.5.2) and Functional Annotation of Prokaryotic Taxa (FAPROTAX, v1.2.4).

Co-occurrence analysis of microbial interaction
Weighted co-occurrence analysis with an ensemble of similarity measures (Bray-Curtis, Pearson, Spearman, and MI) and regression techniques (GBLM) (Faust et al., 2012) were used to generate microbial association networks (Mac Aogain et al., 2021).
Specifically, the input file containing the relative abundance of genera in G1-G4 was calculated by five diverse measures and obtained relationships (positive and negative) between genera in G1-G4.The Gephi and Cytoscape software was used to analyze Betweenness Centrality, Closeness Centrality, and Degree of nodes.Following the result of Banerjee et al. (2018), keystone species have the characteristics of low Betweenness Centrality, high Closeness Centrality, and high Degree.To obtain the scores of the genera of the keystone species, the three indices were homogenized and used Eq.1:

Statistical analysis
An input file, which contained the relative abundance of all genera in G1-G4 and the degradation efficiency of TPH, saturates, and aromatics, calculated the relationship between the relative abundance of genera and the degradation efficiency of TPH, saturates, and aromatics by five diversity measures.Alpha diversity in different samples was compared using the Kruskal-Wallis pairwise test (Liddicoat et al., 2020).We used Permutational multivariate analysis of variance (PERMANOVA) to investigate the effects of time and four transfers on each component of beta diversity.Network metrics such as node degree, stress centrality, and betweenness centrality were calculated and visualized using Cytoscape.Statistical analysis and visualization of bacterial consortium and alpha and beta diversities were performed using the R packages reshape2, ggplot2, dplyr, VennDiagram, igraph, and vegan.The efficiency of crude oil degradation was assessed using Eq. 2

Degradation efficiency
Residual content of control treat Data are expressed as mean ± standard deviation (SD) of three replicates.The Wilcoxon-Mann-Whitney test and Kruskal-Wallis test were performed for statistical analysis of values in different groups to determine the presence of significant differences.

Crude oil biodegradation across four successive transfers
In this study, we enriched a crude oil-degrading bacterial consortium generation 1 (G1) from contaminated soil, achieving a TPH degradation rate of 32.29%.The crude oil degradation capacity dramatically decreased with G1 being continuously transferred.In comparison to G1, G4 consortium only achieved a crude oil degradation rate of 23.14% within the same incubation time (Figure 1A).The saturates hydrocarbons were degraded at 18.72, 22.89, 17.82, and 14.16% by bacterial consortia G1-G4 (Figure 1B), this indicates a decline in the degradation capacity of saturated hydrocarbons after G2.While the G1 consortium was capable of degrading aromatics, G2, G3, and G4 consortia exhibited negligible degradation of aromatics (Figure 1C), suggesting a loss in aromatics degradation capacity after G1.The degradation of resins and asphaltenes showed no significant difference between the bio-treatment and control groups (Supplementary Figures S1A,B).These results illustrate that TPH degradation efficiency declined during successive transfer cultures.

Dynamics of bacterial community composition
The rarefaction curve analysis indicated that the sequencing depth in our study was adequate to capture the diversity of bacterial communities in each sample, ensuring representative coverage (Supplementary Figure S2A).Furthermore, statistical analysis using the Kruskal-Wallis test revealed no significant differences in the Shannon diversity index and evenness index among the G1-G4 consortia (p > 0.1; Supplementary Figures S2B,C).When assessing beta diversity based on the Bray-Curtis distance, PERMANOVA analysis demonstrated no significant difference (p > 0.1) in bacterial community structure between G1-G4 consortia (Supplementary Figure S2D).
Fourteen shared OTUs were observed in the four generations, while 24, 10, 9, and 16 unique OTUs were found, respectively, in G1, G2, G3, and G4 (Figure 2A).There were 5 OTUs in 14 shared OTUs identified as the regional core OTUs, accounting for 20.0 to 51.6% of the bacterial abundance of four generations (Figure 2B).All the 59 unique OTUs were classified to 5 phyla, including Actinobacteriota, Bacteroidota, Chloroflexi, Firmicutes, and Proteobacteria (Figure 2C).Compared to another generations, the unique OTUs of G1 were predominantly concentrated in Firmicutes and Proteobacteria.Those of G2 were mainly centered around Bacteroidota, and G4 primarily exhibited a concentration of Actinobacteriota.
In Figure 2D, the representative sequences of highly abundant bacteria (relative abundances exceeding 1%) in this study were selected to construct the phylogenetic tree.We found that the bacterial consortia of G1-G4 had the characteristic of high diversity.In G1, dominant OTUs were primarily concentrated in the Bacillus, Dietzia, Novosphingobium, Xanthomonas genus, and the Pseudomonadales order with Acinetobacter and Pseudomonas genera.In G2, predominant OTUs were mainly centered around Pseudomonas, Dietzia, Novosphingobium, and Paenibacillus.For G3, dominant OTUs were primarily found in Novosphingobium, Sporosarcina genus, and the Enterobacterales order.In G4, predominant OTUs were focused on the Luteibacter genus and the Enterobacterales order.

Functional prediction analysis by FAPROTAX and PICRUSt2
FAPROTAX and PICRUSt2 software were be used to predict the function of bacterial communities.FAPROTAX was used to predict the biochemical cycle of environmental samples.To display potential function of crude-oil degradation of bacterial communities, C transformation functions were promoted, including chemoheterotrophic, aerobic chemoheterotrophic, and aromatic_ compound degradation processes (Figure 3).Our findings demonstrated that as successive generations progress, the aerobic chemoheterotrophic process of G1 and G3 became significantly more active in G4.The aromatic degradation processes of G1 exhibited pronounced activation in G2, G3, and G4.Notably, no significant differences in chemoorganotrophic were observed across G1 to G4.
Alkanes primarily undergo intricate transformations facilitated by specific enzymatic systems such as alkane hydroxylases, ethanol dehydrogenases, acetaldehyde dehydrogenases, and acetyl-CoA.These processes culminate in their eventual entry into the fatty acid betaoxidation pathway, leading to their complete oxidation into carbon dioxide and water (Rojo, 2009).The degradation of polycyclic aromatic hydrocarbons (PAHs) involves a series of sequential reactions including hydroxylation, dehydrogenation, isomerization, and ring cleavage.This intricate sequence results in the generation of intermediates that become integrated into the tricarboxylic acid (TCA) cycle.Subsequently, these intermediates undergo enzymatic conversions orchestrated by microorganisms, ultimately leading to the production of carbon dioxide and water.Notably, the ring cleavage process constitutes a pivotal step in the degradation of PAHs, involving two distinct pathways (Varjani, 2017): (1) ortho-cleavage catalyzed by catechol 1,2-dioxygenase (C12O) at ortho positions, and (2) metacleavage catalyzed by catechol 2,3-dioxygenase (C23O) at meta positions.Utilizing the PICRUSt2 software, we acquired the abundance of genes encoding these enzymes.As illustrated in the Figure 4, the relative abundance of the alkane hydroxylase gene (K00496) in G1 notably surpassed that in G3 and G4, while both G2 and G3 exceed G4.The acetaldehyde dehydrogenase gene (K00001) in G1 and G3 also exhibited significantly higher abundance than in G4.Furthermore, the acetyl-CoA gene (K01897), along with the C12O gene (K03381) and C23O gene (K07104) involved in aromatic hydrocarbon metabolism, exhibits significantly higher abundance in G1, G2, and G3 compared to G4. Notably, the ethanol dehydrogenase gene (K00128) displays negligible variation between G1 and G4.Consequently, the trend observed across G1 to G4 generations unveils a marked decline in the relative prevalence of genes linked to alkane and aromatic hydrocarbon metabolism, excluding the ethanol dehydrogenase gene.

Functional decline of bacterial consortia
Based on key-score values and the relationship between genera, we constructed Figure 5. Pseudomonas exhibited the highest key-score  D).Their relative abundances were 18, 0.21, 0.03, and 1.77%, respectively (Supplementary Table S1).These findings indicate that the presence of key functional species in the consortium is not solely dependent on their relative abundance.Alternatively, important functional species can exist as either dominant populations with high relative abundance or non-dominant populations with low relative abundance.Table 1 showed the correlation between the genus of top key scores in G1-G4 and the degradation efficiency of TPH, saturates, and aromatics.The correlations ranged from 15.42 to −9.54 for TPH, 11.91 to −9.79 for saturates, and 8.03 to −5.96 for aromatics.
These results suggest a decline in the degradation capacity of TPH, saturated hydrocarbons, and aromatic hydrocarbons by the dominant species in the consortia.In Table 2, the highest positive correlation with TPH and saturated hydrocarbons degradation was observed for the Dietzia genus, with correlation values of 19.48 and 18.38, respectively (p < 0.05).Similarly, the Bacillus genus showed the highest positive correlation with aromatic hydrocarbons degradation, with a correlation value of 16.56 (p < 0.05).Conversely, in Table 3, Serratia exhibited the highest negative correlation with TPH and saturated hydrocarbon degradation, with correlation values of −19.38 and − 15.78, respectively (p < 0.05).Additionally, Psychrobacter showed the highest negative correlation with aromatic  Li et al. 10.3389/fmicb.2023.1270916Frontiers in Microbiology 07 frontiersin.orghydrocarbon degradation, with a correlation value of −8.92 (p < 0.05).These results emphasize the important roles of the Dietzia and Bacillus genera in the degradation of TPH, saturated hydrocarbons, and aromatic hydrocarbons, while highlighting the negative influence of Serratia and Psychrobacter on the degradation process.
The key scores of the Dietzia genus decreased from 3.62E-02 to 2.58E-02 during successive transfer cultures (Table 4).In contrast, the key scores of the Serratia and Psychrobacter genera increased from 0 to 2.87E-02 and 3.72E-02, respectively, after G2 (Table 4).The decline in the degradation of TPH and saturated hydrocarbons can be ascribed to the gradual weakening of the relationship between Dietzia and other genera in the consortia, as well as the progressive strengthening of the relationship between Serratia and Psychrobacter with other genera.These changes in inter-genera relationships within the consortia have significant implications for the overall degradation efficiency of TPH and saturated hydrocarbons.
The Bacillus genus, exhibiting the highest positive correlation (16.56) with aromatic hydrocarbon degradation, exhibit positive correlation with Mycetocola genus only in G1 (Figures 5A-D).This correlation pattern aligns with the observed degradation pattern of aromatic hydrocarbons during successive transfer cultures, suggesting that the relationship between these two genera plays a role in influencing the degradation of aromatic hydrocarbons.

Discussion
In environmental bioremediation, microbial consortia are more effective in pollutant degradation compared to single microorganisms (Gurav et al., 2017).Unlike the stochastic nature of the bottom-up approach, the top-down method facilitates the acquisition of degradation microbial consortia with heightened environmental competitiveness.Nevertheless, microbial consortia enriched through the top-down approach necessitate multiple generations for practical application.In this study, we employed the top-down approach to cultivate a stable and efficient microbial consortium for crude oil degradation.However, degradation efficiency and related functional potential diminish across generations.To obtain more steadfast and competitively adept microbial consortia, a thorough understanding of instability reasons through microbial interactions becomes imperative.To explore microbial interactions, researchers often employ correlation analysis or build regression linear models to construct co-occurrence networks (Matchado et al., 2021).In order to obtain more reliable species relationships, a common approach is to integrate multiple analysis methods, which has been extensively applied in plant research (Duran et al., 2018;Zhang et al., 2018) and soil microbiome study (Mandakovic et al., 2018).However, when studying species relationships in petroleum-degrading microbial communities, current practices have been limited to utilizing a single method to construct co-occurrence networks (Alvarez- Barragan et al., 2022;Xiao et al., 2022;Zhou et al., 2023).To infer more robust species interactions in the investigation of petroleum-degrading microbial communities, we adopt a method similar to that proposed by Mac Aogain et al. ( 2021), integrating five different correlation analysis methods to build a co-occurrence network.This integration approach alleviated the occurrence of spurious correlations among species.
During successive transfers, the abundance of the most critical genera in each generation did not exhibit the highest values, contrary to the conventional belief that species importance within consortia is solely dictated by their relative abundance.In fact, rare microorganisms have been found to exert a significant positive impact on the degradation capacity of consortia, as demonstrated by Delgado-Baquerizo et al. (2016).The Desulfosporosinus genus, despite its scant representation (0.006%) among the total 16S rRNA genes in the microbial community of peatlands, exhibits exceptional efficiency in sulfate reduction (Pester et al., 2010).These findings underscore the crucial role of rare species in providing the necessary genetic resources for the intricate degradation processes that occur within consortia.The constrained proliferation of these rare species is likely attributed to specific environmental conditions that impose unfavorable growth circumstances, as proposed by Jousset et al. (2017).
The strongest positive correlation was observed between the degradation rates of total petroleum hydrocarbons (TPH) and saturated hydrocarbons and the abundance of the Dietzia genus.Numerous studies have reported the hydrocarbon degradation capacity of the Dietzia genus (Wang et al., 2011;Venil et al., 2021).Comparative genomic analysis has revealed that the genome of this FIGURE 3 FAPROTAX analysis was utilized to evaluate the relative abundance of carbon transformation within G1-G4.Li et al. 10.3389/fmicb.2023.1270916Frontiers in Microbiology 08 frontiersin.orggenus contains a higher number of genes associated with lipid transport, metabolism, secondary metabolic production, synthesis, transport, and metabolism compared to other bacterial genomes (Fang et al., 2021).Lipid transport and metabolism play crucial roles in the degradation of components related to crude oil, while certain secondary metabolites, such as rhamnolipids, contribute to the emulsification of crude oil (Wang et al., 2014;Liu et al., 2018).Our study revealed an intriguing finding: the key scores of the Dietzia genus showed a gradual decline during successive transfers, indicating a weakening relationship between Dietzia and other genera in the consortium.However, when the Dietzia genus was co-cultured with other genera, it resulted in enhanced degradation of alkanes (Hu et al., 2020).This suggests that the declining relationship between Dietzia and other genera influences the efficiency of crude oil degradation.In conclusion, the relative abundance of the Dietzia genus in the consortia and its interaction with other genera play a significant role in the degradation of total oil and saturated hydrocarbons.The highest positive correlation was observed between the degradation rate of aromatic hydrocarbons and the abundance of the Bacillus genus, which is well-known for its capability to degrade aromatic hydrocarbons (Das and Mukherjee, 2007;Eskandari et al., 2017;Ghorbannezhad et al., 2022).Interestingly, our analysis revealed a significant positive correlation between the Bacillus genus and the Mycetocola genus, but only in the G1 generation of the bacterial consortia.This correlation aligns with the observed pattern of aromatic hydrocarbon degradation during successive transfers.Mycetocola genus has been reported to counteract the toxic effects of tolaasin I produced by the Pseudomonas genus, which inhibits the growth of both Gram-negative bacteria like Escherichia coli and Gram-positive bacteria like Bacillus subtilis (Rainey et al., 1991;Hermenau et al., 2020;Castaldi et al., 2022).Notably, in the G1-G4 bacterial consortia, the Pseudomonas genus was consistently present with a relative abundance above 7.93%.Therefore, we speculate that the lack of Mycetocola genus in the G2-G4 consortia could have resulted in the persistence of the toxic effect of tolaasin I, thereby limiting the degradation of aromatic hydrocarbons by the Bacillus genus.In summary, our findings suggest that the abundance of the Bacillus genus is positively correlated with the degradation of aromatic hydrocarbons.Furthermore, the presence of the Mycetocola genus appears to be crucial in alleviating the toxic effects of tolaasin I and facilitating the degradation of aromatic hydrocarbons by the Bacillus genus.During successive transfers, we observed that the degradation rates of total petroleum hydrocarbons (TPH) and saturated hydrocarbons exhibited the strongest negative correlation with the relative abundance of the Serratia genus.Similarly, the degradation rates of aromatic hydrocarbons showed the highest negative correlation with the relative abundance of the Psychrobacter genus.It is noteworthy that both the Serratia and Psychrobacter genera have been reported to possess crude oil degradation capabilities (Fagbemi and Sanusi, 2017;Lasek et al., 2017;Semai et al., 2021).By analyzing the relationships between the Serratia and Psychrobacter genera and other bacterial species, we noticed a gradual increase in the key scores of these genera within the bacterial consortia during successive transfers.This suggests a strengthening association between these genera and other members of the consortium.Considering the degradation capabilities of Serratia and Psychrobacter genera, it is reasonable to speculate that the Co-occurrence analysis of bacterial network along successive transfer cultures.Interactions between microbes (nodes) are represented by connecting lines (edges), and the key-scores are reflected by node size.(A-D) Visualization of positive and negative interactions between all taxa on G1, G2, G3, G4.Interactions between microbes are classified as red if the sign of the edge weights between them is positive (positive correlation) and vice versa, the intensity of the colors reflected the strength of the correlations.
TABLE 1 The relationship between the degradation of TPH, saturates, aromatics, and the genus of top key-scores.

Genus
The correlation coefficient with TPH degradation The correlation coefficient with saturates degradation The correlation coefficient with aromatics degradation competition among microbial species within the consortium contributed to a decline in the overall crude oil degradation capacity (Abtahi et al., 2020).
By employing network analysis methods, we can provide insights into the underlying causes of the functional decline in crude oil degradation.Nonetheless, to substantiate our research outcomes, it is

Genus
The correlation coefficient with TPH degradation

Genus
The correlation coefficient with saturates degradation

Genus
The  imperative to conduct subsequent experimental studies focusing on the pertinent microorganisms and genes.These additional investigations will contribute to a more comprehensive and robust understanding of the complex interactions within the microbial community and their implications for crude oil degradation.

Conclusion
The relevant functional potential and degradation capacity of the bacterial consortium for crude oil gradually declined during the successive transfers.To investigate the factors influencing the degradation function of the bacterial consortia, we employed the 16S rRNA amplification technique and conducted bioinformatics statistical analysis.Our results revealed that the relative abundance of key genera within the bacterial consortia was not the sole determinant of their importance.The relative abundance of the Dietzia genus and its interactions with other genera emerged as critical factors influencing the degradation of TPH and saturated hydrocarbons.Furthermore, the decreasing relative abundance of the Bacillus genus and its interaction with the Mycetocola genus were found to impact the degradation of aromatic hydrocarbons.These findings highlight the intricate and dynamic nature of microbial interactions within crude oil degradation processes, underscoring the importance of gaining a comprehensive understanding of microbial community dynamics for the development of effective bioremediation strategies.

FIGURE 1
FIGURE 1Degradation efficiency of crude oil and its two fractions in four successive transfers: (A) TPH degradation efficiency, (B) saturates degradation efficiency, and (C) aromatics degradation efficiency.The abbreviations used are CK for the control treatment and BT for the bacterial consortium treatment.The significance level is indicated by asterisks (*), with p < 0.05.

FIGURE 2
FIGURE 2 Bacterial taxonomic composition.(A) The Venn plot showing the shared and unique OTUs between G1-G4.(B) Core OTUs composition of four generations at the genus level.(C) Unique OTUs composition of four generations at the phylum level.(D) Taxonomic dendrogram showing the highly abundant bacteria of G1-G4.Color ranges identify order within the tree.Heatmap gradient represent the relative abundance of each OTU in G1-G4.(E) The relative abundance of bacterial phylum in the samples.(F) The relative abundance of the 10 most represent bacterial genera in the samples.Within each column, means with different letters are significantly different at p < 0.05.

TABLE 2
The strongest positive relationship between the degradation of TPH, saturates, aromatics, and the genera.
The significance level is indicated by asterisks (*), with p < 0.05.

TABLE 3
The strongest negative relationship between the degradation of TPH, saturates, aromatics, and the genera.The significance level is indicated by asterisks (*), with p < 0.05.