Dysferlin and Other Non-Red Cell Proteins Accumulate in the Red Cell Membrane of Diamond-Blackfan Anemia Patients

Diamond Blackfan Anemia (DBA) is a congenital anemia usually caused by diverse mutations in ribosomal proteins. Although the genetics of DBA are well characterized, the mechanisms that lead to macrocytic anemia remain unclear. We systematically analyzed the proteomes of red blood cell membranes from multiple DBA patients to determine whether abnormalities in protein translation or erythropoiesis contribute to the observed macrocytosis or alterations in the mature red blood cell membrane. In depth proteome analysis of red cell membranes enabled highly reproducible identification and quantitative comparisons of 1100 or more proteins. These comparisons revealed clear differences between red cell membrane proteomes in DBA patients and healthy controls that were consistent across DBA patients with different ribosomal gene mutations. Proteins exhibiting changes in abundance included those known to be increased in DBA such as fetal hemoglobin and a number of proteins not normally found in mature red cell membranes, including proteins involved in the major histocompatibility complex class I pathway. Most striking was the presence of dysferlin in the red blood cell membranes of DBA patients but absent in healthy controls. Immunoblot validation using red cell membranes isolated from additional DBA patients and healthy controls confirmed a distinct membrane protein signature specific to patients with DBA.


Introduction
Diamond Blackfan Anemia (DBA) is a rare, congenital anemia characterized by a paucity of erythroid progenitor cells in the bone marrow, reticulocytopenia, and an increased size (macrocytosis) of the remaining circulating red blood cells (RBCs), with elevated adenosine deaminase activity and increased fetal hemoglobin. [1] In the majority of patients, DBA is caused by a heterozygous mutation in a gene that encodes one of several ribosomal proteins in either the small or large ribosomal subunits. To date, mutations in RPS19, RPS17, RPS24, RPL35A, RPL5, RPL11, RPS7, RPS10, RPS26, RPL15 and RPL26 have been associated with DBA, with RPS19 being most frequently mutated, accounting for 20-25% of known mutations. [2] Currently, mutations in these ribosomal genes account for disease in at least 50-70% of patients with DBA. [2,3] The phenotypic presentation of DBA is highly variable, ranging from 1) asymptomatic, to 2) a macrocytic anemia that improves with corticosteroid treatment, to 3) patients whose anemia is so severe that they are dependent on RBC transfusions. Interestingly, the degree of anemia often varies during the course of disease and is usually not present at birth but develops during the first year of life. Anemia might be initially responsive to corticosteroids but may become so severe that the patient becomes transfusion-dependent later in life. In about 15% of patients the anemia improves spontaneously to the extent that no further treatment is required [4].
While the genetic basis of DBA is well characterized, the mechanism by which impaired ribosome biosynthesis leads to macrocytic RBCs with elevated levels of fetal hemoglobin and adenosine deaminase activity are still unclear. DBA mainly affects the erythroid lineage and does not respond to treatment with erythropoietin. [5] In vitro colony forming assays with bone marrow cells from patients with DBA demonstrated a decreased proliferation and increased apoptosis in erythroid colonies that pinpoints the defect to the erythropoietin-dependent stages of RBC differentiation. [6,7,8] Additionally, in vitro culture studies and studies in animal models of DBA demonstrated that inactivation of p53 improves or even restores RBC production, suggesting a p53 dependent cell cycle checkpoint may be involved in pathogenesis [9,10].
Previous studies have investigated gene expression in purified bone marrow populations isolated from DBA patients in remission, [11] a mouse model using fetal liver erythroid cells with knockdown of RPS19, [12] and a transgenic mouse model with an RPS19 mutation that causes a DBA-like phenotype. [13] Results from these studies highlighted the upregulation of apoptotic genes such as TNFRSF10B and FAS [11] as well as decreased expression of transcription factor MYB [11,12] and its target, c-kit. [12,13] Because DBA is a disease of defective ribosome biogenesis, we hypothesized that some of the RBC abnormalities may be due to faulty translation or protein sorting and degradation during erythroid maturation. Hence, quantitative label-free proteomic analysis of RBC membranes from DBA patients could provide new insights into downstream effects of altered ribosome biogenesis in RBC progenitors that would not be apparent in gene expression analysis. Additionally, any differences in the RBC proteome could lead to new diagnostic or prognostic DBA biomarkers or identify novel therapeutic targets.
In this study, we performed label-free proteomic analysis of RBC membranes from four DBA patients with differing ribosomal gene mutations and compared the results with those from matched healthy controls. A number of proteins showed altered abundance levels in DBA RBC membranes compared with controls, however, none of the major RBC membrane proteins exhibited significant abundance changes. The most dramatic difference was the presence, in DBA RBCs, of proteins not normally observed in RBC membranes, such as dysferlin and members of the MHC class I pathway. Furthermore, most observed protein changes did not appear to correlate with specific ribosomal protein gene mutations. Our results suggest that DBA RBCs have a characteristic pattern of RBC membrane protein composition that differs markedly from the proteome of normal RBC membranes and points towards the activation of a cellular stress response that appears to be relatively specific to DBA.

