Neuropeptidomic analysis of the embryonic Japanese quail diencephalon

Abstract

Background

Endogenous peptides such as neuropeptides are involved in numerous biological processes in the fully developed brain but very little is known about their role in brain development. Japanese quail is a commonly used bird model for studying sexual dimorphic brain development, especially adult male copulatory behavior in relation to manipulations of the embryonic endocrine system. This study uses a label-free liquid chromatography mass spectrometry approach to analyze the influence of age (embryonic days 12 vs 17), sex and embryonic day 3 ethinylestradiol exposure on the expression of multiple endogenous peptides in the developing diencephalon.

Results

We identified a total of 65 peptides whereof 38 were sufficiently present in all groups for statistical analysis. Age was the most defining variable in the data and sex had the least impact. Most identified peptides were more highly expressed in embryonic day 17. The top candidates for EE₂ exposure and sex effects were neuropeptide K (downregulated by EE₂ in males and females), gastrin-releasing peptide (more highly expressed in control and EE₂ exposed males) and gonadotropin-inhibiting hormone related protein 2 (more highly expressed in control males and displaying interaction effects between age and sex). We also report a new potential secretogranin-2 derived neuropeptide and previously unknown phosphorylations in the C-terminal flanking protachykinin 1 neuropeptide.

Conclusions

This study is the first larger study on endogenous peptides in the developing brain and implies a previously unknown role for a number of neuropeptides in middle to late avian embryogenesis. It demonstrates the power of label-free liquid chromatography mass spectrometry to analyze the expression of multiple endogenous peptides and the potential to detect new putative peptide candidates in a developmental model.

Background

Bird models have been instrumental to the overall understanding of neural sex differences and endocrine influences on brain development [1, 2]. Japanese quail is a commonly used bird model for studying sex-specific brain development and behavior, especially the influence of the hormonal milieu during early brain development on sexual behaviour [3–5]. Sex-specific neural development has traditionally been associated with sex hormones produced by the gonads. Exogenous estrogen exposure before embryonic day 12 (ed12) causes demasculinization of the male sexual behaviour in Japanese quail [3, 5], which shows the important role of sex hormones in brain differentiation. Several studies in songbirds have shown that there are exceptions to this classical model [6–9], implicating an intrinsic genetic influence on sex-specific neural development. Our previous studies on gene expression in early chicken [10] and quail embryos (ongoing) also indicate a sex-specific brain development that is at least partly independent of the influence of gonadal hormones. Quail ed12 corresponds roughly to chicken ed14 which is close to the developmental stage when gonad derived plasma testosterone levels peak (~ed13.5) in male embryos and the endocrine hypothalamic-pituitary-gonadal (HPG) axis becomes established [11–13]. Endogenous gonad derived plasma estradiol levels increase gradually in female chicken embryos until ~ed13.5 where they increase markedly [14, 15]. Male chicken embryo plasma estradiol levels are significantly lower than female levels at all stages although the differences become more pronounced after ed13.5 [15]. The establishment of the male HPG axis is believed to be dependent on both hypothalamic gonadotropin-releasing hormone release and opioid peptide regulation [12]. Several endogenous peptides, including opioid peptides, are messengers with sex-specific activity, some of them being stored and secreted in an endocrine-like manner [16]. Opioid peptide dependent pain and analgesia [17] are well known examples of neural sex differences. Pain and analgesia are conditions that are sensitive to neonatal alterations in gonadal steroid hormones [18]. Studies using the opioid receptor antagonist naloxone on the medial preoptic area (POM) in the anterior diencephalon point to a role for opioid peptides in the regulation of both avian and mammalian sexual behavior [19–21].

A group of neuropeptides may have the same pre-cleavage protein precursor in common, meaning that one gene may code for several peptides with different functions. Although genome information from different species has facilitated prediction of putative precursor cleavage sites and peptide/precursor amino acid sequences, little is known about the expression levels of different peptides under different conditions. This is especially true regarding the peptidome content in the developing brain. The sequencing of the chicken genome has enabled prediction and expression analysis of neuropeptides in avian models such as Japanese quail. Neuropeptidomics is the technological approach for detailed analysis of endogenous peptides in the brain. Recently, a novel approach to study a large number of neuropeptides using nL/min flow liquid chromatography (nanoLC) electrospray ionization mass spectrometry was developed to investigate the endogenous neuropeptide content of brain tissue samples from rats and mice [22–24].

