Study on the regulatory mechanism of the earthworm microbial community in vitro and in vivo under cadmium stress

Response patterns of the microbiome during hexavalent chromium remediation by Tagetes erecta L

Chromium pollution, particularly hexavalent chromium [Cr(VI)], may threaten the environment and human health. This study investigated the potential of Tagetes erecta L. (Aztec marigold) for phytoremediation of soil contaminated with Cr(VI), and focused on the effects of varying concentrations of Cr(VI) on both the physicochemical properties of soil and microbiome of Tagetes erecta L. We observed that Tagetes erecta L. showed tolerance to Cr(VI) stress and maintained normal growth under these conditions, as indicated by bioconcentration factors of 0.33-0.53 in shoots and 0.39-0.70 in roots. Meanwhile, the structure and diversity of bacterial communities were significantly affected by Cr(VI) pollution. Specifically, Cr(VI) had a more significant effect on the microbial community structure in the endophytic of Tagetes erecta L. than in the rhizosphere (p < 0.05). The genera Devosia and Methylobacillus were positively correlated with Cr(VI) concentrations. Biomarkers such as Bacilli and Pseudonocardia were identified under the different Cr(VI)-contaminated treatments using LEfSe. In addition, the interaction and stability of the endophytic microbiome were enhanced under Cr(VI) stress. This study explored the interactions between heavy metals, microorganisms, and plants, providing valuable insights for developing in situ bioremediation of Cr(VI)-contaminated soils.

Intensive ammonium fertilizer addition activates iron and carbon conversion coupled cadmium redistribution in a paddy soil under gradient redox conditions

While over-fertilization and nitrogen deposition can lead to the enrichment of nitrogen in soil, its effects on heavy metal fractions under gradient moisture conditions remains unclear. Here, the effect of intensive ammonium (NH4+) addition on the conversion and interaction of cadmium (Cd), iron (Fe) and carbon (C) was studied. At relatively low (30-80 %) water hold capacity (WHC) NH4+ application increased the carbonate bound Cd fraction (F2Cd), while at relatively high (80-100 %) WHC NH4+ application increased the organic matter bound Cd fraction (F4Cd). Iron‑manganese oxide bound Cd fractions (F3Cd) and oxalate-Fe decreased, but DCB-Fe increased in NH4+ treatments, indicating that amorphous Fe was the main carrier of F3Cd. The variations in F1Cd and F4Cd observed under the 100-30-100 % WHC treatment were similar to those observed under low moisture conditions (30-60 % WHC). The C=O/C-H ratio of organic matter in soil decreased under the 30-60 % WHC treatment, but increased under the 80-100 % WHC treatment, which was the dominant factor influencing F4Cd changes. The conversion of NH4+ declined with increasing soil moisture content, and the impact on oxalate-Fe was greater at 30-60 % WHC than at 80-100 % WHC. Correspondingly, genetic analysis showed the effect of NH4+ on Fe and C metabolism at 30-60 % WHC was greater than at 80-100 % WHC. Specifically, NH4+ treatment enhanced the expression of genes encoding extracellular Fe complexation (siderophore) at 30-80 % WHC, while inhibiting genes encoding Fe transmembrane transport at 30-60 % WHC, indicating that siderophores simultaneously facilitated Cd detoxification and Fe complexation. Furthermore, biosynthesis of sesquiterpenoid, steroid, butirosin and neomycin was significantly correlated with F4Cd, while glycosaminoglycan degradation metabolism and assimilatory nitrate reduction was significantly correlated with F2Cd. Overall, this study gives a more comprehensive insight into the effect of NH4+ on activated Fe and C conversion on soil Cd redistribution under gradient moisture conditions.

Exploring the Toxicity of Oxytetracycline in Earthworms (Eisenia fetida) Based on the Integrated Biomarker Response Method

Antibiotic contamination has become a global environmental issue of widespread concern, among which oxytetracycline contamination is very severe. In this study, earthworm (Eisenia fetida) was exposed to oxytetracycline to study its impact on the soil environment. The total protein (TP), catalase (CAT), superoxide dismutase (SOD), peroxidase (POD), malondialdehyde (MDA), glutathione S-transferase (GST), and glutathione peroxidase (GPX) oxidative stress indicators in earthworms were measured, and the integrated biomarker response (IBR) approach was used to evaluate the toxic effect of oxytetracycline on earthworms. A Technique for Order Preference by Similarity to an Ideal Solution (TOPSIS) and a path analysis model were used to explore the physiological and metabolic processes of earthworms after stress occurs. The results showed that SOD, GPX, and GST play important roles in resisting oxytetracycline stress. In addition, stress injury showed a good dose-effect relationship, and long-term stress from pollutants resulted in the most serious damage to the head tissue of earthworms. These results provide a theoretical basis for understanding the toxic effect of oxytetracycline on soil animals, monitoring the pollution status of oxytetracycline in soil, and conducting ecological security risk assessment.

