In-silico analysis of myeloid cells across the animal kingdom reveals neutrophil evolution by colony-stimulating factors
Abstract
Neutrophils constitute the largest population of phagocytic granulocytes in the blood of mammals. The development and function of neutrophils and monocytes is primarily governed by the granulocyte colony-stimulating factor receptor family (CSF3R/CSF3) and macrophage colony-stimulating factor receptor family (CSF1R/IL34/CSF1) respectively. Using various techniques this study considered how the emergence of receptor:ligand pairings shaped the distribution of blood myeloid cell populations. Comparative gene analysis supported the ancestral pairings of CSF1R/IL34 and CSF3R/CSF3, and the emergence of CSF1 later in lineages after the advent of Jawed/Jawless fish. Further analysis suggested that the emergence of CSF3 lead to reorganisation of granulocyte distribution between amphibian and early reptiles. However, the advent of endothermy likely contributed to the dominance of the neutrophil/heterophil in modern-day mammals and birds. In summary, we show that the emergence of CSF3R/CSF3 was a key factor in the subsequent evolution of the modern-day mammalian neutrophil.
Introduction
Phagocytes are key effector immune cells responsible for various biological processes; from orchestrating responses against invading pathogens to maintaining tissue homeostasis and neutrophils are the most abundant population of granulocytic phagocytes present in mammalian blood (Adrover et al., 2019; Hidalgo et al., 2019; Ng et al., 2019; Yvan-Charvet and Ng, 2019; Evrard et al., 2018). Neutrophil and heterophils (a functionally analogous granulocytic phagocyte population present in non-mammals and some mammals; Montali, 1988) arise from a shared pool of haematopoietic stem cells and mitotic myeloid progenitor cells that can also differentiate into monocytes, eosinophils and basophils following exposure to the relevant growth factor (Adrover et al., 2019; Yvan-Charvet and Ng, 2019; Mehta et al., 2014). The development and life cycle of mammalian neutrophils through a continuum of multipotent progenitors to a post-mitotic mature cell has been well described and has been recently reviewed (Hidalgo et al., 2019; Yvan-Charvet and Ng, 2019).
The development and function of myeloid phagocytes is mediated through lineage-specific transcription factors and pleiotropic glycoproteins - termed colony-stimulating factors (CSFs)- acting in concert on myeloid progenitor cells. CSF1, CSF3, and their cognate receptors, are lineage-specific and responsible for the differentiation and function of monocytes/macrophages and neutrophils, respectively (for full nomenclature see Supplementary file 1; Dai et al., 2002; Yoshida et al., 1990; Umeda et al., 1996; Liu et al., 1996; Lieschke et al., 1994). There is a large body of evidence demonstrating the requirement of CSFs for cell development as multiple studies in knockout mice have shown that CSF1R/CSF1 and CSF3R/CSF3 are linked to the development of monocytes and neutrophils in vivo. The loss of CSF3 or CSF3R directly affects neutrophil populations resulting in a severe neutropenia, but not the complete loss of mature neutrophils in the models studied (Liu et al., 1996; Lieschke et al., 1994; Basu et al., 2000). The loss of CSF1 caused reduced cavity development of the bone marrow, loss of some progenitor populations, monocytopenia, and a reduced population of neutrophils in the bone marrow, although interestingly, elevated levels of neutrophils were observed in the periphery (Dai et al., 2002; Yoshida et al., 1990; Hibbs et al., 2007).
Similarly, in humans, single gene mutations have been described in both CSF3 and CSF3R resulting in severe congenital neutropenia (SCN). Furthermore, mutations in ELANE, a gene that encodes for the neutrophil granule serine protease -neutrophil elastase- and acts as a negative regulator of CSF3R signalling can also result in SCN (Piper et al., 2010; Hunter et al., 2003; Garg, 2020; Nayak et al., 2015). In contrast to some animal models, individuals present as children with early onset life-threatening infections because they lack mature neutrophils in the circulation as the neutrophil progenitors in the bone marrow do not progress beyond the myelocyte/promyelocyte stage (Lakshman and Finn, 2001). These studies demonstrate that CSF receptor/ligand pairings are essential for homeostatic neutrophil development and are intrinsically linked with neutrophil function, arguably making them ideal surrogates to study neutrophil evolution. Through multiple methods we examined the emergence of the respective CSF ligand and receptor genes and proteins across the Chordate phylum and demonstrated how CSF1R/CSF1 and CSF3R/CSF3 pairings contributed to the evolutionary adaptations of the mammalian neutrophil.
Results
The neutrophilic/heterophilic granulocyte is the predominant granulocyte in the blood of mammals and aves
The presence of analogous myeloid granulocytes and agranulocytes was examined by comparing available complete blood count (CBC) data, where applicable, for various animal orders and demonstrated the possible distribution of myeloid cells in blood across evolution. The earliest chordates were represented by two groups, the first; jawed fish (4.10 × 108 million years ago [MYA]), which included; coelacanths, elasmobranchs (Whale shark) and chimaeras (Australian Ghost sharks). Jawless fish containing lampreys and hagfish (3.6 × 108 MYA) represented the second group. The next group represented chronologically was the amphibians (3.5 × 108 MYA); containing anurans and gymnophiona. The reptilian class was represented by three different orders, squamata (3.3 × 108 MYA), crocodilia (2.4 × 108 MYA) and testudines (2.1 × 108 MYA). The following group was the closely related avian class (7.0 × 107 MYA), which, contained birds from both paleognaths and neognaths. The final class of chordates evaluated was Mammalia, where all three orders; monotremata (1.1 × 108 MYA), marsupial (6.5 × 107 MYA) and placental (6.25 × 107 MYA) were represented. All the species examined for cell distribution and for the subsequent gene sequence and/or protein homology studies were visualised in a species tree (Figure 1a).

Population comparison of blood myeloid cell subset distribution in chordates demonstrates predominance of the neutrophilic granulocyte in birds and mammals.
(A) Species tree of the animals within Phylum Chordata examined, sub-classified by animal order or class. (B) Meta-analysis of aggregated phylum Chordata complete blood cell counts (excluding all lymphoid cells). Data is visualised as a floating bar and the line represents the mean value and shows absolute counts per cell type per order (i) and composition of myeloid cells per cell type (ii). (N) = Neutrophils, (H) = Heterophils, (A) = Azurophils, (M) = Monocytes.
All classes within the phylum Chordata were evaluated for populations of neutrophils/heterophils; eosinophils and basophils. Azurophils, a specialised population of granulocytes – analogous to both the mammalian neutrophil and monocyte- but unique to reptiles were also included. A meta-analysis was performed on representative aggregated data to show absolute counts per cell type per order and finally, the proportional composition of myeloid cells only per cell type (Figure 1b).
The presence of all myeloid cells was confirmed in at least one representative species from each order examined. Jawed fish, which comprised sharks and rays, were unique in having both heterophils (2.0 × 109/L) and neutrophils (0.045 × 109/L) in the shark. All other orders had either neutrophilic or heterophilic populations only. Heterophils formed a sizeable population of granulocyte cells in Amphibians (the order of Anura; 2.9 × 109/L), Reptiles (the orders; squamata, 6.9 x x109/L; testudines, 2.1 × 109/L; and crocodiles, 3.2 × 109/L) and Aves (12.7 × 109/L). In contrast, neutrophils were the majority population in all three mammalian orders (monotremata, 6.0 × 109/L; marsupalia, 5.3 × 109/L; and placental, 5.2 × 109/L [Figure 1bi]). Eosinophils and basophils were over-represented in the orders of amphibians, reptiles and aves, when compared to mammals, where they constituted minority populations [Figure 1bi].
Analysis of the proportional composition of total myeloid white blood cells (WBC) showed interesting changes in the distribution of granulocytes across the different species classes. Jawed fish lineages, represented by sharks and rays, retained both a majority heterophil population (49.0%) and a minority neutrophil population (1.6%). The cold-blooded orders; anura (29.7%), squamata (55.3%), testudines (41.9%) and crocodilia (43.7%); retained a large population of neutrophils/heterophils; although proportions of the population varied within related animal orders and between the different species classes (Figure 1bii). In contrast, within the warm-blooded classes and orders; aves (84.7%), monotremata (83.8%), marsupial (82.6%) and placental (81.8%) the opposite distribution was observed and the neutrophil/heterophil was the most abundant peripheral blood cell type (Figure 1bii). Interestingly, basophils and eosinophils were well represented in the anuran order (31.5%, and 23.2%, respectively). A similar pattern was observed in the reptilian lineage; however, there were differences between orders as squamata (10.1%) and crocodilia (13.8%) had lower proportions of eosinophils compared to testudines (37.5%). Conversely, basophils were elevated in crocodilia (27.2%) compared to squamata and testudines (9.7% and 0.6% respectively) [Figure 1bii]. Suggesting that within closely related orders that share similar habitats or environments, individual species favoured different configurations of myeloid WBCs. Eosinophils were present as a small population in the mammalian lineages, again variation between the respective classes was observed as monotremata had a smaller proportion of cells (6.7%) compared to the marsupial and placental orders, which had similar values (8.2% and 10.0% respectively). In contrast, compared to all the other classes, basophils seemed to have been largely lost from the mammalian classes of monotremata (0.7%) and marsupials (0.6%), with a minority population observed in placental animals (1.7%) [Figure 1bii].
Monocytes are blood agranulocytic cells that are closely related to neutrophils sharing both growth and survival factors, early development pathways with common progenitors, and in some instances, effector functions in the respective terminally differentiated cells (Evrard et al., 2018; Yáñez et al., 2017; Halene et al., 2010). Interestingly, the monocyte population peaked both in terms of number and proportion within the reptilian class, in the orders; squamata (5.06 × 109/L and 25% respectively) and testudines (1.2 × 109/L and 19.8% respectively). The azurophil – a specialised myeloid cell population- peaked in the crocodilian lineage (1.5 × 109/L and 10.5% respectively), however; it remains unclear whether to classify the azurophil as a distinct cell type, as we have done, or as a subset of monocyte or neutrophil. Aves is the most closely related phylogenetic group to reptiles as they share a recent common ancestor. Although aves has lost the azurophil, they retain a small population of monocytes (4.9%). Interestingly, a similar proportion is observed in the more distantly-related mammalian class, across all three orders (monotremata (8.8%), marsupial (8.5%) and placental (6.9%)), suggesting that the advent of endothermy may have played a role in the distribution of monocytes in warm-blooded animals [Figure 1bii] In summary, we show in phylum Chordata that the majority of peripheral myeloid blood cells are comprised of granulocytes, although the proportion of basophils, eosinophils and neutrophils varies according to the respective orders and lineages. By the advent of mammals and birds however, the neutrophil has become the predominant granulocyte of the blood.
Interrogation of CSF1/CSF1R and CSF3/CSF3R and C/EBP gene family reveals their loss in the early lineages
The development and maturation of a neutrophil from a multipotent progenitor to a lineage committed post-mitotic cell is driven by CSF3 working in concert with a number of transcription factors, including Pu.1, GFI1 and Runx, whose roles have been well established (Yvan-Charvet and Ng, 2019; Halene et al., 2010). One of the most heavily involved family of transcription factors in neutrophil development is the CCAT enhancer binding-protein family (C/EBP), comprising six members (C/EBPα, β, δ, γ, ε, and τ) (Halene et al., 2010; Zhang et al., 1997; Yamanaka et al., 1997; Lekstrom-Himes and Xanthopoulos, 1999). Multiple studies have shown that C/EBPα, C/EBPβ, and C/EBPε are required for different stagesof granulopoiesis and are further detailed in Supplementary file 2. Briefly, under steady state conditions c/ebpα-/- mice are neutropenic as there is a block in the maturation stage between the common myeloid progenitor (CMP) and the granulocyte- monocyte progenitor (GMP) (Smith et al., 1996; Zhang et al., 2004). In contrast, c/ebpε-/- mice lack functionally mature neutrophils as the stage between GMP and terminal differentiation is blocked (Morosetti et al., 1997; Chih et al., 1997). Thus, demonstrating that c/ebpα and c/ebpε have distinct roles in the early (myeloid progenitor progression) and later (maturation) stages of neutrophil development (Halene et al., 2010; Zhang et al., 1997; Yamanaka et al., 1997; Lekstrom-Himes and Xanthopoulos, 1999; Smith et al., 1996; Zhang et al., 2004; Morosetti et al., 1997; Chih et al., 1997). C/ebpβ is required during emergency granulopoiesis as deficient mice could not mobilise a response to systemic fungal infection or cytokine stimulation, even though steady state granulopoiesis remained unaltered (Zhang et al., 2002; Hirai et al., 2006; Screpanti et al., 1995; Tanaka et al., 1995). Having established, the presence of the neutrophil/heterophils across the phylum Chordata, we then evaluated the presence of CSFR/CSF and C/EBP genes in the different animal classes to further understand the emergence of granulopoiesis in phylum Chordata.
Gene data for fifty-nine species was collected from the NCBI Gene and Ensembl databases for the following genes; CSF1R, IL34, CSF1, CSF3R, CSF3, C/EBPα, C/EBPβ and C/EBPε and used to generate a heat map. Interestingly, CSF1, CSF1R, CSF3R and CSF3 are absent from the cartilaginous fishes; the Australian Ghost shark and the Whale shark (Figure 2a). In contrast, the Coelacanth - another early lineage - retained CSF1R, IL34, CSF3R and CSF3 genes. CSF1 was absent from all the jawed fishes although IL34 was present. CSF1R was the only the receptor present in both the lamprey and hagfish, and CSF3R and all other ligands were absent (Figure 2a). The Jawed/Jawless fish accounted for the largest loss of genes within in an animal order as generally; CSF1R, IL34, CSF1, CSF3R, and CSF3 were largely all present in amphibia, reptilia, aves, and mammals. The exception being in the reptilian order of squamata, where IL34 had been lost. There were also examples of gene loss at the species level, CSF1 and CSF3 were absent from the Ostrich; however, they had not been lost overall in from the group (Figure 2a).

