Transcriptomic and metabolomic landscape of the molecular effects of glyphosate commercial formulation on Apis mellifera ligustica and Apis cerana cerana

Synergistic effects between microplastics and glyphosate on honey bee larvae

Microplastic (MPs) pollution has emerged as a global ecological concern, however, the impact of MPs exposure, particularly in conjunction with other pollutants such as glyphosate (GLY) on honey bee remains unknown. This study investigated the effects of exposure to different concentrations of MPs and their combination with GLY on honey bee larvae development, or during the larvae period, regulation of major detoxification, antioxidant and immune genes, and oxidative stress biomarkers. Results revealed that combined exposure to MPs and GLY decreased larvae survivorship and weight, while exposure to MPs alone showed no significant differences. Both MPs and GLY alone downregulated the defensin-1 gene, but only combined exposure with GLY downregulated the hymenoptaecin gene and increased catalase enzyme activity. The data suggest a synergistic effect of MPs and GLY, leading to immunosuppression and reduced larvae survival and weight. These findings highlight potential risks of two prevalent environmental pollutants on honey bee health.

Intake of imidacloprid in lethal and sublethal doses alters gene expression in Apis mellifera bees

Bees are important pollinators for ecosystems and agriculture; however, populations have suffered a decline that may be associated with several factors, including habitat loss, climate change, increased vulnerability to diseases and parasites and use of pesticides. The extensive use of neonicotinoids, including imidacloprid, as agricultural pesticides, leads to their persistence in the environment and accumulation in bees, pollen, nectar, and honey, thereby inducing deleterious effects. Forager honey bees face significant exposure to pesticide residues while searching for resources outside the hive, particularly systemic pesticides like imidacloprid. In this study, 360 Apis mellifera bees, twenty-one days old (supposed to be in the forager phase) previously marked were fed syrup (honey and water, 1:1 m/v) containing a lethal dose (0.081 μg/bee) or sublethal dose (0.00081 μg/bee) of imidacloprid. The syrup was provided in plastic troughs, with 250 μL added per trough onto each plastic Petri dish containing 5 bees (50 μL per bee). The bees were kept in the plastic Petri dishes inside an incubator, and after 1 and 4 h of ingestion, the bees were euthanised and stored in an ultra-freezer (-80 °C) for transcriptome analysis. Following the 1-h ingestion of imidacloprid, 1516 genes (73 from lethal dose; 1509 from sublethal dose) showed differential expression compared to the control, while after 4 h, 758 genes (733 from lethal dose; 25 from sublethal) exhibited differential expression compared to the control. All differentially expressed genes found in the brain tissue transcripts of forager bees were categorised based on gene ontology into functional groups encompassing biological processes, molecular functions, and cellular components. These analyses revealed that sublethal doses might be capable of altering more genes than lethal doses, potentially associated with a phenomenon known as insecticide-induced hormesis. Alterations in genes related to areas such as the immune system, nutritional metabolism, detoxification system, circadian rhythm, odour detection, foraging activity, and memory in bees were present after exposure to the pesticide. These findings underscore the detrimental effects of both lethal and sublethal doses of imidacloprid, thereby providing valuable insights for establishing public policies regarding the use of neonicotinoids, which are directly implicated in the compromised health of Apis mellifera bees.

Changes in the proteome of Apis mellifera acutely exposed to sublethal dosage of glyphosate and imidacloprid

Uncontrolled use of pesticides has caused a dramatic reduction in the number of pollinators, including bees. Studies on the effects of pesticides on bees have reported effects on both metabolic and neurological levels under chronic exposure. In this study, variations in the differential expression of head and thorax-abdomen proteins in Africanized A. mellifera bees treated acutely with sublethal doses of glyphosate and imidacloprid were studied using a proteomic approach. A total of 92 proteins were detected, 49 of which were differentially expressed compared to those in the control group (47 downregulated and 2 upregulated). Protein interaction networks with differential protein expression ratios suggested that acute exposure of A. mellifera to sublethal doses of glyphosate could cause head damage, which is mainly associated with behavior and metabolism. Simultaneously, imidacloprid can cause damage associated with metabolism as well as, neuronal damage, cellular stress, and impairment of the detoxification system. Regarding the thorax-abdomen fractions, glyphosate could lead to cytoskeleton reorganization and a reduction in defense mechanisms, whereas imidacloprid could affect the coordination and impairment of the oxidative stress response.

Impact of chronic exposure to field level glyphosate on the food consumption, survival, gene expression, gut microbiota, and metabolomic profiles of honeybees

Glyphosate (GLY) is among the most widely used pesticides in the world. However, there are a lot of unknowns about chronic exposure to GLY's effects on Honeybee (HB) behavior and physiology. To address this, we carried out five experiments to study the impact of chronic exposure to 5 mg/kg GLY on sugar consumption, survival, gene expression, gut microbiota, and metabolites of HB workers. Our results find a significant decrease in sugar consumption and survival probability of HB after chronic exposure to GLY. Further, genes associated with immune response, energy metabolism, and longevity were conspicuously altered. In addition, a total of seven metabolites were found to be differentially expressed in the metabolomic profiles, mainly related the sucrose metabolism. There was no significant difference in the gut microbiota. Results suggest that chronic exposure to field-level GLY altered the health of HB and the intricate toxic mechanisms. Our data provided insights into the chronic effects of GLY on HB behavior in food intake and health, which represents the field conditions where HB are exposed to pesticides over extended periods.

Current Insights into Sublethal Effects of Pesticides on Insects

The effect of pesticides on insects is often discussed in terms of acute and chronic toxicity, but an important and often overlooked aspect is the impact of sublethal doses on insect physiology and behavior. Pesticides can influence various physiological parameters of insects, including the innate immune system, development, and reproduction, through a combination of direct effects on specific exposed tissues and the modification of behaviors that contribute to health and reproductive success. Such behaviors include mobility, feeding, oviposition, navigation, and the ability to detect pheromones. Pesticides also have a profound effect on insect learning and memory. The precise effects depend on many different factors, including the insect species, age, sex, caste, physiological condition, as well as the type and concentration of the active ingredients and the exposure route. More studies are needed to assess the effects of different active ingredients (and combinations thereof) on a wider range of species to understand how sublethal doses of pesticides can contribute to insect decline. This review reflects our current knowledge about sublethal effects of pesticides on insects and advancements in the development of innovative methods to detect them.

A systemic study of cyenopyrafen in strawberry cultivation system: Efficacy, residue behavior, and impact on honeybees (Apis mellifera L.)