Ethics Statement
Whole blood was collected from individuals after informed written consent was obtained. All experiments were in accordance with the principles expressed in the Declaration of Helsinki with procedures approved by the IRB of the Children's Hospital of Philadelphia and University of Pennsylvania.

RBC Membrane Processing
Approximately 5-10 mL whole blood was collected in K 2 EDTA tubes and all materials were kept at 0-4uC throughout processing to minimize proteolysis. RBC purification and membrane preparation was performed as previously described. [14] Briefly, RBCs were purified by filtering with a leukocyte depletion filter (PlasmodipurH, Accurate Chemical & Scientific Corp) and washing the cells 4 times in PBS with subsequent removal of the top layer of RBC. RBC membranes were prepared by lysing with hypotonic buffer (5 mM sodium phosphate, 1 mM EDTA, 0.1 mM DFP (diisopropylfluorophosphate), pH 8.0) followed by 4-5 washes in the same buffer, until white RBC membranes were obtained. Protein concentration was determined using the Modified Lowry Protein Assay (Thermo Scientific). Aliquots of each sample were stored immediately at 280uC until future use. Pools of control RBC membranes were prepared by combining equal protein amounts of each individual control sample.

Sample Preparation and LC-MS/MS
RBC membrane fractionation was performed by separating 12 mg total protein on 10% Bis-Tris mini gels with MOPS SDS buffer (Invitrogen) until the dye front reached 3 cm. Three parallel gel lanes of the same sample were sliced into 3061 mm fragments and corresponding slices were pooled. Pooled gel slices were reduced, alkylated, digested with trypsin and analyzed by LC-MS/ MS as previously described [14].

Data Processing and Label-Free Analysis
Data was processed through the Rosetta Elucidator software package (Ver. 3.3, Rosetta Biosoftware). Using this software, the raw data was aligned across LC-MS/MS full scans for corresponding gel slices to align identical MS features (discrete m/z signals at a given retention time), which were quantitated based on peak heights of ion intensities, as previously described. [15,16,17] Corresponding peptide sequences were identified using the SEQUEST algorithm (Ver. 28, rev. 13, University of Washington) and the human UniRef 100 database (Ver. June, 2011) that was expanded to include expected environmental contaminants. [17] Proteins were annotated using ProteinTeller scores .0.95 and PeptideTeller scores .0.8 to obtain a protein false discovery rate of ,0.5%. Raw data and search results can be found in the Peptide Atlas data repository with the dataset identifier PASS00334 (www.peptideatlas.org).

Statistical Analysis
Protein lists were generated with a 2 peptide minimum and common protein contaminants (e.g. trypsin, keratins, etc.) were excluded. An intensity cutoff of 1610 4 was applied to peptides to eliminate erroneous protein intensities attributed to background noise. Final protein intensities were calculated by summing the associated peptide intensities followed by normalizing based on the median log10 intensity ratio of DBA/Control. In samples where no protein was detected, to avoid division by zero a fold change threshold of 6100 was applied. 95% confidence intervals were calculated by fitting the standard deviation versus the mean protein intensities of both control datasets to an exponential decay function. [14] This confidence interval was then transposed onto a plot of protein intensities from pair wise DBA versus control to identify significantly changed proteins.

