Skip to main content
  • Research article
  • Open access
  • Published:

A map of gene expression in neutrophil-like cell lines

Abstract

Background

Human neutrophils are central players in innate immunity, a major component of inflammatory responses, and a leading model for cell motility and chemotaxis. However, primary neutrophils are short-lived, limiting their experimental usefulness in the laboratory. Thus, human myeloid cell lines have been characterized for their ability to undergo neutrophil-like differentiation in vitro. The HL-60 cell line and its PLB-985 sub-line are commonly used to model human neutrophil behavior, but how closely gene expression in differentiated cells resembles that of primary neutrophils has remained unclear.

Results

In this study, we compared the effectiveness of differentiation protocols and used RNA sequencing (RNA-seq) to compare the transcriptomes of HL-60 and PLB-985 cells with published data for human and mouse primary neutrophils. Among commonly used differentiation protocols for neutrophil-like cell lines, addition of dimethyl sulfoxide (DMSO) gave the best combination of cell viability and expression of markers for differentiation. However, combining DMSO with the serum-free-supplement Nutridoma resulted in increased chemotactic response, phagocytic activity, oxidative burst and cell surface expression of the neutrophil markers FPR1 and CD11b without a cost in viability. RNA-seq analysis of HL-60 and PLB-985 cells before and after differentiation showed that differentiation broadly increases the similarity in gene expression between the cell lines and primary neutrophils. Furthermore, the gene expression pattern of the differentiated cell lines correlated slightly better with that of human neutrophils than the mouse neutrophil pattern did. Finally, we created a publicly available gene expression database that is searchable by gene name and protein domain content, where users can compare gene expression in HL-60, PLB-985 and primary human and mouse neutrophils.

Conclusions

Our study verifies that a DMSO-based differentiation protocol for HL-60 and PLB-985 cell lines gives superior differentiation and cell viability relative to other common protocols, and indicates that addition of Nutridoma may be preferable for studies of chemotaxis, phagocytosis, or oxidative burst. Our neutrophil gene expression database will be a valuable tool to identify similarities and differences in gene expression between the cell lines and primary neutrophils, to compare expression levels for genes of interest, and to improve the design of tools for genetic perturbations.

Background

Neutrophils are the most abundant immune cell population circulating in the human bloodstream, representing 50–70% of all leukocytes. They are generated from myeloid precursors in the bone marrow, where they undergo several stages of maturation, namely myeloblast, promyelocyte, myelocyte, metamyelocyte, band cell and finally polymorphonuclear neutrophil [1]. From the bloodstream, neutrophils are rapidly recruited to sites of inflammation, migrating up gradients of chemical cues in a process called chemotaxis. Neutrophils are among the world’s most accurate chemotaxing cells, and serve as a leading model for eukaryotic chemotaxis [2]. Once at the inflammation site, neutrophils respond to the chemoattractant signals with cytotoxic and inflammatory responses, including secretion of granules (degranulation), ingestion of microbes and other particles (phagocytosis), and generation of neutrophil extracellular traps that catch and kill extracellular microbes. Neutrophil responses require a fine balance: underactivity leaves the body vulnerable to invasion, while over-activity is implicated in inflammatory and autoimmune diseases, and cancer progression [1, 3].

Despite much interest in understanding the molecular mechanisms underlying human neutrophil behaviors, their use in experimental research is limited by their very short life-span. They begin to undergo apoptosis within 6 to 12 h after isolation [1], making techniques such as genetic perturbation impractical. To address this issue, human myeloid cell lines have been established and characterized for their ability to undergo neutrophil-like differentiation in vitro [4, 5]. The HL-60 cell line, established in 1977 from an acute myeloid leukemia patient, is one of the most commonly used [6]. HL-60 cells are at the myeloblast stage of development but can be induced to differentiate terminally in vitro to either a neutrophil-like or a monocyte-like state. Neutrophilic differentiation causes a decrease in cell size, increased nuclear pyknosis and segmentation, and alterations in gene and protein expression. This in vitro differentiation makes HL-60 cells a good model to study neutrophil behaviors such as chemotaxis and phagocytosis [7,8,9], as well as the myeloid differentiation process itself [10]. A second commonly used cell line, PLB-985, was originally reported as being established in 1985 from a different patient [11]. However, further characterization determined that PLB-985 is actually a sub-line of HL-60 [12]. PLB-985 cells have similar, but slightly different, properties from those of HL-60 cells, and they are used to study processes including neutrophil ROS production [13] and chemotaxis [14,15,16]. These cell lines serve as primary model systems, but even so, they do not fully recapitulate all neutrophil behaviors [5, 17]. Important outstanding questions include to what degree the in vitro differentiated cells recapitulate neutrophil gene expression, and what differences in gene expression exist that may affect particular neutrophil behaviors.

Additionally, it is still not clear which of the commonly used reagents and protocols best induces differentiation of HL-60 and PLB-985 cells into a neutrophil-like state. Dimethyl sulfoxide (DMSO), all-trans retinoic acid (ATRA), and dibutyryl cyclic adenosine monophosphate (dbcAMP) are all commonly used to induce differentiation [5, 11, 18]. While ATRA and dbcAMP are thought to regulate transcription through cognate retinoic acid receptors and cAMP signaling, respectively, the mechanism for DMSO is much less clear. DMSO may act by altering membrane permeability and microviscosity, by affecting synthesis of nucleic acid and proteins, or by inducing changes in chromatin [19,20,21]. Additional reagents have also been reported to improve the efficiency of differentiation. In particular, addition of caffeic acid (CA) [22], replacement of serum with the supplement Nutridoma [23], and stimulation with the cytokine granulocyte colony-stimulating factor (G-CSF) [24] have each been used to promote neutrophilic differentiation in some contexts.

We wanted to address the above outstanding questions by first determining an optimal protocol for efficient neutrophilic differentiation, and then characterizing gene expression in differentiated cells for systematic comparison with primary neutrophils. Our study indicates that DMSO-based differentiation achieves the best balance of neutrophil marker expression with high cell viability and function, and that the use of Nutridoma can further increase cell surface expression of at least some neutrophil markers, including the Formyl Peptide Receptor 1 (FPR1). Our RNA sequencing (RNA-seq) analysis of HL-60 and PLB-985 transcriptomes confirmed that differentiation markedly increases the correlation in gene expression between the cell lines and human neutrophils, slightly exceeding the correlation between mouse and human neutrophil gene expression. Finally, we organized our data, together with previously published data, in a public neutrophil gene expression database that is searchable by gene name and by protein domain content. This database will serve as a resource for the community for easy comparison of gene expression in the cell lines versus human and mouse primary neutrophils, and as a reference for improving gene perturbation strategies in neutrophil models.

Results and discussion

Comparison of commonly used differentiation protocols

Since neutrophilic differentiation is often associated with increased cell death, we set out to determine which differentiation protocol produces the most highly penetrant differentiation while maintaining high cell viability. We compared the most commonly used differentiation agents in PLB-985 cells by measuring their ability to induce expression of neutrophil surface markers and reduce proliferation while maintaining cell viability. Initially, we compared differentiation for 6 days with DMSO, ATRA (with 0.5% N,N-Dimethyl formamide (DMF)) and dbcAMP [11, 25, 26].

To assess differentiation, we measured both early and late differentiation markers. We chose CD11b (Integrin Alpha M, ITGAM) as an early marker, since it begins to appear within 3 days of differentiation with DMSO [27], and FPR1 as a late neutrophil differentiation marker since it is usually only detectable after 5 or 6 days of differentiation. We measured both markers by flow cytometry, using a monoclonal antibody against CD11b and a fluorescent peptide ligand for FPR1 (FLPEP). An isotype control for the CD11b antibody showed no staining (data not shown). While undifferentiated cells did not express CD11b or FPR1, a substantial percentage of cells expressed both markers when treated with any of the three differentiation inducers (Fig. 1a, Additional file 1: a and b). Treatment with either DMSO or ATRA with DMF caused efficient expression of CD11b, while dbcAMP only weakly induced this early differentiation marker (Fig. 1a, Additional file 1: a). FPR1 was induced with partial penetrance by all differentiation agents. DMSO and dbcAMP gave the highest fractions of positive cells, but a sizeable fraction of cells with low or undetectable receptor levels also remained (Fig. 1b, Additional file 1: b, and Table 1). We also analyzed differentiation with ATRA and DMF independently, but their combination resulted in higher levels of both CD11b and FPR1 than either compound on its own (Additional file 1).

