A Strategy for Discovery and Verification of Candidate Biomarkers in Cerebrospinal Fluid of Preclinical Alzheimer’s Disease

Introduction

Alzheimer’s disease (AD) is the most common form of dementia among the elderly. It is characterized by deposition of peptide amyloid-β (Aβ) as amyloid plaques and of protein tau as neurofibrillary tangles (Blennow et al., 2006). Complemented with molecular imaging techniques, cerebrospinal fluid (CSF) markers have been established to reflect AD pathologies (Blennow et al., 2010; Zetterberg, 2015). Emerging clinical cohort studies suggest that subtle cognitive alteration begins years before mild cognitive impairment (Sperling et al., 2011). It is crucial to establish a panel of reliable protein biomarkers for diagnosis of AD at very earliest stage (Hampel et al., 2010). CSF Aβ42 is a well-established marker of Aβ plaque pathology that is highly predictive of AD in cognitively normal individuals (Dubois et al., 2016). Although several peptide and protein biomarkers in CSF have been used for aiding in AD diagnosis (Puchades et al., 2003; Hall et al., 2012), an unambiguous diagnosis of AD at preclinical stage is still lacking. To establish novel and reliable biomarkers capable of early diagnosis in preclinical AD, mass spectrometry (MS) is implemented for unbiased candidate protein biomarker discovery as the principal enabling technology. MS-based label-free quantification holds special promise to discover significant disease-related protein dysregulations for large cohorts with high throughput and sensitivity. The potential biomarkers can be further verified through targeted proteomics with high accuracy and reproducibility (Peterson et al., 2012).

Enzyme-linked immunosorbent assay (ELISA) is widely used for biomarker verification in clinical laboratories, but the dependence on high-quality antibodies restrict its feasibility to quantify novel protein biomarkers (Hüttenhain et al., 2009). Absolute quantification (AQUA) strategy was first proposed to use known concentrations of synthesized heavy isotope-encoded peptides as internal standards to quantify concentrations of native peptides in the clinical sample (Gerber et al., 2003). However, the number of biomarker candidates to be verified would be limited due to the high cost of AQUA. Alternatively, mass difference labeling strategies such as stable isotope labeling by amino acids in cell culture (SILAC) (Ong et al., 2002; Seyfried et al., 2010), combined precursor isotopic labeling and isobaric tagging (cPILOT) (Robinson and Evans, 2012), amino acid-coded tagging (Pan et al., 2003), and mass differential tags for relative and absolute quantification (mTRAQ) (Kang et al., 2010) introduce heavy isotopes metabolically or chemically into target peptides to yield different precursor ion masses, which can be quantified through full MS scans. For example, mTRAQ can react with free N-terminus of a peptide through N-hydroxysuccinimide (NHS) ester group, there is no need to preselect target peptide to incorporate stable isotopes (Ross et al., 2004). mTRAQ is a relatively cost-efficient approach compared to AQUA, but like AQUA (Pannee et al., 2013), the three channels used for mTRAQ limit its ability to construct a standard curve with at most two channels to label standards and the third to label sample, leading to inaccurate estimates when protein amounts span a wide dynamic range in CSF.

The development of a set of five mass-difference reagents, isotopic N,N-dimethyl leucine (iDiLeu) tags, offers greatly reduced cost and significantly improves quantification throughput as compared to AQUA and mTRAQ (Greer et al., 2015). iDiLeu labeling strategy enables construction of a four-point calibration curve to quantify analytes in a single liquid chromatography–mass spectrometry (LC–MS) run. Each of the 5-plex iDiLeu reagents is comprised of an isotopic N,N-dimethyl leucine and an amine-reactive triazine ester moiety, which selectively labels the N-terminus of peptides as well as lysine side chains (Supplementary Figure S1). Mass additions of 141.1154, 144.1313, 147.1409, 150.1631, and 153.1644 Da are efficiently incorporated into peptides by d0, d3, d6, d9, and d12 labels, respectively. Different mass shifts can be discrete in the MS¹ precursor ion scan for quantification, while tandem mass (MS²) yields abundant sequence-specific fragment ions for identification. Four iDiLeu labels are used to set up internal calibration curve by labeling peptide standards with known concentrations and the fifth channel is used to label target peptide of interest in CSF, making absolute quantification achievable based on the intensity of precursor ions in one LC–MS run. Here, we report a comparative global proteomics study in CSF from healthy and preclinical AD individuals to discover candidate biomarkers, which are further validated by targeted quantitative analysis with 5-plex iDiLeu tagging strategy (Figure 1).

FIGURE 1

Figure 1. Workflow for label-free quantification-based discovery proteomics and isotopic labeling-based targeted proteomics.

Materials and Methods

Participants

Twelve enrollees (six individuals in preclinical AD stage and equal number of healthy controls) from Wisconsin Alzheimer’s Disease Research Center (ADRC) participated in this study. The participants comprised of three women in preclinical AD and three gender, family history, and years of education matched healthy controls respectively, as well as three paired men. Table 1 summarizes the demographic characteristics of the participants. Classification of participants as cognitively normal was based on a comprehensive neuropsychological test battery (Almeida et al., 2015). Global cognitive function was evaluated by Mini-Mental State Examination (MMSE), on which a score of ≤24 is deemed an indication of cognition impairment (Folstein et al., 1975). All participants had MMSE scores of ≥28, indicating cognition normal. Classification of preclinical AD was based on evidence of amyloid beta (Aβ) accumulation on ¹¹C Pittsburgh Compound B Positron emission tomography (PET) imaging and hypometabolism on ¹⁸F Fluoro-2-deoxy-glucose (FDG)-PET. The University of Wisconsin Institutional Review Board approved all study procedures. All participants provided signed informed consent.

