Architecture of the neutrophil compartment in mice [scRNA-seq clonal]

In this study, we present NeuMAP, a comprehensive single-cell analysis of neutrophils spanning over 40 anatomical, physiological, and pathological contexts in mice. NeuMAP confirms and expands previous models of neutrophil diversity, revealing the organization of neutrophils into distinct functional hubs under both normal and pathological conditions. Furthermore, we delineate prototypical trajectories as key organizers of granulopoiesis and examine neutrophil dynamics along these trajectories during acute inflammation and cancer. Integrating insights from fate mapping, mutant mouse models, and in vitro experiments, we identified immunological signals guiding neutrophil through these typical trajectories of granulopoiesis, that balance immune protection, tissue homeostasis and repair. Specifically, we identified IFNB, GMCSF, and TGFB as drivers of pro-inflammatory, cancer-associated, and mature neutrophil states, respectively and found that the transcription factor JUNB drives the angiogenic and immunosuppressive function of neutrophils. Additionally, we uncover conserved transcriptional signatures for human neutrophil states, validate their prognostic significance in cancer patients, and introduce a proof-of-concept strategy for exploring the diagnostic potential of neutrophil functional diversity across various diseases using single-cell transcriptomics of blood neutrophils. Overall design: Model of Clonal Hematopoiesis and PCSK9-Induced Hypercholesterolemia: To model TET2 loss-of-function-driven clonal hematopoiesis, we performed an adoptive bone marrow transfer (ADT) without pre-conditioning. Ten-week-old unirradiated JAXBoy (PtprcK302E) recipient mice were intravenously injected with a total of 1.5 × 107 unfractionated CD45.2? Tet2–/– bone marrow (BM) cells, administered as three consecutive daily doses of 5 × 106 cells. Donor cells were harvested from age-matched littermate Tet2–/– mice (8–10 weeks old) by flushing femurs and tibias following euthanasia. To induce hypercholesterolemia, a recombinant AAV vector encoding a gain-of-function form of murine or PCSK9 (pAAV/D377Y-mPCSK9) was delivered via a single tail vein injection. One week later, mice were placed on either a high-cholesterol Western diet (Envigo, TD.88137; 42% calories from fat, 0.2% cholesterol) or a matched control diet for 13 weeks. At endpoint, ADT mice were euthanized and Tet2–/– (CD45.2?) or WT (CD45.1?) neutrophils were isolated from peripheral blood, lung and bone marrow by cell sorting. Live CD45.1+CD11b+-Ly6G+ and CD45.2+CD11b+-Ly6G+ cells were sorted out. After sorting, the neutrophil suspensions were Taged (BD Multiplexing Kit) and pooled in equal proportion. For the generation of single-cell whole-transcriptomes we used a BD Rhapsody system according to the manufacturer instructions. Briefly, cell suspensions from each condition were incubated with Sample Tags (BD) for 20 minutes at RT. Cells were then washed three times and pooled in a single tube. Viability and cell concentration using the Countess III cell counter (Thermofisher). 60,000 cells were loaded into a Rhapsody Single Cell Analysis System cartridge. Cell capture and cDNA synthesis were performed according to manufacturer's instructions; cells were isolated into nanowells by gravity, then cells lysed and mRNAs together with sample tags oligos were released and captured by beads present in the nanowells. Each bead contained a unique oligo named “cell label” to identify each individual bead. All beads present in the cartridge were collected and cDNA synthesis took place in a single reaction. At this point each cDNA and Sample Tag oligo were attached to its corresponding “cell label” oligo. Two separated indexed libraries were prepared for whole transcriptome analysis and Sample Tag demultiplexing following manufacturer's instructions. The average size of the libraries was calculated using the 2100 Bioanalyzer (Agilent) and the concentration was determined using the Qubit fluorometer (Thermofisher). Finally, libraries were combined and sequenced together in a paired-end run (60x42) using a NextSeq 2000 system (Illumina) and a P2 flow cell. Output files were processed with NextSeq 1000/2000 Control Software Suite v1.4.1. FastQ files for each sample were obtained using BCL Convert v3.6.3 software (Illumina).