The aims of this study were to characterize the endogenous peptide content of the diencephalon including the sexually dimorphic area POM in the developing quail brain and to analyze how its expression is related to age, sex and the effect of ethinylestradiol (EE₂) exposure. Quail ed12 diencephalon was chosen as both the stage of development and the tissue are of interest in this context. Ed12 is at the end of the sensitivity period for embryonic estrogen exposure leading to behavioral alterations in male copulatory behavior [5] and likely in the time period when the male endocrine HPG axis becomes established and a more marked increase of female estradiol synthesis occurs. Quail ed17 is somewhat before hatching and a time where the male embryos lost their sensitivity to exogenous estradiol and their sex specific brain and neuroendocrine development has become established. The detection of neuropeptides and the quantification of their expression levels were done using liquid chromatography and ion trap mass spectrometry together with a label-free approach.

Results

Peptide identification

Japanese quail embryonic diencephalon protein samples from 38 individuals were separated by molecular weight, obtaining < 10 kD samples. These samples were then examined using nano LC coupled to an ESI-LTQ MS system or an LTQ-FTICR system. A total of 65 peptides were identified (Table 1) of which 38 (Table 2) were present in a sufficient number of individuals to enable expression profiling (see below). Among the 65 peptides, 61 were derived from neuropeptides or neuropeptide precursor proteins. The four remaining peptides belonged to CRMP-2, FKBP5 (also called PPIA), Stathmin and Thymosin beta.

Table 1 Peptides identified by LTQ and FTICR-MSMS in quail diencephalon.

Table 2 Differentially expressed peptides

Expression profiling

A subset of the LTQ-MS samples were divided into eight groups by sex (male vs female), age (ed12 vs ed17) and EE₂-exposure (tissues collected from ed12 and ed17 males and females) with four to five samples in every group. The samples were ed12 male controls (n = 5), ed12 female controls (n = 4), ed12 male EE₂ exposed (n = 4), ed12 female EE₂ exposed (n = 4), ed17 male controls (n = 4), ed17 female controls (n = 4), ed17 male EE₂ exposed (n = 4) and ed17 female EE₂ exposed (n = 5). The alignment of peptide spectra was performed using the DeCyder MS software. Expression data for 204 peptides, whereof 38 identified, was first analyzed by unsupervised methods (hierarchical clustering, principal component analysis, PCA) and then by a supervised method (Partial Least Squares Discriminant Analysis, PLS/DA). A clear separation between ed12 and ed17 peptide expression was visible in a hierarchical clustering (Figure 1A) of all eight experimental groups with the majority of the identified peptides being more highly expressed in ed17 than in ed12 diencephalon (Figure 1B). A similar pattern was detected with PCA (Figure 1C, additional file 1).

Figure 1

Unsupervised hierarchical clustering and principal component analysis. (A) Visualization of group differences (ed12MC, ed12FC, ed12ME, ed12FE, ed17MC, ed17FC, ed17ME, ed17FE) for 204 peptides, using two-way clustering of the log2 ratio (Xclass-Xtot) between the group median value (X_class) and the total median value for every peptide (X_tot). (B) One-way clustering of 38 identified peptides. Positive values are colored red, negative are green. (C) Principal component analysis using two components (C1, C2). The last letter(s) in the group definitions ed12/17XX stand for M = male, F = female, C = controls, E = ethinylestradiol exposed. # designates the peptide numbers in the data set.

Both analyses indicate a weak separating effect of EE₂ exposure in ed12 samples. PLS/DAs were conducted to explore the maximized separation between groups regarding age, sex and EE₂ exposure. There were no successful models separating all three parameters (sex, age and EE₂) leading to a PLS/DA model for age and EE₂ exposure (Figure 2, additional file 2). The effect of EE₂ exposure was more distinct on ed12 than on ed17 (Additional file 3). A three-dimensional PLS/DA score plot (Figure 2A) shows the separation of age (ed12 and ed17) and EE₂ exposure (controls, C and exposed, E), independently of sex (ed12C, ed12E, ed17C, ed17E), with the second and third component being most effective in separating controls from EE₂ exposure (Figure 2B-D). Next, a virtual importance in the projection (VIP) (Figure 3) was plotted to see how much the identified and unidentified peptides influence the different class separations in the PLS/DA. The preprotachykinin 1 peptide neuropeptide K (PPT NPK, VIP 1.59) was most important for overall separation among the identified peptides, followed by a peptide derived from CRMP-2 (VIP 1.45). These two peptides also display a significantly interaction effect when using moderated F-statistics (see below) (Table 2). Four unidentified peptides (VIP 1.6-1.85) were more relevant than NPK (data not shown). NPK, the top PLS/DA candidate for overall effects on classes based on age and EE₂ exposure, was also the most influential among ed12 samples (see Additional file 3).

Figure 2

PLS/DA on age and EE ₂ exposure using 204 quail peptides. (A) Three-dimensional score plot for all observations using the first three components (explaining ~75% of the data) in a model on age and EE₂ exposure effects. The PLS/DA model uses the dummy variables DA1 (ed12C, red color), DA2 (ed12E, green color), DA3 (ed17C, blue color), and DA4 (ed17E, black color). In (B), the second and third components are combined in a score scatter plot, separating controls (red and blue triangles) from EE₂ individuals (green and black boxes). The influence of identified peptides on the positioning of samples in the score scatter plot is seen in the loading scatter plot (C). The 38 identified peptides in the loading plot are listed in (D). Score/loading axes for DC2 and DC3 are designated t(2)/w*c(2) and t(3)/w*c(3) respectively. ed12/17C designates ed12/17 controls and ed12/17E designates ethinylestradiol exposed.

Figure 3

Virtual importance in the projection plot (VIP). The VIP shows the role of identified peptides in the separation of the four age and EE₂ exposure PLS/DA classes (ed12C, ed12E, ed17C, ed17E). The preprotachykinin-1 neuropeptide K (NPK) is the top candidate among identified peptides for being affected by both age and EE₂ exposure. Values above one are the most relevant for the separation of all classes. Unidentified peptides are not shown. Error bars designate confidence intervals derived by jack knifing.

Multivariate analysis methods such as PCA and PLS/DA are efficient for illustrating general trends and patterns in complex data but more limited in determining the significance/effect for single factors (peptides). An F-test was used to detect statistically significant (adjusted p < 0.05) differences in singular peptide expression between groups regarding the response variables age, EE₂ exposure and sex. The age differences among identified peptides were dominated by a general upregulation between ed12 and ed17 (Figure 1B, Table 2). There were a few exceptions, such as NPK and the cerebellin derived peptides (Table 2). Only few peptides were significantly affected by either sex or EE₂ exposure (Table 2). The top EE₂ candidates were NPK and the CRMP-2 derived peptide (Figure 4). This was also indicated by their localization (peptides #118, 217) in the PLS/DA VIP (Figure 3) and the loading scatter plot (Figure 2C-D) where they are most prominent in separating control samples (ed12C, ed17C) from EE₂ samples (ed12E, ed17E). The CRMP-2 peptide displays EE₂ induced downregulation of expression in males at ed12 and ed17 (in average ~57%), whereas NPK is downregulated in both males and females at ed12 and ed17 (in average 50%). Peptides derived from gastrin-releasing peptide (GRP; APLQPGGSPALTKIYPR) and gonadotropin inhibitory hormone related peptide 2 (GnIH-RP2) were the top candidates for sex specific differences (Figure 4). GRP was more highly expressed (in average 66%) in all male groups compared to their respective female groups independent of EE₂ exposure (Figure 4) and GnIH-RP2 was more highly expressed (in average 52%) in ed12 males and in ed17 control males.

Figure 4

Candidate diencephalon peptides sensitive to sex differences or EE ₂ . Visualization of the differences of candidate peptides (NPK, CRMP-2, GnIH-RP2, and GRP) between all eight groups (ed12MC, ed12FC, ed12ME, ed12FE, ed17MC, ed17FC, ed17ME, ed17FE) as calculated (Xclass-Xtot) for the hierarchical clustering in figure 1. The stars designate significance levels (*** p < 0.001, ** p < 0.01 and * p < 0.05) for the different factors in the model. ^a FDR adjusted p-value significance for F-test. Interaction effects were tested using Wald F-test. The last letters in the group definitions ed12/17XX stand for M = male, F = female, C = controls, E = ethinylestradiol exposed.

Secretogranins

Peptides from three secretogranins (SCG1, SCG2, and SCG5) were identified (Table 1). Two out of six SCG2 peptides were detected in all samples. The two SCG2 derived peptides (TNEIVEEQYTPQSL and SGKLSFLEDE) behave in different manners as seen by F-test statistics (the former being non-significant and the later significant, F-test statistics, adjusted p < 0.05). Both show an age effect, becoming more highly expressed in ed17 but only SGKLSFLEDE shows a significant interaction effect and approaches significance for a treatment effect (t-test p-value 0.062). SGKLSFLEDE is located between two basic amino acid cleavage sites (KR, RR) and corresponds to human SCG2 285-295 (Figure 5). The full peptide probably includes an additional C-terminal methionine (SGKLSFLEDE M) if one considers the full sequence between the cleavage sites. Such a peptide was also discovered (Table 1) but its presence across all samples was insufficient for matching purposes. The amino acid sequence of the secretogranin-1 (chromagranin B) peptide (IHEGEEGEAEEE) is less conserved between birds and mammals and the peptide as such seemingly unaffected by age. The peptide is located in a non-conserved sequence region of secretogranin-1 (data not shown) and no deduction about relevance or functionality was possible.

Figure 5

Japanese quail secretogranins. (A) Overview of avian secretogranin-2 (based on chicken sequence XP_422624) showing its granin protein domain and the position and multiple sequence alignment to other species for two identified peptides: TNEIVEEQYTPQSL (Secretoneurin peptide) and SGKLSFLEDE. Bt is Bos taurus, Ss is Sus scrofa, Hs is Homo sapiens, Mm is Mus musculus, Gg is Gallus gallus, Xl is Xenopus leavis, and Dr is Danio rerio. (B) MS/MS fragmentation graph for SGKLSFLEDE in GPM browser. RI stands for relative intensity. C, Visualization of expression differences between the groups as calculated (Xclass-Xtot) for the hierarchical clustering in figure 1. The stars designate significance levels (*** p < 0.001, ** p < 0.01 and * p < 0.05) for the different factors in the model. ^a FDR adjusted p-value significance for F-test. Interaction effects were tested using Wald F-test. The last letters in the group definitions ed12/17XX stand for M = male, F = female, C = controls, E = ethinylestradiol exposed.

Protachykinins

Among the peptides identified, several were part of the protachykinin 1 precursor (designated PPT). Representative peptide sequences for all peptide forms (substance P [SP], NPK, neurokinin A [NKA] and C-terminal flanking peptide [CTFP]) except neuropeptide gamma were detected. FTICR-MS enabled us to detect two previously unknown phosphorylations on the C-terminal flanking peptide (CTFP), on S6 or S7 (single phosphorylated; Figure 6A) and S10 (double phosphorylated with S6/7; Figure 6B). The detection of the truncated NPK sequence DAGYGQISH (Figure 6C) shows that quail has the same sequence deletion in its Protachykinin gene as seen in chicken compared to mammals. The log2 ion intensity levels for non-phosphorylated and single-phosphorylated CTFP isoforms were relatively similar (log2 ion intensity ~18-19) across all groups whereas the expression of the double-phosphorylated isoform was much lower (log2 ion intensity ~12-13).

Figure 6

The preprotachykinin-1 peptides. SLNSGSSERSIAQNYE is similar to the C-terminal flanking peptide (CTFP) in mammalian protachykinin 1. Two previously unknown serine phosphorylations at positions S6 or S7 and S10 were discovered. (A) FTICR derived MS/MS spectra for single phosphorylated form (Mw 1820) at S6 or S7 and (B) S10 (Mw 1900.3). Amino acid sequence alignment (C) between chicken (continuous sequence in black color) and human (red text) preprotachykinin 1 precursor shows the differences between birds and mammals in the detected tachykinins, especially neuropeptide K (NPK). (D) Visualization of differences between the groups as calculated (Xclass-Xtot) for the hierarchical clustering in figure 1. The stars designate significance levels (*** p < 0.001, ** p < 0.01 and * p < 0.05) for the different factors in the model. ^a FDR adjusted p-value significance for F-test. Interaction effects were tested using Wald F-test. The last letters in the group definitions ed12/17XX stand for M = male, F = female, C = controls, E = ethinylestradiol exposed.

Discussion

There is very little known about neuropeptide expression patterns and regulation in the developing embryonic brain. This is the first peptidomics study on neuropeptide expression in the developing avian brain. A total of 65 peptides from the embryonic quail diencephalon were identified (Table 1). Using the minimum matching criteria of max one missing sample per group resulted in 204 matched peptides, including 38 identified. These were plotted for differences related to age (ed12 and ed17), sex and EE₂ exposure. The general patterns in the data set were explored using unsupervised and supervised multivariate analysis and unsupervised hierarchical clustering. These methods demonstrated a clear age separation between ed12 and ed17 individuals. The fact that most identified peptides were more highly expressed in ed17 compared to ed12 was not representative for all plotted peptides (Figure 1A-B, Table 2). The PCA and hierarchical clustering also revealed a less distinctive but still discernable separation of control and EE₂ exposed individuals, although mainly at ed12. Overall, the data indicates that the ed12 quail diencephalon shows a more distinct EE₂ sensitive and sex peptide expression pattern as compared to ed17 (Figure 1, Additional file 3). If true, it may be correlated with the establishment of the HPG axis occurring at this developmental stage in male quail embryos [12, 13]. No overall sex-specific patterns were detected, which was somewhat unexpected considering the sexual dimorphic nature of the diencephalon and the large sex differences observed in birds on the mRNA level [10, 25, 26]. The mRNA differences are mainly based on gene dosage effects from the Z chromosomes (males ZZ; females ZW). Only one peptide (GRP) in our data was derived from a gene located on the Z chromosome, possibly explaining the absence of large scale sex differences. Interestingly, the sex chromosome Z derived GRP peptide was the top candidate for sex specific peptide expression and was significantly more highly expressed in males compared to females in both controls and EE₂ exposed individuals. The absence of general sex differences was also unexpected in relation to the detection of candidates for the demasculinizating EE₂ effects on the copulatory behavior in adult males following embryonic estrogen exposure [5, 27].

Most of the identified peptides were part of larger neuropeptide precursor proteins. Peptides derived from the same secretogranin-2 precursor displayed expression differences indicating peptide specific regulation. Both peptides were more highly expressed in ed17 but only SGKLSFLEDE was significantly affected over all groups (F-test statistics, adjusted p < 0.001) with a significant interaction effect and near significant EE₂ treatment effect (p-value 0.062). The positioning of SGKLSFLEDE between basic amino acid pairs (KR and RK, Figure 5A) and its expression pattern compared to Secretoneurin indicates that it possesses a specific biological function and makes it a potential new avian neuropeptide candidate. Besides SGKLSFLEDE, two other forms (SGKLSFLE and SGKLSFLEDEM) were also detected, although not in all the samples. The later, SGKLSFLEDEM is probably the complete peptide as it represents the complete sequence between the two cleavage sites in the precursor protein.

Previous studies have demonstrated that opioid peptides are involved in aggressive and sexual behaviors in adult Japanese quail [21, 28] and also in the establishment of the male chicken embryonic HPG axis [12]. We detected a large number of proenkephalin precursor (PENK) derived peptides (Table 1, Additional file 4), although the majority were propeptides located in front of the opioid peptides. None of the identified opioid peptides (MERF or MERSL) differed significantly between groups except with regard to age. The detection of Met-Enkephalin (Met-Enk) and Leu-Enkephalin (Leu-Enk) propeptides indicates that Met-Enk and Leu-Enk also are expressed in the developing diencephalon. Their absence may be explained by the fact that their short sequences make it more difficult to detect them with LC-MS. The presence of the complete Leu-Enk propeptide (VGRPEWWLDYQ) in all groups and it being significantly affected by age and EE₂ exposure plus exhibiting interaction effects (Table 2) may indicate that Leu-Enk also is differentially expressed. Among the peptides derived from another large neuropeptide precursor protein, Preprotachykinin 1, we identified two previously unknown phosphorylations on the C-terminal flanking peptide (CTFP), on S6 or S7 (single phosphorylated; Figure 6A) and S10 (double phosphorylated with S6/7; Figure 6B). Quail SP possesses the same Arg for Lys substitution at the third position compared with humans as previously described in chicken [29], alligator [30] and Burmese python [31] whereas quail NKA is identical to chicken, alligator, Burmese python and human NKA [31]. NKA is the ten amino acid long C-terminal part of NPK (36 amino acids) in mammals. The LC-MS detection and identification of the truncated NPK sequence DAGYGQISH (Figure 6C) shows that quail has the same sequence deletion in its Protachykinin gene compared to mammals as seen in chicken. Antisera against a mammalian NPK fragment which encompasses some of the chicken NPK sequence shows immunoreactivity in the developing chicken embryonic spinal cord which weakens after hatching [32]. The nine amino acid long quail NPK could work in a functionally inert propeptide fashion to NKA in quail considering that NPK is converted into NKA in mammals [33]. It is not clear if the truncated avian NPK has a biological function although it is indicative that it is the top candidate for influencing overall separation in the multivariate datasets with regard to age and EE₂ exposure, especially at ed12 (Figures 2 and 3, Additional file 3) and becomes significantly downregulated in all groups after EE₂ exposure (Figure 4, Table 2). The different age dependent expression patterns seen in our data for NPK and NKA (NPK is unaffected by age whereas NKA is significantly upregulated) can also be interpreted as involvement in separate biological processes during diencephalon development. NPK is known to regulate gonadotropin secretion in the mammalian hypothalamus [34], and our data indicates that this may also be true for avian diencephalon development.

There were only few peptides exhibiting sex differences in expression, the top candidates being the GRP and GnIH-RP2 peptides. The first 17 amino acids of GRP was as previously mentioned the only ChrZ-derived peptide detected in our data and it was generally more highly expressed in males compared to females (in average 66%). GRP is an interesting candidate for sex specific brain development as expression of its receptor has been related to the layering of cells in the chick telencephalon, optical tectum and cerebellum [35]. The GnIH-RP2 peptide is part of the gonadotropin inhibitory hormone (GnIH) precursor protein (UniProt GNIH_COTJA) whose main peptide (GnIH) inhibits gonadotropin release in quail and is strongly expressed in all avian species studied so far [36, 37]. The GnIH-RP2 peptide was significantly differentially expressed between groups with a clear age effect from ed12 to ed17 (Figure 4, Table 2). It also displayed significant sex and interaction effects. It should be noted that the identified sequence (ARSSIQSLLNLSQ) is actually somewhat longer than the calculated GnIH-RP2 sequence (SSIQSLLNLSQRF) and does not posses the characteristic C-terminal RF-amide. The observed peptide may be part of a larger propeptide-GNIH-RP2 peptide sequence that has become degraded at both the N and C-terminal. The GnIH precursor protein contains a neuropeptide cleavage site (RxnR where n = 2, 4 or 6) [38] adjacent to a somewhat longer sequence (SPLARSSIQSLLNLSQ) of the peptide in question (Table 2, Additional file 4). Unlike GnIH itself, there is little known about the role of GnIH-RP peptides. Although this may be the case, the overall expression pattern of this particular peptide implies a biological role in avian embryonic diencephalon development for GnIH-RP2. Both GnIH-RP peptides posses affinity for the GnIH receptor [39], indicating that GnIH-RP2 is involved in the regulation of gonadotropin release from the developing pituitary. NPK and GnIH-RP2 are both potential regulators of gonadotropin regulation during diencephalon development and the male HPG axis formation. The EE₂ effects on GnIH-RP2 expression are not very influential/statistical significant in either type of data analysis (multivariate or univariate) and therefore not conclusive. Still, the fact that EE₂ also seems to modulate GnIH-RP2 expression (Figure 4) is intriguing considering that the NPK expression is EE₂ sensitive. It could indicate that the male quail embryonic window of sensitivity (regarding adult male copulatory behavior) before and around ed12 to exogenous EE₂ exposure is partly influenced by an altered gonadotropin regulation and HPG axis formation.

It is important to note that this study is foremost of hypothesis generating nature. We have analyzed the expression of a large number of neuropeptides (including detailed amino acid sequence information) under different circumstances in the developing avian brain in a manner that is impractical or impossible for the more traditional immunoassay methods. Omics data is valuable when the experimental design includes several variables (in our case age, sex and EE₂ exposure corresponding to three variables or eight experimental groups) but requires at the same time an increasing number of samples and makes the analysis and presentation of the results more difficult. Further studies are needed to better analyze the role of the individual neuropeptides during development using for instance more data points for age or the addition of aromatase inhibitors to better deduce the impact of EE₂. The absence of neuropeptides such as vasotocin was for instance unexpected although this does not provide any conclusions on their presence or role in embryonic brain development. The absence of proof is not the proof of absence. Peptides such as vasotocin contain two cysteine amino acids and can contain a disulfide bridge, a property that makes them more difficult to fragment and identify with MS/MS. Additional peptides could for instance be detectable in less hydrophilic sample preparation buffer settings than the weakly acidic 0.25% acetic acid buffer used in this study. Future peptidomics studies such as this study may provide even more data when other sample preparations strategies or additional mass spectrometry systems besides LTQ-MS are added.

Conclusion

To summarize, we report the first large scale neuropeptidomic analysis of the developing avian diencephalon at ed12 and ed17 in Japanese quail, a common model for studying sex-specific neural development. The peptidome was analyzed in relation to age, sex and EE₂ treatment at ed3. Both age and EE₂ treatment affected the quail peptidome, without any clear discernible influence of gender. A potential new avian neuropeptide candidate was discovered in the secretogranin-2 precursor protein and previously unknown phosphorylation sites in the protachykinin C-terminal flanking peptide identified. GnIH-RP2 and GRP are likely to be differentially expressed in a sex specific manner and NPK in an EE₂ sensitive manner during embryonic brain development in Japanese quail. Our approach shows the strength of labeling-free LC-MS peptidomics and how to better characterize multiple known and unknown endogenous peptides in the developing organism.

Methods

Embryos and sample collection

Fertilized Japanese quail (Coturnix coturnix japonica) eggs were obtained from a local breeder. The eggs were incubated at 37.5°C and 60% relative humidity and turned every 3 h. On day three of incubation, eggs were injected with 300 ng ethinylestradiol (EE₂) each. The EE₂ was dissolved in an emulsion of peanut oil, lecithin, and water [40]. Controls received emulsion only. The holes in the shells were sealed with melted paraffin wax and the eggs were returned to the incubator. The diencephalon was excised as quickly as possible from embryonic brains on ed12 and ed17 and was placed in pre-weighed Eppendorf tubes. All samples were quickly frozen in liquid nitrogen, and stored at -80°C. These samples were then used for neuropeptidomic analysis. A tissue sample from each embryo was also collected for DNA isolation and genetic sexing according to a PCR-based method [41] in which intron sequences of different lengths in the W-linked gene CHD1W (females) and Z-linked gene CHD1Z (both sexes) are amplified. All animal experiments were approved by the local Ethics Committee for Animal Research, Uppsala, Sweden (reference number C 198/6).

Sample preparation

Each frozen tissue sample was weighed and then placed on aluminum foil and put into the Denator Stabilizer instrument (prototype acquired from Denator Biotechnology, Gothenburg, Sweden) for heat denaturation for 45 seconds to protect the samples from post mortem protein degradation [42]. The order for denaturation of the samples was randomized. The tissue samples were then transferred to low-retention Eppendorf tubes and suspended in pre-chilled extraction solution (0.25% acetic acid; 0.2 mg tissue/μl) and homogenized by sonication (Vibra cell 750, Sonics & Materials Inc., Newtown, CT, USA) for 30 seconds [22, 23]. Each sample suspension was then centrifuged at 20,000 g for 30 minutes at 4°C to remove insoluble material. The supernatant was transferred to Microcon 10 kDa cut-off spin columns (YM-10, Millipore, Bedford, MA, USA) and centrifuged for 45 minutes at 14,000 g at 4°C. The resulting peptide filtrate was then frozen at -80°C until analysis.

Experimental design

All samples were run in a complete randomized block fashion with a total of five blocks (Table 3). Each block included replicates of the same sample (Q52) that was run repeatedly (six times, designated 52a-f). This resulted in a total of 42 runs, using 37 unique samples.

Table 3 Sample annotation and run order design

MS analysis

Five μl of each sample was desalted on a Nano-Precolumn (LC Packings, Amsterdam, Netherlands) using Ettan MDLC (GE Healthcare). The peptides were then separated by a 40 minutes long run with a gradient from 3 to 80% acetonitrile in 0.25% acetic acid on a 15 cm long Proxeon nanoESI emitter (75 μm inner diameter, 360 μm outer diameter) packed in house with Reprosil-Pur C18 3 μm resin (Dr. Maisch GmbH, Ammerbuch-Entringen, Germany). At a flow rate of approximately 150 nl/min the peptides were electro sprayed into a linear ion trap mass spectrometer (LTQ, Thermo Electron, San Jose, CA, USA). The spray voltage was 1.8 kV, ES source capillary temperature was 160°C, and 35 units of collision energy were used to obtain peptide fragmentation tandem MS mode (MSMS). To increase the number of identifications, eight μl of different samples were also applied to an Agilent 1100 nanoflow system coupled to a 7-tesla hybrid LTQ Fourier transform ion cyclotron resonance (FTICR) MS (Thermo Electron, Bremen, Germany) equipped with a nanoelectrospray ion source (Proxeon Biosystems, Odense, Denmark) as described previously [43]. Mass spectrometric experiments were performed using unattended data-dependent acquisition in which the mass spectrometer automatically switches between a high resolution survey scan (r 100.000), (accuracy 5 ppm) followed by consecutive ECD and CAD fragmentation (r 25.000), (accuracy 0.02 Da) of the two most abundant multiply charged peptides eluting at this moment from the nano-LC column. Additional experiments were performed where the precursor masses were as previously detected with high mass accuracy (accuracy 0.02 Da), but the fragmentation was performed solely in the LTQ (accuracy 0.5 Da). Highmass accuracy spectra were deisotoped and filtered using the peptide window approach [43]. Detection of phosphorylated peptides was performed by PhosTShunter [44]. Raw MS data associated with the study may be downloaded from ProteomCommons.org Tranche network https://proteomecommons.org/index.jsp using the following hash:

Vui+Q2v4uXpRm2sWc5j8N3aNn+JISmG2RW9gnNMx0hBsse7DGTF3s7e+ci3lpmHUVsTuOKvvQTerSKqe5X3uzYSGbaUAAAAAAAAYWg==

Peptide identification

The raw LTQ-MS data was converted to dta files by Xcalibur 1.4 SR1 and assembled by an in-house developed script to Mascot generic files. These files were searched against downloaded peptide rodent databases from the SwePep endogenous peptide database [45] and a precleaved chicken protein database derived from the NCBI RefSeq database. The pre-cleavage was performed as described by Fälth et al [24]. Estimation of false positives was conducted by searching all spectra against reversed databases [24]. Searches using MASCOT [46] and X! Tandem [47] (included in the Global Proteome Manager, GPM) against the previously mentioned databases were performed while checking for potential modifications (M oxidation, C-terminal amidation, Q/N deamidation, STY phosphorylations and N-terminal acetylation). Fragment mass error tolerance for both X! Tandem and MASCOT was set at 0.7 Da and the parent mass error was set to ± 1.5 Da. Two phosphorylated peptides were detected and characterized with FTICR MS and MASCOT as described above.

LC-MS image analysis

All LTQ-data files were imported into the DeCyder software (version 2.0) for peak detection (PepDetect) and spot matching between samples (PepMatch). PepDetect parameters for peak detection were "Typical peak width: 0.2 min", "Ion-trap mass resolution: 0.7 u", "charge states: 1-10". All files were background corrected using a smooth surface model. Identities from the MASCOT and X! Tandem searches were incorporated into the DeCyder image files. Peaks that were not present after the initial DeCyder spot detection were added manually for all files.

Removal of outliers

MS-data file outliers (see also table 3) were visually detected by the presence of increased protein degradation (Q62), i.e., more peptide peaks, or general low ion intensity (Q10 and Q30), the later by PCA analysis (see below). Removal of sample outliers left 39 MS-data files, of which six were technical replicates (ed12 female control; Q52a-Q52f) to each other, corresponding to 34 unique samples.

Data normalization and analysis

Sample Q62 was removed from the data set before data normalization and analysis due to increased degradation. Ion intensity peptide match data was then exported for further data analysis. The resulting data was log2 transformed. In order to correct for global intensity differences between peptide runs the data was normalized in two steps as previously described [48]. A linear regression was fitted for each individual run to a median run that were constructed of all median peptide values for peptides that were matched in > 50% of all runs. Based on the linear regression equation new values were predicted for each run. A locally-weighted polynomial regression (Lowess) was then fitted for each matched peptide against the run order and the mean value across all runs were added to retain the native intensity dimension. For each matched peptide a proportion of 0.5 neighbors (runs), weighted by their distance to the measurement, were used for controlling the smoothness of the fit.

The data was grouped into eight groups: ed12 male controls (ed12MC, n = 5), ed12 female controls (ed12FC, n = 5), ed12 male EE₂ (ed12ME, n = 5), ed12 female EE₂ (ed12FE, n = 4), ed17 male controls (ed17MC, n = 4)), ed17 female controls (ed17FC, n = 4), ed17 male EE₂ (ed17ME, n = 4), and ed17 female EE₂ (ed17FE, n = 5) (Table 4). A minimum matching cut-off (max one missing file per group) resulted in 204 peptides. All 204 peptides had been manually inspected to control for mismatching. An unsupervised multivariate analysis by principal component analysis (PCA) using the SIMCA-P software (version 12, Umetrics AB, Umeå, Sweden) indicated previously mentioned Q10 (ed12ME) and Q30 (ed12FC) to be outliers which led them to be removed. Characterization of normalized log2 peptide ion intensities was performed by hierarchical clustering in the PermutMatrix program [49] using complete linkage and Euclidian distance as the distance metric. To better visualize age differences, the median value from the log2 ion intensities for each peptide and group (x_class: ed12MC, ed12FC, ed12ME, ed12FE, ed17MC, ed17FC) was calculated and compared to the total average (x_tot) ion intensity for each peptide (x_class-x_tot). A PCA model with two principal components (C1-C2) explained ~42% of the variation in the data (R2X = 0.421; Q2X = 0.298), Eigen values for the components being 11 (C1), 3.33 (C2) (Additional file 1). No patterning was discovered for sex, leaving age and EE₂ (corresponding to four classes in the data; combined male and female ed12C, ed12E, ed17C and ed17E) as parameters for further supervised multivariate analysis. An optimized class separation based on the two variables of age and EE₂ (corresponding to four classes: ed12C, ed12E, ed17C and ed17E) exposure was conducted by means of Partial Least Square Discriminant Analysis (PLS/DA). The age & EE₂ PLS/DA resulted in a model with three components (DC1-3; Eigen values 11, 2.62 and 2.12) explaining ~72% of the variation in the data (R2X = 0.464; Q2X = 0.367 and R2Y = 0.722, Additional file 2). See also additional file 3 for additional information. The most influential peptides for the separation of the four classes were identified by calculating the VIP (variable importance in the projection). Peptides with large VIP, larger than 1, are the most relevant for explaining the separation. Confidence intervals for VIP were derived from jack knifing. Both PCA and PLS/DA models were optimized using Q2 cross validation scores and had the Hotelling's T² tolerance region set at 95% and the DmodX critical distance (Dcrit) set to 0.05. The median log2 ion intensity value for each peptide from the six technical replicates (Q52a-Q52f) was used for both PCA and PLS/DA analysis.

Table 4 Experimental groups and PLS/DA classes for age and EE2 exposure

To test for differences in peptide expression between groups (not including outlier samples Q10, Q30 and Q62) and using the median value of the six technical replicates (Q52a-Q52f)) a linear model was employed estimating the effects of the fixed effects age, sex and treatment. The linear model used was:

where μ was the average peptide expression, A _i represented the age factor (i = 12,17), S _j the sex factor (j = M, F), T _k the treatment factor (k = C, T) and AS _ij , AT _ik , ST _jk and AST _ijk the interaction factors and ε _ijkl was the NID(0, σ²) error component. Marginal sum of squares were calculated for main effects and interaction effects. An F-test was performed to test the significance of the model and the p-values were adjusted for multiple testing [50]. The joint significance of all interaction terms was tested using Wald F-test. Least squares means for statistical significant (p < 0.05) main effects were calculated (see additional file 4). For pair wise t-testing, three levels of statistical significance were used: p < 0.05 (*), p < 0.01 (**) and p < 0.001 (***) (see Additional file 5). The nlme [51] library available in the R [52] software was used for the statistical analysis. Follow up comparative sequence analysis for specific peptides was done using BLAST 2 sequences alignment and ClustalW.