Effects of biochar-derived dissolved organic matter on the gut microbiomes and metabolomics in earthworm Eisenia fetida

The ecological risks of biochar-derived dissolved organic matter (DOM) to soil invertebrates at different organismal levels remains limited. This study comprehensively explored the ecological risks of biochar-derived DOM on earthworm gut through assessments of enzyme activity response, histopathology, gut microbiomes, and metabolomics. Results demonstrated that DOM disturbed the digestive enzymes in earthworm, especially for 10% DOM300 groups. The integrated biomarker response v2 (IBRv2) indicated that the perturbation of earthworm digestive enzymes induced by DOM was both time-dependent and dose-dependent. Pathological observations revealed that 10% DOM300 damaged intestinal epithelium and digestive lumen of earthworms. The significant damage and injury to earthworms caused by DOM300 due to its higher concentrations of heavy metal ions and organic substrates (e.g., toluene, hexane, butanamide, and hexanamide) compared to DOM500 and DOM700. Analysis of 16S rRNA from the gut microbiota showed a significant decrease in genera (Verminephrobacter, Bacillus, and Microbacteriaceae) associated with inflammation, disease, and detoxification processes. Furthermore, 10% DOM300 caused the abnormality of metabolites, such as glutamate, fumaric acid, pyruvate, and citric acid, which were involved in energy metabolism, These findings contributed to improve our understanding of the toxic mechanism of biochar DOM from multiple perspectives.

Biotransformation of roxarsone by earthworms and subsequent risk of soil arsenic release: The role of gut bacteria

The organoarsenical feed additive roxarsone (ROX) is a ubiquitous threat due to the unpredictable levels of arsenic (As) released by soil bacteria. The earthworms representing soil fauna communities provide hotspots for As biotransformation genes (ABGs). Nonetheless, the role of gut bacteria in this regard is unclear. In this study, the changes in As speciation, bacterial ABGs, and communities were analyzed in a ROX-contaminated soil (50 mg/kg As in ROX form) containing the earthworm Eisenia feotida. (RE vs. R treatment). After 56 d, earthworms reduced the levels of both ROX and total As by 59 % and 17 %, respectively. The available As content was 10 % lower in the RE than in R treatment. Under ROX stress, the total ABG abundance was upregulated in both earthworm gut and soil, with synergistic effects observed following RE treatment. Besides, the enrichment of arsM and arsB genes in earthworm gut suggested that gut bacteria may facilitate As removal by enhancing As methylation and transport function in soil. However, the bacteria carrying ABGs were not associated with the ABG abundance in earthworm gut indicating the unique strategies of earthworm gut bacteria compared with soil bacteria due to different microenvironments. Based on a well-fit structural equation model (P = 0.120), we concluded that gut bacteria indirectly contribute to ROX transformation and As detoxification by modifying soil ABGs. The positive findings of earthworm-induced ROX transformation shed light on the role of As biomonitoring and bioremediation in organoarsenical-contaminated environments.

Inoculation of cadmium-tolerant bacteria to regulate microbial activity and key bacterial population in cadmium-contaminated soils during bioremediation