The pesticide application method is one of the important factors affecting its effectiveness and residues, and the risk of pesticides to non-target organisms. To elucidate the effect of application methods on the efficacy and residue of cyenopyrafen, and the toxic effects on pollinators honeybees in strawberry cultivation, the efficacy and residual behavior of cyenopyrafen were investigated using foliar spray and backward leaf spray in field trials. The results showed that the initial deposition of cyenopyrafen using backward leaf spray on target leaves reached 5.06-9.81 mg/kg at the dose of 67.5-101.25 g a.i./ha, which was higher than that using foliar spray (2.62-3.71 mg/kg). The half-lives of cyenopyrafen in leaves for foliar and backward leaf spray was 2.3-3.3 and 5.3-5.9 d, respectively. The residues (10 d) of cyenopyrafen in leaves after backward leaf spray was 1.41-3.02 mg/kg, which was higher than that after foliar spraying (0.25-0.37 mg/kg). It is the main reason for the better efficacy after backward leaf spray. However, the residues (10 d) in strawberry after backward leaf spray and foliar spray was 0.04-0.10 and < 0.01 mg/kg, which were well below the established maximum residue levels of cyenopyrafen in Japan and South Korea for food safety. To further investigate the effects of cyenopyrafen residues after backward leaf spray application on pollinator honeybees, sublethal effects of cyenopyrafen on honeybees were studied. The results indicated a significant inhibition in the detoxification metabolic enzymes of honeybees under continuous exposure of cyenopyrafen (0.54 and 5.4 mg/L) over 8 d. The cyenopyrafen exposure also alters the composition of honeybee gut microbiota, such as increasing the relative abundance of Rhizobiales and decreasing the relative abundance of Acetobacterales. The comprehensive data on cyenopyrafen provide basic theoretical for environmental and ecological risk assessment, while backward leaf spray proved to be effective and safe for strawberry cultivation.

Odorant Receptors Expressing and Antennal Lobes Architecture Are Linked to Caste Dimorphism in Asian Honeybee, Apis cerana (Hymenoptera: Apidae)

Insects heavily rely on the olfactory system for food, mating, and predator evasion. However, the caste-related olfactory differences in Apis cerana, a eusocial insect, remain unclear. To explore the peripheral and primary center of the olfactory system link to the caste dimorphism in A. cerana, transcriptome and immunohistochemistry studies on the odorant receptors (ORs) and architecture of antennal lobes (ALs) were performed on different castes. Through transcriptomesis, we found more olfactory receptor genes in queens and workers than in drones, which were further validated by RT-qPCR, indicating caste dimorphism. Meanwhile, ALs structure, including volume, surface area, and the number of glomeruli, demonstrated a close association with caste dimorphism. Particularly, drones had more macroglomeruli possibly for pheromone recognition. Interestingly, we found that the number of ORs and glomeruli ratio was nearly 1:1. Also, the ORs expression distribution pattern was very similar to the distribution of glomeruli volume. Our results suggest the existence of concurrent plasticity in both the peripheral olfactory system and ALs among different castes of A. cerana, highlighting the role of the olfactory system in labor division in insects.

Glyphosate affects larval gut microbiota and metamorphosis of honey bees with differences between rearing procedures

The honey bee Apis mellifera is a sentinel species of the pollinator community which is exposed to a wide variety of pesticides. In the last half-century, the pesticide most applied worldwide has been the herbicide glyphosate (GLY) used for weed control and with microbiocide effects. After its application in crops, the GLY residues have been detected in flowers visited by honey bees as well as in the stored food of their hives. Therefore, the honey bee brood can ingest the herbicide during larval development. Recent studies proved that GLY has detrimental effects on adult honey bees and other insects associated with the disturbance of their gut microbiota. GLY induces changes in the growth, metabolism and survival of honey bees and stingless bees reared in vitro. However, the effect of GLY on larval microbiota is unknown so far and there are few studies with an in-hive exposure to GLY. For these reasons, this study aims to determine whether GLY induces dysbiosis in honey bee larvae and affects their metamorphosis during the exposure period (pre-defecation) and the post-exposure period. Furthermore, we assessed this herbicide in vitro and in the hive to compare its effects on different rearing procedures. Finally, we tested the pigment BLUE1 as an indirect exposure marker to detect and estimate the in-hive intake concentration of GLY. Our results indicate that the intake of field-relevant concentrations of GLY induced a slowdown in growth with dysbiosis in the larval gut microbiota followed by late effects on their metamorphosis such as teratogenesis and mortality of newly emerged bees. Nevertheless, brood from the same colonies expressed different signs of toxicity depending on the rearing procedure and in a dose-dependent manner.

Sulfoxaflor effects depend on the interaction with other pesticides and Nosema ceranae infection in the honey bee (Apis mellifera)

Honey bees health is compromised by many factors such as the use of agrochemicals in agriculture and the various diseases that can affect them. Multiple studies have shown that these factors can interact, producing a synergistic effect that can compromise the viability of honey bees. This study analyses the interactions between different pesticides and the microsporidium Nosema ceranae and their effect on immune and detoxification gene expression, sugar consumption and mortality in the Iberian western honey bee (Apis mellifera iberiensis). For this purpose, workers were infected with N. ceranae and subjected to a sugar-water diet with field concentrations of the pesticides sulfoxaflor, azoxystrobin and glyphosate. Increased sugar intake and altered immune and cytochrome P450 gene expression were observed in workers exposed to sulfoxaflor and infected with N. ceranae. None of the pesticides affected Nosema spore production in honey bee gut. Of the three pesticides tested (alone or in combination) only sulfoxaflor increased mortality in honey bees. Taken together, our results suggest that the effects of sulfoxaflor were attenuated in contact with other pesticides, and that Nosema infection leads to increase sugar intake in sulfoxaflor-exposed bees. Overall, this underlines the importance of studying the interaction between different stressors to understand their overall impact not only on honey bee but also on wild bees health.

Gene expression in bumble bee larvae differs qualitatively between high and low concentration imidacloprid exposure levels

