Validation and implementation of a method for microarray gene expression profiling of minor B-cell subpopulations in man
- Kim Steve Bergkvist1,
- Mette Nyegaard1, 2,
- Martin Bøgsted1, 3, 11,
- Alexander Schmitz1,
- Julie Støve Bødker1,
- Simon Mylius Rasmussen1,
- Martin Perez-Andres1,
- Steffen Falgreen1,
- Anders Ellern Bilgrau1, 3,
- Malene Krag Kjeldsen1,
- Michael Gaihede4, 11,
- Martin Agge Nørgaard5,
- John Bæch6,
- Marie-Louise Grønholdt7,
- Frank Svendsen Jensen8,
- Preben Johansen9,
- Karen Dybkær1, 11 and
- Hans Erik Johnsen1, 10, 11Email author
© Bergkvist et al.; licensee BioMed Central Ltd. 2014
Received: 25 October 2013
Accepted: 28 January 2014
Published: 31 January 2014
This report describes a method for the generation of global gene expression profiles from low frequent B-cell subsets by using fluorescence-activated cell sorting and RNA amplification. However, some of the differentiating compartments involve a low number of cells and therefore it is important to optimize and validate each step in the procedure.
Normal lymphoid tissues from blood, tonsils, thymus and bone marrow were immunophenotyped by the 8-colour Euroflow panel using multiparametric flow cytometry. Subsets of B-cells containing cell numbers ranging from 800 to 33,000 and with frequencies varying between 0.1 and 10 percent were sorted, subjected to mRNA purification, amplified by the NuGEN protocol and finally analysed by the Affymetrix platform.
Following a step by step strategy, each step in the workflow was validated and the sorting/storage conditions optimized as described in this report. First, an analysis of four cancer cell lines on Affymetrix arrays, using either 100 ng RNA labelled with the Ambion standard protocol or 1 ng RNA amplified and labelled by the NuGEN protocol, revealed a significant correlation of gene expressions (r ≥ 0.9 for all). Comparison of qPCR data in samples with or without amplification for 8 genes showed that a relative difference between six cell lines was preserved (r ≥ 0.9). Second, a comparison of cells sorted into PrepProtect, RNAlater or directly into lysis/binding buffer showed a higher yield of purified mRNA following storage in lysis/binding buffer (p < 0.001). Third, the identity of the B-cell subsets validated by the cluster of differentiation (CD) membrane profile was highly concordant with the transcriptional gene expression (p-values <0.001). Finally, in normal bone marrow and tonsil samples, eight evaluated genes were expressed in accordance with the biology of lymphopoiesis (p-values < 0.001), which enabled the generation of a gene-specific B-cell atlas.
A description of the implementation and validation of commercially available kits in the laboratory has been examined. This included steps for cell sorting, cell lysis/stabilization, RNA isolation, RNA concentration and amplification for microarray analysis. The workflow described in this report will enable the generation of microarray data from minor sorted B-cell subsets.
KeywordsMicroarray gene expression profiling RNA purification B-cell subpopulations Fluorescence activated cell sorting
Haematological malignancies are characterised by a continuous sub-clonal selection that becomes abnormally and non-homogeneously distributed within individual tumours. In most cases, malignant transformation and metastases are clonal since they are derived from single cells that, at least initially, preserve many features of the hierarchical structure of the normal tissue of origin. However, this has not been analysed systematically for malignancies [1–3]. Global gene expression profiling (GEP) and genetic alterations have already resulted in changes in the classification of malignant B-cell disorders [4–8], including identification of the normal cell of origin. Consequently, highly pure normal subpopulations are essential in order to investigate the complex cellular and molecular mechanisms involved in the stepwise B-cell differentiation in normal and malignant conditions [9–11].
The molecular mechanisms that control B-cell lymphopoiesis are regulated by the coordinated activity of a group of so-called master regulatory transcription factors (TF) . TFs can be divided into those that maintain B-cell lineage commitment (such as PAX5) and the TFs involved in the germinal centre (GC) reaction (such as BACH2 and BCL6), and those that promote and facilitate end-stage differentiation, notably IRF4, PRDM1, and XBP1. This ensures the separation of a range of well-defined subsets including pre-BI cells, pre-BII cells, immature (I), naive (N), centrocytes (CC), centroblasts (CB), memory B-cells (M), plasmablasts (PB) and end-stage antibody producing plasma cells (PC) [13–15].
A limiting factor in studies of B-cell subsets, however, is that some of the differentiating compartments involve a low number of cells and it is therefore important to optimize and validate each step in the procedure for global GEP of minor B-cell subsets following fluorescence-activated cell sorting (FACS). In addition, the conventional labelling kits require an input of 50 – 500 ng of total RNA, which corresponds to around 105 to 106 cells, depending on cell type. However, it is practically impossible to obtain such numbers from some B-cell subpopulations as they only constitute between 0.1 and 10% of the lymphoid tissue .
Early human B-cell development has been characterized by GEP, using FACS purified B-cell subsets in bone marrow (BM) [17, 18]. However, whereas the early B-cell differentiation has been characterized, the simultaneous sorting of post-germial centre B-cells in lymphoid tissue has been explored in less detail. The concept behind the present project is that a detailed workflow for the generation of GEP from minor B-cell subsets will allow us to establish a B-cell specific gene atlas. This will increase our knowledge of the B-cell differentiation and, ultimately, to use these gene lists in post-GC disease classification. The aims of the study were to validate and implement a fast and efficient method for isolation and generation of GEP from B-cell subsets in peripheral blood, tonsils, thymus and BM in order to generate a gene-specific B-cell atlas. The strategy omits an immunomagnetic purification step for B-cell enrichment before FACS, as often performed [17–20], which is problematic when low frequent B-cells are sorted.
A flow chart of the established protocol and methods for sorting the B-cell subsets in different tissue [see Additional file 1].
Impact of amplification
Gene specific amplification determined by qPCR
Total RNA from one million cells was extracted from six multiple myeloma (MM) cancer cell lines (CCLs) (MOLP-8, KMS-12-BM, RPMI-8226, OPM-2, LP-1, and KMM-1) using RNAeasy Plus Micro equipment (QIAGEN, Hilden, Germany). Genomic DNA was removed using gDNA Eliminator Spin Columns (QIAGEN, Hilden, Germany). RNA quality was evaluated with an Agilent 2100 bioanalyzer (Agilent Technologies, Inc., Palo Alto, CA) (RIN > 9.8). Total RNA from each CCL was processed in parallel by either directly converting 500 ng to cDNA (non-amplified) with SuperScript III First-Strand Synthesis Supermix (Invitrogen, Paisley, UK) or by amplifying 5 ng with an Ovation Pico WTA system (NuGEN Technologies, Inc., San Carlos, CA), as described by the manufacturer. QPCR assays were performed by comparing amplified cDNA at 25 ng/reaction to non-amplified cDNA derived from SuperScript III at 25 ng/reaction (total RNA equivalents). Commercially available Taqman primer probes sets, previously described in the qPCR Section, were used.
Comparing NuGEN protocol to standard protocol
Total RNA from one million cells of the same four CCLs KMM-1, OPM-2, U2932_M, and SU-DHL-5 was subjected to the NuGEN protocol or to the standard protocol from Ambion (Ambion WT Expression kit, Ambion, Inc., Austin, TX) following the manufactures recommendations. The input of total RNA was 1 ng for the NuGEN protocol whereas the input for the Ambion protocol was 100 ng.
Optimisation of storage buffer
Selection of storage buffer for sorted cells
The storage buffer was examined by sorting 15,000 fresh naive tonsil cells from a single donor directly into 12 separate tubes containing 450 μl of either lysis/binding buffer, RNAlater (Ambion, Austin, TX) or PrepProtect. mRNA was isolated using the μMACS™ technology (Miltenyi Biotech, Bergisch-Gladbach, Germany), allowing isolation on μ Columns and elution with pure water. This technology is referred to as magnetic bead isolation (MBI). MBI purification from sorted cells was performed in triplicates, either directly after cell sorting or after 14 days of storage of cells in the various RNA extraction buffers at 4°C, -20°C or -80°C. Before purification, all samples were equilibrated to RT and cells in RNAlater and PrepProtect was recovered from the storage solution by a 5 minute centrifugation (RT) at 5000 g and re-suspended in 450 μl of lysis/binding buffer, following MBI purification. Elutes were stored at -80°C before examining the yields with the TaqMan pre-developed endogenous control assay PPIA (333763 F) by RT-qPCR.
Ethical statement and tissue preparation
All samples were collected following informed consent in accordance with the research protocol accepted by the Ethics Committee for the North Denmark Region (N-20080062MCH). The cells were either FACS sorted fresh (i.e. samples collected, processed and sorted within the same day), or vital cryopreserved for storage and thawed before sorting.
Isolation of cells from tissues
The isolation strategy of all tissues is described in the following section, including the number of donors (n) for each tissue and the frequency of sorted cells used for generating GEP and establishing a B-cell atlas from different B-cell subsets in tonsils and BM.
Peripheral blood mononuclear cells (PBMNC) (frequency of sorted cells: 4,500 – 130,000)
PBMNC (n = 6 ) were isolated by diluting the peripheral blood sample 1:1 in phosphate-buffered saline (PBS) and placing the sample in a LeucoSep tube (Greiner Bio-One, Frickenhausen, Germany), according to the manufacturer’s instructions. The mononuclear cells (MNCs) were washed once and the red blood cells were lysed by adding 9 ml of Easylyse (DAKO, Glostrup, Denmark) to the pellet and incubating the sample for 15 minutes at room temperature (RT). The MNCs were washed twice in PBS and sorted fresh.
Tonsils (frequency of sorted cells: 5,600 – 20.000)
Tonsils (n = 8) were obtained by routine tonsillectomy as previously described . In brief, tonsils were placed on ice until homogenization in cold RPMI medium 1640 (Gibco, Invitrogen, UK) using a Medimachine (cod. 79200, Dako, Glostrup, Denmark) with a 35 μM sterile medicon. The cell suspension was passed through a 40 μM filter to remove debris and aggregates. The cell suspension was diluted in PBS before MNCs were isolated using Ficoll-Paque Plus (GE Health Care, Uppsala, Sweden) according to the manufacturer’s instructions. The cells were vital cryopreserved.
Thymus (frequency of sorted cells: 1,000 – 33,000)
Thymus (n = 7) was obtained from patients undergoing cardiac surgery and placed on ice immediately after removal and cut into smaller pieces. BM was simultaneously obtained from the same patient and processed as described in the following section. The Thymus tissue was squeezed and chopped with a tweezer in PBS. The cells were washed once in PBS and MNCs were purified using Ficoll-Paque Plus and sorted fresh.
BM (frequency of sorted cells: 800 – 25,000)
BM (n = 7) from sternum was obtained from patients undergoing cardiac surgery by physical scraping and scooping of the sternum and placed on ice immediately after removal. The BM was homogenized in 2 ml PBS using a syringe. The red blood cells were lysed by adding 20 ml of Easylyse (DAKO, Glostrup, Denmark) and incubating the sample for 45 min at RT. The samples were washed once in PBS before passing through a 40 μM filter to remove debris and aggregates and sorted fresh.
Human malignant B-cell lines
CCLs originating from MM MOLP-8, KMS-12-BM, RPMI-8226, LP-1, OPM-2 and diffuse large B-cell lymphoma (DLBCL) SU-DHL-5 were purchased from DSMZ [German Collection of Microorganisms and Cell Cultures, Braunschweig, Germany], whereas the MM CCL KMM-1 was purchased from JCRB [Japanese Collection of Research Bioresources Cell Bank, Japan]. DLBCL CCL U2932_M was generously provided by Jose A Martinez-Climent [Molecular Oncology Laboratory, University of Navarra Pamplona, Spain].
Multiparametric flow cytometry
MNCs were washed twice in PBS and in the final wash 2% BSA was added. The cells were stained with a 6–8 colour panel of monoclonal antibodies (mAbs). Isolating the B-cell subsets from PBMMC, tonsils, thymus and BM had the following five mAbs in common: CD20 clone 2H7 conjugated with pacific blue (eBioscience, San Diego, CA), CD45 clone 2D1 conjugated with anemonia majano cyan (BD Biosciences, San Jose, CA), CD10 clone HI10a conjugated with phycoerythrin/cyanin7 (BD Biosciences, San Jose, CA), CD27 clone L128 conjugated with Allophycocyanin (BD Biosciences, San Diego, CA), CD38 clone HIT2 conjugated with Alexa Flour 700 (ExBio Vestec, Czech Republic). CD19 conjugated with PERCPCy5.5 (BD Biosciences, San Diego, CA) were included in all samples except in tonsils. In BM, CD34 conjugated with phycoerythrin (BD Biosciences, San Diego, CA) were included. In tonsils, CD3 clone SK7 conjugated with fluorescein isothiocyanate (BD Biosciences, San Jose, CA), CD44 clone IM7 conjugated with peridinin chlorophyll protein/cyanin 5.5 (eBiosciences, San Diego, CA), and CXCR4 clone 12G5 conjugated with phycoerythrin (Beckman coulter, Brea, CA) were included. In thymus, sIgM conjugated with FITC (DAKO, Carpinteria, US) and sIgG conjugated with phycoerythrin (SouthernBiotech, Alabama, US) were included.
CD marker combinations used to distinguish between the B-cell subsets [see Additional file 1]
Briefly, all cells but tonsils were incubated for 30 minutes at room temperature (RT) (tonsils on ice) in darkness. Sorting of B-cell subsets was performed using a FACSAria2 cell sorter (BD Biosciences, San Jose, CA) at RT. Compensation was automatically calculated by FACSDiva software using single-stained control samples with the mAbs previously listed. Immediately before acquisition, the cells were filtered through a 35-μm filter (Cell Stainer, BD Biosciences). The purity of the isolated B-cell subsets (>90%) was confirmed by sorting approximately 1,000 cells into PBS and reacquisition of the sorted B-cell subsets. The cells were sorted into 450 μl of lysis/binding buffer (Miltenyi Biotech, Bergisch-Gladbach, Germany), except for the first PBMNC processed, which were sorted into 1 ml of PrepProtect (Miltenyi Biotech, Bergisch-Gladbach, Germany). The sorted B-cell subsets were stored at -20°C.
Global GEP of sorted B-cell subsets
Gene expression analysis of FACS sorted cells was performed using the Gene Chip Human Exon 1.0 ST (Exon) or the Gene chip HG U133 Plus 2.0 (U133) arrays (Affymetrix, Santa Clara, CA). In brief, mRNA was isolated by MBI and the eluted mRNA was concentrated to 5–10 μl by a volume reduction step using a speedVac Concentrator 5310 (Eppendorf, Hamburg, Germany). Five μl was used as input for amplification with the Ovation Pico WTA system (NuGEN Technologies, Inc., San Carlos, CA) following the manufacturer’s protocol. Due to technical challenges, it was not possible to determine the concentration of mRNA derived from five μl of elute. Finally, the samples were hybridised to the Exon or U133 array. This is referred to as the NuGEN protocol.
Reverse transcription and qPCR
Complementary DNA synthesis
Five μl of mRNA from sorted cells or 500 ng of total RNA from CCLs were used in cDNA synthesis using SuperScript III First-Strand Synthesis Supermix (Invitrogen, Paisley, UK). Random hexamers (50 μM) and oligo(dT) (50 μM) primers were used in a final volume of 20 μl.
qPCR was performed on a LightCycler 480 II using the LightCycler 480 Probes Masters PCR mix (Roche Diagnostics, Hvidovre, Denmark) and Taqman gene expression assays (Applied Biosystems, Foster City, CA) in a final reaction volume of 20 μl. Water controls, no-RT samples, and inter-run calibrator samples were included on all plates. Commercially available Taqman primer/probe sets were used (Applied Biosystems, Foster City, CA) for three endogenous control genes: PPIA [4333763 F], GAPDH [4333764 F] and TBP [4333769 F], for three transcription factors involved in the late B-cell lymphopoiesis: IRF4 [Hs01056534_m1], PRDM1 [Hs00153357_m1] and XBP1 [Hs00231936_m1] and two genes related to the oncogenesis of MM: MGST1 [Hs00220393_m1] and WHSC1 alias MMSET [Hs00983716_m1].
All data processing and statistical analysis was performed using Affymetrix GeneChip Command Console Software (AGCC), Partek Genomics Suite version 6.5 (Partek Inc., St. Louis, MO, USA), built-in Excel macros, and the statistical software system R, version 2.15.1. For all tests, p-values below 5% (p < 0.05) were considered statistically significant.
Gene specific amplification determined by qPCR
For a specific gene, x, the difference in Cq values (dCq) between two amplified CCLs, 1 and 2, (Cqx,1 – Cqx,2)amplified, was plotted against the dCq-value for the same two non-amplified cell lines (Cqx,1 – Cqx,2)non-amplified. Pearson’s correlation coefficient between the two dCq values was calculated. A significant positive Pearson’s correlation coefficient indicates preservation of differential gene expression after amplification . In addition, for each gene, the six cell lines were ranked from high to low expression. This ranking was performed for both the non-amplified and the amplified cell lines, and compared using Spearman’s rank correlation. A test for inconsistent ranking was carried out by an exact permutation test.
Selection of storage buffer for sorted cells
The Cq values were analysed by a 2-way analysis of variance (ANOVA) with the two factors storage buffer at three levels (Lysis buffer; PrepProtect; RNAlater) and condition at four levels (Directly; 4°C; -20°C; -80°C).
Data processing and analysis
CEL-files from the two array types U133 and Exon were generated by AGCC [for results of the initial quality control analyses, see Additional file 2]. CEL-files were imported into Partek and RMA normalised. Due to in-homogenous variance, the difference in expression values between the B-cell subsets for each of the endogenous control genes was tested by a Kruskal-Wallis one-way analysis of variance. Furthermore, for biological validation, the expression values of selected genes with well-known function in each B-cell subset of the BM were analysed. Significant expressions between functional classes of B-cell subsets were tested by a two-tailed t-test with unequal variance. Correlations between two measurements were evaluated throughout by the Pearson’s correlation coefficient and denoted by an r. Finally, plots of gene differences versus gene averages, MA plots, were used to examine the concordance between methods for generating global GEP.
Concordance between the pre-defined CD marker and array based transcript expressions
The B-cell subsets were divided into two groups based on the pre-defined positive or negative CD marker. The mean and standard deviations were calculated for the positive and negative groups based on the gene expression value for the CD marker. Discordance was registered if a gene expression value for a B-cell subset was observed in the opposite group as defined by the CD marker. Fisher’s exact test for 2x2 tables was used to test for independence between the groupings based on CD marker and GEP [see Additional file 3].
Freshly sorted samples or cryopreservation before sorting
The ratios of excluded samples were compared for freshly sorted samples as well as for those cryopreserved before sorting. Samples were either excluded due to low amplification yield or unsatisfactory quality control parameters on array. The significance between the ratios was tested by a chi-square test for equality of proportions.
Impact of amplification
Gene specific amplification determined by qPCR
Comparing NuGEN protocol to standard protocol
Amount of starting material and percent present calls on Exon arrays
Amount of aRNA/cDNA after processing
The reproducibility between the two protocols was first examined by the narrow set of previously used genes [see Additional file 6]. The expression of the selected genes generated by the NuGEN protocol was plotted against the expression of the genes generated by the Ambion protocol as presented in Figure 1B. We observed a high degree of correlation between the two protocols regarding the expression from the eight genes (r = 0.950). Next, we measured the reproducibility of the gene expression on a global scale by generating MA plots [see Additional file 7]. We observed a Pearson’s correlation coefficient ranging from r = 0.892 to r = 0.904 for all pair-wise comparisons of protocols, demonstrating an acceptable reproducibility and robustness of the NuGEN protocol.
Optimisation of storage buffer
Selection of storage buffer for sorted cells
A priori, the handling procedures were considered important, including methods for RNA extractions [25–29]. In one previous study, the best method for RNA extraction from flow-sorted cells was MBI , which was selected for the present study. A pilot study revealed that the recovery of relative few cells (50,000-200,000) after centrifugation varied by up to 50% (Unpublished observations). Cells sorted in conventional storage buffers such as PrepProtect or RNAlater include a centrifugation step in order to re-suspend the cells in lysis buffer. Therefore, we wanted to test if the lysis/binding buffer could be used as a storage buffer, thereby omitting the centrifugation step in the procedure.
Performance of the protocol
Amplified yield from FACS tonsil B-cell subsets
The technical variation in the protocol, including flow sorting, freezing, purifying mRNA, up-concentrating RNA and amplifying the samples, was addressed by sorting two B-cell subsets, namely the N and the PB from tonsils in triplicates from a single donor. By sorting different numbers of cells (500, 1,000 and 5,000 PBs and 1,000, 5,000 and 10,000 N cells), we found that the amplified cDNA yield per cell varied depending on cell type, with PBs providing a higher yield compared to N cells (Figure 2B). Based on the amplified cDNA yield, 1,000 PB and 5,000 N was sufficient for generating global GEP, according to the input requirement in the protocol. These samples were only intended to be used to investigate the technical variation in the protocol and not used to generate global GEP.
Tissue-independent yield and QC on microarray data
Tissue independent yield and QC on microarray
Sorted cells (x1000)
Total no of Amp. samples
Amp. yield cDNA μg
PBMNC (F) n = 6
5 F; 1 M
49 ± 13
5.3 ± 0.9
6.5 ± 0.9
7.1 ± 0.6
5.3 ± 1.3
PBMNC (F) n = 3
1 F; 2 M
43 ± 17
6.6 ± 1.0
7.1 ± 1.1
8.6 ± 1.5
9.3 ± 0.3
PBMNC (C) n = 3
6.0 ± 1.5
7.3 ± 1.1
6.4 ± 1.2
9.8 ± 0.5
Tonsils (C) n = 8
6 F; 2 M
17 ± 10
6.8 ± 0.8
8.8 ± 0.9
8.8 ± 0.9
6.7 ± 0.9
9.3 ± 0.6
Thymus (F) n = 7
2 F; 5 M
63 ± 12
7.6 ± 2.2
6.4 ± 3.1
8.4 ± 3.1
7.2 ± 4.1
BM (F) n = 7
7.8 ± 2.2
7.8 ± 1.2
6.0 ± 1.7
6.3 ± 2.2
7.1 ± 2.5
8.0 ± 2.2
Concordance between the predefined CD marker and array-based transcript expressions
A quality control step was included to validate the identity of the sorted B-cell subsets by correlating the pre-defined surface expressed CD markers used for FACS of the B-cell subsets to the CD markers transcript expression levels on microarray [see Additional file 3]. In summary, there was a highly significant (p < 0.001) concordant expression of CD protein markers and transcript expressions. For example, in a total of 44 PBMNC samples concordance was observed for CD20, CD10, CD27, and CD38 gene expression in 43, 43, 42, and 41 samples respectively [see Additional file 3: Table S2].
Global GEP data were generated from sorted B-cell subsets in PBMNC and lymphoid tissues. The initial data analysis was carried out to explore similarities and differences between samples and to determine whether the samples could be grouped into distinct B-cell subsets. The data matrix was subjected to principal component analysis (PCA) and the results are shown in the Additional file 8. It is noted that the preBI and II; immature, naive and memory B-cells; plasmablasts and plasma cells tended to cluster in separate biological relevant groups, as expected. In addition, naive and memory B-cells tended to cluster together, regardless of tissue, as did the centroblasts and centrocyte populations in tonsils, regardless of array type. Of note, the tonsil donors on the Exon and U133 arrays are not identical.
Likewise, in the tonsil data set, eight evaluated genes were expressed in accordance with the biology of lymphopoiesis (p-values < 0.001) [see Additional file 9]. The genes included the formerly evaluated transcription factors PAX5, PRDM1 and IRF4. In addition, four genes (BCL6, AICDA, BACH2 and CXCR4) involved in the GC reaction and the proliferation marker MKI67 were included. Finally, with the implemented and validated protocol, a gene specific B-cell atlas was generated from the BM and tonsil tissue [see Additional files 10 and 11].
The aim of this project was to implement and validate the steps in a protocol for generating global GEP of FACS sorted low frequent B-cell subsets in lymphoid tissue and blood. In the present study, we have validated a protocol for microarray studies of minor subsets and sorted cells from different tissues in the range 800 to 33,000 cells into lysis/binding buffer.
In order to use GEP for the analysis of populations below 100,000 cells, RNA has to be reverse-transcribed and amplified. The amplification technologies have been recently evaluated with respect to reproducibility and sensitivity and it was observed that the NuGEN protocol was the most suitable for amplification of pico amounts of RNA . In addition, a study comparing different amplification protocols, including NuGEN protocol, Message Amp (Ambion), Small Sample Target Labelling Assay Version II (Affymetrix) and BioArray Small Sample Amplification Protocol (Enzo Life Sciences) found that all technologies performed acceptably. The NuGEN protocol, however, resulted in the most sensitive and specific data when 10 ng of RNA was amplified .These studies were the basis for choosing the NuGEN protocol in this study. The amplification by the NuGEN protocol was first evaluated by comparing the expression of eight genes in six non-amplified and amplified CCLs. We noticed a gene-specific amplification, meaning that some sequences or parts of transcripts were amplified better than others. However, differential gene expression was preserved across the CCLs (r = 0.917), which is in line with previous findings . Next we compared the NuGEN protocol to the standard protocol from Ambion by examining the generated global gene expression data on the Exon array. The reproducibility was good, both when the expression was compared globally (r = 0.892 to 0.904) and compared for a narrow set of genes (r = 0.950). Correlation coefficients between methods are usually lower than 0.90 when compared .
By sorting cells directly into lysis/binding buffer, a significantly higher yield was obtained both for direct and postponed time for preparation of mRNA purification. The reduced yields obtained from cell storage in RNAlater are very likely due to the centrifugation step introduced to recover the cells prior to lysis, but do not explain the higher Cq values of approximately 3–4 obtained in PrepProtect. Of notice, reducing the centrifugation force to 3000 g the same Cq values was obtained for cells stored in PrepPretect; however, any inhibition in the enzymatic assay has not been tested, which also may explain the poor yield.
By using the NuGEN amplification method and MBI technology, we obtained sufficient amplified cDNA yields from 1,000 PB and 5,000 N B-cells to ensure successful microarray analysis.
As the present work represents one of the first global transcriptome studies with sorted low frequent B-cell subsets, the reliability of the results obtained has been considered a critical issue. Certainly, contamination of the sorted B-cell subsets with “rosetting” non B-cells cannot be excluded. However, it is evident from the results obtained in the present analysis that such contaminations do not compromise gene expression analysis. First, the CD markers used for identification and sorting are concordantly expressed at the transcript level [see Additional file 5]. Secondly, biological validation of the generated GEP was conducted in the BM and tonsil samples for eight genes and we found them to be expressed in accordance with our knowledge of the B-cell lymphopoiesis.
A protocol for microarray analysis of FACS-sorted low frequent B-cells was implemented and validated to achieve new insights into B-cell differentiation. The protocol consisted of CD panels for sorting early/mature and post-GC B-cells in BM, tonsils, thymus and blood; sorting cells directly into lysis/binding buffer; purifying mRNA by MBI; and concentrating RNA by speedVac concentrator following amplification with NuGEN’s technology. This protocol represents a significant advance over established protocols by allowing the possibility to sort the cells directly into lysis/binding buffer, which also can be used for long-term storage. In addition, complete CD panels for the simultaneous sorting of pre-GC, GC and post-GC B-cells in lymphoid tissues and blood is provided. A gene specific B-cell atlas was generated in BM and tonsils and a future goal is to assign these profiles in post-GC disease classification by “cell of origin”.
Global gene expression profiling
Following fluorescence-activated cell sorting
Peripheral blood mononuclear cells
Cancer cell lines
Diffuse large B-cell lymphoma
Gene Chip Human Exon 1.0 ST
Gene chip HG U133 Plus 2.0
Magnetic bead isolation
Affymetrix GeneChip Command Console Software
Pearson’s correlation coefficient
Analysis of variance
Cluster of differentiation
in vitro transcription
Principal component analysis.
The authors would like to acknowledge Ann-Maria Jensen, Helle Høholt and Helle Stiller, Department of Haematology, Aalborg Hospital Science and Innovation Centre, Aalborg University Hospital, Aalborg, Denmark for their expert technical assistance.
This project has been supported by the European Myeloma Network by means of a travel grant to MPA - an independent organization supporting research in Europe.
KSB, MKK, MPA and HEJ were supported by the EU 6th FP to MSCNET (LSHC-CT-2006-037602), the Danish Cancer Society, the Danish Research Agency (#2101-07-0007), and the KE Jensen Foundation (2006–2008).
- Greaves M, Maley CC: Clonal evolution in cancer. Nature. 2012, 481: 306-313. 10.1038/nature10762.PubMedPubMed CentralView ArticleGoogle Scholar
- Valent P, Bonnet D, De MR, Lapidot T, Copland M, Melo JV: Cancer stem cell definitions and terminology: the devil is in the details. Nat Rev Cancer. 2012, 12: 767-775. 10.1038/nrc3368.PubMedView ArticleGoogle Scholar
- Johnsen HE, Kjeldsen MK, Urup T, Fogd K, Pilgaard L, Boegsted M: Cancer stem cells and the cellular hierarchy in haematological malignancies. Eur J Cancer. 2009, 45 (Suppl 1): 194-201.PubMedView ArticleGoogle Scholar
- Hoffmann MH, Klausen TW, Boegsted M, Larsen SF, Schmitz A, Leinoe EB: Clinical impact of leukemic blast heterogeneity at diagnosis in cytogenetic intermediate-risk acute myeloid leukemia. Cytometry B Clin Cytom. 2012, 82: 123-131.PubMedView ArticleGoogle Scholar
- Chiorazzi N, Ferrarini M: Cellular origin(s) of chronic lymphocytic leukemia: cautionary notes and additional considerations and possibilities. Blood. 2011, 117: 1781-1791. 10.1182/blood-2010-07-155663.PubMedPubMed CentralView ArticleGoogle Scholar
- Rosenwald A, Wright G, Chan WC, Connors JM, Campo E, Fisher RI: The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma. N Engl J Med. 2002, 346: 1937-1947. 10.1056/NEJMoa012914.PubMedView ArticleGoogle Scholar
- Rasmussen T, Jensen L, Johnsen HE: The clonal hierachy in multiple myeloma. Acta Oncol. 2000, 39: 765-770. 10.1080/028418600750063479.PubMedView ArticleGoogle Scholar
- Kuppers R: The biology of Hodgkin’s lymphoma. Nat Rev Cancer. 2009, 9: 15-27. 10.1038/nrc2542.PubMedView ArticleGoogle Scholar
- Kjeldsen MK, Perez-Andres M, Schmitz A, Johansen P, Boegsted M, Nyegaard M: Multiparametric flow cytometry for identification and fluorescence activated cell sorting of five distinct B-cell subpopulations in normal tonsil tissue. Am J Clin Pathol. 2011, 136: 960-969. 10.1309/AJCPDQNP2U5DZHVV.PubMedView ArticleGoogle Scholar
- Shaffer AL, Rosenwald A, Staudt LM: Lymphoid malignancies: the dark side of B-cell differentiation. Nat Rev Immunol. 2002, 2: 920-932. 10.1038/nri953.PubMedView ArticleGoogle Scholar
- Rasmussen T, Haaber J, Dahl IM, Knudsen LM, Kerndrup GB, Lodahl M: Identification of translocation products but not K-RAS mutations in memory B cells from patients with multiple myeloma. Haematologica. 2010, 95: 1730-1737. 10.3324/haematol.2010.024778.PubMedPubMed CentralView ArticleGoogle Scholar
- Matthias P, Rolink AG: Transcriptional networks in developing and mature B cells. Nat Rev Immunol. 2005, 5: 497-508. 10.1038/nri1633.PubMedView ArticleGoogle Scholar
- Morgan GJ, Walker BA, Davies FE: The genetic architecture of multiple myeloma. Nat Rev Cancer. 2012, 12: 335-348. 10.1038/nrc3257.PubMedView ArticleGoogle Scholar
- O’Neil J, Look AT: Mechanisms of transcription factor deregulation in lymphoid cell transformation. Oncogene. 2007, 26: 6838-6849. 10.1038/sj.onc.1210766.PubMedView ArticleGoogle Scholar
- Welner RS, Pelayo R, Kincade PW: Evolving views on the genealogy of B cells. Nat Rev Immunol. 2008, 8: 95-106. 10.1038/nri2234.PubMedView ArticleGoogle Scholar
- Caraux A, Klein B, Paiva B, Bret C, Schmitz A, Fuhler GM: Circulating human B and plasma cells. Age-associated changes in counts and detailed characterization of circulating normal C. Haematologica. 2010, 95: 1016-1020. 10.3324/haematol.2009.018689.PubMedPubMed CentralView ArticleGoogle Scholar
- Hystad ME, Myklebust JH, Bo TH, Sivertsen EA, Rian E, Forfang L: Characterization of early stages of human B cell development by gene expression profiling. J Immunol. 2007, 179: 3662-3671.PubMedView ArticleGoogle Scholar
- van Zelm MC, van der Burg M, de Ridder D, Barendregt BH, de Haas EF, Reinders MJ: Ig gene rearrangement steps are initiated in early human precursor B cell subsets and correlate with specific transcription factor expression. J Immunol. 2005, 175: 5912-5922.PubMedView ArticleGoogle Scholar
- Klein U, Tu Y, Stolovitzky GA, Keller JL, Haddad J, Miljkovic V: Transcriptional analysis of the B cell germinal center reaction. Proc Natl Acad Sci USA. 2003, 100: 2639-2644. 10.1073/pnas.0437996100.PubMedPubMed CentralView ArticleGoogle Scholar
- Nakayama Y, Stabach P, Maher SE, Mahajan MC, Masiar P, Liao C: A limited number of genes are involved in the differentiation of germinal center B cells. J Cell Biochem. 2006, 99: 1308-1325. 10.1002/jcb.20952.PubMedView ArticleGoogle Scholar
- Vermeulen J, Derveaux S, Lefever S, De SE, De PK, Yigit N: RNA pre-amplification enables large-scale RT-qPCR gene-expression studies on limiting sample amounts. BMC Res Notes. 2009, 2: 235-10.1186/1756-0500-2-235.PubMedPubMed CentralView ArticleGoogle Scholar
- Lauss M, Vierlinger K, Weinhaeusel A, Szameit S, Kaserer K, Noehammer C: Comparison of RNA amplification techniques meeting the demands for the expression profiling of clinical cancer samples. Virchows Arch. 2007, 451: 1019-1029. 10.1007/s00428-007-0522-4.PubMedView ArticleGoogle Scholar
- Eberwine J, Yeh H, Miyashiro K, Cao Y, Nair S, Finnell R: Analysis of gene expression in single live neurons. Proc Natl Acad Sci USA. 1992, 89: 3010-3014. 10.1073/pnas.89.7.3010.PubMedPubMed CentralView ArticleGoogle Scholar
- Phillips J, Eberwine JH: Antisense RNA Amplification: A linear amplification method for analyzing the mRNA population from single living cells. Methods. 1996, 10: 283-288. 10.1006/meth.1996.0104.PubMedView ArticleGoogle Scholar
- Deng MY, Wang H, Ward GB, Beckham TR, McKenna TS: Comparison of six RNA extraction methods for the detection of classical swine fever virus by real-time and conventional reverse transcription-PCR. J Vet Diagn Invest. 2005, 17: 574-578. 10.1177/104063870501700609.PubMedView ArticleGoogle Scholar
- Khymenets O, Ortuno J, Fito M, Covas MA, Farre M, de la Torre R: Evaluation of RNA isolation procedures from human blood and its application for gene expression studies (Sod-1, Sod-2). Anal Biochem. 2005, 347: 156-158. 10.1016/j.ab.2005.07.007.PubMedView ArticleGoogle Scholar
- Ruettger A, Neumann S, Wiederanders B, Huber R: Comparison of different methods for preparation and characterization of total RNA from cartilage samples to uncover osteoarthritis in vivo. BMC Res Notes. 2010, 3: 7-10.1186/1756-0500-3-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Rump LV, Asamoah B, Gonzalez-Escalona N: Comparison of commercial RNA extraction kits for preparation of DNA-free total RNA from Salmonella cells. BMC Res Notes. 2010, 3: 211-10.1186/1756-0500-3-211.PubMedPubMed CentralView ArticleGoogle Scholar
- Eldh M, Lotvall J, Malmhall C, Ekstrom K: Importance of RNA isolation methods for analysis of exosomal RNA: evaluation of different methods. Mol Immunol. 2012, 50: 278-286. 10.1016/j.molimm.2012.02.001.PubMedView ArticleGoogle Scholar
- Mack E, Neubauer A, Brendel C: Comparison of RNA yield from small cell populations sorted by flow cytometry applying different isolation procedures. Cytometry A. 2007, 71: 404-409.PubMedView ArticleGoogle Scholar
- Lossos IS, Czerwinski DK, Wechser MA, Levy R: Optimization of quantitative real-time RT-PCR parameters for the study of lymphoid malignancies. Leukemia. 2003, 17: 789-795. 10.1038/sj.leu.2402880.PubMedView ArticleGoogle Scholar
- Oettinger MA, Schatz DG, Gorka C, Baltimore D: RAG-1 and RAG-2, adjacent genes that synergistically activate V(D)J recombination. Science. 1990, 248: 1517-1523. 10.1126/science.2360047.PubMedView ArticleGoogle Scholar
- Bartholdy B, Matthias P: Transcriptional control of B cell development and function. Gene. 2004, 327: 1-23. 10.1016/j.gene.2003.11.008.PubMedView ArticleGoogle Scholar
- Cattoretti G, Shaknovich R, Smith PM, Jack HM, Murty VV, Alobeid B: Stages of germinal center transit are defined by B cell transcription factor coexpression and relative abundance. J Immunol. 2006, 177: 6930-6939.PubMedView ArticleGoogle Scholar
- Lu R, Medina KL, Lancki DW, Singh H: IRF-4,8 orchestrate the pre-B-to-B transition in lymphocyte development. Genes Dev. 2003, 17: 1703-1708. 10.1101/gad.1104803.PubMedPubMed CentralView ArticleGoogle Scholar
- Clement-Ziza M, Gentien D, Lyonnet S, Thiery JP, Besmond C, Decraene C: Evaluation of methods for amplification of picogram amounts of total RNA for whole genome expression profiling. BMC Genomics. 2009, 10: 246-10.1186/1471-2164-10-246.PubMedPubMed CentralView ArticleGoogle Scholar
- Viale A, Li J, Tiesman J, Hester S, Massimi A, Griffin C: Big results from small samples: evaluation of amplification protocols for gene expression profiling. J Biomol Tech. 2007, 18: 150-161.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.