
==== Front
101688624
Sci Immunol
Sci Immunol
Science immunology
2470-9468

39093958
10.1126/sciimmunol.ado1227
ems198347
Article
Recruited atypical Ly6G+ macrophages license alveolar regeneration after lung injury
Ruscitti C. 12
Abinet J. 12
Maréchal P. 12
Meunier M. 12
de meeûs C. 23
Vanneste D. 12
Janssen P. 12
Dourcy M. 24
Thiry M. 5
Bureau F. 26
Schneider C. 7
Machiels B. 24
Hidalgo A. 89
Ginhoux F. 10111213
Dewals B.G. 24
Guiot J. 1415
Schleich F. 1415
Garigliany M-M. 23
Bellahcène A. 16
Radermecker C. 12*
Marichal T. 1217*
1 Laboratory of Immunophysiology, GIGA Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
2 Faculty of Veterinary Medicine, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
3 Department of Pathology, FARAH Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
4 Laboratory of Immunology-Vaccinology, FARAH Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
5 Laboratory of Cellular and Tissular Biology, GIGA Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
6 Laboratory of Cellular and Molecular Immunology, GIGA Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
7 Institute of Physiology, https://ror.org/02crff812 University of Zurich ; Zurich, Switzerland
8 Area of Cell & Developmental Biology, Centro Nacional de Investigaciones Cardiovasculares Carlos III; Madrid, Spain
9 Vascular Biology and Therapeutics Program and Department of Immunobiology, Yale University School of Medicine; New Haven, CT, USA
10 Shanghai Institute of Immunology, Shanghai JiaoTong University School of Medicine; Shanghai, China
11 Inserm U1015, Gustave Roussy, Bâtiment de Médecine Moléculaire ; Villejuif, France
12 https://ror.org/03vmmgg57 Singapore Immunology Network (SIgN) , https://ror.org/036wvzt09 Agency for Science, Technology and Research (A*STAR) ; Singapore, Singapore
13 Translational Immunology Institute, https://ror.org/00xcwps97 SingHealth Duke-NUS Academic Medical Centre ; Singapore, Singapore
14 Laboratory of Pneumology, GIGA Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
15 Department of Respiratory Medicine, CHU University Hospital; Liège, Belgium
16 Metastasis Research Laboratory, GIGA Institute, https://ror.org/00afp2z80 University of Liège ; Liège, Belgium
17 https://ror.org/04qbvw321 Walloon Excellence in Life Sciences and Biotechnology (WELBIO) Department, WEL Research Institute; Wavre, Belgium
* Corresponding authors: c.radermecker@uliege.be; t.marichal@uliege.be
* These authors contributed equally to this work and are co-last authors

02 8 2024
02 8 2024
30 8 2024
07 9 2024
9 98 eado1227eado1227
https://creativecommons.org/licenses/by/4.0/ This work is licensed under a BY 4.0 International license.
The lung is constantly exposed to airborne pathogens and particles that can cause alveolar damage. Hence, appropriate repair responses are essential for gas exchanges and life. Here, we deciphered the spatiotemporal trajectory and function of an atypical population of macrophages after lung injury. Post-influenza A virus (IAV) infection, short-lived monocyte-derived Ly6G-expressing macrophages (Ly6G+ Macs) were recruited to the alveoli of lung perilesional areas. Ly6G+ Macs engulfed immune cells, exhibited a high metabolic potential and clustered with alveolar type 2 (AT2) epithelial cells in zones of active epithelial regeneration. Ly6G+ Macs were partially dependent on GM-CSF and IL-4 receptor signaling and were essential for AT2-dependent alveolar regeneration. Similar macrophages were recruited in other models of injury and in the airspaces of lungs from patients suspected of pneumonia. This study identifies perilesional alveolar Ly6G+ Macs as a spatially-restricted, short-lived macrophage subset promoting epithelial regeneration post-injury, thus representing an attractive therapeutic target for treating lung damage.
==== Body
pmcIntroduction

Severe respiratory viral infections represent a global health issue and a major threat for the healthcare systems as they often require hospitalization, such as seen during annual influenza A virus (IAV) outbreaks or the Covid-19 pandemic. Acute lung infectious episodes are typically associated with excessive lung inflammation, damage, and abnormal tissue repair that can lead to acute respiratory distress syndrome (ARDS), pneumonia and death (1–5). Deciphering the mechanisms eliciting appropriate lung regeneration and host recovery after viral-triggered lung injury is urgently needed to improve clinical management and broaden therapeutic opportunities.

Blood monocytes are heavily recruited to the lungs during the acute inflammatory phase post-infection, thereby contributing to host innate defense mechanisms. When not appropriately regulated, they are also thought to contribute to uncontrolled inflammation via the aberrant release of cytokines, which, in its extreme form, is known as the “cytokine storm”, well described in severe Covid-19 patients (3, 6–9). Recruited monocytes can also differentiate into monocyte-derived macrophages (Mo-Macs) that are either short-lived or can establish long-term residency in particular niches and can have functional consequences for lung immunity (10–15). The idea that recruited Mo-Macs are considered pathogenic after influenza infection (15, 16) or Covid-19 (6, 17), while lung-resident alveolar macrophages (AM) and interstitial macrophages (IM) exert beneficial roles (13, 18–24) is likely oversimplistic. Indeed, recruited Mo-Macs are increasingly recognized as heterogenous and can adopt distinct functional identities that depend on their differentiation trajectory, the diseased tissue microenvironment, the extent and phase of inflammation, and their activation state (9, 13, 14, 18, 25). In this regard, the fate and functions of short-lived Mo-Macs after lung injury remain incompletely described.

Here, we used an in vivo model of IAV-triggered injury to investigate the spatiotemporal trajectory and function of an atypical population of Mo-Macs expressing Ly6G, a marker considered to be restricted to granulocytes. We found that Ly6G+ Macs emerged transiently from Ccr2-dependent monocytes during the early recovery phase post-IAV, populated the alveolar lumen of lung perilesional areas and could promote progenitor AT2 differentiation and alveolar re-epithelization. Our study thus unravels the fate and function of a previously undescribed Ly6G+ Mac population that engages in crosstalk with epithelial cells to promote epithelial repair after viral-triggered lung injury.

Results

Lung Ly6G+ Macs emerge in the early recovery phase post-IAV infection

To investigate the dynamic of myeloid cell responses after lung injury, we performed time-course flow cytometry studies in a clinically relevant mouse model of lung infectious injury following IAV infection (2). Eight to twelve weeks-old C57BL/6 wild-type (WT) mice were infected intranasally (i.n.) with 5 plaque-forming units (PFU) of IAV H1N1 strain PR8/34, which triggered self-limiting disease with a peak in viral RNA at day 5 post-IAV and viral clearance at day 10 post-IAV (fig. S1, A and B). In this model, Ly6G+CD11b+CD64− neutrophils (Neu) increased at day 5 post-IAV and returned to baseline by day 15 (Fig. 1, A and B). A partial loss of Ly6G−CD64+SiglecF+CD11c+ AM was observed between day 5 and day 10 post-IAV, which was restored at day 15 post-IAV, as described (15) (Fig. 1, A and B). The numbers of CD64−Ly6C− monocytes (Ly6C− Mo) remained stable over the course of infection, unlike those of classical CD64−Ly6C+ Mo (Ly6C+ Mo) and inflammatory CD64+Ly6C+ monocytes (iMo) that peaked around day 5 post-IAV (Fig. 1, A and B). Macrophages (Macs) resembling IM (IM-like), defined as F4/80+CD11b+Ly6G−SiglecF−Ly6C−CD64+ cells and likely encompassing resident IM and recruited Mo-Macs, increased over time (Fig. 1, A and B). We also observed emergence, from day 5 onwards, of a distinct population of IAV-triggered Ly6G+CD11b+CD64+ Macs that fell in the classical Ly6G+CD11b+ neutrophil gate but were clearly distinct from neutrophils based on their elevated CD64 expression (Fig. 1, C and D), which we call Ly6G+ Macs hereafter. Ly6G+ Macs were absent in the blood (Fig. 1E), peaked at day 10 post-IAV in the lung and could still be detected at days 15 and 20 post-IAV (Fig. 1D). Morphologically, Ly6G+ Macs analyzed at day 10 post-IAV exhibited a kidney-shaped nucleus, like iMo, and possessed numerous cytoplasmic vacuolated structures and a cell membrane rich in protrusions (Fig. 1F and fig. S1C). Phenotypically, Ly6G+ Macs were F4/80hiSiglecFloCD11clo (fig. S1, D and E) and expressed high levels of the chemokine receptor CXCR4, of type II major histocompatibility complex (MHC-II), of CD101 and of CD319, a regulator of Mac functions (26, 27) (Fig. 1, G and H). However, Ly6G+ Macs exhibited low expression of the neutrophil activation marker CD177 (Fig. 1, G and H).

Next, we performed single cell (sc) RNA-sequencing (scRNA-seq) analyses of lung myeloid cells at day 10 post-IAV. Lung CD45+F4/80+ and/or CD11b+ cells were sorted from five mock-infected and 5 IAV-infected mice and were subjected to sc droplet encapsulation (28), scRNA-seq and quality control filtering. The curated data were integrated with a published dataset of steady-state lung monocytes and IM (29) and projected to global and condition-specific uniform manifold approximation and projection (UMAP) plots (Fig. 2A). Myeloid cells from mock-infected mice mainly comprised clusters annotated as AM (Chil3, Ear1, Fapb1; C6), Ly6C− Mo (Ace, Nr4a1, Fcgr4; C3), Ly6C+ Mo (Ccr2, Ly6c2; C4) and dendritic cells (DC; H2-Ab1, Cd209a, Flt3; C9) (Fig. 2A-C and fig. S2A). Few cycling Macs (Birc5, Top2a, Mki67; C11) were also detected in mock-infected mice, along with few CD206− IM (C1qa, C1qc, H2-Ab1, Cd74, Tmem119; C1) and CD206+ IM (C1qa, C1qc, Mrc1, Maf; C8) (29) (Fig. 2A-C and fig. S2A). Ten days post-IAV, AM (C6) disappeared, a small cluster of IAV-associated AM was present instead (Chil3, Ear1, Ear 2; C12), Neu were recruited (S100a8, S100a9, Mmp9; C7), and IM (C1, C8) and iMo (Ccr2, Ly6c2, Irf7; C5) expanded (Fig. 2A-C and fig. S2A). Of note, clusters C2 and C10 were specifically triggered by IAV. C10 had an elevated content in mitochondrial genes and low numbers of detected genes (fig. S2B) and was therefore annotated as dying Macs. C2 expressed significantly higher levels of cathepsins (Ctsb, Ctsz), galectins (Lgals1, Lgals3), Arginase-1 (Arg1) and osteopontin (Spp1) as compared to the other clusters (Fig. 2D). Intracellular flow cytometry staining for Arginase-1 and osteopontin showed that the combined expression of these two proteins was restricted to Ly6G+ Macs at day 10 post-IAV (Fig. 2, E and F), supporting that C2 corresponded to Ly6G+ Macs identified by flow cytometry. Of note, Ly6G+ Macs expressed high levels of Csf1r and Slamf7 (coding for CD319) but did not express any of the neutrophil-related transcripts Csf3r, S100a8, S100a9, Mmp8, Mmp9, Mpo, Slpi or Cd177 (fig. S2, A and C). Ly6g transcripts were not detectable in Neu nor Ly6G+ Macs (fig. S2, A and C). Ly6G+ Macs (C2) displayed both “M1-like” or “M2-like” signature scores and genes and could not be categorized as such on the basis of their expression profile (fig. S2, D and E). Together, our data show that a phenotypically and transcriptionally distinct subset of Ly6G+ Macs emerges during a specific time window corresponding to early weight recovery post-IAV.

Ly6G+ Macs arise from recruited monocytes and are partially dependent on GM-CSF receptor signaling

We next investigated the origin of Ly6G+ Macs and asked whether they could expand via local proliferation. We observed that the percentage of cells positive for the proliferation Ki67 was very low in Ly6G+ Macs (fig. S3A). Next, we treated mice at day 10 post-IAV intraperitoneally (i.p.) with 5-ethynyl-2′-deoxyuridine (EdU) 4 hours before analysis. While the percentage of EdU+ cells was higher in AM compared to all other lung myeloid cells, indicative of active proliferation, virtually no EdU+ cells were detected in Ly6G+ Macs, ruling out their active proliferation (fig. S3, B and C). Third, we assessed whether Ly6G+ Macs arose from the BM or from local progenitor monocytes (30). We generated chimeric mice in which lethally irradiated, thorax-protected CD45.2 WT mice were reconstituted with BM cells from Ms4a3CreR26LSLtdTomato mice (subsequently referred to as Ms4a3tdTom), in which the progeny of granulocyte monocyte progenitors (GMPs) is constitutively labelled (31). At week 4 after transfer, the percentages of tdTomato+ blood Ly6C+ Mo was around 50%, while the percentages of tdTomato+ lung AM and IM were very low, confirming efficient reconstitution and thorax protection (fig. S3, D-G). At day 10 post-IAV, we found that the percentage of tdTomato+ Ly6G+ Macs was similar to that of Ly6C+ Mo (fig. S3, H and I), consistent with a major contribution of BM-derived GMPs, but not local progenitors, to Ly6G+ Macs.