Neonicotinoid pesticides negatively impact bumble bee health, even at sublethal concentrations. Responses to the neonicotinoid imidacloprid have been studied largely at individual adult and colony levels, focusing mostly on behavioral and physiological effects. Data from developing larvae, whose health is critical for colony success, are deficient, particularly at the molecular level where transcriptomes can reveal disruption of fundamental biological pathways. We investigated gene expression of Bombus impatiens larvae exposed through food provisions to two field-realistic imidacloprid concentrations (0.7 and 7.0 ppb). We hypothesized both concentrations would alter gene expression, but the higher concentration would have greater qualitative and quantitative effects. We found 678 genes differentially expressed under both imidacloprid exposures relative to controls, including mitochondrial activity, development, and DNA replication genes. However, more genes were differentially expressed with higher imidacloprid exposure; uniquely differentially expressed genes included starvation response and cuticle genes. The former may partially result from reduced pollen use, monitored to verify food provision use and provide additional context to results. A smaller differentially expressed set only in lower concentration larvae, included neural development and cell growth genes. Our findings show varying molecular consequences under different field-realistic neonicotinoid concentrations, and that even low concentrations may affect fundamental biological processes.

The effects of glyphosate, pure or in herbicide formulation, on bumble bees and their gut microbial communities

The widespread use of glyphosate-based formulations to eliminate unwanted vegetation has increased concerns regarding their effects on non-target organisms, such as honey bees and their gut microbial communities. These effects have been associated with both glyphosate and co-formulants, but it is still unknown whether they translate to other bee species. In this study, we tested whether glyphosate, pure or in herbicide formulation, can affect the gut microbiota and survival rates of the eastern bumble bee, Bombus impatiens. We performed mark-recapture experiments with bumble bee workers from four different commercial colonies, which were exposed to field relevant concentrations of glyphosate or a glyphosate-based formulation (0.01 mM to 1 mM). After a 5-day period of exposure, we returned the bees to their original colonies, and they were sampled at days 0, 3 and 7 post-exposure to investigate changes in microbial community and microbiota resilience by 16S rRNA amplicon sequencing and quantitative PCR. We found that exposure to glyphosate, pure or in herbicide formulation, reduced the relative abundance of a beneficial bee gut bacterium, Snodgrassella, in bees from two of four colonies when compared to control bees at day 0 post-exposure, but this reduction became non-significant at days 3 and 7 post-exposure, suggesting microbiota resilience. We did not find significant changes in total bacteria between control and exposed bees. Moreover, we observed an overall trend in decreased survival rates in bumble bees exposed to 1 mM herbicide formulation during the 7-day post-exposure period, suggesting a potential negative effect of this formulation on bumble bees.

Risk assessment of predatory lady beetle Propylea japonica's multi-generational exposure to three non-insecticidal agrochemicals

The effects of non-insecticidal agrochemicals on pest natural predators remain largely unexplored except bees and silkworm. The herbicide quizalofop-p-ethyl (QpE), fungicide thiophanate-methyl (TM), and plant growth regulator mepiquat chloride (MC) have been extensively applied as non-insecticidal agrochemicals. Here, we systematically evaluated multiple effects of these 3 non-insecticidal agrochemicals on three generations of Propylea japonica, an important agroforestry predatory beetle, including the effects on its development, reproduction, enterobacteria, and transcriptomic response. The results showed that QpE exhibited a hormetic effect on P. japonica, thus significantly increasing the survival rate of generation 2 (F2) females, generation 3 (F3) females, and F3 males and body weight of F3 males. However, three successive generations exposed to TM and MC had no significant effect on longevity, body weight, survival rate, pre-oviposition period, and fecundity of P. japonica. Additionally, we investigated the effects of MC, TM, and QpE exposure on gene expression and gut bacterial community of F3 P. japonica. Under MC, TM, and QpE exposure, the overwhelming genes of P. japonica (99.90 %, 99.45 %, and 99.7 %) remained unaffected, respectively. Under TM and MC exposure, differentially expressed genes (DEGs) were not significantly enriched in any KEGG pathway, indicating TM and MC did not significantly affect functions of P. japonica, but under QpE exposure, the expression levels of drug metabolism-related genes were down-regulated. Although QpE treatment did not affect gut dominant bacterial community composition, it significantly increased relative abundances of detoxification metabolism-related bacteria such as Wolbachia, Pseudomonas and Burkholderia in P. japonica. However, TM and MC had no significant effect on the gut bacterial community composition and relative abundance in P. japonica. This study revealed for the first time the mechanism by which P. japonica might compensate for gene downregulation-induced detoxification metabolism decline through altering symbiotic bacteria under QpE exposure. Our findings provide reference for the rational application of non-insecticidal agrochemicals.

Fox transcription factor AccGRF1 in response to glyphosate stress in Apis cerana cerana

Glyphosate is an herbicide commonly used in agriculture, and its widespread use has adversely affected the survival of nontarget organisms. Among these organisms, bees in particular are important pollinators, and declining bee populations have severely affected crop yields around the world. However, the molecular mechanism by which glyphosate harms bees remains unclear. In our experiment, we screened and cloned a glyphosate-induced gene in Apis cerana cerana (A. c. cerana) and named glyphosate response factor 1 (AccGRF1). Sequence analysis showed that AccGRF1 contains a winged-helix DNA binding domain, which suggests that it belongs to the Forkhead box (Fox) protein family. qRT-PCR and heterologous expression in Escherichia coli and yeast showed that AccGRF1 can respond to glyphosate and oxidative stress. After AccGRF1 knockdown by means of RNA interference (RNAi), the resistance of A. c. cerana to glyphosate stress improved. The results suggested that AccGRF1 is involved in A. c. cerana glyphosate stress tolerance. This study reveals the functions of Fox transcription factors in response to glyphosate stress and provides molecular insights into the regulation of glyphosate responses in honeybees.

AccsHSP21.7 enhances the antioxidant capacity of Apis cerana cerana

BACKGROUND

The widespread use of glyphosate has many adverse effects on Apis cerana cerana. Due to the incomplete understanding of the molecular mechanisms of glyphosate toxicity, there are no available methods for mitigating the threat of glyphosate to Apis cerana cerana. Small heat shock proteins (sHSPs) play an important role in resisting oxidative stress, but their mechanism of action in Apis cerana cerana remains unclear.

RESULTS

In this experiment, we cloned and identified AccsHSP21.7. Studies have shown that AccsHSP21.7 contains binding motifs for various transcription factors related to oxidative stress. Abiotic stresses induced the expression of AccsHSP21.7. Bacteriostatic testing of a recombinant AccsHSP21.7 protein proved that Escherichia coli overexpressing AccsHSP21.7 showed increased resistance to oxidative stress. Knocking down the AccsHSP21.7 gene caused significant damage to midgut cells, which seriously disrupted the antioxidant system in Apis cerana cerana and greatly increased mortality under glyphosate stress.

CONCLUSION

This study investigated the relationship between antioxidant regulation and the AccsHSP21.7 gene at the molecular level, and the results have guiding significance for the improvement of stress resistance in Apis cerana cerana. © 2023 Society of Chemical Industry.