Fig. 1
figure 1

DMSO-based differentiation protocol gives the best differentiation potential based on expression of differentiation markers and cell viability. PLB-985 cells were differentiated into a neutrophil-like state by culturing in different media (DMSO, ATRA + DMF, or dbcAMP) for 6 days. Undifferentiated cells were also analyzed. Cells were stained with an antibody against CD11b, chosen as an early differentiation marker (a), with FLPEP (a fluorescent ligand of FPR1), chosen as a late differentiation marker (b), or with NucRed Dead 647 Probe, to measure cell death (c). Samples were measured by cytometry and data was analyzed using MATLAB. d Cell growth was monitored at days zero, three, and six of the differentiation protocols, using the trypan blue dye exclusion test. The number of cells were normalized to the initial number of cells. These experiments were performed three times and a representative experiment is shown

Table 1 PLB-985 cell death and expression of FPR1 under several differentiation inducer agents

We measured the loss of cell viability by staining cells with the NucRed Dead 647 reagent, which binds to DNA in cells with compromised plasma membrane integrity. While DMSO caused little if any apparent cell death, ATRA with DMF and dbcAMP each resulted in about 50% inviability in the cell population (Fig. 1c and Table 1). We also monitored cell proliferation by counting viable cells using the trypan blue exclusion method. All three treatments reduced proliferation, with dbcAMP acting the most rapidly (Fig. 1d). However, by day 6, the number of cells in the ATRA with DMF and dbcAMP conditions had markedly decreased, consistent with high rates of cell death. These results also coincide with reports from the initial characterization of the PLB-985 cell line [11].

While our initial results indicated that DMSO is the most suitable agent for inducing neutrophil marker expression without hampering cell viability, we wanted to explore variations of the protocol that might further improve differentiation efficiency. We pursued combinations of DMSO, ATRA [28] and CA [22], as these have been reported to improve differentiation. However, we observed increased cell death and marginal improvement of FPR1 expression at best compared to DMSO alone (Table 1). In an attempt to decrease cell death while still inducing differentiation [29], we limited induction with ATRA to 12 h followed by DMSO differentiation media for 5 days. However, the fraction of cells expressing FPR1 was no higher than for DMSO alone (data not shown).

Replacing serum with Nutridoma during differentiation increases FPR1 surface expression and chemotactic efficiency

We also tested the effects of serum and serum replacement on differentiation. High concentrations of serum have been reported to inhibit differentiation [23]. However, the cell death rate increases markedly under serum-free conditions [23]. Nutridoma, a serum free differentiation supplement, has been reported to improve cell differentiation while providing serum-like support of cell viability [23, 30]. To test how serum concentration and Nutridoma affect differentiation in PLB-985 cells, we analyzed upregulation of FPR1 upon differentiation in media containing DMSO, 0.5% fetal bovine serum (FBS), and 2% Nutridoma. Surprisingly, replacing serum with Nutridoma increased the percentage of cells that express FPR1 by approximately 20% compared to DMSO alone (Table 1 and Fig. 2a), while cell viability remained similar to that of non-differentiated cells (Table 1). Consistent with our FPR1 results, we observed that the presence of Nutridoma in the differentiation media also increased the surface expression of CD11b (Fig. 2b). We also validated the Nutridoma differentiation protocol in HL-60 cells, where we saw similar increases in CD11b and FPR1 surface expression (Additional file 2).

Fig. 2
figure 2

Replacing serum with Nutridoma during differentiation increases FPR1 surface expression and chemotactic efficiency. PLB-985 cells were differentiated into a neutrophil-like state by culturing in media supplemented with 1.3% DMSO and 9% FBS or supplemented with 1.3% DMSO, 2% Nutridoma and 0.5% FBS, for 6 days. Then cells were stained with FLPEP (a) or an antibody against CD11b (b) and measured by cytometry. Data was analyzed using MATLAB. The experiment was repeated three times and a representative experiment is shown. (c) Cells differentiated as in (a) and (b) were plated under agarose and analyzed in an automated chemotaxis assay by time-lapse microscopy with chemoattractant uncaging of Nv-fMLF. A cell directionality parameter measuring the angular bias of cell movement towards the gradient source is shown. This experiment was performed three times. Error bars represent the standard error of the mean. Images and statistics were processed using custom MATLAB software. Asterisk indicates a p-value of less than 0.05 using a t-test. d Cells differentiated as in (a) and (b) were mixed in suspension with pHrodo Green-labeled dead Staphylococcus aureus bioparticles for 2 h at 37 degrees. Phagocytosis of the particles was then analyzed by cytometry. e PLB-985 cells were differentiated into a neutrophil-like state as in (a) and (b). Cells were then incubated with NBT solution and 100 ng/mL of PMA or 1 μM fMLF, at 37 °C for 15 min, and measured by cytometry

Next, we wanted to determine if the increase in CD11b- and FPR1-expressing cells is correlated with a more general increase in the functional maturation of the cell population. With this in mind, we tested three major behaviors associated with neutrophilic differentiation: chemotaxis, phagocytosis and oxidative burst. We used an automated chemotaxis assay to measure the response of PLB-985 cells to a gradient of fMLF, an FPR1 ligand [25]. When cells were differentiated with DMSO supplemented with Nutridoma and 0.5% FBS, we observed a statistically significant increase (approximately 30%) in the directional accuracy of chemotaxis, validating the efficacy of this differentiation media (Fig. 2c). We measured phagocytosis using pHrodo Green labeled dead Staphylococcus aureus bioparticles. The addition of Nutridoma to the differentiation media resulted in a higher percentage of cells phagocytosing pHrodo-labeled particles (Fig. 2d, Additional file 1: c). To quantify oxidative burst, we measured the ability of PLB-985 cells to produce reactive oxygen species using the nitroblue tetrazolium (NBT) reduction assay. In response to stimulation with either PMA or fMLF, a higher fraction of cells differentiated with Nutridoma generated an oxidative burst response, as measured by reduction of NBT (Fig. 2e). Thus, in addition to increasing the surface expression of neutrophil markers, the Nutridoma differentiation protocol results in increased responsiveness for chemotaxis, phagocytosis, and oxidative burst.

To further assess the degree of differentiation, we assessed nuclear size and morphology by imaging Hoechst-stained nuclei. Although differentiation of the cell line did not recapitulate the lobed morphology of primary neutrophils, it resulted in bean shaped or partially lobed nuclei (Additional file 3). Also, as expected, the nuclear size decreased after each of the differentiation protocols. Somewhat surprisingly, this decrease was smaller using the Nutridoma differentiation protocol.

Finally, we wanted to compare the results obtained with the PLB-985 cell line to primary neutrophils as a benchmark for the validation of the differentiated cell line. We isolated highly purified neutrophils from whole human blood by immunomagnetic negative selection, and performed cell surface staining and functional assays. The primary neutrophils showed highly uniform expression of CD11b and FPR1, with CD11b expression exceeding that of the cell line, but FPR1 expression actually being lower than that of the highest expressing PLB-985 cells (Fig. 3a and b). The primary neutrophils efficiently phagocytosed pHrodo-labeled S. aureus particles, but at a level that was similar to that observed for PLB-985 cells differentiated with DMSO and Nutridoma (Fig. 3c). Finally, we assessed chemotaxis. Primary neutrophils showed markedly lower basal motility and slightly reduced stimulated speed when compared to differentiated PLB-985 cells (Fig. 3d and e). However, their chemotactic responses were strong, with a marked increase in speed after stimulation (Fig. 3d and e) and higher average directional accuracy than the cell line (Fig. 3f). Interestingly, we found that the difference in chemotactic directionality could be explained by the heterogeneous expression of FPR1 in the cell line. Since only about 70% of cells differentiated with DMSO and Nutridoma express FPR1 (Fig. 1b, Table 1), we simulated the chemotaxis of a mixed cell population in which 70% chemotax equivalently to primary neutrophils and 30% move with random directionality. This simulated data closely matched the observed directionality of the differentiated PLB-985 cells (Fig. 3g).

