Tviblindi algorithm identifies branching developmental trajectories of human B-cell development and describes abnormalities in RAG-1 and WAS patients.

Detailed knowledge of human B-cell development is crucial for the proper interpretation of inborn errors of immunity and malignant diseases. It is of interest to understand the kinetics of protein expression changes during development, but also to properly interpret the major and possibly alternative developmental trajectories. We have investigated human samples from healthy individuals with the aim of describing all B-cell developmental trajectories. We validated a 30-parameter mass cytometry panel and demonstrated the utility of "vaevictis" visualization of B-cell developmental stages. We used the trajectory inference tool "tviblindi" to exhaustively describe all trajectories leading to all developmental ends discovered in the data. Focusing on Natural Effector B cells, we demonstrated the dynamics of expression of nuclear factors (PAX-5, TdT, Ki-67, Bcl-2), cytokine and chemokine receptors (CD127, CXCR4, CXCR5) in relation to the canonical B-cell developmental stage markers. We observed branching of the memory development, where follicular memory formation was marked by CD73 expression. Lastly, we performed an analysis of two example cases of abnormal B-cell development caused by mutations in RAG-1 and Wiskott-Aldrich syndrome gene in patients with primary immunodeficiency. In conclusion, we developed, validated, and presented a comprehensive set of tools for the investigation of B-cell development in the bone marrow compartment.

Defects in B-lymphopoiesis and B-cell maturation underlie prolonged B-cell depletion in ANCA-associated vasculitis.

OBJECTIVES: B-cell depletion time after rituximab (RTX) treatment is prolonged in antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) compared with other autoimmune diseases. We investigated central and peripheral B-cell development to identify the causes for the defect in B-cell reconstitution after RTX therapy. METHODS: We recruited 91 patients with AAV and performed deep phenotyping of the peripheral and bone marrow B-cell compartment by spectral flow and mass cytometry. B-cell development was studied by in vitro modelling and the role of BAFF receptor by quantitative PCR, western blot analysis and in vitro assays. RESULTS: Treatment-naïve patients with AAV showed low transitional B-cell numbers, suggesting impaired B-lymphopoiesis. We analysed bone marrow of treatment-naïve and RTX-treated patients with AAV and found reduced B-lymphoid precursors. In vitro modelling of B-lymphopoiesis from AAV haematopoietic stem cells showed intact, but slower and reduced immature B-cell development. In a subgroup of patients, after RTX treatment, the presence of transitional B cells did not translate in replenishment of naïve B cells, suggesting an impairment in peripheral B-cell maturation. We found low BAFF-receptor expression on B cells of RTX-treated patients with AAV, resulting in reduced survival in response to BAFF in vitro. CONCLUSIONS: Prolonged depletion of B cells in patients with AAV after RTX therapy indicates a B-cell defect that is unmasked by RTX treatment. Our data indicate that impaired bone marrow B-lymphopoiesis results in a delayed recovery of peripheral B cells that may be further aggravated by a survival defect of B cells. Our findings contribute to the understanding of AAV pathogenesis and may have clinical implications regarding RTX retreatment schedules and immunomonitoring after RTX therapy.

A high-throughput pipeline for validation of antibodies.

Western blotting (WB) is widely used to test antibody specificity, but the assay has low throughput and precision. Here we used preparative gel electrophoresis to develop a capture format for WB. Fractions with soluble, size-separated proteins facilitated parallel readout with antibody arrays, shotgun mass spectrometry (MS) and immunoprecipitation followed by MS (IP-MS). This pipeline provided the means for large-scale implementation of antibody validation concepts proposed by an international working group on antibody validation (IWGAV).

DUX4r, ZNF384r and PAX5-P80R mutated B-cell precursor acute lymphoblastic leukemia frequently undergo monocytic switch.

