Next Article in Journal
First Report on Microplastics Quantification in Poultry Chicken and Potential Human Health Risks in Pakistan
Previous Article in Journal
Whole-Transcriptome Analysis on the Leaves of Rosa chinensis Jacq. under Exposure to Polycyclic Aromatic Hydrocarbons
Previous Article in Special Issue
Assessing the Impact of Microplastic Filaments Contaminated with PAHs on Mytilus coruscus Larvae through Surface Contact
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Metagenomic and Metatranscriptomic Analyses Reveal the Response of Black Soldier Fly (Hermetia illucens) Larvae Intestinal Microbes and Reduction Mechanisms to High Concentrations of Tetracycline

1
Key Laboratory of Agricultural Microbial Enzyme Engineering (Ministry of Agriculture), School of Life Sciences, Henan Agricultural University, Zhengzhou 450046, China
2
Department of Infectious Diseases and Public Health, Jockey Club College of Veterinary Medicine and Life Sciences, City University of Hong Kong, Kowloon, Hong Kong
*
Authors to whom correspondence should be addressed.
Toxics 2023, 11(7), 611; https://doi.org/10.3390/toxics11070611
Submission received: 11 June 2023 / Revised: 10 July 2023 / Accepted: 10 July 2023 / Published: 13 July 2023

Abstract

:
Black soldier fly (Hermetia illucens L.) larvae (BSFL) possess remarkable antibiotic degradation abilities due to their robust intestinal microbiota. However, the response mechanism of BSFL intestinal microbes to the high concentration of antibiotic stress remains unclear. In this study, we investigated the shift in BSFL gut microbiome and the functional genes that respond to 1250 mg/kg of tetracycline via metagenomic and metatranscriptomic analysis, respectively. The bio-physiological phenotypes showed that the survival rate of BSFL was not affected by tetracycline, while the biomass and substrate consumption of BSFL was slightly reduced. Natural BSFL achieved a 20% higher tetracycline degradation rate than the germ-free BSFL after 8 days of rearing. Metagenomic and metatranscriptomic sequencing results revealed the differences between the entire and active microbiome. Metatranscriptomic analysis indicated that Enterococcus, Vagococcus, Providencia, and Paenalcaligenes were the active genera that responded to tetracycline. Furthermore, based on the active functional genes that responded to tetracycline pressure, the response mechanisms of BSFL intestinal microbes were speculated as follows: the Tet family that mediates the expression of efflux pumps expel tetracycline out of the microbes, while tetM and tetW release it from the ribosome. Eventually, tetracycline was degraded by deacetylases and novel enzymes. Overall, this study provides novel insights about the active intestinal microbes and their functional genes in insects responding to the high concentration of antibiotics.

1. Introduction

Antibiotics are widely used in medicine, agriculture, and animal husbandry [1], with China producing more than 200,000 tons of antibiotics yearly [2]. The annual global consumption of antimicrobial drugs is also increasing and is projected to grow by 104,000 tons by 2050 [3]. However, due to the low efficiency of tetracycline biosorption by animals and humans, 50–80% of tetracycline is released into the environment, leading to the spread of multiple drug-resistance genes and significant toxicity into the ecosystem [4]. Physical and chemical techniques, such as electrochemical methods and the adsorption of porous materials, are being used to remove antibiotics in the environment [5]. However, these techniques are costly and tend to cause secondary contamination. In contrast to traditional methods, bioremediation techniques can effectively remove antibiotics [6]. Currently, many studies have been conducted to design strategies that remove antibiotics via efficient degradation of strains, including Klebsiella sp.TR5 [7], Trichosporon mycotoxinivorans XPY-10 [8], Stenotrophomonas maltophilia DT1 [9], Serratia marcescens WW1 [10], Arthrobacter nicotianae OTC-16 [6], and Alcaligenes sp. T17 [11]. However, most high-efficient antibiotic degradation microbes are pathogenic and risk spreading antibiotic resistance genes. Moreover, these strains struggle to consistently coexist with the indigenous microbial community [6,12]. Therefore, novel approaches are needed to remediate antibiotic contamination.
Recently, there has been growing interest in insects due to their potential ability to treat different types of pollutants. This ability is attributed to their wide range of feeding habits, high economic value, easy management, and low risk of antibiotic resistance gene transmission [13]. Among insects, black soldier fly larvae (BSFL), also known as Hermetia illucens (Diptera: Stratiomyidae), have shown high capacity to degrade most antibiotics, which is primarily credited to its powerful intestinal flora [14,15]. For instance, Yang et al. discovered that BSFL effectively degraded 84.22% of ciprofloxacin in poultry manure within 12 days [16]. Moreover, they isolated Klebsiella pneumoniae CIP1 and Proteus mirabilis CIP5, which can efficiently degrade more than 84% of 50 mg/kg ciprofloxacin in 12 days. In addition, Luo et al. achieved an 84.9% reduction in lincomycin in fermentation residue via BSFL treatment within 12 days [17]. Acinetobacter and Sphingobacterium are presumed to be the dominant groups involved in lincomycin degradation. Furthermore, Liu et al. achieved 63.3–82.7% oxytetracycline degradation at an initial concentration of 100–3042.3 mg/kg in 7–8 days after BSFL treatment [5,18]. Moreover, they identified the remarkable alteration of the relative abundances of lgnatzschineria, Enterococcus, and Providencia.
To discover the main functional genes, molecular methods were used to detect the abundance/expression level of related functional genes in previous works. For instance, Luo et al. explored the expression levels of six lincomycin resistance genes (lnuA, lnuB, lnuC, lnuD, lnuE, and lnuF) via qPCR [17]. Similarly, Liu et.al selectively identified 180 or 285 ARGs in the gut microbial of BSFL under the pressure of oxytetracycline using high-throughput qPCR (HT-qPCR) [5,18]. The results clearly showed that the abundances of tetA, tetB, and tetC increased. However, almost all existing studies on the BSFL gut flora and functional genes that respond to tetracycline were based on a combination of MiSeq and qRT-PCR analyses, which are limited in providing a comprehensive understanding of the overall environment in terms of the microbial community and all of its functional genes changes.
Metagenomic analysis is a field of study that examines the genome DNA (gDNA) of environmental microbial communities, revealing the functional genes that are present in these communities [19]. In contrast to metagenomics, metatranscriptomics analyzes the mRNA of microbial communities in the environment to reveal the active micro-organisms and their functional gene resources that respond to changes in environmental factors [20]. However, there is still a lack of studies related to the changes in BSFL gut flora composition and functional genes recorded via metagenomic and metatranscriptomic analysis techniques. Furthermore, previous work mainly focused on the degradation of BSFL to low concentrations of antibiotics. The effects of high concentrations of antibiotics on gut microbes and functional genes of BSFL are unclear.
In this study, we aimed to explore the changes in gut microbial community structure and functional genes at the genomic and transcriptomic levels in response to high tetracycline stress in BSFL via combined comparative metagenomic and metatranscriptomic analyses. The objectives of the study included determining the changes in the physiological characteristics of BSFL and the gut microbial community in response to high tetracycline presence, as well as the changes in the abundances of associated antibiotic resistance genes and degradation genes, thus inferring potential mechanisms of antibiotic degradation used by BSFL. The findings of this study provide a theoretical basis for the application of BSFL to dispose of tetracycline-contained contaminants.

2. Results

2.1. Effect of Tetracycline on the BSFL Growth