Fig. 3
figure 3

Comparison of differentiated PLB-985 cells with primary neutrophils. a Undifferentiated PLB-985 cells, PLB-985 cells differentiated with DMSO and Nutridoma, and primary neutrophils were stained with antibody against CD11b and analyzed by cytometry. Staining of primary neutrophils with an isotype control antibody is also shown. b The same cell samples were stained with FLPEP and analyzed by cytometry. Unstained samples are also shown for comparison. c Phagocytosis of pHrodo Green-labeled dead Staphylococcus aureus bioparticles. Negative controls in which no particles were present are also shown. d-f Differentiated PLB-985 cells and primary neutrophils were plated under agarose and analyzed in an automated chemotaxis assay by time-lapse microscopy with chemoattractant uncaging of Nv-fMLF. Mean cell speed was measured both before (d) and after (e) generation of an fMLF gradient. A cell directionality parameter measuring the angular bias of cell movement towards the gradient source is shown (f). An angular bias of 90 degrees indicates perfect directionality, and zero degrees indicates random orientation. This experiment was performed three times. Error bars represent the standard error of the mean. Asterisk indicates a p-value of less than 0.05 and double asterisks indicate a p-value of less than 0.01 using a t-test. g Histograms of the instantaneous directionality of individual cells is shown for the same experiments analyzed in (f). Here an angle of zero indicates optimal directionality, and 180 degrees indicates movement in the opposite direction. A third curve (yellow) indicates simulated data for a population of neutrophils in which only 70% of cells express the receptor FPR1

The presence of Nutridoma, rather than the absence of serum, is responsible for increased FPR1 surface expression and chemotactic efficiency

We wanted to determine whether the increased FPR1 expression was due to the reduced serum levels or a positive effect from the Nutridoma supplement. To test this, we titrated the FBS concentration in differentiation media containing fixed concentrations of DMSO and Nutridoma. We found that most of the increase in FPR1 expression was due to the presence of Nutridoma, and FPR1 expression depended only weakly on the concentration of FBS. However, 0.5% FBS gave a slightly higher population of FPR1 positive cells (Additional file 4: a), consistent with previous reports [23].

Although the molecular components of Nutridoma are proprietary, we tested the cytokine G-CSF as a candidate active component. G-CSF is a primary regulator of myeloid cell differentiation in humans, and it has been shown to promote neutrophilic differentiation in vitro in combination with ATRA [24, 31]. To investigate whether G-CSF could be the factor in Nutridoma causing the cells to upregulate FPR1 expression, we differentiated PLB-985 cells using DMSO combined with either Nutridoma or G-CSF [32]. However, based on CD11b and FPR1 expression, G-CSF was not only unable to substitute for Nutridoma, but it actually inhibited DMSO-mediated differentiation (Additional file 4: b and c). This result was surprising, given the role of G-CSF in neutrophil differentiation in vivo. However, G-CSF acts on multiple cell types in coordination with other stimuli, and it is possible that in this context it is promoting formation of an alternate or intermediate cell state. Based on the above results, we have adopted an optimized differentiation protocol with DMSO in media supplemented with 0.5% FBS and 2% Nutridoma.

Transcriptional profiles of HL-60 and PLB-985 cells

Having validated an optimal differentiation protocol using DMSO and Nutridoma based on expression of surface markers, we next wanted to compare global gene expression in myeloid cell lines to that of primary neutrophils. While another group had previously compared gene expression in HL-60 cells differentiated with either DMSO or ATRA using microarrays [33], no analysis of HL-60 or PLB-985 cell differentiation had yet been completed with the depth and precision that is now possible with RNA-seq. To this end, we extracted total RNA from the cell lines before and after differentiation (as well as at two intermediate time points for PLB-985 cells) and processed it for RNA-seq (Fig. 4a). To account for measurement noise and day-to-day variability, we repeated the RNA-Seq experiment on two independent days separated by more than 2 months. For comparison, we collated raw datasets for primary human neutrophils [34,35,36], undifferentiated HL-60 cells [37], and primary mouse neutrophils [38, 39] from previously published studies [40]. We re-analyzed these data sets along-side our own data using identical parameters and normalization (see Methods for details, and Additional file 5 and Additional file 6 for the full data set).

Fig. 4
figure 4

Transcriptional profiles of HL-60 and PLB-985 cells before and after differentiation and comparison to that of primary neutrophils. a Experimental design: HL-60 and PLB-985 cells were differentiated by culturing in media supplemented with 1.3% DMSO and 9% FBS or 1.3% DMSO, 2% Nutridoma and 0.5% FBS. Total RNA of cells harvested at the indicated time points was isolated. PolyA enrichment and RNA-seq were performed by Applied Biological Materials Inc. through their Total RNA Sequencing service. RNA-seq data was subjected to our analysis pipeline as described in Methods. b Heat map of Spearman correlations between transcriptional profiles of undifferentiated and differentiated HL-60 and PLB-985 cells and human and mouse primary neutrophils, organized by unbiased hierarchical clustering. Color bar indicates correlation strength. Each row/column represents an independent sample. For the cell line data, the labels d2, d4, and d6 indicate 2, 4, and 6 days of differentiation, respectively. For the Thomas et al. data, the labels s1, s2, s3, and s4 indicate sample numbers, which correspond to the sample numbers in Table 2

To compare human and mouse neutrophil transcriptomes, we restricted our analysis to the set of genes mapped as homologs between the two organisms in NCBI’s Homologene database [40]. We then performed unbiased hierarchical clustering of the different transcriptomes using the Spearman rank correlation coefficient as a similarity metric (Fig. 4b) and a principal component analysis (Additional file 7). As expected, in both analyses the human primary neutrophil datasets clustered together, as did the mouse primary neutrophil datasets. The cell line datasets also clustered together, with clear subclusters for undifferentiated and differentiated subsets. The latter lie closer to the human primary neutrophil datasets. As one point of caution, the datasets reanalyzed in this study used different protocols for neutrophil isolation, and the presence of some immature cells in the Taylor et al. study, which used density centrifugation of bone marrow cells, cannot be excluded. A comparison of the neutrophil isolation methods and RNA processing strategies is presented in Table 2.

Table 2 Overview of processing of cells and RNA for RNA-seq data sets analyzed in this study

Differentiation results in a neutrophil-like gene expression pattern

The differentiation of the cell lines with either DMSO or DMSO and Nutridoma resulted in large-scale changes in gene expression (Fig. 5a) and increased the correlation of their gene expression profiles with that of primary human neutrophils (Fig. 4b). The correlation strength increased steadily throughout differentiation, with the biggest change occurring in the first 2 days (Figs. 4b and 5b). To compare the global change in transcription during differentiation, we plotted density-colored scatter plots comparing primary neutrophil data for protein-coding genes to that of undifferentiated or differentiated (with DMSO and Nutridoma) PLB-985 cells (Fig. 5c). The plots revealed a global increase in correlation strength, with the bulk of points falling in a tighter more elongated cloud along the diagonal after differentiation. We saw the same effect when HL-60 cell datasets were compared to primary neutrophil data (Additional file 8).

Fig. 5
figure 5

PLB-985 cell differentiation results in neutrophil-like gene expression patterns. a Histogram of log10 fold expression changes between replicate measurements of undifferentiated PLB-985 cells (ND), and between undifferentiated cells and those differentiated with either 1.3% DMSO and 9% FBS (DMSO) or DMSO + 2% Nutridoma + 0.5% FBS (Nutri). b Spearman’s correlation coefficients for the similarity between transcriptional profiles of primary human neutrophils and PLB-985 cells at 0, 2, 4 and 6 days post differentiation with DMSO + 2% Nutridoma + 0.5% FBS. c Density-colored scatter plots of expression values (normalized FPKM on a log10 scale) for protein-coding genes for primary human neutrophils versus undifferentiated PLB-985 cells (left) or those differentiated with DMSO + Nutridoma (right)