TABLE 1

Table 1. Characteristics of the current study participants.

CSF Protein Digestion

One mL CSF samples were collected through lumbar puncture. After 3 kDa molecular weight cut off separation, the concentrated protein was subjected to Sigma LC2 IgY column and 14 most abundant proteins were depleted to reduce sample complexity. 95–99% depletion efficiency was achieved and the unbound protein was measured by 660 nm protein assay kit (Thermo Scientific Pierce) at absorbance of 660 nm in Tecan Ultra 384 multi-detection microplate reader (Männedorf, Switzerland). Equal amounts of CSF proteins were reduced by adding dithiothreitol to a final concentration of 5 mM and incubated at room temperature for 1 h. 15 mM of iodoacetamide was used for alkylation of cysteines by incubating for 30 min at room temperature in the dark. After quenching unreacted iodoacetamide with 5 mM dithiothreitol, the protein mixture was diluted with 50 mM Tris buffer (pH 8) to a final urea concentration of 1 M prior to digestion with trypsin at a protein:enzyme ratio of 50:1 at 37°C for 16 h. The digestion reaction was quenched by acidification with 10% trifluoroacetic acid (TFA) to pH 3, followed by desalting with Bond Elut OMIX C18 pipette tips (Agilent Technologies, Santa Clara, CA, United States). For targeted proteomics, another one mL CSF samples were dried down in vacuum via centrifugation with a SpeedVac concentrator (Thermo Scientific, Waltham, MA, United States) and resuspended in 100 μL of lysis buffer, which contained 8 M urea, 50 mM tris base (adjust pH to 8 with hydrochloric acid), 5 mM CaCl₂, 20 mM NaCl, and 1 tablet of EDTA-free protease inhibitor cocktail. CSF proteins were quantified and digested as above.

Synthesis and Activation of iDiLeu Reagents

Synthesis of 5-plex iDiLeu isobaric tags was previously described in detail (Xiang et al., 2010; Greer et al., 2015). iDiLeu tags stored at -20°C are very stable until activation for labeling. The activation buffer consisted of 4-(4,6-dimethoxy-1,3,5-triazin-2-yl)-4-methylmorpholinium tetrafluoroborate (DMTMM) and N-methylmorpholine (NMM) at 0.6x molar ratios of iDiLeu tags dissolved in anhydrous N,N-dimethylformamide (DMF). Activation of iDiLeu tags was conducted by vortexing at room temperature for 1 h. The mixture was centrifuged at 14,000 ×g for 1 min and the supernatant was transferred immediately for peptide labeling.

Peptides Labeling by 5-Plex iDiLeu Tags

A stock solution of synthesized neurosecretory VGF peptide (THLGEALAPLSK) was prepared to be 1 μg/μL. Five aliquots were set up for labeling with d0, d3, d6, d9, d12 iDiLeu reagent, respectively. In each aliquot, 20 μg of peptide standards were resuspended in 20 μL of 0.5 M triethylammonium bicarbonate (TEAB) buffer. To ensure sufficient labeling, peptide standards ware labeled with a 20x w/w excess of iDiLeu reagents. 60 μL of activated iDiLeu reagents were added to each peptide standards aliquot to make the organic:aqueous ratio of 75%. The reaction was vortexing for 2 h under room temperature. 50% hydroxylamine was added to final 0.25% to quench the reaction and incubated for another 10 min. Labeled peptides were vacuum-dried. Two hundred μg of apoE protein standard was digested based on the digestion procedure described above. Each of five aliquots containing 20 μg of apoE protein standard digests was labeled with the five-channel iDiLeu reagents individually.

For correction factor, 2.5 μg of peptide standards were reconstituted in strong cation exchange (SCX) reconstitution buffer. By-products of iDiLeu labeling reaction were removed from SCX SpinTips (Protea Biosciences, Morgantown, WV, United States) according to the manufacturer’s protocol. The eluate was dried in vacuo and desalted with ZipTip C18 pipette tips (Merck Millipore, Darmstadt, Germany).

Cerebrospinal fluid protein digests were reconstituted in 20 μL of 0.5 M TEAB and labeled with d0 reagent separately according to the peptide standards labeling procedure described above. The other four channels labeled peptide standards were spiked into each d0-labeled CSF sample in a ratio of 1:10:50:100. The combined sample was cleaned up with SCX SpinTips and desalted with Bond Elut OMIX C18 pipette tips. All the labeled samples were then dried in vacuo and reconstituted in 25 μL of 3% acetonitrile (ACN), 0.1% formic acid (FA) in water.

LC–MS/MS Acquisition

Samples were analyzed using a Waters nanoAcquity ultra-performance liquid chromatography (UPLC) system (Milford, MA, United States) coupled to a Thermo Scientific Q-Exactive Orbitrap mass spectrometer (San Jose, CA, United States). Peptide mixture was loaded onto a fabricated column with an integrated emitter tip, which packed with Bridged Ethylene Hybrid C18 particles (75 μm × 150 mm, 1.7 μm, 130 Å). Mobile phase A was 0.1% FA in H₂O and mobile phase B was 0.1% FA in ACN. Peptides were separated with a gradient elution of 4 to 35% B over 90 min at a flow rate of 300 nL/min. Full MS scan was acquired in profile mode ranging from m/z 350 to 1800 at a resolution of 70 K. Automatic gain control (AGC) target was 1 × 10⁶, and maximum injection time was 120 ms. Tandem mass spectra were acquired at centroid mode. The top 15 most abundant precursor ions were selected for higher-energy collisional dissociation (HCD) fragmentation with a dynamic exclusion for 40 s with a 10 ppm tolerance. Data-dependent acquisition parameters were set as resolution power of 17.5 K, isolation window of 2.0 Th, normalized collision energy (NCE) of 27, the maximum injection time of 150 ms, AGC target of 1 × 10⁵, and fixed first mass of m/z 110. For targeted proteomics, an inclusion list containing all the precursor ions of peptide of interest was constructed. Each sample was acquired in technical triplicates.

