key: cord-0807213-yyie09ib authors: Li, Yang; Lai, Dan-yun; Tao, Sheng-ce title: SARS-CoV-2 Spike linear epitope scanning via a peptide microarray through sera profiling date: 2021-07-10 journal: STAR Protoc DOI: 10.1016/j.xpro.2021.100707 sha: afc5077b699d0c9fc340b24e9e57325e1e0f68f6 doc_id: 807213 cord_uid: yyie09ib Host humoral immunological response plays an essential role in protection against pathogens. Identification- of B cell epitopes on antigens is required for accurate diagnosis and vaccine development. To map SARS-CoV-2 Spike linear epitopes, we developed a protocol of profiling sera from COVID-19 patients via a peptide microarray designed according to Spike protein. The protocol is also applicable for other antigens or sample types. This protocol is rapid, high-throughput and the cost is acceptable while it needs specialized microarray facilities. Lead Contact *Corresponding Author, e-mail: taosc@sjtu.edu.cn (S.-C. Tao) Summary Host humoral immunological response plays an essential role in protection against pathogens. Identification-of B cell epitopes on antigens is required for accurate diagnosis and vaccine development. To map SARS-CoV-2 Spike linear epitopes, we developed a protocol of profiling sera from COVID-19 patients via a peptide microarray designed according to Spike protein. The protocol is also applicable for other antigens or sample types. This protocol is rapid, high-throughput and the cost is acceptable while it needs specialized microarray facilities. For complete details on the use and execution of this protocol, please refer to (Li et al., 2020 (Li et al., , 2021a (Li et al., , 2021b . Note: 1) 50 μL positive control are good for around 200 slides, i. e., about 2,800 (200x14) sera samples. 2) Mixing of a set of samples instead of selecting only one as a positive control is to avoid bias. Microarray layout design and positive probe preparation Microarray layout design is to assign a position in an array for each probe and would be used to set printing parameters. The probes include not only the peptides of target antigen(s) but also some negative and positive controls which are used to justify the results of the microarray experiment. Incubation buffer Dissolve 0.1 g BSA into 10 mL Wash buffer or Mix 1 mL 10% (wt/vol) BSA (stored at -20°C) into 9 mL Wash buffer. This solution should be freshly prepared before use. Printing buffer J o u r n a l P r e -p r o o f Mix 2.5 mL glycerol and 1 mL 10 x PBS buffer into 6.5 mL dH2O to a final concentration of 25% (vol/vol) glycerol. This solution can be stored at 4 °C for at least 6 months. Step-by-Step Method Details Conjugation of peptides with BSA To efficiently immobilize the peptide on the substrate slide, the peptides are conjugated to a large carrier protein. Here we used BSA as the carrier. Other protein can be used. However, the protein may bring non-specific binding so as to generate high background or false positive signals. One strategy to lower this effect is that the carrier protein is added into the Blocking buffer and serum Incubation buffer. 1. Reconstitute the lyophilized peptides with double distilled water (ddH2O) to 4 mg/mL. Notably, for hydrophobic peptides which is difficult to be dissolved with water, DMSO (20 μL) could be added to 1 mg peptide followed by adding 230 μL ddH2O. However, a few peptides that still cannot be completely dissolved would have a compromised concentration. 2. BSA activation by Sulfo-SMCC. a. Dissolve 10 mg BSA with 2.5 mL 1 x PBS (phosphate buffer saline). b. Dissolve 2 mg Sulfo-SMCC with 100 μL ddH2O and pipette up and down to mix to generate homogenous solution. c. Add 97.7 μL dissolved Sulfo-SMCC into the BSA solution, gently shake and incubate for 30 min at room temperature (22-27℃). The molar ratio of BSA to Sulfo-SMCC is 1: 30. d. Dialyze in 1 x PBS at 4℃ to remove the excess Sulfo-SMCC. The dialysis buffer (1 x PBS) should be changed at least twice after 3-4 h of dialysis each. 3. Peptide conjugation. a. Add 250 μL activated BSA to 1 mg reconstituted peptides (1:1 in w/w). b. Incubate at room temperature (22-27℃) for 4 h on a rotator. Optional: Conjugation efficiency can be estimated by electrophoresis separation (Figure 1A ). It is highly recommended particularly for those peptides which do not dissolve completely. High quality microarray is essential for high quality data acquisition. This part provides detailed procedures as well as considerations to fabricate peptide microarray. For the Spike peptides, we created an array of 32 x 20 features with triplicates for each probe. 4. Load the dissolved peptides and controls to a 384-well plate a. Dilute the BSA-peptide conjugates with Printing buffer (1 x PBS with 25% glycerol) to 0.3 mg/mL. Note: The concentration refers to the peptides. The optimal concentration of peptides for serum detection is 0.1-1 mg/mL according to our experience. However, the optimal concentration that largely depends on the sample types or antibody concentrations, if required, it needs to be determined in a case-by-case manner. b. According to the designed layout of the microarray, add the probes into the designated wells of the 384-well plate. For the Spike peptides, we created an array of 32 x 20 features with triplicates for each probe ( Figure 1B) . The minimal volume is 10 μL per well, while the maximal volume is 50 μL. The minimal and maximal volume would be different for alternative types of 384-well plates. Note: the probes should be handled on ice. c. Centrifuge the plate at 800 x g for 1 min at 4℃. Pause Point: the prepared plate could be stored at -80℃ for long term (at least one year). 5. Microarray printing a. Use the printer (Super Marathon) to print the probes on PATH substrate slides according to the manufacture's instruction (https://www.arrayjet.co.uk/instruments/printers-in-detail). b. Incubate the microarrays at 4℃ with humidity less than 50% overnight (over 12 h). c. Put the slides into the dedicated slide box (5 or 1 slide(s) per box), seal the box with Parafilm and store the box at -80℃. Pause Point: the microarrays can be stored at -80℃ for several years. 6. GAL (GenePix Array List) file generation a. Create an Excel file and list the position of the probes on 384-well plate ( Table 1) . b. Create GAL file using the printer software by combining of this file with the printing program. The GAL file is for extracting the data from the microarray (Step 17). Optional: Microarray quality evaluation. The steps for this experiment are similar with the following part of Serum Profiling (step 9-17) except that 1) the Blocking buffer were replaced by a commercial blocking buffer (QuickBlock blocking buffer) and 2) the microarray was incubated with an Anti-BSA antibody (1:1000) followed by incubation with a Cy3 labeled Anti-Mouse secondary antibody (1:1000). Critical: It is better to take 1.5 minutes or longer to finish this process to avoid trailing effect (a phenomenon that a big tail attached to the otherwise circular feature). e. Keep the slides at room temperature (22-27℃) and shake at 20-30 rpm (revolutions per minute) on a shaker for 3 h. d. Dry the slides by centrifugation. As list in Materials and Equipment table, we used a specialized slide centrifuge SlideWasher™ 8 with default parameters. Alternatively, a normal staining rack with a compatible centrifuge can be used (500 x g for 1-2 min). The dried slides can be stored at 4℃ in a dry environment for at least one month without compromised signal intensity for serum profiling. 9. Sera sample preparation. Centrifuge the tube at 12,000 x g at 4℃ for 20 min. c. Add 200 μL diluted sample into one well that contains one subarray. For each slide, the last two blocks or wells are added with Incubation buffer and positive sera sample as blank control and positive control, respectively, i. e., block # 13 for blank control, and block #14 for positive control. A blank control is used to monitor unexpected contaminations. Note: 1) The samples can be added with a multichannel pipette. 2) Immediately record the sera number, the block position and the slide number. d. Incubate the slide at room temperature (22-27℃) for 2 h. It is recommended to keep them in a box to prevent evaporation. 11. Pre-wash and wash. a. Remove the solution from each well with a multichannel pipette, add 200 μL Wash buffer (1 x PBS with 0.1% tween-20) to each well, and pipette up and down for 15-20 times. b. Repeat the above step twice. c. Uninstall the 16-well chamber. Critical: Do not leave the surface of the slide dry. d. Place the slide in a container with ~20 mL Wash buffer. Keep the container at room temperature (22-27℃) and shake for 90-100 rpm on a shaker for 10 min. e. Repeat the above step twice. a. Dilute the secondary antibody (1: 1000) with Incubation buffer. b. Remove the Wash buffer from the container and add 15 mL diluted secondary antibodies. Depending on the isotypes of antibody of interest, multiple isotypes could be detected simultaneously, only if the isotype specific fluorescence labeled secondary antibodies are available and compatible with the microarray scanner. Here, we detected the IgG and IgM isotypes using Cy3 (532 nm excitation) and Alexa Fluor 647 (635 nm excitation) labeled isotype specific antibody, respectively ( Figure 1C ). Note: if only one slide is used, a smaller container would be suggested with a smaller incubation volume. c. Keep the container at room temperature (22-27℃) and shake for 20-30 rpm on a shaker for 1 h. Keep the slides shielded from light. 13. Wash and dry a. Remove the solution from the container and add 30 mL Wash buffer. Keep the container at room temperature (22-27℃) and shake for 90-100 rpm on a shaker for 10 min. b. Repeat the above step twice. c. Wash twice with ddH2O. d. Dry the slide by centrifugation. Note: From step 15 on, the slides should be kept under dark. 14. Microarray scanning (link to Troubleshooting, problem 2-4) a. Open the Scanner Genepix 4200A and warm up for 15 min. b. Scan the slide with 532 nm (IgG) and 635 nm (IgM) laser, successively. Recommended scanning parameters: 95% power, PMT 200-300. However, the optimal scanning intensity should be adjusted according to the overall signal intensity. c. After scanning is completed, save the two images as one single TIF file. Critical: 1) All sides which belongs to the same experiment should apply the same parameters. 2) Avoid scraping or touching the surface of the microarrays in any way during the whole process. Quantitative signal intensity (indicating the antibody response) to each peptide for one sample could be acquired by transforming the color intensity of the microarray image to digital signals, which is microarray data extraction. Then, normalization is performed to enable comparations of the signal intensities among the samples. a. Use the software Genepix 6.0 to align the features. The software enables automatic alignment of the GAL file (step 8) to the features of the microarray image, however, it is necessary to double check all aligned features and adjust the incorrectly aligned features. It is a time-consuming process. b. Extract the data of all features to generate a GPR file (raw data) for each slide. a. Calculate signal intensity for each feature: subtract Foreground Median (F median) with Background Median (B median). IgG (532 nm channel) and IgM (635 nm channel) are calculated separately. b. Calculate the average of the triplicates of each probe in the same subarray. 17. Perform normalization among slides. Use the data of the positive sera sample (the data of block #14 for each slide) to calculate the normalization factor (Factor N) for each slide. IgG (532 nm channel) and IgM (635 nm channel) are calculated separately. a. With the data of block #14 of all slides, calculate the overall average signal intensity for each peptide (not for the control probes) to create a reference data set. b. Exclude the peptides with low signal intensity (< 300). c. Calculate the Pearson correlation coefficient between the reference data set and the data of block #14 for each slide. If the coefficient is less than 0.98, exclude some spots far away from the diagonal line. d. Perform linear fitting and calculate the Factor N for each slide using the following function in Excel software: LINEST (Y, X , FALSE), where the Y is the reference data set, the X is the data of block #14 from the slide , and the FLASE means to make the constant term equal to 0. e. For each slide, all the signal intensities are normalized by multiplying the corresponding Factor N. Due to potential non-specific binding, sera from healthy people were set as negative controls to define the basal signal intensity for each peptide, so response frequency could be specifically depicted for each peptide in the patient group. For the map, higher frequency value implies widespread epitopes of the target antigen. The peptides with higher response frequency imply epitopes that elicit stronger antibody responses. 18. Set threshold for each peptide a. Use serum samples from healthy people to define the baseline and calculate the threshold value as Mean value + 2 x SD (standard deviation). b. Calculate positive rate or frequency response for each peptide for the patient group. 19. Generate linear epitope map. For each peptide, provide the total patient sera number and positive number as well as the peptide sequence to the webpage-based tool ImmunomeBrowser (Dhanda et al., 2018) . An epitope map would be generated ( Figure 1E ). For the map, higher frequency value implies widespread epitopes of the target antigen. For SARS-CoV-2 Spike, several hot areas in rich of linear epitopes were indicated in the map. The BSA-peptide conjugation generates an obvious shift in gel electrophoresis when compared with BSA and Sulfo-SMCC activated BSA ( Figure 1A ). The quality of the peptide microarray can be evaluated by Anti-BSA antibody ( Figure 1B) . A qualified peptide microarray exhibits bright, uniform, organized and regular features. For serum derived from infected patients, the microarray image is expected to exhibit high differences among the peptides and normally most of the peptides are negative. The background is low and uniform (Figure 1C) . High correlations could be achieved between the repeats of the positive sera samples analyzed on different slides. The normalization factor (Factor N) could be easily calculated J o u r n a l P r e -p r o o f by linear fitting (Figure 1D ). An example (SARS-CoV-2 Spike) of linear epitope map generated by ImmunomeBrowser is shown ( Figure 1E ). Theoretically, this protocol cannot profile conformational epitopes of the target antigens. It may also fail to detect some linear epitopes that were not covered by set of 12-aa peptides on the microarray. Feature alignment (step 17) could dramatically affect the signals of the features, so apply of the same standard for all features during this process is important to achieve comparable and highly reliable quantitative data. Problem 1: The peptides with highly hydrophobic properties might fail to be synthesized (step 2 in "Before You Begin" section). One possible solution is that extend or shorten the peptides to reduce the number of hydrophobic amino acids. However, we did not try in our experiment. Problem 2: False positive signals of the peptides were observed for the bank control sample (# 13 block) on some slide(s) (step 16). The False positive on the control slide may be caused by one of the following reasons: 1) Contaminations from other serum samples are brought when loading serum samples; 2) Contaminations from other serum samples are brought when performing the pre-wash step. 3) The Blocking buffer, Incubation buffer or Wash buffer is contaminated. Discard the current buffer and prepare new buffer, then re-analyze the samples. Carefully operate during the incubation and wash steps. Spike protein was incubated with Anti-BSA followed by Cy3 labeled secondary antibody. Each spot indicates one probe and each probe is printed in triplicate. The peptides that failed to be synthesized or conjugated were replaced by Printing buffer. C) A representative image of serum profiling with the peptide microarray as same as in B. Green, the 532 nm channel, IgG response; Red, the 635 nm channel, IgM response. D) Calculation of Factor N for slide . Each spot indicates one peptide. The Factor N for the slide is 1.4094, which would be used for normalization. E) An Example of the linear epitope map for SARS-CoV-2 Spike generated by ImmunomeBrowser. The map is reprinted with permission from Li et al., 2020. J o u r n a l P r e -p r o o f J o u r n a l P r e -p r o o f ImmunomeBrowser: A tool to aggregate and visualize complex and heterogeneous epitopes in reference proteins Linear epitopes of SARS-CoV-2 spike protein elicit neutralizing antibodies in COVID-19 patients Linear epitope landscape of the SARS-CoV-2 Spike protein constructed from 1,051 COVID-19 patients Systematic evaluation of SARS-CoV-2 spike protein derived peptides for diagnosis of COVID-19 patients Very high background of the microarray (step 16). Potential Solution: Increase the blocking time or concentration of the BSA in Blocking buffer (step 9); Decrease the concentration of the serum samples (step 12); decrease the power of the scanner (step 16); Change the substrate slides. Many merged features in the microarray (step 16). Use a bigger spot-spot distance when printing (step 7); lower the concentration of the printed samples (step 12); extend the equilibrium time (step 9). Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Sheng-ce Tao (taosc@sjtu.edu.cn). This study did not generate new unique reagents. This study did not generate/analyze datasets/code.