The perennial ryegrass Lolium perenne can be used in conjunction with cadmium (Cd)-tolerant bacteria such as Cdq4-2 (Enterococcus spp.) for bioremediation of Cd-contaminated soil. In this study, a theoretical basis was provided to increase the efficiency of L. perenne remediation of Cd-contaminated soil using microorganisms to maintain the stability of the soil microbiome. The experimental design involved three treatment groups: CK (soil without Cd addition) as the control, 20 mg·kg-1 Cd-contaminated soil, and 20 mg·kg-1 Cd-contaminated soil + Cdq4-2, all planted with L. perenne. The soil was collected on day 60 to determine the soil microbial activity and bacterial community structure and to analyze the correlation between soil variables, the bacterial community, available Cd content in the soil, Cd accumulation, and L. perenne growth. The soil microbial activity and bacterial community diversity decreased under Cd stress, and the soil microbial community composition was changed; while inoculation with Cdq4-2 significantly increased soil basal respiration and the activities of urease, invertase, and fluorescein diacetate (FDA) hydrolase by 83.65%, 79.72%, 19.88%, and 96.15% respectively; and the stability of the community structure was also enhanced. The Actinobacteriota biomass, the amount of available Cd, and the above- and belowground Cd content of L. perenne were significantly negatively correlated with the total phosphorus, total potassium, and pH. The activity of urease, invertase, and FDA hydrolase were significantly positively correlated with the biomasses of Acidobacteriota and L. perenne and significantly negatively correlated with the Chloroflexi biomass. Further, the available soil Cd content and the above- and belowground Cd levels of L. perenne were significantly positively correlated with the Actinobacteriota biomass and significantly negatively correlated with the Gemmatimonadetes biomass. Overall, inoculating Cd-tolerant bacteria improved the microbial activity, diversity, and abundance, and changed the microbial community composition, facilitating the remediation of Cd-contaminated soil by L. perenne.

Deciphering survival strategies: Oxidative stress and microbial interplay in Eisenia fetida under tetracycline contamination

Soil contamination resulting from residual antibiotics presents a pressing need to understand the survival mechanisms of soil organisms in polluted environments. This study focused on Eisenia fetida, and tetracycline stress experiments were conducted in a controlled environment using sterile artificial soil. The stress concentration ranged from 0 to 600 mg/kg, and stress cycles lasted either 10 or 30 days. The objective of this study was to assess the effects of oxidative stress and the changes in microbial communities both within and outside the earthworms. A comprehensive mathematical model was developed to elucidate the responses of organisms following exposure to stress utilizing factor analysis, grey relational analysis, and hierarchical entropy weight analysis. Under tetracycline stress, the initial stage (1-3 days) exhibited a coordinated regulation of oxidative stress and microbial communities in the soil with the assistance of CAT and GPX enzymes. The subsequent stage (4-5 days) further emphasized the influence of soil microbial communities. A notable "feedback regulation" of soil microbial communities on oxidative stress was observed during the third stage (6-8 days). Earthworms maintained a metabolic balance in the fourth stage (9-10 days). In the long term, the stress-induced a self-detoxification mechanism within soil microbial communities, which collaborated with GPX to respond to oxidative stress.

Response behavior of antibiotic resistance genes and human pathogens to slope gradient and position: An environmental risk analysis in sloping cultivated land

Soils, especially in farmlands, are key media for the transmission of antibiotic resistance genes (ARGs) and their hosts from the environment to humans. Sloping farmland is an important agricultural resource, but there lack of studies on the fate and risk of ARGs in sloping land. Also, the behavior and drivers of ARGs in response to slope gradient and position are unclear. Here, metagenomics was used to investigate the profiles of ARGs, mobile genetic elements, and microbial communities in soils from lands of five slope gradients (5°, 10°, 15°, 20°, and 25°) with two slope positions (uphill and downhill). Results showed that while the abundance (except 15°) and diversity (except 20°) of ARGs increased as the slope gradient increased, the diversity of ARGs with health risk, especially the high-risk ones, decreased. For slope positions, abundant and diverse ARGs were more likely to accumulate at downhill. Furthermore, 52 bacterial genera and 12 human pathogenic bacteria (HPB) species were identified as the potential hosts for ARGs with high risk, and abundant HPB species were also detected in the soils with low gradients at downhill. Moreover, the structural equation model analysis revealed that the slope gradient and the slope position have both direct and indirect effects on the abundance of ARGs. Further correlation analysis revealed that the slope gradient has a positive effect (p < 0.05) on nitrite nitrogen in the soils. Also, the slope position has a negative effect (p < 0.05) on total phosphorus and microbial nitrogen, while positively affected (p < 0.05) on particulate nitrogen and microbial carbon, which were the key factors driving the behavior of ARGs. Overall, this study provided comprehensive information on ARGs with health risks and their potential pathogenic hosts in sloping farmland. It can be important for controlling antibiotic resistance transmission and be consistent with the One Health framework.

The Relationship between Acid Production and the Microbial Community of Newly Produced Coal Gangue in the Early Oxidation Stage