To assess the effect of tetracycline on the growth and substrate utilization of BSFL, a high concentration of tetracycline (1250 mg/kg) was added to the larvae’s diet. The results showed that tetracycline did not affect the survival of BSFL during the 8-day feeding time (Figure 1a). In addition, after 8 days of rearing, the BSFL experienced a modest decrease in biomass and substrate consumption rate of 1.20 g/10 larva and 63.2%, respectively, compared to the control group without antibiotic supplementation (1.53 g/10 larva and 73.5%) (Figure 1b,c). Meanwhile, the natural BSFL group achieved a tetracycline degradation rate of 32.21%, which was far higher than those of CK (9.92%) and the germ-free model groups (11.51%) (Figure 1d). These results indicated that BSFL has excellent tolerance of the high concentration of tetracycline, and its intestinal flora plays an essential role in the tetracycline degradation process.

2.2. Overview of Metagenomic and Metatranscriptomic Raw Data

Metagenomic and metatranscriptomic sequencing methods were used to analyze the microbial community and functional genes involved in the process of tetracycline degradation. Raw reads of 113,499,814, 121,370,872, 115,167,960, and 117,241,914 were obtained for the CK-MG, TC-MG, CK-MT, and TC-MT groups, respectively. The sequencing data generated 41,586, 43,057, 9526, and 2500 contigs from 78.55%, 85.77%, 92.77%, and 80.61% clean reads, respectively. Additionally, the ratio of clean bases to raw bases for each group was 77.91%, 85.09%, 92.20%, and 80.20%, respectively. Clean reads Q20 and Q30 represent the percentage of bases with a Phred value of >20 or 30. The Q20 and Q30 statistics in all samples exceeded 97.4% and 94.0%, respectively (Table 1). The N50 and N90 for MG data were ≥4314 bp and 650 bp, respectively, while for MT data, they were ≥1078 and 565 bp, respectively. All of the above results implied that the MG and MT data were of the high sequencing quality required for further analysis.

2.3. Composition and Diversity of Microbial Communities in BSFL Gut

Metagenomic sequencing data revealed that the microbial community structure was largely changed at the phylum level in TC group compared to CK group. Proteobacteria was the most dominant phylum, having a relative abundance of 75.98% and 52.65% in the CK and TC groups, respectively, followed by Bacteroidetes and Firmicutes (Figure 2a). Notably, the relative abundance of Bacteroidetes increased from 5.16 to 30.84% under the stress of tetracycline. At the genus level, the most abundant microbes remained the same, but their abundances showed considerable distinctions between the TC and CK groups. Ignatzschineria and Providencia were the most abundant genera in the CK and TC groups, accounting for 38.78% and 14.98%, respectively (Figure 2b). Meanwhile, the addition of tetracycline increased the percentage of Moheibacter and Sphingobacterium from 0.10 to 13.07% and from 1.35 to 5.89%, respectively. These results indicated that the indigenous intestinal microbiome of BSFL was influenced by tetracycline.
The results of metatranscriptomic analysis were markedly different from the metagenomic results. Firmicutes was the most dominant phylum, having a richness of 71.20% and 71.10% in CK and TC groups, respectively (Figure 2c), followed by Proteobacteria and Actinobacteria. The relative abundance of Actinobacteria was decreased from 9.76% in the CK group to 6.70% in the TC group. At the genus level, the relative abundances of the dominant microbes, including Enterococcus, Vagococcus, Providencia, and Paenalcaligenes, were increased with TC treatment (Figure 2d). The above results demonstrate the differences between the total and active microbes in BSFL intestine.

2.4. Antimicrobial Resistance (AMR) Gene Analysis

Functional genes of BSFL intestinal microbes at genetic and expressed levels were revealed via annotation of resistance genes in the CARD database. At the DNA level, CK and TC groups contained 707 and 711 AMR genes, respectively. The overall effect of tetracycline on the types and abundances of antibiotic resistance genes in BSFL gut microbes was small (Figure 3a–d). Specifically, the abundances of genes associated with ATP-binding cassette (ABC) antibiotic efflux pump and major facilitator superfamily (MFS) antibiotic efflux pump were slightly decreased, whereas the abundances of genes related to tetracycline-resistant ribosomal protection protein and tetracycline inactivation enzyme were slightly elevated by 0.22% and 0.30%, respectively (Figure 3a,b). The abundance of the tetA(48), which was the most abundant tetracycline resistance gene in both groups, decreased from 44.58% in the CK group to 41.06% in the TC group (Figure 3c,d). Moreover, the abundance of tetB(48) was decreased from 3.82 to 0.89%, whereas the abundances of tetB(60) and tet(X) were elevated from 6.85% and 0.60% to 9.61% and 3.89%, respectively, after tetracycline treatment. These findings suggest that tetracycline has a minor influence on the abundance of tetracycline resistance genes in BSFL gut microbes.
In addition, changes in antibiotic resistance genes at the transcriptional level were analyzed. The CARD database showed that 153 and 118 AMR genes were detected at the RNA level in the CK and TC groups, respectively. Moreover, the genes’ expression of ABC antibiotic efflux pump, MFS antibiotic efflux pump, and tetracycline-resistant ribosomal protection protein were upregulated from 7.63%, 1.82%, and 0.95% to 8.00%, 3.60%, and 2.31%, respectively (Figure 4a,b). However, for tetracycline inactivation enzymes, the tet(X) family was not annotated. These results indicated that there must be other enzymes in the BSFL gut that presented as tetracycline-degrading enzymes. Furthermore, the analyses of the tetracycline resistance genes revealed that tet(31), tetA(60), tetS, and tetB(46) were the most dominant genes in CK group, with expression levels decreasing to 0.19%, 10.64%, 23.81%, and 10.05% in TC group. Further, the abundances of tet(33), tet(A), tet (B), tet(L), tetA(60), tetR, tetM, tetW, and tet(L) were increased in the TC group compared to the CK group (Figure 4c,d). The above results suggested that tetracycline has effects on the expression of the genes of the overall BSFL gut microbes. Further, these results also indicate that genes that respond to antibiotic stress in the BSFL gut flora vary greatly in terms of the copy numbers and expression levels.

2.5. Differentially Expressed Genes in Metatranscriptomic Library

To further investigate the influence of tetracycline on the expression of all genes in the BSFL intestinal microbe, the differentially expressed genes between two groups were analyzed. It showed that 1159 genes were upregulated, whereas 4980 genes were downregulated, under tetracycline stress (Figure 5a). KEGG enrichment analysis showed that among the top 15 functional gene enriched pathways, D-Glutamine and D-glutamate metabolism underwent significant upregulation (p = 0.03), with the gene ratio reaching 0.77 (Figure 5b), which is the only ratio reported that is related to tetracycline resistance [21].

2.6. Potential Tetracycline Degradation Genes