Confocal Microscopy
The procedure for RBC fixing and permeabilization was adapted from a study by Franco et al and performed at room temperature. [18] Washed RBCs were fixed in 1 mL of freshly prepared 0.5% acrolein (VWR) in PBS for 5 min followed by quenching with 4 washes of cold 0.1 M glycine/PBS. Cells were permeabilized with equal volume of 0.1% Triton-X 100/PBS for 5 minutes followed by 3 washes with blocking buffer (0.1 M glycine/ PBS with 0.2% cold water fish gelatin (Sigma-Aldrich)) and incubated in the blocking buffer for a minimum of 1 hour. RBCs were incubated with a mouse anti-dysferlin antibody and rabbit anti-a-spectrin antibody, [19] for a minimum of 40 minutes and subsequently washed 3 times with blocking buffer. RBCs were then incubated for 40 minutes with a goat anti-mouse antibody conjugated with Alexa FluorH 488 and a goat anti-rabbit antibody conjugated with Alexa FluorH 633 (Molecular Probes) followed by 3 washes. Cells were mounted onto poly-l-lysine coated slides with Aqua-Mount (Thermo-Scientific). Images were obtained with a Leica TCS SP5 II scanning laser confocal microscope (Leica Microsystems) equipped with AOBS and HyD detectors.

Evaluation of DBA Patient Proteomes
A total of four transfusion-independent DBA patients were chosen for proteomic analysis (Table 1). We focused on the analysis of RBC membranes, i.e. ''white ghosts'', comprising mainly membrane and membrane skeleton proteins. RBC membrane preparations from two female DBA patients (D1 and D2) were compared to two RBC membrane samples (C1 and C2) each pooled from three age, sex and ethnically matched healthy control donors (C1a-C1c and C2a-C2c). Similarly, two male DBA patients (D3 and D4) were compared to two pools (C3 and C4) of two matched healthy control donors (C3a-C3b and C4a-C4b). Control pools were used to average natural variations in protein levels while at the same time keeping the scope of the proteomics analyses within practical limits. For each DBA patient, two blood samples drawn at least three months apart were prepared separately for replicate analyses, with the exception of D2 in which two aliquots of the same blood draw were processed in parallel because only one blood sample was available. Pair wise comparisons were performed using normalized average protein intensities from DBA patient replicates and appropriate matched control pools. An overview of the experimental scheme is shown in Figure 1. The comparison between D1 and D2 vs. controls C1 and C2 resulted in 1304 high confidence protein identifications, and the comparison between D3 and D4 vs. C3 and C4 identified 1116 proteins. Significant changes in protein abundance were determined by calculating 95% confidence intervals from the control pools, which were imposed onto protein intensity plots for each DBA patient and their appropriate control pools (Figures S1 and S2 in File S1). Significantly changed proteins were identified with log10 intensities outside of the 95% confidence intervals that were defined using the healthy control datasets. Fold changes for those significantly changed proteins were then calculated from the ratio of log10 protein intensities of DBA patients versus matched controls.