Investigating the effects of glyphosate on the bumblebee proteome and microbiota

Glyphosate is one of the most widely used herbicides globally. It acts by inhibiting an enzyme in an aromatic amino acid synthesis pathway specific to plants and microbes, leading to the view that it poses no risk to other organisms. However, there is growing concern that glyphosate is associated with health effects in humans and an ever-increasing body of evidence that suggests potential deleterious effects on other animals including pollinating insects such as bees. Although pesticides have long been considered a factor in the decline of wild bee populations, most research on bees has focussed on demonstrating and understanding the effects of insecticides. To assess whether glyphosate poses a risk to bees, we characterised changes in survival, behaviour, sucrose solution consumption, the digestive tract proteome, and the microbiota in the bumblebee Bombus terrestris after chronic exposure to field relevant doses of technical grade glyphosate or the glyphosate-based formulation, RoundUp Optima+®. Regardless of source, there were changes in response to glyphosate exposure in important cellular and physiological processes in the digestive tract of B. terrestris, with proteins associated with oxidative stress regulation, metabolism, cellular adhesion, the extracellular matrix, and various signalling pathways altered. Interestingly, proteins associated with endocytosis, oxidative phosphorylation, the TCA cycle, and carbohydrate, lipid, and amino acid metabolism were differentially altered depending on whether the exposure source was glyphosate alone or RoundUp Optima+®. In addition, there were alterations to the digestive tract microbiota of bees depending on the glyphosate source No impacts on survival, behaviour, or food consumption were observed. Our research provides insights into the potential mode of action and consequences of glyphosate exposure at the molecular, cellular and organismal level in bumblebees and highlights issues with the current honeybee-centric risk assessment of pesticides and their formulations, where the impact of co-formulants on non-target organisms are generally overlooked.

Parasite and Pesticide Impacts on the Bumblebee (Bombus terrestris) Haemolymph Proteome

Pesticides pose a potential threat to bee health, especially in combination with other stressors, such as parasites. However, pesticide risk assessment tests pesticides in isolation from other stresses, i.e., on otherwise healthy bees. Through molecular analysis, the specific impacts of a pesticide or its interaction with another stressor can be elucidated. Molecular mass profiling by MALDI BeeTyping® was used on bee haemolymph to explore the signature of pesticidal and parasitic stressor impacts. This approach was complemented by bottom-up proteomics to investigate the modulation of the haemoproteome. We tested acute oral doses of three pesticides—glyphosate, Amistar and sulfoxaflor—on the bumblebee Bombus terrestris, alongside the gut parasite Crithidia bombi. We found no impact of any pesticide on parasite intensity and no impact of sulfoxaflor or glyphosate on survival or weight change. Amistar caused weight loss and 19–41% mortality. Haemoproteome analysis showed various protein dysregulations. The major pathways dysregulated were those involved in insect defences and immune responses, with Amistar having the strongest impact on these dysregulated pathways. Our results show that even when no response can be seen at a whole organism level, MALDI BeeTyping® can detect effects. Mass spectrometry analysis of bee haemolymph provides a pertinent tool to evaluate stressor impacts on bee health, even at the level of individuals.

The sulfoximine insecticide sulfoxaflor exposure reduces the survival status and disrupts the intestinal metabolism of the honeybee Apis mellifera

Honeybees (Apis mellifera) are indispensable pollinators in agricultural production, biodiversity conservation, and nutrients provision. The abundance and diversity of honeybees have been rapidly diminishing, possibly related to the extensive use of insecticides in ecosystems. Sulfoxaflor is a novel sulfoximine insecticide that, like neonicotinoids, acts as a competitive modulator of nicotinic acetylcholine receptors (nAChR) in insects. However, few studies have addressed the negative effects of sulfoxaflor on honeybees at environmentally relevant concentrations. In the present study, adult workers were fed a 50% (w/v) of sugar solution containing different concentrations (0, 0.05, 0.5 and 2.0 mg/L) of sulfoxaflor for two weeks consecutively. The survival rates, food intake, and body weight of the honeybees significantly decreased after continuous exposure at higher doses (0.5 and 2.0 mg/L) of sulfoxaflor when compared with the control. The change in the metabolites in the honeybee gut was determined using high-throughput non-targeted metabolomics on day 14 after sulfoxaflor treatment. The results revealed that 24 and 105 metabolites changed after exposure to 0.5 and 2.0 mg/L sulfoxaflor, respectively, compared with that of the control groups. A total of 12 changed compounds including pregenolone and glutathione were detected as potential biomarkers, which were eventually found to be enriched in pathways of the steroid hormone biosynthesis (p = 0.0001) and glutathione metabolism (p = 0.021). These findings provide a new perspective on the physiological influence of sulfoxaflor stress in honeybees.

Integrative Analysis of Transcriptome and Metabolome to Illuminate the Protective Effects of Didymin against Acute Hepatic Injury

Based on the multiomics analysis, this study is aimed at investigating the underlying mechanism of didymin against acute liver injury (ALI). The mice were administrated with didymin for 2 weeks, followed by injection with lipopolysaccharide (LPS) plus D-galactosamine (D-Gal) to induce ALI. The pathological examination revealed that didymin significantly ameliorated LPS/D-Gal-induced hepatic damage. Also, it markedly reduced proinflammatory cytokines release by inhibiting the TLR4/NF-κB pathway activation, alleviating inflammatory injury. A transcriptome analysis proved 2680 differently expressed genes (DEGs) between the model and didymin groups and suggested that the PI3K/Akt and metabolic pathways might be the most relevant targets. Meanwhile, the metabolome analysis revealed 67 differently expressed metabolites (DEMs) between the didymin and model groups that were mainly clustered into the glycerophospholipid metabolism, which was consistent with the transcriptome study. Importantly, a comprehensive analysis of both the omics indicated a strong correlation between the DEGs and DEMs, and an in-depth study demonstrated that didymin alleviated metabolic disorder and hepatocyte injury likely by inhibiting the glycerophospholipid metabolism pathway through the regulation of PLA2G4B, LPCAT3, and CEPT1 expression. In conclusion, this study demonstrates that didymin can ameliorate LPS/D-Gal-induced ALI by inhibiting the glycerophospholipid metabolism and PI3K/Akt and TLR4/NF-κB pathways.

20-Hydroxyecdysone and Receptor Interplay in the Regulation of Hemolymph Glucose Level in Honeybee (Apis mellifera) Larvae