Analysis of chordate CSF1/CSF1R, CSF3/CSF3R, and neutrophil-related transcription factors reveals their loss in Jawed/Jawless fish lineages.
(A) Heat map of CSF1R/IL34/CSF1, CSF3R/CSF3, and neutrophil-related transcription factors in selected members of Phylum Chordata. (B) Syntenic maps of CSF1R/IL34/CSF1 and CSF3R/CSF3 in selected Jawed fish compared to human.
During granulocytic differentiation, there is temporal expression of C/EBP gene family members according to the maturation stage of the neutrophil. C/EBPα is expressed highest in early progenitors and decreases through differentiation to low levels of expression in mature granulocytes. The inverse expression pattern is seen with CEBPε, where it is highly expressed in mature neutrophils and there’s minimal expression in progenitors. C/EBPβ follows a different pattern with expression increasing in progenitor populations and levels being maintained through to maturation (Scott et al., 1992). Given the importance of the three members to mammalian neutrophil development we considered the conservation of these members across phylum Chordata. Interestingly, within the early Jawed/Jawless lineages there had been local loss of C/EBPα, C/EBPε or both from some members, however all fish examined retained C/EBPβ and the Coelacanth had all three genes, indicative of their presence in a shared common ancestor (Figure 2a). The analysis showed very strongly supported evolutionary conservation of the C/EBP gene family members in tetrapod lineages, as almost every group had a minimal combination of C/EBPα or C/EBPβ and C/EBPε (the exception being the aves class where CEBP/ε was absent). Interestingly, within the aves class both ostriches and turkeys lacked all the gene family members examined. Taken together, these results suggested that there were local factors determining the retention and loss of genes as well as a strong level of redundancy among the C/EBP gene family. A high degree of functional redundancy between family members has been well described and is supported by various knockout studies in mice. These show that at different stages of granulopoiesis, C/EBP gene family members - particularly C/EBPβ- can in part compensate functionally for each other (Jones et al., 2002; Akagi et al., 2010).
Interestingly, within the early lineages of the jawed fishes, IL34 was observed in two species and the Coelacanth was the only species to have CSF1R, CSF3R and CSF3 present. As these three distinct lineages share a common ancestor, this would indicate that there had been local gene loss in both the Whale shark and Australian Ghost shark. CSF1 was absent from all lineages and it was unclear as to whether this absence was due to loss or that CSF1 had not yet evolved. To further address this, we used a syntenic approach (Engström et al., 2007), that takes advantage of a process where genes that have co-evolved together physically co-localise within the loci, and using this to manually map out the orthologous gene locations in the Jawed fish. Syntenic maps were generated for the CSF1R (CSF1R, IL34, CSF1) and CSF3R (CSF3R and CSF3) family of genes using the human orthologue as a reference point (Figure 2b).
In humans, the CSF1R gene is located downstream of the HMGBX3 gene and upstream of PGDFRB (a gene paralogue of CSF1R) and CDX1. A similar orientation is observed in the Coelacanth, where CSF1R is downstream of HMGBX3 and adjacent to PDGFRB and CDX1. The orientation is inversed in the Australian Ghost shark, where PDGFRB and CDX1 co-localise together upstream but HMGBX3 and CSF1R have been lost (Figure 2b). Interestingly in the Whale shark, CSF1R and all the flanking genes are absent, suggesting that this section in its entirety may have been lost (Figure 2b). Human IL34 co-localises with FUK, COG4, and SF3B3 upstream, and the MTSS1L and VAC14 genes immediately downstream. Both the Coelacanth and the Whale shark have partially retained the syntenic combination but not in the same location. IL34 is situated immediately downstream of COG4 and SF3B3, whereas MTSS1L and VAC14 are located elsewhere on the Coelacanth chromosome (Figure 2b). In the Whale shark, IL34 is immediately adjacent to MTSS1L and these genes are both downstream of FUK and COG4, which have co-localised next to each other (Figure 2b). The gene arrangement of the Australian ghost shark most closely resembles the human, as FUK, COG4, SF3B3, MTSS1L and VAC14 all co-localise to the same region on the chromosome and the absence of IL34, suggest it has been loss in the process of a local gene rearrangement (Figure 2b). The CSF1 human orthologue has formed a contiguous block with AHCYL1 and STRIP1 and is flanked upstream by AMPD2 and EPS8L3. Interestingly although CSF1 is absent from all members, the flanking genes have largely been retained. Similar to the human arrangement, in the Coelacanth; AMPD2 and EPS8L3 co-localised together, while AHCYL1 and STRIP1 are located in close proximity in a different location (Figure 2b). In both the Australian Ghost shark and Whale Sharks, AHCYL1, and AMPD2 are located together, STRIP1 is located elsewhere in the Australian Ghost shark and has been lost entirely from the Whale shark (Figure 2b). These lineages arose early in evolution but have very similar synteny structure to the later emerging human chromosome, suggesting that the CSF1 gene entered this location in an ancestor that emerged after the jawed fish.
A similar pattern emerged in the CSF3 family of proteins when analysed. In humans, the CSF3R gene is located downstream of STK40, LSM10, OSCP1, and MRPS15. The Coelacanth is the only species to retain CSF3R and that is located immediately adjacent to MRPS15, however STK40, LSM10, and OSCP1 have been lost. In contrast, although both the Whale shark and Australian Ghost shark have presumably lost the CSF3R gene, they have retained the four other flanking genes, either in one location as in the Australian Ghost shark, or distributed along the chromosome, as in the Whale Shark (Figure 2b). Human CSF3 co-localises downstream of GSDMA and PSMD3 and is adjacent to MED24. Again, a similar arrangement is observed in the Coelacanth, with CSF3 flanked by PSMD3 and MED24 upstream and downstream respectively. As before, the Whale shark and Australian Ghost shark are similar in their gene arrangements as GSDMA and PSMD3 have been lost, and the only gene retained in the Australian Ghost shark is MED24 (Figure 2b). Intriguingly, in the Whale shark, MED24 co-localises with an IL6-like gene, which could be a functional paralogue of CSF3 (Figure 2b). These results suggest there is more than one receptor-ligand family involved in the development and maturation of heterophils and neutrophils in the early lineages. Taken together, these data suggest that CSF3R/CSF3 and CSF1R were present in a common ancestor to early lineages and there have been some local losses in selected members. CSF1 appears to have evolved independently of its cognate receptor, after the emergence of Jawed lineages and prior to the advent of the tetrapod lineages.
Analysis of chordate orthologous protein homology further supports the ancestral pairing of CSF1R/IL34 and CSF3R/CSF3 in early lineages
The syntenic analysis established the presence of the CSF1R and CSF3R gene families in Chordates. Interestingly, CSF1R/IL34 and CSF3R/CSF3 had already evolved by the emergence of Chordates, as evidenced by their existence in some of the early lineages of Jawed/Jawless fish. However, the gene data in isolation did not provide a complete understanding and additional analysis was needed. Orthologous proteins are considered to have the same function in different species and therefore, it’s broadly assumed that the proteins will largely be conserved at the primary and structural levels (Forslund et al., 2011). To further elucidate the evolutionary process, we compared the shared sequence similarity of the orthologous CSF1R and CSF3R protein families, as this data is widely available. The protein sequences for the following proteins; CSF1R, IL34, CSF1, CSF3R, and CSF3 for fifty-nine species were collated from the NCBI and Ensembl databases. The shared sequence similarity for each individual sequence was generated using the NCBI BLAST engine tool and the relevant human orthologue submitted as the query. The resulting data was then plotted in groups as identified by animal order and visualised in a bar chart (Figure 3a).