As both differentiated human cell lines and primary mouse neutrophils serve as models for human neutrophils, we asked which gene expression pattern correlates more strongly with that of human neutrophils. While both expression patterns correlated more strongly with the human neutrophil data than the undifferentiated cell line data did, we noted that the differentiated cell line data actually correlated more strongly than the mouse neutrophil data (Fig. 4b). These results should be taken with some caution, as there may be systematic differences due to different experimental protocols or imperfect mapping of homologs. However, based on global gene expression, our study suggests that the differentiated cell lines resemble human neutrophils as much, or even more so, than mouse neutrophils do. We also noted that while we did see systematic differences between datasets collected by different groups, these differences were generally small relative to the differences between cell types or differentiation states. Indeed, every individual human neutrophil dataset correlated better with the differentiated cell lines than it did with the mouse neutrophil datasets (Fig. 4b).

Nutridoma affects gene expression during differentiation subtly, but consistently in HL-60 and PLB-985 cells

Since the presence of Nutridoma increased cell surface expression of FPR1 and CD11b during differentiation (Table 1, Fig. 2a, b), we investigated whether Nutridoma broadly increases neutrophil-like differentiation at the level of transcription. First, we compared the expression of FPR1 in our RNA-seq data to surface expression data at different time points during differentiation (Additional file 9: a). Consistent with our earlier observations, the addition of Nutridoma increased the FPR1 surface expression level markedly between day 4 and day 6. However, while Nutridoma also increased FPR1 expression at the mRNA level, this change was considerably smaller than what we observed at the cell surface protein level (Additional file 9: a). These results suggest that Nutridoma may influence the cells at the protein level to make them more functionally mature, such as an effect on translation or vesicular trafficking of receptors to the plasma membrane.

Next, we compared global gene expression patterns in cells differentiated with DMSO or with the combination of DMSO and Nutridoma. Although each protocol resulted in gene expression patterns that correlated approximately equally well with the primary human neutrophil data (Fig. 4b), we wanted to check if there were reproducible differences in gene expression between the two conditions. To do this, we compared the log10 fold-changes in expression between the two protocols and compared these changes from each of our two replicates of the RNA-Seq experiments. Even though the differences in gene expression were typically modest, they did correlate between the two replicates (Additional file 9: b, c). A further analysis of gene annotations using the DAVID functional annotation tool [41], revealed that the set of genes upregulated during differentiation in the presence of Nutridoma (relative to differentiation with DMSO only) are enriched for several annotations including the Gene Ontology (GO) Biological Process annotation of “immune response”, the GO Cellular Compartment annotation “plasma membrane”, and STAT5B transcription factor binding sites (Additional file 10).

Together our observations indicate that, although the correlation in gene expression between primary neutrophils and PLB-985 cells differentiated with or without Nutridoma is roughly equivalent, the presence of this supplement in the differentiation media subtly, but reproducibly affects gene expression at the mRNA level. In addition, Nutridoma also causes non-transcriptional changes which increase the cell surface expression of the neutrophil markers FPR1 and CD11b.

SNP analysis confirms that PLB-985 is genetically identical to HL-60

Although the PLB-985 cell line was originally described as a distinct cell line established from an acute myeloid leukemia patient in 1985 [11], Drexler and colleagues more recently reported that this cell line is actually a subclone of the HL-60 cell line [12]. The description of PLB-985 as a misidentified cell line can be found in the National Center for Biotechnology Information (NCBI) BioSample database (https://www.ncbi.nlm.nih.gov/biosample/3151776) or in the ICLAC Database of Cross-contaminated or Misidentified Cell Lines (http://iclac.org/databases/cross-contaminations/). The RNA-seq data that we generated for both cell lines allowed us to investigate this reclassification. We used the Genome Analysis Toolkit (GATK) pipeline to call variant genotypes for each sample from our RNA-Seq raw data and from the primary neutrophil data which could serve as a negative control [42]. As expected, we found identical genotypes for replicate samples, and distinct genotypes for primary neutrophil data from different donors. Consistent with previous reports, we also found identical genotypes for the PLB-985 and HL-60 cell lines (Additional file 11). However, even though they are genetically identical, our RNA-Seq data indicates that there are reproducible gene expression differences between PLB-985 and HL-60 cells.

We conclude that PLB-985 is a sub-line of HL-60 with some differences in gene expression which may make it better suited for some experimental investigations.

A searchable database of neutrophil and neutrophil-like cell line gene expression

A major goal of our study was to create an accessible database of gene expression for neutrophils and neutrophil-like cell lines. Such a tool would allow researchers to assess the expression of gene families of interest, or compare gene expression in cell lines to that in primary neutrophils to verify similarities or identify key differences. To highlight how our compiled datasets can be used, we analyzed expression levels of four sets of genes known to be important for neutrophil functions: cell surface receptors, proteins involved in reactive oxygen production, guanine nucleotide exchange factors (GEFs) for Rho family GTPases, and adhesion receptors (Fig. 6a-d). In each case, differentiation broadly increased the similarity in gene expression of these gene sets between the cell lines and primary neutrophils, including upregulation of chemoattractant receptors, oxidative burst-related genes, adhesion molecules, and other neutrophil-associated genes (Fig. 6a-d). In contrast, related oxidase components and adhesion receptors not associated with neutrophils remained lowly expressed (Fig. 6b, d). However, key differences between the cell lines and primary neutrophils are also apparent. For example, the formyl peptide receptor FPR1 is highly expressed in the cell lines, but expression of the IL-8 cytokine receptor CXCR1 is lacking (Fig. 6a). Similarly, many GEFs including PREX1, ARHGEF2 (also known as GEF-H1), VAV1, and ARHGEF6 (also known as alpha-PIX) are expressed at levels close to or equal to that found in neutrophils. However, DOCK5 and PLEKHG3 are expressed at much lower levels.

Fig. 6
figure 6

An easily searchable online database of neutrophil and neutrophil-like cell gene expression. a-d Gene expression data is shown (normalized FPKM data on a log10 scale) for the indicated genes in undifferentiated PLB-985 cells, PLB-985 cells differentiated with the DMSO only protocol, PLB-985 cells differentiated with the Nutridoma + DMSO protocol, and primary human neutrophils (collated from 3 published studies). Shown is data for select receptors (a), genes related to the production of reactive oxygen species (b), select guanine nucleotide exchange factors (GEFs) for Rho Family GTPases (c) and adhesion molecules (d). e Our collated RNA-seq data is available in an easily searchable form on our lab website: http://collinslab.ucdavis.edu/neutrophilgeneexpression/. Data can be searched by gene name, or by the presence of specific PFAM domains of interest. Data can also be sorted by expression level in any of the samples. In this example, we identify the top 10 genes containing PKINASE domain, ordered by expression level in primary human neutrophils

Finally, to make our RNA-seq results easily accessible to other researchers, we published a comprehensive, interactive database on our website (http://collinslab.ucdavis.edu/neutrophilgeneexpression/) (Fig. 6e). This database is searchable by gene id, symbol, or protein domain and can be sorted by expression level in any of the cell types. In addition to our own RNA-seq data, we included published RNA-seq data for primary human and Mus musculus (mouse) neutrophils, and links to mouse and Dictyostelium discoideum homologs. D. discoideum is an amoeboid protozoan and a valuable model organism for the study of chemotaxis and phagocytosis. This complete merged data set is also available as a table in Additional file 12.

Conclusions

We first compared established differentiation protocols for myeloid cell lines, with a focus on maximizing cell surface expression of neutrophil markers without compromising cell viability. Our study verifies that a DMSO-based differentiation protocol for HL-60 and PLB-985 cell lines gives superior differentiation and cell viability relative to other common protocols, and indicates that addition of Nutridoma may be preferable for studies of chemotaxis, phagocytosis and oxidative burst. Our RNA-seq data allowed us to generate a map of gene expression in differentiated HL-60 and PLB-985 cells which could be combined with primary human and mouse neutrophil data from the literature. We annotated this map, creating a database that is searchable by gene name and by protein domain content to allow rapid interrogation of gene expression in neutrophils and neutrophil-like cell lines. Our neutrophil gene expression database will be a valuable tool to identify similarities and differences in gene expression between the cell lines and primary neutrophils, to compare expression levels for genes of interest, and to improve the design of tools for genetic perturbations.

Methods

Reagents and antibodies

DMSO, ATRA, CA, dbcAMP, PMA, fMLF and NBT were purchased from Sigma; Nutridoma-CS from Roche, Trypan Blue solution from Corning, G-CSF from Applied Biological Materials (abm), FLPEP (catalog # F1314); Sytox Blue (catalog #S11348), TO-PRO3 Ready Flow dead cell stain (catalog #R37170), and NucRed Dead 647 (catalog #R37113) from Life Technologies; and Anti-CD11b-APC (Clone ICRF44, catalog # 301309), isotype control mouse IgG1κ (Clone MOPC-21, catalog # 400107) and Fc Receptor Blocking Solution (catalog # 422302) from Biolegend.

Cell culture and differentiation

HL-60 and PLB-985 cells were obtained as a gift from the laboratory of Dr. Orion Weiner. Generation of the stable cell lines expressing fusions of histone H2B to mTurquoise or mCherry were described elsewhere [25]. The PLB-985 cell line expressing H2B-mCherry was used for all the experiments described throughout the manuscript. Cells were cultured in complete media [RPMI 1640 with HEPES and GlutaMAX (Gibco, catalog # 72400–047) supplemented with 9% fetal bovine serum (FBS) as well as streptomycin (100 mg/mL) and penicillin (100 U/mL) (P/S)] at 37 °C in a humidified atmosphere with 5% CO2. Cell cultures were passaged two to three times per week, maintaining cell densities between 105 and 2 × 106 cells per ml. The cells were differentiated into a neutrophil-like state by culturing at an initial density of 2 × 105 using one of seven media recipes: a) RPMI 1640 complete media supplemented with 1.3% DMSO for 6 days; b) RPMI 1640 media supplemented with 5% FBS, 0.5% DMF and 2 uM ATRA (renewed on day 3 of the 6-day differentiation period); c) RPMI 1640 complete media supplemented with 750 uM dbcAMP for 6 days; d) RPMI 1640 complete media supplemented 1.3% DMSO and 2 uM ATRA for 6 days; e) RPMI 1640 complete media supplemented with 1.3% DMSO, 2 uM ATRA and 52 uM CA for 6 days; f) RPMI 1640 media supplemented with different FBS concentrations (0.5, 2, 5 or 10%), 1.3% DMSO and 2% Nutridoma for 6 days; g) RPMI 1640 media supplemented with 0.5% FBS, 1.3% DMSO and 30 ng/ml G-CSF for 6 days.

Human primary neutrophil isolation

Ethical approval for the study of neutrophils from adult healthy controls was granted by the Institutional Review Board (IRB) from the University of California, Davis. All participants gave written, informed consent. Blood was collected by finger prick using a microlet lancing device to isolate small drops of blood. Neutrophils were isolated using negative selection with the EasySep Direct Human Neutrophil isolation kit (StemCell), following the manufacturer’s instructions.

Nuclear area calculation

Fluorescence images were acquired using a Nikon Ti-E inverted microscope with a 20× (0.75 NA) Plan Apochromat objective with an optovar for an additional 1.5× magnification. DNA was detected using mCherry tagged histone H2B for cell lines and Hoechst staining for primary human neutrophils. All cell types were stained with Sytox Blue to determine viability. Viability scoring, nuclear segmentation and calculation of nuclear areas were carried out in MATLAB using custom software. Briefly, nuclei were segmented using an intensity threshold. A minimum area threshold of 100 pixels (~44uM) was implemented to filter out debris. Dead cells (nuclei positive for Sytox stain) were also excluded.

Immunolabeling and cytometry

Differentiation of PLB-985 cells into a neutrophil-like state was assessed by measuring the levels of CD11b and the formyl peptide receptor FPR1 on the cell surface. For detection of CD11b, cells were harvested, washed with phosphate-buffered saline (PBS) and stained with primary antibody antiCD11b-APC (clone ICRF44) along with Fc Receptor Blocking Solution (for 45 min at 4 °C in FACS buffer (0.5% bovine serum albumin and 0.05% sodium azide in PBS). After fluorescent labeling, the samples were washed and acquired with the BD FACSCanto II flow cytometer (BD Biosciences, Franklin Lakes, NJ). Relative levels of FPR1 were assessed by binding of N-formyl-norleucyl-leucyl-phenylalanyl-norleucyl-tyrosyl-lysine-fluorescein (FLPEP) (Life Technologies product number F1314), a fluorescent ligand of the receptor. Cells were mixed with ice cold media containing FLPEP to give a final concentration of 10 nM, placed on ice for 10 min, and acquired with the BD FACSCanto II flow cytometer. For one condition (DMSO+ATRA+CA), cellular autofluorescence was high in the fluorescein channel, and it was corrected by subtracting the mean autofluorescence from the unstained control samples. We measured cell death using NucRed Dead 647 Probe (Life Technologies, product number R37113) or TO-PRO3 (R37170) according to the manufacturer’s protocol. FACS data was analyzed using MATLAB (Mathworks, Inc.). To determine the percentage of cells expressing the target, we set thresholds using an appropriate negative control (isotype immunoglobulins) for CD11b or non-stained cells for FLPEP and NucRed Dead 647.

Cell growth

To measure cell proliferation, cells were counted at day three and six of the differentiation protocols, using the trypan blue dye exclusion test. The number of cells were normalized against the initial number of cells.

Chemotaxis assay

The chemotaxis assay was performed as previously described [25]. Briefly, differentiated PLB-985 cells were plated in 96-well plates under low melting temperature agarose and imaged by time-lapse microscopy with chemoattractant uncaging. A caged UV sensitive derivative N-nitroveratryl derivative (Nv-fMLF) of N-formyl-methioninealanine-phenylalanine was used, and UV photo-uncaging was performed to generate chemical gradients. Image processing and statistical analyses of chemotaxis were performed using custom MATLAB software.

Phagocytosis assay

Cells were mixed in suspension with pHrodo Green-labeled dead Staphylococcus aureus bioparticles for 2 h at 37 degrees. The pHrodo Green dye markedly increases fluorescence upon reaching an acidic environment such as that of the phagosome, which greatly limits the contribution of non-internalized particles to the fluorescence signal. Phagocytosis of the particles was analyzed by cytometry with the BD FACSCanto II flow cytometer.

NBT reduction assay

Differentiated cells were harvested and resuspended in Leibovitz’s L-15 media (ThermoFisher, product number 11415–064) supplemented with FBS. Cells were then mixed with ice cold media containing 2X NBT solution (Sigma, product number N5514) and 200 ng/ml of phorbol 12-myristate 13- acetate (PMA) or 2 μM of fMLF in L-15 media and incubated at 37 °C for 15 min. After incubation, deposition of purple formazan granules (from the reduction of NBT) was acquired with the BD FACSCanto II flow cytometer.

RNA-seq analysis pipeline

Total RNA was isolated using the RNeasy Plus Mini Kit (Qiagen) according to the manufacturer’s protocol. PolyA enrichment and RNA-seq were performed by Applied Biological Materials Inc. through their Total RNA Sequencing service. The samples were subjected to polyA enrichment followed by fragmentation, first and second strand synthesis, adenylation of 3′ ends, adapter ligation, DNA fragment enrichment, and real-time PCR quantification. For each sample, over thirty million 1 × 75-base single-end reads were produced on the Illumina NextSeq 500 platform. Human raw RNA-seq reads were aligned to the GRCh38.p5 version of the human genome using the GENCODE GRCh38 (release 24) genome annotations and the STAR two-pass algorithm (STAR version 2.4) [43]. Mouse RNA-seq reads were aligned to the GRCm38 version of the mouse genome using the GENCODE M8 annotations. Gene expression was quantified using the open-source software Cufflinks (version 2.2.1) [44]. Gene expression was normalized using the “uniform genes” strategy of Glusman et al. [45]. Briefly, a set of “uniform genes” was defined as all genes with expression levels between the 50th and 90th percentile in all samples. A scaling factor was computed for each dataset as the geometric mean of the expression levels of these “uniform genes.” Each dataset was then normalized by dividing by its own scaling factor and multiplying by the mean scaling factor of the human datasets. Normalization, comparison and correlation of the different transcriptomes, and hierarchical clustering of the data, was done using MATLAB. Log base 10 expression values were calculated after normalization, and genes with zero read counts were assigned a value of − 4 to avoid infinite values. Hierarchical clustering was done using one minus the Spearman rank correlation coefficient as the distance metric and average linkage. Pairwise Spearman rank correlation coefficients were displayed in the cluster heatmap. Principal component analysis was performed using MATLAB’s “pca” function with the default options. To avoid a large influence from noise in lowly expressed genes, all log10 expression below − 2 were set to − 2 for the principal component analysis.

For our final data tables, we averaged the data over replicate samples using the geometric mean. For the primary human and mouse neutrophil datasets, we first averaged the replicate measurement performed by each laboratory group (note that the Thomas et al. [36] and Wright et al. [35] datasets are from the same laboratory), and then averaged over the different laboratories to give equal weight to each lab to account for potential lab-to-lab systematic differences.

For the density-colored scatter plots in Fig. 5a, the averaged gene expression profiles were used. Only protein-coding genes were included in the plot, based on NCBI’s “Gene type” designations (18,976 genes) [40]. Data for noncoding RNAs were much more variable, perhaps in part due to the polyA enrichment step of our experimental processing.

For the cluster of gene expression profiles in Fig. 4b and the principal component analysis in Additional file 7, only genes mapped between human and mouse as homologs in NCBI’s Homologene database [40] were included.

Variant calling for genotyping was performed using the GATK pipeline and the GATK best practices guidelines [42]. After variant calling, the genotype data was filtered to keep only loci with a read depth of at least 10 in all samples and a quality score of greater than 100. Genotypes were then compared to determine the fraction of loci with identical genotypes for each pair of samples.

Generation of the database

Normalized RNA-seq data (prepared as described above) was compiled into a database using MySQL and published on our website (collinslab.ucdavis.edu). Dictyostelium discoidum homologs were downloaded from InParanoid (http://inparanoid.sbc.su.se/cgi-bin/index.cgi).

Accession numbers

HL-60 and PLB-985 undifferentiated and differentiated RNA-seq data were deposited at the NCBI Gene Expression Omnibus (GEO) database (Accession: GSE103706). Previously published data for primary human neutrophils from the Wright et al. [GEO: GSE40548], Jiang et al. [GEO: GSE66895], and Thomas et al. [GEO: GSE70068] studies, for undifferentiated HL-60 cells from the Cusanovich et al. study [GEO: GSE68103], and for primary mouse neutrophils from the Taylor et al. [GEO: GSE55090] and Coffelt et al. [GSE55633] studies were downloaded from NCBI’s GEO database.

Abbreviations

ATRA:

All-trans retinoic acid

CA:

Caffeic acid

CD11b:

Integrin Alpha M, ITGAM

dbcAMP:

Dibutyryl cyclic adenosine monophosphate

DMF:

N,N-Dimethyl formamide

DMSO:

Dimethyl sulfoxide

FACS:

Fluorescence-activated cell sorting

FBS:

Fetal bovine serum

FLPEP:

N-formyl-norleucyl-leucyl-phenylalanyl-norleucyl-tyrosyl-lysine-fluorescein

fMLF:

N-formyl-methionyl-leucyl phenylalanine

FPR1:

Formyl Peptide Receptor 1

GATK:

Genome Analysis Toolkit

G-CSF:

Granulocyte colony-stimulating factor

GEF:

Guanine nucleotide exchange factor

MFI:

Mean fluorescence intensity

NBT:

Nitro Blue tetrazolium

NCBI:

National Center for Biotechnology Information

P/S:

Streptomycin and penicillin

PBS:

Phosphate-buffered saline

PMA:

Phorbol 12-myristate 13-acetate

RNA-seq:

RNA sequencing

SNPs:

Single-nucleotide polymorphisms

References

  1. Kolaczkowska E, Kubes P. Neutrophil recruitment and function in health and inflammation. Nat. Rev. Immunol. 2013;13:159–75.

    Article  PubMed  CAS  Google Scholar 

  2. Rappel WJ, Loomis WF. Eukaryotic chemotaxis. Wiley Interdiscip Rev Syst Biol Med. 2009;1:141–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Mantovani A, Cassatella MA, Costantini C, Jaillon S. Neutrophils in the activation and regulation of innate and adaptive immunity. Nat Rev Immunol. 2011;11:519–31.

    Article  PubMed  CAS  Google Scholar 

  4. Collins SJ, Gallo RC, Gallagher RE. Continuous growth and differentiation of human myeloid leukaemic cells in suspension culture. Nature. 1977;270:347–9.

    Article  PubMed  CAS  Google Scholar 

  5. Gupta D, Shah HP, Malu K, Berliner N, Gaines P, Gupta D, et al. Differentiation and Characterization of Myeloid Cells. Curr. Protoc. Immunol. Hoboken, Wiley; 2014. p. 22F.5.1-22F.5.28

  6. Collins SJ, Ruscetti FW, Gallagher RE, Gallo RC. Terminal differentiation of human promyelocytic leukemia cells induced by dimethyl sulfoxide and other polar compounds. Proc Natl Acad Sci U S A. 1978;75:2458–62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Hauert AB, Martinelli S, Marone C, Niggli V. Differentiated HL-60 cells are a valid model system for the analysis of human neutrophil migration and chemotaxis. Int J Biochem Cell Biol. 2002;34:838–54.

    Article  PubMed  CAS  Google Scholar 

  8. Fleck RA, Romero-Steiner S, Nahm MH. Use of HL-60 cell line to measure opsonic capacity of pneumococcal antibodies. Clin Vaccine Immunol. 2005;12:19–27.

    Article  CAS  Google Scholar 

  9. Skoge M, Wong E, Hamza B, Bae A, Martel J, Kataria R, et al. A worldwide competition to compare the speed and chemotactic accuracy of neutrophil-like cells. PLoS One. 2016;11:1–19.

    Article  CAS  Google Scholar 

  10. Birnie GD. The HL60 cell line: a model system for studying human myeloid cell differentiation. Br J Cancer Suppl. 1988;9:41–5.

    PubMed  PubMed Central  CAS  Google Scholar 

  11. Tucker KA, Lilly MB, Heck L, Rado TA. Characterization of a new human diploid myeloid leukemia cell line (PLB-985) with granulocytic and monocytic differentiating capacity. Blood. 1987;70:372–8.

    PubMed  CAS  Google Scholar 

  12. Drexler HG, Dirks WG, Matsuo Y, MacLeod RAF. False leukemia–lymphoma cell lines: an update on over 500 cell lines. Leukemia. 2003;17:416–26.

    Article  PubMed  CAS  Google Scholar 

  13. Ashkenazi A, Marks RS. Luminol-dependent chemiluminescence of human phagocyte cell lines: comparison between DMSO differentiated PLB 985 and HL 60 cells. Luminescence. 2009;24:171–7.

    Article  PubMed  CAS  Google Scholar 

  14. Yang HW, Collins SR, Meyer T. Locally excitable Cdc42 signals steer cells during chemotaxis. Nat Cell Biol. 2015;18:191–201.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Servant G, Weiner OD, Neptune ER, Sedat JW, Bourne HR. Dynamics of a chemoattractant receptor in living neutrophils during chemotaxis. Mol Biol Cell. 1999;10:1163–78.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Majumdar R, Tavakoli Tameh A, Parent CA. Exosomes mediate LTB4 release during neutrophil chemotaxis. PLoS Biol. 2016;14:1–28.

    Article  CAS  Google Scholar 

  17. Yaseen R, Blodkamp S, Lüthje P, Reuner F, Völlger L, Naim HY, et al. Antimicrobial activity of HL-60 cells compared to primary blood-derived neutrophils against Staphylococcus aureus. J. Negat. Results Biomed. 2017;16:2.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Collins SJ. The HL-60 promyelocytic leukemia cell line: proliferation, differentiation, and cellular oncogene expression. Blood. 1987;70:1233–44.

    PubMed  CAS  Google Scholar 

  19. Friend C, Scher W, Holland JG, T Sato. Hemoglobin synthesis in murine virus-induced leukemic cells in vitro: stimulation of erythroid differentiation by dimethyl sulfoxide. Proc Natl Acad Sci. 1971;68:378–82.

  20. Siebenlist U, Bressler P, Kelly K. Two distinct mechanisms of transcriptional control operate on c-myc during differentiation of HL60 cells. Mol Cell Biol. 1988;8:867–74.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Jorgenson KF, Antoun GR, Zipf TF. Chromatin structural analysis of the 5′ end and contiguous flanking region of the myeloperoxidase gene. Blood. 1991;77:159–64.

    PubMed  CAS  Google Scholar 

  22. Veselska R, Zitterbart K, Auer J, Neradil J. Differentiation of HL-60 myeloid leukemia cells induced by all-trans retinoic acid is enhanced in combination with caffeic acid. Int J Mol Med. 2004;14:305–10.

  23. Pedruzzi E, Fay M, Elbim C, Gaudry M, Gougerot-Pocidalo MA. Differentiation of PLB-985 myeloid cells into mature neutrophils, shown by degranulation of terminally differentiated compartments in response to N-formyl peptide and priming of superoxide anion production by granulocyte-macrophage colony-stimulating fact. Br J Haematol. 2002;117:719–26.

    Article  PubMed  CAS  Google Scholar 

  24. Nakajima H, Ihle JN. Granulocyte colony-stimulating factor regulates myeloid differentiation through CCAAT/enhancer-binding protein epsilon. Blood. 2001;98:897–905.

    Article  PubMed  CAS  Google Scholar 

  25. Collins SR, Yang HW, Bonger KM, Guignet EG, Wandless TJ, Meyer T. Using light to shape chemical gradients for parallel and automated analysis of chemotaxis. Mol Syst Biol. 2015;11:804.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Marin-Esteban V, Turbica I, Dufour G, Semiramoth N, Gleizes A, Gorges R, et al. Afa/Dr diffusely adhering Escherichia coli strain C1845 induces neutrophil extracellular traps that kill bacteria and damage human enterocyte-like cells. Infect Immun. 2012;80:1891–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Chang HH, Oh PY, Ingber DE, Huang S. Multistable and multistep dynamics in neutrophil differentiation. BMC Cell Biol. 2006;7:11.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Le Cabec V, Calafat J, Borregaard N. Sorting of the specific granule protein, NGAL, during granulocytic maturation of HL-60 cells. Blood. 1997;89:2113–21.

    PubMed  CAS  Google Scholar 

  29. Meyer PA, Kleinschnitz C. Retinoic Acid Induced Differentiation and Commitment in HL-60 Cells. Environ Health Perspect. 1990;88:179–82.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Yu L, Zhen L, Dinauer MC. Biosynthesis of the phagocyte NADPH oxidase cytochrome b 558. J Biol Chem. 1997;272:27288–94.

    Article  PubMed  CAS  Google Scholar 

  31. Bunce CM, French PJ, Durham J, Stockley RA, Michell RH, Brown G. Indomethacin potentiates the induction of HL60 differentiation to neutrophils, by retinoic acid and granulocyte colony-stimulating factor, and to monocytes, by vitamin D3. Leukemia. 1994;8:595–604.

    PubMed  CAS  Google Scholar 

  32. Fleck RA, Athwal H, Bygraves JA, Hockley DJ, Feavers IM, Stacey GN. Optimization of nb-4 and hl-60 differentiation for use in opsonophagocytosis assays. In Vitro Cell Dev Biol Anim. 2003;39:235–42.

    Article  PubMed  Google Scholar 

  33. Huang S, Eichler G, Bar-Yam Y, Ingber DE. Cell fates as high-dimensional attractor states of a complex gene regulatory network. Phys Rev Lett. 2005;94:1–4.

    Google Scholar 

  34. Jiang K, Zhu L, Buck MJ, Chen Y, Carrier B, Liu T, et al. Disease-associated single-nucleotide polymorphisms from noncoding regions in juvenile idiopathic arthritis are located within or adjacent to functional genomic elements of human neutrophils and CD4+ T cells. Arthritis Rheumatol. 2015;67:1966–77.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Wright HL, Thomas HB, Moots RJ, Edwards SW. RNA-Seq reveals activation of both common and cytokine-specific pathways following neutrophil priming. PLoS One. 2013;8:e58598.

  36. Thomas HB, Moots RJ, Edwards SW, Wright HL. Whose gene is it anyway? The effect of preparation purity on neutrophil transcriptome studies. PLoS One. 2015;10:1–15.

    Google Scholar 

  37. Cusanovich DA, Daza R, Adey A, Pliner HA, Christiansen L, Gunderson KL, et al. Multiplex single-cell profiling of chromatin accessibility by combinatorial cellular indexing. Science. 2015;348:910–4.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Taylor A, Tang W, Bruscia EM, Zhang PX, Lin A, Gaines P, et al. SRF is required for neutrophil migration in response to inflammation. Blood. 2014;123:3027–36.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Coffelt SB, Kersten K, Doornebal CW, Weiden J, Vrijland K, Hau C-S, et al. IL-17-producing γδ T cells and neutrophils conspire to promote breast cancer metastasis. Nature. 2015;522:345–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. NCBI Resource Coordinators. Database resources of the National Center for biotechnology information. Nucleic Acids Res. 2017;45:D12–7.

  41. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4:44–57.

    Article  CAS  Google Scholar 

  42. Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, del Angel G, Levy-Moonshine A, et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinforma. Hoboken, NJ, USA: John Wiley & Sons, Inc.; 2013. p. 11.10.1–11.10.33.

  43. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    Article  PubMed  CAS  Google Scholar 

  44. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

  45. Glusman G, Caballero J, Robinson M, Kutlu B, Hood L. Optimal Scaling of Digital Transcriptomes. PLoS One. 2013;8:e77885.

Download references

Acknowledgements

We are grateful to colleagues who generously provided reagents. We thank George Bell, Dean Natwick and Michael Pargett for stimulating discussion and critical assessment of the manuscript.

Funding

This work was funded by a Kimmel Scholar Award, an NIH Director’s New Innovator Award (DP2 HD094656), a UC Davis Hellman Fellowship, and a grant from the UC Davis Academic Senate Committee on Research to SRC. BLRG was funded by a Floyd and Mary Schwall Fellowship and an NIH IMSD Fellowship (R25 GM056765). These funding bodies had no role in the design of the study, collection, analysis, and interpretation of data, or in writing the manuscript.

Availability of data and materials

The datasets generated and analyzed during the current study are available in the NCBI GEO database under the accession number GSE103706.

Author information

Authors and Affiliations

Authors

Contributions

ER participated in the experimental work, data compilation, analysis and writing of the manuscript. BLRG participated in data analysis, generation of the database, and writing of the manuscript. SRC participated in the data compilation, analysis and writing of the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Sean R. Collins.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for the study of neutrophils from healthy adult volunteers was granted by the Institutional Review Board (IRB) from the University of California, Davis (IORG0000251). All participants gave written, informed consent.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Comparison of cell surface marker expression with different differentiation protocols. PLB-985 cells were differentiated into a neutrophil-like state by culturing in media supplemented with 1.3% DMSO and 2 μM ATRA for 4 days, 1.3% DMSO and 2 μM ATRA for 6 days, 0.5% DMF for 6 days, 1 μM ATRA for 6 days, 2 μM ATRA for 6 days, or 0.5% DMF and 2 μM ATRA for 6 days. Undifferentiated cells were also analyzed. Cells were stained with an antibody against CD11b, chosen as an early differentiation marker (a), or with FLPEP (a fluorescent ligand of FPR1), as a late differentiation marker (b). (c) PLB-985 cells differentiated for 6 days with either 1.3% DMSO, 1.3% DMSO and 2% Nutridoma-CS, 0.5% DMF and 2 μM ATRA, or dbcAMP were mixed in suspension with pHrodo Green-labeled dead Staphylococcus aureus bioparticles for 2 hours at 37 degrees. Phagocytosis of the particles was then analyzed by cytometry and data was analyzed using MATLAB. (PDF 162 kb)

Additional file 2:

The effect of Nutridoma on surface marker expression is also seen in HL-60 cells. HL-60 cells were differentiated into a neutrophil-like state by culturing in media supplemented with 1.3% DMSO and 9% FBS or supplemented with 1.3% DMSO, 2% Nutridoma and 0.5% FBS, for 6 days. Cells were stained with an antibody against CD11b (a) or the fluorescent FPR1 ligand FLPEP and measured by cytometry. Undifferentiated cells were also analyzed. Data was analyzed using MATLAB. (PDF 134 kb)

Additional file 3:

Differentiation of PLB985 cells reduces nuclear area and alters nuclear morphology. (a) Nuclear areas were measured by fluorescence imaging of either histone H2B-mCherry for PLB-985 cells or Hoechst-labeled DNA for primary neutrophils. A histogram of nuclear areas for undifferentiated PLB-985 cells (blue; n = 6203), PLB-985 cells differentiated with DMSO alone (red; n = 4338), PLB-985 cells differentiated with DMSO + Nutridoma CS (yellow; n = 4393), PLB-985 cells differentiated with DMF + ATRA (purple; n = 2386) and human primary neutrophils (green; n = 2264). The X axis represents nuclear area in square pixels, where 1 pixel = 0.4389 μm. Staining with Sytox Blue was used to exclude dead cells from the analysis. (b-f) Representative images of nuclei where scale bar equals 20 μm. DNA is shown in red and Sytox Blue dye (a membrane impermeable DNA-binding dye) is shown in green. (PDF 1850 kb)

Additional file 4:

(a) PLB-985 cells were differentiated into a neutrophil-like state by culturing in media supplemented with 1.3% DMSO, 2% Nutridoma and the indicated percentage of FBS. Cells were stained with FLPEP and measured by cytometry. (b-c) PLB-985 cells were differentiated into a neutrophil-like state by culturing in media supplemented with 1.3% DMSO 0.5% FBS and either 2% Nutridoma or 30 ng/mL G-CSF for 6 days. Cells were stained with anti-CD11b (b) or with FLPEP (c) and measured by cytometry. Cytometry data was analyzed using MATLAB in all the experiments. Data for the undifferentiated cells and those differentiated with DMSO or DMSO + Nutridoma is the same as that shown in Fig. 2. (PDF 158 kb)

Additional file 5:

Full tables of processed gene expression values measured by RNA-Seq in this study. This Excel file contains two sheets. The first sheet contains FPKM gene expression values generated by Cufflinks for all samples analyzed by RNA-Seq in this study. The second sheet contains the log10-transformed normalized expression values. These values were computed from the FPKM values as described in the Methods section. (XLSX 15648 kb)

Additional file 6:

Full tables of reanalyzed gene expression data for primary neutrophils and HL-60 cells from previously published studies. This Excel file contains four sheets. The first sheet contains FPKM gene expression values generated by Cufflinks for all primary human neutrophil and HL-60 samples reanalyzed in this study. The second sheet contains the corresponding log10-transformed normalized expression values. The third sheet contains FPKM gene expression values for all primary mouse neutrophil samples reanalyzed in this study, and the fourth sheet contains the corresponding log10-transformed normalized values. (XLSX 21707 kb)

Additional file 7:

Principal component analysis of the gene expression data analyzed in this study. Principal component analysis was carried out using MATLAB’s “pca” function with the default options for all data sets included in this study (the same data sets that are included in Fig. 4b). Only log10 expression values for genes with both human and mouse homologs were included for the analysis. To avoid a large influence from noise in lowly expressed genes, all log10 expression below − 2 were set to − 2 for this analysis. Shown is a scatter plot of the first two principal components. Together, these two components explained 66% of the variance in the data set. (PDF 357 kb)

Additional file 8:

Similarity of gene expression between HL-60 cells and primary human neutrophils. Density-colored scatter plots of expression values (normalized FPKM on a log10 scale) for protein-coding genes for primary human neutrophils versus undifferentiated HL-60 cells (left) or those differentiated with DMSO + Nutridoma (right). (PDF 364 kb)

Additional file 9:

Nutridoma affects gene expression during differentiation subtly, but consistently in HL-60 and PLB-985 cells. (a) PLB-985 cells were differentiated into a neutrophil-like state by culturing in media supplemented with 1.3% DMSO, 2% Nutridoma and 0.5% FBS. At the indicated days post-differentiation, cells were harvested and processed for RNA-seq analysis or stained with FLPEP. Shown are the expression level of FPR1 at the mRNA level (from our RNA-seq data) and at the surface of the cells (FLPEP staining). Values for undifferentiated cells and for cell differentiated with 1.3% DMSO and 9% FBS are shown for comparison. Values were normalized by the values for the sample differentiated with only DMSO. (b) Density-colored scatter plots of expression values (normalized FPKM on a log10 scale) for protein-coding genes for HL-60 (left) or PLB-985 cells (right) cells differentiated with DMSO + Nutridoma versus those differentiated with DMSO. These data represent the averages of two replicate experiments. (c) Density-colored scatter plots of log10 fold-changes between the two protocols (subtracting DMSO values from DMSO + Nutridoma) comparing each of our two replicates of the RNA-Seq experiments for HL-60 (left) or PLB-985 (right) cells. To avoid noise from lowly expressed genes, only protein-coding genes with a mean log10 normalized expression value of − 1 or higher were included. The Pearson’s (R) and Spearman’s correlation coefficients are shown. (PDF 417 kb)

Additional file 10:

Summary of enriched gene annotations for the set of genes upregulated in the presence of Nutridoma. This Excel file contains two sheets. The first contains a set of genes determined to be upregulated in differentiation with Nutridoma relative to differentiation with DMSO alone. This gene set was selected based on a mean log10 normalized expression value of at least − 1 across our differentiation samples, and a minimum 3-fold upregulation in each Nutridoma differentiation sample compared to its corresponding DMSO differentiation sample (for both PLB-985 and HL-60 cells). The table contains NCBI GeneIDs, Ensembl GeneIDs, gene names, and the mean log10 normalized expression values in Nutridoma and DMSO differentiation samples. Sheet 2 contains a summary of annotations enriched in for this set of upregulated genes. Annotation enrichment was calculated using the DAVID online gene annotation analysis tool [41]. (XLSX 34 kb)

Additional file 11:

SNP analysis confirms that PLB-985 is genetically identical to HL-60. Variant calling for genotyping was performed using the GATK pipeline and the GATK best practices guidelines for each sample from our RNA-Seq raw data and from the primary human neutrophil RNA-Seq data. The heatmap indicates the fraction of loci with identical genotypes for each pair of samples (see Methods for more details). (PDF 377 kb)

Additional file 12:

Table of cell line and primary neutrophil gene expression data. This comma separated value file contains the final averaged log10 normalized gene expression values for undifferentiated and differentiated cell lines, as well as primary human and mouse neutrophils. This file contains all of the data included in our online searchable neutrophil gene expression database. (CSV 2970 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rincón, E., Rocha-Gregg, B.L. & Collins, S.R. A map of gene expression in neutrophil-like cell lines. BMC Genomics 19, 573 (2018). https://doi.org/10.1186/s12864-018-4957-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-4957-6

Keywords