RESUMEN
Recent studies using advanced techniques such as single cell RNA sequencing (scRNAseq), high parameter flow cytometry, and proteomics reveal that neutrophils are more heterogeneous than previously appreciated. Unique subsets have been identified in the context of bacterial and parasitic infections, cancer, and tissue injury and repair. The characteristics of infiltrating neutrophils differ depending on the nature of the inflammation-inciting stimulus, the stage of the inflammatory response, as well as the tissue microenvironment in which they accumulate. We previously described a new subpopulation of immature Ly6Glow neutrophils that accumulate in the peritoneal cavity 3 days following intraperitoneal (i.p.) administration of the fungal cell wall extract, zymosan. These neutrophils express markers of alternative activation and possess neuroprotective/regenerative properties. In addition to inducing neurite outgrowth of explanted neurons, they enhance neuronal survival and axon regeneration in vivo following traumatic injury to the optic nerve or spinal cord. In contrast, the majority of neutrophils that accumulate in the peritoneal fluid 4 hours following i.p. zymosan injection (4h NΦ) have features of conventional, mature Ly6Ghi neutrophils and lack neuroprotective or neuroregenerative properties. In the current study, we expand upon on our previously published observations by performing a granular, in-depth analysis of these i.p. zymosan-modulated neutrophil populations using scRNAseq and high parameter flow cytometry. We also analyze cell lysates of each neutrophil population by liquid chromatography/mass spectrometry. Circulating blood neutrophils, harvested from naive mice, are analyzed in parallel as a control. When samples were pooled from all three groups, scRNAseq revealed 11 distinct neutrophil clusters. Pathway analyses demonstrated that 3d NΦ upregulate genes involved in tissue development and wound healing, while 4h NΦ upregulate genes involved in cytokine production and perpetuation of the immune response. Proteomics analysis revealed that 3d NΦ and 4h NΦ also express distinct protein signatures. Adding to our earlier findings, 3d NΦ expressed a number of neuroprotective/neuroregenerative candidate proteins that may contribute to their biological functions. Collectively, the data generated by the current study add to the growing literature on neutrophil heterogeneity and functional sub-specialization and might provide new insights in elucidating the mechanisms of action of pro-regenerative, neuroprotective neutrophil subsets.