Recently, we described B-cell precursor acute lymphoblastic leukemia (BCP-ALL) subtype with early switch to the monocytic lineage and loss of the B-cell immunophenotype, including CD19 expression. Thus far, the genetic background has remained unknown. Among 726 children consecutively diagnosed with BCP-ALL, 8% patients experienced switch detectable by flow cytometry (FC). Using exome and RNA sequencing, switch was found to positively correlate with three different genetic subtypes: PAX5-P80R mutation (5 cases with switch out of 5), rearranged DUX4 (DUX4r; 30 cases of 41) and rearranged ZNF384 (ZNF384r; 4 cases of 10). Expression profiles or phenotypic patterns correlated with genotypes, but within each genotype they could not identify cases who subsequently switched. If switching was not taken into account, the B-cell-oriented FC assessment underestimated the minimal residual disease level. For patients with PAX5-P80R, a discordance between FC-determined and PCR-determined MRD was found on day 15, resulting from a rapid loss of the B-cell phenotype. Discordance on day 33 was observed in all the DUX4r, PAX5-P80R and ZNF384r subtypes. Importantly, despite the substantial phenotypic changes, possibly even challenging the appropriateness of BCP-ALL therapy, the monocytic switch was not associated with a higher incidence of relapse and poorer prognosis in patients undergoing standard ALL treatment.

MetaMass, a tool for meta-analysis of subcellular proteomics data.

We report a tool for the analysis of subcellular proteomics data, called MetaMass, based on the use of standardized lists of subcellular markers. We analyzed data from 11 studies using MetaMass, mapping the subcellular location of 5,970 proteins. Our analysis revealed large variations in the performance of subcellular fractionation protocols as well as systematic biases in protein annotation databases. The Excel and R versions of MetaMass should enhance transparency and reproducibility in subcellular proteomics.

ERG deletions in childhood acute lymphoblastic leukemia with DUX4 rearrangements are mostly polyclonal, prognostically relevant and their detection rate strongly depends on screening method sensitivity.

ERG-deletions occur recurrently in acute lymphoblastic leukemia, especially in the DUX4-rearranged subtype. The ERG-deletion was shown to positively impact prognosis of patients with IKZF1-deletion and its presence precludes assignment into IKZF1 plus group, a novel high-risk category on AIEOP-BFM ALL trials. We analyzed the impact of different methods on ERG-deletion detection rate, evaluated ERG-deletion as a potential marker for DUX4-rearranged leukemia, studied its associations with molecular and clinical characteristics within this leukemia subtype, and analyzed its clonality. Using single-nucleotide-polymorphism array, genomic polymerase chain reaction (PCR) and amplicon-sequencing we found ERG-deletion in 34% (16 of 47), 66% (33 of 50) and 78% (39 of 50) of DUX4-rearranged leukemia, respectively. False negativity of ERG-deletion by single-nucleotide-polymorphism array caused IKZF1 plus misclassification in 5 patients. No ERG-deletion was found outside the DUX4-rearranged cases. Within DUX4-rearranged leukemia, the ERG-deletion was associated with higher total number of copy-number aberrations, and, importantly, the ERG-deletion positivity by PCR was associated with better outcome [5-year event-free survival (EFS), ERG-deletion-positive 93% vs. ERG-deletion-negative 68%, P=0.022; 5-year overall survival (OS), ERG-deletion-positive 97% vs. ERG-deletion-negative 75%, P=0.029]. Ultra-deep amplicon-sequencing revealed distinct co-existing ERG-deletions in 22 of 24 patients. In conclusion, our data demonstrate inadequate sensitivity of single-nucleotide-polymorphism array for ERG-deletion detection, unacceptable for proper IKZF1 plus classification. Even using more sensitive methods (PCR/amplicon-sequencing) for its detection, ERG-deletion is absent in 22-34% of DUX4-rearranged leukemia and does not represent an adequately sensitive marker of this leukemia subtype. Importantly, the ERG-deletion potentially stratifies the DUX4-rearranged leukemia into biologically/clinically distinct subsets. Frequent polyclonal pattern of ERG-deletions shows that late origin of this lesion is more common than has been previously described.

Genomic landscape of pediatric B-other acute lymphoblastic leukemia in a consecutive European cohort.

Novel biological subtypes and clinically important genetic aberrations (druggable lesions, prognostic factors) have been described in B-other acute lymphoblastic leukemia (ALL) during the last decade; however, due to a lack of studies on unselected cohorts, their population frequency and mutual associations still have to be established. We studied 110 consecutively diagnosed and uniformly treated childhood B-other patients using single nucleotide polymorphism arrays and whole exome/transcriptome sequencing. The frequency of DUX4-rearranged, BCR-ABL1-like, ZNF384-rearranged, ETV6-RUNX1-like, iAMP21 and MEF2D-rearranged subtypes was 27%, 15%, 5%, 5%, 4%, and 2%, respectively; 43% of cases were not classified into any of these subtypes (B-rest). We found worse early response to treatment in DUX4-rearranged leukemia and a strong association of ZNF384-rearranged leukemia with B-myeloid immunophenotype. Of the druggable lesions, JAK/STAT-class and RAS/RAF/MAPK-class aberrations were found in 21% and 43% of patients, respectively; an ABL-class aberration was found in one patient. A recently described negative prognostic factor, IKZF1plus , was found in 14% of patients and was enriched in (but not exclusive for) BCR-ABL1-like subtype. PAX5 fusions (including 4 novel), intragenic amplifications and P80R mutations were mutually exclusive and only occurred in the B-rest subset, altogether accounting for 20% of the B-other group. PAX5 P80R was associated with a specific gene expression signature, potentially defining a novel leukemia subtype. Our study shows unbiased European population-based frequencies of novel ALL subtypes, recurrent (cyto)genetic aberrations and their mutual associations. This study also strengthens and widens the current knowledge of B-other ALL and provides an objective basis for optimization of current genetic diagnostics.

High-resolution Antibody Array Analysis of Childhood Acute Leukemia Cells.

Acute leukemia is a disease pathologically manifested at both genomic and proteomic levels. Molecular genetic technologies are currently widely used in clinical research. In contrast, sensitive and high-throughput proteomic techniques for performing protein analyses in patient samples are still lacking. Here, we used a technology based on size exclusion chromatography followed by immunoprecipitation of target proteins with an antibody bead array (Size Exclusion Chromatography-Microsphere-based Affinity Proteomics, SEC-MAP) to detect hundreds of proteins from a single sample. In addition, we developed semi-automatic bioinformatics tools to adapt this technology for high-content proteomic screening of pediatric acute leukemia patients.To confirm the utility of SEC-MAP in leukemia immunophenotyping, we tested 31 leukemia diagnostic markers in parallel by SEC-MAP and flow cytometry. We identified 28 antibodies suitable for both techniques. Eighteen of them provided excellent quantitative correlation between SEC-MAP and flow cytometry (p< 0.05). Next, SEC-MAP was applied to examine 57 diagnostic samples from patients with acute leukemia. In this assay, we used 632 different antibodies and detected 501 targets. Of those, 47 targets were differentially expressed between at least two of the three acute leukemia subgroups. The CD markers correlated with immunophenotypic categories as expected. From non-CD markers, we found DBN1, PAX5, or PTK2 overexpressed in B-cell precursor acute lymphoblastic leukemias, LAT, SH2D1A, or STAT5A overexpressed in T-cell acute lymphoblastic leukemias, and HCK, GLUD1, or SYK overexpressed in acute myeloid leukemias. In addition, OPAL1 overexpression corresponded to ETV6-RUNX1 chromosomal translocation.In summary, we demonstrated that SEC-MAP technology is a powerful tool for detecting hundreds of proteins in clinical samples obtained from pediatric acute leukemia patients. It provides information about protein size and reveals differences in protein expression between particular leukemia subgroups. Forty-seven of SEC-MAP identified targets were validated by other conventional method in this study.

ETV6/RUNX1-like acute lymphoblastic leukemia: A novel B-cell precursor leukemia subtype associated with the CD27/CD44 immunophenotype.

We have shown previously that ETV6/RUNX1-positive acute lymphoblastic leukemia (ALL) is distinguishable from other ALL subtypes by CD27pos /CD44low-neg immunophenotype. During diagnostic immunophenotyping of 573 childhood B-cell precursor ALL (BCP-ALL), we identified eight cases with this immunophenotype among "B-other ALL" (BCP-ALL cases negative for routinely tested chromosomal/genetic aberrations). We aimed to elucidate whether these cases belong to the recently described ETV6/RUNX1-like ALL defined by the ETV6/RUNX1-specific gene expression profile (GEP), harboring concurrent ETV6 and IKZF1 lesions. We performed comprehensive genomic analysis using single nucleotide polymorphism arrays, whole exome and transcriptome sequencing and GEP on microarrays. In unsupervised hierarchical clustering based on GEP, five out of seven analyzed CD27pos /CD44low-neg B-other cases clustered with ETV6/RUNX1-positive ALL and were thus classified as ETV6/RUNX1-like ALL. The two cases clustering outside ETV6/RUNX1-positive ALL harbored a P2RY8/CRLF2 fusion with activating JAK2 mutations and a TCF3/ZNF384 fusion, respectively, assigning them to other ALL subtypes. All five ETV6/RUNX1-like cases harbored ETV6 deletions; uniform intragenic ARPP21 deletions and various IKZF1 lesions were each found in three ETV6/RUNX1-like cases. The frequency of ETV6 and ARPP21 deletions was significantly higher in ETV6/RUNX1-like ALL compared with a reference cohort of 42 B-other ALL. In conclusion, we show that ETV6/RUNX1-like ALL is associated with CD27pos /CD44low-neg immunophenotype and identify ARPP21 deletions to contribute to its specific genomic profile enriched for ETV6 and IKZF1 lesions. In conjunction with previously published data, our study identifies the ETV6 lesion as the only common genetic aberration and thus the most likely key driver of ETV6/RUNX1-like ALL.

Lymphocyte enrichment using CD81-targeted immunoaffinity matrix.

In mass cytometry, the isolation of pure lymphocytes is very important to obtain reproducible results and to shorten the time spent on data acquisition. To prepare highly purified cell suspensions of peripheral blood lymphocytes for further analysis on mass cytometer, we used the new CD81+ immune affinity chromatography cell isolation approach. Using 21 metal conjugated antibodies in a single tube we were able to identify all basic cell subsets and compare their relative abundance in final products obtained by density gradient (Ficoll-Paque) and immune affinity chromatography (CD81+ T-catch™) isolation approach. We show that T-catch isolation approach results in purer final product than Ficoll-Paque (P values 0.0156), with fewer platelets bound to target cells. As a result acquisition time of 105 nucleated cells was 3.5 shorter. We then applied unsupervised high dimensional analysis viSNE algorithm to compare the two isolation protocols, which allowed us to evaluate the contribution of unsupervised analysis over supervised manual gating. ViSNE algorithm effectively characterized almost all supervised cell subsets. Moreover, viSNE uncovered previously overseen cell subsets and showed inaccuracies in Maxpar™ Human peripheral blood phenotyping panel kit recommended gating strategy. These findings emphasize the use of unsupervised analysis tools in parallel with conventional gating strategy to mine the complete information from a set of samples. They also stress the importance of the impurity removal to sensitively detect rare cell populations in unsupervised analysis. © 2016 International Society for Advancement of Cytometry.

Characterization of leukemias with ETV6-ABL1 fusion.

To characterize the incidence, clinical features and genetics of ETV6-ABL1 leukemias, representing targetable kinase-activating lesions, we analyzed 44 new and published cases of ETV6-ABL1-positive hematologic malignancies [22 cases of acute lymphoblastic leukemia (13 children, 9 adults) and 22 myeloid malignancies (18 myeloproliferative neoplasms, 4 acute myeloid leukemias)]. The presence of the ETV6-ABL1 fusion was ascertained by cytogenetics, fluorescence in-situ hybridization, reverse transcriptase-polymerase chain reaction and RNA sequencing. Genomic and gene expression profiling was performed by single nucleotide polymorphism and expression arrays. Systematic screening of more than 4,500 cases revealed that in acute lymphoblastic leukemia ETV6-ABL1 is rare in childhood (0.17% cases) and slightly more common in adults (0.38%). There is no systematic screening of myeloproliferative neoplasms; however, the number of ETV6-ABL1-positive cases and the relative incidence of acute lymphoblastic leukemia and myeloproliferative neoplasms suggest that in adulthood ETV6-ABL1 is more common in BCR-ABL1-negative chronic myeloid leukemia-like myeloproliferations than in acute lymphoblastic leukemia. The genomic profile of ETV6-ABL1 acute lymphoblastic leukemia resembled that of BCR-ABL1 and BCR-ABL1-like cases with 80% of patients having concurrent CDKN2A/B and IKZF1 deletions. In the gene expression profiling all the ETV6-ABL1-positive samples clustered in close vicinity to BCR-ABL1 cases. All but one of the cases of ETV6-ABL1 acute lymphoblastic leukemia were classified as BCR-ABL1-like by a standardized assay. Over 60% of patients died, irrespectively of the disease or age subgroup examined. In conclusion, ETV6-ABL1 fusion occurs in both lymphoid and myeloid leukemias; the genomic profile and clinical behavior resemble BCR-ABL1-positive malignancies, including the unfavorable prognosis, particularly of acute leukemias. The poor outcome suggests that treatment with tyrosine kinase inhibitors should be considered for patients with this fusion.

Non-apoptotic FAS signaling controls mTOR activation and extrafollicular maturation in human B cells.

Defective FAS (CD95/Apo-1/TNFRSF6) signaling causes autoimmune lymphoproliferative syndrome (ALPS). Hypergammaglobulinemia is a common feature in ALPS with FAS mutations (ALPS-FAS), but paradoxically, fewer conventional memory cells differentiate from FAS-expressing germinal center (GC) B cells. Resistance to FAS-induced apoptosis does not explain this phenotype. We tested the hypothesis that defective non-apoptotic FAS signaling may contribute to impaired B cell differentiation in ALPS. We analyzed secondary lymphoid organs of patients with ALPS-FAS and found low numbers of memory B cells, fewer GC B cells, and an expanded extrafollicular (EF) B cell response. Enhanced mTOR activity has been shown to favor EF versus GC fate decision, and we found enhanced PI3K/mTOR and BCR signaling in ALPS-FAS splenic B cells. Modeling initial T-dependent B cell activation with CD40L in vitro, we showed that FAS competent cells with transient FAS ligation showed specifically decreased mTOR axis activation without apoptosis. Mechanistically, transient FAS engagement with involvement of caspase-8 induced nuclear exclusion of PTEN, leading to mTOR inhibition. In addition, FASL-dependent PTEN nuclear exclusion and mTOR modulation were defective in patients with ALPS-FAS. In the early phase of activation, FAS stimulation promoted expression of genes related to GC initiation at the expense of processes related to the EF response. Hence, our data suggest that non-apoptotic FAS signaling acts as molecular switch between EF versus GC fate decisions via regulation of the mTOR axis and transcription. The defect of this modulatory circuit may explain the observed hypergammaglobulinemia and low memory B cell numbers in ALPS.

Single-cell protein profiling defines cell populations associated with triple-negative breast cancer aggressiveness.

Triple-negative breast cancer (TNBC) is an aggressive and complex subtype of breast cancer that lacks targeted therapy. TNBC manifests characteristic, extensive intratumoral heterogeneity that promotes disease progression and influences drug response. Single-cell techniques in combination with next-generation computation provide an unprecedented opportunity to identify molecular events with therapeutic potential. Here, we describe the generation of a comprehensive mass cytometry panel for multiparametric detection of 23 phenotypic markers and 13 signaling molecules. This single-cell proteomic approach allowed us to explore the landscape of TNBC heterogeneity, with particular emphasis on the tumor microenvironment. We prospectively profiled freshly resected tumors from 26 TNBC patients. These tumors contained phenotypically distinct subpopulations of cancer and stromal cells that were associated with the patient's clinical status at the time of surgery. We further classified the epithelial-mesenchymal plasticity of tumor cells, and molecularly defined phenotypically diverse populations of tumor-associated stroma. Furthermore, in a retrospective tissue-microarray TNBC cohort, we showed that the level of CD97 at the time of surgery has prognostic potential.

An automated analysis of highly complex flow cytometry-based proteomic data.

The combination of color-coded microspheres as carriers and flow cytometry as a detection platform provides new opportunities for multiplexed measurement of biomolecules. Here, we developed a software tool capable of automated gating of color-coded microspheres, automatic extraction of statistics from all subsets and validation, normalization, and cross-sample analysis. The approach presented in this article enabled us to harness the power of high-content cellular proteomics. In size exclusion chromatography-resolved microsphere-based affinity proteomics (Size-MAP), antibody-coupled microspheres are used to measure biotinylated proteins that have been separated by size exclusion chromatography. The captured proteins are labeled with streptavidin phycoerythrin and detected by multicolor flow cytometry. When the results from multiple size exclusion chromatography fractions are combined, binding is detected as discrete reactivity peaks (entities). The information obtained might be approximated to a multiplexed western blot. We used a microsphere set with >1,000 subsets, presenting an approach to extract biologically relevant information. The R-project environment was used to sequentially recognize subsets in two-dimensional space and gate them. The aim was to extract the median streptavidin phycoerythrin fluorescence intensity for all 1,000+ microsphere subsets from a series of 96 measured samples. The resulting text files were subjected to algorithms that identified entities across the 24 fractions. Thus, the original 24 data points for each antibody were compressed to 1-4 integrated values representing the areas of individual antibody reactivity peaks. Finally, we provide experimental data on cellular protein changes induced by treatment of leukemia cells with imatinib mesylate. The approach presented here exemplifies how large-scale flow cytometry data analysis can be efficiently processed to employ flow cytometry as a high-content proteomics method.

The clinical value of circulating free tumor DNA in testicular germ cell tumor patients.

BACKGROUND: Testicular germ cell tumors (TGCT) are unique malignancies of young adult men; their biology is, however, underexplored and there has not been much progress in their treatment for decades. Circulating free tumor DNA (cfDNA) analysis represents a promising way of discovering novel diagnostic and treatment options. OBJECTIVE: The study evaluates the clinical value of cfDNA detection in TGCT patients. DESIGN AND METHODS: Total cfDNA concentration and ratio of its 2 main fragments (180 and 360 bp) were evaluated by spectrophotometry, capillary electrophoresis and qPCR in peripheral blood plasma of 96 TGCT patients (173 samples) and 31 normal controls. Non-parametric tests were used for statistical analyses. RESULTS: The total cfDNA concentration was significantly higher in TGCT than in controls (P < 0.0001), with the highest levels at disease progression, but with no clear threshold between malignant and normal samples. Patients with positive tumor markers had higher cfDNA concentrations than those with negative markers (P = 0.01). Longer 360 bp cfDNA fragments were found in 58% of TGCT patients including almost all samples from relapse or disease progression but no normal controls (P < 0.0001). CONCLUSION: Total cfDNA levels are significantly increased in TGCT patients but without a clear threshold separating normal and tumor samples, thus total cfDNA amount itself is not a sensitive enough marker to identify or monitor TGCT. Longer cfDNA fragments have been found exclusively in a proportion of tumors and predominantly at disease progression, representing a novel potential marker for TGCT monitoring that would deserve further exploration.

Minimal residual disease in BCR::ABL1-positive acute lymphoblastic leukemia: different significance in typical ALL and in CML-like disease.

Recently, we defined "CML-like" subtype of BCR::ABL1-positive acute lymphoblastic leukemia (ALL), resembling lymphoid blast crisis of chronic myeloid leukemia (CML). Here we retrospectively analyzed prognostic relevance of minimal residual disease (MRD) and other features in 147 children with BCR::ABL1-positive ALL (diagnosed I/2000-IV/2021, treated according to EsPhALL (n = 133) or other (n = 14) protocols), using DNA-based monitoring of BCR::ABL1 genomic breakpoint and clonal immunoglobulin/T-cell receptor gene rearrangements. Although overall prognosis of CML-like (n = 48) and typical ALL (n = 99) was similar (5-year-EFS 60% and 49%, respectively; 5-year-OS 75% and 73%, respectively), typical ALL presented more relapses while CML-like patients more often died in the first remission. Prognostic role of MRD was significant in the typical ALL (p = 0.0005 in multivariate analysis for EFS). In contrast, in CML-like patients MRD was not significant (p values > 0.2) and inapplicable for therapy adjustment. Moreover, in the typical ALL, risk-prediction could be further improved by considering initial hyperleukocytosis. Early distinguishing typical BCR::ABL1-positive ALL and CML-like patients is essential to enable optimal treatment approach in upcoming protocols. For the typical ALL, tyrosine-kinase inhibitors and concurrent chemotherapy with risk-directed intensity should be recommended; in the CML-like disease, no relevant prognostic feature applicable for therapy tailoring was found so far.

Multiplexed immuno-precipitation with 1725 commercially available antibodies to cellular proteins.

Antibody array analysis of complex samples requires capture reagents with exceptional specificity. The frequency of antibodies with label-based detection may be as low as 5%. Here, however, we show that as many as 25% of commercially available antibodies are useful when biotinylated cellular proteins are fractionated by size exclusion chromatography (SEC) first. A microsphere multiplex with 1725 antibodies to cellular proteins was added to 24 SEC fractions, labelled with streptavidin and analyzed by flow cytometry (microsphere-based affinity proteomics, MAP) The SEC-MAP approach resolved different targets captured by each antibody as reactivity peaks across the separation range of the SEC column (10-670kDa). Complex reactivity profiles demonstrated that most antibodies bound more than one target. However, specific binding was readily detected as reactivity peaks common for different antibodies to the same protein. We optimized sample preparation and found that amine-reactive biotin rarely inhibited antibody binding when the biotin to lysine ratio was kept below 1:1 during labelling. Moreover, several epitopes that were inaccessible to antibodies in native proteins were unmasked after heat denaturation with 0.1% of SDS. The SEC-MAP format should allow researchers to build multiplexed assays with antibodies purchased for use in e.g. Western blotting.