Next Article in Journal
Inhibition of Microtubule Affinity Regulating Kinase 4 by Metformin: Exploring the Neuroprotective Potential of Antidiabetic Drug through Spectroscopic and Computational Approaches
Next Article in Special Issue
MALDI-MSI: A Powerful Approach to Understand Primary Pancreatic Ductal Adenocarcinoma and Metastases
Previous Article in Journal
Anticancer Inhibitors
Previous Article in Special Issue
Forensic Discrimination of Differentially Sourced Animal Blood Using a Bottom-Up Proteomics Based MALDI MS Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mass Spectrometry-Based Peptide Profiling of Haemolymph from Pterostichus melas Exposed to Pendimethalin Herbicide

1
Department of Chemistry and Chemical Technologies, University of Calabria, 87036 Arcavacata di Rende, Italy
2
Department of Biology, Ecology and Earth Science, University of Calabria, 87036 Arcavacata di Rende, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2022, 27(14), 4645; https://doi.org/10.3390/molecules27144645
Submission received: 27 June 2022 / Revised: 16 July 2022 / Accepted: 19 July 2022 / Published: 21 July 2022

Abstract

:
Pendimethalin-based herbicides are used worldwide for pre-emergence selective control of annual grasses and weeds in croplands. The endurance of herbicides residues in the environment has an impact on the soil biodiversity and fertility, also affecting non-target species, including terrestrial invertebrates. Carabid beetles are known as natural pest control agents in the soil food web of agroecosystems, and feed on invertebrates and weed seeds. Here, a mass spectrometry untargeted profiling of haemolymph is used to investigate Pterostichus melas metabolic response after to pendimethalin-based herbicide exposure. Mass spectrometric data are examined with statistical approaches, such as principal component analysis, for possible correlation with biological effects. Those signals with high correlation are submitted to tandem mass spectrometry to identify the associated biomarker. The time course exposure showed many interesting findings, including a significant downregulation of related to immune and defense peptides (M-lycotoxin-Ls4a, Peptide hormone 1, Paralytic peptide 2, and Serine protease inhibitor 2). Overall, the observed peptide deregulations concur with the general mechanism of uptake and elimination of toxicants reported for Arthropods.

Graphical Abstract

1. Introduction

Dinitroaniline herbicides, including trifluralin, ethalfluralin, oryzalin, and pendimethalin, have been used worldwide for pre-emergence selective control of annual grasses and weeds in croplands. Pendimethalin ([N-(1-ethylpropyl)-2,6-dinitro-3,4-xylidine], PND) acts on the polymerization of microtubules, interfering with cell division (mitosis) and, thus, inhibiting the development of roots and shoots in seedlings [1]. However, the endurance of PND residues in the environment have an impact on the growth the soil bacterial [2] and fungal [3], consequently altering the soil biodiversity and fertility [4]. Its sublethal effects have been documented in several non-target species including terrestrial invertebrates, such as springtails, earthworms, pill bugs [5], wasps [6], and beetles [7], performing ecosystem services as decomposers, parasitoids, and pest predators, respectively. Carabid beetles are known as natural pest control agents in the soil food web of agroecosystems [8] to feed on invertebrates (aphids, beetles, lepidopterans, slugs, and dipterans) [9] and weed seeds [10]. However, the routine application of agrochemicals causes the exposure of all beneficial insects inhabiting agricultural landscapes, including carabids in soil trophic web, through direct contact with the sprayed products or by eating contaminated food. Currently, laboratory, field, and semi-field studies have been carried out to assess the detrimental effects of pesticides at the organism level [11], and the resulting impact on diversity and abundance of carabid species [12]. Some studies have also supplied evidence that the exposure to herbicides causes mortality or sublethal effects in carabids [13,14]. Metabolomics provide interesting information for the small molecules, and strongly linked to the phenotype [15]. One of the most established protocols in metabolomics is the metabolic profiling biological fluids, which are easy to collect and store, and contain a wide range of molecules that can be strongly affected by physiological disorder including nutrition or disease [16]. Insect haemolymph plasma is a complex mixture of inorganic ions and cations, amino acids, proteins (storage and lipid transport proteins), lipids, carbohydrates, and their degradation products [17]. It is primarily responsible for supplying nutrients and transferring metabolic wastes to maintain homeostasis. The haemolymph plays a central role in the host immune defense processes and provides constant carriage of haemocytes involved in cellular immune responses, i.e., pathogen recognition, phagocytosis, melanization, nodulation, and encapsulation, and humoral effectors including antimicrobial peptides, enzymes that regulates coagulation and melanization, and reactive oxygen (ROS) and nitrogen (RNS) species [18]. Therefore, this fluid represents an opportunity to access to metabolic pool and to collect evidence of metabolic disorders along time-periods [19]. Mass spectrometry profiling (MS) is a technique commonly used to characterize complex matrices and biological fluids [20,21]. Matrix assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) has become increasingly valuable for insect taxonomy [22,23], and also for monitoring significant haemolymph molecular changes during developmental stages [24] or infections [25]. This technique is used to analyze peptides and small proteins located in cells and tissues [26,27] or released into the haemolymph, such as hormones, to regulate diuresis, heartbeat, ecdysis [28], and antimicrobial peptides involved in immune responses [29]. In addition, it is applied to evaluate alterations in protein distribution and expression effects under the agrochemical exposure [30,31,32]. A high resolution of MS data and accurate structure characterization can lead to study a complete set of low molecular weight species in a biological sample drawing a specific fingerprint, time related and ambient associated [33]. This approach has the advantage to highlight both known and unknown metabolites. Identification of metabolites, however, is necessary to draw biological conclusions from untargeted MS data. This step is generally performed by searching the experimental MS/MS data through databases available to the public for free. MS-untargeted profiling of haemolymph is a useful tool to investigate the metabolic response to pendimethalin-based herbicide exposure. P melas, a eurytopic and thermophilus clay soil species inhabiting pastures, open forests, and forest edges in Central and Southern Europe, was selected as a case study. The carabid beetle Pterostichus melas italicus Dejean, 1828 (Coleoptera, Carabidae) is a generalist predator of pests including aphids, lepidopterans, and dipterans [34]. The exposure effect is tested in vivo on P. melas males and females in the reproductive phase of their life cycle to simulate the field exposure in their main period of activity. A recommended field rate (4 L per ha, 330 g/L−1 of active ingredient) is used to evaluate the variability of responses in 21 days, corresponding to the half-life of PND [1,4]. The MS data (all signals) are examined with statistical approaches such as principal component analysis (PCA) for possible correlation with the observed biological effects. Those signals with high correlation are submitted to MS/MS experiments to identify the associated biomarker. Finally, relevant protein–protein interactions network analysis is performed by the STRING [35] database (Search Tool for the Retrieval of Interacting Genes/Proteins).

2. Results and Discussion

For metabolic profiling, haemolymph samples from both males and females of P. melas control and PND-treated were collected at 2, 7, and 21 days after the initial exposure to PND.

2.1. MALDI MS Metabolites Identification

