Proteomic analysis of the hepatotoxicity of Microcystis aeruginosa in adult zebrafish (Danio rerio) and its potential mechanisms

Bloom-forming planktonic Microcystis and benthic Oscillatoria-induced oxidative stress and inflammatory responses in juvenile silver carp and bighead carp

As global warming and water eutrophication, the multiple proliferation of harmful cyanobacteria can form algal blooms and cause serious ecological problems. In recent years, the large-scale and persistent cyanobacterial blooms occur frequently worldwide and have attracted widespread attention due to the harmful impacts. Among these harmful bloom-forming cyanobacteria, the ecological and toxicological impacts of planktonic cyanobacteria have been extensively studied. However, research on the ecological risks and adverse effects of harmful benthic cyanobacteria is lagging. Filter-feeding fish could suffer from more toxic stimuli than other fish due to their special feeding habits. To investigate and compare the complex toxic effects of different kinds of harmful cyanobacteria on fish, three different-sized (i.e. small, medium, and large) juvenile silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) were exposed to cyanobacterial blooms-related density (1 × 106 cells/mL) of Microcystis aeruginosa (i.e. generating microcystins) and Oscillatoria sp. (i.e. generating cylindrospermopsin) for 3 d, after which biomarkers of oxidative stress and inflammation in the liver of fish were detected. The silver carp and bighead carp can effectively ingest Microcystis cells but cannot effectively ingest Oscillatoria cells through the measurement of the levels of cyanotoxins. Both Microcystis and Oscillatoria cells can induce different levels of oxidative stress and inflammatory responses in the liver of these juvenile filter-feeding fish via altering the biochemical parameters of the antioxidant system (e.g. superoxide dismutase activity) and immune system (e.g. interleukin-1β level). Therefore, our research identified potential data gaps that how the different types of cyanobacteria induce toxic effects in the liver of juvenile filter-feeding fish in a short time. This study contributes to a better understanding of the short-term adverse effects of different cyanobacterial species on juvenile fish, suggesting that the benthic toxic cyanobacteria-induced ecological and health risks require further attention.

Environmental signatures and fish proteomics: A multidisciplinary study to identify the major stressors in estuaries located in French agricultural watersheds

Watersheds and estuaries are impacted by multiple anthropogenic stressors that affect their biodiversity and functioning. Assessing their ecological quality has consequently remained challenging for scientists and stakeholders. In this paper, we propose a multidisciplinary approach to identify the stressors in seven small French estuaries located in agricultural watersheds. We collected data from landscape (geography, hydrobiology) to estuary (pollutant chemistry) and fish individual scales (environmental signatures, proteomics). This integrative approach focused on the whole hydrosystems, from river basins to estuaries. To characterize each watershed, we attempted to determine the land use considering geographic indicators (agricultural and urbanised surfaces) and landscape patterns (hedges density and riparian vegetation). Juveniles of European flounder (Platichthys flesus) were captured in September, after an average residence of five summer months in the estuary. Analyses of water, sediments and biota allowed to determine the concentrations of dissolved inorganic nitrogen species, pesticides and trace elements in the systems. Environmental signatures were also measured in flounder tissues. These environmental parameters were used to establish a typology of the watersheds. Furthermore, data from proteomics on fish liver were combined with environmental signatures to determine the responses of fish to stressors in their environments. Differential protein abundances highlighted a dysregulation related to the detoxification of xenobiotics (mainly pesticides) in agricultural watersheds, characterized by intensive cereal and vegetable crops and high livestock. Omics also revealed a dysregulation of proteins associated with the response to hypoxia and heat stress in some estuaries. Furthermore, we highlighted a dysregulation of proteins involved in urea cycle, immunity and metabolism of fatty acids in several systems. Finally, the combination of environmental and molecular signatures appears to be a relevant method to identify the major stressors operating within hydrosystems.

Harmful planktonic Microcystis and benthic Oscillatoria-induced toxicological effects on the Asian clam (Corbicula fluminea): A survey on histopathology, behavior, oxidative stress, apoptosis and inflammation