The COG database was utilized to annotate the functional genes to analyze potential tetracycline degradation-related genes. The metagenomic data revealed that the diversity and abundances of functional genes in the overall BSFL intestinal microbiota remained unchanged under tetracycline pressure (Figure 6a). In contrast, metatranscriptomic data exhibited that genes belonging to the functional category of “translation, ribosomal structure and biogenesis” increased from 19.02% (CK) to 25.66% (TC) (Figure 6b). This result indicated that tetracycline has an obvious effect on the translation and expression of genes, and the upregulated genes are likely to be related to the degradation of tetracycline. Further, the above results implied that the functional genes of BSFL gut exhibited a positive transcriptional response to tetracycline, but did not cause any changes in their gene copy numbers.
The distribution of functional genes based on the CAZy database was similar to that of the COG analysis. There was a remarkable difference in the abundances of annotated carbohydrates between the genomic and transcriptomic levels. Further, while functional gene distribution was not distinctly different between the two groups at the genomic level, it showed quite differences at the transcriptomic level (Figure 7a,b). At the DNA level, glycosyl transferases (GTs, 37.68–41.64%) and glycoside hydrolases (GHs, 30.08–31.14%) were the most two dominant groups, followed by carbohydrate esterases (CEs, 15.57–17.31%) and carbohydrate-binding modules (CBMs, 7.54–8.71%) (Figure 7a). In contrast, the gene abundances of both CEs and CBMs were increased at the transcriptional level, from 12.51 to 14.26% and from 23.91 to 35.55%, respectively, while GTs, GHs, and AAs decreased in the TC group in comparison to the CK group (Figure 7b). These results illustrate that tetracycline affects not only the abundances of these genes, but also their expression level.
Tet (X), horseradish peroxidase, manganese peroxidase, glutathion S-transferase, lignin peroxidase, and laccase are known as the tetracycline degradation enzymes, and their corresponding genes were identified via metagenomic and metatranscriptomic libraries, respectively. Genes involved in antibiotic degradation, specifically AA&CE, were annotated and are presented in Table 2. Notably, AA1 (laccase), AA2 (manganese peroxidase), AA7 (esterase), and AA10 (monooxygenases), which have been reported to be involved in tetracycline degradation, were identified in the metagenomic database, but were not identified in the metatranscriptomic database. This result suggested that the efficient tetracycline degradation ability of the BSFL gut microflora was not caused by known tetracycline-degrading enzymes. Therefore, we speculated that novel and efficient tetracycline-degrading enzymes may exist in the BSFL gut microflora. Moreover, CEs that were in charge of degrading other antibiotics presented distinctive differences between metagenomic and metatranscriptomic levels. The expression of genes corresponding to CE1, CE3, CE4, CE6, CE7, CE8, and CE11 was upregulated after tetracycline treatment (Table 2), which indicated that they may be involved in tetracycline degradation.

2.7. Possible Mechanisms of BSFL Intestinal Microbiome Response to High Tetracycline Concentration Stress

Under high levels of TC stress, the BSFL’s active gut microbial community structure and expression of related functional genes were changed. Based on the results presented above, mechanisms of the BSFL intestinal microbiota response to the high tetracycline concentrations were supposed, as shown in Figure 8. Firstly, Tet family genes, including tet(33), tet(A), tet(B), tet(L), and tetA(60), act as efflux pumps to remove tetracycline from the cell. In addition, tetM and tetW may dissociate tetracycline that has already bound to the ribosome, thus promoting the diffusion of tetracycline in microbes, as well as accelerating subsequent degradation. Notably, other known tetracycline-degrading enzyme genes, including Tet(X), horseradish peroxidase, manganese peroxidase, lignin peroxidase, glutathione-S transferase, and laccase, were not detected at the transcriptomic level. Deacetylase (CE3, CE4, CE7, CE8, and CE11) was identified as the only type of tetracycline-degrading enzyme used in this study, which assisted in the hydrolysis of tetracycline and, ultimately, its removal. Therefore, it is concluded that there may be novel and highly efficient tetracycline-degrading enzymes present in the BSFL’s intestinal microbes.

3. Discussion

BSFL has been utilized for the degradation of various antibiotics [17,22,23]. In this study, we investigated the ability of BSFL to withstand and degrade a high concentration of exogenous tetracycline (1250 mg/kg). The results showed that BSFL maintained a 100% survival rate for 8 days after exposure to the high concentration of tetracycline (Figure 1a), indicating their strong ability to resist this antibiotic. Although the biomass and substrate consumption rate decreased by only 0.33 g/10 larva and 10%, respectively (Figure 1b,c), the BSFL intestinal microbiome still demonstrated a significant ability to degrade tetracycline. Intestinal flora assisted BSFL in achieving a tetracycline degradation rate of 32.21%, which was higher than that of the germ-free group (11.51%) (Figure 1d). This result contrasts with those of a previous study by Cai et al., which showed that normal BSFL achieved a 95.5–96.9% degradation rate at low initial concentrations of 20–80 mg/kg tetracycline, while sterile BSFL reduced tetracycline content by ~60% [22]. The difference between each study’s results could be due to the fact that the high tetracycline concentration inhibits the activity of some micro-organisms, thus increasing the time required to achieve a higher degradation rate. Moreover, Liu et al. found that BSFL could degrade substrates containing oxytetracycline (initial concentration of 434.4 mg/kg) by 65.9% within 8 days [5]. These results suggest that BSFL can effectively degrade high concentrations of tetracycline, with the assistance of BSFL intestinal flora, without affecting their own growth and survival.
The composition of the gut microbiota at the phylum level revealed that Proteobacteria, Bacteroidetes, and Firmicutes were the dominant microbes at the genetic level (Figure 2a), which is consistent with a previous study that found that these three phyla account for more than 99% of the BSFL’s gut microbial community [22]. Moreover, Yun discovered that Proteobacteria and Firmicutes constitute a major part of the gut microbiota of 218 different insect species [24]. After feeding the tetracycline-containing substrate for 8 days, the relative abundances of Proteobacteria and Firmicutes in the BSFL gut flora decreased by 23.33% and 3.39%, respectively, while that of Bacteroidetes increased by 25.68%. The results of this study were in accordance with a previous study’s gDNA-based Miseq sequencing results [22]. This consistency indicates that gDNA-based analysis tools differ only slightly in terms of the relative abundances of different phyla, but not in terms of the composition and abundance changes in the microbial community.
Interestingly, at the transcriptional level, the microbial structure presented differently to that at the genomic level, suggesting differences in the composition of the active and entire microbiomes. Furthermore, the microbial community did not show huge changes between two groups. In both the CK and TC groups, Firmicutes was the dominant phylum, followed by Proteobacteria, Actinobacteria, Arthropoda, and Bacteroides (Figure 2c). Previous studies indicated that Firmicutes and Proteobacteria carry a variety of tetracycline resistance genes [25]. The relative abundance of Actinobacteria was decreased, while that of Bacteroides was increased. This result is consistent with a study that shows a strong correlation between Bacteroides and the degradation of tetracycline in disposed pig manure by housefly larvae [26]. Ji et al. also discovered that the relative abundance of Bacteroides in the human intestine rose from 1.68% to 5.70% under tetracycline exposure [27].
At the genus level, the dominant microbes in BSFL gut also presented differently based on gDNA and mRNA analysis (Figure 2b,d). The relative abundances of active genera, including Enterococcus, Providencia, Paenalcaligenes, and Carnobacterium, was increased via tetracycline treatment (Figure 2d), which indicated that they may correlate with the tetracycline-degrading ability of BSFL. Enterococcus is frequently isolated from tetracycline-contaminated environments and considered to be a reservoir of tetracycline and macrolide–lincosamide–streptogramin (MLS) resistance genes [28]. Providencia is one of the core bacteria in the BSFL gut and believed to be responsible for the tetracycline-modifying gene tetX [29]. Paenalcaligenes were also reported to strongly respond to antibiotic exposure [30]. Interestingly, Carnobacterium was not previously known to possess resistance to tetracycline [31]. This shift may be due to gene horizontal transfer, which allowed it to acquire the tetracycline resistance genes.
Similar to the results of structural analysis of the BSFL gut microbiota, their functional genes differed distinctly between the transcriptomic and genomic levels (Figure 3 and Figure 4). Metatranscriptomic analysis provided a comprehensive understanding of the whole gene expression and regulation of BSFL gut microbiome in response to tetracycline. Firstly, the KEGG annotation results showed a significant upregulation of genes associated with D-Glutamine and D-glutamate metabolism under tetracycline stress (Figure 5b). Zhao et al. reported that glutamine modulates the permeability of bacterial cell membranes, which, in turn, affects the amount of antibiotics that enter the bacteria [21]. Furthermore, the CARD database comparison results illustrated the expression differences of antimicrobial resistance genes before and after tetracycline treatment. Genes that express ABC-binding antibiotic efflux pump, major facilitator superfamily antibiotic efflux pump, and tetracycline-resistant ribosomal protein families were also upregulated (Figure 4a,b). The efflux pump is one of the main means through which bacteria develop drug resistance [32]. Notably, Tet family genes, including tet(33), tet(A), tet(B), tet(L), tetA(60), tetR, tetM, and tetW, showed distinct differences at expressed levels (Figure 4c,d). Among these genes, tet(33), tet(A), tet(B), tet(L), and tetA(60) are mainly involved in the expression of the efflux pump [33,34]. TetR is a repressor of the tetracycline resistance element, the N-terminal region of which forms a helix-turn-helix structure that binds DNA. Binding of tetracycline to tetR reduces the repressor’s affinity for the tetA promoter–operator sites [35]. tetM was reported to induce the translocation of ribosome-bound tetracycline to the free state and shorten the half-life of the drug in most gram-positive bacteria [36]. TetW was found to be widespread in the BSFL gut community in the previous study, which blocked the action of tetracycline through ribosome non-covalent modifications [22,37]. Meanwhile, we were surprised to find that tet(X), which has tetracycline degradation ability, was not expressed in either the CK or TC group, which is in contrast to the metagenomic results. Therefore, we speculated that the BSFL’s intestinal microbiota selectively express certain genes to maintain a balance among them in order to resist the external high concentration of tetracycline stimulation.
Drug resistance genes are not the only factors that aid the host in resisting tetracycline stress, as some proteins with hydrolytic and redox functions may also be involved in tetracycline degradation. The differentially expressed genes annotated in the COG database that illustrated active functional genes under the “Translation, ribosomal structure and biogenesis” category were elevated with tetracycline treatment (Figure 6b). This observation is in line with a previous study by Liu et. al., which used a knockout approach to find that translation elongation factor 4 (LepA) affects the ribosome biogenesis of cells under the stress of tetracycline [38]. Results of the CAZy database annotation showed that genes that enable tetracycline-degrading enzymes that were present under the AA module were detected at the metagenomic level (Figure 7). Specifically, AA1 gene family that annotated for laccase, and AA2 gene family that annotated for manganese peroxidase and lignin peroxidase, were identified. All three types of protein were previously reported to possess the ability to degrade tetracycline [12,39,40,41,42]. However, no regulatory changes in the genes under the AA module occurred at the RNA-seq level. This result suggests that there may be other metabolic pathways involved in tetracycline degradation.
In this study, seven annotated gene families (CE1, 3–4, 6–8, 11) were found to be upregulated after tetracycline treatment (Table 2). CEs were annotated as esterase, which were reported to degrade β-lactam and macrolide antibiotics by hydrolyzing ester bonds [43]. This finding suggests that a correlation exists between different antibiotic resistance levels. Meanwhile, deacetylases classified as CE3, CE4, CE7, CE8, and CE11 deserve our attention due to their roles in the hydrolysis of acetamides [44], which may be involved in tetracycline degradation via deamination. Furthermore, 171 unknown function genes were upregulated in our study and may be involved in tetracycline degradation. Therefore, the subsequent study will be focused on exploring the unknown efficient tetracycline-degrading enzymes.