The kinetics of Ly6G+ Mac emergence post-IAV was comparable but delayed compared to that of Ly6C+ Mo and iMo, consistent with the idea that recruited Ly6C+ Mo could give rise to Ly6G+ Macs. Supporting this, Slingshot trajectory analyses of the scRNA-seq data identified two main trajectories starting from Ly6C+ Mo, transiting through iMo to give rise to either IM-like cells or Ly6G+ Macs (32) (Fig. 3A). Genes that exhibited the same pattern of downregulation along pseudotime in each trajectory encompassed the classical monocyte-associated genes Ccr2 and Ly6c2 (Fig. 3B). We also found, in both trajectories, a time-restricted upregulation of interferon-stimulated genes (Ifi209, Ifitm3, Ifi47, Isg15) that likely corresponded to transitioning iMo (Fig. 3B). Finally, trajectory-specific genes were gradually and specifically upregulated along pseudotime in the IM-like (e.g., C1qa, C1qc, C1qb, Mrc1, Cd74, H2-Ab1, H2-Eb1) or the Ly6G+ Mac (e.g., Arg1, Spp1, Ccl2, Ccl7, Ctsb) trajectories (Fig. 3B).

Next, we infected the monocyte fate-mapper mice Ms4a3tdTom (31) and Cx3cr1GFP (33) with IAV and found that, at day 10 post-IAV, Ly6G+ Macs were Tomato+ and GFPhi, respectively, consistent with their GMP-derived monocytic origin (Fig. 3, C and D). To address the dependency of Ly6G+ Mac on Ccr2-dependent BM-derived Ly6C+ Mo (34), we generated BM competitive chimeras in which lethally irradiated CD45.1/CD45.2 WT mice were engrafted with a 1:1 mix of CD45.2 Ccr2−/− and CD45.2 Ms4a3tdTom BM cells. At week 4 after reconstitution, most blood Ly6C+ Mo were of donor Tomato+ origin, as expected (fig. S3, J and K) (34). At day 10 post-IAV, the majority of Ly6G+ Macs was also of donor Ms4a3tdTom origin, indicating their dependency on Ccr2 (Fig. 3, E and F).

We next sought to assess the fate and lifespan of Ly6G+ Macs. The abundance of Ly6G+ Macs during a limited time window post-IAV suggested that they might be short-lived. In line with this, trajectory analyses in Fig. 3A suggested that Ly6G+ Macs (C2) gave rise to Macs with low RNA content (C10) (fig. S2B). We performed EdU pulse experiments at day 7 post-IAV and found that Ly6G+ Macs staining positive for EdU at day 10 post-IAV were completely cleared from the lung at day 17 (Fig. 3G). Annexin V/Propidium iodide (PI) staining at day 10 post-IAV further supported that Ly6C+ Mo transitioning to iMo and differentiating into Ly6G+ Macs became progressively more sensitive to death, with a substantial portion of Ly6G+ Macs being either early or late apoptotic (Annexin V+/PI+/−), or necrotic (Annexin V−/PI+) (Fig. 3H and fig. S3L). Hence, our data show that Ly6G+ Macs represent a short-lived Mac subset arising from Ccr2-dependent Ly6C+ Mo.

The Ly6G signal on Ly6G+ Macs, considered to be neutrophil-specific, was unexpected and required careful validation. First, we verified that the Ly6G fluorescence intensity was virtually absent in unstained or isotype antibody (Ab)-stained CD11b+ cells from IAV-infected WT mice (fig. S4A). Second, we evaluated whether Ly6G+ Mac precursors, namely iMo, could intrinsically upregulate Ly6G on their surface when isolated from the lungs of IAV-infected mice at day 10 post-IAV. We found that iMo from IAV-infected mice and stimulated ex vivo with GM-CSF, and to a less extent with M-CSF, upregulated Ly6G protein on their surface (Fig. 3, I and J). Importantly, we also found that lung iMo isolated from IAV-infected Ly6gCreERT2R26LSLtdTomato mice (Ly6gtdTom) treated with tamoxifen and stimulated with GM-CSF ex vivo became tdTomato+, indicative of active Ly6g gene transcription in monocytic cells (Fig. 3, K and L). These data demonstrate that Ly6G can be actively expressed by Mo-Macs.

Given the ability of GM-CSF to trigger Ly6G expression on lung iMo isolated from IAV-infected mice, we assessed the dependency of Ly6G+ Macs on GM-CSF receptor signaling in vivo. We generated BM competitive chimeras in which thorax-protected, lethally irradiated CD45.1/CD45.2 WT mice were engrafted with a 1:1 mix of CD45.1 Csf2ra-/- and CD45.2 Csf2ra+/+ BM donor cells. At week 4 after reconstitution, blood Neu and Mo of donor origin arose equally from CD45.1 Csf2ra−/− and CD45.2 Csf2ra+/+ BM cells (fig. S4, B and C). At day 10 post-IAV, we found that CD45.2 Csf2ra+/+ BM cells had a competitive advantage over CD45.1 Csf2ra−/− BM cells to become Ly6G+ Macs, which was not observed among most other lung myeloid cells (Fig. 3M and fig. S4C-E), and the percentage of Arg-1+ cells was lower in Csf2ra−/− Ly6G+ Macs as compared to Csf2ra−/− Ly6G+ Macs (Fig. 3, N and O), indicating their partial dependency on GM-CSF receptor signaling.

Ly6G+ Macs exhibit distinct ultrastructural features associated with elevated metabolic and phagocytic abilities

Gene Set Enrichment Analyses (GSEA) of the transcriptomic profile of Ly6G+ Macs (C2) compared to all other clusters identified a response to interferon-γ and cytokines, chemotactic and viral processes, an active metabolic state, a highly developed endomembrane system and elevated phagocytic abilities in Ly6G+ Macs (Fig. 4A). We analyzed FACS-sorted Ly6G+ Macs by transmission electron microscopy at day 10 post-IAV and found that they exhibited a kidney-shaped nucleus, microvilli-rich membrane, secretory granules and a cytoplasm rich in rugous endoplasmic reticulum (RER), Golgi apparatus, lysosomes and autophagy vacuoles, distinguishing them from Neu, Ly6C+ Mo and IM-like cells (Fig. 4B). The morphology of Ly6G+ Macs was reminiscent to that of an atypical population of monocytes, called SatM monocytes, which arise from particular GMPs during the fibrotic phase in lungs post-bleomycin, contribute to fibrosis and were regulated by CCAAT/enhancer binding protein β (C/EBPβ) (35). We conducted a single-cell regulatory network inference and clustering (SCENIC) analysis (36) and found that C/EBPβ activity was lower in Ly6G+ Macs compared to Ly6C− Mo, Ly6C+ Mo, iMo and Neu (fig. S5A). Moreover, we generated a SatM signature score based on the genes unregulated in SatM monocytes (35), mapped such score to our scRNA-seq data and found that Ly6G+ Macs displayed a lower SatM score compared to IM, Ly6C+ Mo, Ly6C- Mo or iMo (fig. S5, B and C), supporting that Ly6G+ Macs are not dependent on C/EBPβ and are transcriptionally distinct from SatM monocytes (35).

Next, we characterized the metabolic profile of Ly6G+ Macs using a metabolic flux assay and found that the extracellular acidification rate (ECAR) was higher in Ly6G+ Mac compared to IM-like cells, both at baseline and under stress, supporting that the glycolytic pathway was highly active in Ly6G+ Macs (Fig. 4C). Moreover, while the basal mitochondrial oxygen consumption rate (OCR) was similar between Ly6G+ Mac and IM-like cells, the OCR under stress conditions was higher in Ly6G+ Macs compared to IM-like cells (Fig. 4D), supporting that they possessed a high metabolic potential (Fig. 4E).

We administered infected mice with fluorescently-labeled E. coli particles intratracheally (i.t.) at day 10 post-IAV and confirmed that Ly6G+ Macs were highly phagocytic compared to Neu, Ly6C+ Mo, IM-like cells and AM (Fig. 4, F and G). Hence, we asked whether Ly6G+ Macs could perform efferocytosis in vivo. To this end, we infected BM competitive chimeras in which lethally irradiated CD45.2 WT mice were engrafted with a 1:1 mix of Cx3cr1GFP/+ and Ms4a3tdTom BM cells. Of note, at day 10 post-IAV, 60% of Cx3cr1GFP/+ Ly6G+ Macs were also tdTomato+ (Fig. 4, H and I and fig. S6, A-C), demonstrating that they were highly potent in engulfing myeloid cells in vivo. To assess whether Ly6G molecules could be transferred from Neu to Ly6G+ Macs during efferocytosis, we infected BM competitive chimeras in which lethally irradiated CD45.1/CD45.2 WT mice were engrafted with a 1:1 mix of Cx3cr1GFP/+ cells and Ly6g−/− (i.e., homozygous Ly6gCreERT2 mice) (37) or Ly6g+/+ BM cells. At day 10 post-IAV, we found that the levels of Ly6G on Cx3cr1GFP+ Ly6G+ Macs from Cx3cr1GFP+: Ly6g−/− BM chimeric mice, in which half of the Neu were Ly6g−/−, were similar to those from Cx3cr1GFP+: Ly6g+/+ BM chimeric mice, supporting no evidence for a Ly6G transfer from Neu to Ly6G+ Macs (fig. S6, D and E). Together, these data show that IAV-triggered Ly6G+ Macs are characterized by metabolic, morphological and efferocytic properties distinct from other lung myeloid cells.

Ly6G+ Macs populate the alveolar lumen of perilesional regenerating areas

Next, we investigated the localization and the spatial organization of Ly6G+ Macs. First, we performed confocal microscopy staining of lung sections of infected Cx3cr1GFP mice at day 10 post-IAV. By defining Ly6G+ Macs as cells double positive for Ly6G and GFP, we found that they were located in the alveolar lumen (Fig. 5A), which was also confirmed by in situ electron microscopy (Fig. 5B).

To further investigate the spatial distribution of Ly6G+ Macs and the molecular signatures of Ly6G+ Macs-rich areas, we performed spatial transcriptomic analyses using GeoMx Digital Spatial Profiler (DSP), which allows whole-genome transcript analyses within regions of interest (ROIs). Lung sections from 2 mock- and 4 IAV-infected mice were collected at day 10 post-IAV, stained with anti-CD68 and anti-Ly6G antibodies, and ROIs were selected in control lungs (4 ROIs), extralesional zones (4 ROIs), intralesional zones (5 ROIs) and zones rich in Ly6G+CD68+ cells that were mostly located in the periphery of consolidated areas (perilesional, 11 ROIs) (Fig. 5C and fig. S7A). Unsupervised principal component (PC) analysis showed that perilesional ROIs were separated from the other regions (Fig. 5D). Volcano plot representation of the differentially expressed (DE) genes between conditions and the heatmap of the 522 significantly upregulated genes in perilesional areas compared to the other areas supported that perilesional areas were transcriptionally very active (fig. S7, B and C). General cellular deconvolution indicated that perilesional zones were also enriched in tissue Macs as compared to the other regions (fig. S7D). We next mapped cell signature scores of lung myeloid cell populations analyzed by scRNA-seq to the ROIs and confirmed that perilesional zones were enriched in Ly6G+ Mac compared to the other zones (Fig. 5E and fig. S7E). GSEA analyses indicated that perilesional areas were enriched in biological responses related to cytoskeleton activity, epithelial cell migration and elevated metabolic activity compared to intralesional areas, consistent with intense remodeling activities (Fig. 5F).

Next, we took advantage of a publicly available scRNA-seq dataset of alveolar epithelial cell states present during alveolar regeneration after bleomycin-induced lung injury (38) and containing type 1 and type 2 alveolar epithelial cells (AT1 and AT2, respectively), as well as transitional states appearing during AT2 to AT1 differentiation, called primed AT2 and damage-associated transient progenitors (DATPs) (38). By mapping the signature scores of such transitional epithelial cell states to the ROIs, we found that Ly6G+ Mac-rich perilesional zones were enriched in primed AT2 and DATPs compared to control and intralesional zones (fig. S7, F and G). Accordingly, the Ly6G+ Mac score correlated positively with those of primed AT2 and DATPs (Fig. 5G). We also confirmed by confocal microscopy that Ly6G+ Macs were particularly abundant in the periphery of consolidated areas and clustered with AT2 cells, while intralesional consolidated areas, which exhibited low levels of staining for AT1 and AT2, contained few Ly6G+ Macs (Fig. 5H and fig. S8). Altogether, these data are consistent with perilesional areas serving as the site of active epithelial regeneration post IAV, and that Ly6G+ Macs, which cluster in such areas, contribute to this process.

Ly6G+ Macs promote alveolar epithelial regeneration through IL-4R signaling

To formally assess the function of Ly6G+ Macs in vivo, we aimed to generate a transgenic mouse strain in which Ly6G+ Mac differentiation was impaired. Thus, we applied the SCENIC algorithm to our scRNA-seq data to map gene regulatory networks and predict transcription factor activities in Ly6G+ Macs (36). Of note, c-Maf and MafB exhibited a high regulon activity in Ly6G+ Macs and IM-like cells, as described (39), but not in other lung myeloid cells (Fig. 6A and fig. S9). Elevated c-Maf and MafB protein levels were also detected in Ly6G+ Macs at day 10 post-IAV by flow cytometry (Fig. 6, B and C), and Maf and Mafb transcript levels were elevated in lung Ly6G+ Mac-rich perilesional areas of IAV-infected mice (Fig. 6D). We generated mice with myeloid-restricted Maf and Mafb deficiency by crossing Maf and Mafb floxed mice (Maf/Mafbfl/fl) with mice constitutively expressing Cre recombinase under the control of the lysozyme M promoter (Lyz2Cre), called Maf/MafbMyeloKO mice hereafter. At day 10 post-IAV, Maf/MafbMyeloKO mice showed a virtual absence of Ly6G+ Macs, while numbers of Neu, AM, IM-like cells, Ly6C+ Mo and iMo were similar and numbers of Ly6C− Mo were higher compared to control mice (Fig. 6E). Hence, we employed this model to address the consequences of Ly6G+ Mac deficiency on viral control, morbidity and alveolar epithelial repair following IAV infection.

We assessed the levels of lung mRNA coding for the non-structural influenza protein NS1 post-IAV and found that they were not significantly different between Maf/MafbMyeloKO and controls and returned to baseline at day 10 post-IAV (Fig. 6F), supporting that Ly6G+ Macs did not substantially influence host viral control. However, Maf/MafbMyeloKO mice lost more weight post-IAV compared to controls (Fig. 6G), suggestive of a more severe IAV-induced pathology.

Histopathological analyses of lung sections at day 20 post-IAV indicated broader lesional areas in Maf/MafbMyeloKO mice compared to controls, as well as more pronounced dysplastic repair and bronchiolization of the alveoli, based on quantification of mucus area in lung lesional areas (Fig. 6H-J). These results suggested that, in the absence of Ly6G+ Macs, the classical pathway of alveolar epithelial regeneration involving progenitor AT2 expansion and differentiation towards AT1 (40, 41) was defective and compensated by dysplastic repair. Next, we evaluated the numbers of AT1, AT2 and regenerating AT2 (regAT2) at day 20 post-IAV in Maf/MafbMyeloKO and control mice by flow cytometry and observed a significant decrease in the numbers of AT2 and regAT2 in Maf/MafbMyeloKO mice compared to controls (Fig. 6, K and L), confirming that AT2 were less able to expand and differentiate into AT1 in the absence of Ly6G+ Macs. Of note, i.t. transfer of Ly6G+ Macs isolated from lungs of WT mice at day 10 post-IAV into IAV-infected Maf/MafbMyeloKO mice partially improved weight recovery and restored numbers of AT2 to the levels observed in IAV-infected control mice (Fig. 6, M and N). These results suggest that Ly6G+ Macs are key players of euplastic epithelial regeneration after IAV-induced lung injury.

To determine whether Ly6G+ Mac-rich perilesional areas were imprinted by a type 2 reparative environment (19, 42, 43), we mapped a type 2 signature score based on genes involved in IL-4 receptor downstream signaling pathways to the DSP spatial transcriptomic data and found that perilesional areas exhibited the highest type 2 score as compared to the other ROIs (Fig. 7A). Hence, we asked whether IL-4 receptor signaling, whose activation is known to induce a repair phenotype in Macs (19, 42, 43), was involved in Ly6G+ Mac identity and function. First, we found that Ly6G+ Macs expressed high levels of the IL-4 receptor α chain (IL-4Rα) (Fig. 7, B and C). Next, we generated BM competitive chimeras in which lethally irradiated CD45.1/CD45.2 WT mice were engrafted with a 1:1 mix of CD45.2 Il4ra−/− and CD45.2 Ms4a3tdTom BM cells. At week 4 after reconstitution, efficient BM reconstitution was confirmed in the blood (fig. S10, A and B). At day 10 post-IAV, we found that Ly6G+ Macs of donor Ms4a3tdTom origin exhibited a competitive advantage over those of donor Il4ra−/− origin, which was not observed among other lung myeloid cells (Fig. 7D and fig. S10, C and D), and the remaining Il4ra−/− Ly6G+ Macs were impaired in their ability to express Arg-1 (Fig. 7, E and F). Finally, WT chimeric mice fully reconstituted with Il4ra−/− or WT BM cells were generated (Il4ra−/− BM -> WT or WT BM -> WT, respectively) and infected with IAV. We found, like in Maf/MafbMyeloKO mice, that Il4ra−/− BM -> WT mice had an impaired recovery post-IAV compared to WT BM -> WT mice (Fig. 7G). These data suggest that Ly6G+ Macs exert their function via IL-4R-dependent pathways, at least in part.

Finally, we asked whether Ly6G+ Macs could directly influence AT2 fate and whether cell-cell contacts were needed. To this end, we performed a scratch assay in vitro using the MLE-12 mouse AT2 cell line and evaluated the confluence of AT2 cells 12 hours post-scratch in the presence or absence of Ly6G+ Macs isolated from infected lung at day 10 post-IAV. Co-culture with Ly6G+ Macs, but not with Neu, IM-like cells or iMo was associated with an increase in cell confluence (Fig. 7H), indicating that Ly6G+ Macs can directly and specifically promote wound healing in vitro. A similar scratch assay was also performed using conditioned medium (CM) from FACS-sorted Ly6G+ Macs that were cultured overnight with or without the type 2 cytokines IL-4 and IL-13. In this setting, CM from IL-4/13-pulsed Ly6G+ Macs could promote wound healing compared to control medium (containing only IL-4 and IL-13) or CM from unpulsed Ly6G+ Macs (Fig. 7, I and J). We performed proteome profiling on such CM and found that Ly6G+ Macs were highly potent in secreting soluble factors, among which were chemokines (CCL5, CXCL16, CCL12, CXCL10), cytokines (TNF-α, IL-10, IL-1α) and osteopontin, some of which were increased upon IL-4R activation (Fig. 7K). Some of the molecules detected in the CM of Ly6G+ Macs had their transcript levels significantly upregulated in Ly6G+ Macs (C2) as compared to other clusters (fig. S11). Altogether, our data demonstrate that Ly6G+ Macs can release soluble factors upon IL-4 receptor triggering that act directly on AT2 to promote epithelial regeneration (fig. S12).

Ly6G+ Macs belong to a conserved host response to injury across organs, triggers and species

We evaluated whether Ly6G+ Macs were specifically recruited in the IAV model or were also triggered in other models of injury. First, we used a model of non-infectious lung injury based on bleomycin (bleo) instillation and performed time-course flow cytometry analyses. Ly6G+ Macs expressing high levels of Arg-1 and CXCR4 were mostly present between day 7 and day 14 post-bleo, which correlated with signs of epithelial damage, as reflected by the decrease in numbers of AT1 and AT2 (Fig. 8A-E). We also found similar Ly6G+ Macs peaking at days 1 and 2 post-treatment in an acute model of acetaminophen-induced liver injury, which correlated with the release of alanine aminotransaminase (ALT) in plasma (Fig. 8F-I). Our data thus suggested that Ly6G+ Macs are a component of a conserved response to tissue damage, regardless of the organ or trigger.

Finally, we asked whether Macs sharing a similar transcriptomic signature were also present in the broncho-alveolar lavage fluid (BALF) of diseased humans. We performed scRNA-seq analyses of BALF cells from 7 patients with a suspicion of pneumonia and manually annotated the cell clusters based on the most upregulated genes (Fig. 8, J and K and fig. S13). Next, we mapped a Ly6G+ Mac score based on orthologous genes in humans to the BALF cells and found that cells exhibiting the highest Ly6G+ Mac score belonged to the same cluster C9 identified as Mo-Macs based on their high expression of monocyte genes and their low expression of AM-associated genes (Fig. 8, L and M). SCENIC analyses (36) predicted higher MAF and MAFB activities in the Mo-Mac cluster compared to other clusters (Fig. 8N), further supporting that the airspace of human pneumonia lungs contains Mo-Macs that are transcriptionally similar to mouse Ly6G+ Macs.

Discussion

Restoration of gas exchanges after lung injury is critical for life and relies on appropriate regulation of inflammation and regeneration of the damaged alveoli. While recent progress has been made in understanding the epithelial-intrinsic mechanisms underlying alveolar regeneration post-injury (3, 38, 41, 44–47), an important gap resides in our understanding of the innate immune-epithelial crosstalk taking place to promote epithelial repair and host recovery. While recruited Mo-Macs are often seen as culprits and drivers of disease progression in different contexts, such as in Covid-19, interstitial fibrosis or lung cancer (13, 17, 18, 48), advances in single cell and spatial technologies have enabled new opportunities to investigate the spatiotemporal regulation of Mo-Mac responses in-depth. Here, using such approaches combined with lineage tracing, BM chimeras, gene targeting, multi-parameter flow cytometry and imaging, our work identifies a previously undescribed atypical population of short-lived recruited Ly6G+ Macs that critically contributes to alveolar epithelial regeneration post-injury in mice.

To our knowledge, there is no report of such Ly6G+ Macs in the literature. Of note, Ly6G is largely considered as a neutrophil-specific marker (49), and, in many studies employing cytometry, anti-Ly6G antibodies are included to gate out “neutrophils” before gating on Macs. Hence, Ly6G+ Macs might have been previously overlooked and considered as part of the neutrophil compartment. Here, we provide evidence that the Ly6G signal is specific and that lung iMo, the precursors of Ly6G+ Macs, can actively express Ly6G gene and protein upon GM-CSF stimulation ex vivo. Of note, GM-CSF is mainly produced by AT2 cells (41, 50), which are located in the vicinity of Ly6G+ Macs, and Ly6G+ Macs were dependent on GM-CSF receptor signaling for their generation and Arg-1 expression post-IAV infection in vivo. Further supporting an intrinsic upregulation of Ly6G, we found no evidence of Ly6G protein transfer from neutrophils to Ly6G+ Macs. In addition to Ly6G expression, Ly6G+ Macs were phenotypically and transcriptionally distinct from neutrophils and exhibited key developmental, phenotypic and transcriptomic macrophage features including their dependency on Ccr2, elevated Cx3cr1 expression, and their dependency on the transcription factors c-Maf and MafB. However, Ly6g does not appear as a target gene of c-Maf and MafB in the ChiP atlas, suggesting alternative gene regulatory mechanisms.

We found that Ly6G+ Macs were transiently recruited to the alveolar spaces of particular lung areas and remained phenotypically and transcriptionally distinct from tissue-resident AM. Our data thus suggest that the local microenvironment of Ly6G+ Macs, which can shape Mac identity (13, 51–53), is dynamically regulated and distinct from that of AM. First, we found that Ly6G+ Macs originate from BM-derived monocytes recruited to the lung in a Ccr2-dependent manner. Second, we showed that such inflammatory monocytes differentiating into Ly6G+ Macs could phagocytose GMP-derived myeloid cells, such as neutrophils, in vivo. Interestingly, such a process has been shown to trigger a metabolic rewiring that is associated with Arg-1 activity and aerobic respiration and is important for the resolution of inflammation and tissue repair (54, 55). Third, we provided evidence that GM-CSF and type 2 cytokine signaling through the IL-4 receptor are involved in generation of Ly6G+ Macs and their function. While the cellular source of type 2 cytokines remains unknown, a peak of T helper type 2 (Th2) and type 2 innate lymphoid cells (ILC2s) has been reported in lungs of IAV-infected mice around day 10 post-infection (56), a time point that coincides with the peak of Ly6G+ Macs. Moreover, Ly6G+ Macs can release the Th2-attracting and ILC2-activating signals CCL22 and IL-33, respectively, consistent with the idea that Ly6G+ Macs can contribute to the type 2 milieu that promotes their repair phenotype (57–59). Fourth, we found that Ly6G+ Macs were spatially-restricted to perilesional areas, zones that were enriched in transitional epithelial cell states involved in AT2-mediated alveolar regeneration (38). Perilesional areas were also sites of intense cytoskeleton activity, aerobic respiration, extracellular matrix deposition and cell migration, all of which are involved in active alveolar epithelial regeneration (3, 9, 41). Our findings are in line with a previous report that identified damaged zones of IAV-infected lungs that were in the periphery of consolidated areas and were sites of active tissue regeneration and AT2 cell proliferation and differentiation (47).