Coal gangue is a solid waste formed during coal production, and the acid mine drainage it generates during open-pit storage severely pollutes the ecological environment of mining areas. Microorganisms play a crucial catalytic role in acidification, and their species and gene functions change during the oxidation process of coal gangue. In this study, the changes in microbial community structure were investigated during the initial acidification process for newly produced gangue exposed to moisture by monitoring the changes in pH, EC, sulfate ion concentration, and the iron oxidation rate of gangue leaching solutions. Moreover, the composition and functional abundance of microbial communities on the surface of the gangue were analyzed with rainfall simulation experiments and 16S rRNA sequencing. The study yielded the following findings: (1) The critical period for newly produced gangue oxidation spanned from 0~15 d after its exposure to water; the pH of leaching solutions decreased from 4.65 to 4.09 during this time, and the concentration and oxidation rate of iron in the leaching solutions remained at low levels, indicating that iron oxidation was not the main driver for acidification during this stage. (2) When the gangue was kept dry, Burkholderia spp. dominated the gangue microbial community. When the gangue was exposed to moisture, the rate of acidification accelerated, and Pseudomonas replaced Burkholderia as the dominant genus in the community. (3) In terms of gene function, the microbial community of the acidified gangue had stronger nitrogen cycling functions, and an increase in the abundance of microorganisms related to the sulfur cycle occurred after day 15 of the experiment. The microbial community in the acidified gangue had more stress resistance than the community of the newly formed gangue, but its potential to decompose environmental pollutants decreased.

Enhancement of carbon sequestration via MIL-100(Fe)@PUS in bacterial-algal symbiosis treating municipal wastewater

Bacterial-algal symbiosis (BAS) is a promising carbon neutrality technology to treat municipal wastewater. However, there are still non-trivial CO₂ emissions in BAS due to the slow diffusion and biosorption of CO₂. Aiming to reduce CO₂ emissions, the inoculation ratio of aerobic sludge to algae was further optimized at 4:1 on the base of favorable carbon conversion. MIL-100(Fe) served as CO₂ adsorbents was immobilized on polyurethane sponge (PUS) to increase the interaction with microbes. When MIL-100(Fe)@PUS was added to BAS in the treatment of municipal wastewater, zero CO₂ emission was achieved and the carbon sequestration efficiency was increased from 79.9% to 89.0%. Most genes related to metabolic function were derived from Proteobacteria and Chlorophyta. The mechanism of enhanced carbon sequestration in BAS could be attributed to both enrichment of algae (Chlorella and Micractinium) and increased abundance of functional genes related to PS I, PS II and Calvin cycle in photosynthesis.

Mapping and Identifying Candidate Genes Enabling Cadmium Accumulation in Brassica napus Revealed by Combined BSA-Seq and RNA-Seq Analysis

Rapeseed has the ability to absorb cadmium in the roots and transfer it to aboveground organs, making it a potential species for remediating soil cadmium (Cd) pollution. However, the genetic and molecular mechanisms underlying this phenomenon in rapeseed are still unclear. In this study, a 'cadmium-enriched' parent, 'P1', with high cadmium transport and accumulation in the shoot (cadmium root: shoot transfer ratio of 153.75%), and a low-cadmium-accumulation parent, 'P2', (with a cadmium transfer ratio of 48.72%) were assessed for Cd concentration using inductively coupled plasma mass spectrometry (ICP-MS). An F₂ genetic population was constructed by crossing 'P1' with 'P2' to map QTL intervals and underlying genes associated with cadmium enrichment. Fifty extremely cadmium-enriched F₂ individuals and fifty extremely low-accumulation F₂ individuals were selected based on cadmium content and cadmium transfer ratio and used for bulk segregant analysis (BSA) in combination with whole genome resequencing. This generated a total of 3,660,999 SNPs and 787,034 InDels between these two segregated phenotypic groups. Based on the delta SNP index (the difference in SNP frequency between the two bulked pools), nine candidate Quantitative trait loci (QTLs) from five chromosomes were identified, and four intervals were validated. RNA sequencing of 'P1' and 'P2' in response to cadmium was also performed and identified 3502 differentially expressed genes (DEGs) between 'P1' and 'P2' under Cd treatment. Finally, 32 candidate DEGs were identified within 9 significant mapping intervals, including genes encoding a glutathione S-transferase (GST), a molecular chaperone (DnaJ), and a phosphoglycerate kinase (PGK), among others. These genes are strong candidates for playing an active role in helping rapeseed cope with cadmium stress. Therefore, this study not only sheds new light on the molecular mechanisms of Cd accumulation in rapeseed but could also be useful for rapeseed breeding programs targeting this trait.