Cyanobacterial blooms are worldwide distributed and threaten aquatic ecosystems and public health. The current studies mainly focus on the adverse impacts of planktonic cyanobacteria or pure cyanotoxins, while the benthic cyanobacteria-induced ecotoxic effects are relatively lacking. The cyanobacterial cell-induced toxic effects on aquatic organisms might be more serious and complex than the pure cyanotoxins and crude extracts of cyanobacteria. This study explored the chronic effects of toxin-producing planktonic Microcystis aeruginosa (producing microcystin) and benthic Oscillatoria sp. (producing cylindrospermopsin) on the behaviors, tissue structures, oxidative stress, apoptosis, and inflammation of the Asian clams (Corbicula fluminea) under 28-d exposure. The data showed that both M. aeruginosa and Oscillatoria sp. can decrease the behaviors associated with the feeding activity and induce tissue damage (i.e. gill and digestive gland) in clams. Furthermore, two kinds of cyanobacteria can alter the antioxidant enzyme activities and increase antioxidant, lipid oxidation product, and neurotransmitter degrading enzyme levels in clams. Moreover, two kinds of cyanobacteria can activate apoptosis-related enzyme activities and enhance the proinflammatory cytokine levels of clams. In addition, two kinds of cyanobacteria can disturb the transcript levels of genes linked with oxidative stress, apoptosis, and inflammation. These results suggested harmful cyanobacteria can threaten the survival and health of clams, while the benthic cyanobacteria-induced adverse effects deserve more attention. Our finding also indicated that it is necessary to focus on the entire algal cell-induced ecotoxicity when concerning the ecological impacts of cyanobacterial blooms.

Can a shift in dominant species of Microcystis alter growth and reproduction of waterfleas?

The bloom-forming species Microcystis wesenbergii and M. aeruginosa occur in many lakes globally, and may exhibit alternating blooms both spatially and temporally. As environmental changes increase, cyanobacteria bloom in more and more lakes and are often dominated by M. wesenbergii. The adverse impact of M. aeruginosa on co-existing organisms including zooplanktonic species has been well-studied, whereas studies of M. wesenbergii are limited. To compare effects of these two species on zooplankton, we explored effects of exudates from different strains of microcystin-producing M. aeruginosa (Ma905 and Ma526) and non-microcystin-producing M. wesenbergii (Mw908 and Mw929), on reproduction by the model zooplankter Daphnia magna in both chronic and acute exposure experiments. Specifically, we tested physiological, biochemical, molecular and transcriptomic characteristics of D. magna exposed to Microcystis exudates. We observed that body length and egg and offspring number of the daphnid increased in all treatments. Among the four strains tested, Ma526 enhanced the size of the first brood, as well as total egg and offspring number. Microcystis exudates stimulated expression of specific genes that induced ecdysone, juvenile hormone, triacylglycerol and vitellogenin biosynthesis, which, in turn, enhanced egg and offspring production of D. magna. Even though all strains of Microcystis affected growth and reproduction, large numbers of downregulated genes involving many essential pathways indicated that the Ma905 strain might contemporaneously induce damage in D. magna. Our study highlights the necessity of including M. wesenbergii into the ecological risk evaluation of cyanobacteria blooms, and emphasizes that consequences to zooplankton may not be clear-cut when assessments are based upon production of microcystins alone.

Tea polyphenols inhibit blooms caused by eukaryotic and prokaryotic algae

With changes in global climate, blooms are becoming more frequent and difficult to control. Therefore, the selection of algal suppressor agents with effective inhibition and environmental safety is of paramount importance. One of the main treatment strategies is to inhibit the release of harmful algal toxins. Tea polyphenols (TP) are natural products that have been widely used in medicine, the environment, and other fields due to their antibacterial and antioxidant properties. To investigate their potential application in the treatment of algal blooms, TP were applied to three different microalgae. TP exhibited strong inhibitory effects towards all three microalgae. They stimulate the accumulation of ROS in algal cells, leading to lipid peroxidation and subsequent damage to the cell membrane, resulting in the rupture and necrosis of Cyclotella sp. and Chlorella vulgaris cells. Remarkably, it was observed that lower concentrations of TP exhibited the ability to induce apoptosis in M. aeruginosa cells without causing any structural damage. This outcome is particularly significant as it reduces the potential risk of microcystin release resulting from cell rupture. Overall, blooms dominated by different algae can be treated by adjusting the concentration of TP, a new algal suppressor, indicating strong potential treatment applications.

Is a lower-toxicity strain of Microcystis aeruginosa really less toxic?