Data Analysis

In discovery proteomics analysis, protein identification was performed using MaxQuant (1.5.6.5) against SwissProt human database with 1% false discovery rate (FDR) at protein level. The first search peptide tolerance for precursor and product ion were 20 ppm and 0.02 Da, respectively. The maximum missed cleavages per peptide was 2. Fixed modification was set as carbamidomethylation of cysteine residues (+57.0215 Da). Oxidation of methionine (+15.9949 Da) was selected as variable modifications. The intensity obtained from label-free quantification in MaxQuant was uploaded into Perseus for advanced downstream analysis. Firstly, the reverse protein sequences created during database search, the most common contaminants in the MS analyses, and proteins identified only by modified sequences were removed. Secondly, proteins identified based on only one peptide and proteins for which no unique peptide was identified were removed. To ensure accurate quantification, proteins identified in each sample along three technical replicates were normalized to the summed intensity, followed by determination of the ratio between preclinical AD and healthy control as well as statistical analysis of Student’s t-test. Gene ontology (GO) analysis of the significantly regulated proteins in preclinical AD patients was performed by DAVID, which is a tool to provide enrichment analysis of the most relevant GO terms associated with a given list (Huang et al., 2009).

In targeted proteomics analysis, protein identification was performed on Peaks Studio 7 software (Bioinformatics Solutions, Inc., Waterloo, ON, Canada). The data refinement was applied to correct precursor mass by default. All the raw files were searched against UniProt Homo sapiens reviewed database with trypsin as digestion enzyme. The error tolerance for precursor mass was 25 ppm using monoisotopic mass and 0.02 Da for fragment ion. The maximum missed cleavages per peptide was two, allowed to be cleaved at both ends of the peptides. Fixed modification was set as carbamidomethylation of cysteine residues (+57.0215 Da). iDiLeu labels (+141.1154 Da for d0, +144.1313 Da for d3, +147.1409 Da for d6, +150.1631 Da for d9, and +153.1644 Da for d12) of peptide N-termini and lysine residues, as well as oxidation of methionine (+15.9949 Da) were selected as variable modifications. Estimation of FDR was enabled. Peptides with FDR < 1% were considered to be unambiguous identification. Peak areas generated by Genesis peak detection algorithm in Thermo Xcalibur 2.2 software were used for quantification. The precursor ion integration tolerance was 15 ppm. Retention time of extracted ion chromatogram of 5-plex iDiLeu-labeled peptides was required to be within 2 min. Isotopic interference correction factors were applied to each sample to rectify the raw values.

Results and Discussion

Label-Free Quantification for Discovery Proteomics

To discover candidate protein biomarkers for early diagnosis of AD in preclinical stage, label-free based discovery proteomics displayed in Figure 1 was employed for measurement of protein changes between healthy and preclinical AD individuals as a faster and cleaner technique. Given sex differences in the risk of AD (Mielke et al., 2014), the study was designed so that gender could be taken into account in all the analyses. As shown in Figure 2A, 708 and 690 proteins were identified in healthy and preclinical AD for women with 1% protein-level FDR, respectively. The identified proteins in men were comparable to women (659 and 653 proteins in each group, Figure 2B). Taken collectively, 732 and 704 proteins with more than one unique peptide were identified in healthy controls and preclinical AD patients.

FIGURE 2

Figure 2. Venn diagram illustration of proteins identified in healthy controls (gray) and preclinical AD patients (orange) for women (A) and men (B), respectively. Volcano plots of quantified proteins in women (C) and men (D).

All the identified proteins were submitted to Perseus for advanced downstream analysis (Tyanova et al., 2016). To ensure accurate quantification, proteins identified based on only one peptide were removed. Proteins identified in sex-matched three pairs of healthy and preclinical AD subjects along with all the technical replicates were utilized for label-free quantification. As shown in volcano plots (Figures 2C,D), 79 and 98 proteins were significantly altered in preclinical AD for women and men in comparison to corresponding healthy controls, respectively (Supplementary Tables S1A,B). They were sequentially used for gene ontology (GO) enrichment analysis to understand what biological processes were statistically over-represented. The top 10 biological processes with Fisher Exact p-value < 0.05 in preclinical AD for women and men were shown in Figure 3. Negative regulation of endopeptidase activity, platelet degranulation, and cell adhesion were the three most prominent biological processes in both women and men, all of which were related to AD pathology. Many endopeptidases are capable of proteolyzing Aβ, especially neprilysin, which exhibits the most potent Aβ-degrading activity (Shirotani et al., 2001). Down-regulated neprilysin in the hippocampus and cerebral cortex has been reported in AD development (Shirotani et al., 2001). Platelet degranulation played an important role in platelet-mediated amyloid beta (Aβ) oligomerization (Donner et al., 2016). It was suggested that Aβ stimulated aberrant neuritic growth by activation of cell adhesion signaling pathways (Cotman et al., 1998; Grace et al., 2002).

FIGURE 3