By disrupting myeloid-specific cMaf and MafB-dependent pathways, we obtained Maf/MafbMyleoKO mice in which Ly6G+ Macs were no longer able to differentiate post-IAV, thus representing a valuable tool to address their functions. We identified Ly6G+ Macs as essential actors to license optimal alveolar epithelial regeneration requiring the differentiation of progenitor AT2 cells towards AT1 cells (41). In the absence of Ly6G+ Macs, AT2-to-AT1 transitioning cells were virtually absent and a more pronounced dysplastic alveolar repair associated with a bronchiolization of the alveoli was observed. Such epithelial phenotype was associated with exacerbated morbidity and is reminiscent of what is observed in severe forms of respiratory viral infections (3, 41). Additionally, local adoptive transfer of Ly6G+ Macs in IAV-infected Maf/MafbMyeloKO mice improved weight recovery and restored AT2 numbers to levels seen in IAV-infected control mice, suggesting that Ly6G+ Macs could support AT2 expansion post-injury. We further dissected the underlying mechanisms ex vivo and found that Ly6G+ Macs could directly promote wound healing of murine AT2 cells through IL-4 receptor-mediated release of soluble factors. Among these factors, osteopontin is expressed by Ly6G+ Macs and is a ligand for the receptor CD44 (60). Of note, CD44hi AT2 cells represent a subset of AT2 cells with stem cell properties (61), consistent with the hypothesis that the osteopontin-CD44 axis might trigger alveolar regeneration, even though osteopontin release was not potentiated by type 2 cytokine stimulation of Ly6G+ Macs under the experimental conditions tested. The elevated Arg-1 expression by Ly6G+ Macs could also influence AT2 cells, either via the local deprivation of L-Arginine or the generation of ornithine and polyamines (62, 63). Finally, Ly6G+ Macs can also secrete pro-inflammatory cytokines, such as IL-1 and TNF-α, which have been shown to support alveolar regeneration (38, 64).

In the last part of the work, we provided evidence that Macs similar to Ly6G+ Macs are part of a host response to injury that is independent of the organ or the initial trigger and that is conserved across species. Indeed, we found that Ly6G+ Macs were also recruited in a bleomycin-induced model of non-infectious lung injury and in a model of acute acetaminophen-induced liver injury in mice. Even though Ly6G+ Mac numbers were lower in these models compared to the IAV model, it is noteworthy that the peak of Ly6G+ Macs correlated with the presence of damage, as attested by the drop in AT1 and AT2 in the lungs or the release of the hepatic enzyme ALT in the liver. These data are consistent with the idea that they may contribute to tissue repair in these models as well, even though it remains to be addressed experimentally. Interestingly, we also reported the presence of transcriptionally similar Macs in the BALF of pneumonia patients by performing scRNA-seq analyses of BALF cells. Of note, cells exhibiting an elevated Ly6G+ Mac score belonged to a cluster identified as Mo-Macs, and SCENIC analyses suggested that such Mo-Macs displayed high MAF and MAFB activities, reminiscent of Ly6G+ Macs in mice. Speculating that such Mo-Macs also depend on GM-CSF and exhibit similar reparative functions in human lungs, our findings provide a rationale to investigate the benefits of inhaled GM-CSF, beyond the restoration of the AM niche (65), to improve epithelial regeneration after severe viral-induced disorders. Indeed, one could speculate that the fate or functions of Ly6G+ Macs are modified and become dysregulated in uncontrolled forms of respiratory viral infections or in chronic fibrotic diseases. Of note, SPP1, which encodes osteopontin, has been linked to fibrosis and is often used as a proxy for “pro-fibrotic”, pathogenic Macs (48, 66–68). Our data support that Spp1+ Mo-Macs, like Ly6G+ Macs, can also exert beneficial roles, while other Spp1+ Mo-Macs can become dysregulated, persistent and pathogenic, like in chronic Covid-19 or idiopathic pulmonary fibrosis. Understanding what drives beneficial or pathological responses of SPP1+ Macs represents an interesting avenue for future research.

No therapeutic options exist to promote lung regeneration so far. By characterizing in-depth a short-lived atypical Mac population that licenses alveolar regeneration post-injury, our findings could serve as a basis to devise novel myeloid-centered regenerative strategies for medically relevant conditions such as severe or chronic respiratory viral infections or acute respiratory distress syndrome.

Materials and Methods

Study design

In this study, we investigated the spatio-temporal distribution, transcriptional regulation, fate, identity and function of Ly6G+ Macs in an infectious model of lung injury. To this end, we employed flow cytometry, microscopy, single-cell and spatial transcriptomic approaches, bone marrow (BM) chimeras, monocyte fate-mapping and gene targeting. In most of the mouse experiments, 4 to 10 mice per group per time point were used to identify differences between groups with at least 80% power and 5% significance level. In some experiments, no statistical methods were used to pre-determine sample sizes but our sample sizes are similar to those reported in previous publications (39, 69, 74, 75). No statistical methods to pre-determine sample size were used for the analyses of human BALF cells from pneumonia patients. Data from independent experiments were pooled for analysis in each data panel, unless otherwise indicated. No data were excluded from the analyses and all attempts at replication were successful and gave similar results. Histopathological examination of lung sections was blinded. Allocation of animals into experimental groups was done randomly at the start of the experiments. The specific numbers of mice, the number of experimental replicates and the statistical tests performed are all included in the respective figure legends.

Mice