Some well-known hazards of blooming cyanobacteria are caused by toxic metabolites such as microcystins (MCs), though many other bioactive chemicals of unknown toxicity are present in their exudates. It is also unclear whether toxicity of cyanobacterial cells depends on growth phases in the life cycle. In this study, we compared toxicity to Daphnia magna of Microcystis aeruginosa - a common cyanobacterial species - exudates (MaE) from two MC-producing strains over both exponential growth and stationary phases in acute and chronic experiments. Specifically, we assessed mitochondrial dysfunction, oxidative stress and lipid peroxidation, and filtering activity and heartbeat rate of Daphnia exposed to MaE. All MaE treatments induced common characteristics of Microcystis toxicity including disorder in the mitochondrial membrane and aberrant heart rate. MaE from cells at stationary growth phase were more toxic than those at exponential phase. Surprisingly, the MC-lower strain had higher toxicity than MC-higher one. Microcystis at different stage of blooms may differentially affect waterfleas owing to variable MaE-induced physiological dysfunction, abundance and grazing rate. Our study suggested that Microcystis strains with lower microcystin-producing ability might release other detrimental chemicals and should not be ignored in harmful bloom monitoring.

Spatiotemporal patterns and threshold of chlorophyll-a in Lake Taihu based on microcystins

Chlorophyll-a (Chl-a) is considered as an indicator of phytoplankton biomass dynamically reflecting the growth of algae. Therefore, determination of Chl-a threshold is of vital importance to the health of aquatic ecosystems and drinking water security. This research is aimed to investigate the spatial and temporal distributions of Chl-a and microcystin (MC) concentrations using Geographic Information System (GIS) and identify the Chl-a threshold in Lake Taihu based on available guideline values of MCs. Nearly, the same characteristics of spatiotemporal variation of Chl-a and MCs were observed in Lake Taihu. Overall, the lakewide distributions of Chl-a and MCs were highly variable over time and space. The Chl-a concentration in the winter and spring was relatively low, and gradually increasing in summer and autumn, with the maximum concentration observed in August. But the maximum MCs concentration appeared in October, 2 months lagging behind the Chl-a. The highest annual average Chl-a and MCs concentrations were observed in Zhushan Bay, Meiliang Bay, and Gonghu Bay in northwest of Lake Taihu, following by West Zone and Center Zone. Dongtaihu Bay, East Zone, and South Zone always present good water quality. Referring to the guideline value of MCs, the Chl-a threshold was determined as 10-15 mg·m^-3 based on the linear regression correlation between Chl-a and MCs. The establishment of Chl-a threshold is useful for eutrophication control, water quality management, and drinking water utilities in developing water safety plans.

Mechanism for the Potential Inhibition Effect of Microcystin-LR Disinfectant By-Products on Protein Phosphatase 2A

The secondary contamination of microcystin disinfection by-products (MC-DBPs) is of concern due to the residual structure similar to their original toxin. Based on identification and preparation, the potential inhibition effect of typical MCLR-DBPs (associated with the oxidation of Adda⁵) on PP2A was confirmed in the sequence of MCLR > P1 > P4 > P3 ≈ P2 > P7 ≈ P6 ≈ P5 > P8. To elucidate the molecular mechanism underlying the inhibition effect, the interaction models for typical MCLR-DBPs and PP2A were constructed using a modeling-based-on-ligand-similarity approach, and the candidate interaction parameters between typical MCLR-DBPs and PP2A were obtained by molecular docking. By analyzing the correlation between inhibition data and candidate interaction parameters, the key interaction parameters were filtered as hydrogen bonds "Adda⁵"←Asn₁₁₇, "Adda⁵"←His₁₁₈, MeAsp³←Arg₈₉, Arg⁴←Arg₂₁₄, Arg⁴→Pro₂₁₃; ionic bonds Glu⁶-Arg₈₉, Asp₈₅-Mn₁²⁺, Asp₅₇-Mn₂²⁺; and metal bonds Glu⁶-Mn₁²⁺, Glu⁶-Mn₂²⁺. With the gradual intensification of chlorination, Adda⁵ was destroyed to varying degrees. The key interactions changed correspondingly, resulting in the discrepant inhibition effects of typical MCLR-DBPs on PP2A.

Evaluation of phoxim toxicity on aquatic and zebrafish intestinal microbiota by metagenomics and 16S rRNA gene sequencing analysis