The hormone 20-hydroxyecdysone (20E) and the ecdysone receptors (ECR and USP) play critical roles in the growth and metabolism of insects, including honeybees. In this study, we investigated the effect of 20E on the growth and development of honeybee larvae by rearing them in vitro and found reduced food consumption and small-sized pupae with increasing levels of 20E. A liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based analysis of widely targeted metabolomics was used to examine the changes in the metabolites after an exogenous 20E application to honeybee larvae and the underlying mechanisms. A total of 374 different metabolites were detected between the control group and the 20E treatment group, covering 12 subclasses. The most significant changes occurred in 7-day-old larvae, where some monosaccharides, such as D-Glucose and UDP-galactose, were significantly upregulated. In addition, some metabolic pathways, such as glycolysis/gluconeogenesis and galactose metabolism, were affected by the 20E treatment, suggesting that the 20E treatment disrupts the metabolic homeostasis of honeybee larvae hemolymph and that the response of honeybee larvae to the 20E treatment is dynamic and contains many complex pathways. Many genes involved in carbohydrate metabolism, including genes of the glycolysis and glycogen synthesis pathways, were downregulated during molting and pupation after the 20E treatment. In contrast, the expression levels of the genes related to gluconeogenesis and glycogenolysis were significantly increased, which directly or indirectly upregulated glucose levels in the hemolymph, whereas RNA interference with the 20E receptor EcR-USP had an opposite effect to that of the 20E treatment. Taken together, 20E plays a critical role in the changes in carbohydrate metabolism during metamorphosis.

Glyphosate-Based Herbicide Causes Cellular Alterations to Gut Epithelium of the Neotropical Stingless Bee Melipona quadrifasciata quadrifasciata (Hymenoptera: Meliponini)

Glyphosate-based herbicides (GBH) are the best-selling pesticides in Brazil, with hundreds of thousands of tons sold per year. There is no study investigating morphological alterations caused by GBH on the epithelium of the gut in bees. Here, we aimed to demonstrate effects of chronic ingestion of GBH in the midgut digestive cells of the Brazilian stingless bee Melipona quadrifasciata quadrifasciata Lepeletier 1836. We kept forager workers of M. quadrifasciata in laboratory conditions and fed on food contaminated with three different concentrations of GBH for 10 days, after which the midgut digestive cell structure and ultrastructure were analyzed. The presence of GBH in food did not affect food consumption, indicating that M. quadrifasciata bees do not reject food contaminated with GBH. As digestive cells of the midgut release apocrine secretion as a detoxication mechanism, we expected that the ingestion of food contaminated with GBH in the present study affect the height of midgut digestive cells. However, such reduction did not occur, probably because of the low-test concentrations. Although there were differences in digestive cell ultrastructure, ingestion of GBH impaired midgut digestive cell cohesion by disorganizing the smooth septate junctions between cells, which may probably be caused by the adjuvant "polyethoxylated tallow amine" present in the GBH. Previous studies demonstrated that GBH increase bees' sensibility to intestine infections, based on the present results we hypothesized that the loss of cell cohesion in the midgut epithelium favors pathogenic microbial infections and harms food absorption, increasing bees' mortality.

Review on effects of some insecticides on honey bee health

Insecticides, one of the main agrochemicals, are useful for controlling pests; however, the indiscriminate use of insecticides has led to negative effects on nontarget insects, especially honey bees, which are essential for pollination services. Different classes of insecticides, such as neonicotinoids, pyrethroids, chlorantraniliprole, spinosad, flupyradifurone and sulfoxaflor, not only negatively affect honey bee growth and development but also decrease their foraging activity and pollination services by influencing their olfactory sensation, memory, navigation back to the nest, flight ability, and dance circuits. Honey bees resist the harmful effects of insecticides by coordinating the expression of genes related to immunity, metabolism, and detoxification pathways. To our knowledge, more research has been conducted on the effects of neonicotinoids on honey bee health than those of other insecticides. In this review, we summarize the current knowledge regarding the effects of some insecticides, especially neonicotinoids, on honey bee health. Possible strategies to increase the positive impacts of insecticides on agriculture and reduce their negative effects on honey bees are also discussed.

Identification of the cuticle protein AccCPR2 gene in Apis cerana cerana and its response to environmental stress

Cuticular proteins (CPs) are known to play important roles in insect development and defence responses. The loss of CP genes can lead to changes in insect morphology and sensitivity to the external environment. In this study, we identified the AccCPR2 gene, which belongs to the CPR family (including the R&R consensus motif) of CPs, and explored its function in the response of Apis cerana cerana to adverse external stresses. Our results demonstrated that AccCPR2 was highly expressed in the late pupal stage and epidermis, and the expression of AccCPR2 may be induced or inhibited under different stressors. RNA interference experiments showed that knockdown of AccCPR2 reduced the activity of antioxidant enzymes, led to the accumulation of oxidative damage and suppressed the expression of several antioxidant genes. In addition, knockdown of AccCPR2 also reduced the pesticide resistance of A. cerana cerana. The overexpression of AccCPR2 in a prokaryotic system further confirmed its role in resistance to various stresses. In summary, AccCPR2 may play pivotal roles in the normal development and environmental stress response of A. cerana cerana. This study also enriched the theoretical knowledge of the resistance biology of bees.

Metabolomic analysis of honey bee (Apis mellifera L.) response to glyphosate exposure

Glyphosate is among the world's most commonly used herbicides in agriculture and weed control. The use of this agrochemical has unintended consequences on non-target organisms, such as honey bees (Apis mellifera L.), the Earth's most prominent insect pollinator. However, detailed understanding of the biological effects in bees in response to sub-lethal glyphosate exposure is still limited. In this study, 1H NMR-based metabolomics was performed to investigate whether oral exposure to an environmentally realistic concentration (7.12 mg L-1) of glyphosate affects the regulation of honey bee metabolites in 2, 5, and 10 days. On Day 2 of glyphosate exposure, the honey bees showed significant downregulation of several essential amino acids, including leucine, lysine, valine, and isoleucine. This phenomenon indicates that glyphosate causes an obvious metabolic perturbation when the honey bees are subjected to the initial caging process. The mid-term (Day 5) results showed negligible metabolite-level perturbation, which indicated the low glyphosate impact on active honeybees. However, the long-term (Day 10) data showed evident separation between the control and experimental groups in the principal component analysis (PCA). This separation is the result of the combinatorial changes of essential amino acids such as threonine, histidine, and methionine, while the non-essential amino acids glutamine and proline as well as the carbohydrate sucrose were all downregulated. In summary, our study demonstrates that although no significant behavioral differences were observed in honey bees under sub-lethal doses of glyphosate, metabolomic level perturbation can be observed under short-term exposure when met with other environmental stressors or long-term exposure.

Population Structure and Genetic Diversity of Chinese Honeybee (Apis Cerana Cerana) in Central China

Central China has a rich terrain with a temperate monsoon climate and varied natural environments for the Chinese honeybee (Apis cerana cerana). However, little comprehensive research on population genetic diversity has been done in this area. A population survey of the structure and genetic diversity of Apis cerana cerana in this area is deeply needed for understanding adaptation to variable environments and providing more references for the protection of honeybee biodiversity. In this study, we present a dataset of 72 populations of Chinese honeybees collected from nine sites by whole genome sequencing in Central China. We obtained 2,790,214,878 clean reads with an average covering a depth of 22×. A total of 27,361,052 single nucleotide polymorphisms (SNPs) were obtained by mapping to the reference genome with an average mapping rate of 93.03%. Genetic evolution analysis was presented via the population structure and genetic diversity based on the datasets of SNPs. It showed that Apis cerana cerana in plains exhibited higher genetic diversity than in mountain areas. The mantel test between Apis cerana cerana groups revealed that some physical obstacles, especially the overurbanization of the plains, contributed to the differentiation. This study is conducive to elucidating the evolution of Apis cerana in different environments and provides a theoretical basis for investigating and protecting the Chinese honeybee.

Mushrooms Do Produce Flavonoids: Metabolite Profiling and Transcriptome Analysis of Flavonoid Synthesis in the Medicinal Mushroom Sanghuangporus baumii

Mushrooms produce a large number of medicinal bioactive metabolites with antioxidant, anticancer, antiaging, and other biological activities. However, whether they produce flavonoids and, if so, how they synthesize them remains a matter of some debate. In the present study, we combined flavonoid-targeted metabolomics and transcriptome analysis to explore the flavonoid synthesis in the medicinal mushroom Sanghuangporus baumii. The S. baumii synthesized 81 flavonoids on a chemically defined medium. The multiple classes of flavonoids present were consistent with the biosynthetic routes in plants. However, paradoxically, most of the genes that encode enzymes involved in the flavonoid biosynthetic pathway are missing from S. baumii. Only four genes related to flavonoid synthesis were found in S. baumii, among which phenylalanine ammonia lyase gene (PAL) is a key gene regulating flavonoid synthesis, and overexpression of SbPAL increases the accumulation of flavonoids. These results suggest that the flavonoid synthesis pathway in S. baumii is different from that in known plants, and the missing genes may be replaced by genes from the same superfamilies but are only distantly related. Thus, this study provides a novel method to produce flavonoids by metabolic engineering using mushrooms.

The effects of beta-cypermethrin, chlorbenzuron, chlorothalonil, and pendimethalin on Apis mellifera ligustica and Apis cerana cerana larvae reared in vitro

BACKGROUND

Declines in bee populations and diversity have drawn international attention. The long-term use of chemical pesticides has affected bee behavior and physiology. This study aimed to investigate the effects of chronic exposure to four commonly used chemical pesticides (beta-cypermethrin, chlorbenzuron, chlorothalonil and pendimethalin) on the growth of Apis mellifera ligustica and Apis cerana cerana larvae reared in vitro.

RESULTS

Pesticide type and concentration were the main factors affecting honeybee fitness. Beta-cypermethrin and chlorbenzuron had chronic toxic effects on bee larvae. They reduced the fitness of A. m. ligustica and A. c. cerana even at low doses of 323.5 ng g^-1 for beta-cypermethrin and 62.6 ng g^-1 for chlorbenzuron in bee bread. The effects were positively associated with the dietary amounts of pesticides. By contrast, chlorothalonil and pendimethalin exposure did not affect bee larvae despite changes in enzyme activities. Caution is still needed with chlorothalonil, which led to a decrease in harvest adult bee numbers at a high dose (6937.2 ng g^-1 ). Furthermore, a difference in pesticide resistance was observed, suggesting that A. m. ligustica may tolerate toxic effects better than A. c. cerana.

CONCLUSION

This study sheds new light on chronic toxicity in bee larvae exposed to residues in bee bread. The results could guide the scientific and rational use of chemical pesticides to reduce the potential risks to A. m. ligustica and A. c. cerana. © 2021 Society of Chemical Industry.

Distinct molecular impact patterns of abamectin on Apis mellifera ligustica and Apis cerana cerana

The effects of insecticides on bee health are a topic of intensive research. Although abamectin is toxic to bees, the molecular impact of abamectin needs to be clarified. Here, we found that Apis cerana cerana exhibited a higher mortality rate when exposed to abamectin than Apis mellifera ligustica. In addition, A. cerana cerana had markedly higher numbers of differentially expressed genes (DEGs), differentially expressed proteins (DEPs) and differentially expressed metabolites (DEMs) than A. mellifera ligustica during exposure to abamectin. These results indicate that abamectin exposure exerts stronger effects on A. cerana cerana than on A. mellifera ligustica. In addition, six DEGs, two DEPs and two DEMs overlapped between the two bee species under abamectin exposure; however, some genes or proteins from the zinc finger protein, superoxide dismutase and peroxiredoxin families and the energy metabolism pathway were only unregulated in A. cerana cerana, which indicates a significant difference in the impact of abamectin on the two bee species. Despite these differences, several of the same gene families, such as heat shock proteins, cytochrome P450, odorant-binding proteins and cuticle proteins, and pathways, including the carbohydrate metabolism, immune system, lipid metabolism, amino acid metabolism, sensory system, locomotion and development pathways, were influenced by abamectin exposure in both A. cerana cerana and A. mellifera ligustica. Together, our results indicate that abamectin causes adverse effects on bees and thus poses a risk to bee populations and that abamectin exposure affects A. cerana cerana more strongly than A. mellifera ligustica. These findings improve our understanding of the behavioural and physiological effects of abamectin on bees.

Integrated metabolomics and transcriptomics analysis reveals new biomarkers and mechanistic insights on atrazine exposures in MCF‑7 cells

Atrazine (ATZ) is a widely used herbicide worldwide and is a long-suspected endocrine-disrupting chemical. However, most endocrine-disrupting toxicity studies on ATZ have been based on animal models and those investigating inner mechanisms have only focused on a few genes. Therefore, the possible link between ATZ and endocrine-disrupting toxicity is still unclear. In this study, multi-omics and molecular biology techniques were used to elucidate the possible molecular mechanisms underlying the effect of ATZ exposure on MCF-7 proliferation at environmentally relevant concentrations. Our study is the first report on ATZ-induced one carbon pool by folate metabolic disorder in MCF-7 cells. A concentration of 1 μM ATZ yielded the highest cell viability and was selected for further mechanistic studies. A total of 34 significantly changed metabolites were identified based on metabolomic analysis, including vitamins, amino acids, fatty acids, and corresponding derivatives. Folate and pyridoxal have potential as biomarkers of ATZ exposure. One carbon pool by folate metabolic pathway was identified based on metabolic pathway analysis of the significantly altered pathways. Moreover, FTCD and MTHFD related to this pathway were further identified based on transcriptomic analysis and protein assays. Folate and different forms of 5,6,7,8-tetrahydrofolate, which participate in purine synthesis and associate with methyl groups (SOPC, arachidonic acid, and L-tryptophan) in one carbon pool by the folate metabolic pathway, potentially promote MCF-7 cell proliferation. These findings on the key metabolites and regulation of the related differentially expressed genes in folate metabolism will shed light on the mechanism of MCF-7 cell proliferation after ATZ exposure. Overall, this study provides new insights into the mechanistic understanding of toxicity caused by endocrine-disrupting chemicals.

Stereoselective toxicity mechanism of neonicotinoid dinotefuran in honeybees: New perspective from a spatial metabolomics study

Development of stereoisomeric neonicotinoid pesticides with lower toxicity is key to preventing global population declines of honeybees, whereas little is known about the in situ metabolic regulation of honeybees in response to stereoisomeric pesticides. Herein, we demonstrate an integrated mass spectrometry imaging (MSI) and untargeted metabolomics method to disclose disturbed metabolic expression levels and spatial differentiation in honeybees (Apis cerana) associated with stereoisomeric dinotefuran. This method affords a metabolic network mapping capability regarding a wide range of metabolites involved in multiple metabolic pathways in honeybees. Metabolomics results indicate more metabolic pathways of honeybees can be significantly affected by S-(+)-dinotefuran than R-(-)-dinotefuran, such as tricarboxylic acid (TCA) cycle, glyoxylate and dicarboxylate metabolism, and various amino acid metabolisms. MSI results demonstrate the cross-regulation and spatial differentiation of crucial metabolites involved in the TCA cycle, purine, glycolysis, and amino acid metabolisms within honeybees. Taken together, the integrated MSI and metabolomics results indicated the higher toxicity of S-(+)-dinotefuran arises from metabolic pathway disturbance and its inhibitory role in the energy metabolism, resulting in significantly reduced degradation rates of detoxification mechanisms. From the view of spatial metabolomics, our findings provide novel perspectives for the development and applications of pure chiral agrochemicals.

Glyphosate induces immune dysregulation in honey bees

BACKGROUND

Similar to many other animals, the honey bee Apis mellifera relies on a beneficial gut microbiota for regulation of immune homeostasis. Honey bees exposed to agrochemicals, such as the herbicide glyphosate or antibiotics, usually exhibit dysbiosis and increased susceptibility to bacterial infection. Considering the relevance of the microbiota-immunity axis for host health, we hypothesized that glyphosate exposure could potentially affect other components of the honey bee physiology, such as the immune system.

RESULTS

In this study, we investigated whether glyphosate, besides affecting the gut microbiota, could compromise two components of honey bee innate immunity: the expression of genes encoding antimicrobial peptides (humoral immunity) and the melanization pathway (cellular immunity). We also compared the effects of glyphosate on the bee immune system with those of tylosin, an antibiotic commonly used in beekeeping. We found that both glyphosate and tylosin decreased the expression of some antimicrobial peptides, such as apidaecin, defensin and hymenoptaecin, in exposed honey bees, but only glyphosate was able to inhibit melanization in the bee hemolymph.

CONCLUSIONS

Exposure of honey bees to glyphosate or tylosin can reduce the abundance of beneficial gut bacteria and lead to immune dysregulation.

Effects of glyphosate exposure on honeybees

Honeybees show an important pollination ability and play vital roles in improving crop yields and increasing plant genetic diversity, thereby generating tremendous economic benefits for humans. However, honeybee survival is affected by a number of biological and abiotic stresses, including the effects of fungi, bacteria, viruses, parasites, and especially agrochemicals. Glyphosate, a broad-spectrum herbicide that is primarily used for weed control in agriculture, has been reported to have lethal and sublethal effects on honeybees. Here, we summarize recent advances in research on the effects of glyphosate on honeybees, including effects on their behaviors, growth and development, metabolic processes, and immune defense, providing a detailed reference for studying the mechanism of action of pesticides. Furthermore, we provide possible directions for future research on glyphosate toxicity to honeybees.

Transcriptomic analysis of salivary gland and proteomic analysis of oral secretion in Helicoverpa armigera under cotton plant leaves, gossypol, and tannin stresses

Gossypol and tannin are involved in important chemical defense processes in cotton plants. In this study, we used transcriptomics and proteomics to explore the changes in salivary gland functional genes and oral secretion (OS) proteins after feeding with artificial diet (containing gossypols and tannins) and cotton plant leaves. We found that dietary cotton plant leaves, gossypols and tannins exerted adverse impacts on the genes that regulated the functions of peptidase, GTPase, glycosyl hydrolases in the salivary glands of the Helicoverpa armigera (H. armigera). However, GST, UGT, hydrolases, and lipase genes were up-regulated to participate in the detoxification and digestive of H. armigera. The oral secretory proteins of H. armigera were significantly inhibited under the stress of gossypol and tannin, such as enzyme activity, but some proteins (such as PZC71358.1) were up-regulated and involved in immune and digestive functions. The combined analysis of transcriptomics and metabolomics showed a weak correlation, and the genes and proteins involved were mainly in digestive enzyme activities. Our work clarifies the deleterious physiological impacts of gossypols and tannins on H. armigera and reveals the mechanism by which H. armigera effectively mitigate the phytotoxic effects through detoxification and immune systems.

Glyphosate inhibits melanization and increases susceptibility to infection in insects

Melanin, a black-brown pigment found throughout all kingdoms of life, has diverse biological functions including UV protection, thermoregulation, oxidant scavenging, arthropod immunity, and microbial virulence. Given melanin's broad roles in the biosphere, particularly in insect immune defenses, it is important to understand how exposure to ubiquitous environmental contaminants affects melanization. Glyphosate-the most widely used herbicide globally-inhibits melanin production, which could have wide-ranging implications in the health of many organisms, including insects. Here, we demonstrate that glyphosate has deleterious effects on insect health in 2 evolutionary distant species, Galleria mellonella (Lepidoptera: Pyralidae) and Anopheles gambiae (Diptera: Culicidae), suggesting a broad effect in insects. Glyphosate reduced survival of G. mellonella caterpillars following infection with the fungus Cryptococcus neoformans and decreased the size of melanized nodules formed in hemolymph, which normally help eliminate infection. Glyphosate also increased the burden of the malaria-causing parasite Plasmodium falciparum in A. gambiae mosquitoes, altered uninfected mosquito survival, and perturbed the microbial composition of adult mosquito midguts. Our results show that glyphosate's mechanism of melanin inhibition involves antioxidant synergy and disruption of the reaction oxidation-reduction balance. Overall, these findings suggest that glyphosate's environmental accumulation could render insects more susceptible to microbial pathogens due to melanin inhibition, immune impairment, and perturbations in microbiota composition, potentially contributing to declines in insect populations.

Is glyphosate toxic to bees? A meta-analytical review

Glyphosate (GLY) is an herbicide widely used in agriculture. First considered as non-toxic or slightly toxic to bees, GLY and its different formulations have shown, more recently, to affect negatively the survival, development and behavior of these insects, even when used in doses and concentrations recommended by the manufacturer. Thus, the results of research on the toxicity of GLY to bees are often conflicting, which makes a meta-analysis interesting for data integration, generating a statistically reliable result. Therefore, this study aimed to evaluate the GLY effects on mortality of bees through a meta-analysis. For this, a search was carried out in the databases Web of Science, CAPES (Coordination for the Improvement of Higher Education Personnel - Brazil), Scopus, and PubMed. Papers that evaluated the effect of GLY on bee mortality published between 1945 and October 2020, were considered. After obtaining the data, R software was used to perform the meta-analytical tests. Sixteen papers on mortality were selected with 34 data sets. Most of the sets demonstrated differences between the control and experimental groups, showing that the treatments with GLY caused higher mortality of bees. The results considering the methodology used (ingestion or contact), the phase of the biological cycle (adults or larvae), and the dose (ecologically relevant dose and recommended by the manufacturer) were different when compared with their respective control groups. Therefore, GLY can be considered toxic to bees. It is important to emphasize that this meta-analysis identified that papers assessing the toxicity of GLY to bees are still scarce, for both lethal and sublethal effects, mainly for stingless and solitary bee species.

Impact of Chronic Exposure to Sublethal Doses of Glyphosate on Honey Bee Immunity, Gut Microbiota and Infection by Pathogens

Glyphosate is the most used pesticide around the world. Although different studies have evidenced its negative effect on honey bees, including detrimental impacts on behavior, cognitive, sensory and developmental abilities, its use continues to grow. Recent studies have shown that it also alters the composition of the honey bee gut microbiota. In this study we explored the impact of chronic exposure to sublethal doses of glyphosate on the honey bee gut microbiota and its effects on the immune response, infection by Nosema ceranae and Deformed wing virus (DWV) and honey bee survival. Glyphosate combined with N. ceranae infection altered the structure and composition of the honey bee gut microbiota, for example by decreasing the relative abundance of the core members Snodgrassella alvi and Lactobacillus apis. Glyphosate increased the expression of some immune genes, possibly representing a physiological response to mitigate its negative effects. However, this response was not sufficient to maintain honey bee health, as glyphosate promoted the replication of DWV and decreased the expression of vitellogenin, which were accompanied by a reduced life span. Infection by N. ceranae also alters honey bee immunity although no synergistic effect with glyphosate was observed. These results corroborate previous findings suggesting deleterious effects of widespread use of glyphosate on honey bee health, and they contribute to elucidate the physiological mechanisms underlying a global decline of pollination services.

Cloning and expression studies on glutathione S-transferase like-gene in honey bee for its role in oxidative stress

Glutathione S-transferases (GSTs) constitute an important multifunctional enzyme family that plays vital roles in cellular detoxification and protecting organisms against oxidative stress caused by reactive oxygen species (ROS). In this study, we isolated a GST-like gene from Apis cerana cerana (AccGSTL) and investigated its antioxidant functions under stress conditions. We found that AccGSTL belongs to the Sigma class of GSTs. Real-time quantitative PCR and western blotting analyses showed that the mRNA and protein levels of AccGSTL were altered in response to oxidative stress caused by various external stimuli. In addition, a heterologous expression analysis showed that AccGSTL overexpression in Escherichia coli (E. coli) cells enhanced resistance to oxidative stress. After AccGSTL silencing with RNA interference (RNAi) technology, the expression of some antioxidant genes was inhibited, and the enzymatic activities of POD, CAT, and SOD were decreased. In conclusion, these data suggest that AccGSTL may be involved in antioxidant defense under adverse conditions in A. cerana cerana.

Integrated Analysis of the Transcriptome and Metabolome Reveals Genes Involved in Terpenoid and Flavonoid Biosynthesis in the Loblolly Pine (Pinus taeda L.)

Loblolly pine (Pinus taeda L.) is an important tree for afforestation with substantial economic and ecological value. Many metabolites with pharmacological activities are present in the tissues of P. taeda. However, the biosynthesis regulatory mechanisms of these metabolites are poorly understood. In the present study, transcriptome and metabolome analyses were performed on five tissues of P. taeda. A total of 40.4 million clean reads were obtained and assembled into 108,663 unigenes. These were compared with five databases, revealing 39,576 annotated unigenes. A total of 13,491 differentially expressed genes (DEGs) were observed in 10 comparison groups. Of these, 487 unigenes exhibited significantly different expressions in specific tissues of P. taeda. The DEGs were explored using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes metabolic pathway analysis. We identified 343 and 173 candidate unigenes related to the biosynthesis of terpenoids and flavonoids, respectively. These included 62 R2R3-MYB, 30 MYB, 15 WRKY, seven bHLH, seven ERF, six ZIP, five AP2, and one WD40 genes that acted as regulators in flavonoid and/or terpenoid biosynthesis. Additionally, metabolomics analysis detected 528 metabolites, among which 168 were flavonoids. A total of 493 differentially accumulated metabolites (DAMs) were obtained in 10 comparison groups. The 3,7-Di-O-methyl quercetin was differentially accumulated in all the comparison groups. The combined transcriptome and metabolome analyses revealed 219 DEGs that were significantly correlated with 45 DAMs. Our study provides valuable genomic and metabolome information for understanding P. taeda at the molecular level, providing a foundation for the further development of P. taeda-related pharmaceutical industry.