Shared sequence similarity analysis of Chordate CSFR/CSF protein homology further supports the ancestral pairing of CSF1R/IL34 and CSF3R/CSF3 in early lineages.
(A) Representative plots of % shared sequence similarity for CSF1R, CSF3R, IL34, CSF1, and CSF3 in the respective sub-groups of Phylum Chordata, data is visualised as a floating bar and the line represents the mean value. Baseline shared similarity is represented by the red line. (B) Graphical representation of % shared sequence similarity in sub-groups of Phylum Chordata versus time.
The level of shared sequence similarity observed in the receptors, and therefore conservation of the primary protein sequence, varied across the different animal orders. The highest levels of shared sequence conservation for both CSF1R and CSF3R were observed in the placental mammals (CSF1R; 86.7%, CSF3R; 81.7), which would be anticipated as the human is a member of this group. The lowest level of shared sequence similarity was observed in the Jawed/Jawless fish (CSF1R; 43.1%, CSF3R; 33.0%) (Figure 3a). Interestingly, the CSF1R and CSF3R protein sequence values from early mammals, monotremata (CSF1R; 60.4%, CSF3R; 46.9%) and marsupalia (CSF1R; 60.4%, CSF3R; 51.0%), ranged between the earlier orders; aves (CSF1R; 51.5%, CSF3R; 39.4%), testudines (CSF1R; 53.4%, CSF3R; 41.2%), crocodilia (CSF1R; 52.6%, CSF3R; 37.6%), squamata (CSF1R; 47.1%, CSF3R; 37.2%), amphibia (CSF1R; 47.2%, CSF3R; 35.7%) and placental mammals (Figure 3a). While there is not a clear consensus as to what the minimum percentage of shared sequence similarity needed to correlate with conservation of function is, it is notable that in this dataset for both CSF1R and CSF3R - independently of each other - the baseline value of shared sequence similarity was approximately 40% (Figure 3a).
The syntenic analysis also demonstrated that of the three ligands, IL34 emerged the earliest and this was reflected in the protein data. As observed with the receptors, the highest level of shared sequence similarity is in the placental mammal group (83.3%), and the lowest in the jawed/jawless fish (38.8%) (Figure 3a). The IL34 protein sequence values from early mammals; monotremata (50.0%) and marsupalia (54.8%), grouped very closely with the earlier orders of; aves (48.5%), crocodilia (46.1%), testudines (43.6%), amphibia (38.0%) and the jawed/jawless fishes (Figure 3a). IL34 was absent from the members of squamata examined. CSF3 is the next best conserved ligand, which again agreed with the synteny data. The highest level of shared sequence similarity was in the placental mammal group (84.2%), and the lowest in the jawed/jawless fish (30.7%) (Figure 3a). The CSF3 protein sequence values were spread among the orders from monotremata (57.2%) and marsupalia (61.1%) to aves (41.0%), crocodilia (35.4%), testudines (39.7%), squamata (30.4%) and amphibia (32.9.1%) (Figure 3a). Finally, CSF1, which is completely absent in jawed/jawless fish, had similar values to CSF3, in terms of shared sequence similarity across the respective groups of; placental mammals (79.4%), monotremata (39.4%), marsupalia (46.5), aves (31.3%), crocodilia (28.0%), testudines (34.5%), squamata (32.6%) and amphibia (28.4%). Intriguingly, the 40% baseline is applicable to the ligand data. In IL34, which is the oldest ligand, the shared sequence similarity for the majority of the groups was above 40%. Whereas for both CSF3 and CSF1, which emerged later, in the majority of the groups the shared sequence similarity was under 40% (Figure 3a).
The analysis of orthologous CSF1R and CSF3R protein families illustrated that the mammalian proteins- largely within the placental mammals- had changed considerably compared to the other lineages and therefore were excluded from subsequent analysis. To further interrogate how the respective CSF1R and CSF3R families co-evolved in the early lineages, the shared sequence similarities for each order were plotted against time (Figure 3b). Interestingly, the trajectories for CSF1R and IL34 and CSF3R and CSF3, largely tracked to each other. This suggested that they evolved at a similar pace across the same period of time and is indicative of evolutionary pressure restricting changes to cognate receptors and ligands. As expected, CSF1 did not have the same restrictive pattern in the earlier lineages as it emerged later (Figure 3b). These results supported the early emergence of CSF3R, CSF1R, and IL34. Interestingly, although CSF3 developed later than CSF3R, the two have co-evolved in step together. In contrast, while CSF1 would appear to be the principal ligand of CSF1R in mammals, the data supported the ancestral pairing of IL34 and CSF1R in early lineages. However, it should also be noted that IL34 has two other receptors, Sydecan (CD138) (Segaliny et al., 2015) and protein-tyrosine phosphatase zeta (Nandi et al., 2013) and these interactions could also account for the earlier emergence of IL34.
The emergence of CSF3R/CSF3 and onset of endothermy likely influenced the distribution of neutrophils in Chordates during evolution
The mammalian neutrophil shares many developmental and functional properties with counterparts present in other animal orders such as phagocytosis, oxidative burst, degranulation and cell motility (Genovese et al., 2013; Styrt, 1989). However; a notable species-specific difference is the distribution of the neutrophil/heterophil within peripheral blood, which suggests that evolutionary pressures might also be involved. To further address this, we reconstructed a set of related timescales plotting the percentage distribution of the chordate myeloid WBCs described in Figure 1 and the percentage of protein shared sequence similarity of chordate CSFR/CSF described in Figure 3 versus time. To simplify the model, the timescales were plotted chronologically, where applicable, and the emergence point of the earliest member within an animal order was used. The timescales were then divided into two distinct periods; the first focussed on the early lineages of jawed/jawless fish, amphibia and early reptilia and the second focussed on all reptilia, aves and mammalia (Figure 4).