The haemolymph from males and females of P. melas was analyzed by MS, in both control and PND-treated specimens. MS spectra revealed different ion species with different intensity reflecting the individual variability (sex or response to treatment). The identities of some of the detected molecules were assigned through proteomics-generated data sets. A total of twenty-three peptides are identified in a single spectrum (Table 1). Although, the indicated taxonomy species are relatively distant from Pterostichus melas, 12 peptides turned out to be neuropeptides belonging to known families including pyrokinin, FMRFamide, NPF, allatotropin, NP, antidiuretic factor, peptide hormone 1, and peptide tarsal-less AA (Table S1). NeuroPIpred predictive tool was used to recognize the neuropeptides among those identified (https://webs.iiitd.edu.in/raghava/neuropipred/).
The search was performed using SVM threshold: 0.0 and Natural model. A score of 2.7 assigned by the Machine Learning Algorithm allow us to label 12 sequences as putative neuropeptides. Figure 1 shows MS/MS spectrum of pyrokinin-related peptide DGAETPGAAASLWFGPRV-Amide (m/z 1800.94). This signal dominates female chemical profiles in both control and treated samples; conversely, it isn’t detected in control as well treated male samples. Pyrokinins (m/z 1800) are a family of well-studied neuropeptides with myotropic, pheromonotropic, and melanotropic roles in several insects [36], involved in myo-stimulatory activity on visceral muscle [37], including oviduct muscle [38], in pheromone biosynthesis [39] and induction or inhibition of diapause hormone (DH) [40]. FMRFamide (m/z 1201.54) peptides are pleiotropic and thus participate in many processes including the stimulation of various muscles [41,42] and the promotion of the stress-induced sleep [43].
They regulate oogenesis and oviposition in the hemipteran Rhodnius prolixus [44] and heart physiology in Diptera [45]. In Tenebrio molitor and Zophobas atratus, they have myotropic effects in the regulation of contractile activity of the heart, ejaculatory duct, oviduct and the hindgut [46]. Allatotropins (m/z 1367.77) and allatostatins such as allatostatin-3 (ALL3_RHOPR, m/z 1846.97), identified in the haemolymph of P. melas, have inhibitory or stimulatory effects, respectively, on juvenile hormone biosynthesis in the corpora allata [47]. The role of these peptides has also been demonstrated in regulating several other functions including release of digestive enzymes in the midgut [48], myostimulation influencing the physiology of hearts in the dorsal vassel, gut, and reproductive system. The allatostatin effects have been considered on the egg movements in the oviducts and food transit in tenebrionid beetles [49]. Some neuro-derivate factors in insects are involved in the maintenance of homeostasis. There are numerous families of diuretic and anti–diuretic peptides involved in controlling active transport processes of the gut including the creation of urine in the Malpighian tubules and the uptake of water and salts in the hindgut [50]. Antidiuretic factor “ADFa peptide” (ADFA_TENMO, m/z 1542.80) has not been detected outside coleoptera; it is probably the fragment of larger proteins that are not normally liberated and part of regulatory functions [51]. Antimicrobial peptides (AMP) are small proteins having a broad range of activity against bacteria, fungi and viruses as an essential factor for innate immune response [52]. Bombolitin (m/z 1862.19), first isolated from Apoidea (Hymenoptera) venoms, is involved in the mast cell degranulation also showing as bactericide and fungicide properties [53]. Vespulakinin 1 (m/z 1960.11) described in Vespula maculifroms is a vasoactive peptide [54]. Phospholipases A2 (PLA2, m/z 2637.59) have been described in various hymenopteran species and act as neuro- and myotoxins by hydrolyzing membrane phospholipids of motor nerves and subsequently leading to cell lysis [55]. Metalnikowin IIB (m/z 2040.10) is a proline-rich antimicrobial peptide [56], first identified in Palomena prasina [57]. Cecropins (m/z 4083.40) first isolated from the haemolymph of pupae in giant silk moths, Hyalophora cecropia, were recognized as a part of immune response in Lepidoptera, Coleoptera, and Diptera [58]. Cecropins are released by fat bodies and haemocytes in haemolymph and have broad range of antimicrobial activity against bacteria and fungi [59]. Cupiennin-1c (m/z 3770.22) is a cytolytic peptide belonging to a family of antimicrobial peptides isolated from Cupiennius salei, displaying lytic activity towards bacteria, trypanosomes, plasmodia, human blood cells, and cancer cells [60]. Paralytic peptide (m/z 2476.16), first in lepidopteran species, belongs to the cytokine family. It delays larval growth, induces local paralysis around the wound reducing the loss of haemolymph, regulates the immune responses [61], and boosts cell proliferation and morphological variation [62]. Theraphotoxins including U3-theraphotoxin-Hhn1r (m/z 3671.70) and dermonecrotic toxins such as LgSicTox-beta-LOXN4 (m/z 4019.92) have been identified first in spiders to show highly selective inhibitor activity on cell ion channels and cellular lysis, respectively. Serine protease inhibitors (serpins, m/z 3750.72) are functionally distinct, and structurally conserved proteins present in all higher eukaryotes. They are involved in a wide array of physiological functions, including development, innate immune response, reproduction, and host-pathogen interactions [63]. Cell signaling pathways include small peptides such as peptide tarsal-less AA (m/z 3722.89) that coordinate multiple biological processes of cell division, survival, migration and differentiation, giving rise to the embryonic formation of the organs or regulating the innate immune responses [64]. The peptide tarsal-less AA has been studied in Drosophyla to be involved in leg morphogenesis [65], while in the silkworm Bombix mori, it participates in the development of the trachea, silk glands, and Malpighian tubules and is involved in the response to infections of pathogens [66].

2.2. Bioinformatic Analysis

Neuropeptides are responsible for the regulation of several processes, such as homeostasis, development, metabolism, and reproduction, and are produced from larger precursor proteins called prepropeptides [67,68]. Some of the neuropeptides detected in the haemolymph sample perform the duty for multiple functions arranging, as transmitters and co-transmitters with other neuropeptides, several biological process. Protein association network analysis for dynamic/transient interactions of components in signaling and regulatory pathways may be useful to evaluate the hexapoda’s physiological state. Notwithstanding, the large phylogenetic distance across these huge and various classes, neuropeptides display clear similarity among insect taxa [67,68]. In fact, the structure of the genes encoding both neuropeptides and receptors are highly conserved during evolutionary history [69].
Bioinformatic data-analysis is performed using ID of prepropeptide for the analogous reference species Drosophila melanogaster, which is the only other higher insect whose complete genome sequence is available in the public domain. Today at least 50 genes encoding precursors of neuropeptides, peptide hormones, and protein hormones are reported in D. melanogaster. Protein association network analysis for both male and female samples was performed by STRING database (https://string-db.org/). For female control samples, the network is composed of 17 nodes (prepropeptides) and 13 edges (interactions); expected number of edges 1 and PPI enrichment p-value: 3.38 × 10−10. To create this network, a value of 3 for the MCL clustering coefficient was chosen. Two protein clusters are observed, the main protein network showed the set composed by 6 nodes (AstA, Dh44, FMRFa, Hug, Ilp1, sNPF) (Figure 2A) is referred to “Neuropeptide signaling pathway, and insulin-like superfamily”. Accordingly, two networks are observed for male (Figure 2B), and the principal network is composed by 5 nodes (AstA, Dh44, FMRFa, Ilp1, sNPF). The central node in the networks related to both male and female is the neuropeptide sNPF (Figure 2A,B). This feature can’t be explained attributing a unifying global function to sNPF, but underlies that this neuropeptide accomplishes several distinct roles, and it has distributed functions. In fact, Neuropeptide F (sNPF) (m/z 1213) has been suggested to have multiple functional roles in foraging, feeding, alcohol sensitivity, stress, aggression, reproduction, learning function, and circadian rhythm [70,71]. The NPF role in male fertility and female oocyte maturation is agree with the specimens analyzed that were in the reproductive phase of their life cycle.

2.3. Chemical Component Profile of Control and PND-Treated Groups

The chemical component profile of haemolymph specimens (Figure 3) is affected by the individual variability (sex, hormonal variation, and response to treatment). The visual comparative analysis of the spectra, over time after PND exposition, shows changes of the overall profile with consistent modifications in specific regions. Male profiles appear more complex and the feature seems to be related to the time and the treatment. The change overtime of control male profiles is probably related to the insect growth. Otherwise, female profiles (control and treated samples) all seem similar, less complex and dominated by the signal of pyrokinin (m/z 1800). This peptide is absent in male profiles (Figure 3). The detection of pyrokinin in haemolymph of the females may likely related to its myotropic action on the ovipositor contraction during egg laying, as observed in cricket Gryllus bimaculatus [72].
To analyze the relationship between spectra and PND treatment, a pre-processing step is adopted. Nine replicates are aligned and average, in order to obtain a unique representative profile for each specimen used in subsequent chemometric analysis. The data matrix, obtained using the signals of the m/z values selected on the basis of a threshold greater than 400 in terms of signal-to-noise in at least four of the analyzed samples (43 variables), was subjected to principal component analysis (PCA). PCA is a chemometric tool widely used in data exploratory analysis in order to have an overview of data.
PCA transforms the original variables, using an orthogonal linear transformation, to a new set of uncorrelated variables known as principal components (PCs). Representation of the principal component scores and loadings in a bidimensional plot can be used as powerful visualization tool, pointing out patterns hidden in the data set and finding possible correlations between variables. The scores of the samples and the loadings of the variables on the two first principal components are plotted (Figure 4). The information retained is 39.39% of the total variance. The plot of scores obtained by PCA highlights that the female control samples are very similar to each other, while the male control ones are placed in points far from each other, especially for the control sample of the treatment at 7 days (Figure 4A). In any case, all the male control samples are different from female control samples mostly along the first principal components. The difference between these samples can be attributed to the variables placed on the far left and the far right of the loading plot, corresponding to compounds present at higher concentration in male and female samples, respectively (Figure 4B). In particular, the female samples have higher concentrations of peptide tarsal-less AA (m/z 3722.89, Table 1) and pyrokinin (m/z 1800.94). On the other hand, higher concentrations of phospholipase A1 verutoxin-1 (Fragment) (m/z 2637.59), vespulakinin-1 (m/z 1960.11), serine protease inhibitor 2 (m/z 3750.72), bombolitin-3 (m/z 1862.19), M-lycotoxin-Ls4a (m/z 2211.28), dermonecrotic toxin LgSicTox-beta-LOXN4 (m/z 4019.92), and peptide hormone 1 (m/z 2010.94) were found in the male samples (Table 1, Figure 4A,B).
The comparison between the females from the control group and treated samples underlines that there is no significant difference between samples belonging to the group exposed at 7 days and 21 days, whereas a shifting can be observed along PC1 going from control sample to treated one at 2 days from the initial exposure, which indicates an increase of the signal for the m/z values on the left of the plot of loadings. Therefore, the exposure at 2 days involves the concentration increase of phospholipase A1 verutoxin-1 (Fragment) (m/z 2637.59), vespulakinin-1 (m/z 1960.11), serine protease inhibitor 2 (m/z 3750.72), bombolitin-3 (m/z 1862.19), M-lycotoxin-Ls4a (m/z 2211.28), dermonecrotic toxin LgSicTox-beta-LOXN4 (m/z 4019.92), and Peptide hormone 1 (m/z 2010.94), as well as of the unknown compounds placed on the left of the plot. The greatest difference is observed in males at 7 days of exposure compared to the control ones. In particular, the two samples have very different values for both first and second principal component demonstrating that the treatment involves a decrease of the signal for the variables on the left of the plot of loadings (high negative loading values on PC1) and, at the same time, an increase for the variables at the bottom of the loading plot (high negative loading values on PC2, M-lycotoxin-Ls4a of m/z 2211.28, and Cecropin of m/z 4083.40, among the identified compounds). The samples after 21 days from the initial exposure showed a shifting similar but in opposite direction to that observed for the female samples at 2 days of exposure. Finally, a different behavior is detected for the samples belonging to the treatment at 2 days, for which a clear separation occurs exclusively on the second principal component. In this case, the position of the treated sample at higher score values on PC2 means that the treatment determines a decrease of the signal for the m/z values being at the bottom of the loading plot, among which there are the aforementioned peptides M-lycotoxin-Ls4a (m/z 2211.28), and Cecropin (m/z 4083.40).
In order to better evaluate the effect of treatment on the fingerprint profiles obtained by mass spectrometry analysis, another data matrix was constructed using the ratio of the signals of the PND-treated specimens to the signals of the control specimens for each selected m/z value. The application of PCA to this matrix allows to obtain a score plot (Figure 5) in which the variation of the MS profile respect to PND-treatment is clearly shown. Treatment certainly involves different consequences for males and females. In fact, for females the most important variation is substantially along the first principal component, while for males the variations concern only PC2. Therefore, the variables responsible for the modifications of profiles due to treatment are different for males and females. However, a trend common to males and females can be observed. In fact, for both types of samples the greatest variation is noticed passing from the first to the second treatment, while the third treatment are between the first two. This indicates that the involved variables first are subjected to a great variation and then return to assume intermediate values. In particular, for females, a significant decrease of signal ratio passing from the first to the second treatment is observed for the variables placed on the far right of the loading plot shown in Figure 5 (1217.18; 1800.94; 1822.60; 1862.19; 1871.69; 1960.11; 2040.10; 2288.06; 2476.16; 2637.59; 3588.24; 3685.98; 4019.92). On the other hand, for males the variables that show a substantial decrease are the following: 3650.00, 4026.30, 4078.35, and 4233.57. A total of 13 and 4 differentially regulated metabolites were found in female and male, respectively. There are interesting dysregulation patterns in neuropeptides with significant changes observed resulting from PND exposure. However, for males, the treatment at 21 days involves also the increase of signal values for the variables at the bottom of the loadings plot (high negative loading values on PC2), that is 3649.04, 3611.07, 4133.56, and 2822.56. PCA elaboration highlighted differences in the responses of P. melas according to gender.

2.4. Probing Peptide Changes after PND Treatment

Figure 6 A,B depicts the significant expression changes of the identified peptides in haemolymph from females and males of P. melas as a ratio to the control after each exposure duration. These ratios are expressed as a log2 ratio for easier viewing (positive results show upregulation, negative show downregulation). There are interesting and significant peptide changes resulting from PND treatment. There are also interesting general trends observed between the two groups. Some peptides show a constant pattern of up-downregulation over 21 days, others show a hyperarousal pattern. In these cases, the peptide is increased (Metalnikowin-2B, VDKPDYRPRPWPRNMI, and M-lycotoxin-Ls4a, IASHLAFEKLSKLGSKHTML, Figure 6A) at one time point and returns to baseline. A significant downregulation of Paralytic peptide 2 (ENFAGGCTPGYQRTADGRCKPTF, Figure 6B) is observed in male sample after 2 days of exposure suggesting that PND is able to suppress the PP-dependent induction analogously to a pharmacologic p38 MAPK inhibitor [61]. In female samples the downregulation of Paralytic peptide 2 is observed only after 21 days of treatment. In particular, in females of P. melas is observed the upregulation of PP2 in first two days, followed by an important decrease leading to significative downregulation after 21 days of treatment. The temporal shift in upregulation between the groups indicate that distinct pathways are involved in their regulation. The observed trend can be explained considering that PP2 participates to the physiological processes, such as embryonic morphogenesis and larval growth rates. In particular, the observed upregulation of pyrokin and PP2 neuropeptides after 2 days of exposure suggests the activation of signaling pathways to increase muscle contractions in the ovipositor activity and embryonic morphogenesis (Figure 6, panel A).
The downregulation of PSK (Peptide hormone 1, SDLTWTYQSPGDPTNSKN) in males, that become significant only after 21 days, is probably correlated to captivity condition and food copiousness. In fact, the gut peptide PSK acts antagonistically to the hunger signal provided by the adipokinetic hormone (AKH) inhibiting the pTRPγ channel, which is activated under conditions of food shortage [73].
In response to PDN exposure, significant changes occur in males, resulting in downregulation of serpin (EISCEPGTTFQDKCNTCRCGKDGKSAAGCTLKACPQ, m/z 3750.72) at the time points 2 and 7 days, turning to baseline after 21 days. In females, serpin is in baseline at the time points 2 and 7 days, but significantly upregulated after 21 days. Serpins carry out various physiological functions in insects, including development, digestion, host-pathogen interactions, and innate immune response. Serpin plays a key role in immune pathways, its involvement in the regulation of the Toll pathway and prophenoloxidase activation is well documented [74]. Its downregulation was observed in Plutella xylostella after fungal and viral infection, and treatment with mycotoxin [74].
The effect of PND exposure is less important in females than in males, suggesting than females are less sensitive. The observed differences in dysregulation of peptides between males and females highlight that they probably activate different mechanisms of defense for uptake or elimination of toxicants. These data agree with the general mechanism of uptake and elimination of toxicants reported for Arthropods [75]. The females lower their physiological exposure, minimizing any toxic effects of PND in the tissues, by depositing toxicants into their eggs. In males, which do not have this physiological advantage, a significant downregulation of related to immune and defense peptides (M-lycotoxin-Ls4a, Peptide hormone 1, Paralytic peptide 2, and Serine protease inhibitor 2) is observed.

3. Materials and Methods

3.1. Sample Collection and Treatment

Adults of P. melas were collected in an organic olive grove (39°59′27.56′′ N, 16°15′32.64′′ E, 1202 m a.s.l. San Marco Argentano, Calabria, Southern Italy) in October 2019 by using in vivo pitfall traps (plastic jars 9 cm in diameter) containing fruit as an attractant. In the laboratory, beetles were identified by using a dichotomous key, separated by sex and kept in 5 L plastic boxes that were filled to a depth of 6 cm with soil from the capture site. Males and females were separately housed and maintained to be acclimatized for 1 month at 60% relative humidity (rh), under a natural photoperiod, and at room temperature. They were fed mealworms and fruit (organic apples) ad libitum.
A commercial formulation of pendimethalin (PND; Activus EC, product n° HRB00858-39; active ingredient 330 g/L−1) was tested in males and females exposed for 21 days at a recommended field rate (4 L per ha, for cereal and vegetable crops), considering that the half-life of pendimethalin ranges from 24.4 to 34.4 days in acidic sandy soil [1,4]. The experimental design included 6 × 2 plastic boxes (180.5 cm2) filled with the clean sandy soil (pH 5 approximately) from the capture site. Exposure was carried out by spraying the PND solution (7.2 µL of Activus in 14 mL of distilled water) with a pipette on the soil surface of 6 boxes for the treated groups (3 boxes for male and 3 boxes for female, each box containing 10 insects) to simulate the field exposure by contact with the contaminated soil. The control groups, comprising 6 boxes of (3 boxes for male and 3 boxes for female, each box contain 10 insects), were sprayed with distilled water. Males and females, kept separately, were introduced 15 min after the PND solution (or water) was sprayed. Thereafter, 4 boxes (2 control and 2 treated) were used after to 2, 7, and 21 days, respectively. A single application of PND-based herbicide was carried out at 0 day, for treated groups.

3.2. Haemolymph Collection

Males and females from both control (n = 9) and treated (n = 9) groups were randomly chosen at 2, 7 and 21 days after the initial exposure and anaesthetized in a cold chamber at 0 °C for 3 min. Haemolymph was collected by puncturing cold anaesthetized beetles ventrally at the pro-mesothorax articulation with a 29-gauge needle. A pool of 15 μL of haemolymph was collected from three specimens, promptly diluted 1:10 into cold ammonium bicarbonate saline solution (NH4HCO3, 50 mM; Sigma-Aldrich, Darmstadt, Germania) and centrifuged at 1700 rpm for 10 min at 4 °C. Cell-free haemolymph obtained as supernatant was collected and stored at −20 °C before chemical analyses.

3.3. Sample Preparation

An aliquot of each sample (10 µL) was added with 90 µL of solution CHCl3/CH3OH diluted 1:3 (v/v), yielding a lipophilic supernatant fraction and a pellet. The supernatant fractions were analyzed by mass spectrometry. A total of 10 µL of lipophilic supernatant fractions were completely air-dried at room temperature and solubilized in 10 μL of α-cyano-4-hydroxycinnamic acid as MALDI matrix (α-CHCA, 5 mg/mL; H2O/CH3CN, 50:50, v:v; 0.3% TFA). Each sample was spotted onto a 384-well insert Opt-TOF TM stainless steel plate (AB SCIEX, Darmstadt, Germany) and a rapid drying protocol [76,77] was adopted to reduce the inhomogeneous co-crystallization of the analyte with the matrix.

3.4. MALDI MS Analysis

Mass spectrometry analyses (MS) were performed by 5800 MALDI TOF-TOF Analyzer (AB SCIEX, Darmstadt, Germany) equipped with an Nd: YLF Laser with λ = 345 nm wavelength of <500 ps pulse length and to 1000 Hz repetition rate. MS analysis was performed in positive reflectron mode, mass accuracy of 10 ppm and a mass range 900–5000 Da were set to record the untargeted fingerprint profile of low molecular weight species. The peptide mass standards kit (Calibration Mixture 1, AB SCIEX) was used to calibrate the MALDI TOF/TOF mass spectrometer. Spectra with signal-to-noise below 200 were automatically discarded by the instrument and at least 4500 laser shots were typically accumulated with a 400 Hz laser pulse. Each haemolymph sample was spotted three times and three spectra were acquired for each spot in order to obtain a set of nine data point (9 row data). Each spectrum results from the accumulation of 4500 laser shots and the laser was set to continuous random movement to obtain a uniform crystal ablation. The corresponding set of nine data points (9 spectra) was merged and averaged to obtain a unique representative profile for each specimen. All spectra were processed using Data Explorer version 4.11 (AB Sciex) and all nine data points were merged and averaged to obtain a data matrix useful for statistical analysis. A comparative quantification analysis was performed. Collision-induced dissociation (CID)-MS/MS analysis were performed to characterize the selected ion species. MS/MS spectra were performed at a collision energy of 1 kV, and ambient air was used as the collision gas (10−6 Torr) acquiring up to 6000 laser shots e averaging a pulse rate of 1000 Hz.

3.5. Database Proteomics and Targeting Predictions

The MS/MS data were processed to produce pick lists for searching in the SwissProt database using the MASCOT search software (© Matrix Science 2021, http://www.matrixscience.com). The mass tolerance of the parent and fragments for MS/MS data search was set at 20 ppm and 0.20 Da, respectively. During database search the following query were considered “Metazoa (Animals)” taxonomy and choosing “noCleave”. A peaklist of 50 ions of 10% higher intensity than the noise level was generically used for database searching. Methionine oxidation was included in the variable modifications. Although several MS/MS spectra showed intense and well-resolved ion signals, all spectra were manually checked to validate MASCOT results. Neuropeptide homologs are determined by Domain Enhanced Lookup Time Accelerated Basic Local Alignment Search Tool (DELTA-BLAST), limiting the searches to homologs in Drosophila melanogaster. Network analysis was performed by STRING (Search Tool for the Retrieval of Interacting Genes) software (v. 11) (http://stringdb.org/). NeuroPIpred was used to evaluate potential neuropeptides in P. melas (https://webs.iiitd.edu.in/raghava/neuropipred/).

3.6. Statistical Analyses

The MS based untargeted strategy gives an impartial approach on every peak recorded and assures a comprehensive collection of data, without selection criteria. To evaluate possible differences among males and females from both control and PND-treated groups, principal component analysis (PCA) was performed by Statistica 8.0 software package (StatSoft 2007 Edition, Tulsa, OK, USA). The use of PCA allows investigating complex samples and novel metabolic pathways [78]. The analysis included exact mass and IMS clustering. To ensure the quality of the results after data processing, only the signals present in all three biological replicates were considered for differential abundance analysis.
Peptides with significant Student’s t-test (p-value < 0.05, two-tailed) results were considered differentially expressed. Peptides were considered upregulated if the Log2 was greater than 0.5 and downregulated if the Log2 was lower than −0.5.

4. Conclusions

Exposure of P. melas to PND based herbicide altered the metabolic profile of haemolymph. This result indicates that PND interferes with the epigenetic mechanisms of P. melas, and may sabotage the immune system and slow down its development. PCA elaboration highlighted important differences between males and females’ responses after exposure to PND indicating evidence for a sex-specific detoxification response. The metabolomic response highlights the effects of herbicides on non-target animal species, suggesting that their ecological role might be compromised.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27144645/s1, Table S1: Predicted neuropeptides by NeuroPIpred predictive tool (https://webs.iiitd.edu.in/raghava/neuropipred/).

Author Contributions

Conceptualization, F.T., D.A., A.G. and A.N.; methodology, F.T., M.L.V., D.A., A.G., A.T. and A.N.; software, D.A., A.T. and A.N.; validation, D.A., A.T. and A.N.; formal analysis, D.A., A.T. and A.N.; investigation, D.A., A.T. and A.N.; data curation, D.A., A.T. and A.N.; writing—original draft preparation, A.G., D.A., A.T. and A.N.; writing—review and editing, A.G., A.T. and A.N., visualization, A.T. and A.N.; supervision, A.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

There are no conflicts of interest to declare.

Sample Availability

Samples are available from the authors.

References

  1. Vighi, M.; Matthies, M.; Solomon, K.R. Critical Assessment of Pendimethalin in Terms of Persistence, Bioaccumulation, Toxicity, and Potential for Long-Range Transport. J. Toxicol. Environ. Health Part B 2017, 20, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Strandberg, M.; Scott-Fordsmand, J.J. Effects of Pendimethalin at Lower Trophic Levels—A Review. Ecotoxicol. Environ. Saf. 2004, 57, 190–201. [Google Scholar] [CrossRef] [PubMed]
  3. Roca, E.; D’Errico, E.; Izzo, A.; Strumia, S.; Esposito, A.; Fiorentino, A. In Vitro Saprotrophic Basidiomycetes Tolerance to Pendimethalin. Int. Biodeterior. Biodegrad. 2009, 63, 182–186. [Google Scholar] [CrossRef]
  4. Kočárek, M.; Artikov, H.; Voříšek, K.; Borůvka, L. Pendimethalin Degradation in Soil and Its Interaction with Soil Microorganisms. Soil Water Res. 2016, 11, 213–219. [Google Scholar] [CrossRef] [Green Version]
  5. Belden, J.B.; Phillips, T.A.; Clark, B.W.; Coats, J.R. Toxicity of Pendimethalin to Nontarget Soil Organisms. Bull. Environ. Contam. Toxicol. 2005, 74, 769–776. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Oliver, J.B.; Reding, M.E.; Moyseenko, J.J.; Klein, M.G.; Mannion, C.M.; Bishop, B. Survival of Adult Tiphia vernalis (Hymenoptera: Tiphiidae) After Insecticide, Fungicide, and Herbicide Exposure in Laboratory Bioassays. J. Econ. Entomol. 2006, 99, 288–294. [Google Scholar] [CrossRef]
  7. Giglio, A.; Cavaliere, F.; Giulianini, P.G.; Kurtz, J.; Vommaro, M.L.; Brandmayr, P. Continuous Agrochemical Treatments in Agroecosystems Can Modify the Effects of Pendimethalin-Based Herbicide Exposure on Immunocompetence of a Beneficial Ground Beetle. Diversity 2019, 11, 241. [Google Scholar] [CrossRef] [Green Version]
  8. De Heij, S.E.; Willenborg, C.J. Connected Carabids: Network Interactions and Their Impact on Biocontrol by Carabid Beetles. BioScience 2020, 70, 490–500. [Google Scholar] [CrossRef]
  9. Giglio, A.; Brandmayr, P.; Talarico, F.; Giulianini, P.G. Effects of alternative and specialised diet on development and survival of larvae and pupae in Carabus (Chaetocarabus) lefebvrei (Coleoptera: Carabidae). Entomol. Gen. 2012, 33, 263–271. [Google Scholar] [CrossRef]
  10. Foffová, H.; Ćavar Zeljković, S.; Honěk, A.; Martinková, Z.; Tarkowski, P.; Saska, P. Which seed properties determine the preferences of carabid beetle seed predators? Insects 2020, 11, 757. [Google Scholar] [CrossRef]
  11. Giglio, A.; Cavaliere, F.; Giulianini, P.G.; Mazzei, A.; Talarico, F.; Vommaro, M.L.; Brandmayr, P. Impact of Agrochemicals on Non-Target Species: Calathus Fuscipes Goeze 1777 (Coleoptera: Carabidae) as Model. Ecotoxicol. Environ. Saf. 2017, 142, 522–529. [Google Scholar] [CrossRef] [PubMed]
  12. Holland, J.M.; Luff, M.L. The Effects of Agricultural Practices on Carabidae in Temperate Agroecosystems. Integr. Pest Manag. Rev. 2000, 5, 109–129. [Google Scholar] [CrossRef]
  13. Cavaliere, F.; Brandmayr, P.; Giglio, A. DNA Damage in Haemocytes of Harpalus (Pseudophonus) Rufipes (De Geer, 1774) (Coleoptera, Carabidae) as an Indicator of Sublethal Effects of Exposure to Herbicides. Ecol. Indic. 2019, 98, 88–91. [Google Scholar] [CrossRef]
  14. Prosser, R.S.; Anderson, J.C.; Hanson, M.L.; Solomon, K.R.; Sibley, P.K. Indirect Effects of Herbicides on Biota in Terrestrial Edge-of-Field Habitats: A Critical Review of the Literature. Agric. Ecosyst. Environ. 2016, 232, 59–72. [Google Scholar] [CrossRef]
  15. Alonso, A.; Marsal, S.; Julià, A. Analytical Methods in Untargeted Metabolomics: State of the Art in 2015. Front. Bioeng. Biotechnol. 2015, 3, 23. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, X.; Li, Q.; Xu, Z.; Dou, J. Mass Spectrometry-Based Metabolomics in Health and Medical Science: A Systematic Review. RSC Adv. 2020, 10, 3092–3104. [Google Scholar] [CrossRef] [Green Version]
  17. Chapman, R.F.; Simpson, S.J.; Douglas, A.E. The Insects: Structure and Function, 5th ed.; Cambridge University Press: New York, NY, USA, 2013. [Google Scholar]
  18. Rosales, C. Cellular and Molecular Mechanisms of Insect Immunity. In Insect Physiology and Ecology; Shields, V.D.C., Ed.; InTech: Rijeka, Croatia, 2017. [Google Scholar] [CrossRef] [Green Version]
  19. Aliferis, K.A.; Copley, T.; Jabaji, S. Gas Chromatography–Mass Spectrometry Metabolite Profiling of Worker Honey Bee (Apis mellifera L.) Hemolymph for the Study of Nosema Ceranae Infection. J. Insect Physiol. 2012, 58, 1349–1359. [Google Scholar] [CrossRef]
  20. Aiello, D.; Giambona, A.; Leto, F.; Passarello, C.; Damiani, G.; Maggio, A.; Siciliano, C.; Napoli, A. Human Coelomic Fluid Investigation: A MS-Based Analytical Approach to Prenatal Screening. Sci. Rep. 2018, 8, 10973. [Google Scholar] [CrossRef] [Green Version]
  21. Imbrogno, S.; Aiello, D.; Filice, M.; Leo, S.; Mazza, R.; Cerra, M.C.; Napoli, A. MS-Based Proteomic Analysis of Cardiac Response to Hypoxia in the Goldfish (Carassius auratus). Sci. Rep. 2019, 9, 18953. [Google Scholar] [CrossRef] [Green Version]
  22. Chavy, A.; Nabet, C.; Normand, A.C.; Kocher, A.; Ginouves, M.; Prévot, G.; Vasconcelos dos Santos, T.; Demar, M.; Piarroux, R.; de Thoisy, B. Identification of French Guiana Sand Flies Using MALDI-TOF Mass Spectrometry with a New Mass Spectra Library. PLoS Negl. Trop. Dis. 2019, 13, e0007031. [Google Scholar] [CrossRef]
  23. Vega-Rúa, A.; Pagès, N.; Fontaine, A.; Nuccio, C.; Hery, L.; Goindin, D.; Gustave, J.; Almeras, L. Improvement of Mosquito Identification by MALDI-TOF MS Biotyping Using Protein Signatures from Two Body Parts. Parasites Vectors 2018, 11, 574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hou, Y.; Zou, Y.; Wang, F.; Gong, J.; Zhong, X.; Xia, Q.; Zhao, P. Comparative Analysis of Proteome Maps of Silkworm Hemolymph during Different Developmental Stages. Proteome Sci. 2010, 8, 45. [Google Scholar] [CrossRef] [Green Version]
  25. Arafah, K.; Voisin, S.N.; Masson, V.; Alaux, C.; Le Conte, Y.; Bocquet, M.; Bulet, P. MALDI–MS Profiling to Address Honey Bee Health Status under Bacterial Challenge through Computational Modeling. Proteomics 2019, 19, 1900268. [Google Scholar] [CrossRef]
  26. Dani, F.R.; Francese, S.; Mastrobuoni, G.; Felicioli, A.; Caputo, B.; Simard, F.; Pieraccini, G.; Moneti, G.; Coluzzi, M.; della Torre, A.; et al. Exploring Proteins in Anopheles Gambiae Male and Female Antennae through MALDI Mass Spectrometry Profiling. PLoS ONE 2008, 3, e2822. [Google Scholar] [CrossRef] [PubMed]
  27. Verhaert, P.D.; Conaway, M.C.P.; Pekar, T.M.; Miller, K. Neuropeptide Imaging on an LTQ with VMALDI Source: The Complete ‘All-in-One’ Peptidome Analysis. Int. J. Mass Spectrom. 2007, 260, 177–184. [Google Scholar] [CrossRef]
  28. Wegener, C.; Reinl, T.; Jänsch, L.; Predel, R. Direct Mass Spectrometric Peptide Profiling and Fragmentation of Larval Peptide Hormone Release Sites in Drosophila Melanogaster Reveals Tagma-Specific Peptide Expression and Differential Processing: Peptidomics of Hormone Release Sites in Drosophila. J. Neurochem. 2006, 96, 1362–1374. [Google Scholar] [CrossRef] [PubMed]
  29. Chowdhury, T.; Mandal, S.M.; Dutta, S.; Ghosh, A.K. Identification of a Novel Proline-rich Antimicrobial Protein from the Hemolymph of Antheraea Mylitta. Arch. Insect Biochem. Physiol. 2021, 106, e21771. [Google Scholar] [CrossRef]
  30. Catae, A.F.; da Silva Menegasso, A.R.; Pratavieira, M.; Palma, M.S.; Malaspina, O.; Roat, T.C. MALDI-Imaging Analyses of Honeybee Brains Exposed to a Neonicotinoid Insecticide: MALDI-Imaging of Honeybee Brains. Pest Manag. Sci. 2019, 75, 607–615. [Google Scholar] [CrossRef]
  31. Roat, T.C.; dos Santos-Pinto, J.R.A.; dos Santos, L.D.; Santos, K.S.; Malaspina, O.; Palma, M.S. Modification of the Brain Proteome of Africanized Honeybees (Apis mellifera) Exposed to a Sub-lethal Doses of the Insecticide Fipronil. Ecotoxicology 2014, 23, 1659–1670. [Google Scholar] [CrossRef]
  32. Thornton, B.J.; Elthon, T.E.; Cerny, R.L.; Siegfried, B.D. Proteomic Analysis of Atrazine Exposure in Drosophila melanogaster (Diptera: Drosophilidae). Chemosphere 2010, 81, 235–241. [Google Scholar] [CrossRef]
  33. Aiello, D.; Lucà, F.; Siciliano, C.; Frati, P.; Fineschi, V.; Rongo, R.; Napoli, A. Analytical Strategy for MS-Based Thanatochemistry to Estimate Postmortem Interval. J. Proteome Res. 2021, 20, 2607–2617. [Google Scholar] [CrossRef] [PubMed]
  34. Sunderland, K.D. Invertebrate Pest Control by Carabids. In The Agroecology of Carabid Beetles; Holland, J.M., Ed.; Intercept Limited: Andover, UK, 2002; pp. 165–214. [Google Scholar]
  35. Szklarczyk, D.; Franceschini, A.; Kuhn, M.; Simonovic, M.; Roth, A.; Minguez, P.; Doerks, T.; Stark, M.; Muller, J.; Bork, P.; et al. The STRING Database in 2011: Functional Interaction Networks of Proteins, Globally Integrated and Scored. Nucleic Acids Res. 2011, 39, D561–D568. [Google Scholar] [CrossRef] [PubMed]
  36. Shiomi, K. Pyrokinin. In Handbook of Hormones; Elsevier: Amsterdam, The Netherlands, 2021; pp. 729–730. [Google Scholar] [CrossRef]
  37. Lajevardi, A.; Paluzzi, J.-P.V. Receptor Characterization and Functional Activity of Pyrokinins on the Hindgut in the Adult Mosquito, Aedes Aegypti. Front. Physiol. 2020, 11, 490. [Google Scholar] [CrossRef]
  38. Holman, G.M.; Cook, B.J.; Nachman, R.J. Primary Structure and Synthesis of a Blocked Myotropic Neuropeptide Isolated from the Cockroach, Leucophaea Maderae. Comp. Biochem. Physiol. Part C Comp. Pharmacol. 1986, 85, 219–224. [Google Scholar] [CrossRef]
  39. Raina, A.K.; Menn, J.J. Pheromone Biosynthesis Activating Neuropeptide: From Discovery to Current Status. Arch. Insect Biochem. Physiol. 1993, 22, 141–151. [Google Scholar] [CrossRef]
  40. Zhang, Q.; Nachman, R.J.; Kaczmarek, K.; Zabrocki, J.; Denlinger, D.L. Disruption of Insect Diapause Using Agonists and an Antagonist of Diapause Hormone. Proc. Natl. Acad. Sci. USA 2011, 108, 16922–16926. [Google Scholar] [CrossRef] [Green Version]
  41. Sedra, L.; Lange, A.B. The Female Reproductive System of the Kissing Bug, Rhodnius prolixus: Arrangements of Muscles, Distribution and Myoactivity of Two Endogenous FMRFamide-like Peptides. Peptides 2014, 53, 140–147. [Google Scholar] [CrossRef] [PubMed]
  42. Suggs, J.M.; Jones, T.H.; Murphree, C.S.; Hillyer, J.F. CCAP and FMRFamide-like Peptides Accelerate the Contraction Rate of the Antennal Accessory Pulsatile Organs (Auxiliary Hearts) of Mosquitoes. J. Exp. Biol. 2016, 219, 2388–2395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Lenz, O.; Xiong, J.; Nelson, M.D.; Raizen, D.M.; Williams, J.A. FMRFamide Signaling Promotes Stress-Induced Sleep in Drosophila. Brain Behav. Immun. 2015, 47, 141–148. [Google Scholar] [CrossRef] [Green Version]
  44. Sedra, L.; Lange, A.B. Cloning and Expression of Long Neuropeptide F and the Role of FMRFamide-like Peptides in Regulating Egg Production in the Chagas Vector, Rhodnius prolixus. Peptides 2016, 82, 1–11. [Google Scholar] [CrossRef]
  45. Hillyer, J.F. Insect Heart Rhythmicity Is Modulated by Evolutionarily Conserved Neuropeptides and Neurotransmitters. Curr. Opin. Insect Sci. 2018, 29, 41–48. [Google Scholar] [CrossRef] [PubMed]
  46. Marciniak, P.; Witek, W.; Szymczak, M.; Pacholska-Bogalska, J.; Chowański, S.; Kuczer, M.; Rosiński, G. FMRFamide-Related Peptides Signaling Is Involved in the Regulation of Muscle Contractions in Two Tenebrionid Beetles. Front. Physiol. 2020, 11, 456. [Google Scholar] [CrossRef] [PubMed]
  47. Abou El Asrar, R.; Cools, D.; Vanden Broeck, J. Role of Peptide Hormones in Insect Gut Physiology. Curr. Opin. Insect Sci. 2020, 41, 71–78. [Google Scholar] [CrossRef]
  48. Rudwall, A.J.; Sliwowska, J.; Nassel, D.R. Allatotropin-like Neuropeptide in the Cockroach Abdominal Nervous System: Myotropic Actions, Sexually Dimorphic Distribution and Colocalization with Serotonin. J. Comp. Neurol. 2000, 428, 159–173. [Google Scholar] [CrossRef]
  49. Lubawy, J.; Marciniak, P.; Kuczer, M.; Rosiński, G. Myotropic Activity of Allatostatins in Tenebrionid Beetles. Neuropeptides 2018, 70, 26–36. [Google Scholar] [CrossRef] [PubMed]
  50. Lubawy, J.; Urbański, A.; Colinet, H.; Pflüger, H.-J.; Marciniak, P. Role of the Insect Neuroendocrine System in the Response to Cold Stress. Front. Physiol. 2020, 11, 376. [Google Scholar] [CrossRef]
  51. Li, B.; Predel, R.; Neupert, S.; Hauser, F.; Tanaka, Y.; Cazzamali, G.; Williamson, M.; Arakane, Y.; Verleyen, P.; Schoofs, L.; et al. Genomics, Transcriptomics, and Peptidomics of Neuropeptides and Protein Hormones in the Red Flour Beetle Tribolium castaneum. Genome Res. 2008, 18, 113–122. [Google Scholar] [CrossRef] [Green Version]
  52. Wu, Q.; Patočka, J.; Kuča, K. Insect Antimicrobial Peptides, a Mini Review. Toxins 2018, 10, 461. [Google Scholar] [CrossRef]
  53. Stöcklin, R.; Favreau, P.; Thai, R.; Pflugfelder, J.; Bulet, P.; Mebs, D. Structural Identification by Mass Spectrometry of a Novel Antimicrobial Peptide from the Venom of the Solitary Bee Osmia rufa (Hymenoptera: Megachilidae). Toxicon 2010, 55, 20–27. [Google Scholar] [CrossRef]
  54. Grosch, J.; Hilger, C.; Bilò, M.B.; Kler, S.; Schiener, M.; Dittmar, G.; Bernardin, F.; Lesur, A.; Ollert, M.; Schmidt-Weber, C.B.; et al. Shedding Light on the Venom Proteomes of the Allergy-Relevant Hymenoptera Polistes dominula (European Paper Wasp) and Vespula Spp. (Yellow Jacket). Toxins 2020, 12, 323. [Google Scholar] [CrossRef]
  55. Perez-Riverol, A.; Lasa, A.M.; dos Santos-Pinto, J.R.A.; Palma, M.S. Insect Venom Phospholipases A1 and A2: Roles in the Envenoming Process and Allergy. Insect Biochem. Mol. Biol. 2019, 105, 10–24. [Google Scholar] [CrossRef] [PubMed]
  56. Ebbensgaard, A.; Mordhorst, H.; Overgaard, M.T.; Nielsen, C.G.; Aarestrup, F.M.; Hansen, E.B. Comparative Evaluation of the Antimicrobial Activity of Different Antimicrobial Peptides against a Range of Pathogenic Bacteria. PLoS ONE 2015, 10, e0144611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Chernysh, S.; Cociancich, S.; Briand, J.-P.; Hetru, C.; Bulet, P. The Inducible Antibacterial Peptides of the Hemipteran Insect Palomena Prasina: Identification of a Unique Family of Prolinerich Peptides and of a Novel Insect Defensin. J. Insect Physiol. 1996, 42, 81–89. [Google Scholar] [CrossRef]
  58. Brady, D.; Grapputo, A.; Romoli, O.; Sandrelli, F. Insect Cecropins, Antimicrobial Peptides with Potential Therapeutic Applications. Int. J. Mol. Sci. 2019, 20, 5862. [Google Scholar] [CrossRef] [Green Version]
  59. Yakovlev, A.Y.; Nesin, A.P.; Simonenko, N.P.; Gordya, N.A.; Tulin, D.V.; Kruglikova, A.A.; Chernysh, S.I. Fat Body and Hemocyte Contribution to the Antimicrobial Peptide Synthesis in Calliphora vicina R.-D. (Diptera: Calliphoridae) Larvae. In Vitro Cell. Dev. Biol. Anim. 2017, 53, 33–42. [Google Scholar] [CrossRef] [PubMed]
  60. Kuhn-Nentwig, L.; Willems, J.; Seebeck, T.; Shalaby, T.; Kaiser, M.; Nentwig, W. Cupiennin 1a Exhibits a Remarkably Broad, Non-Stereospecific Cytolytic Activity on Bacteria, Protozoan Parasites, Insects, and Human Cancer Cells. Amino Acid 2011, 40, 69–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Ishii, K.; Hamamoto, H.; Kamimura, M.; Nakamura, Y.; Noda, H.; Imamura, K.; Mita, K.; Sekimizu, K. Insect Cytokine Paralytic Peptide (PP) Induces Cellular and Humoral Immune Responses in the Silkworm Bombyx Mori. J. Biol. Chem. 2010, 285, 28635–28642. [Google Scholar] [CrossRef] [Green Version]
  62. Song, L.; Wang, F.; Dong, Z.; Hua, X.; Xia, Q. Label-Free Quantitative Phosphoproteomic Profiling of Cellular Response Induced by an Insect Cytokine Paralytic Peptide. J. Proteom. 2017, 154, 49–58. [Google Scholar] [CrossRef]
  63. Shakeel, M.; Xu, X.; Mandal, S.; Jin, F. Role of Serine Protease Inhibitors in Insect-host-pathogen Interactions. Arch. Insect Biochem. Physiol. 2019, 102, e21556. [Google Scholar] [CrossRef]
  64. Zhu, M.; Hu, X.; Cao, G.; Xue, R.; Gong, C. Functions and Impact of Tal-like Genes in Animals with Regard to Applied Aspects. Appl. Microbiol. Biotechnol. 2018, 102, 6841–6845. [Google Scholar] [CrossRef]
  65. Pueyo, J.I.; Couso, J.P. Tarsal-Less Peptides Control Notch Signalling through the Shavenbaby Transcription Factor. Dev. Biol. 2011, 355, 183–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Cao, G.; Gong, Y.; Hu, X.; Zhu, M.; Liang, Z.; Huang, L.; Yu, L.; Xu, J.; Li, K.; Zar, M.S.; et al. Identification of Tarsal-Less Peptides from the Silkworm Bombyx Mori. Appl. Microbiol. Biotechnol. 2018, 102, 1809–1822. [Google Scholar] [CrossRef] [PubMed]
  67. Kim, Y.; Bark, S.; Hook, V.; Bandeira, N. NeuroPedia: Neuropeptide Database and Spectral Library. Bioinformatics 2011, 27, 2772–2773. [Google Scholar] [CrossRef] [Green Version]
  68. Yeoh, J.G.C.; Pandit, A.A.; Zandawala, M.; Nässel, D.R.; Davies, S.-A.; Dow, J.A.T. DINeR: Database for Insect Neuropeptide Research. Insect Biochem. Mol. Biol. 2017, 86, 9–19. [Google Scholar] [CrossRef]
  69. Veenstra, J.A. Two Lys-Vasopressin-like Peptides, EFLamide, and Other Phasmid Neuropeptides. Gen. Comp. Endocrinol. 2019, 278, 3–11. [Google Scholar] [CrossRef] [PubMed]
  70. Cui, H.; Zhao, Z. Structure and Function of Neuropeptide F in Insects. J. Integr. Agric. 2020, 19, 1429–1438. [Google Scholar] [CrossRef]
  71. Nässel, D.R.; Wegener, C. A Comparative Review of Short and Long Neuropeptide F Signaling in Invertebrates: Any Similarities to Vertebrate Neuropeptide Y Signaling? Peptides 2011, 32, 1335–1355. [Google Scholar] [CrossRef]
  72. Abdoun, K.; Mesnier-Sabin, M.; Baudry-Partiaoglou, N.; Nicolas, P.; Cohen, P. Separation of Oviposition-Stimulating Peptides and Myotropic Factors from Head Extracts of Galleria mellonella L.: Comparative Effects of Myotropic and Non-Myotropic Factors on Egg Laying. J. Comp. Physiol. B 1995, 165, 102–109. [Google Scholar] [CrossRef]
  73. Wicher, D.; Derst, C.; Gautier, H.; Lapied, B.; Heinemann, S.H.; Agricola, H.J. The satiety signaling neuropeptide perisulfakinin inhibits the activity of central neurons promoting general activity. Front. Cell. Neurosci. 2007, 1, 3. [Google Scholar] [CrossRef] [Green Version]
  74. Chen, J.; Cui, D.; Ullah, H.; Hao, K.; Tu, X.; Zhang, Z. Serpin7 Controls Egg Diapause of Migratory Locust (Locusta migratoria) by Regulating Polyphenol Oxidase. FEBS Open Bio 2020, 10, 707–717. [Google Scholar] [CrossRef] [Green Version]
  75. McClellan-Green, P.; Romano, J.; Oberdörster, E. Does Gender Really Matter in Contaminant Exposure? A Case Study Using Invertebrate Models. Environ. Res. 2007, 104, 183–191. [Google Scholar] [CrossRef] [PubMed]
  76. Aiello, D.; Siciliano, C.; Mazzotti, F.; Di Donna, L.; Athanassopoulos, C.M.; Napoli, A. A Rapid MALDI MS/MS Based Method for Assessing Saffron (Crocus sativus L.) Adulteration. Food Chem. 2020, 307, 125527. [Google Scholar] [CrossRef] [PubMed]
  77. Napoli, A.; Athanassopoulos, C.M.; Moschidis, P.; Aiello, D.; Di Donna, L.; Mazzotti, F.; Sindona, G. Solid Phase Isobaric Mass Tag Reagent for Simultaneous Protein Identification and Assay. Anal. Chem. 2010, 82, 5552–5560. [Google Scholar] [CrossRef] [PubMed]
  78. Naccarato, A.; Furia, E.; Sindona, G.; Tagarelli, A. Multivariate Class Modeling Techniques Applied to Multielement Analysis for the Verification of the Geographical Origin of Chili Pepper. Food Chem. 2016, 206, 217–222. [Google Scholar] [CrossRef] [PubMed]
Figure 1. MALDI MS/MS of the m/z 1800.94.
Figure 1. MALDI MS/MS of the m/z 1800.94.
Molecules 27 04645 g001
Figure 2. Protein–protein interaction network by STRING software version 11, https://string-db.org/ (panel (A), female; panel (B), male).
Figure 2. Protein–protein interaction network by STRING software version 11, https://string-db.org/ (panel (A), female; panel (B), male).
Molecules 27 04645 g002
Figure 3. Untargeted molecular profiles of control (C) and treated (T) male insects (panel (A)) and female insects (panel (B)), after 2 (I), 7 (II), and 21 (III) days of PND treatment.
Figure 3. Untargeted molecular profiles of control (C) and treated (T) male insects (panel (A)) and female insects (panel (B)), after 2 (I), 7 (II), and 21 (III) days of PND treatment.
Molecules 27 04645 g003
Figure 4. Scores (A) and loadings (B) plots obtained by Principal Component Analysis applied to mass spectrometry dataset.
Figure 4. Scores (A) and loadings (B) plots obtained by Principal Component Analysis applied to mass spectrometry dataset.
Molecules 27 04645 g004
Figure 5. Scores (A) and loadings (B) plots obtained by Principal Component Analysis applied to the ratio of the signals of treated specimens to the signals of the control specimens.
Figure 5. Scores (A) and loadings (B) plots obtained by Principal Component Analysis applied to the ratio of the signals of treated specimens to the signals of the control specimens.
Molecules 27 04645 g005
Figure 6. The ratio for the intensity after 2, 7, and 21 days of PND exposure relative to the control is depicted on the y-axis as a log2 ratio to clearly see up/downregulation of peptides relative to the x-axis. Panel (A,B) show up/downregulation of peptides for female and male groups, respectively. The error bars denote the standard error of the mean (n = 3).
Figure 6. The ratio for the intensity after 2, 7, and 21 days of PND exposure relative to the control is depicted on the y-axis as a log2 ratio to clearly see up/downregulation of peptides relative to the x-axis. Panel (A,B) show up/downregulation of peptides for female and male groups, respectively. The error bars denote the standard error of the mean (n = 3).
Molecules 27 04645 g006
Table 1. Identified proteins by MALDI MS/MS and Mascot software by Matrix Science (www.matrixscience.com). The MS/MS data were processed using a mass tolerance of 20 ppm and 0.2 Da for the precursor and fragment ions, respectively. a According to “UniProtKB” (http://www.uniprot.org/).
Table 1. Identified proteins by MALDI MS/MS and Mascot software by Matrix Science (www.matrixscience.com). The MS/MS data were processed using a mass tolerance of 20 ppm and 0.2 Da for the precursor and fragment ions, respectively. a According to “UniProtKB” (http://www.uniprot.org/).
IDaName aSequenceExactm/z
1.FAR14_SARBUFMRFamide-14DPHHDFMRF1201.52131201.54
2.NPF1_LEPDENeuropeptide NPF-1ARGPQLRLRF1213.72821213.75
3.ALLTR_ACRHIAllatotropin-related peptideGFKNVALSTARGF1367.74351367.77
ALLTR_BANDIAllatotropin-related peptideGFKNVALSTARGF
ALLTR_EUSSEAllatotropin-related peptideGFKNVALSTARGF
ALLTR_NEZVIAllatotropin-related peptideGFKNVALSTARGF
ALLTR_ONCFAAllatotropin-related peptideGFKNVALSTARGF
ALLTR_PENRUAllatotropin-related peptideGFKNVALSTARGF
ALLTR_PERAMAllatotropin-related peptideGFKNVALSTARGF
ALLTR_PYRAPAllatotropin-related peptideGFKNVALSTARGF
4.ADFA_TENMOAntidiuretic factorVVNTPGHAVSYHVY1542.77051542.80
5.TXS6D_CUPSAShort cationic peptide-6dINKYREWKNKKN1620.89741620.92
6.PPK_SCHGRPyrokininDGAETPGAAASLWFGPRV—Amide1800.90321800.94
7.ALL3_RHOPRAllatostatin-3QVSLKYPEGKMYSFGL1846.94131846.97
8.BOL3_BOMPEBombolitin-3IKIMDILAKLGKVLAHV1862.16651862.19
9.BRK_VESMCVespulakinin-1TATTRRRGRPPGFSPFR1960.07411960.11
10.LYC1_LYCSIM-lycotoxin-Ls3aGKLQAFLAKMKEIAAQTL1961.12571961.16
11.PH1_PERAMPeptide hormone 1SDLTWTYQSPGDPTNSKN2010.90452010.94
12.MK2B_PALPRMetalnikowin-2BVDKPDYRPRPWPRNMI2040.06012040.10
13.LYC40_LYCSIM-lycotoxin-Ls4aIASHLAFEKLSKLGSKHTML2211.23232211.28
14.PAP2_SPOEXParalytic peptide 2ENFAGGCTPGYQRTADGRCKPTF2476.11382476.16
15.PA11_VESVEPhospholipase A1 verutoxin-1 (Fragment)GLLPKVKLVPEQISFILSTRENR2637.54552637.59
16.TXC6D_CUPSACupiennin-6dFINTIKLLIEKYREWKNKQSS2638.47212638.52
17.HN423_CYRHAU3-theraphotoxin-Hhn1rDCAGYMRECKEKLCCSGYVCSSRWKWCVLPAP3671.62223671.70
18.MSPI2_MELSASerine protease inhibitor 2EISCEPGTTFQDKCNTCRCGKDGKSAAGCTLKACPQ3750.64763750.72
19.TXC1C_CUPSACupiennin-1cGFGSLFKFLAKKVAKTVAKQAAKQGAKYIANKQTE3770.14843770.22
20.TALAA_DROMEPeptide tarsal-less AALDPTGTYRRPRDTQDSRQKRRQDCLDPTGQY3722.81693722.89
21.BX4_LOXGADermonecrotic toxin LgSicTox-beta-LOXN4ADSRKPDDRYDMSGNDALGDVKLATYEDNPWETFK4019.83574019.92
22.CEC_CALVICecropinGWLKKIGKKIGRVGQHTRDATIQGLAVAQQAANVAATAR4083.31674083.40
23.DIUH1_TENMODiuretic hormone 1SPTISITAPIDVLRKTWEQERARKQMVKNREFLNSLN4369.35664369.44
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aiello, D.; Giglio, A.; Talarico, F.; Vommaro, M.L.; Tagarelli, A.; Napoli, A. Mass Spectrometry-Based Peptide Profiling of Haemolymph from Pterostichus melas Exposed to Pendimethalin Herbicide. Molecules 2022, 27, 4645. https://doi.org/10.3390/molecules27144645

AMA Style

Aiello D, Giglio A, Talarico F, Vommaro ML, Tagarelli A, Napoli A. Mass Spectrometry-Based Peptide Profiling of Haemolymph from Pterostichus melas Exposed to Pendimethalin Herbicide. Molecules. 2022; 27(14):4645. https://doi.org/10.3390/molecules27144645

Chicago/Turabian Style

Aiello, Donatella, Anita Giglio, Federica Talarico, Maria Luigia Vommaro, Antonio Tagarelli, and Anna Napoli. 2022. "Mass Spectrometry-Based Peptide Profiling of Haemolymph from Pterostichus melas Exposed to Pendimethalin Herbicide" Molecules 27, no. 14: 4645. https://doi.org/10.3390/molecules27144645

APA Style

Aiello, D., Giglio, A., Talarico, F., Vommaro, M. L., Tagarelli, A., & Napoli, A. (2022). Mass Spectrometry-Based Peptide Profiling of Haemolymph from Pterostichus melas Exposed to Pendimethalin Herbicide. Molecules, 27(14), 4645. https://doi.org/10.3390/molecules27144645

Article Metrics

Back to TopTop