Phoxim is one of the main organophosphorus pesticides used in agricultural production. However, little information is known about how it affects the aquatic microbial community and the intestinal microbiota of fish. Herein, we utilized shotgun metagenomics and 16S rRNA gene sequencing to reveal the aquatic eco-risk of phoxim. Seven days of phoxim exposure significantly changed the composition of aquatic microbial community, obliterated the interactions between microorganisms, and thus reduced the complexity and stability of the microbial community. During long-time exposure (i.e., 14 days), most of the ecological functions were restored due to the redundancy of the microbial community. However, phoxim exposure promoted the dissemination of elfamycin resistance gene. The zebrafish gut microbial community also recovered from a temporary ecological disorder of aquatic microbiota, but phoxim continually affected zebrafish growth and swimming behavior. Overall, our results demonstrated that phoxim exposure significantly changed the structure and function of the microbial community and displayed a negative impact on freshwater ecosystems in a short exposure time.

Algicide capacity of Paucibacter aquatile DH15 on Microcystis aeruginosa by attachment and non-attachment effects

The excessive proliferation of Microcystis aeruginosa can lead to ecological damage, economic losses, and threaten animal and human health. For controlling Microcystis blooms, microorganism-based methods have attracted much attention from researchers because of their eco-friendliness and species-specificity. Herein, we first found that a Paucibacter strain exhibits algicidal activity against M. aeruginosa and microcystin degradation capability. The algicidal activity of DH15 (2.1 × 10⁴ CFU/ml) against M. aeruginosa (2 × 10⁶ cells/ml) was 94.9% within 36 h of exposure. DH15 also degraded microcystin (1.6 mg/L) up to 62.5% after 72 h. We demonstrated that the algicidal activity of DH15 against M. aeruginosa can be mediated by physical attachment and indirect attack: (1) Both washed cells and cell-free supernatant could kill M. aeruginosa efficiently; (2) Treatment with DH15 cell-free supernatants caused oxidative stress, altered the fatty acid profile, and damaged photosynthetic system, carbohydrate, and protein metabolism in M. aeruginosa. The combination of direct and indirect attacks supported that strain DH15 exerts high algicidal activity against M. aeruginosa. The expression of most key genes responsible for photosynthesis, antioxidant activity, microcystin synthesis, and other metabolic pathways in M. aeruginosa was downregulated. Strain DH15, with its microcystin degradation capacity, can overcome the trade-off between controlling Microcystis blooms and increasing microcystin concentration. Our findings suggest that strain DH15 possesses great potential to control outbreaks of Microcystis blooms.

Salt-alkalization may potentially promote Microcystis aeruginosa blooms and the production of microcystin-LR

The development of saline-alkali lands has contributed to the increasing discharge of alkaline salt-laden wastewater, which poses a threat to aquatic organisms. However, the comprehensive effect of alkaline salt on Microcystis aeruginosa, a harmful cyanobacterium, remains unclear. In this study, the growth, physiology, cell ultrastructure and production of microcystin-LR (MC-LR) in Microcystis aeruginosa exposed to four levels of alkaline salt stress were evaluated. The growth of Microcystis aeruginosa was stimulated at an electrical conductivity (EC) of 2.5 mS/cm compared to the control, as supported by the increased cell density, photosynthetic pigment and protein contents. Microcystis aeruginosa could tolerate a certain level of alkaline salt (i.e., EC of 5 mS/cm) via increasing photosynthetic pigment contents to protect cells from alkaline salt stress, but the antioxidant defence system and cell ultrastructure were not affected. When EC increased to 7.5 mS/cm, alkaline salt caused oxidative stress and toxicity in Microcystis aeruginosa, as evidenced by analysis of the integrated biomarker response (IBR). Furthermore, the photosynthetic pigment and protein contents decreased, and cell apoptosis associated with ultrastructural changes was observed. Therefore, we propose that EC of 7.5 mS/cm is a threshold for growth of Microcystis aeruginosa. Additionally, the intracellular MC-LR content was stimulated by alkaline salt, and the highest value was observed at EC of 2.5 mS/cm. The extracellular MC-LR content increased with the increasing alkaline salt concentration. When EC was 7.5 mS/cm, the extracellular MC-LR content was significantly higher than in the control and was associated with the upregulated mcyH gene. This study recommends that more attention should be paid to the risk of Microcystis aeruginosa bloom and microcystin-LR pollution in lakes located in salinization regions.

Proteomic analysis of zebrafish brain damage induced by Microcystis aeruginosa bloom