Changes in chordate blood granulocyte distribution from Jawed/Jawless lineages through to placental mammals are multi-factorial and likely driven by the emergence of CSF3 and onset of endothermy.
(A) Graphical representation of % population distribution of myeloid white blood cells versus time (i) and % shared sequence similarity of CSF1R and CSF3R protein families versus time for Jawed/Jawless fish, Amphibia and the reptilian order of Squamata (ii). (B) Graphical representation of % population distribution of myeloid white blood cells versus time (i) and % shared sequence similarity of CSF1R and CSF3R protein families versus time for the Reptilian orders of squamata, testudines and crocodilia, Aves, and the Mammalian orders of monotremata, marsupalia and Placental (ii). Error bars represent standard deviation.
The synteny studies demonstrated that CSF1R and IL34 were already in existence prior to the appearance of the jawed/jawless fish lineages. However, CSF3R and CSF3 were only present in the Coelacanth and although an IL6-like paralogue was observed in the Whale shark, a cognate receptor was not identified. These results suggested there was a greater level of diversity in the growth factors responsible for WBCs in these lineages and is reflected in the distribution of WBCs in Jawed/jawless fish, where heterophils, neutrophils, monocytes, eosinophils and basophils are all present. While heterophils and monocytes largely dominated, there was a more even distribution of minority populations of neutrophils, basophils and eosinophils (Figure 4ai). By the advent of amphibia, distinct neutrophil and heterophil populations had seemingly been lost in favour of a single neutrophil or heterophil population. The amphibian order was the only lineage to demonstrate a largely even distribution of neutrophils, eosinophils, basophils, and monocytes, although there was a slight decline compared to the levels observed in Jawed/Jawless fish (Figure 4ai).
At some stage during tetrapod evolution after the emergence of CSF3 and CSF1 haematopoiesis largely transitioned to the tissue-specific compartment of the bone marrow, which likely had implications for WBC distribution in subsequent orders. In early reptilia, represented by the squamata order of lizards, the neutrophil becomes the dominant granulocyte and there is concomitant reduction in basophils and eosinophils. Interestingly, there is an increase in monocyte populations, which coincides with the emergence of the monocyte/macrophage specific growth factor CSF1 (Figure 4aii). These changes occur in ectothermic species suggesting that this is more in response to cell-intrinsic factors, rather than external factors such as environment (Figure 4a).
The reptilian lineage comprises three orders; squamata, crocodilia and testudines, who emerged over a large timescale. Accordingly, there are many intra-order differences such as the environments in which members live; that is in the water versus on land or the presence or absence of limbs or exoskeletons. Interestingly, the neutrophils seemingly peaked in squamata as proportions steadily declined with the lowest levels observed in testudines (Figure 4bi). In contrast, eosinophil levels, which remained low across the squamata and crocodilia classes, peaked in testudines. However, across this period in time; CSF1R/CSF1 and CSF3R/CSF3 had become established as blood myeloid cell-specific factors and were unlikely to be driving the changes in WBC distribution (Figure 4bi). Notably, as the neutrophil populations declined, both the basophil and eosinophil populations increased, with basophils peaking in the crocodilia lineage before a dramatic decline in testudines. Suggesting, environmental pressure on the immune response of certain orders favoured either the basophil or the eosinophil at the expense of neutrophils and the other granulocyte (Figure 4bi). The decline was not only restricted to neutrophil population, as the monocyte population had a similar pattern. Monocytes were at their highest levels in squamata before a rapid decline in crocodilia, however the population expands again in testudines. Notably, as the azurophil emerged in crocodilia, monocytes declined, suggesting that they effectively took over the role of monocytes in crocodilia (Figure 4bi).
Interestingly, by the advent of testudines, neutrophil/heterophils were no longer the dominant granulocyte, although granulocytes overall still represented the largest proportion of myeloid WBCs (Figure 4bi). However, as the orders of aves, monotremata and marsupalia emerge there is a noticeable change. A limitation of visualising data in this form is the implication of the distinct sequential evolution of different animal orders. However; it’s likely there would have been a degree of overlap between the emergence of distinct orders. Therefore, it is noticeable that there was a rapid increase in the both avian heterophil population and monotremata neutrophil population compared to testudines (Figure 4bi). The increase in neutrophils was mirrored by an equally rapid decrease in eosinophils in both the aves and mammalian orders (Figure 4bi). Neutrophils and eosinophils also share early common progenitors as part of the development pathway in the bone marrow, therefore neutrophil production appeared to be dominating the cell machinery. Interestingly, the switch to the predominance of the neutrophil over the eosinophil, coincides with the emergence of endothermy, strongly suggesting that external factors are behind these changes in myeloid WBC distribution as there are no noticeable changes in the distribution of the CSF1R/CSF1 and CSF3R/CSF3 ligand/receptor pairings. However, the impact is filtering down to the gene/protein level as there are changes in the shared sequence similarity of CSF3, as it increases more rapidly between Testudines and Monotremata, than in the period between Squamata and Testudines. Interestingly, these changes were not reflected in the CSF3R trajectory (Figure 4bii). Presumably, the arrival of endothermy resulted in the appearance of novel pathogens for which a neutrophil-mediated response was more appropriate than an eosinophil-mediated one, for example, neutrophils have an array of killing mechanisms such as superior phagocytosis capabilities and preferential production of intracellular respiratory burst that favour intracellular killing of smaller targets (Borregaard and Cowland, 1997; Kovács et al., 2014; Hatano et al., 2009). In contrast, eosinophils favour extracellular killing by deploying granules such major basic protein (MBP) and eosinophil cationic protein (ECP) that are toxic to parasites such as helminths (Shamri et al., 2011). These results suggest that in response to the emergence of endothermy and presumed associated novel pathogens, the mammalian host system has selected for the neutrophil as the predominant granulocyte of the blood. The monocyte population had re-emerged as a single population by the appearance of aves and the distribution remained consistent in both the avian and mammalian orders (Figure 4bi). Intriguingly, the greatest period of change for myeloid WBC distribution is between testudines and monotremata. However, the equivalent period for protein homology happens much later between the mammalian orders of marsupalia and placentalia (Figure 4bii) and is common across all CSFR/CSF pairings. Therefore, as WBC function and distribution had previously been established in early lineages, this suggested that another external factor was responsible for the rapid change, such as the emergence of internalized pregnancy (Figure 4bii).
Discussion
The mammalian neutrophil is a highly specialised cell that acts as a first responder to insults against a host immune system as well as acting in an equally important sentinel role (Hidalgo et al., 2019; Ng et al., 2019). They are functionally conserved across phylum Chordata and constitute the largest population of myeloid cells in the blood of birds and mammals. In humans, for example, up to one billion neutrophils per kilogram of body weight are produced in the bone marrow each day (Yvan-Charvet and Ng, 2019). The immune response has evolved in such a way as to be able to respond efficiently to a variety of threats, and it is interesting that a resource such as the bone marrow should be expended on the production and maintenance of the relatively short-lived neutrophil at the expense of the eosinophil and basophil. The timescale modelling demonstrates that prior to the advent of tetrapoda lineages, the neutrophil was in a pool of different WBC populations at the disposal of the jawed/jawless fish. However, the appearance of CSF3, altered the distribution and with each emerging animal order, a different granulocyte was favoured, presumably, for cell-specific adaptations in response to environmental challenges. Thus, suggesting that CSFR1R/CSF1 and CSF3R/CSF3 signalling conferred adaptations on the neutrophil that proved evolutionarily advantageous. The work discussed here focusses on the mammalian neutrophil within the context of blood, further work is required to define similar aspects as regards to other tissue resident populations such as macrophages or granulocytes.
The bone marrow is an essential site for granulopoiesis and myelopoiesis. As multipotent haematopoietic stem cells (HSCs) progress through different vascular niches within the bone marrow they sequentially lose their potential to form other lineages in response to environmental cues from bone marrow dwelling macrophages and stromal endothelial cells (Ng et al., 2019; Yvan-Charvet and Ng, 2019; Greenbaum and Link, 2011). The emergence of CSF1 in tetrapod lineages likely led to the appearance of bone and bone marrow. A consequence of this was the gradual re-organisation of haematopoiesis away from existing haematopoietic tissues, such as the eosinophil rich Leydig organ of sharks to tissue-specific compartments in the bone marrow. CSF1, in conjunction with another factor, receptor activator of nuclear factor-κB ligand (RANKL), co-ordinate the reabsorption of old bone through haematopoietically derived osteoclasts to allow the generation of new bone. CSF1 and RANKL (which emerged at a similar point in evolution) orchestrate bone-remodelling through their respective receptors, CSF1R and RANK (Kim and Kim, 2016). This suggests that by the time a bone structure had evolved in amphibia, the bone marrow had become the principle site of haematopoiesis, although there are some exceptions within various anurans (de Abreu Manso et al., 2009).
In contrast to peripheral blood, CSF3 is expressed on a number of cells in the bone marrow including; neutrophils, monocytes, B cells, myeloid progenitors and HSCs (Greenbaum and Link, 2011; Petit et al., 2002; Semerad et al., 2005; Lévesque et al., 2003). Although CSF3 can likely act directly on HSCs through their receptor, it is believed to indirectly mobilise HSCs through a monocytic intermediary that secretes CSF3, which leads to suppression of the CXCR4:CXCL12 axis, alteration of the bone marrow niche, and the subsequent release of HSCs (Greenbaum and Link, 2011; Petit et al., 2002; Semerad et al., 2005; Lévesque et al., 2003). In a similar process, CSF3 can suppress B cell lymphopoiesis by again targeting CXCL12 and suppressing other B cell tropic factors or stromal cells that favour the lymphoid niche (Day et al., 2015). Thus, the emergence of CSF3R/CSF3 conferred the adaptation, or advantage, of control of the biological machinery i.e. CSF3 provides a mechanism through which haematopoiesis can be shaped and deployed in favour of maximal neutrophil generation. This broadly aligns with the neutrophil starting to dominate peripheral populations in the transition from amphibia to the lizards of early squamata following the appearance of bone.
Mammalian neutrophil production occurs in the haematopoietic cords present within the venous sinuses and the daily output is approximately 1.7 × 109/kg (Summers et al., 2010). As our data demonstrate high levels of neutrophil output are common across warm-blooded chordates and neutrophil/heterophils are the predominant granulocyte in birds and mammals, in what could be considered an example of convergent evolution. One of the fundamental requirements of the host innate immune response is to be able to respond rapidly to perceived threats that can be present at any site in the body, which requires the cellular arm to be constantly present and available. Theoretically this can be achieved in the steady state by having either low volumes of long-lived cells or high volumes of short-lived cells circulating in the periphery. As a consequence, there are potential biological trade-offs when considering each setting, firstly; the longer a cell survives in the periphery, the more effort is required by the host to maintain its survival in terms of providing appropriate cues and growth factors. Secondly, while a short-lived cell does not need as much host input for survival, a high turnover is required in order to ensure it doesn’t compete for growth factors with other cell types or cause damage by being retained beyond its usefulness. The latter setting fits with the observed neutrophil life cycle and could be considered an evolutionary adaptation. The CSF3R/CSF3 signalling pathway is essential for generating high neutrophil numbers without adverse effects to the host.
The neutrophil has a short circulatory half-life in the steady state of approximately one day (Ng et al., 2019; Yvan-Charvet and Ng, 2019), although this does remain controversial as some estimates of the circulatory lifespan are as much as five days (Pillay et al., 2010). By contrast, a mature eosinophil has a short circulatory lifespan of approximately 18 hr before relocating to the tissue, where it will survive for a further 2–5 days (Park and Bochner, 2010). Similarly, intermediate and non-classical monocytes can survive in the periphery for four and seven days respectively (Patel et al., 2017). In human studies, CSF3 has been shown to ‘effectively’ shorten the lifespan of neutrophil myeloblasts and promyelocytes by decreasing their time spent in cell cycle and accelerating their progress to maturation (Lord et al., 1989). Although, it is unclear what the exact effect is on the lifespan of post-mitotic neutrophils, studies have shown that the addition of exogenous CSF3 can delay apoptosis in mature neutrophils (van Raam et al., 2008; Zimbelman et al., 2002). This suggests that under certain conditions; such as infection, CSF3 can extend the lifespan of a mature neutrophil. Interestingly, classical monocytes also express CSF3R and have a short circulating lifespan of approximately 24 hr (Patel et al., 2017) These studies suggest that under steady state conditions and below a certain threshold, CSF3 has an as-yet-undocumented role either directly or indirectly in maintaining a short neutrophil lifespan and thus allowing the efficient turnover required to sustain high levels of granulopoiesis.
CSF3 orchestrates the life cycle of neutrophils in the bone marrow microenvironment by marshalling neutrophil progenitors through different development stages through to maturity, which is indicative of the inductive model of differentiation (Stanley, 2009). Accordingly, in response to this, CSF3R is expressed though every life stage, although at differing levels across the neutrophil populations, with the highest levels observed on mature cells where it is expressed at between two and three-fold more than on progenitors (Demetri and Griffin, 1991). CSF3R/CSF3 signalling has a dual role in controlling the distribution of neutrophils by both retaining a population of mature neutrophils in the bone marrow as a reservoir and facilitating the egress of other neutrophils into the periphery. Thus, a major evolutionary adaptation that CSF3R/CSF3 has conferred to the neutrophil is the ability to move, both on a population-wide and individual cell level.
There is abundant production of neutrophils daily in the bone marrow; however, there are far fewer neutrophils in circulation in the blood than are produced during granulopoiesis as the total neutrophil population is effectively stored in the bone marrow or marginated in intravascular pools within the spleen and liver (Ussov et al., 1995). Neutrophils can transit between sites in response to CSF3 signalling and it is estimated that 49% of cells are present in the circulating pool and the remaining 51% are marginated in discrete vascular pools (Summers et al., 2010; Athens et al., 1961). In the event of an infection, neutrophils can be mobilised from the marginated pools and bone marrow in response to CSF3-induced production of mobilising signals such as CXCL1 (Köhler et al., 2011). The co-ordinated egress of neutrophils from the bone marrow is achieved by CSF3 interacting with the CXCR4/CXCL12 and CXCR2/CXCL2 signalling pathways. CSF3 disrupts the CXCR4-CXCL12 retention axis by reducing CXCL12 release from endothelial stromal cells and reducing CXCR4 expression on neutrophils, thus allowing the movement of mature neutrophils through the venous sinuses to the periphery. CSF3-induced expression of CXCR2 on neutrophils then causes their migration to the vasculature along a CXCL2 chemotactic gradient (Adrover et al., 2019; Eash et al., 2010; Eash et al., 2009).
As neutrophils enter the bloodstream, they need to be able to migrate easily around the circulatory system under high flow conditions. This requires the neutrophil’s physical form to have high deformability and flexibility as it encounters the different diameters of the vasculature. The neutrophil nucleus is functionally adapted to this role because of its multi-nucleated structure and the distinct protein composition of the nuclear envelope, features that are widely conserved across mammalian species (Manley et al., 2018). CSF3 in co-ordination with C/EBPε and the ETS factors; Pu.1 and GA-binding protein (GABP) are responsible for the transcriptional control of the essential neutrophil nuclear structural proteins; Lamin A, Lamin C and Lamin B receptor (LBR) (Malu et al., 2016; Cohen et al., 2008). In comparison to other cell nuclear protein compositions, neutrophils have a low proportion of Lamin A and Lamin C, which is believed to make the nucleus more flexible for easier transit (Manley et al., 2018; Malu et al., 2016; Cohen et al., 2008). In contrast, the levels of LBR are increased, which is required for nuclear lobulation and subsequent neutrophil maturation (Manley et al., 2018; Malu et al., 2016). However, recent work has also demonstrated that band cells (immature neutrophils) are capable of just as efficient migration with incomplete nuclear segmentation as mature neutrophils (van Grinsven et al., 2019). Taken together, these factors support the argument that CSF3 signalling and C/EBPε are intrinsic to the evolution of mammalian neutrophil cell motility.
The evolution of CSF3R/CSF3 has been essential to the development of the neutrophil/heterophil in chordates and its own existence as the principle neutrophil growth factor is evolutionarily advantageous. Some members of the jawed/jawless lineages are unique in that they have populations of both neutrophils and heterophils, whereas later linages favour either neutrophils or heterophils. Our analysis shows that CSF3R/CSF3 emerged before the advent of jawed/jawless fish as both are present in the Coelacanth and absent from the other lineages. Intriguingly, an IL6-like gene was present in the same syntenic location of the whale shark, suggesting the possibility that heterophils and neutrophils were independently controlled by an IL6-like protein and CSF3 in the early lineages. IL6 is an important pleiotropic pro-inflammatory cytokine that plays key roles in infection, inflammation and haematopoiesis (Tanaka et al., 2014). Although CSF3R and IIL6R, which are functional paralogues, diverged from each other many millions of years ago, there is still functional redundancy at the cell level in modern-day mammals, as neutrophils are present - though at vastly reduced numbers - in CSF3-/- and CSF3R-/- mice (Liu et al., 1996; Lieschke et al., 1994). IL-6 can act on neutrophil progenitors and immature neutrophils thus supporting neutrophil development in the bone marrow of CSF3-deficient mice (Basu et al., 2000). However, mature neutrophils are refractory to IL-6 signalling as the expression of the IL-6R subunit gp130 is lost during maturation (Wilkinson et al., 2018). From amphibia onwards, the protein analysis here suggests that CSF3R and CSF3 co-evolved closely together and in contrast to CSF1R/CSF1/IL34, CSF3 is likely the only ligand of CSF3R.
Using various in-silico approaches, these findings have demonstrated how essential CSF3R and CSF3 (and to a lesser extent CSF1R and CSF1) are to the predominance of mammalian neutrophils in blood. Through the course of evolution, CSF3R/CSF3 signalling has accrued many properties that are responsible for the survival of the neutrophil and its ability to function, such as cell motility and mobility. Interestingly, although CSF3 is present in jawed/jawless fish lineages, it’s not until the emergence of tetrapoda that CSF3R and CSF3 begin to dominate haematopoiesis. Our current studies have been limited to chordates in which the data is publicly available and as the datasets continue to be expanded and updated, our hypothesis may need to change to reflect that. However, we would argue that this approach is a valid method for exploring the origins of the neutrophil granulocyte and would further be supported by the physical isolation and characterisation of neutrophils, associated genes and proteins in the different animal orders within Phylum Chordata.
We also envisage exploring the phagocytic granulocyte of two model species. The amphioxus is considered the basal chordate and a macrophage-like population has been identified, although a neutrophil/heterophil has yet to be described within the rudimentary circulatory system (Han et al., 2010). The lungfish, a primitive airbreathing fish, is unique among jawed/jawless lineages as it lives in freshwater and can survive on the land for up to one year. Accordingly, the lungfish has many adaptations and is considered to be a one of the closest living relatives to tetrapods making it a good model species for further study (Takezaki and Nishihara, 2017). Given how essential motility and mobility are to neutrophil function and development, it would be useful to discern when it emerged in evolution by identifying if equivalent cells and functional gene/protein orthologues are present in either species. These comparative studies would answer fundamental questions about the origin of the neutrophilic phagocyte.
Materials and methods
Species selection
Request a detailed protocolNinety-four animals from all the major classes were selected, thirty-five of which, were used for the calculation of haematological parameters and the remaining fifty-nine for the bioinformatics-based studies. Analyses were performed on non- mammalian lineages including; jawed/jawless fish, amphibia, non-avian reptiles – both non-crocodilian and crocodilian-, aves, and the mammalian lineages; monotremata, marsupialia and placentalia. Urochordates (tunicata), and cephalochordate (Amphoxi) were excluded from analysis as there was insufficient coverage or insufficient annotation of sequence data in the Pubmed Gene and Ensembl databases. Similarly, the lungfish (Dipnoi) was also excluded from this analysis because of insufficient coverage of sequence data. Finally, teleost lineages were excluded from the analysis owing to having gone through three rounds of genome wide duplication, in contrast to all other chordates who have only undergone two rounds (Glasauer and Neuhauss, 2014). A species tree was generated using the NCBI Common Taxonmy browser common tree tool (Sayers et al., 2009; Benson et al., 2009) and visualised using the Interactive Tree of Life (iTOL) web browser tool (Letunic and Bork, 2019).
Comparative analysis of haematological parameters and CSF/CSFR gene presence and synteny in chordates
Request a detailed protocolComplete blood count (CBC) data for thirty-eight species representing each animal order or class were collated from existing literature to perform a meta-analysis of myeloid WBC in phylum Chordata. Where available, a representative data for each gender within in species class or order were used. Only counts for sub-adult or adults that were calculated to the standardised concentration (109/L) or equivalent were used and all lymphoid data was excluded (references for CBC data are listed in Supplementary file 3). The proportional composition for each myeloid subset; neutrophil, eosinophil, basophil and monocytes as part of the overall myeloid population was also calculated.
Gene sequence for fifty-nine species - representing different animal orders and classes - were retrieved from the NCBI Gene databases (Pruitt et al., 2007) in most instances or from annotated entries in Ensembl (v95) (Yates et al., 2020) for the following genes; CSF1R, IL34, CSF1, CSF3R, CSF3, C/EBPα, C/EBPβ and C/EBPε (gene identification numbers for species examined are listed in Supplementary file 4). A heatmap visualising gene presence or absence was generated using maplotlib/seaborn library for Python package 3.7.7.
Syntenic maps were generated manually for selected genes in the Jawed fish lineage. To generate the maps, syntenic blocks were identified as a section of the human chromosome containing the gene of interest flanked by x number of genes. The blocks were then manually compared to similar regions in the chromosomes in Jawed fish to identify orthologous genes. Between three and five genes were chosen per syntenic block to function as anchor points of reference. Anchor points were selected based on their situational proximity, either upstream or downstream, to gene of interest and were present in all species examined. Multiple genes were identified as anchor points to mitigate for the random loss of genes during the process of evolutionary gene rearrangement.
Generation of percentage shared sequence similarity plots and timescales
Request a detailed protocolAs with the gene sequences, protein sequences for the identified species were retrieved from the NCBI Protein databases. Presumed orthologous sequences were screened using the NCBI Basic Local Alignment Search tool (BLAST) [RRID:SCR_004870] to generate a percentage score of sequence similarity (protein accession numbers for species examined are listed in Supplementary file 5). The parameters to determine shared sequence similarity were as follows; a sequence was deemed to be homologous where coverage of the protein sequence was equal to or greater than 40% of the total protein sequence examined and the e value was between 0 × 10−20 - 0. Percentage scores were averaged per order. Shared sequence similarity was also plotted against animal order or class. The same data was also used to generate two sets of timescales for each order/member for a given receptor/ligand family, where either the calculated mean shared sequence similarity value or haematological parameters were plotted on a timescale based on the emergence of the earliest known ancestor of that order versus time (million years ago [mya]).
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
Leukokinetic studies. IV. the total blood, circulating and marginal granulocyte pools and the granulocyte turnover rate in normal subjectsJournal of Clinical Investigation 40:989–995.https://doi.org/10.1172/JCI104338
-
Localization of hematopoietic cells in the bullfrog (Lithobates catesbeianus)Cell and Tissue Research 337:301–312.https://doi.org/10.1007/s00441-009-0803-0
-
CXCR2 and CXCR4 antagonistically regulate neutrophil trafficking from murine bone marrowJournal of Clinical Investigation 120:2423–2431.https://doi.org/10.1172/JCI41649
-
Domain architecture conservation in orthologsBMC Bioinformatics 12:326.https://doi.org/10.1186/1471-2105-12-326
-
Inducible expression of a disease-associatedThe Journal of Biological Chemistry 295:7492–7500.https://doi.org/10.1074/jbc.RA120.012366
-
The avian heterophilDevelopmental & Comparative Immunology 41:334–340.https://doi.org/10.1016/j.dci.2013.03.021
-
Whole-genome duplication in teleost fishes and its evolutionary consequencesMolecular Genetics and Genomics 289:1045–1060.https://doi.org/10.1007/s00438-014-0889-2
-
The primitive immune system of amphioxus provides insights into the ancestral structure of the vertebrate immune systemDevelopmental & Comparative Immunology 34:791–796.https://doi.org/10.1016/j.dci.2010.03.009
-
C/EBPbeta is required for 'emergency' granulopoiesisNature Immunology 7:732–739.https://doi.org/10.1038/ni1354
-
Signaling pathways in osteoclast differentiationChonnam Medical Journal 52:12–17.https://doi.org/10.4068/cmj.2016.52.1.12
-
Neutrophil disorders and their managementJournal of Clinical Pathology 54:7–19.https://doi.org/10.1136/jcp.54.1.7
-
Interactive tree of life (iTOL) v4: recent updates and new developmentsNucleic Acids Research 47:W256–W259.https://doi.org/10.1093/nar/gkz239
-
Disruption of the CXCR4/CXCL12 chemotactic interaction during hematopoietic stem cell mobilization induced by GCSF or cyclophosphamideJournal of Clinical Investigation 111:187–196.https://doi.org/10.1172/JCI15994
-
Systems approach to phagocyte production and activation: neutrophils and monocytesAdvances in Experimental Medicine and Biology 844:99–113.https://doi.org/10.1007/978-1-4939-2095-2_6
-
Comparative pathology of inflammation in the higher vertebrates (reptiles, birds and mammals)Journal of Comparative Pathology 99:1–26.https://doi.org/10.1016/0021-9975(88)90101-6
-
Receptor-type Protein-tyrosine Phosphatase ζ Is a Functional Receptor for Interleukin-34Journal of Biological Chemistry 288:21972–21986.https://doi.org/10.1074/jbc.M112.442731
-
Pathogenesis of ELANE-mutant severe neutropenia revealed by induced pluripotent stem cellsJournal of Clinical Investigation 125:3103–3116.https://doi.org/10.1172/JCI80924
-
Heterogeneity of neutrophilsNature Reviews. Immunology 19:255–265.https://doi.org/10.1038/s41577-019-0141-8
-
Eosinophil survival and apoptosis in health and diseaseAllergy, Asthma and Immunology Research 2:87–101.https://doi.org/10.4168/aair.2010.2.2.87
-
The fate and lifespan of human monocyte subsets in steady state and systemic inflammationJournal of Experimental Medicine 214:1913–1923.https://doi.org/10.1084/jem.20170355
-
Database resources of the national center for biotechnology informationNucleic Acids Research 37:D5–D15.https://doi.org/10.1093/nar/gkn741
-
Syndecan-1 regulates the biological activities of interleukin-34Biochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1853:1010–1021.https://doi.org/10.1016/j.bbamcr.2015.01.023
-
Eosinophils in innate immunity: an evolving storyCell and Tissue Research 343:57–83.https://doi.org/10.1007/s00441-010-1049-6
-
Lineage commitment: cytokines instruct, at last!Cell Stem Cell 5:234–236.https://doi.org/10.1016/j.stem.2009.08.015
-
Species variation in neutrophil biochemistry and functionJournal of Leukocyte Biology 46:63–74.https://doi.org/10.1002/jlb.46.1.63
-
Neutrophil kinetics in health and diseaseTrends in Immunology 31:318–324.https://doi.org/10.1016/j.it.2010.05.006
-
Support for lungfish as the closest relative of tetrapods by using slowly evolving Ray-finned fish as the outgroupGenome Biology and Evolution 9:evw288.https://doi.org/10.1093/gbe/evw288
-
IL-6 in inflammation, immunity, and diseaseCold Spring Harbor Perspectives in Biology 6:a016295.https://doi.org/10.1101/cshperspect.a016295
-
Effects of macrophage colony-stimulating factor on macrophages and their related cell populations in the osteopetrosis mouse defective in production of functional macrophage colony-stimulating factor proteinThe American Journal of Pathology 149:559–574.
-
Granulocyte margination in bone marrow: comparison with margination in the spleen and liverScandinavian Journal of Clinical and Laboratory Investigation 55:87–96.https://doi.org/10.3109/00365519509075382
-
Granulocytes are unresponsive to IL-6 due to an absence of gp130The Journal of Immunology 200:3547–3555.https://doi.org/10.4049/jimmunol.1701191
-
Granulopoiesis and neutrophil homeostasis: a metabolic, daily balancing actTrends in Immunology 40:598–612.https://doi.org/10.1016/j.it.2019.05.004
Decision letter
-
Patricia J WittkoppSenior Editor; University of Michigan, United States
-
Antonis RokasReviewing Editor; Vanderbilt University, United States
-
Leo KoendermanReviewer; UMC Utrecht, Netherlands
-
Clare PridansReviewer; University of Edinburgh, United Kingdom
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
In this study, Pinheiro and colleagues examine how neutrophils and other myeloid cells evolved. Using a wide range of relational taxonomic data, the authors show how pairings of genes in the granulocyte colony-stimulating factor receptor family (CSF3R/CSF3) and in the macrophage colony-stimulating factor receptor family (CSF1R/IL34/CSF1) evolved and contributed to granulocyte adaptations. This work sheds light into the evolution of mammalian neutrophils.
Decision letter after peer review:
Thank you for submitting your article "Analysis of receptor-ligand pairings and distribution of myeloid subpopulations across the animal kingdom reveals neutrophil evolution was facilitated by colony-stimulating factors" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Patricia Wittkopp as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Leo Koenderman (Reviewer #1); Clare Pridans (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). Specifically, we are asking editors to accept without delay manuscripts, like yours, that they judge can stand as eLife papers without additional data, even if they feel that they would make the manuscript stronger. Thus the revisions requested below only address clarity and presentation.
Summary:
Pinheiro and colleagues have described a fascinating view on the evolution of neutrophils and other myeloid cells. The authors used a wide range of relational taxonomic data to show how CSF1/CSF1R and CSF3/CSF3R pairings evolved and contributed to granulocyte adaptations. This is a very original and potentially important piece of work that sheds light into the evolution of mammalian neutrophils.
Essential revisions:
1) First of all Figure 4, is missing. This is an important figure as the text (as it is now) is difficult to understand.
2) At several locations in the article the authors imply that G-CSF is inducing differentiation fitting with an inductive model (e.g. Introduction). At the same time the authors rightly mention the presence of mature neutrophils in G-CSF-/- mice (as well as mature eosinophils in IL5R-/- mice) more pointing at a stochastic model. This latter model assumes that expression of CSF-R's is more random, and only committed progenitors expressing these receptors will proliferate rather than differentiate in response to these CSF. Please provide sufficient arguments for the inductive model or change part of the interpretations when a stochastic model is more likely.
3) In the whole article data are provided on numbers in peripheral blood. Only a minority of myeloid cells reside in the blood, the majority is in the tissues. The situation with neutrophils is uncertain. Please discuss.
4) The part on C-EBP transcription factors is difficult to follow. Please help the reader understand why they are so important (based on KO strategies) while there is no clear picture in evolution as the genes are sometimes present, sometimes not. Some species have many some only one. Simply stating redundancy in the system does not really fit the knock-out studies.
5) The part described in the subsection “The emergence of CSF3R/CSF3 and onset of endothermy likely influenced the distribution of neutrophils in Chordates during evolution”/Supplementary Figure 1 is not adding much to the article. It is only human with no evolutionary perspective. Consider removing.
6) Please provide some more insight into the issue of eosinophils versus neutrophils. Now it is implied that the co-evolution with endothermia is relevant. Many articles suggest that eosinophils are more specialized in killing large targets (extracellular killing/e.g. parasites) vs. neutrophils small targets (intracellular killing/e.g., bacteria). Can the authors provide their ideas about the functional difference of the cells in the evolutionary perspective.
7) Discussion. It is stated that neutrophils comprise the largest population of myeloid cells in mammals. This needs supportive evidence, as macrophages are thought to be the largest population at least in the tissues.
8) Discussion. Although the issue of the lamins is well taken formal proof that the segmented nuclear morphology of neutrophils is important for movement and trans-cellular migration is yet to be determined (e.g. van Grinsven et al., 2019 ).
9) Introduction. Young children with SCN often have mutations in the ELANE gene rather than the GSF-R gene. Can the authors discuss how ELANE fits with the model they are presenting?
10) Please provide the definitions of neutrophils and heterophils as they can be present as different cells in the same species.
11) Please provide a supplementary file containing all the references used for Figure 1H (complete blood count data; CBC). This would be a useful source of data for researchers interested in other blood cell types.
12) Regarding the CBC data – the authors should mention in the text if all the samples were obtained from adults. While we appreciate that n values are low for some species, do you obtain the same result if you analyse males and females separately? This may be worth mentioning given that neutrophil numbers have been reported to be higher in women.
13) Please provide a supplementary file containing all the NCBI gene and Ensembl accession numbers for each gene, in each species (Figure 2A).
14) The authors may want to mention that there are other receptors for IL-34 which may explain its expression (in fish, Figure 2A) in the absence of Csf1r.
15) Please provide a supplementary file containing all the NCBI protein accession numbers used for Figure 3A.
16) Please include isotype controls on histogram in Supplementary Figure 1A, C and D.
17) Please include full gating strategy for Supplementary Figure 1A.
18) Why was 72h chosen for the mobility assays (Supplementary Figure 1B)? At this point, monocytes cultured in CSF1 would begin differentiating into macrophages, and this may affect their mobility.
19) Supplementary Figure 1C – please include the antibodies in the Lin cocktail for flow cytometry in the figure legend.
20) Please mention in text and figure legend that human blood was used (there is no mention of it within text).
21) Was a dead cell exclusion dye used for flow cytometry of human blood and neutrophils? And did you look at FSC-A v FSC-H to exclude doublets? If not, how can you exclude the possibility that the Cxcr4 hi neutrophils are not dying or doublets?
22) Since this is a large selection of taxa groups, can specification (of a subset) be divided into more detail?
23) It would be helpful to check the sequence conservation of the receptors across these taxonomic families and see whether there are any minor evolution instances where they mutated. If the receptors have mutated, do they have a particular residue that mutated?
https://doi.org/10.7554/eLife.60214.sa1Author response
Essential revisions:
1) First of all Figure 4, is missing. This is an important figure as the text (as it is now) is difficult to understand.
We thank the reviewer for bringing this to our attention. This was a problem with transfer from BioXriv and eLife (pre-print copy had Figure 4 included). We have now resolved the issue with editorial team and Figure 4 is present in its entirety.
2) At several locations in the article the authors imply that G-CSF is inducing differentiation fitting with an inductive model (e.g. Introduction). At the same time the authors rightly mention the presence of mature neutrophils in G-CSF-/- mice (as well as mature eosinophils in IL5R-/- mice) more pointing at a stochastic model. This latter model assumes that expression of CSF-R's is more random, and only committed progenitors expressing these receptors will proliferate rather than differentiate in response to these CSF. Please provide sufficient arguments for the inductive model or change part of the interpretations when a stochastic model is more likely.
We thank the reviewers for raising these interesting points. As highlighted, functional redundancy is a prominent feature of myelopoiesis, as in the knock-out models identified, IL-6 and CSF2 can both fully or partially function as surrogates for either CSF3, IL5 or both [1, 2] (all references at end of this document). Although this is indicative that redundancy favours the stochastic model, we would argue in homeostatic conditions CSF3 favours an inductive model, as it is exclusively associated with neutrophil development and function. Competitive repopulation assays in CSF3R null mice have shown that CSF3 is responsible for almost all basal granulopoiesis including the regulation, production and maintenance of committed myeloid progenitors [3]. Similarly, in vitro CSF3 treatment of colony forming units -granulocyte, monocyte (CFU-GM), which are capable of producing monocyte or neutrophil populations, resulted in the uniform production of neutrophils, which is supportive of the inductive model [4, 5]. These studies also support a similar model where CSF3 is considered an instructive cytokine as opposed to permissive [6]. The bioavailability of CSF3R could be considered another indicator of the inductive model, as although it is also expressed on haematopoietic stem cells (HSCs), higher levels of expression are restricted to specific neutrophil progenitors increasing through to maturity [7, 8]. Furthermore, in basal conditions, CSF3 is present at higher concentrations in blood serum than CSF2, therefore the CSF3/CSF3R axis “effectively” outcompetes the CSF2/CSF2Ra/CSF2rb axis [8]. Within this context, it could be argued that the elevated presence of CSF3/CSF3R as the default or principle mechanism for neutrophil generation, in mammals at least, strongly supports an inductive model of development. We have now added a new sentence to reflect this hypothesis (Discussion).
3) In the whole article data are provided on numbers in peripheral blood. Only a minority of myeloid cells reside in the blood, the majority is in the tissues. The situation with neutrophils is uncertain. Please discuss.
We thank the reviewers for highlighting this important point. We agree the majority of myeloid cells are expected to reside within tissue. However, the focus of our work was to examine the blood microenvironment and evolution of bona fide blood myeloid cells. Correlating the evolution of CSFs with heterogeneity of myeloid subsets. We provocatively highlight that adaption of granulocyte like behaviour (migration between tight junctions) may correlate with arrival of lungs, blood vessels and acceleration of microbial infection diversity. More work will now be needed to define evolution of other cells outside the vasculature, such as tissue resident macrophages and bone fide subpopulations of granulocytes, such as neutrophils. We have now added text to reflect this in the Discussion.
4) The part on C-EBP transcription factors is difficult to follow. Please help the reader understand why they are so important (based on KO strategies) while there is no clear picture in evolution as the genes are sometimes present, sometimes not. Some species have many some only one. Simply stating redundancy in the system does not really fit the knock-out studies.
We thank the reviewer for their comments and agree the section needs further clarity. Three members of the C/EBP gene family; C/EBPε, C/EBPβ and C/EBPε have been shown to be required for neutrophil development. C/EBPβ, is the earliest expressed member during neutrophil development and is present on HSCs, common myeloid progenitors (CMPs) and granulocyte monocyte progenitors (GMPs) [9]. C/EBPβ is necessary for the transition from the CMP to GMP and when absent, granulopoiesis is impeded and does not progress beyond the CMP stage [9-13]. C/EBPβ is the next member to be expressed and is present on GMPs through to maturation and is required for emergency and/or cytokine induced granulopoiesis. C/EBPβ-deficient mice failed to mobilise granulocytes in response to systemic fungal infection or cytokine stimuli [9, 14-17]. Interestingly, although C/EBPβ is not required for steady-state haematopoiesis, it has been shown in models to functionally compensate for the loss of C/EBPβ [18]. Finally, C/EBPε is required for the terminal differentiation of granulocytes and is the last of the three to be expressed and is only present on neutrophil-committed promyelocytes until maturation [19-24]. We have now expanded this section of work by including a table (new Supplementary file 2) that further details the C/EBP-related knock-out studies to help the reader understand the relevance of these transcription factors to neutrophil development. We agree that system redundancy is not necessarily an appropriate measure and have re-written the section to reflect that the C/EBP gene family is very strongly evolutionarily conserved across Phylum Chordata, indicative of their importance to neutrophil development (subsection “Interrogation of CSF1/CSF1R and CSF3/CSF3R and C/EBP gene family reveals their loss in the early lineages”).
5) The part described in the subsection “The emergence of CSF3R/CSF3 and onset of endothermy likely influenced the distribution of neutrophils in Chordates during evolution”/Supplementary Figure 1 is not adding much to the article. It is only human with no evolutionary perspective. Consider removing.
We thank the reviewers for their comments and agree that as the work currently stands it does not add much to the in silico work presented. We have decided to remove this work, as more lab work is required to complete the preliminary studies discussed and is out of scope of current in silico studies. We are now seeking grant funding to complete this work.
6) Please provide some more insight into the issue of eosinophils versus neutrophils. Now it is implied that the co-evolution with endothermia is relevant. Many articles suggest that eosinophils are more specialized in killing large targets (extracellular killing/e.g. parasites) vs. neutrophils small targets (intracellular killing/e.g., bacteria). Can the authors provide their ideas about the functional difference of the cells in the evolutionary perspective.
We thank the reviewers for their comments. We have considered how eosinophils and neutrophils are specialized towards intra- and extracellular killing respectively Firstly, neutrophils have an arsenal of primary and secondary granules such as; defensins, cathelicidins and lysozymes, which have known antimicrobial properties [25]. In contrast, although eosinophils have a similar set up of primary and secondary granules, eosinophil cationic protein (ECP) and major basic protein (MBP) are both key for extracellular killing, as they are toxic to large targets such as helminth parasites and hemoflagellates as well bacteria [26]. While both eosinophils and neutrophils use respiratory bursts as a killing function, the mechanisms differ. Eosinophils are known to have more O2- species than neutrophils, and this is believed to be required for the extracellular targeting of the respiratory burst, a feature unique to eosinophilic granulocytes and implicated in host-pathogen defences against parasites [27-29]. In contrast, neutrophils favour intracellular respiratory bursts for the destruction of internalised pathogens [27-29]. Interestingly, although eosinophils can phagocytose extracellular S. aureus and E. coli, they are unable to kill them as efficiently as neutrophils, who are by far the more efficient cell type [30].
Our analysis suggests that as early as the emergence of amphibians, neutrophils and eosinophils had likely acquired their respective functional properties, however, the population proportions are relatively similar. Therefore, it could be posited that the re-organisation of blood granulocyte populations is linked to the environment and likely corresponds with the type and prevalence of pathogens, i.e. pathogen exposure shapes the host myeloid cell pool and selects for the most appropriate cell type. In some ectothermic lineages, where presumably parasites and viral or bacterial pathogens are more prevalent, there is a near equal distribution of eosinophils and neutrophils. However, it could be hypothesised that following the transition from ectothermy to endothermy, there would have been changes to existing pathogens such as loss of pathogenicity or function as they failed to cross over and infect the newly emerging endothermic species, as evidenced by the inability of some viruses and fungi to infect modern day mammals [31-33]. Ultimately, this “thermic switch” would lead to the emergence of novel pathogenic strains that could survive in endothermic hosts (or both). As neutrophils and heterophils dominate the blood of land-dwelling mammals and birds respectively, this again supports an argument that in response to external selection pressures, the host species are redistributing their myeloid blood cell pool and selecting for neutrophil predominance. We have commented further on this within the manuscript (subsection “The emergence of CSF3R/CSF3 and onset of endothermy likely influenced the distribution of neutrophils in Chordates during evolution”).
7) Discussion. It is stated that neutrophils comprise the largest population of myeloid cells in mammals. This needs supportive evidence, as macrophages are thought to be the largest population at least in the tissues.
We thank the reviewer for their comment. We have adjusted our statement in the manuscript to reflect that it applied to the blood microenvironment only (Discussion), where neutrophils are the largest population of myeloid cells. Also see point 3 above.
8) Discussion. Although the issue of the lamins is well taken formal proof that the segmented nuclear morphology of neutrophils is important for movement and trans-cellular migration is yet to be determined (e.g. van Grinsven et al., 2019 ).
We thank the reviewers for bringing this important and highly relevant article to our attention. It is now evident that immature neutrophils (band cells) can efficiently migrate with incomplete segmentation, suggesting at least, that it is not always a requirement for movement and trans-cellular migration in all neutrophil populations. We have amended the text of the manuscript to reflect this (Discussion).
9) Introduction. Young children with SCN often have mutations in the ELANE gene rather than the GSF-R gene. Can the authors discuss how ELANE fits with the model they are presenting?
We thank the reviewers for their interesting comments, and we have further considered the role of ELANE in the model. The ELANE gene encodes for the neutrophil granule serine protease- neutrophil elastase- and is a negative regulator of CSF3R [34-37]. in vitro studies have shown the ELANE can downregulate CSF3R expression on the surface of mature neutrophils and can abrogate proliferative signals provided by CSF3 for expansion of the relevant committed myeloid progenitors [34, 35]. Similarly, it has also been shown that ELANE can antagonise the in vitro activity of CSF3 [34, 35]. Therefore, ELANE is heavily implicated in the regulation of neutrophil development via CSF3R/CSF3 signalling. During SCN the arrest of neutrophil development at the promyelocyte stage and concomitant apoptosis of those precursors results in acute neutropenia. There are several theories as to how aberrant ELANE expression could be implicated in this. Some studies have shown that improper ELANE expression results in endoplasmic reticulum stress (ER stress) which can lead to subcellular mislocalization and subsequent cell death of progenitors [36, 37]. Alternatively, other groups posit that disrupted or truncated CSF3R/CSF3 signalling as mediated by the aberrant expression of mutated ELANE could be a possible mechanism for the subsequent loss of neutrophils [34, 35]. A recent study that used an inducible expression system in mice, showed that a specific ELANE mutation diminished enzymatic activity which resulted in impaired granulocytic differentiation and reduced expression of key neutrophil associated genes such as GFi1, C/EBPδ, C/EBPε and CSF3R [37]. Taken together, these studies suggest there could be multiple mechanisms through which ELANE can dysregulate granulopoiesis and is likely mutation specific. ELANE fits into our model as the disruption caused be improperly functioning ELANE to granulopoiesis is significant as it directly impacts on CSF3R/CSF3 signalling and through this has downstream consequences for the development of neutrophils from committed progenitors. We have added a line to reflect that mutations in the ELANE gene are also involved in SCN pathogenesis (Introduction).
10) Please provide the definitions of neutrophils and heterophils as they can be present as different cells in the same species.
We thank the reviewer for their comments. Heterophils are described as a phagocytic leukocyte that is functionally analogous to the mammalian neutrophil. Heterophils can be distinguished from neutrophils using Romanowsky staining techniques as they do not stain neutral because of the presence of red -orange eosinophilic cytoplasmic granules [38, 39] A definition has been added to manuscript (Introduction).
11) Please provide a supplementary file containing all the references used for Figure 1B (complete blood count data; CBC). This would be a useful source of data for researchers interested in other blood cell types.
We thank the reviewer for their comments. A supplementary file (Supplementary file 3) has been added that contains all the references included in the meta-analysis.
12) Regarding the CBC data – the authors should mention in the text if all the samples were obtained from adults. While we appreciate that n values are low for some species, do you obtain the same result if you analyse males and females separately? This may be worth mentioning given that neutrophil numbers have been reported to be higher in women.
We thank the reviewer for their comments and have clarified within the Materials and methods section, that datasets used for the metanalysis used sub-adult and adult data and were pooled with a minimum of male and female set where the data was available. It was not possible to take this approach for every species considered and within the context of this current study it is not possible to comment definitively on if there are gender-specific differences observed.
13) Please provide a supplementary file containing all the NCBI gene and Ensembl accession numbers for each gene, in each species (Figure 2A).
We thank the reviewer for their comments. A supplementary file (Supplementary file 4) has been added that contains all the ID numbers for each gene, in each species examined.
14) The authors may want to mention that there are other receptors for IL-34 which may explain its expression (in fish, Figure 2A) in the absence of Csf1r.
We thank for the reviewer for bringing these receptors to our attention. We have now included a line referencing the IL-34 receptors, Protein tyrosine phosphatase zeta and CD138 (Sydecan), but have cited this in the general discussion regarding the emergence and co-evolution of CSF1R and IL34 (subsection “Analysis of Chordate orthologous protein homology further supports the ancestral pairing of CSF1R/IL34 and CSF3R/CSF3 in early lineages”).
15) Please provide a supplementary file containing all the NCBI protein accession numbers used for Figure 3A.
We thank the reviewer for their comments. A supplementary file (Supplementary file 5) has been added that contains all the ID numbers for each protein, in each species examined.
16) Please include isotype controls on histogram in Supplementary Figure 1A, C and D.
We thank for the reviewer for their comments, this section of work has now been removed.
17) Please include full gating strategy for Supplementary Figure 1A.
We thank for the reviewer for their comments, this section of work has now been removed.
18) Why was 72h chosen for the mobility assays (Supplementary Figure 1B)? At this point, monocytes cultured in CSF1 would begin differentiating into macrophages, and this may affect their mobility.
We thank for the reviewer for their comments, this section of work has now been removed.
19) Supplementary Figure 1C – please include the antibodies in the Lin cocktail for flow cytometry in the figure legend.
We thank for the reviewer for their comments, this section of work has now been removed.
20) Please mention in text and figure legend that human blood was used (there is no mention of it within text).
We thank for the reviewer for their comments, this section of work has now been removed.
21) Was a dead cell exclusion dye used for flow cytometry of human blood and neutrophils? And did you look at FSC-A v FSC-H to exclude doublets? If not, how can you exclude the possibility that the Cxcr4 hi neutrophils are not dying or doublets?
We thank for the reviewer for their comments, this section of work has now been removed.
22) Since this is a large selection of taxa groups, can specification (of a subset) be divided into more detail?
We thank the reviewers for their comments. Where applicable larger taxonomic groups have been subdivided into smaller groups such as reptilia being divided in squamata, crocodilia and testudines. Similarly, mammalia has also been subdivided into the respective groups of monotremata, marsupalia and placentalia. For the purpose of these studies, the cell population of interest does not further subset as determined by CSF3R.
23) It would be helpful to check the sequence conservation of the receptors across these taxonomic families and see whether there are any minor evolution instances where they mutated. If the receptors have mutated, do they have a particular residue that mutated?
We thank the reviewers for this good suggestion. However, we feel that this is beyond the scope of the original focus of the study. The work here addresses changes on a species or sub-species wide scale and how this might have contributed overall to the functional properties of mammalian neutrophils. Although the suggested idea is very interesting, we feel that focussing on changes at the residue level would not really fit in with the wider macro level view and might confuse the overall message.
References
1) Walker, F., et al., IL6/sIL6R complex contributes to emergency granulopoietic responses in G-CSF– and GM-CSF–deficient mice. Blood, 2008. 111(8): p. 3978-3985.2)Kopf, M., et al., IL-5-deficient mice have a developmental defect in CD5+ B-1 cells and lack eosinophilia but have normal antibody and cytotoxic T cell responses. Immunity, 1996. 4(1): p. 15-24.3) Richards, M.K., et al., Pivotal role of granulocyte colony-stimulating factor in the development of progenitors in the common myeloid pathway. Blood, 2003. 102(10): p. 3562-3568.4) Sonoda, Y., et al., Analysis in serum-free culture of the targets of recombinant human hemopoietic growth factors: interleukin 3 and granulocyte/macrophage-colony-stimulating factor are specific for early developmental stages. Proceedings of the National Academy of Sciences, 1988. 85(12): p. 4360.5) Avalos, B.R., Molecular analysis of the granulocyte colony-stimulating factor receptor. Blood, 1996. 88(3): p. 761-777.6) Stanley, E.R., Lineage Commitment: Cytokines Instruct, At Last! Cell Stem Cell, 2009. 5(3): p. 234-236.7) Demetri, G.D. and J.D. Griffin, Granulocyte colony-stimulating factor and its receptor. Blood, 1991. 78(11): p. 2791-808.8) Lee, K.Y., et al., Varying expression levels of colony stimulating factor receptors in disease states and different leukocytes. Experimental and Molecular Medicine, 2000. 32(4): p. 210-215.9) Scott, L.M., et al., A novel temporal expression pattern of three C/EBP family members in differentiating myelomonocytic cells. Blood, 1992. 80(7): p. 1725-35.10) Smith, L.T., et al., PU.1 (Spi-1) and C/EBP α regulate the granulocyte colony-stimulating factor receptor promoter in myeloid cells. Blood, 1996. 88(4): p. 1234-47.11) Ma, O., et al., Granulopoiesis requires increased C/EBPα compared to monopoiesis, correlated with elevated Cebpa in immature G-CSF receptor versus M-CSF receptor expressing cells. PLoS One, 2014. 9(4): p. e95784.12) Zhang, P., et al., Enhancement of hematopoietic stem cell repopulating capacity and self-renewal in the absence of the transcription factor C/EBP α. Immunity, 2004. 21(6): p. 853-63.13) Ford, A.M., et al., Regulation of the myeloperoxidase enhancer binding proteins Pu1, C-EBP α, -β, and -δ during granulocyte-lineage specification. Proc Natl Acad Sci U S A, 1996. 93(20): p. 10838-43.14) Zhang, P., et al., Induction of granulocytic differentiation by 2 pathways. Blood, 2002. 99(12): p. 4406-12.15) Hirai, H., et al., C/EBPbeta is required for 'emergency' granulopoiesis. Nat Immunol, 2006. 7(7): p. 732-9.16) Screpanti, I., et al., Lymphoproliferative disorder and imbalanced T-helper response in C/EBP β-deficient mice. EMBO J, 1995. 14(9): p. 1932-41.17) Tanaka, T., et al., Targeted disruption of the NF-IL6 gene discloses its essential role in bacteria killing and tumor cytotoxicity by macrophages. Cell, 1995. 80(2): p. 353-61.18) Jones, L.C., et al., Expression of C/EBPbeta from the C/ebpalpha gene locus is sufficient for normal hematopoiesis in vivo. Blood, 2002. 99(6): p. 2032-6.19) Yamanaka, R., et al., Impaired granulopoiesis, myelodysplasia, and early lethality in CCAAT/enhancer binding protein epsilon-deficient mice. Proc Natl Acad Sci U S A, 1997. 94(24): p. 13187-92.20) Morosetti, R., et al., A novel, myeloid transcription factor, C/EBP epsilon, is upregulated during granulocytic, but not monocytic, differentiation. Blood, 1997. 90(7): p. 2591-600.21) Verbeek, W., et al., C/EBPepsilon -/- mice: increased rate of myeloid proliferation and apoptosis. Leukemia, 2001. 15(1): p. 103-11.22) Chih, D.Y., et al., Modulation of mRNA expression of a novel human myeloid-selective CCAAT/enhancer binding protein gene (C/EBP epsilon). Blood, 1997. 90(8): p. 2987-94.23) Lekstrom-Himes, J. and K.G. Xanthopoulos, CCAAT/enhancer binding protein epsilon is critical for effective neutrophil-mediated response to inflammatory challenge. Blood, 1999. 93(9): p. 3096-105.24) Gombart, A.F., et al., Aberrant expression of neutrophil and macrophage-related genes in a murine model for human neutrophil-specific granule deficiency. J Leukoc Biol, 2005. 78(5): p. 1153-65.25) Borregaard, N. and J.B. Cowland, Granules of the Human Neutrophilic Polymorphonuclear Leukocyte. Blood, 1997. 89(10): p. 3503-3521.26) Shamri, R., J.J. Xenakis, and L.A. Spencer, Eosinophils in innate immunity: an evolving story. Cell and tissue research, 2011. 343(1): p. 57-83.27) Bolscher, B.G.J.M., et al., NADPH:O2 oxidoreductase of human eosinophils in the cell-free system. FEBS Letters, 1990. 268(1): p. 269-273.28) Lacy, P., et al., Divergence of Mechanisms Regulating Respiratory Burst in Blood and Sputum Eosinophils and Neutrophils from Atopic Subjects. The Journal of Immunology, 2003. 170(5): p. 2670.29) Kovács, I., et al., Comparison of proton channel, phagocyte oxidase, and respiratory burst levels between human eosinophil and neutrophil granulocytes. Free Radical Research, 2014. 48(10): p. 1190-1199.30) Hatano, Y., et al., Phagocytosis of heat-killed Staphylococcus aureus by eosinophils: comparison with neutrophils. APMIS, 2009. 117(2): p. 115-23.31) Robert, V.A. and A. Casadevall, Vertebrate Endothermy Restricts Most Fungi as Potential Pathogens. The Journal of Infectious Diseases, 2009. 200(10): p. 1623-1626.32) Chinchar, V.G., Ranaviruses (family Iridoviridae): emerging cold-blooded killers. Archives of Virology, 2002. 147(3): p. 447-470.33) Wirth, W., et al., Ranaviruses and reptiles. PeerJ, 2018. 6: p. e6083-e6083.34) Hunter, M.G., et al., Proteolytic cleavage of granulocyte colony-stimulating factor and its receptor by neutrophil elastase induces growth inhibition and decreased cell surface expression of the granulocyte colony-stimulating factor receptor. Am J Hematol, 2003. 74(3): p. 149-55.35) Piper, M.G., et al., Neutrophil elastase downmodulates native G-CSFR expression and granulocyte-macrophage colony formation. J Inflamm (Lond), 2010. 7(1): p. 5.36) Nayak, R.C., et al., Pathogenesis of ELANE-mutant severe neutropenia revealed by induced pluripotent stem cells. J Clin Invest, 2015. 125(8): p. 3103-16.37) Garg, B., et al., Inducible expression of a disease-associated. J Biol Chem, 2020. 295(21): p. 7492-7500.38) Montali, R.J., Comparative pathology of inflammation in the higher vertebrates (reptiles, birds and mammals). Journal of Comparative Pathology, 1988. 99(1): p. 1-26.39) Horobin, R.W., How Romanowsky stains work and why they remain valuable — including a proposed universal Romanowsky staining mechanism and a rational troubleshooting scheme. Biotechnic and Histochemistry, 2011. 86(1): p. 36-51.
https://doi.org/10.7554/eLife.60214.sa2Article and author information
Author details
Funding
Medical Research Council (MR/M003159/1)
- Kevin J Woollard
Kidney Research UK (RP_019_20160303)
- Kevin J Woollard
Kidney Research UK (RP_002_20170914)
- Kevin J Woollard
British Heart Foundation (PG/18/41/33813)
- Kevin J Woollard
The authors declare that there was no direct funding for this work. Grants from MRC (MR/M003159/1), Kidney Research UK (RP_019_20160303, RP_002_20170914) and BHF (PG/18/41/33813) support the Woollard lab. KJW is now an employee for AstraZeneca (BioPharmaceuticals R&D, Cambridge, UK). All of this work was performed at Imperial College London. No funding or support was received from AstraZeneca.
Senior Editor
- Patricia J Wittkopp, University of Michigan, United States
Reviewing Editor
- Antonis Rokas, Vanderbilt University, United States
Reviewers
- Leo Koenderman, UMC Utrecht, Netherlands
- Clare Pridans, University of Edinburgh, United Kingdom
Publication history
- Received: June 21, 2020
- Accepted: November 24, 2020
- Accepted Manuscript published: November 25, 2020 (version 1)
- Version of Record published: December 4, 2020 (version 2)
Copyright
© 2020, Pinheiro et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,270
- Page views
-
- 128
- Downloads
-
- 5
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Ecology
- Evolutionary Biology
Most phytophagous insect species exhibit a limited diet breadth and specialize on a few or a single host plant. In contrast, some species display a remarkably large diet breadth, with host plants spanning several families and many species. It is unclear, however, whether this phylogenetic generalism is supported by a generic metabolic use of common host chemical compounds (‘metabolic generalism’) or alternatively by distinct uses of diet-specific compounds (‘multi-host metabolic specialism’)? Here, we simultaneously investigated the metabolomes of fruit diets and of individuals of a generalist phytophagous species, Drosophila suzukii, that developed on them. The direct comparison of metabolomes of diets and consumers enabled us to disentangle the metabolic fate of common and rarer dietary compounds. We showed that the consumption of biochemically dissimilar diets resulted in a canalized, generic response from generalist individuals, consistent with the metabolic generalism hypothesis. We also showed that many diet-specific metabolites, such as those related to the particular color, odor, or taste of diets, were not metabolized, and rather accumulated in consumer individuals, even when probably detrimental to fitness. As a result, while individuals were mostly similar across diets, the detection of their particular diet was straightforward. Our study thus supports the view that dietary generalism may emerge from a passive, opportunistic use of various resources, contrary to more widespread views of an active role of adaptation in this process. Such a passive stance towards dietary chemicals, probably costly in the short term, might favor the later evolution of new diet specializations.
-
- Evolutionary Biology
- Genetics and Genomics
The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.