Figure 3. GO analysis of significantly regulated proteins in women (A) and men (B) with preclinical AD. GO term with Fisher Exact p-value smaller than 0.05 was considered to be a strongly enriched biological process.

Many of the differentially expressed proteins have been demonstrated to have similar up- or down-regulation patterns in earlier studies of MCI or AD dementia, including elevated amyloid beta precursor protein in women with preclinical AD (p-value = 0.006), which could be further proteolytically cleaved into neurotoxic Aβ42 peptide (Zlokovic, 2005; Sleegers et al., 2006). Secretogranin II, a precursor for neuropeptide secretoneurin, widely distributed in endocrine and nervous tissues, was decreased (p-value = 0.00003) in men with preclinical AD as compared to healthy subjects (Li et al., 2008). Interestingly, some proteins displayed opposite regulation trends between women and men in preclinical AD. Neurosecretory protein VGF (VGF), synthesized by neurons where it promotes dendritic growth and survival of cortical neurons (Sato et al., 2012), was significantly down-regulated in men with preclinical AD but increased in women. Apolipoprotein E (apoE) has three common isoforms, all of which differentially modulate Aβ aggregation and clearance (Bu, 2009). ApoE was up-regulated in women with preclinical AD patients but slightly decreased in men. Biomarker candidates of VGF and apoE determined from discovery proteomics need to be further validated through isotopic labeling-based targeted proteomics.

Quantitative Performance of 5-Plex iDiLeu Labeling Strategy

To demonstrate the strong performance of 5-plex iDiLeu labeling strategy, a mixture of 5-plex iDiLeu-labeled peptide standards were employed for characterization and quantification. The starting material, L-leucine or isotopic L-leucine (L-leucine-1-¹³C, ¹⁵N) used for iDiLeu synthesis are commercially available, but the slight impurities of these reagents can diminish quantification accuracy. Furthermore, when 5-plex iDiLeu reagents only label the N-terminus of a peptide with +2 charge state, the m/z of five monoisotopic precursors will be 1.5 Th apart, which will inevitably introduce isotopic interference. Therefore, correction factors need to be firstly determined for each of the 5-plex iDiLeu tags to ensure accurate quantification (Supplementary Table S2). A series of equations based on i-Tracker were set up to rectify raw values with correction factors (Shadforth et al., 2005; Greer et al., 2015). Each channel has its own correction equation, where the peak area of both monoisotopic precursor and any possible interferences from neighboring isotopic precursor were taken into consideration.

Peptide standards were labeled with 5-plex iDiLeu reagents. The combination of 5-plex iDiLeu labeled peptides was subjected to LC–MS/MS analysis in triplicates at the ratio of 1:1:1:1:1. Extracted ion chromatograms of precursor ions were used for quantification of 5-plex iDiLeu-labeled peptides. The median ratio of d0:d3:d6:d9:d12 was 0.89:1.08:0.87:0.85:1.00 as shown in the boxplots (Supplementary Figure S2). Medians displayed high accuracy with <15% error, which confirms quantification based on MS¹-precursor level was reliable in high-resolution and accurate mass (HRAM) measurement by Q-Exactive Orbitrap mass spectrometer.

To determine the linear concentration range of this iDiLeu labeling strategy, standards of VGF peptide and tryptic apoE peptides were individually labeled with 5-plex iDiLeu tags and combined at the ratio of 1:10:50:100. The highest concentration stock of VGF peptide ranged from 30 to 3000 fmol/μL. This stock solution was serially diluted twice to prepare concentrations of 3–300 and 0.3–30 fmol/μL. Figure 4 displayed peak area of iDiLeu-labeled VGF peptide at the MS¹ level as a function of concentration across LC–MS/MS analysis in triplicates by normalization to d9 channel. Each calibration curve revealed a linear response comprising R² around 0.999 spanning two orders of magnitude. The higher the concentration, the lower the coefficient of variation (CV). The CV was below 0.8% in the concentration of 30–3000 fmol/μL but increased to be 8% among the triplicates in 0.3–30 fmol/μL. The limits of quantification for VGF peptide and tryptic apoE peptides were all around 0.1 fmol/μL. The high linearity and low limit of quantification confirmed that iDiLeu tagging is a robust strategy for absolute quantification.

FIGURE 4

Figure 4. Calibration curves of 5-plex iDiLeu labeled VGF peptide in triplicates. The normalized peak area of precursor ions plotted as a function of peptide concentrations at 0.3–30, 3–300, 30–3000 fmol/μL.

Absolute Quantification of Biomarkers in Preclinical Alzheimer’s Disease

The iDiLeu labeling strategy has been evaluated with high accuracy and linearity, allowing for verification of biomarker candidates in CSF samples. We set up the targeted proteomics strategy to perform quantification in both healthy and preclinical AD subjects using 5-plex iDiLeu reagents (Figure 1). Each CSF protein digests was labeled with d0 channel, respectively. Peptide standards were composed of peptides derived from proteins VGF and apoE. The standards were labeled with d3, d6, d9, d12, individually, followed by combination at the ratio of 1:10:50:100. The four channels labeled peptide standards were spiked into each CSF sample with equal amount at the same ratio. The final concentration of peptide standard was ranged from 0.04 to 4 pmol/μL for proteins VGF and apoE.

