A rapid MALDI-TOF mass spectrometry workflow for Drosophila melanogaster differential neuropeptidomics
Molecular Brain volume 6, Article number: 60 (2013)
Neuropeptides are a diverse category of signaling molecules in the nervous system regulating a variety of processes including food intake, social behavior, circadian rhythms, learning, and memory. Both the identification and functional characterization of specific neuropeptides are ongoing fields of research. Matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) analysis of nervous tissues from a variety of organisms allows direct detection and identification of neuropeptides. Here, we demonstrate an analysis workflow that allows for the detection of differences in specific neuropeptides amongst a variety of neuropeptides being simultaneously measured. For sample preparation, we describe a straight-forward and rapid (minutes) method where individual adult Drosophila melanogaster brains are analyzed. Using a MATLAB-based data analysis workflow, also compatible with MALDI-TOF mass spectra obtained from other sample preparations and instrumentation, we demonstrate how changes in neuropeptides levels can be detected with this method.
Over fifty isotopically resolved ion signals in the peptide mass range are reproducibly observed across experiments. MALDI-TOF MS profile spectra were used to statistically identify distinct relative differences in organ-wide endogenous levels of detected neuropeptides between biological conditions. In particular, three distinct levels of a particular neuropeptide, pigment dispersing factor, were detected by comparing groups of preprocessed spectra obtained from individual brains across three different D. melanogaster strains, each of which express different amounts of this neuropeptide. Using the same sample preparation, MALDI-TOF/TOF tandem mass spectrometry confirmed that at least 14 ion signals observed across experiments are indeed neuropeptides. Among the identified neuropeptides were three products of the neuropeptide-like precursor 1 gene previously not identified in the literature.
Using MALDI-TOF MS and preprocessing/statistical analysis, changes in relative levels of a particular neuropeptide in D. melanogaster tissue can be statistically detected amongst a variety of neuropeptides. While the data analysis methods should be compatible with other sample preparations, the presented sample preparation method was sufficient to identify previously unconfirmed D. melanogaster neuropeptides.
Neuropeptides are a large and diverse class of signaling molecules that affect numerous processes, including behavior, development, heart rate, metabolism, and reproduction [1, 2]. These peptides, mostly exerting their role by acting upon G-protein coupled receptors , can function as classical hormones, localized neurohormones , at muscles and glands, and synaptically, where they can modify the postsynaptic response to classical, fast-acting neurotransmitters . Studies of neuropeptide function that cross a wide variety of aspects of behavior and development have been particularly productive in the model organism Drosophila melanogaster[3, 6–14], which benefits from exceptional genetic manipulation tools developed for the study of the molecular mechanisms of development and behavior. In model insects such as D. melanogaster, Apis mellifera (honeybee), and Tribolium castaneum (red flour beetle), 30–40 genes have been consistently identified as encoding neuropeptides [8, 15, 16], with each gene product potentially producing multiple different mature neuropeptides. To become active, neuropeptides often require multiple post-translational modifications, such as proteolysis and amidation, which are difficult to infer from a genome and necessitate that putative neuropeptides be directly identified in organisms, often using mass spectrometry-based methods. Bioinformatics studies have predicted as many as 156 neuropeptides encoded by 33–119 putative neuropeptide genes in D. melanogaster, and a total of 76 neuropeptides from 21 genes have been detected experimentally [3, 17]. The sensitivity of MS-based methods has allowed for detection and identification of neuropeptides from specific nervous system regions and cellular populations across the developmental lifespan of D. melanogaster[1, 6, 7, 9, 18, 19] permitting precise temporal and spatial localization to be ascribed to various neuropeptides.
While great strides have been made towards comprehensive identification of D. melanogaster neuropeptides, functional characterization is lacking for many. For example, the majority of the peptides derived from the D. melanogaster gene neuropeptide-like precursor 1 (NPLP1) remain “orphaned” without an identified receptor and/or physiological function . Quantitative neuropeptidomics provides a discovery tool for ascertaining functional significance of neuropeptides, with goals of monitoring and quantifying changes in levels of multiple neuropeptides in response to experimental perturbations such as those eliciting complex behavioral responses. For example, isotope labeling followed by UPLC-ESI-QTOF has been used to quantify ~50 of known Apis mellifera brain peptides in the context of foraging, revealing molecular connections between the regulation of food intake in individual insects and this social behavior, as well as distinctions between nectar and pollen gathering . Isotopic labeling from extracts using MALDI-TOF MS combined with direct tissue MALDI imaging has been used to provide complementary information regarding changes in the expression of an array of neuropeptides during feeding in both the brain and pericardial organ of the crab Cancer borealis. A label-free LC-Orbitrap approach was employed to analyze extracts from hypothalamus and striatum from rats, using higher-energy collision dissociation and electron transfer dissociation fragmentation to identify more than 1700 endogenous peptides, revealing upregulation of orexigenic and anorexigeneic neuropeptides in animals fed on a high-fat/high-sucrose diet . Direct access to quantitative neuropeptidomics techniques, however, is often limited to laboratories equipped with a considerable array of specialized instrumentation, reagents, and personnel, preventing these methods from being more routinely utilized by those studying D. melanogaster development and behavior. Thus, we sought to develop a rapid method for performing differential expression neuropeptidomics studies utilizing D. melanogaster that does not require specialized reagents or advanced MS instrumentation. Furthermore, we wanted to present a data analysis workflow utilizing software that could preprocess and statistically analyze MS data regardless of instrument manufacturer.
Here we present a D. melanogaster sample preparation method that, when analyzed with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS), reliably detects an abundance of ions in the peptide mass range, 14 of which we subsequently confirmed by MALDI-TOF/TOF tandem mass spectrometry (MS/MS) to be D. melanogaster neuropeptides. Amongst the neuropeptides we identified by MS/MS fragmentation were three peptides derived from the NPLP1 gene not identified previously in the literature. Utilizing a MATLAB-based spectra preprocessing workflow, we demonstrate the ability to statistically detect differences in the expression of a specific neuropeptide, amongst all the ions we simultaneously observe, without isotopic labeling using MALDI-TOF MS.
Results and discussion
Straight-forward on-target peptide extraction provided adequate signal quality for MALDI-TOF MS profiling as well as targeted MALDI-TOF/TOF MS/MS
We set out to develop a sample preparation strategy for comparing neuropeptidomes from D. melanogaster that: could be performed in minutes, thus preserving labile biomolecules; could detect a large number of ions simultaneously, ideally with abundant enough signal to confidently identify using MALDI-TOF/TOF MS/MS; did not require extensive utilization of specialized reagents or equipment beyond a standard benchtop MALDI-TOF MS (at least for detection); and that utilized, ideally, only a single fly brain as an individual sample for statistical comparisons. The overall sample preparation we used consisted of dissection of individual D. melanogaster brains followed by their direct placement onto a steel MALDI target, an on-target wash, and matrix application. As a dissection medium, a modified dissection saline consisting of 7.5 g/L NaCl, 0.2 g/L KCl, 0.2 g/L CaCl2, and 0.1 g/L NaHCO3 in MilliQ water (pH 7.2)  was chosen as it yields a high number of detected ion signals and a relatively low baseline compared with dissection solutions with higher concentrations of salts. Solutions of fructose , Tris, and ammonium bicarbonate were evaluated as wash steps at various concentrations, with 100 mM ammonium bicarbonate producing the most abundant number of ions detected in the peptide mass range. Included in the detectable peaks with this sample preparation was m/z 1972.0, which we believed and later confirmed was the neuropeptide pigment-dispersing factor (PDF). This relatively low abundance peptide was used in subsequent experiments (see below) to demonstrate that known differences in relative levels of neuropeptides could be detected with this technique. With the ammonium bicarbonate wash, however, excess ammonium bicarbonate on the MALDI target after drying from the wash and subsequent matrix deposition was occasionally observed, which can interfere with homogenous crystallization of MALDI matrix and decrease the quality of acquired spectra (~10% of spectra). Alternatively, we found washing by dipping the brain in ammonium bicarbonate after dissection, but prior to placement on the MALDI target, as an alternative to the on-target wash. Data shown in the profiling experiments comparing flies with varying of PDF were acquired from samples prepared with the on-target wash, with spectra only acquired from samples that were not contaminated by excess ammonium bicarbonate.
Finally, various concentrations of 2,5-dihydroxybenzoic acid (DHB) and α-cyano-4-hydroxycinnamic acid (CHCA) were tested, ranging from 10–50 mg/mL and 5–10 mg/mL, respectively, as a matrix for MALDI-MS analysis, with 10 mg/mL CHCA providing the most reliable and highest quality spectra in terms of number of peaks with signal-to-noise (S/N) ratios greater than 6. During MS acquisition we noted that the entire crystalline matrix surface across a given spotted sample on the target did not yield homogeneous spectra. Specifically, both the surface of the brain itself and the edges of the MALDI spot yielded poor S/N ratios, leaving a “halo” region of high S/N spectra around the tissue (Additional file 1: Figure S1), consistent with what has long been observed in MALDI analyses of biological peptides from tissue samples . We attempted to improve homogeneity by trituration of the matrix solution; by homogenizing the brain on-target using a pipette tip or pressing with a cover slip; and by using a microcentrifuge tube homogenizer, but all of these procedures resulted in reduced S/N spectra. Overall, minimal mechanical perturbation of the brain was found to be important for achieving optimal quality spectra from single brains. Spectra of homogenized samples could be improved using a reversed phase “ZipTip,” but this required ~20 brains and added an additional step. As a result of the lack of MALDI spot homogeneity, we cannot be sure that peptides from all regions of the brain are extracted with identical efficiency. While this might prevent analysis of specific neuropeptides using this method, our primary goals of being able to detect simultaneously a variety of neuropeptides from various regions of the brain (as well as determine distinct differences in the levels of these neuropeptides when making experimental comparisons, with particular focus on PDF, as described below) were achieved.
Overall, raw spectra acquired with the final method yielded spectra with 37 ± 9.6 (mean ± S.D.) isotopically resolved peak distributions of S/N greater than 6 within the m/z 900–4000 range prior to any preprocessing (Figure 1A). Of importance was the ion signal at m/z 1972.0, confirmed in experiments below to correspond to the monoisotopic [M + H]+ of the neuropeptide PDF. Using the yellow white (yw, henceforth referred to as WT or wild-type flies) fly strain as a control strain, the PDF signal from m/z 1972.0 was often detectable in spectra from individual fly brains, but often only barely distinguishable from noise (Figure 1A inset). To demonstrate this sample preparation technique could be used in a differential neuropeptide profiling experiment to statistically identify differences in neuropeptide-derived ion signals, including neuropeptide-derived ion signals with low S/N ratios like PDF, flies expressing varying levels of PDF were obtained and analyzed. In particular, the pdf01 fly strain  (henceforth referred to as “PDF-null” flies) was used as a mutant strain lacking any expression of mature PDF and flies overexpressing PDF throughout the adult nervous system (referred subsequently in the text as “PDF overexpressing” flies) were generated using the GAL4-UAS binary expression system  by crossing the pan-neuronal elav-GAL4 driver line with a UAS-Drm-pdf line (see Methods section). Indeed, in spectra from individual flies, the ion signal corresponding to PDF was never observed in samples from PDF-null flies (Figure 1B), and was almost always observed with S/N > 6 in spectra from PDF overexpressing flies (Figure 1C). From this, we determined an experimental design that would permit changes in PDF to be statistically identified to validate this sample preparation method as a means for differential profiling of neuropeptides.
Preprocessing of spectra permitted statistical identification of distinct detection levels of ions corresponding to neuropeptide PDF
The most accurate methods of MS-based quantification generally involve the use of isotopologue standards . Relative quantitation of neuropeptide standards can be achieved over a thousand-fold concentration range on a MALDI-TOF mass spectrometer using isotopic labeling after careful selection of appropriate data acquisition parameters . Our method could be adapted to isotope dilution using any of the strategies demonstrated to quantify neuropeptides in single neurons or neuron clusters . However, in the current study, we show it is also possible to detect several distinct levels of expression of a neuropeptide in a label-free profiling approach using preprocessed MALDI-TOF MS spectra acquired from individual D. melanogaster brains. To achieve this, MALDI-TOF MS spectra from brains of WT flies were compared to PDF-null flies and PDF overexpressing flies. All brains were dissected from flies within a two hour window centered at Zeitgeber time 2 after entrainment to a 12:12-hour light–dark cycle at 25°C. Within an experiment, fly brains were analyzed individually, with multiple fly brains per genotype analyzed to permit statistical analysis of differences observed. Two full experiments were performed in order to further analyze the experimental reproducibility. In one experiment where the relative expression levels of PDF were compared, spectra were acquired from 10 WT, 7 PDF-null, and 9 PDF overexpressing D. melanogaster brains. In a second experiment, spectra were acquired from 9 WT, 7 PDF-null, and 9 PDF overexpressing D. melanogaster brains. Differences in sample numbers between experiments occurred due to spectra not being acquired from certain prepared samples judged to have poor crystallization, potentially from excess ammonium bicarbonate.
In order to use mass spectra obtained from individual fly brains for the purposes of differential neuropeptide expression analysis, a spectra preprocessing workflow  was employed that includes spectrum denoising, baseline correction, and normalization (see Methods section for full description of preprocessing). Peaks bins were then chosen from peaks identified in a total average spectrum, which was calculated from all spectra (after preprocessing) acquired across the three conditions (Figure 2A). Deisotoping criteria was applied in order to identify isotopically resolved peak distributions from the peaks detected in the total average spectrum. The inset of Figure 2A highlights the peak bins that were considered to be a single isotopomer distribution with a monoisotopic peak at m/z 1972.0 (i.e. PDF). While over 300 peaks were identified in the total average spectra of the two replicate experiments described here, after applying deisotoping criteria, exactly 57 isotopically resolved distributions were identified in either experiment, with 52 ion masses observed common to both replicate experiments (Additional file 2: Table S1).
Peak bins corresponding to isotopically resolved ion signals were used to query across all preprocessed spectra (Figure 2B). From this, a value for every ion signal isotopically resolved in the total average spectrum was assigned in the individual sample spectra regardless of whether that signal was detectable in any particular sample spectrum. Identifying ion signals in a total average spectrum, as opposed to individual spectra, may reduce the sensitivity of feature detection. For example, ion signals only detectable in a particular set of samples might be “averaged out” in the total average spectrum. However, by using a common set of peak bins across all spectra, some value can be established for every ion signal for the purpose of statistical comparisons, avoiding a “missing value” problem for peaks not otherwise “detectable” in a given spectrum. Summed deisotoped intensities were then compared with the non-parametric Kruskal-Wallis analysis of variance (ANOVA) test (α = 0.01). We adjusted for multiple comparisons using the straight-forward and conservative Bonferroni correction.
In the experiment comparing spectra from individual brains of 10 WT, 7 PDF-null, and 9 PDF overexpressing D. melanogaster, there was sufficient power to detect three significant differences between groups, including the distributions with monoisotopic peaks at m/z 1972.0 (Bonferroni-adjusted p-value = 0.0017, confirmed by MALDI-TOF/TOF MS/MS to be the [M + H]+ of PDF), m/ z 2010.0 (Bonferroni-adjusted p-value = 0.0014, inferred to be the [M + K]+ of PDF), and m/z 1203.6 (Bonferroni-adjusted p-value = 0.0099, m/z currently unassigned). Post-hoc analysis of the PDF species confirmed the difference was significant between all three groups, with levels increasing in the logical order (Tukey’s least significant difference procedure, α = 0.05). Figure 2C and D illustrate differences in ion signals corresponding to PDF compared to ion signal that did not vary significantly, m/z 1653.9 (Raw p-value = 0.22, subsequently identified as the [M + H]+ of IPNamide). Post-hoc analysis of m/z 1203.6 suggested detection of this isotope distribution was significantly reduced in the PDF overexpressing condition. In the second full profiling experiment performed, spectra were acquired from 9 WT, 7 PDF-null, and 9 PDF overexpressing flies, with only two isotope distributions being significantly different, the [M + H]+ of PDF (Bonferroni-adjusted p-value = 0.0097) and the [M + K]+ of PDF (Bonferroni-adjusted p-value = 0.0043). The significance of the difference in the unidentified m/z 1203.6 was not replicated in this experiment (Bonferroni-adjusted p-value = 0.1142, see Additional file 2: Table S1 for full results of statistics from both experiments). Post-hoc analysis of the PDF isotope distributions in this experiment again confirmed that PDF was detected at distinct levels across the three conditions in the expected order.
To further evaluate the general reproducibility of this method, the correlation between intensities of isotopically resolved ion signals detected in both experiments was examined. Focusing specifically on the fourteen ion signals later confirmed by MS/MS to be neuropeptides (see next subsection), the correlation between replicates of the mean intensities of these signals within the WT condition was high, with R2 = 0.969 (Pearson correlation, Figure 3). In the other two experimental conditions, the correlation between experiments was still generally high, with R2 = 0.894 for the PDF-null condition and R2 = 0.871 for the PDF overexpressing condition (Additional file 2: Table S2). Expanding this analysis out to all 52 isotopically resolved signals observed in both experiments, the correlation was generally high in each condition, with R2 = 0.957 for the WT condition, R2 = 0.848 for the PDF-null condition, and R2 = 0.914 for the PDF overexpressing condition (Additional file 2: Table S3).
Neuropeptide PDF was potentially detected as both an [M + H]+ at m/z 1972.0 and an [M + K]+ at m/z 2010.0 (as well as a [M + Na]+ at m/z 1994.0, which was not full resolved in the total average spectrum but is clearly present in the PDF overexpressing condition, marked by an asterisk in Figure 2). While both the [M + H]+ and [M + K]+ were found to be reproducibly significantly different between the three genotypes examined, the ratio of the relative intensities of these two forms of PDF ions were not necessarily consistent across an experiment. In particular, we evaluated the correlation between the intensities of the [M + H]+ and [M + K]+ of neuropeptide PDF within experiments (Additional file 2: Table S4). In the second experiment, there was a high correlation between the [M + H]+ and [M + K]+ of neuropeptide PDF in both the WT and PDF overexpressing conditions (R2 = 0.957 and 0.951, respectively). There was no correlation between the [M + H]+ and [M + K]+ signals in the PDF-null samples from the second experiment (R2 = 0.022) as would be expected given there should be no actual PDF signal in those samples. In the first experiment, however, while there remained a strong correlation between the [M + H]+ and [M + K]+ signals in the WT samples (R2 = 0.802), there was very low correlation between the [M + H]+ and [M + K]+ within the PDF overexpressing samples (R2 = 0.110), suggesting that perhaps residual potassium was unevenly distributed amongst samples in this experiment. Despite this inconsistency in the relative abundance of the potassium adduct of PDF compared with the [M + H]+, the differences in PDF levels between genotypes examined in this experiment were substantial enough to be detected whether either ionized form of PDF was considered. However, as this is a potentially confounding source of variability, it is important to consider whether detected changes in ion signal may be due shifts to different adduct ions, potentially from biases introduced during sample preparation such as inconsistent washing with ammonium bicarbonate.
MS/MS analyses confirm many ions detected are neuropeptide including three novel identifications
After utilizing MALDI-TOF MS profiling to identify differences in particular ions, the next logical step would be to identify what those ions are, preferably without additional sample preparation. Indeed, MALDI-TOF/TOF MS/MS data obtained using this sample preparation permitted identification of multiple neuropeptides, including PDF (Figure 4, Table 1). This was reassuring given that one concern might have been that the ions we observed were not in fact neuropeptides, but rather, for example, proteolytic fragments from abundant proteins. The rapidity of this sample preparation, preserving labile biomolecules, may be one reason we do not have this problem. In-source/post-source decay could also yield ions that obscure neuropeptide detection. Thus far, however, none of the 14 molecular ions analyzed by MS/MS appear to result from the breakdown of larger molecules during MS analysis.
Among identifications obtained were three previously predicted but unconfirmed peptides originating from the neuropeptide-like precursor 1 gene (NPLP1, Figure 5). Included in these identifications were two variants of the predicted NPLP1-3, without and without a C-terminal GAE (with the peptide lacking GAE being C-terminally amidated), and a peptide corresponding to the predicted NPLP1-2 but with a C-terminal lysine intact (NIATMARLQSAPSTHRDPK, or following previous convention, DPK peptide for short). It is possible the unamidated, glycine-extended NPLP1-3 (NVAAVARYNSQHGHIQRAGAE) is a precursor to the truncated and amidated variant (NVAAVARYNSQHGHIQRAa, or QRAamide), which may be the functional form of this peptide [31, 32]. Four other peptides derived from NPLP1 have been detected previously [1, 7, 33]. We also identify three of these, IPNamide, MTYamide, and VQQ, using MS/MS and tentatively assigned the fourth, APK peptide (measured monoisotopic m/z 1423.814/theoretical 1423.827). While the VQQ peptide of NPLP1 (NPLP1-4) has been identified as a ligand for receptor guanylate cyclase Gyc76C, serving a role in modulating the innate immune IMD pathway in response to salt stress , precise functions for the remaining NPLP1-derived peptides are unknown. The ability to monitor and detect changes in these peptides in response to various D. melanogaster experimental paradigms will hopefully provide insight into their potential significance.
Curiously, m/z 925.481, determined by MALDI-TOF/TOF MS/MS to be PDNFMRFamide (monoisotopic [M + H]+ = 925.435) was of a slightly higher mass error (50 ppm) compared with other identified peaks (average mass error = 4 ppm). Subsequent analysis with MALDI-Fourier transform ion cyclotron resonance (FTICR) MS, which has substantially higher resolving power compared with MALDI-TOF MS, confirmed this peak was actually a convolution of PDNFMRFamide and a peak tentatively assigned as Drostatin-3 (Ast-A3, monoisotopic [M + H]+ = 925.489, Figure 6), hence the comparatively larger m/z error. Including Drostatin-3, 30 additional molecular ions observed in at least some MALDI-TOF profiling spectra have been tentatively assigned by mass matching to be neuropeptides (Table 2). Of these 30 molecular ions, 19 were abundant enough to be detected (using our stringent deisotoping criteria) in both experiments presented here. This implies that perhaps at least 63% of the molecular ions detected in our MALDI-TOF MS profiling experiments are neuropeptides [(14 neuropeptides detected in both experiments identified with MS/MS + 19 additional peaks detected in both experiments tentatively assigned to neuropeptides)/52 isotopically resolved ion signals detected in both experiments], suggesting this method is highly specific for this type of biomolecule.
We acknowledge the rigor of isotopic labeling-based approaches to quantitative mass spectrometry and are aware of confounds inherent to MALDI MS-based quantification (i.e. differences in ionization efficiency of analytes, ion suppression effects, etc.). These notwithstanding, there are numerous examples of label-free MALDI-TOF MS-based methods providing informative semi-quantitative results [34–36]. Indeed, MALDI-TOF MS neuropeptide peak detection alone has been sufficient to distinguish particular cell types and tissues in Drosophila[6, 19]. Given the benefits of MALDI-TOF MS in terms of relative instrument expense and maintenance, as well as the ease of sample preparation and data acquisition, MALDI-TOF MS, even in instances when isotopic labeling is prohibitive, can serve as a valuable discovery tool, particularly when discovery of relatively more pronounced differences is an acceptable achievement (as opposed to absolute quantification). Here we present a rapid, label-free MALDI-TOF MS-based method and data analysis workflow that permits detection of differences in specific neuropeptides amongst a panel being monitored, using individual D. melanogaster brains as sample points. The described MATLAB-based preprocessing workflow and statistical analysis is compatible with other MALDI-TOF MS sample preparation techniques, including those previously described by other groups that have obtained spectra of excellent quality using other D. melanogaster tissues, including more specific D. melanogaster brain regions such as the antennal lobe  and individual cells . The relatively straight-forward sample preparation method described here was sufficient to enable both detection of distinct levels of neuropeptide expression, as well as identify previously unconfirmed neuropeptides. Similar to other discovery-based methodologies, we strongly recommend differences in ion intensities detected with the described technique be validated by an independent method (such as more quantitative MS and/or immuno-based methods). However, we believe this validation effort is worth the additional time given the relative ease of the initial discovery procedure.
D. melanogaster were reared on standard medium and raised under 12:12-hour light–dark conditions at 25°C. Flies were dissected between one to three hours after lights-on (two hour window centered at Zeitgeber time 2) when PDF expression levels are high . The pan-neuronal elav-GAL4 driver line was Bloomington stock #8760 (Bloomington Drosophila Stock Center at Indiana University, Bloomington, IN, USA).
Generation of UAS-Drm-pdf transgenic flies
Full-length D. melanogaster pdf- cDNA (Drm-pdf- cDNA) was kindly provided by Jeffrey C. Hall  and initially cloned into the pBluescript II SK (+/-) vector. To generate the UAS-Drm-pdf construct, the Drm-pdf- cDNA was then subcloned into the appropriate sites of the polylinker of the pUAST vector . The pUAST vector contains a P-element for which the transposase gene has been replaced by the sequences of the GAL4-specific UAS, the hsp70 TATA-box, the mini-white gene and the SV40 polyadenylation signal. The construct has been confirmed by direct sequencing using vector specific primers. Transgenic flies were generated by germline transformation following standard protocols. Briefly, the pUAST-Drm-pdf-cDNA construct and the transposase gene-containing helper plastmid pUChs ΠΔ2-3  were co-injected into Drosophila w1118 embryos using standard injection protocols. Two homozygous transgenic fly lines with different chromosomal localization of the construct (line no. 77: III chromosome) were obtained.
D. melanogaster brain dissection and on-target extraction
Flies were dissected in a modified insect dissection saline (NaCl 7.5 g/L, KCl 0.2 g/L, CaCl2 0.2 g/L, NaHCO3 0.1 g/L; pH 7.2)  and a single dissected fly brain was transferred with non-locking forceps (Dumont Tweezers #5, 11 cm, 0.025 × 0.005 mm tip, World Precision Instruments) to a stainless steel MALDI target. Excess dissecting saline was removed during the transfer of the brain. While the brain was still on the forceps, the forceps were touched about 2 mm from the fly brain in the void space between the two arms of the forceps, using a KimWipe™ paper (Kimberly-Clark Worldwide, Inc.), which wicked away excess liquid. An on-target wash of 1.0 μL of 100 mM ammonium bicarbonate was performed using a pipettor to add and aspirate the solution. 0.5 μL of 10 mg/mL CHCA in 50% (v/v) acetonitrile, 0.1% (v/v) formic acid was then directly pipetted onto the brain and allowed to dry before MALDI-MS analysis.
MALDI-TOF MS analysis of single dissected fly brains
Mass spectra were acquired on a microflex MALDI-TOF mass spectrometer (Bruker Daltonics Inc., Billerica, MA) equipped with a 337 nm N2 laser. Positive ion mass spectra were acquired from 500 m/z – 4000 m/z in reflectron mode. The acceleration voltage was set at 20 kV and the pulsed-ion extraction was set at 200 ns. One thousand laser shots were acquired for each spectrum. External mass calibration was achieved using a standard peptide mixture of Angiotensin I and II, Substance P, Renin Substrate, and ACTH (Bruker Daltonics Inc.). The externally calibrated mass accuracy of the instrument was approximately 100 parts-per-million (ppm) at m/z 1500.
MALDI-TOF/TOF MS/MS analysis for identification of detected peptides from single brain on-target extraction sample preparation
Fragmentation spectra were acquired in LIFT mode on an autoflex III and an ultraflex III MALDI-TOF/TOF mass spectrometer (Bruker Daltonics Inc., Billerica, MA). MS spectra were acquired in positive ion and reflectron modes. For MS/MS analysis, the source acceleration voltage was set to 8.0 kV and the reflectron voltage was set to 29.5 kV. Mass spectra were acquired with approximately 3000 laser shots summed in 200 to 400 shot increments. External mass calibration was achieved using a standard peptide mixture of Angiotensin I and II, Substance P, Renin Substrate, and ACTH (Bruker Daltonics Inc.). The external calibration mass accuracy of the instrument was approximately 20 ppm in MS mode and <400 ppm in LIFT (MS/MS) mode. MS/MS spectra were not internally calibrated. All spectra were processed with FlexAnalysis software (Bruker Daltonics Inc.).
Individual fly brain MALDI-TOF MS profiling data processing and statistical analysis
For preprocessing and statistical comparison of the relative abundance of PDF from MALDI-TOF profile spectra, all spectra obtained were first realigned by internal calibration using a mass list of neuropeptides we identified with MALDI-TOF/TOF MS/MS that were typically observed in our profile spectra. Internally calibrated spectra were exported to a plain text, two-column (m/z sampling points and corresponding intensities) ASCII format so they could be loaded into MATLAB 2013a (Mathworks, Natick, MA, USA). The msresample function in the MATLAB Bioinformatics Toolbox was used to resample internally calibrated spectra to a uniformly-spaced common set of m/z axis values.
Spectra from the two repeat experiments were analyzed separately. Spectra were denoised  in MATLAB using the undecimated discrete wavelet transform (UDWT) found in the Rice Wavelet Toolbox (http://www.dsp.ece.rice.edu/software/rwt.shtml) with a Daubechies’ scaling filter of length 8, soft thresholding applied, and the thresholding of low pass components enabled. Spectra were then baseline subtracted using the msbackadj function in the MATLAB Bioinformatics Toolbox with the default settings for this function. A total average spectrum of the denoised/baseline-subtracted spectra across all analyzed samples (including across conditions) within an experimental repeat was calculated. The total average spectrum was then normalized to its greatest value (i.e. the base peak). The mspeaks function from the MATLAB Bioinformatics Toolbox was then used to identify peaks from the total average spectrum, with the HeightFilter option of the mspeaks function set to “0.01”, so that only local maxima of the total average spectrum with intensity greater than 1% of the total average spectrum’s base peak are identified as peaks. When applying the mspeaks function, the left and right m/ z locations of the full width at half height (FWHH) limits for each identified peak in the total average spectrum was specified to be returned. These values were used to establish peak bins.
Individual denoised and baseline-subtracted spectra were additionally normalized by dividing each point in the spectrum by the sum of all intensities in the spectrum after denoising and baseline subtraction (i.e. a total ion current normalization, although we avoid the term here to specify total area under the spectrum is taken after denoising and baseline-subtraction). Peak bins from the total average spectrum were used to query across each spectrum being analyzed, with the maximum value in each peak bin returned for each spectrum. This reduced the set of all spectra being processed into a single n × m matrix, the entries of one dimension representing the n peaks found by mspeaks and the entries of the second dimension corresponding to the m spectra being analyzed.
Intensity values taken from certain peak bins were further added together by combining peak bins that corresponded to isotope distributions of the same singly charged ion. Criteria for combining peak bins corresponding to the same isotope distribution required that at least three peaks be detected with centroided m/z’s differing by 1 ± 0.03, with the peak corresponding to the monoisotopic peak in a distribution being the most abundant ion mass in the total average spectrum for distributions less than 1700 m/z and the second peak in a distribution always having to be more abundant than the third and fourth peak in the total average spectrum. Series of peaks meeting this criteria were assigned as “isotopically resolved”.
Non-parametric statistical analysis (Kruskal-Wallis one-way ANOVA) was performed in MATLAB, comparing the sum of the maximum intensities of isotopically resolved distributions in each spectrum across the three fly strains being analyzed. The p-value obtained from Kruskal-Wallis one-way ANOVA for each isotopically resolved distribution was adjusted using a simple Bonferroni correction based on the total number of isotope distributions being statistically compared so that a significance level of α = 0.01 could be used despite multiple comparisons. Tukey’s least significant difference procedure (α = 0.05) was used as a post-hoc analysis to evaluate differences between conditions of isotopically resolved peaks found to be significantly different in the ANOVA analysis.
MALDI-TOF MS/MS data processing for neuropeptide identification
MALDI-TOF/TOF MS/MS fragmentation spectra were exported to the FlexAnalysis software for batch preprocessing consisting of top-hat baseline subtraction; smoothing with four, width 0.15 m/z, Savitzky-Golay cycles; and peak picking using the SNAP algorithm with averagine molecular composition. All MS/MS spectra were combined and exported from FlexAnalysis as a single Mascot (Matrix Science Inc., London, England) generic file (*.mgf). The mgf file was submitted to an in-house Mascot server (version 2.2.07) for putative peptide identification using a 0.15 Da cutoff for precursor and a 0.5 Da cutoff for MS/MS peaks. Spectra were searched against the SwissProt 56.0 database with taxonomy specified as Drosophila (including 5357 protein sequences). Parameters for the Mascot search included enzyme specified as “none”, and variable modifications (C-terminal amidation, N-terminal pyroglutamic acid modification, methionine oxidation, and tyrosine sulfation) were considered. FlexAnalysis preprocessed MS/MS data were also moved to the BioTools software program (Bruker Daltonics Inc.) for manual confirmation of Mascot’s peptide assignments. A putative ID was considered confirmed when at least 3 consecutive b- or y-ions were observed, and in addition, the majority of MS/MS peaks were assigned. Although Mascot is tuned for protein-, not peptide-level identification (its significance scores at the peptide level are conservative), further confirmation came from the Mascot scoring algorithm in the form of a peptide score, peptide rank, and expectation value. For example, 13 of the 14 manually confirmed peptide identifications were also the highest ranking peptide from the Mascot search, and six of the 14 manually assigned peptides had Mascot Scores in the statistically significant range for protein identification.
MALDI-Fourier transform ion cyclotron resonance (FTICR)-MS was performed on an Apex Qe ultra 7 Tesla MALDI-FT-ICR mass spectrometer (Bruker Daltonics Inc., Billerica, MA). Mass spectra were collected in positive ion mode. The external mass accuracy of the instrument was approximately 20 ppm. After internal calibration, mass accuracy ranged from 0–2 ppm with a mean value of 0.5 ppm. Mass calibration was achieved using a standard peptide mixture of Angiotensin I and II, Substance P, Renin Substrate, and ACTH (Bruker Daltonics Inc.). Spectra were processed with DataAnalysis software (Bruker Daltonics Inc.).
Tandem mass spectrometry
Matrix-assisted laser desorption/ionization
Fourier transform ion cyclotron resonance
Neuropeptide pigment-dispersing factor
Neuropeptide-like precursor 1
Baggerman G, Boonen K, Verleyen P, De Loof A, Schoofs L: Peptidomic analysis of the larval Drosophila melanogaster central nervous system by two-dimensional capillary liquid chromatography quadrupole time-of-flight mass spectrometry. J Mass Spectrom. 2005, 40 (2): 250-260. 10.1002/jms.744.
Altstein M, Nassel DR: Neuropeptide signaling in insects. Adv Exp Med Biol. 2010, 692: 155-165. 10.1007/978-1-4419-6902-6_8.
Clynen E, Reumer A, Baggerman G, Mertens I, Schoofs L: Neuropeptide biology in Drosophila. Adv Exp Med Biol. 2010, 692: 192-210. 10.1007/978-1-4419-6902-6_10.
Strand FL: New vistas for melanocortins. Finally, an explanation for their pleiotropic functions. Ann N Y Acad Sci. 1999, 897: 1-16. 10.1111/j.1749-6632.1999.tb07874.x.
Merighi A, Salio C, Ferrini F, Lossi L: Neuromodulatory function of neuropeptides in the normal CNS. J Chem Neuroanat. 2011, 42 (4): 276-287. 10.1016/j.jchemneu.2011.02.001.
Neupert S, Johard HA, Nassel DR, Predel R: Single-cell peptidomics of drosophila melanogaster neurons identified by Gal4-driven fluorescence. Anal Chem. 2007, 79 (10): 3690-3694. 10.1021/ac062411p.
Baggerman G, Cerstiaens A, De Loof A, Schoofs L: Peptidomics of the larval Drosophila melanogaster central nervous system. The Journal of biological chemistry. 2002, 277 (43): 40368-40374. 10.1074/jbc.M206257200.
Schoofs L, Baggerman G: Peptidomics in Drosophila melanogaster. Brief Funct Genomic Proteomic. 2003, 2 (2): 114-120. 10.1093/bfgp/2.2.114.
Predel R, Wegener C, Russell WK, Tichy SE, Russell DH, Nachman RJ: Peptidomics of CNS-associated neurohemal systems of adult Drosophila melanogaster: a mass spectrometric survey of peptides from individual flies. J Comp Neurol. 2004, 474 (3): 379-392. 10.1002/cne.20145.
Hewes RS, Taghert PH: Neuropeptides and neuropeptide receptors in the Drosophila melanogaster genome. Genome Res. 2001, 11 (6): 1126-1142. 10.1101/gr.169901.
Vanden Broeck J: Neuropeptides and their precursors in the fruitfly, Drosophila melanogaster. Peptides. 2001, 22 (2): 241-254. 10.1016/S0196-9781(00)00376-4.
Pyza E, Meinertzhagen IA: The regulation of circadian rhythms in the fly’s visual system: involvement of FMRFamide-like neuropeptides and their relationship to pigment dispersing factor in Musca domestica and Drosophila melanogaster. Neuropeptides. 2003, 37 (5): 277-289. 10.1016/j.npep.2003.06.001.
Wang C, Zhang J, Tobe SS, Bendena WG: Defining the contribution of select neuropeptides and their receptors in regulating sesquiterpenoid biosynthesis by Drosophila melanogaster ring gland/corpus allatum through RNAi analysis. Gen Comp Endocrinol. 2012, 176 (3): 347-353. 10.1016/j.ygcen.2011.12.039.
Renn SC, Park JH, Rosbash M, Hall JC, Taghert PH: A pdf neuropeptide gene mutation and ablation of PDF neurons each cause severe abnormalities of behavioral circadian rhythms in Drosophila. Cell. 1999, 99 (7): 791-802. 10.1016/S0092-8674(00)81676-1.
Hummon AB, Richmond TA, Verleyen P, Baggerman G, Huybrechts J, Ewing MA, Vierstraete E, Rodriguez-Zas SL, Schoofs L, Robinson GE, et al: From the genome to the proteome: uncovering peptides in the Apis brain. Science. 2006, 314 (5799): 647-649. 10.1126/science.1124128.
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 (1): 113-122.
De Loof A: Ecdysteroids, juvenile hormone and insect neuropeptides: Recent successes and remaining major challenges. Gen Comp Endocrinol. 2008, 155 (1): 3-13. 10.1016/j.ygcen.2007.07.001.
Carlsson MA, Diesner M, Schachtner J, Nassel DR: Multiple neuropeptides in the Drosophila antennal lobe suggest complex modulatory circuits. J Comp Neurol. 2010, 518 (16): 3359-3380. 10.1002/cne.22405.
Yew JY, Wang Y, Barteneva N, Dikler S, Kutz-Naber KK, Li L, Kravitz EA: Analysis of Neuropeptide expression and localization in adult drosophila melanogaster central nervous system by affinity cell-capture mass spectrometry. J Proteome Res. 2009, 8 (3): 1271-1284. 10.1021/pr800601x.
Overend G, Cabrero P, Guo AX, Sebastian S, Cundall M, Armstrong H, Mertens I, Schoofs L, Dow JA, Davies SA: The receptor guanylate cyclase Gyc76C and a peptide ligand, NPLP1-VQQ, modulate the innate immune IMD pathway in response to salt stress. Peptides. 2012, 34 (1): 209-218. 10.1016/j.peptides.2011.08.019.
Brockmann A, Annangudi SP, Richmond TA, Ament SA, Xie F, Southey BR, Rodriguez-Zas SR, Robinson GE, Sweedler JV: Quantitative peptidomics reveal brain peptide signatures of behavior. Proc Natl Acad Sci USA. 2009, 106 (7): 2383-2388. 10.1073/pnas.0813021106.
Chen R, Hui L, Cape SS, Wang J, Li L: Comparative neuropeptidomic analysis of food intake via a multi-faceted mass spectrometric approach. ACS Chem Neurosci. 2010, 1 (3): 204-214. 10.1021/cn900028s.
Frese CK, Boender AJ, Mohammed S, Heck AJR, Adan RAH, Altelaar AFM: Profiling of diet-induced Neuropeptide changes in Rat brain by quantitative mass spectrometry. Anal Chem. 2013, 85 (9): 4594-4604. 10.1021/ac400232y.
Christie AE, Stemmler EA, Peguero B, Messinger DI, Provencher HL, Scheerlinck P, Hsu YW, Guiney ME, de la Iglesia HO, Dickinson PS: Identification, physiological actions, and distribution of VYRKPPFNGSIFamide (Val1)-SIFamide) in the stomatogastric nervous system of the American lobster Homarus americanus. J Comp Neurol. 2006, 496 (3): 406-421. 10.1002/cne.20932.
Garden RW, Sweedler JV: Heterogeneity within MALDI samples as revealed by mass spectrometric imaging. Anal Chem. 2000, 72 (1): 30-36. 10.1021/ac9908997.
Brand AH, Perrimon N: Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Development. 1993, 118 (2): 401-415.
Hou X, Xie F, Sweedler JV: Relative quantitation of neuropeptides over a thousand-fold concentration range. J Am Soc Mass Spectrom. 2012, 23 (12): 2083-2093. 10.1007/s13361-012-0481-0.
Rubakhin SS, Sweedler JV: Quantitative measurements of cell-cell signaling peptides with single-cell MALDI MS. Anal Chem. 2008, 80 (18): 7128-7136. 10.1021/ac8010389.
Morris JS, Coombes KR, Koomen J, Baggerly KA, Kobayashi R: Feature extraction and quantification for mass spectrometry in biomedical applications using the mean spectrum. Bioinformatics. 2005, 21 (9): 1764-1775. 10.1093/bioinformatics/bti254.
Wegener C, Reinl T, Jansch 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. J Neurochem. 2006, 96 (5): 1362-1374. 10.1111/j.1471-4159.2005.03634.x.
Czyzyk TA, Ning Y, Hsu MS, Peng B, Mains RE, Eipper BA, Pintar JE: Deletion of peptide amidation enzymatic activity leads to edema and embryonic lethality in the mouse. Dev Biol. 2005, 287 (2): 301-313. 10.1016/j.ydbio.2005.09.001.
Kolhekar AS, Roberts MS, Jiang N, Johnson RC, Mains RE, Eipper BA, Taghert PH: Neuropeptide amidation in Drosophila: separate genes encode the two enzymes catalyzing amidation. J Neurosci. 1997, 17 (4): 1363-1376.
Verleyen P, Baggerman G, Wiehart U, Schoeters E, Van Lommel A, De Loof A, Schoofs L: Expression of a novel neuropeptide, NVGTLARDFQLPIPNamide, in the larval and adult brain of Drosophila melanogaster. J Neurochem. 2004, 88 (2): 311-319.
Romanova EV, Lee JE, Kelleher NL, Sweedler JV, Gulley JM: Comparative peptidomics analysis of neural adaptations in rats repeatedly exposed to amphetamine. J Neurochem. 2012, 123 (2): 276-287. 10.1111/j.1471-4159.2012.07912.x.
Hanrieder J, Wicher G, Bergquist J, Andersson M, Fex-Svenningsen A: MALDI mass spectrometry based molecular phenotyping of CNS glial cells for prediction in mammalian brain tissue. Anal Bioanal Chem. 2011, 401 (1): 135-147. 10.1007/s00216-011-5043-y.
Hettick JM, Kashon ML, Simpson JP, Siegel PD, Mazurek GH, Weissman DN: Proteomic profiling of intact mycobacteria by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Anal Chem. 2004, 76 (19): 5769-5776. 10.1021/ac049410m.
Park JH, Helfrich-Forster C, Lee G, Liu L, Rosbash M, Hall JC: Differential regulation of circadian pacemaker output by separate clock genes in Drosophila. Proc Natl Acad Sci U S A. 2000, 97 (7): 3608-3613. 10.1073/pnas.97.7.3608.
Park JH, Hall JC: Isolation and chronobiological analysis of a neuropeptide pigment-dispersing factor gene in Drosophila melanogaster. J Biol Rhythms. 1998, 13 (3): 219-228. 10.1177/074873098129000066.
Rio DC, Rubin GM: Transformation of cultured Drosophila melanogaster cells with a dominant selectable marker. Mol Cell Biol. 1985, 5 (8): 1833-1838.
Coombes KR, Tsavachidis S, Morris JS, Baggerly KA, Hung MC, Kuerer HM: Improved peak detection and quantification of mass spectrometry data acquired from surface-enhanced laser desorption and ionization by denoising spectra with the undecimated discrete wavelet transform. Proteomics. 2005, 5 (16): 4107-4117. 10.1002/pmic.200401261.
This work was supported by the Department of Defense (contract W81XWH-04-0158) and Howard Hughes Medical Institute to M.R.; Amyotrophic Lateral Sclerosis Association (grant 1856) to J.N.A.; and Y.A.H. was supported in part by the National Institutes of Health (T32NS7292). We would like to thank Sergei Dikler and Jane-Marie Kowaslki from Bruker Daltonics, Inc. for their support and insight. We would like to thank Yuhua Shang for initial help with experiments and dissection.
The authors declare that they have no competing interests.
JNA and MR designed the experiment; KJB and YAH performed the experiments, with further contributions by JQ and AS; GG contributed the UAS-Drm-pdf line. PJK and MLE assisted with instrumentation. JNA, KJB, and JPS analyzed the data, prepared the figures, and wrote the manuscript. All authors have read and given approval to the final version of the manuscript.
Joseph P Salisbury, Kristin J Boggio contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: On-target extraction provides spectra with greater signal-to-noise and more peaks from the region surrounding the tissue, as opposed to acquiring spectra directly from the tissue. A) Acquiring spectra from directly over the deposited D. melanogaster brain (shown at the center of the crosshair encircled in red) did not provide quality spectra reliably. Rather, the region outside the red circle, which made up the visible matrix spot encircled approximately in orange, was where the best signal was obtained. B) Shows the same regions encircled with the crosshairs positioned over an area representative of a region that provides high and varied ion signal in the peptide mass range. (TIFF 4 MB)
Additional file 2: Table S1: This file contains supporting material, including the following tables. A complete list of isotopically resolved ion masses detected in either experiment (with neuropeptide assignments when possible), intensity means with standard deviations calculated within each group for each of these ion masses, and the adjusted p-values from Kruskal-Wallis ANOVA. Table S2. Correlation between experiments of fourteen neuropeptides confirmed by MS/MS fragmentation. Table S3. Correlation between experiments of all isotopically resolved ion signals observed in both experiments. Table S4. Correlation between [M + H]+ and [M + K]+ of neuropeptide PDF within experimental conditions for both repeat experiments. Format: XLSX (Excel Spreadsheet); Size: 45 KB. (XLSX 45 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Salisbury, J.P., Boggio, K.J., Hsu, YW.A. et al. A rapid MALDI-TOF mass spectrometry workflow for Drosophila melanogaster differential neuropeptidomics. Mol Brain 6, 60 (2013). https://doi.org/10.1186/1756-6606-6-60