Comparative Analysis of Significantly Changed Proteins
To further improve the reliability of protein abundance changes, we focused on proteins identified and quantitated by at least three peptides and that were found to be significantly changed in two or more DBA patients based on the 95% confidence intervals described above. Using these criteria, a total of 30 proteins were increased in the RBC membranes of DBA patients ( Table 2). Six proteins, dysferlin (DYSF), transporter 1 ATP binding cassette sub-family B (TAP1), major vault protein (MVP), major histocompatibility complex class I A (HLAA) and serine palmitoyltransferase 1 and 2 (SPTLC1, SPTLC2) were found to be present in all four DBA patients but were not detected in the control pools. Additional members of the MHC class I pathway such as b-2-microglobulin (B2M), antigen peptide transporter 2 (TAP2), MHC class I antigen (HLAB), and tapasin were also significantly elevated in DBA patients. C-1-tetrahydrofolate synthase (MTHFD1), gamma-adducin (ADD3), and epoxide hydrolase (EPHX1) were also significantly elevated in three of the four DBA patients tested. Fetal hemoglobin (HbG1 and HbG2), a known marker of DBA RBCs, was elevated in three of the four DBA patients tested. Interestingly, prior DBA studies focused on total hemoglobin, while the hemoglobin fraction analyzed here is that tightly associated with the membrane. [20] Prior studies suggest that the membrane-associated fraction of hemoglobin reflects the composition of the soluble fraction, although this remains to be verified [21,22]. Table 3 lists the 32 proteins that were significantly decreased in two or more DBA patients. Consistent with the above fHb results, the adult forms of membrane-bound hemoglobin (HbB, HbD, and HbA) were decreased in all four DBA patients compared to controls. Complement receptor 1 (CR1) was not detected in any of the DBA patients but was present in all of the control pools. Several other proteins were decreased in three of the four DBA patients including phosphoribosyl pyrophosphate synthase-associated protein 2 (PRPSAP2), class IVb beta tubulin (TUBB2C), tubulin alpha-ubiquitous chain (TUBA1B), and intercellular adhesion molecule 4 (ICAM4).
The above comparisons evaluated how individual proteins changed relative to that same protein in the controls. To assess abundance levels across proteins, relative abundance was estimated for each protein by summing the top three peptide intensities to minimize bias due to protein size and variations in peptide signal intensities in the mass spectrometer. [23] Figure 2B shows the relative protein intensities for DBA patients of major RBC proteins as well as a subset of proteins from Tables 2 and 3 that were significantly changed in three or more DBA patients. The proteins differentially expressed in DBA patients were of intermediate intensities (10 4 -10 8 ) compared to the most abundant RBC membrane proteins i.e. spectrin and band 3, which showed intensities of 10 9 -10 10 . Therefore, the majority of proteins that were significantly changed in DBA patients were of medium to high abundance and well above the limit of detection for low abundance proteins (approximately 10 4 signal intensity).
To determine whether proteins with altered abundance levels belonged to specific pathways or processes, these proteins were further characterized using gene ontology functional annotation with DAVID (http://david.abcc.ncifcrf.gov/home) ( Table S1 in File S1). Proteins that were significantly increased in DBA patients involved processes related to antigen presentation via the major histocompatibility class I pathway, cell cycle processing, and immune response. Significantly decreased proteins were categorized into oxygen transport and response to oxidative stress.

Immunoblot Validation
Immunoblot validation was performed for proteins of interest where suitable antibodies were available. RBC ghost preparations from 18 healthy controls and nine transfusion-independent DBA patients (Table 1) were analyzed. Figure 3A is a representative immunoblot analysis of six controls and eight DBA patients, four of whom were used in the LC-MS/MS analysis. Dysferlin was present, although at variable levels, in all white ghost preparations from DBA patients but was not detectable in any of the 18 healthy control RBC membranes tested. Dysferlin was barely detectable in DBA patient D6 with the conditions shown in Figure 3A, however, subsequent analysis using more extensive reaction conditions showed a clear band, indicating dysferlin was present albeit at a much lower level than in RBC membranes from other patients (data not shown). ADD3 was also consistently increased in all DBA samples tested compared to normal controls, although to a lesser  extent, which is consistent with the relatively low fold changes in the proteomic data (ranging from 2.3-6.4 fold). MVP was largely absent in control white ghost samples, but variable protein levels were found in five of the DBA patient samples ( Figure 3A). The levels of SPTLC1 were increased in DBA patients, with the majority of controls showing little to no SPTLC1. A distinct increase in SPTLC2 levels can be seen in all of the DBA patients tested compared to the low levels seen in the controls. Additionally, HLA-ABC was found to be present in variable amounts in the white ghost fractions of six DBA patients as well as two control samples. ICAM4 can be seen with variable abundance in the control samples but is present at low to moderate levels in six DBA patients and cannot be detected in two of the DBA patients. CR1 was largely absent in DBA patients whereas moderate levels were detectable in all of the controls, which is consistent with the proteomic data.

Dysferlin in the RBC Membrane
The most striking result was the identification of dysferlin in the RBC membrane of DBA patients since it has never been described to be present in the RBC membrane. Dysferlin has been well studied because of its involvement in forms of muscular dystrophy Figure 2. Fold change and relative abundance of significantly changed proteins with comparison to major RBC membrane proteins. A) Fold change relative to controls shows minimal variability for normal, major RBC membrane proteins and more variability between patients for the proteins that were significantly changed in three or more DBA patients. Bracketed proteins were not detected and apparent fold change is based on a fold change threshold of 6100. B) Magnitude of protein abundance in RBCs were calculated for each protein based on the sum of the top three peptide intensities. Significantly changed proteins ranged from low to high abundance as compared to major RBC membrane proteins that were present in high abundance. doi:10.1371/journal.pone.0085504.g002 but it has never been reported to be present in the RBC membrane in previous proteome analyses. [24,25,26,27] Consistent with these earlier results, we did not detect dysferlin in any of the control donor RBC membranes. We used confocal microscopy to evaluate the location of dysferlin in the RBCs of DBA patients. Figure 3B shows representative images of dysferlin in RBC membranes of DBA patient D1, as seen with immunostaining of dysferlin (green) relative to spectrin as a control (red). Consistent with the proteome analyses, confocal microscopy of RBC from control C2 showed no detectable dysferlin. Interestingly, the levels of dysferlin in the RBCs from individual DBA patients were variable, with some cells displaying intense fluorescence, while other cells showed minimal or undetectable staining. To further evaluate dysferlin location in RBC membranes of DBA patients, we extracted spectrin and actin from ghosts and analyzed this extract and inside-out vesicles using dysferlin immunoblots (Fig. 3C). Dysferlin was present in the inside-out vesicles but was not detectable in the cytoskeletal extract which is consistent with the fact that dysferlin is a transmembrane protein.
In addition to increased depth of proteome analysis, separation of proteins by 1D gel electrophoresis prior to LC-MS/MS provided interesting information on the size distribution of proteins by comparison of spectral counts across gel slices for individual proteins. Figure 3D shows the relative abundance of dysferlin at different molecular weights for patients D1 and D2. In both cases, the majority of the protein was detected at its expected molecular weight, in the fractions corresponding to 250-200 kDa. However, for patient D2, some dysferlin was detected at lower and higher molecular weights (slices 14-16 and 1-3, respectively), which was indicative of proteolysis and covalent crosslinking of the molecule, possibly to itself, as it may be present as a homodimer. [28] The results from the immunoblot correlated with the relative abundance and size distribution of dysferlin in the LC-MS/MS data (Fig. 3).
To further investigate whether dysferlin on the RBC is a general marker of stress hematopoiesis we investigated other blood disorders for the presence of dysferlin on the RBC membrane with western blot analysis. We screened 19 patients with macrocytosis associated with various diagnoses including paroxysmal nocturnal hemoglobinuria, dyskeratosis congenita, acquired aplastic anemia, congenital neutropenia, Shwachman Diamond syndrome, diserythropoietic anemia (Gata1 mutation), and chemotherapy-related macrocytosis (Table 4). Interestingly, of all the samples tested, dysferlin was not in any other inherited bone marrow failure syndrome tested, however it was found on the RBC of four of the eight acquired aplastic anemia patients. Western blot analysis of RBC membrane preparations show low to moderate levels of dysferlin in these acquired aplastic anemia patients compared to DBA patients ( Figure S3 in File S1).