Cyanobacterial blooms constitute a global ecological problem that can seriously threaten human health. One of the most common bloom-forming cyanobacteria in freshwater is Microcystis aeruginosa, whose secretion of toxic substances (microcystins, MCs) have strong liver toxicity and endanger the health of exposed people through contaminated aquatic products and drinking water. However, few studies on the neurotoxicity of M. aeruginosa to zebrafish have simulated the process of an actual cyanobacterial bloom. In this study, we used the zebrafish (Danio rerio) as an effective model organism to study the acute neurotoxicity of M. aeruginosa, and to clarify its principal mechanism of action. A total of 82 upregulated and 26 downregulated proteins were detected by quantitative proteomics analysis in zebrafish brain after exposure to M. aeruginosa. Intriguingly, these proteins with changed expression were related to Synaptic vesicle cycle and terpenoid skeleton biosynthesis pathway, such as ACAT, STX1A, and V-ATPase. The obtained results uniformly indicated that the neurotoxicity of M. aeruginosa seriously damaged the neurotransmitter conduction in the nervous system and brain information storage and transmission of zebrafish and makes it more susceptible to neurological diseases. Our study provides a new perspective on the neurotoxicity risk of cyanobacterial blooms.

Effects of microcystin-producing and non-microcystin-producing Microcystis on the behavior and life history traits of Chironomus pallidivittatus

Species of the genus Microcystis are among the most notorious cyanobacteria in eutrophic lakes worldwide, with ability present adverse effects on many aquatic organisms. In the surface sediments, Microcystis can be ingested by benthic macroinvertebrates such as Chironomus. However, the potential negative effects of Microcystis on Chironomus life history traits remain unclear. In the present study, we investigated the effect of different Microcystis diets on specific behaviors (burrowing activity, locomotion ability) and life history traits of Chironomus pallidivittatus (Diptera, Chironomidae). We also studied the interactive effects of microcystin-producing M. aeruginosa and temperature (15, 20, and 25 °C) stress on chironomid larvae. The results showed that the inhibitory effect on the cumulative emergence and burrowing activity of larvae was more severe when they were fed M. aeruginosa among the three Microcystis diets groups. Locomotion ability (i.e., locomotor distance and velocity) and adult dry weight decreased significantly in the group fed M. aeruginosa. Locomotion was significantly inhibited and mortality increased when the larvae were fed a mixture of M. aeruginosa and M. wesenbergii, which may have been the result of additive or synergistic effect of the toxins. Under the stress of lower temperature, C. pallidivittatus larvae exhibited weaker locomotion and growth ability, and the emerging adults were mostly male. At both the lower and higher temperature conditions, M. aeruginosa cause cumulative emergence decreased, and sex ratio imbalance, which inhibited the reproduction of larvae from the population perspective. The fourth-instar larvae showed better adaption to Microcystis than did the other instars. This study thus highlights the adverse effects of microcystin-producing M. aeruginosa on Chironomus. It also provides a novel perspective on how environmental factors may influence the behavior and life history traits of chironomid larvae, and how they may respond to cyanobacterial blooms and global warming.

Development of an adverse outcome pathway for nanoplastic toxicity in Daphnia pulex using proteomics

Nanoplastics are a growing environmental and public health concern. However, the toxic mechanisms of nanoplastics are poorly understood. Here, we evaluated the effects of spherical polystyrene nanoplastics on reproduction of Daphnia pulex and analyzed the proteome of whole animals followed by molecular and biochemical analyses for the development of an adverse outcome pathway (AOP) for these contaminants of emerging concern. Animals were exposed to polystyrene nanoplastics (0, 0.1, 0.5, 1 and 2 mg/L) via water for 21 days. Nanoplastics negatively impacted cumulative offspring production. A total of 327 differentially expressed proteins (DEPs) were identified in response to nanoplastics which were further validated from gene expression and enzyme activity data. Based on these results, we propose an AOP for nanoplastics, including radical oxygen species production and oxidative stress as the molecular initiating event (MIE); followed by changes in specific signaling pathways (Jak-STAT, mTOR and FoxO) and in the metabolism of glutathione, protein, lipids, and molting proteins; with an end result of growth inhibition and decrease reproductive output. This study serves as a foundation for the development of a mechanistic understanding of nanoplastic toxicity in crustaceans and perhaps other aquatic organisms.

Potassium regulates the growth and toxin biosynthesis of Microcystis aeruginosa