iDiLeu labeled peptides were fragmented under NCE of 27. Figure 5 displayed representative tandem mass spectrum of d0-labeled peptide AATVGSLAGQPLQER in CSF sample. iDiLeu d0 tag was fragmented to produce dimethylated immonium reporter ion at m/z of 114.1. Abundant b- and y- product ions have been matched to specific peptide sequence within an FDR < 1%, which is equivalent to a -10lgP score of 15. This score is derived from the p-value that indicates the statistical significance of the peptide-spectrum match (Zhang et al., 2012). Peptide AATVGSLAGQPLQER of apoE had -10lgP score of 58.13, which presented as a high-quality spectrum with modification of iDiLeu tag. The plot of normalized peak area of each peptide against concentration displayed high linearity in both preclinical AD and control samples, which ensured the accuracy of biomarker quantification (Figure 6). A standard deviation within 3% was demonstrated across the triplicates run in each sample.

FIGURE 5

Figure 5. Example of tandem mass spectrum for iDiLeu d0-labeled apoE peptide in CSF.

FIGURE 6

Figure 6. Calibration curves of proteins VGF and apoE constructed in both healthy (A,B) and preclinical AD (C,D) subjects.

The amounts of VGF were quantified as 55 and 21 ng/mL in CSF of women and men with preclinical AD, respectively. ApoE has higher abundance than VGF, which were 1.38 and 0.9 μg/mL in preclinical AD for women and men, respectively (Figure 7). The isotopic labeling based targeted proteomics suggested that both VGF and apoE were down-regulated in men with preclinical AD but showing opposite regulation pattern in women, which is consistent with that of label-free based discovery proteomics. The average ages for women were younger than men, which may partially contribute to the opposite phenomenon as aging is the greatest risk factor for AD. Moreover, women had significantly higher glucose metabolism in the brain areas involved in the pathological process of AD, while men showed consistent level of cognitive impairment (Perneczky et al., 2007). Sex difference in cortical cytoarchitecture, such as brain size, neuron account, and neuron density, may also lead to brain cognitive reserve difference in women and men (Rabinowicz et al., 2002).

FIGURE 7

Figure 7. The determination of absolute amounts of proteins VGF and apoE in CSF of women (A,B) and men (C,D).

There is no clear consensus regarding whether VGF and apoE are positively or negatively associated with AD dementia (Hesse et al., 2000; Carrette et al., 2003; Liu et al., 2013; Busse et al., 2015; Hölttä et al., 2015). Statistical analysis demonstrated that VGF was significantly down-regulated in men with preclinical AD, possibly due to its regulation of hippocampal function and memory, in part through modulation of brain-derived neurotrophic factor (BDNF) secretion and signaling (Bozdagi et al., 2008). The increased expression of VGF in women may be correlated to immune manifestations with negative effects on brain function and neuronal repair processes, as reported that VGF was expressed with a higher percentage of peripheral CD3+ T cells from AD patients compared to age-matched healthy controls (Akiyama et al., 2000; Busse et al., 2015).

ApoE is a lipoprotein containing 299 amino acids, regulates lipid homeostasis by mediating lipid transport between various cells in the brain. The delivering function of Aβ peptide is thought to initiate toxic events, leading to synaptic dysfunction and neurodegeneration in AD (Liu et al., 2013). Analysis of apoE alleles in AD demonstrated that there was a significant association of 𝜀4 allele with an earlier age of AD onset, resulting from domain interactions in the presence of arginine at residue 112 (Kim et al., 2009; Wolf et al., 2013). ApoE may function as Aβ chaperones for receptor-mediated clearance of soluble extracellular Aβ but may also lead to intraneuronal Aβ accumulation (Lane and Farlow, 2005). The absolute amount and particular isoform of apoE may influence both Aβ fibrillogenesis and metabolism to differing extents. It has been reported that increased cytokines and decreased cholesterol in the brain could reduce apoE levels in AD (Gueguen et al., 2001). However, higher levels of apoE have also been reported in CSF of AD patients compared with age-matched controls (Lindh et al., 1997; Sihlbom et al., 2008). According to our results, 1.38 and 0.9 μg/mL of apoE could be defined as thresholds for women and men with preclinical AD, respectively, but a large cohort study would be necessary to determine the threshold of multiple biomarkers to aid pre-clinical diagnosis.

Conclusion

In summary, we reported a large number of proteins identified in human CSF. A panel of potential protein biomarkers was produced to differentiate healthy controls and preclinical AD patients in this pilot study. Many of the candidate biomarkers have shown similar trends in patients with MCI or AD, including amyloid precursor protein, apolipoprotein E, and several neuroendocrine proteins. More interestingly, the differential expression of proteins in women and men with preclinical AD might provide us with a novel perspective to reveal the pathological mechanism of AD, but a larger cohort study is necessary to further elucidate the correlation between AD and gender differences. For the first time, the novel iDiLeu-enabled standard curve approach has been established for absolute quantification of candidate protein biomarkers in each CSF sample with high accuracy and reproducibility. Taken together, we established a comprehensive strategy to discover and verify candidate biomarkers by integrating label-free based discovery proteomics with isotopic labeling based targeted proteomics. This strategy could be applied for biomarker discovery and verification across a broad spectrum of diseases.

Author Contributions

XZ, JW, and LL designed the research. XZ and JW performed the experiments and analyzed the data. HZ, CC, and OO collected and provided the CSF specimens. XZ, HZ, and LL prepared the manuscript and all authors provided feedback.

Funding