A method to study the effects of combined stress of cadmium and microplastics on the acute toxicity of Eisenia fetida

The compound pollutants formed by microplastics and cadmium present a significant potential threat to the soil-based ecosystem, and it is urgent to carry out relevant ecotoxicological studies. However, the lack of appropriate test methods and scientific mathematical analysis models has restricted the progress of research. Based on an orthogonal test design, a ternary combined stress test was performed to study the effect of microplastics and cadmium on earthworms. This study used the particle size and concentration of microplastics as well as the concentration of cadmium as test factors. Using the improved factor analysis model and the TOPSIS (Technique for Order Preference by Similarity to Ideal Solution) method, a new model was constructed according to the response surface methodology to analyze the acute toxic effects on earthworms under the combined stress of microplastics and cadmium. In addition, the model was tested in a soil-polluted environment. The results show that the model can perfectly integrate the spatiotemporal cross effects of the concentration and time of the applied stress, and the scientific data analysis process ensures the efficient development of ecotoxicological research in the actual compound pollution environment. Moreover, the results of the filter paper test and soil test showed that the equivalent toxicity ratio of cadmium concentration, microplastic concentration, and microplastic particle size to earthworms as 26:35:39 and 23:36:41, respectively. In terms of the interaction effect, a certain positive interaction was observed between the cadmium concentration and that of the microplastics and their particle size, while a negative interaction was observed between the concentration of microplastics and their particle size. This research provides a test basis and model reference for early monitoring of the health of contaminated soils and assessments of ecological safety and security.

Ecological adaptation of earthworms for coping with plant polyphenols, heavy metals, and microplastics in the soil: A review

In recent years, soil pollution by massive accumulation of heavy metals (HMs), microplastics, and refractory hydrocarbon chemicals has become an emerging and global concern, drawing worldwide attention. These pollutants influence soil diversity by hindering the reproduction, abundance, thereby affecting aboveground productivity. The scientific community has recently emphasized the contribution of earthworms to heavy metal accumulation, microplastic degradation, and the decomposition of organic matter in the soil, which helps maintain the soil structure. This review paper aimed to compile scientific facts on how earthworms cope with the effect of HMs, microplastics, and plant polyphenols so that vermiremediation could be widely applied for well-being of the soil ecosystem by environmentalists. Earthworms have special surface-active metabolites in their guts called drilodefensins that help them defend themselves against the oxidative action of plant polyphenols. They also combat the effects of toxic microplastics, and other oxidative compounds by elevating the antioxidant activities of their enzymes and converting them into harmless compounds or useful nutrients. Moreover, earthworms also act as biofilters, bioindicators, bioaccumulators, and transformers of oxidative polyphenols, microplastics, toxic HMs, and other pollutant hydrocarbons. Microorganisms (fungi and bacteria) in earthworms' gut of also assist in the fixation, accumulation, and transformation of these toxicants to prevent their effects. As a potential organism for application in ecotoxicology, it is recommended to propagate earthworms in agricultural fields; isolate, and culture enormously in industry, and inoculate earthworms in the polluted soil, thereby abate toxicity and minimizing the health effect caused by these pollutants as well enhance the productivity of crops.

New insights into the effects of chlorpyrifos on soil microbes: Carbon and nitrogen cycle related microbes in wheat/maize rotation agricultural field

Chlorpyrifos, a broad-spectrum organophosphorus insecticide, has been widely detected worldwide and is a potential neurotoxin and endocrine disruptor. Besides, chlorpyrifos has been proven that have a negative effect on soil microbes. In the present study, chlorpyrifos formulation (LORSBAN®, 45% emulsifiable concentrate) was applied in an agricultural field at the recommended dose (R dose, 270.0 and 337.5 g a.i. ha^-1 for wheat and maize respectively) and double recommended (DR) dose. Chlorpyrifos residue level and effect on soil microbes related to soil carbon and nitrogen cycle function were analyzed. Results showed that the half-lives of chlorpyrifos in wheat and maize field soil were 7.23-8.23 and 1.45-1.77 d, respectively. Application of chlorpyrifos at even DR dose did not result in unacceptable residual chlorpyrifos, where the final residual chlorpyrifos in wheat/maize (leaf, stem, and grain) was meet the requirement of the maximum residual limit (0.5 mg kg^-1 for wheat and 0.05 mg kg^-1 for maize) in China. Chlorpyrifos enhanced the activity of β-glucosidase by increasing the relative abundance of Sphingosinicella and promoted the carbon cycle in wheat field. The changes of cbbLR and cbbLG gene abundance also confirmed that chlorpyrifos could affect the import and export of soil carbon pool. The effect of chlorpyrifos on soil N cycle was determined by changes in the abundance of the bacterial genus Gemmatimonas, which is associated with denitrification. Further analysis of N-cycle functional genes and urease activity showed that chlorpyrifos inhibited nitrogen fixation in wheat field, but promoted nitrogen fixation in maize field. In general, bacterial abundance, urease, and AOA-amoA gene could be early warning markers of chlorpyrifos contamination. The results demonstrated the negative effects of chlorpyrifos on soil microbes especially on soil C and N cycle in actual agricultural field. It provides new insights about chlorpyrifos environmental pollution and its effect on soil ecosystems.

Comparison of gut microbiota between immigrant and native populations of the Silver-eared Mesia (Leiothrix argentauris) living in mining area

The complex gut bacterial communities have a major impact on organismal health. However, knowledge of the effects of habitat change on the gut microbiota of wild birds is limited. In this study, we characterized the gut microbiota of two different subspecies of the Silver-eared Mesia (Leiothrix argentauris), the native subspecies (L. a. rubrogularis) and immigrant subspecies (L. a. vernayi), using 16S rRNA gene high-throughput sequencing. These two subspecies live in a trace metal-contaminated area, and L. a. vernayi was trafficked. They are an excellent system for studying how the gut microbiome of wild animal changes when they move to new habitats. We hypothesized that the immigrant subspecies would develop the same adaptations as the native subspecies in response to habitat changes. The results showed that there were no significant differences in the composition, diversity, or functional metabolism of gut microbiota between native and immigrant subspecies under the combined action of similar influencing factors (the p values of all analyses of variance >0.05). In addition, the composition and functional metabolism of gut microbiota in two subspecies showed adaptation against trace metal damage. Linear discriminant analysis effect size (LEfSe) analysis revealed that Massilia in the intestinal microbiota of immigrant subspecies was significantly higher than that of native subspecies, suggesting that immigrant subspecies suffered habitat change. Finally, we found that these two subspecies living in the mining area had an extremely high proportion of pathogenic bacteria in their gut microbiota (about 90%), much higher than in other species (about 50%) living in wild environment. Our results revealed the adaptation of intestinal microbiota of immigrant Silver-eared Mesias under heavy metals stress, which would provide guidance for biodiversity conservation and pollution management in mining area.

Evolution of the earthworm (Eisenia fetida) microbial community in vitro and in vivo under tetracycline stress

Bioremediation of contaminated soil has received increasing attention, and the effects of antibiotic residues on the soil ecological environment are a current research hotspot. Earthworms are the first choice of soil organisms to indicate the degree of soil pollution, and their detoxification mechanism after antibiotic stress must be further explored. Taking Eisenia fetida as the research object, an antibiotic (tetracycline) stress test was carried out in sterile artificial soil. The stress concentrations were set at 0, 0.3, 3, 30, 300 and 600 mg/kg. The ECO method was used to cultivate microbes in earthworms and soil. The carbon source utilization intensity algorithm developed by our team was used for data statistics, and a factor analysis model was constructed to explore the succession process of microbes of earthworms in vivo and in vitro under tetracycline stress. The results showed that there were four processes in the evolution of microbes under short-term tetracycline stress: at 1-3 days, the microbes in worms played a leading role; at 4-5 days, the microbes in the worms and the soil microbes jointly resisted TET stress; after 6-8 days of stress, the microbes in worms still played the main role, but their role was weakened; and after 9-10 days, soil microbes played a leading role, and tolerant microbes appeared. Under long-term stress, the microbes of earthworms in vivo and in vitro were obvious different, and there may be no regulatory relationship. And the factor analysis model is suitable for the analyse of the changes in microbial communities in vivo and in vitro under TET stress. The research results provide a reference method and model basis for the bioremediation of antibiotic-contaminated soil and the study of earthworm detoxification mechanisms, and help agricultural development.