All experiments, unless otherwise specified, were performed on age-matched 8–12-wk-old male and female mice on the C57BL/6 background. Details about the transgenic strains can be found in the Supplementary Materials. Mice were housed under specific pathogen-free conditions and maintained in a 12-h light–dark cycle with food and water ad libitum. All animal experiments described in this study were carried out in an animal biosafety level 3 containment unit. Experiments were reviewed and approved by the Institutional Animal Care and Use Committee of the University of Liège (ethical approval #2276). The ‘Guide for the Care and Use of Laboratory Animals,’ prepared by the Institute of Laboratory Animal Resources, National Research Council, and published by the National Academy Press, as well as European and local legislations, was followed carefully. Accordingly, the temperature and relative humidity were 21°C and 45-60%, respectively

In vivo models of injury

The mouse-adapted influenza strain A/Puerto Rico/8/34 (H1N1; PR8) was kindly provided by F. Trottein (Institut Pasteur, France). The viral stock suspension (108 Plaque Forming Units [PFU] ml–1) was diluted and 5 PFU were administered intranasally (i.n.) to isoflurane-anesthetized mice in 50 μl of PBS (Thermo Fisher). Control groups received an equal volume of PBS i.n. for mock infection.

For bleomycin-induced lung injury, isoflurane-anesthetized mice were treated intratracheally (i.t.) with a single instillation of 0.06 IU of bleomycin (Bio-Connect) in a volume of 50μl PBS. Control animals received 50μl PBS alone.

For acetaminophen-induced liver injury, mice were fasted during 15 hours with free access to water and were injected intraperitoneally (i.p.) with 300mg kg-1 of acetaminophen (Sigma) in saline solution (NaCl 0.9%). Free access to food was allowed after treatment.

Reagents and antibodies

A complete list of the reagents and antibodies used in this manuscript can be found in Tables S2 and S3, respectively.

Flow cytometry

Staining reactions were performed in the dark at 4°C for 30 minutes with 2% v/v of Fc block (BD Biosciences) to avoid nonspecific binding. For intracellular staining, extracellular-stained cells were fixed and permeabilized with the Foxp3/Transcription factor Staining Buffer Set (Thermo Fisher). For EdU staining, extracellular-stained cells were permeabilized and stained using Click-iT EdU Alexa Fluor 488 Flow Cytometry Assay Kit (Thermo Fisher), according to the manufacturer’s instructions.

Cell viability was assessed using 7-AAD (BD Bioscience) or Fixable Viability Dye eFluor™ 780 (Thermo Fisher). Cell suspensions was analysed with a FACSCANTO II or a LSRFortessa (BD Biosciences). Results were analyzed using FlowJo software (Tree Star). For scRNA-seq, transmission electron microscopy, cytological examination and ex vivo experiments, lung myeloid cells were sorted using a FACSAria III (BD Biosciences) or a Sony MA900.

In vivo treatments

For EdU incorporation experiments shown in fig. S3, B and C, mice were injected i.p. at day 10 post-IAV with 1 mg EdU (Santa Cruz Biotechnology) in 200 μl PBS 4 hours before sacrifice. For experiments addressing the lifespan of Ly6G+ Macs (Fig. 3G), 1mg EdU in 200 μl PBS was injected i.p. twice 5 hours apart at day 7 post-IAV, and EdU incorporation was evaluated in blood leucocytes at day 8 post-IAV. The incorporation of EdU in lung myeloid cells was evaluated at days 10, 14 and 17 post-IAV. Assessment of phagocytic activity was performed as previously described (69). Briefly, isoflurane-anesthetized mice were instilled i.t. with 2.108 pHrodo™ Green E. coli BioParticules (Thermo Fisher) in 100 μl PBS. Lungs were harvested 3 hours later for flow cytometry analyses.

Generation of BM (competitive) chimeras

CD45.2, CD45.1 or CD45.1/CD45.2 WT mice were anesthetized by i.p. injection of 200 μl PBS containing ketamine (Nimatek, Dechra, 75 mg kg-1) and xylazine (Rompun, Bayer, 10 mg kg-1). When mentioned, the thoracic cavity was protected with a 0.6-cm-thick lead cover. Mice were irradiated with two consecutive doses of 6 Gy 15 minutes apart. Once recovered from the anaesthesia, mice were reconstituted by intravenous (i.v.) administration of 2.106 BM cells from Ms4a3tdtom or Il4ra−/− mice, for full chimeras. For mixed BM chimeras, mice were reconstituted i.v. with 2.106 BM cells consisting of a 1:1 mix of BM cells obtained from the following mice: CD45.1 WT, Ms4a3tdtom, Ccr2−/−, CD45.1 Csf2ra−/−, CD45.2 Csf2ra+/+, Cx3cr1GFP+, Il4ra−/−, or homozygous Ly6gCreERT2 mice (also called Ly6g−/− mice). From the day of irradiation, mice were treated for 4 weeks with 0.05 mg ml−1 of enrofloxacin (Baytril, Bayer) in drinking water. Chimerism was assessed by flow cytometry in the blood 4 weeks after irradiation.

scRNA-sequencing and analyses

For mouse experiments, lung myeloid cells were FACS-sorted as living singlet CD45+, F4/80+ and/or CD11b+ cells from lung single-cell suspensions pooled from 5 mock-infected and IAV-infected C57BL/6 male WT mice at day 10 post-IAV. For each sample, an aliquot of Trypan blue-treated cells was examined under the microscope for counting, viability and aggregate assessment following FACS sorting. Viability was above 90% for all samples and no aggregates were observed. Cell preparations were centrifuged and pellets were resuspended in calcium- and magnesium-free PBS containing 0.4 mg ml−1 UltraPure BSA (Thermo Fisher Scientific). The 10X Genomics platform (Single Cell 3’ Solution) was used. For library preparation, approximately 2,000 (Mock group) and 6,000 (IAV group) cells were loaded into the Chromium Controller, in which they were partitioned, their polyA RNAs captured and barcoded using Chromium Single Cell 3’ GEM, Library & Gel Bead Kit v3 (10X Genomics). The cDNAs were amplified and libraries compatible with Illumina sequencers were generated using Chromium Single Cell 3’ GEM, Library & Gel Bead Kit v3 (10X Genomics).

For human BALF cell analyses, chromium Fixed RNA Profiling for multiplexed samples (10X Genomics) was used for scRNA-seq analysis of human BALF cells, allowing the storage of fixed cells and enabling analysis of multiple samples in one single GEM reaction. Fresh samples were directly fixed in a 4% formaldehyde solution after collection for storage at -80°C. For GEM creation, the Multiplex-compatible Chromium Next GEM Single Cell Fixed RNA Human Transcriptome Probe Kit including a Probe Barcode that permits sample multiplexing and subsequent demultiplexing was used.

Details about scRNA-seq analyses can be found in the Supplementary Materials.

Transmission electron microscopy

FACS-sorted myeloid cell populations or lung tissues from IAV-infected mice at day 10 post-IAV were fixed in 2.5% glutaraldehyde (diluted in Sorensen’s buffer: 0.1 M Na2HPO4/NaH2PO4 buffer, pH 7.4) for 1h at 4 °C and postfixed for 30 min in 2% OsO4 (diluted in 0.1 M Sorensen’s Buffer). After dehydration in graded ethanol, samples were embedded in Epon resin. Ultrathin sections obtained with a Reichert Ultracut S ultramicrotome (Reichert Technologies) were contrasted with 2% uranyl acetate and 4% lead citrate. For ultrastructural analyses, random fields of cells were examined under a Jeol TEM JEM-1400 Transmission Electron Microscope at 80 kV, and photographed using an 11-megapixel camera system (Quemesa, Olympus).

Extracellular flux analysis

Oxygen consumption rate (OCR) was measured using Seahorse XF Cell Mito Stress Test (Agilent) according to manufacturer’s recommendations and as described previously (70, 71). Briefly, Neu, IM-like cells and Ly6G+ Macs were FACS-sorted at day 10 post-IAV and seeded (10.104, 7.104 and 8.104 cells/well, respectively) in XFp mini-plates (Agilent) pre-coated with CellTak. Cells were kept in unbuffered serum-free DMEM supplemented with pyruvate (1mM), glutamine (2mM), glucose (10mM), at pH 7.4, 37 °C and ambient CO2 for 1h before the assay. Analysis was performed using the XFp analyser (Seahorse Bioscience) as per manufacturer’s instructions. Additional details can be found in the Supplementary Materials.

Spatial transcriptomic analyses using Digital Spatial Profiling (DSP)

Five-μm-thick formalin-fixed, paraffin-embedded (FFPE) sections were prepared using the protocol from NanoString Technologies. Briefly, 2 tissue slides, each containing 1 mock and 2 IAV samples harvested 10 days post-IAV, were analyzed. Slides were first stained with antibodies against CD68, Ly6G (clone 1A8), and DNA was visualized with 500 nM Syto83. Mouse Whole Transcriptome Atlas probes targeting more than 19,000 targets were hybridized, and slides were loaded on the GeoMx DSP. Briefly, entire slides were imaged at x20 magnification, and Regions of Interest (ROIs) were chosen based on serial Hematoxylin & Eosin sections and on morphological markers to select lesional, perilesional and extralesional areas. ROIs were exposed to ultraviolet light, releasing the indexing oligos and collecting them in a 96-well plate for subsequent processing and sequencing, as described (72). Raw count, third quartile (Q3)− normalized count data of target genes from ROIs were provided by the vendor, which were used as input to downstream analyses. Additional details can be found in the Supplementary Materials.

Immunofluorescence

Immunofluorescence staining of mouse lungs were performed as previously described (39). Briefly, lungs from WT or Cx3cr1GFP+ mice were perfused with 5 ml PBS through the right ventricle then with 5 ml paraformaldehyde (PFA) 4% (Thermo Fisher) in PBS, and lungs were collected. Lungs were fixed for 4 h in 4% PAF at 4 °C, then cryoprotected overnight in 30% sucrose (VWR) in PBS at 4 °C, followed by embedding in optimal cutting temperature compound (OCT) (VWR) and stored at −80 °C. Seven-μm-thick sections were cut and left in a methanol 100% (Merck) bath at −20 °C for 20 minutes prior to be stained. Additional details can be found in the Supplementary Materials. All images were acquired on an LSM 980 with Airyscan 2 inverted confocal microscope (Zeiss) using a LD C-Apochromat ×40/1.1 W objective and Zen Black software. Additional details can be found in the Supplementary Materials.

Ex vivo experiments

For ex vivo stimulation and co-culture experiments, single cell suspensions isolated from IAV-infected lungs of WT or Ly6gtdTom mice at day 10 post-IAV were enriched in CD11b+ cells by a magnetic-activated cell sorting (MACS) using CD11b MicroBeads (Myltenyi). Cells were then stained and FACS-sorted using the gating strategy shown in Fig. 1A. After sorting, cells were counted, spinned down, and either directly added to the co-culture with MLE-12 cells, or seeded in 96 wells at a concentration 5.104 cells/well in complete RPMI (ThermoFischer), containing 1mM sodium pyruvate, 1% vol/vol MEM non-essential amino acids, 50 U ml−1 Penicillin-Streptomycin and 10% vol/vol FBS. For stimulation experiments, recombinant mouse GM-CSF (20 ng ml-1, Peprotech), mouse M-CSF (20 ng ml-1, Peprotech), mouse IL-4 (20 ng ml-1, Peprotech) or mouse IL-13 (20 ng ml-1, Peprotech) were added. When required, Cre-ERT2 activation was achieved by adding 0.02 mg ml-1 of 4-hydroxytamoxifen (Sigma-Aldrich). After 18 hours of culture, cell supernatants were collected (conditioned medium, CM) and cells were harvested for flow cytometry phenotyping.

Murine lung epithelial (MLE)-12 cells (ATCC, CVCL_3751) were used. Scratch Wound Assay were performed using IncucyteS3 (Sartorius). MLE-12 cells were seeded in 96-well (Sartorius) at density of 4.104 cells/well and incubated 24 hours in DMEM/F12 medium. An open wound area was created in the cell monolayer using the IncuCyte ® Wound Maker tool, and subsequently co-cultured with FCAS-sorted cells or incubated with CM from unpulsed or IL-4/IL-13-pulsed Ly6G+ Macs.

Additional details can be found in the Supplementary Materials.

Adoptive transfer of Ly6G+ Macs in vivo

Ly6G+ Macs were isolated from the lungs of CD45.2 WT mice at day 10 poist-IAV. Lung single cell suspensions were first enriched in CD11b+ cells by MACS using CD11b MicroBeads (Miltenyi Biotec) and were FACS-sorted using a Sony MA900. Four hundred thousand (4 x 105) Ly6G+ Macs were resuspended in 50 μl sterile PBS and were instilled i.t. in lightly isoflurane-anesthetized Maf/MafbMyeloKO mice at days 8, 11, 13, 15 post-IAV. Control Maf/MafbMyeloKO mice and WT mice received 50 μl PBS as vehicle.

Human BALFs

The use of human BALF cells was approved in 2022 by the Ethics Reviewing Board of the University Hospital of Liege, Belgium (ref. 2022/159). The characteristics of the patients are summarized in Table S1. Human BALFs were fixed directly after collection for storage and scRNA-seq analyses.

Statistical analysis

Graphs were prepared with GraphPad Prism 9 (GraphPad software) or R Bioconductor (3.5.1) (73). Data distribution was assumed to be normal when parametric tests were performed. Data from independent experiments were pooled for analysis in each data panel, unless otherwise indicated. No data were excluded from the analyses. Statistical analyses were performed with Prism 9 (GraphPad software), and with R Bioconductor (3.5.1) (73) and Seurat (76) for scRNA-seq data, respectively. The statistical analyses performed for each experiment are indicated in the respective figure legends. We considered a P value lower than 0.05 to be significant (*, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001; ns, not significant).

Additional sections and details about Materials and Methods can be found in the Supplemental Materials.

Supplementary Material

Supplementary Materials

Acknowledgments

We thank all members of the Immunophysiology laboratory (GIGA Institute, Liège, Belgium) for discussions; N. Maloujahmoum (Metastasis Research Laboratory, GIGA Institute) for expert technical assistance; S. Ormenese, R. Stefan, A. Hego, G. Lefevre and C. Vanwinge from the GIGA In Vitro Imaging Platform; P. Drion, G. Lambert, L., B. Remy and all staff members from the GIGA Mouse facility and Transgenics Platform; W. Coppieters, L. Karim, M. Deckers, A. Mayer, A. Lavergne and members from the GIGA Genomics Platform; R. Fares, I. Sbai and A. Lio for their excellent administrative support.

Funding

Supported by:

FRFS-Welbio Advanced Grant WELBIO-CR-2022A-10 (TM)

ERC Starting Grant ERC StG 2018 IM-ID: 801823 (TM)

Baillet Latour Biomedical Fund (TM)

Research Project of the F.R.S.-FNRS T015021F (TM)

Research Concerted Action of the Fédération Wallonie-Bruxelles (TM)

Data and materials availability

Ms4a3Cre mice are used under the terms of a material transfer agreement (MTA) signed between Singapore Immunology Network (SIgN), A*STAR Research Entities (Singapore) and the University of Liege (Belgium). Ly6gCreERT2 mice are used under the terms of an MTA signed between Centro Nacional de Investigaciones Cardiovasculares Carlos III (CNIC, Madrid, Spain), and the University of Liege (Belgium). Mouse and human single-cell RNA-seq data have been deposited at GEO and are available under GEO accession GSE244765. All original codes have been deposited at Zenodo and are available via this link: https://zenodo.org/records/11354523. Requests for Mabfl/fl mice should be addressed to t.marichal@uliege.be and will be shared via an MTA.

One-Sentence Summary

Short-lived perilesional Ly6G+ macrophages license alveolar epithelial regeneration after viral-triggered lung injury.

Fig. 1 Ly6G+ Macs culminate during the early recovery phase post-IAV infection.

(A) Representative flow cytometry gating strategy showing live CD45+CD11b+Ly6G+CD64− neutrophils (Neu), CD45+CD11b+Ly6G+CD64+ macrophages (Ly6G+ Mac), CD45+Ly6G−CD11c+SiglecF+ alveolar macrophages (AM), CD45+Ly6G−SiglecF−F4/80+CD11b+Ly6C+CD64− monocytes (Ly6C+ Mo), CD45+Ly6G−SiglecF−F4/80+CD11b+Ly6C−CD64− monocytes (Ly6C− Mo), CD45+Ly6G−SiglecF−F4/80+CD11b+Ly6C+CD64+ inflammatory monocytes (iMo) and CD45+Ly6G−SiglecF−F4/80+CD11b+Ly6C−CD64+ IM-like cells in lungs of C57BL/6 wild-type (WT) mice at day 10 post-IAV. (B) Time course of absolute numbers of Neu, AM, Ly6C− Mo, Ly6C+ Mo, iMo and IM-like cells quantified by flow cytometry at days 0, 5, 10, 15 and 20 post-IAV in WT mice. (C) Representative contour plots of CD64 and CD11b expression within lung CD45+CD11b+Ly6G+ cells at day 10 p.i. in mock-infected or IAV-infected WT mice. (D) Time course of absolute numbers of Ly6G+ Macs quantified by flow cytometry, as in (B). (E) Percentage of Neu and Ly6G+ Macs within Ly6G+CD11b+ cells quantified by flow cytometry in the blood and lungs of WT mice at day 10 post-IAV. (F) Photographs of Neu, Ly6G+ Macs, IM-like cells and iMo sorted by FACS from IAV-infected WT mice at day 10 p.i.. Pictures are representative of 1 of 3 independent sorting experiments, each giving similar results. (G) Representative histograms of CXCR4, MHC-II, CD101, CD319 and CD177 expression in the indicated myeloid cell populations, quantified by flow cytometry at day 10 post-IAV in WT mice. (H) Quantification of expression of the indicated markers, as in (G). (B,D) Data show mean (centerline) ± SEM (colored area) and are pooled from 2-3 independent experiments (n=6 mice per time point); (E,H) mean + SEM and are pooled from 2 independent experiments (n=5-6 mice). (B,D,H) P values were calculated using a one-way ANOVA with Dunnett’s post hoc tests. *, P<0.05; ***, P<0.001; ****, P<0.0001. FMO, fluorescence minus one; ns, not significant; p.i., post-infection. (F) Scale bars: 5 μm.

Fig. 2 Ly6G+ Macs are transcriptionally distinct from other lung myeloid cells at day 10 post-IAV.

(A) UMAP plots of scRNA-seq data depicting the transcriptional identity of FACS-sorted lung live CD45+F4/80+ and/or CD11b+ cells from mock- or IAV-infected WT mice 10 days p.i. (pooled from 5 mice per conditions), merged with a published dataset of steady-state lung monocytes and IMs (29). (B) Frequency of each cluster within each experimental condition, as in (A). (C) Heatmap depicting the single cell expression of the most upregulated genes within each cluster. (D) Expression of the indicated genes within each cluster, as depicted by violin plots (height: expression; width: abundance of cells). (E) Representative histograms of intracellular Arg-1 and osteopontin expression in the indicated lung myeloid cell populations, quantified by flow cytometry at day 10 post-IAV in WT mice. (F) Quantification of Arg-1 and osteopontin expression, as in (E). (F) Data show mean + SEM and are pooled from 2 independent experiments (n=6 mice). P values were calculated using (D) a Wilcoxon rank sum test and compare C2 vs. all other clusters or (F) a one-way ANOVA with Dunnett’s post hoc tests. *, P<0.05; **, P<0.01; ***, P<0.001; ****, P<0.0001. FMO, fluorescence minus one; p.i., post-infection.

Fig. 3 IAV-triggered Ly6G+ Macs are recruited from classical monocytes and are short-lived.

(A) UMAP plot depicting the transcriptional identity and cell trajectories (top), and pseudotime trajectory values (below) of lung Ly6C+ Mo, iMo, Ly6G+ Mac, dying Mac, CD206− IM and CD206+ IM, as in Fig. 2A, evaluated by Slingshot trajectory analyses. (B) Heatmap plot depicting the differentially expressed genes along pseudotime evaluated by tradeSeq in the trajectory starting from Ly6C+ Mo and ending either in IM or in Ly6G+ Mac. (C) Representative histograms of tdTomato (left) and GFP (right) expression in the indicated myeloid cell populations, quantified by flow cytometry at day 10 post-IAV in Ms4a3tdTom and Cx3Cr1GFP mice, respectively. (D) Quantification of tdTomato+ cells (left) and GFP expression (right), as in (C). (E) Representative tdTomato and CD45.1 contour plots and (F) bar graph showing % of host, donor Ccr2−/− and donor Ms4a3tdTom+ chimerism in the indicated cell populations from lethally-irradiated CD45.1/CD45.2 mice reconstituted with a 1:1 mix of CD45.2 Ccr2−/− and Ms4a3tdTom+ BM cells, infected with IAV 4 weeks later and evaluated at day 10 post-IAV. (G) Time course of absolute numbers of EdU+ Ly6G+ Macs and EdU+ IM-like cells quantified by flow cytometry at days 7, 10, 14 and 17 post-IAV in EdU-pulsed WT mice at day 7 post-IAV. (H) Pie chart representation of the mean frequency of Annexin V and PI negative and/or positive fractions within lung Ly6C+ Mo, iMo and Ly6G+ Macs, quantified at day 10 post-IAV in WT mice. (I) Representative Ly6G and FSC contour plots and (J) bar graph showing % of Ly6G+ cells within lung iMo sorted from WT mice at day 10 post-IAV and cultured 18 hours ex vivo with vehicle, M-CSF or GM-CSF. (K) Representative confocal microscopy pictures and (L) representative flow cytometry histograms of tdTomato expression within lung iMo sorted from Ly6gtdTom mice at day 10 post-IAV and treated ex vivo with tamoxifen and GM-CSF or vehicle for 18 hours. (M) Bar graph showing donor Csf2ra−/− chimerism relative to donor Csf2ra+/+ chimerism in the indicated cell populations from thorax-protected, lethally-irradiated CD45.1/CD45.2 mice reconstituted with a 1:1 mix of CD45.1 Csf2ra−/− and CD45.2 Csf2ra+/+ BM cells, infected with IAV 4 weeks later and evaluated at day 10 post-IAV. (N) Representative histograms and (O) quantification of Arg-1+ cells (%) in Ly6G+ Macs from donor Csf2ra+/+ and Csf2ra−/− BM cells, as in (M). Data show (D,F,J,M,O) mean + SEM and (D,F) are representative of 1 of 3 independent experiments (n=3-4 mice), (J) are pooled from 3 independent sorting experiments, each dot representing one biological replicate, (M,O) are pooled from 2 independent experiments (n=10 mice); (G) mean (centerline) ± SEM (colored area) and are pooled from 2 independent experiments (n=6 mice per time point). P values compare CD45.2 donor Ccr2−/− chimerism in (F). P values were calculated using (D) a one-way ANOVA with Dunnett’s post hoc tests, (F) a two-way ANOVA with Tukey’s post hoc tests, (G,J) a one-way ANOVA with Tukey’s post hoc tests, (M) a two-way ANOVA with Sidak’s post hoc tests, (O) a two-tailed Student’s t test. *, P<0.05; ***, P<0.001; ****, P<0.0001. FMO, fluorescence minus one; ns, not significant.

Fig. 4 Ly6G+ Macs exhibit an atypical morphology and possess distinct metabolic, phagocytic and efferocytic capabilities.

(A) GSEA analyses of Ly6G+ Mac (C2) profile compared to other clusters using KEGG, Cellular Components and Biological Process gene sets. The Normalized Enrichment Score (NES), False Discovery Rate (FDR) and the size of the gene set are shown for each process. (B) Representative transmission electron microscopy pictures of Neu, Ly6G+ Mac, IM-like cells and iMo FACS-sorted from lungs of WT mice at day 10 post-IAV. (C) Extracellular acidification rate (ECAR) of FACS-sorted Ly6G+ Mac, Neu and IM-like cells, as in (B), quantified at baseline and under stress over time using a Seahorse assay. (D) Oxygen consumption rate (OCR) of Ly6G+ Mac, Neu and IM-like cells, as in (C). (E) ECAR and OCR of Ly6G+ Macs, Neu and IM-like cells, as in (C,D). (F) Representative histograms of E. coli-FITC signal in the indicated myeloid cell populations, quantified by flow cytometry at day 10 post-IAV and 3 hours after i.t. injection of E. coli-FITC particles. Non-injected mice were used as controls (grey line). (G) Quantification of E. coli-FITC+ cells, as in (F). (H) Representative tdTomato and CD11b contour plots of and (I) bar graph showing % of tdTomato+ cells in the indicated Cx3cr1GFP+ donor cell populations from lethally-irradiated CD45.2 WT mice reconstituted with a 1:1 mix of CD45.2 Cx3cr1GFP+ and Ms4a3tdTom+ BM cells, infected with IAV 4 weeks later and evaluated at day 10 post-IAV. Data show (C,D) mean ± SEM and are representative of 1 of 3 independent experiments, each giving similar results; (G,I) mean + SEM and are pooled from 2 independent experiments (n=6 mice). P values (C,D) compare Ly6G+ Macs vs. IM-like cells or Neu and were calculated using a two-way ANOVA with Bonferroni’s post hoc tests; (G,I) were calculated using a one-way ANOVA with Dunnett’s post hoc tests. **, P<0.01; ****, P<0.0001. (B) Scale bars: 2 μm.

Fig. 5 Ly6G+ Macs populate the alveoli of perilesional regenerating areas.

(A) Representative confocal microscopy picture of lung sections from Cx3cr1GFP mice at day 10 post-IAV, with Ly6G+ Macs identified as Ly6G+Cx3cr1GFP+ cells. (B) Representative in situ electron microscopy picture of Ly6G+ Macs in the vicinity of AT2 and AT1 cells, identified on lung sections from WT mice at day 10 post-IAV. (C) Representative examples of regions of interest (ROIs) selected on lung sections from mock- or IAV-infected WT mice at day 10 post-IAV stained with anti-Ly6G and anti-CD68 antibodies. (D) Unsupervised Principal Component (PC) analysis of the ROIs analyzed by DSP. Percentages indicate the variability explained by each component. (E) Ly6G+ Mac signature score within control, extralesional, perilesional and intralesional ROIs, as depicted by violin plots (height: scores; width: abundance of cells). (F) GSEA analysis of perilesional ROIs compared to intralesional ROIs using Cellular Components, Molecular Function and Biological Process gene sets. The Normalized Enrichment Score (NES), False Discovery Rate (FDR) and the size of the gene set are shown for each process. (G) Correlation of Ly6G+ Mac score with primed AT2 (top) and DATP (bottom) scores of the ROIs. (H) Representative picture of perilesional Ly6G+ Macs (i.e., Ly6G+MHC-II+ cells), pSPC+ AT2 and podoplanin+ AT1 identified by confocal microscopy on lung sections from WT mice at day 10 post-IAV. (A,B,H) Pictures are representative of one of 6 mice, each giving similar results. (E) P values were calculated using a one-way ANOVA with Tukey’s post hoc tests. (G) The correlation analysis used was parametric Pearson correlation coefficient. **, P<0.01; ***, P<0.001. Scale bars: (A) 15, (B) 5, (C) 100, (D) 10 μm.

Fig. 6 C-Maf/MafB-dependent Ly6G+ Macs promote euplastic alveolar epithelial regeneration.

(A) Heatmap depicting predicted activities of c-Maf and MafB across lung myeloid cells post-IAV, evaluated by SCENIC analysis of the scRNA-seq data, as in Fig. 2A. (B) Representative histograms of intracellular c-Maf and MafB expression in the indicated lung myeloid cell populations, quantified at day 10 post-IAV. (C) Quantification of expression of intracellular c-Maf and MafB, as in (B). (D) Expression of Maf and Mafb within control, extralesional, perilesional and intralesional ROIs, as depicted by violin plots (height: normalized counts; width: abundance of cells). (E) Absolute numbers of the indicated lung myeloid cell populations, quantified at day 10 post-IAV in control and Maf/MafbMyeloKO mice. (F) Time course of relative lung NS1 RNA expression, assessed by RT-qPCR at days 0, 7 and 10 post-IAV in control and Maf/MafbMyeloKO mice. (G) Time course of weight, expressed as the % of the original weight at day 0 and assessed at days 0, 7, 10, 15, 20 post-IAV in control and Maf/MafbMyeloKO mice. (H-I) Representative (H) Hematoxylin & Eosin and (I) Periodic Acid Schiff (PAS, bottom) pictures of lung sections of control and Maf/MafbMyeloKO mice at day 20 post-IAV. Pictures are representative of 1 of 7 mice analyzed. (J) Percentage of PAS+ cells in lung lesional areas of control and Maf/MafbMyeloKO mice at day 20 post-IAV. (K) Representative pSPC and Pdpn contour plots of CD45-CD31-EpCam+ cells in mock- or IAV-infected mice at day 20 post-IAV. (L) Absolute numbers of pSPC+Pdpn− AT2, pSPC−Pdpn+ AT1 and pSPC+Pdpn+ regenerating AT2 (reg AT2), quantified as in (K). (M) Time course of weight, expressed as the % of the original weight at day 8 and assessed at days 8, 11, 13, 15 and 20 post-IAV in control mice, in Maf/MafbMyeloKO mice instilled i.t. at days 8, 11, 13 and 15 post-IAV with PBS or with 3 x 105 Ly6G+ Macs isolated from WT mice at day 10 post-IAV. (N) Absolute numbers of pSPC+Pdpn- AT2, pSPC-Pdpn+ AT1 and pSPC+Pdpn+ regenerating AT2 (reg AT2), quantified at day 20 post-IAV, as in (M). (C,E,J,L,N) Data show mean + SEM and are pooled from 2 independent experiments (C:n=5 mice; E,N: n=6 mice/group; J: n=7 mice/group; L: n=8 mice/group). (F,G,M) Data show mean (centerline) ± SEM (colored area) and are pooled from 2-3 independent experiments (F: n=10 mice/group; G: n=7 mice/group; M: n=6 mice/group). P values were calculated using (A) a Wilcoxon rank sum test, (C) a one-way ANOVA with Dunnett’s post hoc tests, (D,L) a one-way ANOVA with Tukey’s post hoc tests, (E,F,G,M) a two-way ANOVA with Sidak’s post hoc tests, (J) a two-tailed Student’s t test, *, P<0.05; **, P<0.01; ***, P<0.001; ****, P<0.0001. FMO, fluorescence minus one; ns, not significant; p.i., post-infection. Scale bars: (H, top) 2 mm, (H, bottom; I) 200 μm.

Fig. 7 IL-4R-dependent Ly6G+ Macs release soluble factors that improve alveolar regeneration from AT2 cells.

(A) Type 2 signature score within control, extralesional, perilesional and intralesional ROIs, as depicted by violin plots (height: scores; width: abundance of cells). (B) Representative histograms of IL-4R expression in the indicated lung myeloid cell populations, quantified at day 10 post-IAV. (C) Quantification of IL-4R expression, as in (H). (D) Bar graph showing donor Il4ra−/− chimerism relative to donor Ms4a3tdTom+ chimerism in the indicated cell populations from lethally-irradiated CD45.1 mice reconstituted with a 1:1 mix of CD45.2 Il4ra−/− and Ms4a3tdTom+ BM cells, infected with IAV 4 weeks later and evaluated at day 10 post-IAV. (E) Representative histograms and (F) quantification of Arg-1+ cells (%) in Ly6G+ Macs from donor Il4ra−/− or Ms4a3tdTom+ BM cells, as in (F). (G) Time course of weight, expressed as the % of the original weight at day 0 and assessed at days 0, 7, 10, 15, 20 post-IAV in lethally-irradiated CD45.1/CD45.2 WT mice reconstituted with CD45.2 Il4ra−/− BM cells (Il4ra−/− BM -> WT) or CD45.2 WT BM cells (WT BM -> WT) and infected with IAV 4 weeks later. (H) Cell confluence of AT2 cells (MLE-12) quantified 12 hours after a standardized scratch by live cell analysis when AT2 cells were co-cultured in the presence of Neu, IM-like cells, iMo or Ly6G+ Macs isolated from the lungs of WT mice at day 10 post-IAV. (I) Representative picture of cell confluence of AT2 cells at 0 and 12 hours post-scratch when AT2 cells were cultured with IL-4/13 or with conditioned medium (CM) of unpulsed or IL-4/13-pulsed Ly6G+ Macs isolated from the lungs of WT mice at day 10 post-IAV. (J) Bar graph of cell confluence of AT2 cells quantified 0, 12 and 18 hours after scratch, as in (K). (K) Heatmap showing the proteome profiling of CM of vehicle and IL-4/13-treated Ly6G+ Macs isolated from the lungs of WT mice at day 10 post-IAV. (C,D,F,H,J) Data show mean + SEM and (C,D,F) are pooled from 2 independent experiments (C: n=6 mice; D,F: n= 8 mice); (H,J) are pooled from 3 independent sorting experiments. (G) Data show mean (centerline) ± SEM (colored area) and are pooled from 2 independent experiments (n=6 mice per group). P values were calculated using (C) a one-way ANOVA with Dunnett’s post hoc tests, (D,G) a two-way ANOVA with Sidak’s post hoc tests, (F) a two-tailed Student’s t test, (H) a one-way or (J) a two-way ANOVA with Tukey’s post hoc tests. *, P<0.05; **, P<0.01; ****, P<0.0001. FMO, fluorescence minus one; ns, not significant; p.i., post-infection.

Fig. 8 Ly6G+ Macs are triggered by other insults and have a human counterpart.

(A) Time course of weight, expressed as the % of the original weight at day 0 and assessed at days 0, 3, 7, 10, 14 and 18 post-injection in C57BL/6 WT mice instilled i.t. with bleomycin (bleo). (B) Time course of absolute numbers of lung Ly6G+ Macs quantified at days 0, 5, 10, 14 and 18 post-bleo in WT mice. (C) Representative histograms of intracellular Arg-1 and CXCR4 expression in the indicated lung myeloid cell populations, quantified at day 14 post-bleo. (D) Quantification of intracellular Arg-1 and CXCR4 expression, as in (C). (E) Time course of absolute numbers of pSPC-Pdpn+ AT1 and pSPC+Pdpn- AT2, quantified at days 0, 7, 10, 14 and 18 post-bleo. (F) Amount of alanine aminotransferase (ALT) in the plasma of WT mice injected i.p. with acetaminophen at days 0, 1, 2, 3 and 4 post-injection. (G) Time course of absolute numbers of liver Ly6G+ Macs quantified at days 0, 1, 2, 3 and 4 post-acetominophen. (H) Representative histograms of intracellular Arg-1 and CXCR4 expression in liver neutrophils (Neu) and Ly6G+ Macs, quantified at day 1 post-acetominophen. (I) Quantification of intracellular Arg-1 and CXCR4 expression, as in (H). (J) UMAP plot depicting the transcriptional identity of human BALF cells collected from 7 patients suspected of pneumonia and analyzed by scRNA-seq. Annotations of cell clusters are shown. (K) Representation of each patient within each cluster, shown as frequency. (L) UMAP feature plot, as in (J), according to the Ly6G+ Mac signature score. The score level in cluster C9 is shown for each patient. (M) Ly6G+ Mac signature score of single cells within each cluster, as depicted by violin plots (height: score; width: abundance of cells). (N) Heatmap depicting predicted activities of MAFB and MAF across BALF cell populations, evaluated by SCENIC analysis of the scRNA-seq data shown in (J). (A,B,E,F,G) Data show mean (centerline) ± SEM (colored area) and are pooled from 2 independent experiments (n=5-6 mice). (D,I) Data show mean + SEM and are pooled from 2 independent experiments (n=5-7 mice). P values were calculated using (A) a two-way ANOVA, (B,D,E,F,G) a one-way ANOVA with Dunnett’s post hoc tests, (I) a two-tailed Student’s t test or (M,N) a Wilcoxon rank sum test. (M) P values compare C9 vs. all other clusters. *, P<0.05; **, P<0.01; ***, P<0.001; ****, P<0.0001. FMO, fluorescence minus one.

Author contributions

T.M. and C.Ra. conceived, supervised and secured funding for the project; C.Ru., C.Ra. and T.M. designed the experiments; C.Ru. did most of the experiments, compiled the data and prepared the figures. J.A. performed all the bioinformatic analyses. P.M., M.M., D.V. and C.Ra. helped with in vivo-related experiments. C. de m. and M-M. G. were implicated in experiments related to the histopathological analyses of myeloid-restricted Maf and Mafb-deficient mice. P.J., J.G. and F.S. were involved in human BALF collection and processing. M.D. and B.M. performed viral titrations and helped with viral RNA quantifications. M.T. performed electron microscopy experiments and analyses. F.B., C.S., A.H., F.G and B.G.D. were implicated in experiments involving BM chimeric mice. A.B. performed extracellular flux analyses. C.Ru, C.Ra and T.M. wrote the manuscript. All authors provided feedback on the original and revised manuscript.

Competing interests

Authors declare that they have no competing interests.
==== Refs
References and notes

1 Kalil AC Thomas PG Influenza virus-related critical illness: pathophysiology and epidemiology Crit Care 2019 23 258 31324202
2 Herold S Becker C Ridge KM Budinger GRS Influenza virus-induced lung injury: pathogenesis and implications for treatment Eur Respir J 2015 45 1463 1478 25792631
3 Wei X Narasimhan H Zhu B Sun J Host Recovery from Respiratory Viral Infection Annu Rev Immunol 2023 41 277 300 36716750
4 Burden of disease - Global Influenza Programme https://www.who.int/teams/global-influenza-programme/surveillance-and-monitoring/burden-of-disease
5 Iuliano AD Roguski KM Chang HH Muscatello DJ Palekar R Tempia S Cohen C Gran JM Schanzer D Cowling BJ Wu P Global Seasonal Influenza-associated Mortality Collaborator Network, Estimates of global seasonal influenza-associated respiratory mortality: a modelling study Lancet 2018 391 1285 1300 29248255
6 Bost P Giladi A Liu Y Bendjelal Y Xu G David E Blecher-Gonen R Cohen M Medaglia C Li H Deczkowska A Host-Viral Infection Maps Reveal Signatures of Severe COVID-19 Patients Cell 2020 181 1475 1488 e12 32479746
7 Merad M Martin JC Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophages Nat Rev Immunol 2020 20 355 362 32376901
8 Guilliams M Mildner A Yona S Developmental and Functional Heterogeneity of Monocytes Immunity 2018 49 595 613 30332628
9 Ruscitti C Radermecker C Marichal T Journey of monocytes and macrophages upon influenza A virus infection Curr Opin Virol 2024 66 101409 38564993
10 Rodriguez-Rodriguez L Gillet L Machiels B Shaping of the alveolar landscape by respiratory infections and long-term consequences for lung immunity Front Immunol 2023 14
11 Machiels B Dourcy M Xiao X Javaux J Mesnil C Sabatel C Desmecht D Lallemand F Martinive P Hammad H Guilliams M A gammaherpesvirus provides protection against allergic asthma by inducing the replacement of resident alveolar macrophages with regulatory monocytes Nat Immunol 2017 18 1310 1320 29035391
12 Aegerter H Kulikauskaite J Crotta S Patel H Kelly G Hessel EM Mack M Beinke S Wack, Influenza-induced monocyte-derived alveolar macrophages confer prolonged antibacterial protection Nat Immunol 2020 21 145 157 31932810
13 Aegerter H Lambrecht BN Jakubzick CV Biology of lung macrophages in health and disease Immunity 2022 55 1564 1580 36103853
14 Guilliams M Svedberg FR Does tissue imprinting restrict macrophage plasticity? Nat Immunol 2021 22 118 127 33462453
15 Li F Piattini F Pohlmeier L Feng Q Rehrauer H Kopf M Monocyte-derived alveolar macrophages autonomously determine severe outcome of respiratory viral infection Sci Immunol 2022 7 eabj5761 35776802
16 Lin KL Suzuki Y Nakano H Ramsburg E Gunn MD CCR2+ monocyte-derived dendritic cells and exudate macrophages produce influenza-induced pulmonary immune pathology and mortality J Immunol 2008 180 2562 2572 18250467
17 Chen ST Park MD Del Valle DM Buckup M Tabachnikova A Thompson RC Simons NW Mouskas K Lee B Geanon D D’Souza D A shift in lung macrophage composition is associated with COVID-19 severity and recovery Sci Transl Med 2022 14 eabn5168 36103512
18 Park MD Silvin A Ginhoux F Merad M Macrophages in health and disease Cell 2022 185 4259 4279 36368305
19 Wynn TA Vannella KM Macrophages in Tissue Repair, Regeneration, and Fibrosis Immunity 2016 44 450 62 26982353
20 Nobs SP Kopf M Tissue-resident macrophages: guardians of organ homeostasis Trends Immunol 2021 42 495 507 33972166
21 Ginhoux F Guilliams M Tissue-Resident Macrophage Ontogeny and Homeostasis Immunity 2016 44 439 49 26982352
22 Huang F-F Barnes PF Feng Y Donis R Chroneos ZC Idell S Allen T Perez DR Whitsett JA Dunussi-Joannopoulos K Shams H GM-CSF in the lung protects against lethal influenza infection Am J Respir Crit Care Med 2011 184 259 268 21474645
23 Schneider C Nobs SP Heer AK Kurrer M Klinke G van Rooijen N Vogel J Kopf M Alveolar macrophages are essential for protection from respiratory failure and associated morbidity following influenza virus infection PLoS Pathog 2014 10 e1004053 24699679
24 Ural BB Yeung ST Damani-Yokota P Devlin JC de Vries M Vera-Licona P Samji T Sawai CM Jang G Perez OA Pham Q Identification of a nerve-associated, lung-resident interstitial macrophage subset with distinct localization and immunoregulatory properties Sci Immunol 2020 5 eaax8756 32220976
25 Bain CC MacDonald AS The impact of the lung environment on macrophage development, activation and function: diversity in the face of adversity Mucosal Immunol 2022 15 223 234 35017701
26 Chen J Zhong M-C Guo H Davidson D Mishel S Lu Y Rhee I Pérez-Quintero L-A Zhang S Cruz-Munoz M-E Wu N SLAMF7 is critical for phagocytosis of haematopoietic tumour cells via Mac-1 integrin Nature 2017 544 493 497 28424516
27 Wu Y Wang Q Li M Lao J Tang H Ming S Wu M Gong S Li L Liu L Huang X SLAMF7 regulates the inflammatory response in macrophages during polymicrobial sepsis J Clin Invest 2023 133 e150224 36749634
28 Zheng GX Terry JM Belgrader P Ryvkin P Bent ZW Wilson R Ziraldo SB Wheeler TD McDermott GP Zhu J Gregory MT Massively parallel digital transcriptional profiling of single cells Nat Commun 2017 8 14049 28091601
29 Vanneste D Bai Q Hasan S Peng W Pirottin D Schyns J Maréchal P Ruscitti C Meunier M Liu Z Legrand C Fievez L MafB-restricted local monocyte proliferation precedes lung interstitial macrophage differentiation Nat Immunol 2023 24 827 840 36928411
30 Sabatel C Radermecker C Fievez L Paulissen G Chakarov S Fernandes C Olivier S Toussaint M Pirottin D Xiao X Quatresooz P Exposure to Bacterial CpG DNA Protects from Airway Allergic Inflammation by Expanding Regulatory Lung Interstitial Macrophages Immunity 2017 46 457 473 28329706
31 Liu Z Gu Y Chakarov S Bleriot C Kwok I Chen X Shin A Huang W Dress RJ Dutertre C-A Schlitzer A Fate Mapping via Ms4a3-Expression History Traces Monocyte-Derived Cells Cell 2019 178 1509 1525 e19 31491389
32 Street K Risso D Fletcher RB Das D Ngai J Yosef N Purdom E Dudoit S Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics BMC Genomics 2018 19 477 29914354
33 Yona S Kim KW Wolf Y Mildner A Varol D Breker M Strauss-Ayali D Viukov S Guilliams M Misharin A Hume DA Fate mapping reveals origins and dynamics of monocytes and tissue macrophages under homeostasis Immunity 2013 38 79 91 23273845
34 Serbina NV Pamer EG Monocyte emigration from bone marrow during bacterial infection requires signals mediated by chemokine receptor CCR2 Nat Immunol 2006 7 311 7 16462739
35 Satoh T Nakagawa K Sugihara F Kuwahara R Ashihara M Yamane F Minowa Y Fukushima K Ebina I Yoshioka Y Kumanogoh A Identification of an atypical monocyte and committed progenitor involved in fibrosis Nature 2017 541 96 101 28002407
36 Aibar S González-Blas CB Moerman T Huynh-Thu VA Imrichova H Hulselmans G Rambow F Marine J-C Geurts P Aerts J van den Oord J SCENIC: single-cell regulatory network inference and clustering Nat Methods 2017 14 1083 1086 28991892
37 Ballesteros I Rubio-Ponce A Genua M Lusito E Kwok I Fernández-Calvo G Khoyratty TE van Grinsven E González-Hernández S Nicolás-Ávila JÁ Vicanolo T Co-option of Neutrophil Fates by Tissue Environments Cell 2020 183 1282 1297 e18 33098771
38 Choi J Park J-E Tsagkogeorga G Yanagita M Koo B-K Han N Lee J-H Inflammatory Signals Induce AT2 Cell-Derived Damage-Associated Transient Progenitors that Mediate Alveolar Regeneration Cell Stem Cell 2020 27 366 382 e7 32750316
39 Vanneste D Bai Q Hasan S Peng W Pirottin D Schyns J Maréchal P Ruscitti C Meunier M Liu Z Legrand C MafB-restricted local monocyte proliferation precedes lung interstitial macrophage differentiation Nat Immunol 2023 24 827 840 36928411
40 Adamson IY Bowden DH The type 2 cell as progenitor of alveolar epithelial regeneration. A cytodynamic study in mice after exposure to oxygen Lab Invest 1974 30 35 42 4812806
41 Planer JD Morrisey EE After the Storm: Regeneration, Repair, and Reestablishment of Homeostasis Between the Alveolar Epithelium and Innate Immune System Following Viral Lung Injury Annu Rev Pathol 2023 18 337 359 36270292
42 Minutti CM Jackson-Jones LH García-Fojeda B Knipper JA Sutherland TE Logan N Ringqvist E Guillamat-Prats R Ferenbach DA Artigas A Stamme C Local amplifiers of IL-4Rα-mediated macrophage activation promote repair in lung and liver Science 2017 356 1076 1080 28495878
43 Bosurgi L Cao YG Cabeza-Cabrerizo M Tucci A Hughes LD Kong Y Weinstein JS Licona-Limon P Schmid ET Pelorosso F Gagliani N Macrophage function in tissue repair and remodeling requires IL-4 or IL-13 with apoptotic cells Science 2017 356 1072 1076 28495875
44 Basil MC Katzen J Engler AE Guo M Herriges MJ Kathiriya JJ Windmueller R Ysasi B Zacharias WJ Chapman HA Kotton DN The Cellular and Physiological Basis for Lung Repair and Regeneration: Past, Present, and Future Cell Stem Cell 2020 26 482 502 32243808
45 Zacharias WJ Frank DB Zepp JA Morley MP Alkhaleel FA Kong J Zhou S Cantu E Morrisey EE Regeneration of the lung alveolus by an evolutionarily conserved epithelial progenitor Nature 2018 555 251 255 29489752
46 Nabhan AN Brownfield DG Harbury PB Krasnow MA Desai TJ Single-cell Wnt signaling niches maintain stemness of alveolar type 2 cells Science 2018 359 1118 1123 29420258
47 Liberti DC Kremp MM Liberti WA Penkala IJ Li S Zhou S Morrisey EE Alveolar epithelial cell fate is maintained in a spatially restricted manner to promote lung regeneration after acute injury Cell Rep 2021 35 109092 33979629
48 Misharin AV Morales-Nebreda L Reyfman PA Cuda CM Walter JM McQuattie-Pimentel AC Chen C-I Anekalla KR Joshi N Williams KJN Abdala-Valencia H Monocyte-derived alveolar macrophages drive lung fibrosis and persist in the lung over the life span J Exp Med 2017 214 2387 2404 28694385
49 Lee PY Wang J-X Parisini E Dascher CC Nigrovic PA Ly6 family proteins in neutrophil biology J Leukoc Biol 2013 94 585 594 23543767
50 Cakarova L Marsh LM Wilhelm J Mayer K Grimminger F Seeger W Lohmeyer J Herold S Macrophage tumor necrosis factor-alpha induces epithelial expression of granulocyte-macrophage colony-stimulating factor: impact on alveolar epithelial repair Am J Respir Crit Care Med 2009 180 521 532 19590023
51 Lavin Y Winter D Blecher-Gonen R David E Keren-Shaul H Merad M Jung S Amit I Tissue-resident macrophage enhancer landscapes are shaped by the local microenvironment Cell 2014 159 1312 1326 25480296
52 van de Laar L Saelens W De Prijck S Martens L Scott CL Van Isterdael G Hoffmann E Beyaert R Saeys Y Lambrecht BN Guilliams M Yolk Sac Macrophages, Fetal Liver, and Adult Monocytes Can Colonize an Empty Niche and Develop into Functional Tissue-Resident Macrophages Immunity 2016 44 755 768 26992565
53 Loos P Baiwir J Maquet C Javaux J Sandor R Lallemand F Marichal T Machiels B Gillet L Dampening type 2 properties of group 2 innate lymphoid cells by a gammaherpesvirus infection reprograms alveolar macrophages Sci Immunol 2023 8 eabl9041 36827420
54 Schilperoort M Ngai D Sukka SR Avrampou K Shi H Tabas I The role of efferocytosis-fueled macrophage metabolism in the resolution of inflammation Immunol Rev 2023 319 65 80 37158427
55 Zhang S Weinberg S DeBerge M Gainullina A Schipma M Kinchen JM Ben-Sahra I Gius DR Yvan-Charvet L Chandel NS Schumacker PT Efferocytosis Fuels Requirements of Fatty Acid Oxidation and the Electron Transport Chain to Polarize Macrophages for Tissue Repair Cell Metab 2019 29 443 456 e5 30595481
56 Li BWS de Bruijn MJW Lukkes M van Nimwegen M Bergen IM KleinJan A GeurtsvanKessel CH Andeweg A Rimmelzwaan GF Hendriks RW T cells and ILC2s are major effector cells in influenza-induced exacerbation of allergic airway inflammation in mice Eur J Immunol 2019 49 144 156 29762870
57 Monticelli LA Sonnenberg GF Abt MC Alenghat T Ziegler CGK Doering TA Angelosanto JM Laidlaw BJ Yang CY Sathaliyawala T Kubota M Innate lymphoid cells promote lung-tissue homeostasis after infection with influenza virus Nat Immunol 2011 12 1045 1054 21946417
58 Gieseck RL Wilson MS Wynn TA Type 2 immunity in tissue repair and fibrosis Nat Rev Immunol 2018 18 62 76 28853443
59 Lechner AJ Driver IH Lee J Conroy CM Nagle A Locksley RM Rock JR Recruited Monocytes and Type 2 Immunity Promote Lung Regeneration following Pneumonectomy Cell Stem Cell 2017 21 120 134 e7 28506464
60 Weber GF Ashkar S Glimcher MJ Cantor H Receptor-ligand interaction between CD44 and osteopontin (Eta-1) Science 1996 271 509 512 8560266
61 Chen Q Kumar VS Finn J Jiang D Liang J Zhao Y-Y Liu Y CD44high alveolar type II cells show stem cell properties during steady-state alveolar homeostasis Am J Physiol Lung Cell Mol Physiol 2017 313 L41 L51 28473330
62 Burrack KS Morrison TE The role of myeloid cell activation and arginine metabolism in the pathogenesis of virus-induced diseases Front Immunol 2014 5 428 25250029
63 Crompton RA Williams H Campbell L Hui Kheng L Saville C Ansell DM Reid A Wong J Vardy LA Hardman MJ Cruickshank SM An Epidermal-Specific Role for Arginase1 during Cutaneous Wound Repair J Invest Dermatol 2022 142 1206 1216 e8 34710388
64 Katsura H Kobayashi Y Tata PR Hogan BLM IL-1 and TNFα Contribute to the Inflammatory Niche to Enhance Alveolar Regeneration Stem Cell Reports 2019 12 657 666 30930244
65 Bosteels C Van Damme KFA De Leeuw E Declercq J Maes B Bosteels V Hoste L Naesens L Debeuf N Deckers J Cole B Loss of GM-CSF-dependent instruction of alveolar macrophages in COVID-19 provides a rationale for inhaled GM-CSF treatment Cell Rep Med 2022 3 100833 36459994
66 Reyfman PA Walter JM Joshi N Anekalla KR McQuattie-Pimentel AC Chiu S Fernandez R Akbarpour M Chen C-I Ren Z Verma R Single-Cell Transcriptomic Analysis of Human Lung Provides Insights into the Pathobiology of Pulmonary Fibrosis Am J Respir Crit Care Med 2019 199 1517 1536 30554520
67 Morse C Tabib T Sembrat J Buschur KL Bittar HT Valenzi E Jiang Y Kass DJ Gibson K Chen W Mora A Proliferating SPP1/MERTK-expressing macrophages in idiopathic pulmonary fibrosis Eur Respir J 2019 54
68 Adams TS Schupp JC Poli S Ayaub EA Neumark N Ahangari F Chu SG Raby BA DeIuliis G Januszyk M Duan Q Single-cell RNA-seq reveals ectopic and aberrant lung-resident cell populations in idiopathic pulmonary fibrosis Sci Adv 2020 6 eaba1983 32832599
69 Schyns J Bai Q Ruscitti C Radermecker C De Schepper S Chakarov S Farnir F Pirottin D Ginhoux F Boeckxstaens G Bureau F Non-classical tissue monocytes and two functionally distinct populations of interstitial macrophages populate the mouse lung Nat Commun 2019 10 3964 31481690
70 Nokin M-J Durieux F Bellier J Peulen O Uchida K Spiegel DA Cochrane JR Hutton C Castronovo V Bellahcène A Hormetic potential of methylglyoxal, a side-product of glycolysis, in switching tumours from growth to death Sci Rep 2017 7 11722 28916747
71 Rademaker G Hennequière V Brohée L Nokin M-J Lovinfosse P Durieux F Gofflot S Bellier J Costanza B Herfs M Peiffer R Bettendorff L Myoferlin controls mitochondrial structure and activity in pancreatic ductal adenocarcinoma, and affects tumor aggressiveness Oncogene 2018 37 4398 4412 29720728
72 Dinnon KH Leist SR Okuda K Dang H Fritch EJ Gully KL De la Cruz G Evangelista MD Asakura T Gilmore RC Hawkins P SARS-CoV-2 infection produces chronic pulmonary epithelial and immune cell dysfunction with fibrosis in mice Sci Transl Med 2022 14 eabo5070 35857635
73 Huber W Carey VJ Gentleman R Anders S Carlson M Carvalho BS Bravo HC Davis S Gatto L Girke T Gottardo R Orchestrating high-throughput genomic analysis with Bioconductor Nat Methods 2015 12 115 21 25633503
74 Radermecker C Sabatel C Vanwinge C Ruscitti C Maréchal P Perin F Schyns J Rocks N Toussaint M Cataldo D Johnston SL Locally instructed CXCR4hi neutrophils trigger environment-driven allergic asthma through the release of neutrophil extracellular traps Nat Immunol 2019 20 1444 1455 31591573
75 Scott CL Zheng F De Baetselier P Martens L Saeys Y De Prijck S Lippens S Abels C Schoonooghe S Raes G Devoogdt N Bone marrow-derived monocytes give rise to self-renewing and fully differentiated Kupffer cells Nat Commun 2016 7 10321 26813785
76 Stuart T Butler A Hoffman P Hafemeister C Papalexi E Mauck WM Hao Y Stoeckius M Smibert P Satija R Comprehensive Integration of Single-Cell Data Cell 2019 177 1888 1902 e21 31178118
77 Jung S Aliberti J Graemmel P Sunshine MJ Kreutzberg GW Sher A Littman DR Analysis of fractalkine receptor CX(3)CR1 function by targeted deletion and green fluorescent protein reporter gene insertion Mol Cell Biol 2000 20 4106 14 10805752
78 Boring L Gosling J Chensue SW Kunkel SL Farese RV Broxmeyer HE Charo IF Impaired monocyte migration and reduced type 1 (Th1) cytokine responses in C-C chemokine receptor 2 knockout mice J Clin Invest 1997 100 2552 2561 9366570
79 Madisen L Zwingman TA Sunkin SM Oh SW Zariwala HA Gu H Ng LL Palmiter RD Hawrylycz MJ Jones AR Lein ES A robust and high-throughput Cre reporting and characterization system for the whole mouse brain Nat Neurosci 2010 13 133 140 20023653
80 Wende H Lechner SG Cheret C Bourane S Kolanczyk ME Pattyn A Reuter K Munier FL Carroll P Lewin GR Birchmeier C The transcription factor c-Maf controls touch receptor development and function Science 2012 335 1373 1376 22345400
81 Clausen BE Burkhardt C Reith W Renkawitz R Forster I Conditional gene targeting in macrophages and granulocytes using LysMcre mice Transgenic research 1999 8 265 77 10621974
82 Schneider C Nobs SP Heer AK Hirsch E Penninger J Siggs OM Kopf M Frontline Science: Coincidental null mutation of Csf2rα in a colony of PI3Kγ-/- mice causes alveolar macrophage deficiency and fatal respiratory viral infection J Leukoc Biol 2017 101 367 376 27468760
83 Andreata F Blériot C Di Lucia P De Simone G Fumagalli V Ficht X Beccaria CG Kuka M Ginhoux F Iannacone M Isolation of mouse Kupffer cells for phenotypic and functional studies STAR Protoc 2021 2 100831 34585164
84 Rieger AM Nelson KL Konowalchuk JD Barreda DR Modified annexin V/propidium iodide apoptosis assay for accurate assessment of cell death J Vis Exp 2011 2597 21540825
85 Van den Berge K Roux de Bézieux H Roux H Street K Saelens W Cannoodt R Saeys Y Dudoit S Clement L Trajectory-based differential expression analysis for single-cell sequencing data Nat Commun 2020 11 1201 32139671
86 Gu Z Eils R Schlesner M Complex heatmaps reveal patterns and correlations in multidimensional genomic data Bioinformatics 2016 32 2847 2849 27207943
87 Subramanian A Tamayo P Mootha VK Mukherjee S Ebert BL Gillette MA Paulovich A Pomeroy SL Golub TR Lander ES Mesirov JP Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles Proc Natl Acad Sci U S A 2005 102 15545 50 16199517
88 Danaher P Kim Y Nelson B Griswold M Yang Z Piazza E Beechem JM Advances in mixed cell deconvolution enable quantification of cell types in spatial transcriptomic data Nat Commun 2022 13 385 35046414
89 Foroutan M Bhuva DD Lyu R Horan K Cursons J Davis MJ Single sample scoring of molecular phenotypes BMC Bioinformatics 2018 19 404 30400809