This research was protectsupported in part by the National Institutes of Health (NIH) (Grant Nos. R01AG052324, 1P50AG033514, R01DK071801, and P41GM108538). Support for this research was also provided by the University of Wisconsin-Madison, Office of the Vice Chancellor for Research and Graduate Education with funding from the Wisconsin Alumni Research Foundation. The Orbitrap instruments were purchased through the support of an NIH shared instrument grant (NIH-NCRR S10RR029531) and Office of the Vice Chancellor for Research and Graduate Education at the University of Wisconsin-Madison. LL acknowledges a Vilas Distinguished Achievement Professorship and the Charles Melbourne Johnson Distinguished Chair Professorship with funding provided by the Wisconsin Alumni Research Foundation and University of Wisconsin-Madison School of Pharmacy. HZ was a Wallenberg Academy Fellow.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2018.00483/full#supplementary-material

References

Akiyama, H., Barger, S., Barnum, S., Bradt, B., Bauer, J., Cole, G. M., et al. (2000). Inflammation and Alzheimer’s disease. Neurobiol. Aging 21, 383–421. doi: 10.1016/S0197-4580(00)00124-X

CrossRef Full Text | Google Scholar

Almeida, R. P., Schultz, S. A., Austin, B. P., Boots, E. A., Dowling, N. M., Gleason, C. E., et al. (2015). Effect of cognitive reserve on age-related changes in cerebrospinal fluid biomarkers of Alzheimer disease. JAMA Neurol. 72, 699–706. doi: 10.1001/jamaneurol.2015.0098

PubMed Abstract | CrossRef Full Text | Google Scholar

Blennow, K., de Leon, M. J., and Zetterberg, H. (2006). Alzheimer’s disease. Lancet 368, 387–403. doi: 10.1016/S0140-6736(06)69113-7

CrossRef Full Text | Google Scholar

Blennow, K., Hampel, H., Weiner, M., and Zetterberg, H. (2010). Cerebrospinal fluid and plasma biomarkers in Alzheimer disease. Nat. Rev. Neurol. 6, 131–144. doi: 10.1038/nrneurol.2010.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bozdagi, O., Rich, E., Tronel, S., Sadahiro, M., Patterson, K., Shapiro, M. L., et al. (2008). The neurotrophin-inducible gene vgf regulates hippocampal function and behavior through a brain-derived neurotrophic factor-dependent mechanism. J. Neurosci. 28, 9857–9869. doi: 10.1523/JNEUROSCI.3145-08.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bu, G. (2009). Apolipoprotein e and its receptors in Alzheimer’s disease: pathways, pathogenesis and therapy. Nat. Rev. Neurosci. 10, 333–344. doi: 10.1038/nrn2620

PubMed Abstract | CrossRef Full Text | Google Scholar

Busse, S., Steiner, J., Glorius, S., Dobrowolny, H., Greiner-Bohl, S., Mawrin, C., et al. (2015). VGF expression by T lymphocytes in patients with Alzheimer’s disease. Oncotarget 6, 14843–14851. doi: 10.18632/oncotarget.3569

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrette, O., Demalte, I., Scherl, A., Yalkinoglu, O., Corthals, G., Burkhard, P., et al. (2003). A panel of cerebrospinal fluid potential biomarkers for the diagnosis of Alzheimer’s disease. Proteomics 3, 1486–1494. doi: 10.1002/pmic.200300470

PubMed Abstract | CrossRef Full Text | Google Scholar

Cotman, C. W., Hailer, N. P., Pfister, K. K., Soltesz, I., and Schachner, M. (1998). Cell adhesion molecules in neural plasticity and pathology: similar mechanisms, distinct organizations? Prog. Neurobiol. 55, 659–669. doi: 10.1016/S0301-0082(98)00025-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Donner, L., Fälker, K., Gremer, L., Klinker, S., Pagani, G., Ljungberg, L. U., et al. (2016). Platelets contribute to amyloid-b aggregation in cerebral vessels through integrin αIIbβ3-induced outside-in signaling and clusterin release. Sci. Signal. 9, 1–17. doi: 10.1126/scisignal.aaf6240

CrossRef Full Text | Google Scholar

Dubois, B., Hampel, H., Feldman, H. H., Scheltens, P., Aisen, P., Andrieu, S., et al. (2016). Preclinical Alzheimer’s disease: definition, natural history, and diagnostic criteria. Alzheimers Dement. 12, 292–323. doi: 10.1016/j.jalz.2016.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Folstein, M. F., Folstein, S. E., and McHugh, P. R. (1975). “Mini-mental state”: a practical method for grading the cognitive state of patients for the clinician. J. Psychiatr. Res. 12, 189–198.

Google Scholar

Gerber, S. A., Rush, J., Stemman, O., Kirschner, M. W., and Gygi, S. P. (2003). Absolute quantification of proteins and phosphoproteins from cell lysates by tandem MS. Proc. Natl. Acad. Sci. U.S.A. 100, 6940–6945. doi: 10.1073/pnas.0832254100

PubMed Abstract | CrossRef Full Text | Google Scholar

Grace, E. A., Rabiner, C. A., and Busciglio, J. (2002). Characterization of neuronal dystrophy induced by fibrillar amyloid β: implications for Alzheimer’s disease. Neuroscience 114, 265–273. doi: 10.1016/S0306-4522(02)00241-5

CrossRef Full Text | Google Scholar

Greer, T., Lietz, C. B., Xiang, F., and Li, L. (2015). Novel isotopic N,N-dimethyl leucine (iDiLeu) reagents enable absolute quantification of peptides and proteins using a standard curve approach. J. Am. Soc. Mass Spectrom. 26, 107–119. doi: 10.1007/s13361-014-1012-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Gueguen, Y., Bertrand, P., Ferrari, L., Batt, A. M., and Siest, G. (2001). Control of apolipoprotein E secretion by 25-hydroxycholesterol and proinflammatory cytokines in the human astrocytoma cell line CCF-STTG1. Cell Biol. Toxicol. 17, 191–199. doi: 10.1023/A:1011996515328

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, S., Öhrfelt, A., Constantinescu, R., Andreasson, U., Surova, Y., Bostrom, F., et al. (2012). Accuracy of a panel of 5 cerebrospinal fluid biomarkers in the differential diagnosis of patients with dementia and/or Parkinsonian disorders. Arch. Neurol. 69, 1445–1452. doi: 10.1001/archneurol.2012.1654

PubMed Abstract | CrossRef Full Text | Google Scholar

Hampel, H., Frank, R., Broich, K., Teipel, S. J., Katz, R. G., Hardy, J., et al. (2010). Biomarkers for alzheimer’s disease: academic, industry and regulatory perspectives. Nat. Rev. Drug Discov. 9, 560–574. doi: 10.1038/nrd3115

PubMed Abstract | CrossRef Full Text | Google Scholar

Hesse, C., Larsson, H., Fredman, P., Minthon, L., Andreasen, N., Davidsson, P., et al. (2000). Measurement of apolipoprotein E (apoE) in cerebrospinal fluid. Neurochem. Res. 25, 511–517. doi: 10.1023/A:1007516210548

CrossRef Full Text | Google Scholar

Hölttä, M., Minthon, L., Hansson, O., Holmén-Larsson, J., Pike, I., Ward, M., et al. (2015). An integrated workflow for multiplex CSF proteomics and peptidomics-identification of candidate cerebrospinal fluid biomarkers of alzheimer’s disease. J. Proteome Res. 14, 654–663. doi: 10.1021/pr501076j

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Hüttenhain, R., Malmström, J., Picotti, P., and Aebersold, R. (2009). Perspectives of targeted mass spectrometry for protein biomarker verification. Curr. Opin. Chem. Biol. 13, 518–525. doi: 10.1016/j.cbpa.2009.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, U. B., Yeom, J., Kim, H., and Lee, C. (2010). Quantitative analysis of mTRAQ-labeled proteome using full MS scans. J. Proteome Res. 9, 3750–3758. doi: 10.1021/pr9011014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J., Basak, J. M., and Holtzman, D. M. (2009). The role of apolipoprotein E in alzheimer’s disease. Neuron 63, 287–303. doi: 10.1016/j.neuron.2009.06.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Lane, R. M., and Farlow, M. R. (2005). Lipid homeostasis and apolipoprotein E in the development and progression of Alzheimer’s disease. J. Lipid Res. 46, 949–968. doi: 10.1194/jlr.M400486-JLR200

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Hung, A. C., and Porter, A. G. (2008). Secretogranin II: a key AP-1-regulated protein that mediates neuronal differentiation and protection from nitric oxide-induced apoptosis of neuroblastoma cells. Cell Death Differ. 15, 879–888. doi: 10.1038/cdd.2008.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindh, M., Blomberg, M., Jensen, M., Basun, H., Lannfelt, L., Engvall, B., et al. (1997). Cerebrospinal fluid apolipoprotein E (apoE) levels in Alzheimer’s disease patients are increased at follow up and show a correlation with levels of tau protein. Neurosci. Lett. 229, 85–88. doi: 10.1016/S0304-3940(97)00429-1

CrossRef Full Text | Google Scholar

Liu, C.-C., Kanekiyo, T., Xu, H., and Bu, G. (2013). Apolipoprotein E and Alzheimer disease: risk, mechanisms and therapy. Nat. Rev. Neurol. 9:106. doi: 10.1038/nrneurol.2012.263

PubMed Abstract | CrossRef Full Text | Google Scholar

Mielke, M. M., Vemuri, P., and Rocca, W. A. (2014). Clinical epidemiology of Alzheimer’s disease: assessing sex and gender differences. Clin. Epidemiol. 6, 37–48. doi: 10.2147/CLEP.S37929

PubMed Abstract | CrossRef Full Text | Google Scholar

Ong, S.-E., Blagoev, B., Kratchmarova, I., Kristensen, D. B., Steen, H., Pandey, A., et al. (2002). Stable isotope labeling by amino acids in cell culture, SILAC, as a simple and accurate approach to expression proteomics. Mol. Cell. Proteomics 1, 376–386. doi: 10.1074/mcp.M200025-MCP200

CrossRef Full Text | Google Scholar

Pan, S., Gu, S., Bradbury, E. M., and Chen, X. (2003). Single peptide-based protein identification in human proteome through MALDI-TOF MS coupled with amino acids coded mass tagging. Anal. Chem. 75, 1316–1324. doi: 10.1021/ac020482s

PubMed Abstract | CrossRef Full Text | Google Scholar

Pannee, J., Portelius, E., Oppermann, M., Atkins, A., Hornshaw, M., Zegers, I., et al. (2013). A selected reaction monitoring (SRM)-based method for absolute quantification of Abeta38, Abeta40, and Abeta42 in cerebrospinal fluid of Alzheimer’s disease patients and healthy controls. J. Alzheimers Dis. 33, 1021–1032. doi: 10.3233/JAD-2012-121471

PubMed Abstract | CrossRef Full Text | Google Scholar

Perneczky, R., Drzezga, A., Diehl-Schmid, J., Li, Y., and Kurz, A. (2007). Gender differences in brain reserve. J. Neurol. 254, 1395. doi: 10.1007/s00415-007-0558-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Peterson, A. C., Russell, J. D., Bailey, D. J., Westphall, M. S., and Coon, J. J. (2012). Parallel reaction monitoring for high resolution and high mass accuracy quantitative, targeted proteomics. Mol. Cell. Proteomics 11, 1475–1488. doi: 10.1074/mcp.O112.020131