Discussion
RBCs in DBA are distinguished from normal RBC by macrocytosis, increased fetal hemoglobin and elevated adenosine deaminase activity. Here we used 1D SDS gel separation followed by LC-MS/MS of RBC ghosts to investigate changes in relative protein levels of over one thousand RBC membrane proteins to identify a unique protein signature that further distinguishes DBA patients from normal controls. To confirm that protein changes were not due to proteins from contaminating cell types, we searched our datasets for proteins that were recently listed in a review of 50 blood cell-type specific proteins. Of the 41 proteins specific to monocytes, T cells, platelets, and neutrophils, only 6 were identified in our proteomic datasets (plasma membrane calcium transporting ATPase and tyrosine-protein kinase SYK from monocytes, JAM-A, Ras-related protein Rab-6B, and Syntaxin-11 from platelets, and cathelicidin antimicrobial peptide precursor from neutrophils.) Additionally, we did not detect the transferrin receptor, a known reticulocyte marker. Furthermore, contaminant proteins were not found to be significantly changed, indicating they were uniformly distributed across different sample preparations.
The DBA patients analyzed in this study harbor different ribosomal protein mutations, with patients D3 and D4 on prednisone treatment. Patients D1 and D2 were without treatment at the time of this study, however they still present with macrocytic anemia (Table 1) and have received transfusions in the past. Comparison of results with the ribosomal mutations and patient hematological status did not reveal any obvious trends in protein abundance as a function of treatment or mutation. Only one protein, TAP2, was significantly increased in patients D3 and D4 only, indicating that the majority of proteins were not induced from steroid treatment alone. Interestingly, two patients with RPL5 mutations showed very high levels of dysferlin compared to the other DBA patients. However, more patient samples with differing ribosomal protein mutations are needed for definitive analysis.
The most surprising result is the identification of dysferlin on the RBC membrane of DBA patients because it has never been reported to be present on the RBC and it does not have a welldefined role in hematopoiesis. Dysferlin is a 230 kDa protein with a type two single-pass membrane spanning domain and seven calcium-binding C2 domains. It is normally highly expressed in muscle and mutations in this protein can lead to adult-onset muscular dystrophies, including Limb Girdle Muscular Dystrophy type 2B, [29] Myoshi Myopathy, [30] and Distal Anterior Compartment Myopathy. [31] Dysferlin is also found at lower expression levels in monocytes, endothelium, brain, kidney, and pancreas [32].
The presence of dysferlin in the RBC membrane appears to be relatively specific to DBA. Immunostaining for dysferlin in RBC ghosts showed variable levels of dysferlin present in the RBCs of all transfusion-independent DBA patients tested but it was not found on RBCs of transfusion dependent DBA patients (data not shown), indicating that dysferlin is an integral RBC membrane protein and does not derive from white blood cell contamination or from transfer to circulating RBCs by monocytes. Based on the accepted copy numbers of a-spectrin and band 3 per red cell (2.42610 5 and 1.2610 6 , respectively) [33] we can estimate the copy numbers of dysferlin in the DBA patients using the protein abundance values shown in Figure 2. These calculations yield approximate copy numbers of 1000, 4000, 2, and 100 for dysferlin in patients D1-D4, respectively. Screening other macrocytic blood disorders revealed the presence of dysferlin in the RBC of four of eight acquired aplastic anemia patients but not in any of the other 11 patients with various bone marrow failure disorders. Interestingly, one of these aplastic anemia patients has slightly higher eADA activity and a copy number neutral 5q loss of heterozygosity in the region of the RPS14 gene, raising the possibility of a functional connection between RP haploinsufficiency and dysferlin accumulation in RBCs. At this point, the potential roles of dysferlin on RBCs or in RBC maturation are mainly speculative, whether these are specific to the DBA pathology and impaired ribosome biogenesis, or whether the expression of dysferlin allows DBA erythropoiesis to produce circulating RBCs and rescue RBC production in the context of impaired ribosome biogenesis remains to be determined. The majority of studies on dysferlin function have been in the context of muscular dystrophy. In muscle, dysferlin was found to accumulate at the site of injury and is involved in membrane repair. [34] The normal function of dysferlin in monocytes is not well characterized, however, dysferlin deficient monocytes displayed increased phagocytosis. [35] A recent report showed that dysferlin interacts with components of the focal adhesion complex, implicating a role in cellular adhesion [36].
Several MHC class I proteins and MVP were present on the RBC membranes of DBA patients but not of healthy controls. MHC class I proteins are found on nucleated cells and are therefore not normally present on the surface of mature RBCs. However, several studies have found MHC class I molecules on the RBC of patients with systemic lupus erythematosus, [37,38] mononucleosis, [39] HIV, and interferon-a treatment. [40] Furthermore, a decrease in CR1 is also associated with systemic lupus erythematosus, further suggesting a link between DBA and an increased immune response. [41] Interestingly, a link between p53 and MVP has been described in which MVP was upregulated during p53 over-expression in human diploid fibroblasts and p53 was shown to directly bind to the MVP promoter. [42] Therefore, the increased activation of p53 in DBA erythroid progenitors [7] may explain the increase of MVP in DBA RBCs.
The results presented herein demonstrate that comparative proteome analysis of RBC membranes sensitively compares the abundance of a thousand RBC proteins and this unbiased approach revealed an unexpected expression of red cell membrane proteins on DBA RBC that would not have been discovered otherwise. The mechanism leading to the differences in the RBC proteome in DBA and the function of these proteins in the differentiation and maturation of DBA RBCs will be the subject of further investigation, which may be relevant to understanding DBA pathogenesis, diagnosis and treatment.

Supporting Information
File S1 Combined file of supporting figures and tables. Figure S1: Derivation of 95% confidence intervals for two sets of donor pools using normalized protein intensities. A) Plot of log10 standard deviation and mean protein intensities with exponential decay fit for donor pools C1 and C2 for comparison of patients D1 and D2, r 2 = 0.8665. B) Plot of log10 standard deviation and mean protein intensities with exponential decay fit for donor pools C3 and C4 for comparison of patients D3 and D4, r 2 = 0.8810. Figure  S2: Normalized log10 intensity plot of DBA patients D1-4 versus appropriate controls pools with a minimum detectable signal threshold of 1610 4 . Figure S3: Western blot analysis of RBC membrane preparations showing dysferlin and actin as a loading control for dysferlin-positive acquired aplastic anemia patients (labeled using their ages from Table 4 in the manuscript) for comparison to DBA patients D1-D4.