Potassium (K⁺) is the most abundant cation in phytoplankton cells, but its impact on Microcystis aeruginosa (M. aeruginosa) has not been fully documented. This study presents evidence of how K⁺ availability affects the growth, oxidative stress and microcystin (MC) production of M. aeruginosa. The iTRAQ-based proteomic analysis revealed that during K⁺ deficiency, serious oxidative damage occurred and the photosynthesis-associated and ABC transporter-related proteins in M. aeruginosa were substantially downregulated. In the absence of K⁺, a 69.26% reduction in cell density was shown, and both the photosynthesis and iron uptake were depressed, which triggered a declined production of ATP and expression of MC synthetases genes (mcyA, B and D), and MC exporters (mcyH). Through the impairment of both the MC biosynthesis and MC transportation out of cells, K⁺ depletion caused an 85.89% reduction of extracellular MC content at the end of the study. However, with increasing in the available K⁺ concentrations, photosynthesis efficiency, the expression of ABC-transporter proteins, and the transcription of mcy genes displayed slight differences compared with those in the control group. This work represents evidence that K⁺ availability can regulate the physiological metabolic activity of M. aeruginosa and K⁺ deficiency leads to depressed growth and MC production in M. aeruginosa.

New Insights in Saccharomyces cerevisiae Response to the Cyanotoxin Microcystin-LR, Revealed by Proteomics and Gene Expression

Microcystins (MCs) are hepatotoxins produced by some cyanobacteria. They are cyclic peptides that inhibit the serine/threonine protein phosphatases (PPs) PP1 and PP2A, especially PP2A. The inhibition of PP2A triggers a series of molecular events, which are responsible for most MC cytotoxic and genotoxic effects on animal cells. It is also known that MCs induce oxidative stress in cells due to the production of reactive oxygen species (ROS). However, a complete characterization of the toxic effects of MCs is still not accomplished. This study aimed to clarify additional molecular mechanisms involved in MC-LR toxicity, using Saccharomyces cerevisiae as eukaryotic model organism. First, a shotgun proteomic analysis of S. cerevisiae VL3 cells response to 1 nM, 10 nM, 100 nM, and 1 μM MC-LR was undertaken and compared to the control (cells not exposed to MC-LR). This analysis revealed a high number of proteins differentially expressed related with gene translation and DNA replication stress; oxidative stress; cell cycle regulation and carbohydrate metabolism. Inference of genotoxic effects of S. cerevisiae VL3 cells exposed to different concentrations of MC-LR were evaluated by analyzing the expression of genes Apn1, Apn2, Rad27, Ntg1, and Ntg2 (from the Base Excision Repair (BER) DNA repair system) using the Real-Time RT-qPCR technique. These genes displayed alterations after exposure to MC-LR, particularly the Apn1/Apn2/Rad27, pointing out effects of MC-LR in the Base Excision Repair system (BER). Overall, this study supports the role of oxidative stress and DNA replication stress as important molecular mechanisms of MC-LR toxicity. Moreover, this study showed that even at low-concentration, MC-LR can induce significant changes in the yeast proteome and in gene expression.

Are we forgetting the "proteomics" in multi-omics ecotoxicology?

Proteomics plays a significant role in discerning the effects of chemical exposures in animal taxa. Multi-omics applications have become more pervasive in toxicology, however questions remain about whether proteomics is being utilized by the community to its full potential - are we placing too much stock in transcriptomics and other omics approaches for developing adverse outcome pathways? Proteins are more relevant than transcripts because they are direct mediators of the resulting phenotype. There is also rarely perfect stoichiometry between transcript and protein abundance and transcript abundance may not accurately predict physiologic response. Proteins direct all levels of phenotype: structural proteins dictate physical form, enzymes catalyze biochemical reactions, and proteins act as signaling proteins, antibodies, transporters, ion pumps, and transcription factors to control gene expression. Molecular initiating events (MIEs) of AOPs predominantly occur at the level of the protein (e.g. ligand-receptor binding) and proteomics can elucidate novel MIEs and mapping KEs in AOPs. This critical review highlights the need for proteomics in multi-omics studies in environmental toxicology and outlines steps required for inclusion and wider acceptance in chemical risk assessment. We also present case studies of multi-omics approaches that utilize proteomics and discuss some of the challenges and opportunities for proteomics in comparative ecotoxicology. Our intention is not to minimize the importance of other omics technologies, as each has strengths and limitations, but rather to encourage researchers to consider proteomics-based methods in multi-omics studies and AOP development.