PubMed Abstract | CrossRef Full Text | Google Scholar

Puchades, M., Hansson, S. F., Nilsson, C. L., Andreasen, N., Blennow, K., and Davidsson, P. (2003). Proteomic studies of potential cerebrospinal fluid protein markers for Alzheimer’s disease. Mol. Brain Res. 118, 140–146. doi: 10.1016/j.molbrainres.2003.08.005

CrossRef Full Text | Google Scholar

Rabinowicz, T., Petetot, J. M.-C., Gartside, P. S., Sheyn, D., Sheyn, T., and De, C.-M. (2002). Structure of the cerebral cortex in men and women. J. Neuropathol. Exp. Neurol. 61, 46–57. doi: 10.1093/jnen/61.1.46

CrossRef Full Text | Google Scholar

Robinson, R. A. S., and Evans, A. R. (2012). Enhanced sample multiplexing for nitrotyrosine-modified proteins using combined precursor isotopic labeling and isobaric tagging. Anal. Chem. 84, 4677–4686. doi: 10.1021/ac202000v

PubMed Abstract | CrossRef Full Text | Google Scholar

Ross, P. L., Huang, Y. N., Marchese, J. N., Williamson, B., Parker, K., Hattan, S., et al. (2004). Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol. Cell. Proteomics 3, 1154–1169. doi: 10.1074/mcp.M400129-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, H., Fukutani, Y., Yamamoto, Y., Tatara, E., Takemoto, M., Shimamura, K., et al. (2012). Thalamus-derived molecules promote survival and dendritic growth of developing cortical neurons. J. Neurosci. 32, 15388–15402. doi: 10.1523/JNEUROSCI.0293-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Seyfried, N. T., Gozal, Y. M., Dammer, E. B., Xia, Q., Duong, D. M., Cheng, D., et al. (2010). Multiplex SILAC analysis of a cellular TDP-43 proteinopathy model reveals protein inclusions associated with SUMOylation and diverse polyubiquitin chains. Mol. Cell. Proteomics 9, 705–718. doi: 10.1074/mcp.M800390-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Shadforth, I. P., Dunkley, T. P., Lilley, K. S., and Bessant, C. (2005). i-Tracker: for quantitative proteomics using iTRAQ. BMC Genomics 6:145. doi: 10.1186/1471-2164-6-145

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirotani, K., Tsubuki, S., Iwata, N., Takaki, Y., Harigaya, W., Maruyama, K., et al. (2001). Neprilysin degrades both amyloid β peptides 1-40 and 1-42 most rapidly and efficiently among thiorphan- and phosphoramidon-sensitive endopeptidases. J. Biol. Chem. 276, 21895–21901. doi: 10.1074/jbc.M008511200

PubMed Abstract | CrossRef Full Text | Google Scholar

Sihlbom, C., Davidsson, P., Sjögren, M., Wahlund, L.-O., and Nilsson, C. L. (2008). Structural and quantitative comparison of cerebrospinal fluid glycoproteins in Alzheimer’s disease patients and healthy individuals. Neurochem. Res. 33, 1332–1340. doi: 10.1007/s11064-008-9588-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sleegers, K., Brouwers, N., Gijselinck, I., Theuns, J., Goossens, D., Wauters, J., et al. (2006). APP duplication is sufficient to cause early onset Alzheimer’s dementia with cerebral amyloid angiopathy. Brain 129, 2977–2983. doi: 10.1093/brain/awl203

PubMed Abstract | CrossRef Full Text | Google Scholar

Sperling, R. A., Aisen, P. S., Beckett, L. A., Bennett, D. A., Craft, S., Fagan, A. M., et al. (2011). Toward defining the preclinical stages of Alzheimer’s disease: recommendations from the national institute on aging-alzheimer’s association workgroups on diagnostic guidelines for alzheimer’s disease. Alzheimers Dement. 7, 280–292. doi: 10.1016/j.jalz.2011.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyanova, S., Temu, T., Sinitcyn, P., Carlson, A., Hein, M. Y., Geiger, T., et al. (2016). The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 13, 731–740. doi: 10.1038/nmeth.3901

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolf, A. B., Caselli, R. J., Reiman, E. M., and Valla, J. (2013). APOE and neuroenergetics: an emerging paradigm in Alzheimer’s disease. Neurobiol. Aging 34, 1007–1017. doi: 10.1016/j.neurobiolaging.2012.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, F., Ye, H., Chen, R., Fu, Q., and Li, L. (2010). N,N-Dimethyl leucines as novel isobaric tandem mass tags for quantitative proteomics and peptidomics. Anal. Chem. 82, 2817–2825. doi: 10.1021/ac902778d

PubMed Abstract | CrossRef Full Text | Google Scholar

Zetterberg, H. (2015). Cerebrospinal fluid biomarkers for Alzheimer’s disease: current limitations and recent developments. Curr. Opin. Psychiatry 28, 402–409. doi: 10.1097/YCO.0000000000000179

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Xin, L., Shan, B., Chen, W., Xie, M., Yuen, D., et al. (2012). PEAKS DB: de novo sequencing assisted database search for sensitive and accurate peptide identification. Mol. Cell. Proteomics 11, M111.010587. doi: 10.1074/mcp.M111.010587

PubMed Abstract | CrossRef Full Text | Google Scholar

Zlokovic, B. V. (2005). Neurovascular mechanisms of Alzheimer’s neurodegeneration. Trends Neurosci. 28, 202–208. doi: 10.1016/j.tins.2005.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar