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.

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.

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.

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.

A novel class of ZNF384 aberrations in acute leukemia.

Fusion of the ZNF384 gene as the 3' partner to several different 5' partner genes occurs recurrently in B-cell precursor acute lymphoblastic and mixed phenotype B/myeloid leukemia. These canonical fusions (ZNF384r) contain the complete ZNF384 coding sequence and are associated with a specific gene expression signature. Cases with this signature, but without canonical ZNF384 fusions (ZNF384r-like cases), have been described previously. Although some have been shown to harbor ZNF362 fusions, the primary aberrations remain unknown in a major proportion. We studied 3 patients with the ZNF384r signature and unknown primary genetic background and identified a previously unknown class of genetic aberration affecting the last exon of ZNF384 and resulting in disruption of the C-terminal portion of the ZNF384 protein. Importantly, in 2 cases, the ZNF384 aberration, indel, was missed during the bioinformatic analysis but revealed by the manual, targeted reanalysis. Two cases with the novel aberrations had a mixed (B/myeloid) immunophenotype commonly associated with canonical ZNF384 fusions. In conclusion, we present leukemia cases with a novel class of ZNF384 aberrations that phenocopy leukemia with ZNF384r. Therefore, we show that part of the so-called ZNF384r-like cases represent the same genetic subtype as leukemia with canonical ZNF384 fusions.

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.

Molecular Basis of Cisplatin Resistance in Testicular Germ Cell Tumors.

The emergence of cisplatin (CDDP) resistance is the main cause of treatment failure and death in patients with testicular germ cell tumors (TGCT), but its biologic background is poorly understood. To study the molecular basis of CDDP resistance in TGCT we prepared and sequenced CDDP-exposed TGCT cell lines as well as 31 primary patients' samples. Long-term exposure to CDDP increased the CDDP resistance 10 times in the NCCIT cell line, while no major resistance was achieved in Tera-2. Development of CDDP resistance was accompanied by changes in the cell cycle (increase in G1 and decrease in S-fraction), increased number of acquired mutations, of which 3 were present within ATRX gene, as well as changes in gene expression pattern. Copy number variation analysis showed, apart from obligatory gain of 12p, several other large-scale gains (chr 1, 17, 20, 21) and losses (chr X), with additional more CNVs found in CDDP-resistant cells (e.g., further losses on chr 1, 4, 18, and gain on chr 8). In the patients' samples, those who developed CDDP resistance and died of TGCT (2/31) showed high numbers of acquired aberrations, both SNPs and CNVs, and harbored mutations in genes potentially relevant to TGCT development (e.g., TRERF1, TFAP2C in one patient, MAP2K1 and NSD1 in another one). Among all primary tumor samples, the most commonly mutated gene was NSD1, affected in 9/31 patients. This gene encoding histone methyl transferase was also downregulated and identified among the 50 most differentially expressed genes in CDDP-resistant NCCIT cell line. Interestingly, 2/31 TGCT patients harbored mutations in the ATRX gene encoding a chromatin modifier that has been shown to have a critical function in sexual differentiation. Our research newly highlights its probable involvement also in testicular tumors. Both findings support the emerging role of altered epigenetic gene regulation in TGCT and CDDP resistance development.

The EuroFlow PID Orientation Tube for Flow Cytometric Diagnostic Screening of Primary Immunodeficiencies of the Lymphoid System.

In the rapidly evolving field of primary immunodeficiencies (PID), the EuroFlow consortium decided to develop a PID orientation and screening tube that facilitates fast, standardized, and validated immunophenotypic diagnosis of lymphoid PID, and allows full exchange of data between centers. Our aim was to develop a tool that would be universal for all lymphoid PIDs and offer high sensitivity to identify a lymphoid PID (without a need for specificity to diagnose particular PID) and to guide and prioritize further diagnostic modalities and clinical management. The tube composition has been defined in a stepwise manner through several cycles of design-testing-evaluation-redesign in a multicenter setting. Equally important appeared to be the standardized pre-analytical procedures (sample preparation and instrument setup), analytical procedures (immunostaining and data acquisition), the software analysis (a multidimensional view based on a reference database in Infinicyt software), and data interpretation. This standardized EuroFlow concept has been tested on 250 healthy controls and 99 PID patients with defined genetic defects. In addition, an application of new EuroFlow software tools with multidimensional pattern recognition was designed with inclusion of maturation pathways in multidimensional patterns (APS plots). The major advantage of the EuroFlow approach is that data can be fully exchanged between different laboratories in any country of the world, which is especially of interest for the PID field, with generally low numbers of cases per center.

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.

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).

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.

Common Variable Immunodeficiency patients with a phenotypic profile of immunosenescence present with thrombocytopenia.

Common variable immunodeficiency (CVID) is a heterogeneous group of diseases. Our aim was to define sub-groups of CVID patients with similar phenotypes and clinical characteristics. Using eight-color flow cytometry, we analyzed both B- and T-cell phenotypes in a cohort of 88 CVID patients and 48 healthy donors. A hierarchical clustering of probability binning "bins" yielded a separate cluster of 22 CVID patients with an abnormal phenotype. We showed coordinated proportional changes in naïve CD4+ T-cells (decreased), intermediate CD27- CD28+ CD4+ T-cells (increased) and CD21low B-cells (increased) that were stable for over three years. Moreover, the lymphocytes' immunophenotype in this patient cluster exhibited features of profound immunosenescence and chronic activation. Thrombocytopenia was only found in this cluster (36% of cases, manifested as Immune Thrombocytopenia (ITP) or Evans syndrome). Clinical complications more frequently found in these patients include lung fibrosis (in 59% of cases) and bronchiectasis (55%). The degree of severity of these symptoms corresponded to more deviation from normal levels with respect to CD21low B-cells, naïve CD4+ and CD27− CD28+ CD4+ T-cells. Next-generation sequencing did not reveal any common genetic background. We delineate a subgroup of CVID patients with activated and immunosenescent immunophenotype of lymphocytes and distinct set of clinical complications without common genetic background.

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.

Appearance of cytomegalovirus-specific T-cells predicts fast resolution of viremia post hematopoietic stem cell transplantation.

BACKGROUND: Cytomegalovirus (CMV) specific T-cells are known to provide long-term control of CMV reactivation, which is a frequent complication of hematopoietic stem cell transplantation. We have studied 58 pediatric patients after hematopoietic stem cell transplantation who suffered from CMV reactivation to reveal which functional T cell subset is best correlating with successful reactivation resolution and which protects from reactivation episode. METHODS: Detection of 30 combinatorial subsets of four types of response to ex vivo CMV stimulation (IFNγ secretion, IL-2 secretion, CD40L upregulation and degranulation) that were detectable on either CD8+ or CD4+ T cells through flow cytometry intracellular cytokine staining was used. RESULTS: We found that the presence of CD8+ dual positive (IFNγ+ and IL-2+) cells is the most accurate functional parameter that can predict fast resolution of CMV reactivation. Next, we show that the presence of CD8+ dual positive (IFNγ+ and IL-2+) and CD8+ IFNγ+ cells provides a protective effect (a hazard risk of 0.28 (confidence interval 0.18 - 0.43) and 0.45 (CI 0.27 - 0.75), respectively) and the presence of corticotherapy increases the risk of reactivation (HR 2.47 (CI 1.82-3.36)). Thus, a patient without corticotherapy and with both of the critical T cell subsets present has a cumulative 19.6 times lower risk of developing CMV reactivation than a patient on corticotherapy and without CD8+ dual positive (IFNγ+ and IL-2+) or CD8+ IFNγ+ cells. CONCLUSIONS: We have established parameters of CMV specific functional response ex vivo that can be used in assisting clinical management of patients with CMV reactivation. © 2015 International Clinical Cytometry Society.

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.