4. Conclusions

This study employed a combined comparative metagenomic and metatranscriptomic approach to investigate alterations in the BSFL intestinal microbiome under a high concentration of tetracycline pressure. BSFL achieved an impressive 32.21% tetracycline reduction in 8 days under the 1250 mg/kg tetracycline concentration. Through microbiome sequencing, we discovered that tetracycline does indeed affect the overall ecological community’s structure and abundances in the gut, though the composition of active genera at the transcriptional level remained largely unchanged. Through the annotation of the functional genes in several databases, such as eggNOG, CAZy, and CARD, the abundance changes in potential active tetracycline-degrading genes in the BSFL gut microbiome were examined. Further, the superior tetracycline resistance mechanisms of BSFL gut microbes have been confirmed. The findings of this study provide a valuable reference for future investigations into insects’ in situ bioremediations of antibiotic contamination. Our future research will be focused on exploring the novel and highly efficient tetracycline degradation enzymes.

5. Materials and Methods

5.1. BSFL Husbandry and Substrate

The BSFL utilized in this study were acquired from the culturing center of Henan Agriculture University. The 3rd instar larvae were reared at a moisture content of 70% and a temperature of 30 °C based on the protocol established in our previous study [45]. Next, they were divided into two groups, including the CK (fed soybean meal only) and TC (fed soybean meal containing a final concentration of 1250 mg tetracycline per kg dry weight) groups. Deionized water was added to the substrates daily to maintain a moisture content close to 70%. The tetracycline used in this study was purchased from MCE (10 g, CAS No.: 60-54-8; MedChemExpress, Monmouth Junction, USA).
The survival rate, biomass, and substrate consumption rate were determined as follows:
  • BSFL’s survival rate (%) = (N2/N1) × 100%, where N1 and N2 are the initial and final number of BSFL in each group, respectively. Ten larvae were randomly selected and weighed on the 8th day of rearing to determine their biomass.
  • BSFL’s substrate consumption rate (%) = [(W1 − W2)/W1] × 100%, where W1 and W2 denote the initial and final dry weights of the wheat bran substrate during 8 days of rearing, respectively.

5.2. Germ-Free and Sterile BSFL Model Construction

The effect of tetracycline on the growth of larvae and the ability of non-sterile BSFL to degrade tetracycline was investigated as follows: A germ-free BSFL model was constructed according to our previous research [45]. Fresh BSFL eggs weighing approximately 0.4 g were washed with 1 mL of 2.7% NaClO and 1 mL of Sporgon solution (Beijing Mingyangkehua Bio-Technology Co. Ltd., Beijing, China) and rinsed with sterile water to obtain germ-free BSFL. The sterilized eggs were inoculated with brain heart infusion (BHI) medium and further incubated at 37 °C for 24 h. The cultures were subsequently spread on BHI agar plates to verify the sterilization effect.
Autoclaved wheat bran was moistened with 100 mL of distilled H2O containing tetracycline to achieve a final concentration of up to 1250 mg/kg of dry weight. Approximately 50 germ-free or sterile BSFL were introduced into each flask via the following treatments: (i) tetracycline-spiked substrate and non-sterile larvae, (ii) tetracycline-spiked substrate and sterile larvae, and (iii) tetracycline-spiked substrate and no larvae. The experiments were carried out over 8 days, during which a 1-gram sample of the substrate was collected daily and stored at −20 °C to measure the concentration of tetracycline.

5.3. Tetracycline Determination

Tetracycline concentration was determined using the HPLC method. Samples underwent freeze-drying, crushing using a ceramic mortar, and screening through a 100-mesh sieve. Subsequently, 1 g of sample from each group was mixed with 100 μL of the internal standard solution (200 mg/L tetracycline) in a centrifuge tube. Next, 13 mL of phosphate buffer (pH = 3) and 2 mL of methanol solution were added into the tubes. After 1 min of vertexing, the tubes were sonicated for 20 min and centrifuged at 7500 rpm for 1 min. These steps were repeated 3 times to collect the supernatant, which was subsequently filtered via a 0.22-micrometer filter [22]. A high-performance liquid chromatograph (Agilent Technologies) equipped with an octadecyl X-Terra® C18 hybrid silica column (250 mm × 4.6 mm, 5 μm) was used to analyze the concentration of tetracycline. A series of tetracycline concentrations (0, 0.1, 0.5, 1, 2.5, 5 mg/L) were diluted from the standard solution (500 mg/L tetracycline) and used to draw a standard curve. The aqueous phase was composed of sodium EDTA (0.025 mol/L), calcium chloride (0.035 mol/L), and sodium acetate (0.075 mol/L) with a 7.0 pH value; methanol:acetonitrile (75:25, v/v) was used as the organic phase. The UV detection wavelength and flow rate were 355 nm and 0.8 mL/min, respectively, with the HPLC column temperature set at 27℃. The injection volume was 50 μL [6].

5.4. DNA and RNA Extraction, Library Construction, Sequencing

E.Z.N.A Soil Kit (Omega Bio-Tek, Norcross, GA) was used to extract gDNA from BSFL intestinal microbiota following the manufacturer’s instructions. The RNA of BSFL intestinal microbiota were extracted using our previously developed thiocyanate-high EDTA method [46]. A NanoDrop ND-2000 spectrophotometer (Thermo Fisher, MA, USA) and agarose gel electrophoresis were first used to determine the quantity and quality of gDNA/RNA. Six replicates of the CK and CR groups were mixed as one sample and then sent to Biozeron Ltd. (Biozeron, Shanghai, China) for the subsequent library construction and meta-sequencing via the Illumina Novaseq 6000 sequencing platform (Illumina, San Diego, CA, USA).

5.5. Data Sequencing Analysis

Raw data were extracted via the following steps to obtain a clean read dataset: (i) Trimmomatic (version0.35, available at https://github.com/usadellab/Trimmomatic) (accessed on 20 January 2022) was used to trim raw reads; (ii) reads with low-quality scores (<20) and short lengths (≤75 bp) were trimmed; (iii) Kraken 2 software annotated the species’ taxonomic data from clean reads; (iv) For RNA samples, the removal of rRNA contamination was achieved by comparing the data to the SILVA LSU (23 S/28 S) and SILVA SSU (16 S/18 S) databases via SortMeRNA software (version 2.1) (http://bioinfo.lifl.fr/RNA/sortmerna/, accessed on 20 January 2022).
Obtained clean reads were then assembled de novo using Megahit (http://github.com/voutcn/megahit, accessed on 22 January 2022). The scaffolds were then broken into contigs without gaps. Contigs of >500 bp in length were used for gene prediction. METAProdigal software (http://prodigal.ornl.gov/, accessed on 23 January 2022) was used to annotate the open reading frames (ORFs) of each sample, and CD-HIT software (http://www.bioinformatics.org/cd-hit/, accessed on 23 January 2022) was used to construct a non-redundant gene catalog with parameters set at 90% coverage and 95% consistency [47]. Next, all clean reads were further analyzed for taxonomic and functional annotations using various databases. The BLAST (BLAST Version 2.9.0+, http://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 26 January 2022) search of each unigene was conducted against the NCBI-NR, Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/, accessed on 26 January 2022), Non-supervised Orthologous Groups (eggNOG, http://www.ncbi.nlm.nih.gov/COG/, accessed on 26 January 2022), Comprehensive Antibiotic Resistance Database (CARD, https://card.mcmaster.ca/, accessed on 26 January 2022), Carbohydrate-Active EnZymes (CAZy, http://www.cazy.org/, accessed on 26 January 2022) and other databases to predict gene functions and metabolic pathways.

5.6. Data Availability

The datasets containing metagenomic and metatranscriptomic data were deposited into the NCBI Sequence Read Archive (SRA) using accession numbers PRJNA908247 and PRJNA908248, respectively.

5.7. Statistical Analysis

Statistical data were analyzed via GraphPad Prism software (8.3.0). One-way ANOVA analysis of variance was used to calculate the characterization differences in all groups. P < 0.05 was considered statistically significant.

6. Patents

Two patents were generated during the process of carrying out this work. The patent numbers are as follows: CN 02210699597.8 and CN202210699598.2.

Author Contributions

Y.P.: conceptualization, methodology, formal analysis, writing—original draft, writing—review and editing, project administration, funding acquisition, and supervision. M.S.: data curation and formal analysis. J.Z.: investigation and validation. A.L.: data curation and methodology. H.C.: project administration, funding acquisition, and supervision. X.K.: project administration, and supervision. H.N.: formal analysis, investigation, software, writing—original draft, writing—review and editing, and supervision. S.Y.: methodology, resources, visualization, writing—review and editing, supervision, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Key Research and Development Program of China (No. 2021YFD1700900, 2022YFC3401700), the Young Elite Scientists Sponsorship Program of the Henan Association for Science and Technology (No. 2022HYTP041), and the Key Technologies R&D Program of Henan Province (No. 232102320103).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article. The datasets of metagenomic and metatranscriptomic were deposited into the NCBI Sequence Read Archive (SRA), with accession numbers of PRJNA908247 and PRJNA908248, respectively.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kang, Y.; Li, Q.; Xia, D.; Shen, M.; Mei, L.; Hu, J. Short-term thermophilic treatment cannot remove tetracycline resistance genes in pig manures but exhibits controlling effects on their accumulation and spread in soil. J. Hazard. Mater. 2017, 340, 213–220. [Google Scholar] [CrossRef] [PubMed]
  2. Lyu, J.; Yang, L.; Zhang, L.; Ye, B.; Wang, L. Antibiotics in soil and water in China–a systematic review and source analysis. Environ. Pollut. 2020, 266, 115147. [Google Scholar] [CrossRef] [PubMed]
  3. Ben, Y.; Fu, C.; Hu, M.; Liu, L.; Wong, M.H.; Zheng, C. Human health risk assessment of antibiotic resistance associated with antibiotic residues in the environment: A review. Environ. Res. 2019, 169, 483–493. [Google Scholar] [CrossRef]
  4. Li, W.; Zhang, G. Detection and various environmental factors of antibiotic resistance gene horizontal transfer. Environ. Res. 2022, 212, 113267. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, C.; Yao, H.; Wang, C. Black soldier fly larvae can effectively degrade oxytetracycline bacterial residue by means of the gut bacterial community. Front. Microbiol. 2021, 12, 663972. [Google Scholar] [CrossRef] [PubMed]
  6. Shi, Y.; Lin, H.; Ma, J.; Zhu, R.; Sun, W.; Lin, X.; Zhang, J.; Zheng, H.; Zhang, X. Degradation of tetracycline antibiotics by Arthrobacter nicotianae OTC-16. J. Hazard. Mater. 2021, 403, 123996. [Google Scholar] [CrossRef] [PubMed]
  7. Kang, Y. Tetracycline degradation by Klebsiella sp. strain TR5: Proposed degradation pathway and possible genes involved. Chemosphere 2020, 253, 126729. [Google Scholar]
  8. Huang, X.; Zhang, X.; Feng, F.; Xu, X. Biodegradation of tetracycline by the yeast strain Trichosporon mycotoxinivorans XPY-10. Prep. Biochem. Biotechnol. 2016, 46, 15–22. [Google Scholar] [CrossRef]
  9. Leng, Y.; Bao, J.; Chang, G.; Zheng, H.; Li, X.; Du, J.; Snow, D.; Li, X. Biotransformation of tetracycline by a novel bacterial strain Stenotrophomonas maltophilia DT1. J. Hazard. Mater. 2016, 318, 125–133. [Google Scholar] [CrossRef]
  10. Bhatt, P.; Jeon, C.-H.; Kim, W. Tetracycline bioremediation using the novel Serratia marcescens strain WW1 isolated from a wastewater treatment plant. Chemosphere 2022, 298, 134344. [Google Scholar] [CrossRef]
  11. Chen, X.; Shen, W.; Chen, J.; Zhu, Y.; Chen, C.; Xie, S. Tetracycline biotransformation by a novel bacterial strain Alcaligenes sp. T17. Sci. Total Environ. 2022, 832, 155130. [Google Scholar] [CrossRef] [PubMed]
  12. Wen, X.; Jia, Y.; Li, J. Degradation of tetracycline and oxytetracycline by crude lignin peroxidase prepared from Phanerochaete chrysosporium—A white rot fungus. Chemosphere 2009, 75, 1003–1007. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, T.; Awasthi, M.K.; Chen, H.; Duan, Y.; Awasthi, S.K.; Zhang, Z. Performance of black soldier fly larvae (Diptera: Stratiomyidae) for manure composting and production of cleaner compost. J. Environ. Manag. 2019, 251, 109593. [Google Scholar] [CrossRef] [PubMed]
  14. Wu, N.; Wang, X.; Mao, Z.; Liang, J.; Liu, X.; Xu, X. Bioconversion of chicken meat and bone meal by black soldier fly larvae: Effects of straw addition on the quality and microbial profile of larval frass. J. Environ. Manag. 2022, 307, 114579. [Google Scholar] [CrossRef]
  15. Raksasat, R.; Lim, J.W.; Kiatkittipong, W.; Kiatkittipong, K.; Ho, Y.C.; Lam, M.K.; Font-Palma, C.; Zaid, H.F.Z.; Cheng, C.K. A review of organic waste enrichment for inducing palatability of black soldier fly larvae: Wastes to valuable resources. Environ. Pollut. 2020, 267, 115488. [Google Scholar] [CrossRef]
  16. Yang, C.; Ma, S.; Li, F.; Zheng, L.; Tomberlin, J.K.; Yu, Z.; Zhang, J.; Yu, C.; Fan, M.; Cai, M. Characteristics and mechanisms of ciprofloxacin degradation by black soldier fly larvae combined with associated intestinal microorganisms. Sci. Total Environ. 2021, 811, 151371. [Google Scholar] [CrossRef]
  17. Luo, X.; Yang, Q.; Lin, Y.; Tang, Z.; Tomberlin, J.K.; Liu, W.; Huang, Y. Black soldier fly larvae effectively degrade lincomycin from pharmaceutical industry wastes. J. Environ. Manag. 2022, 307, 114539. [Google Scholar] [CrossRef]
  18. Liu, C.; Yao, H.; Chapman, S.J.; Su, J.; Wang, C. Changes in gut bacterial communities and the incidence of antibiotic resistance genes during degradation of antibiotics by black soldier fly larvae. Environ. Int. 2020, 142, 105834. [Google Scholar] [CrossRef]
  19. Xiao, K.Q.; Li, L.G.; Ma, L.P.; Zhang, S.Y.; Bao, P.; Zhang, T.; Zhu, Y.G. Metagenomic analysis revealed highly diverse microbial arsenic metabolism genes in paddy soils with low-arsenic contents. Environ. Pollut. 2016, 211, 1–8. [Google Scholar] [CrossRef]
  20. Menezes, A.D.; Clipson, N.; Doyle, E. Comparative metatranscriptomics reveals widespread community responses during phenanthrene degradation in soil. Environ. Microbiol. 2012, 14, 2577–2588. [Google Scholar] [CrossRef]
  21. Zhao, X.-L.; Chen, Z.-G.; Yang, T.-C.; Jiang, M.; Wang, J.; Cheng, Z.-X.; Yang, M.-J.; Zhu, J.-X.; Zhang, T.-T.; Li, H. Glutamine promotes antibiotic uptake to kill multidrug-resistant uropathogenic bacteria. Sci. Transl. Med. 2021, 13, eabj0716. [Google Scholar] [CrossRef] [PubMed]
  22. Cai, M.; Ma, S.; Hu, R.; Tomberlin, J.K.; Yu, C.; Huang, Y.; Zhan, S.; Li, W.; Zheng, L.; Yu, Z.; et al. Systematic characterization and proposed pathway of tetracycline degradation in solid waste treatment by Hermetia illucens with intestinal microbiota. Environ. Pollut. 2018, 242, 634–642. [Google Scholar] [CrossRef]
  23. Mei, H.; Li, C.; Li, X.; Hu, B.; Lu, L.; Tomberlin, J.K.; Hu, W. Characteristics of tylosin and enrofloxacin degradation in swine manure digested by black soldier fly (Hermetia illucens L.) larvae. Environ. Pollut. 2022, 293, 118495. [Google Scholar] [CrossRef]
  24. Yun, J.H.; Roh, S.W.; Whon, T.W.; Jung, M.J.; Bae, J.W. Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl. Environ. Microbiol. 2014, 80, 5254–5264. [Google Scholar] [CrossRef]
  25. Peng, D.; Feng, Y.Z.; Wang, Y.M.; Lin, X.G. Enrichment of the antibiotic resistance gene tet(L) in an alkaline soil fertilized with plant derived organic manure. Front. Microbiol. 2018, 9, 1140. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, Z.J.; Shen, J.G.; Wang, H.; Liu, M.; Wu, L.H.; Ping, F.; He, Q.; Li, H.Y.; Zheng, C.F.; Xu, X.H. Attenuation of veterinary antibiotics in full-scale vermicomposting of swine manure via the housefly larvae (Musca domestica). Sci. Rep. 2013, 4, 6844. [Google Scholar] [CrossRef] [PubMed]
  27. Jung, J.Y.; Ahn, Y.; Khare, S.; Gokulan, K.; Piñeiro, S.A.; Cerniglia, C.E. An in vitro study to assess the impact of tetracycline on the human intestinal microbiome. Anaerobe 2018, 49, 85–94. [Google Scholar] [CrossRef]
  28. Cauwerts, K.; Decostere, A.; Graef, E.M.D.; Haesebrouck, F.; Pasmans, F. High prevalence of tetracycline resistance in Enterococcus isolates from broilers carrying the erm(B) gene. Avian Pathol. 2007, 36, 395–399. [Google Scholar] [CrossRef]
  29. Niu, S.-H.; Liu, S.; Deng, W.-K.; Wu, R.-T.; Cai, Y.-F.; Liao, X.-D.; Xing, S.-C. A sustainable and economic strategy to reduce risk antibiotic resistance genes during poultry manure bioconversion by black soldier fly Hermetia illucens larvae: Larval density adjustment. Ecotoxicol. Environ. Saf. 2022, 232, 113294. [Google Scholar] [CrossRef]
  30. Li, K.; Cao, R.; Mo, S.; Yao, R.; Wu, J. Swine manure composting with compound microbial inoculants: Removal of antibiotic resistance genes and their associations with microbial community. Front. Microbiol. 2020, 11, 592592. [Google Scholar] [CrossRef]
  31. Khaira, M.B.; Yusuf, M.B.; Khan, F. Insights to antimicrobial resistance: Heavy metals can inhibit antibiotic resistance in bacteria isolated from wastewater. Environ. Monit. Assess. 2022, 194, 252. [Google Scholar] [CrossRef] [PubMed]
  32. Webber, M.A.; Piddock, L.J.V. The importance of efflux pumps in bacterial antibiotic resistance. J. Antimicrob. Chemother. 2003, 51, 9–11. [Google Scholar] [CrossRef] [PubMed]
  33. Molale, L.G.; Bezuidenhout, C.C. Antibiotic resistance, efflux pump genes and virulence determinants in Enterococcus spp. from surface water systems. Environ. Sci. Pollut. Res. 2016, 23, 21501–21510. [Google Scholar] [CrossRef] [PubMed]
  34. Jordi Vila, J.L.M. Clinical impact of the over-expression of efflux pump in nonfermentative Gram-negative bacilli, development of efflux pump inhibitors. Curr. Drug Targets 2008, 9, 797–807. [Google Scholar] [CrossRef] [PubMed]
  35. Ss, A.; Assg, A.; Jvp, B. Evaluation of the inhibitory effect of caffeic acid and gallic acid on tetR and tetM efflux pumps mediating tetracycline resistance in Streptococcus sp., using computational approach. J. King Saud Univ. Sci. 2020, 32, 904–909. [Google Scholar]
  36. Scott, L.J. Eravacycline: A review in complicated intra-abdominal infections. Drugs 2019, 79, 315–324. [Google Scholar] [CrossRef]
  37. Tian, B.; Fadhil, N.H.; Powell, J.E.; Kwong, W.K.; Moran, N.A. Long-term exposure to antibiotics has caused accumulation of resistance determinants in the gut microbiota of honeybees. MBio 2012, 3, e00312–e00377. [Google Scholar] [CrossRef]
  38. Liu, B.; Chen, C. Translation elongation factor 4 (LepA) contributes to tetracycline susceptibility by stalling elongating ribosomes. Antimicrob. Agents Chemother. 2018, 62, 10-1128. [Google Scholar] [CrossRef]
  39. Yang, J.; Lin, Y.; Yang, X.; Ng, T.B.; Ye, X.; Lin, J. Degradation of tetracycline by immobilized laccase and the proposed transformation pathway. J. Hazard. Mater. 2017, 322, 525–531. [Google Scholar] [CrossRef]
  40. Wen, X.; Zeng, Z.; Du, C.; Huang, D.; Zeng, G.; Xiao, R.; Lai, C.; Xu, P.; Zhang, C.; Wan, J. Immobilized laccase on bentonite-derived mesoporous materials for removal of tetracycline. Chemosphere 2019, 222, 865–871. [Google Scholar] [CrossRef]
  41. Wen, X.; Jia, Y.; Li, J. Enzymatic degradation of tetracycline and oxytetracycline by crude manganese peroxidase prepared from Phanerochaete chrysosporium. J. Hazard. Mater. 2010, 177, 924–928. [Google Scholar] [CrossRef] [PubMed]
  42. Sun, X.; Leng, Y.; Wan, D.; Chang, F.; Huang, Y.; Li, Z.; Xiong, W.; Wang, J. Transformation of tetracycline by manganese peroxidase from Phanerochaete chrysosporium. Molecules 2021, 26, 6803. [Google Scholar] [CrossRef] [PubMed]
  43. Jeon, J.H.; Kim, S.J.; Lee, H.S.; Cha, S.S.; Lee, J.H.; Yoon, S.H.; Koo, B.S.; Lee, C.M.; Choi, S.H.; Lee, S.H. Novel metagenome-derived carboxylesterase that hydrolyzes β-lactam antibiotics. Appl. Environ. Microbiol. 2011, 77, 7830–7836. [Google Scholar] [CrossRef] [PubMed]
  44. Zhou, G.; He, Y.; Zhang, H. Screening and 16S rRNA Analysis of the Bacteria of Producing Chitin Deacetylase. In Proceedings of the 4th International Conference on Bioinformatics and Biomedical Engineering, Chengdu, China, 18–20 June 2010; pp. 1–4. [Google Scholar] [CrossRef]
  45. Pei, Y.; Zhao, S.; Chen, X.; Zhang, J.; Ni, H.; Sun, M.; Lin, H.; Liu, X.; Chen, H.; Yang, S. Bacillus velezensis EEAM 10B strengthens nutrient metabolic process in black soldier fly larvae (Hermetia illucens) via changing gut microbiome and metabolic pathways. Front. Nutr. 2022, 9, 880488. [Google Scholar] [CrossRef]
  46. Pei, Y.; Mamtimin, T.; Ji, J.; Khan, A.; Kakade, A.; Zhou, T.; Yu, Z.; Zain, H.; Yang, W.; Ling, Z.; et al. The guanidine thiocyanate-high EDTA method for total microbial RNA extraction from severely heavy metal-contaminated soils. Microb. Biotechnol. 2021, 14, 465–478. [Google Scholar] [CrossRef]
  47. Hyatt, D.; LoCascio, P.F.; Hauser, L.J.; Uberbacher, E.C. Gene and translation initiation site prediction in metagenomic sequences. Bioinformatics 2012, 28, 2223–2230. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The bio-physiological parameters’ (a) survival rate, (b) Biomass (g/10 larva), (c) substate consumption rate of BSFL in the CK and TC groups, and (d) tetracycline degradation efficiency in the CK, germ-free, and natural BSFL groups. Data are presented as mean ± standard deviation. Values with different letters mean significant differences at p < 0.05, as determined via Tukey’s test.
Figure 1. The bio-physiological parameters’ (a) survival rate, (b) Biomass (g/10 larva), (c) substate consumption rate of BSFL in the CK and TC groups, and (d) tetracycline degradation efficiency in the CK, germ-free, and natural BSFL groups. Data are presented as mean ± standard deviation. Values with different letters mean significant differences at p < 0.05, as determined via Tukey’s test.
Toxics 11 00611 g001
Figure 2. BSFL intestinal microbial composition at (a) phylum and (b) genus levels based on metagenomic sequencing data; (c) phylum and (d) genus levels based on metatranscriptomic sequencing data in CK and TC groups.
Figure 2. BSFL intestinal microbial composition at (a) phylum and (b) genus levels based on metagenomic sequencing data; (c) phylum and (d) genus levels based on metatranscriptomic sequencing data in CK and TC groups.
Toxics 11 00611 g002
Figure 3. Abundances of antibiotic resistance gene (AMR) types in the (a) CK and (b) TC groups. Abundances of Tet family specific genes in the (c) CK and (d) TC groups based on metagenomic sequencing data.
Figure 3. Abundances of antibiotic resistance gene (AMR) types in the (a) CK and (b) TC groups. Abundances of Tet family specific genes in the (c) CK and (d) TC groups based on metagenomic sequencing data.
Toxics 11 00611 g003
Figure 4. Abundances of antibiotic resistance gene (AMR) types in the (a) CK and (b) TC groups. Abundances of Tet family specific genes in the (c) CK and (d) TC groups based on metatranscriptomic sequencing data.
Figure 4. Abundances of antibiotic resistance gene (AMR) types in the (a) CK and (b) TC groups. Abundances of Tet family specific genes in the (c) CK and (d) TC groups based on metatranscriptomic sequencing data.
Toxics 11 00611 g004
Figure 5. Differential gene expression analysis between the TC and CK groups. (a) Differences in gene expression profiles; (b) the top 15 KEGG enrichment terms of differentially expressed genes.
Figure 5. Differential gene expression analysis between the TC and CK groups. (a) Differences in gene expression profiles; (b) the top 15 KEGG enrichment terms of differentially expressed genes.
Toxics 11 00611 g005
Figure 6. Distribution of COG functional categories in the CK and TC groups based on (a) metagenomic and (b) metatranscriptomic sequencing data.
Figure 6. Distribution of COG functional categories in the CK and TC groups based on (a) metagenomic and (b) metatranscriptomic sequencing data.
Toxics 11 00611 g006
Figure 7. The relative abundances of genes that encode different CAZy classes (GTs: Glycosyl Transferases; GHs: Glycoside Hydrolases; CEs: Carbohydrate Esterases; AAs: Auxiliary Activities; CBMs: Carbohydrate-Binding Modules; PLs: Polysaccharide) in the CAZy database at level 1 in the CK and TC groups based on (a) metagenomic and (b) metatranscriptomic sequencing data.
Figure 7. The relative abundances of genes that encode different CAZy classes (GTs: Glycosyl Transferases; GHs: Glycoside Hydrolases; CEs: Carbohydrate Esterases; AAs: Auxiliary Activities; CBMs: Carbohydrate-Binding Modules; PLs: Polysaccharide) in the CAZy database at level 1 in the CK and TC groups based on (a) metagenomic and (b) metatranscriptomic sequencing data.
Toxics 11 00611 g007
Figure 8. Speculated tetracycline response mechanisms of BSFL’s intestinal microbiota. (i) Efflux pump, () Dissociate TC from ribosome, () Deacetylase, and () Novel tetracycline degradation enzymes.
Figure 8. Speculated tetracycline response mechanisms of BSFL’s intestinal microbiota. (i) Efflux pump, () Dissociate TC from ribosome, () Deacetylase, and () Novel tetracycline degradation enzymes.
Toxics 11 00611 g008
Table 1. Metagenomics and metatranscriptomics sequencing data.
Table 1. Metagenomics and metatranscriptomics sequencing data.
SamplesCK-MGTC-MGCK-MTTC-MT
Raw Reads113499814121370872115167960117241914
Raw Bases17024972100182056308001727519400017586287100
Clean Reads8915436810409831210684055294504924
Clean Bases13264807930154909764191592737462614103693095
Contigs41,58643,05795262500
Clean Reads in Raw Reads (%)78.55%85.77%92.77%80.61%
Clean Bases in Raw Bases (%)77.91%85.09%92.20%80.20%
Q2097.63%97.48%98.25%98.73%
Q3095.23%94.89%94.06%95.60%
Table 2. Abundance values of partial carbohydrate-active enzymes (AA&CE) in the CK and TC groups.
Table 2. Abundance values of partial carbohydrate-active enzymes (AA&CE) in the CK and TC groups.
FamilyCK-MGTC-MGCK-MTTC-MTDescription
AA10.502221.76443--Laccase/p-diphenol:oxygen oxidoreductase/ferroxidase (EC 1.10.3.2); ferroxidase (EC 1.10.3.-); Laccase-like multicopper oxidase (EC 1.10.3.-)
AA22.227359.84698--Manganese peroxidase (EC 1.11.1.13); versatile peroxidase (EC 1.11.1.16); lignin peroxidase (EC 1.11.1.14); peroxidase (EC 1.11.1.-)
AA3334.669310.651177.314153.316Cellobiose dehydrogenase (EC 1.1.99.18); glucose 1-oxidase (EC 1.1.3.4); aryl alcohol oxidase (EC 1.1.3.7); alcohol oxidase (EC 1.1.3.13); pyranose oxidase (EC 1.1.3.10)
AA4189.943294.33--Vanillyl-alcohol oxidase (EC 1.1.3.38)
AA6199.316380.5124541.034518.021,4-benzoquinone reductase (EC. 1.6.5.6)
AA765.8805125.302--Glucooligosaccharide oxidase (EC 1.1.3.-); chitooligosaccharide oxidase (EC 1.1.3.-)
AA1075.5844146.551--AA10 (formerly CBM33) proteins are copper-dependent lytic polysaccharide mono-oxygenases (LPMOs); some proteins have been shown to act on chitin, while others act on cellulose.
AA122.682675.62463--The pyrroloquinoline quinone-dependent oxidoreductase activity was demonstrated as present for the CC1G_09525 protein of Coprinopsis cinerea.
CE11464.21816.33864.206910.166Acetyl xylan esterase (EC 3.1.1.72); cinnamoyl esterase (EC 3.1.1.-); feruloyl esterase (EC 3.1.1.73); carboxylesterase (EC 3.1.1.1); S-formylglutathione hydrolase (EC 3.1.2.12)
CE237.894530.416--Acetyl xylan esterase (EC 3.1.1.72).
CE3180.225156.92325.027629.0947Acetyl xylan esterase (EC 3.1.1.72).
CE4501.018768.3341443.961964.17Acetyl xylan esterase (EC 3.1.1.72); chitin deacetylase (EC 3.5.1.41); chito-oligosaccharide deacetylase (EC 3.5.1.-); peptidoglycan GlcNAc deacetylase (EC 3.5.1.-); peptidoglycan N-acetylmuramic acid deacetylase (EC 3.5.1.-).
CE50.965993.93697--Acetyl xylan esterase (EC 3.1.1.72); cutinase (EC 3.1.1.74)
CE639.713777.6012133.101948.039Acetyl xylan esterase (EC 3.1.1.72).
CE761.6056177.51169.404781.6384Acetyl xylan esterase (EC 3.1.1.72); cephalosporin-C deacetylase (EC 3.1.1.41).
CE843.02780.070918.128725.6374Pectin methylesterase (EC 3.1.1.11).
CE9294.834212.118474.72395.33N-acetylglucosamine 6-phosphate deacetylase (EC 3.5.1.25); N-acetylgalactosamine-6-phosphate deacetylase (EC 3.5.1.80).
CE10596.7321009.64225.639135.508Arylesterase (EC 3.1.1.-); carboxyl esterase (EC 3.1.1.3); acetylcholinesterase (EC 3.1.1.7); cholinesterase (EC 3.1.1.8); sterol esterase (EC 3.1.1.13); brefeldin A esterase (EC 3.1.1.-).
CE11376.726322.36228.196301.047UDP-3-0-acyl N-acetylglucosamine deacetylase (EC 3.5.1.-).
CE1298.2928251.826203.141122.558Pectin acetylesterase (EC 3.1.1.-); rhamnogalacturonan acetylesterase (EC 3.1.1.-); acetyl xylan esterase (EC 3.1.1.72)
CE1446.9864308.344--N-acetyl-1-D-myo-inosityl-2-amino-2-deoxy-alpha-D-glucopyranoside deacetylase (EC 3.5.1.89); diacetylchitobiose deacetylase (EC 3.5.1.-); mycothiol S-conjugate amidase (EC 3.5.1.-)
CE155.1411325.4286--4-O-methyl-glucuronoyl methylesterase (EC 3.1.1.-)
CE1622.568551.4155--Acetylesterase (EC 3.1.1.6) is active on various carbohydrate acetyl esters
Others20,856.926,276.421,728.525,421-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pei, Y.; Sun, M.; Zhang, J.; Lei, A.; Chen, H.; Kang, X.; Ni, H.; Yang, S. Comparative Metagenomic and Metatranscriptomic Analyses Reveal the Response of Black Soldier Fly (Hermetia illucens) Larvae Intestinal Microbes and Reduction Mechanisms to High Concentrations of Tetracycline. Toxics 2023, 11, 611. https://doi.org/10.3390/toxics11070611

AMA Style

Pei Y, Sun M, Zhang J, Lei A, Chen H, Kang X, Ni H, Yang S. Comparative Metagenomic and Metatranscriptomic Analyses Reveal the Response of Black Soldier Fly (Hermetia illucens) Larvae Intestinal Microbes and Reduction Mechanisms to High Concentrations of Tetracycline. Toxics. 2023; 11(7):611. https://doi.org/10.3390/toxics11070611

Chicago/Turabian Style

Pei, Yaxin, Mengxiao Sun, Jiran Zhang, Aojie Lei, Hongge Chen, Xiangtao Kang, Hongyuhang Ni, and Sen Yang. 2023. "Comparative Metagenomic and Metatranscriptomic Analyses Reveal the Response of Black Soldier Fly (Hermetia illucens) Larvae Intestinal Microbes and Reduction Mechanisms to High Concentrations of Tetracycline" Toxics 11, no. 7: 611. https://doi.org/10.3390/toxics11070611

APA Style

Pei, Y., Sun, M., Zhang, J., Lei, A., Chen, H., Kang, X., Ni, H., & Yang, S. (2023). Comparative Metagenomic and Metatranscriptomic Analyses Reveal the Response of Black Soldier Fly (Hermetia illucens) Larvae Intestinal Microbes and Reduction Mechanisms to High Concentrations of Tetracycline. Toxics, 11(7), 611. https://doi.org/10.3390/toxics11070611

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop