Introduction

The respiratory system comprises the nasal cavity, airways, and lungs, and has the largest surface area exposed to the external environment. It is thus constantly challenged by various stimuli, ranging from odorants, pollutants, toxins, allergens, to viruses, bacteria, and fungi. Prompt detection and appropriate responses to harmful inhalants are critical to maintaining not only the pivotal gas-exchange function of the system but also the overall health or even vitality of an individual. How the host responds to the external challenges and what damages can be inflicted upon the host depend on the types of harmful stimulus, dosage, host genetic background and immune history (Yunis et al. 2023).

Over the past decades, viral infection has been salient risks to human health. The severe acute respiratory syndrome coronavirus (SARS-CoV-1), and H1N1 influenza A virus and SARS-CoV2 have caused pandemics (Peiris 2003; Donaldson et al. 2009; Chen et al. 2020). However, it remains elusive how the host and infectious viruses interact and how the host regulates its responses to the invasion. Both under- and over-reactions by the host appear to engender serious consequences.

Recent studies have shown that several types of epithelial cells expressing the chemosensory receptors play important roles in detecting and responding to the external stimuli in the respiratory system (Carey and Lee 2019;Tizzano and Finger 2013). Among them is a rare type of microvillous cells, tuft cells (Billipp et al. 2021), which have also been found in a number of tissues, including the gastrointestinal tract, pancreas, respiratory airways and lungs, gallbladder and thymus (Bezencon et al. 2008; Miller et al. 2018; Montoro et al. 2018; Nevalainen 1977; Saqui-Salces et al. 2011; DelGiorno et al. 2020). These cells can sense noxious stimuli and invading pathogens, such as allergens, bacteria, protists and helminths (Gerbe et al. 2016; Howitt et al. 2016; Lei et al. 2018; Luo et al. 2019; von Moltke et al. 2016). Tuft cells are known to express canonical chemosensory signaling pathways, including the sweet and umami taste receptor subunit Tas1r3 (Howitt et al. 2020), bitter taste receptors Tas2rs (Luo et al. 2019; Imai et al. 2020), and olfactory receptors Vmn2r26 (Xiong et al. 2022), as well as their downstream signaling proteins including the heterotrimeric G protein subunits α-gustducin, Gβ1 and Gγ13, phospholipases, and the transient receptor potentialion channel Trpm5 (Howitt et al. 2016; Doyle et al. 2023). Even the signal transduction of the Sucnr1 receptor for succinic acid and free fatty acid receptor 2 Ffar2 for propionate, characteristic metabolites of certain microbes, is mediated by α-gustducin and / or Trpm5 signaling proteins to monitor changes in gut microbiota and airway infection. And ablation of either α-gustducin or Trpm5 diminishes the intestinal and tracheal epithelia’s responses to the infection of the parasitic helminths and microbial colonization (Lei et al. 2018; Nadjsombati et al. 2018; Schneider et al. 2018; Keshavarz et al. 2022; Perniss et al. 2023).

Once stimulated, tuft cells can release a number of output signaling molecules, including IL-25, prostaglandins (Kotas et al. 2022), cysteinyl leukotrienes (Bankova et al. 2018; Bankova and Boyce 2018; McGinty et al. 2020), acetylcholine (Hollenhorst et al. 2020;Krasteva et al. 2011), and other substances, to stimulate type 2 innate lymphoid cells, nerve terminals or adjacent cells, initiating innate and adaptive immune responses, leading to the elimination of irritants and pathogens, and eventually to tissue repair and remodeling (Finger et al. 2003; Tizzano et al. 2010; Saunders et al. 2014). However, tuft cells are heterogeneous in terms of gene expression patterns, signal transduction pathways and output effectors (Banerjee et al. 2018). For example, some tuft cells express subsets of Tas2r receptors, or cysteinyl leukotriene receptor 3 (Cysltr3) while others do not (Luo et al. 2019; Imai et al. 2020; Bankova et al. 2018). It seems important to determine ligand profiles, output signals and physiological roles for tuft cells residing in different tissues to fully understand their protective functions against various pathogens and at different stages of disease progression.

Although they are normally only found in the nasal cavity and trachea, rarely found distal to the lung hilus (Krasteva et al. 2011; Tizzano et al. 2010; Saunders et al. 2014), tuft cells can be ectopically formed in the distal lung following severe viral infection or chemical damage (Rane et al. 2019; Barr et al. 2022; Roach et al. 2022; Melms et al. 2021). These ectopic tuft cells of molecular diversities emerge around 12 days post viral infection, suggesting that they may not play any role in the initial response to the infection or during viral clearance (Rane et al. 2019; Barr et al. 2022). However, it is unclear whether each subtype of these diverse tuft cells makes unique contributions to the subsequent inflammation resolution, tissue repair and remodeling. In this study, we conditionally nullified the expression of a heterotrimeric G protein subunit, Gγ13, in the choline acetyltransferase (ChAT)-expressing tuft cells, which resulted in the generation of fewer tuft cells, but conferred even severer injury, slower recovery and higher fatality following the H1N1 influenza infection. Our data indicate that subsets of ectopic tuft cells with or without Gγ13 expression play different important roles in the inflammation resolution and tissue repair.

Results

Severe injury induces the generation of taste signaling proteins-expressing ectopic tuft cells

To verify the previously reported ectopic generation of lung tuft cells upon severe injury, we intranasally inoculated adult mice with a sublethal dose of H1N1 viruses, and observed the body-weight changes over a course of 25 days post infection (dpi) (Figure 1A). The mice were then killed and their lungs were dissected out for analyses (Figure 1A). The mice showed bodyweight loss over the days post infection with a peak loss of up to 25% of their original bodyweight around 8 to 9 dpi, and their lungs displayed lesion areas, of which the histological analysis showed massive alveolar damage and severe epithelial dysplasia with altered tissue structures as well as an increased number of CD45+ immune cells in comparison with the uninfected control (Figure 1, B, C and D). Analysis of the H1N1 injured-lungs of the ChAT-Cre: Ai9 mice, which express tdTomato proteins in the ChAT-expressing cells, showed that about 78% ChAT-expressing cells also expressed the tuft cell marker Dclk1 whereas about 60% and 12% of ChAT-expressing cells expressed the other two taste signaling proteins Gα-gustducin and Plcβ2, respectively (Figure 1, E and F). In addition, chemically induced severe lung injury with bleomycin, lipopolysaccharide and powder of house dust mite extracts also led to the generation of dysplastic tuft cells, but to a lesser extent (Figure 1— figure supplement 1), confirming the severe injury-induced formation of ectopic tuft cells in the lung parenchyma (Rane et al. 2019; Barr et al. 2022; Huang et al. 2022).

H1N1 infection causes severe damage to mouse lung.

(A) Schematic of H1N1 intranasal inoculation procedure. (B) Bodyweight changes of uninfected control and H1N1-infected mice over 25 days post infection. Data are presented as means ± SD (n=6). Unpaired two-tailed student t tests were performed. (C) Representative images of lungs at 0 and 25 dpi, and their H&E tissue sections (Bars: 5 mm in the whole lung images and 100 µm in the lung section images). (D) Immunostaining images of lung sections of 0 and 25 dpi with an antibody to CD45. (E) Immunostaining of the ChAT-Cre: Ai9 lung sections of 0 and 25 dpi with the antibody to the tuft cell marker Dclk1 (green), and ChAT (Ai9, red tdTomato). (F) Immunostaining of the ChAT-Cre: Ai9 lung sections of 25 dpi with antibodies to Dclk1 (green), Gα-gustducin (green) or Plcβ2 (green) and ChAT (Ai9, red tdTomato). Scale bars in (D, E, F): 50 µm. ***P<0.001.

Figure 1—figure supplement 1. Chemically induced severe injury also engenders ectopic formation of lung tuft cells.

Ectopic lung tuft cells respond to bitter tasting compounds

To assess how the severe lung injury may affect the bitter taste receptor expression, we quantified and comparatively analyzed Tas2r expression levels by performing quantitative reverse transcription-PCR on the WT lung tissues before and after H1N1 infection. The results showed that the expression levels of most of murine 35 Tas2rs remained unchanged while Tas2r105, Tas2r108, Tas2r118, Tas2r137 and Tas2r138 were upregulated, and Tas2r135 and Tas2r143 were downregulated (Figure 2A), which is largely consistent with the single cell RNAseq data (Figure 2—figure supplement 1) (Barr et al. 2022), suggesting that the upregulated Tas2rs are indeed expressed in these ectopic tuft cells.

Expression and functional responses of bitter taste receptors.

(A) Quantitative RT-PCR analysis of bitter taste receptor expression in the lung tissues at 0 dpi versus 25 dpi. Expression of Tas2r105, Tas2r108, Tas2r118, Tas2r137 and Tas2r138 was upregulated whereas Tas2r135 and Tas2r143 were downregulated at 25 dpi. Data are presented as means ± SD (n=3). Unpaired two-tailed student t tests were performed. (B-K) Calcium imaging and response traces of FACs-sorted tuft cells isolated from H1N1-infected lungs at 25 dpi. (B, D, F, H, and J) Left panels are grey images of tuft cells at rest in the HBSS buffer whereas right panels are grey images of the same cells after being stimulated by denatonium benzoate (D.B.) or quinine. Activated cells appeared brighter. (C, E, G, I, K) traces of calcium responses of tuft cells on the left. Some tuft cells responded to D.B, which was inhibited by AITC (B, C), gallein (F, G), or U73122 (J, K) whereas others responded to quinine, which was also inhibited by AITC (D, E) or gallein (H, I). Scale bar: 50 µm. **P<0.01, ***P<0.001.

Figure 2—figure supplement 1. Expression of taste2 receptors in tuft cell clusters.

Figure 2—figure supplement 2. FACs sorting of ectopic tuft cells.

To determine whether the taste signal transduction pathways present in these ectopic lung tuft cells are functional, we isolated these tuft cells from the ChAT-Cre: Ai9 mice at 25 dpi by fluorescence-activated cell sorting (FACs) (Figure 2—figure supplement 2) and assessed their intracellular calcium responses to two common bitter taste substances denatonium benzoate (D.B.) and quinine. It is known that D.B. activates mouse bitter taste receptor Tas2r105 whereas quinine can stimulate multiple Tas2rs, including 3 of the 5 upregulated bitter taste receptors: Tas2r105, Tas2r108 and Tas2r137 (Lossow et al. 2016). The results showed that subsets of tuft cells responded to these two compounds by increasing the intracellular calcium contractions. Application of the bitter taste inhibitor allyl isothiocyanate (AITC), G protein βγ moiety inhibitor gallein or Plcβ2 inhibitor U73122 was able to completely inhibit these cells’ calcium responses to the bitter compounds (Figure 2, B-K), indicating that the Tas2rs and their downstream signaling proteins are functional in subsets of tuft cells, and may play some physiological roles during the transition from inflammatory response to resolution.

Conditional nullification of the taste signaling protein Gγ13 exacerbates the H1N1-inflicted disease severity and suppresses ectopic tuft cell expansion

To determine what role the taste signaling proteins may play in the H1N1 infection-triggered lung injury, we inoculated H1N1 to the whole-body Trpm5 knockout mice (Trpm5−∕−) and the conditional Gγ13 knockout mice (ChAT-Cre: Gng13flox∕flox, i.e., Gng13-cKO) as well as WT control mice. Results indicated that while WT and Trpm5−∕− mice showed similar bodyweight loss/recovery curves with nearly no mortality, Gng13-cKO mice exhibited significantly more bodyweight loss with a maximum of 30% instead of 20% of the original bodyweight, and a slower recovery process with a shifted maximum bodyweight loss from 9 dpi to 10 dpi, and a significant increase in the mortality rate to 37.5%, the majority of which occurred between 8 and 12 dpi (Figure 3, A and B). Furthermore, the H1N1-infected lungs of WT and Trpm5−∕− mice at 25 dpi displayed similar injured surfaces areas and damaged tissue volumes whereas those of Gng13-cKO showed significantly greater injured areas and volumes (Figure 3, C and D, Figure 3—figure supplement 1). And the ratios of wet-to-dry lung masses of both WT and Gng13-cKO mice at 25 dpi were significantly greater than their respective ones at 0 dpi whereas the ratio of the Gng13-cKO mice at 25 dpi was even greater than that of WT mice at 25 dpi (Figure 3E), supporting the notion that Gng13-cKO mice experienced severer damage following H1N1 infection comparing with WT control.

H1N1 infection inflicts differential lung injury and tuft cell dysplasia on WT versus mutant mice.

(A) Gng13-cKO mice showed Figure 3 continued on next page greater bodyweight loss than WT or Trpm5−∕− mice over days post infection. Data are presented as means ± SD (n=6). Unpaired two-tailed ttests were performed. (B) Kaplan-Meier survival curves of WT, Trpm5−∕− and Gng13-cKO mice following H1N1 inoculation. A significant number of Gng13-cKO mice died between 9 and 12 dpi (i.e., 2, 4, 2 mice and 1 mouse died at 9, 10, 11, 12 dpi, respectively), reducing the overall survival rate to 62.5%, significantly lower than those of Trpm5−∕− or WT mice, of which nearly all survived. Data are presented as means ± SEM (n= 44, 17 and 14 for WT, Trpm5−∕− and Gng13-cKO mice, respectively). The curves were determined by a log-rank test. (C, D) Images and statistical analysis of the injured areas on the lungs of WT, Trpm5−∕− and Gng13-cKO mice at 25 dpi. The injured areas that are marked by dashed lines were significantly greater in the Gng13-cKO mice than in Trpm5−∕− or WT mice while no significant difference was found between the latter two. Data are presented as means ± SD (n=5), and unpaired two-tailed student t tests were performed. (E) Comparative analysis of wet to dry weight ratios of WT and Gng13-cKO lungs at 0 and 25 dpi. While no difference in the ratio was found between WT and Gng13-cKO mice at 0 dpi, the ratios of these mice at 25 dpi were significantly greater than their corresponding ones at 0 dpi. However, the Gng13-cKO ratio was even greater than WT ratio at 25 dpi. Data are presented as means ± SD (n=6), and unpaired two-tailed student t tests were performed. (F-M) Identification of tuft cells in the Krt5-expressing tissue areas of H1N1-injured lungs of WT and Gng13-cKO mice at 14, 20, 25 and 50 dpi using antibodies to Krt5 and to the tuft cell marker Dclk1. The densities of tuft cells were significantly higher in WT than in Gng13-cKO lungs at 20, 25 and 50 dpi, but not at an earlier time point of 14 dpi. Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001.

Figure 3—figure supplement 1. Statistical analysis of the injured tissue volumes on the lungs of WT, Trpm5−∕− and Gng13-cKO mice at 25 dpi.

Figure 3—figure supplement 2. Ectopic generation of tuft cells in WT, Trpm5−∕− and Gng13-cKO mice over days post H1N1 infection.

To investigate how the abolishment of taste signaling proteins affects the formation of ectopic lung tuft cells in response to H1N1 infection, we performed immunohistochemistry on lung tissue sections, using an anti-Krt5 antibody to demarcate H1N1-damaged areas and an anti-Dclk1 anti-body to identify tuft cells. We found as early as 14 dpi a small number, i.e., ∼6 and ∼4 tuft cells per mm2 of the Krt5-expressing tissue area in WT and Gng13-cKO lungs, respectively, and no significant differences in the tuft cell density between these two genotypes were found at this time point. However, at 20 and 25 dpi, the densities of tuft cells in WT were increased to 120 and 175 cells per mm2, respectively, which are significantly more than the corresponding densities of Gng13-cKO, i.e., 17 and 57 cells / mm2, respectively; and at 50 dpi, the tuft cell density in WT was decreased to 102 cells / mm2, which, however, was still significantly higher than 56 cells / mm2 of Gng13-cKO (Figure 3, F-M). In contrast, the densities of tuft cells in Trpm5−∕− mice were nearly identical to those of WT at both 14 and 25 dpi, and no significant difference was found between WT and Trpm5−∕− mice (Figure 3—figure supplement 2, A-E). Thus, although tuft cells first appeared at 14 dpi in all WT, Trpm5−∕− and Gng13-cKO mice, and no statistic difference in their densities was found among these mice, the subsequent densities were increased in all WT, Trpm5−∕− and Gng13-cKO mice; and the increase was much less in the Gng13-cKO mice than in WT or Trpm5−∕−, and consequently, the densities in the Gng13-cKO mice were significantly less than in WT or Trpm5−∕− mice at 20, 25 and 50 dpi (Figure 3—figure supplement 2F).

To determine how many ectopic tuft cells in the infected WT lung express Gγ13 and whether any remaining ectopic tuft cells in the infected Gng13-cKO lung express Gγ13, double immunostaining was carried out with antibodies to Dclk1 and to Gγ13 on the lung tissue sections of 25 dpi of both WT and Gng13-cKO mice. The results showed that about 28.6% of Dclk1+ ectopic tuft cells express Gγ13 in WT sections while none of the remaining Dclk1+ tuft cells in the Gng13-cKO sections was Gγ13+ (Figure 3—figure supplement 2, G and H), indicating an effective conditional ablation of Gγ13 expression. Furthermore, reanalysis of the previously published single tuft cell RNAseq dataset GSE197163 showed that about 57% of the Trpm5-GFP+ tuft cells were Gγ13+ in WT mice, some of which also expressed Alox5, a key enzyme to the biosynthesis of pro-resolving mediators (Figure 3—figure supplement 2, I and J)(Barr et al. 2022, Halade et al. 2022). This result is consistent with the previous report indicating the heterogeneity of dysplastic tuft cells in the severely injured lungs, and this study ablated the Gγ13+ tuft cells including some expressing the pro-resolving enzyme Alox5.

Gγ13 mutant mice display stronger inflammation and augmented pyroptosis following H1N1 infection

To determine whether Gγ13 conditional knockout affects the inflammation resolution following H1N1 infection, we performed immunohistochemistry on the lung sections with an antibody to the immune cell marker CD64. The results showed that at both 20 and 25 dpi, the Gng13-cKO mice had significantly more CD64+ immune cells in the injured areas than WT, although both types of mice indicated a trend of inflammation resolution with a reduction in the number of immune cells from 20 to 25 dpi (Figure 4, A and B).

Stronger inflammatory response and upregulated expression of pyroptotic genes in the H1N1-infected Gng13-cKO lungs.

(A, B) Immunostaining of WT and Gng13-cKO lung tissue sections with an antibody to the immune cell marker CD64 indicates significantly more immune cells in the mutant lung than in WT at both 20 and 25 dpi. (C) qRT-PCR indicates that expression levels of gasdermin D and E (Gsdmd, Gsdme) and other proteins of the pyroptosis pathway: Tlr4, Nlpr3, Asc, Casp1 and Il-1b were significantly upregulated in both WT and Gng13-cKO lungs at 25 dpi comparing with those at 0 dpi. Further, the expression levels of these genes in Gng13-cKO were significantly higher than in WT at 25 dpi. Data are presented as means ± SD (n=3). Unpaired two-tailed student t tests were performed. (D, E) Western blot analysis shows that the protein levels of NLRP3, full-length Gasdermin D (Gsdmd-FL), N-terminal fragment of Gasdermin D (Gsdmd-N), pro-caspase 1 (pro-Casp1), cleaved caspase 1(Cleaved Casp-1 (p20)), and IL-1β in the Gng13-cKO lungs were significantly higher than the corresponding levels of WT at both 20 and 25 dpi (β-actin as an internal control). And over this phase, the expression levels of these proteins were reduced at 25 dpi in comparison with those at 20 dpi in both genotypes. Data are presented as means ± SD (n=4), and unpaired two-tailed t tests were performed. *P<0.05, *P<0.01, ***P<0.001, ****P<0.0001.

Figure 4—figure supplement 1. qRT-PCR analysis of gasdermin gene expression in WT and Gng13-cKO lungs.

We also carried out quantitative reverse transcription-PCR (qRT-PCR) to assess the expression of pyroptosis-related gasdermin genes, and found that Gsdmd and Gsdme were the two most abundantly expressed ones (Figure 4—figure supplement 1). Further analysis showed that before H1N1 infection, the expression levels of Gsdmd and Gsdme were similar between WT and Gng13-cKO mice at 0dpi; and at 25 dpi, their expression levels were significantly higher than those at 0 dpi in both WT and Gng13-cKO mice; furthermore, Gsdmd and Gsdme expression levels were significantly higher in the Gng13-cKO mice than in WT at 25 dpi (Figure 4C). We then performed additional qRT-PCR to examine the expression patterns of other genes involved in the pyroptosis, Tlr4, Nlpr3, Asc, Casp1 and Il-1b, and found a similar pattern for all these genes, i.e., both WT and Gng13-cKO mice showed more expression at 25 dpi than at 0 dpi, whereas Gng13-cKO exhibited even higher expression levels than WT at 25 dpi (Figure 4C).

To validate their expression at the protein level, we isolated proteins from the injured lung areas and performed Western blot analysis with antibodies to NLRP3, full-length Gsdmd, N-terminal fragment of Gsdmd, pro-Casp-1, cleaved Casp-1, IL-1β and β-actin as an internal reference (Figure 4D). Quantitative analysis showed that at both 20 and 25 dpi, the protein levels of all these pyroptosis pathway proteins were greater in the Gng13-cKO mice than in WT mice, and comparison within the same genotype showed that the levels of these proteins at 20 dpi were higher than at 25 dpi in both WT and Gng13-cKO mice (Figure 4E).

To determine which cells underwent pyroptosis, we performed immunohistochemistry on the injured lung sections with antibodies to the two most highly expressed gasdermins, Gsdmd and Gsdme. Results showed that at 25 dpi, there were about 180 and 77 cells/mm2 in Gng13-cKO mice expressing Gsdmd and Gsdme, respectively, which are significantly higher than the corresponding densities of 88 and 50 cells/mm2 in WT (Figure 5, A and B; Figure 5—figure supplement 1, A and B). Colocalization of immunostaining with antibodies to Gsdmd, the macrophage marker F4/80 and the epithelial cell marker EpCAM, indicated that the majority of Gsdmd+ cells, i.e., about 76.7% and 80.6% in WT and Gng13-cKO mice, respectively, at 25 dpi were also F4/80+, and the remaining Gsdmd+ cells, i.e. 21.9% and 17.8% of Gsdmd+ cells in WT and Gng13-cKO mice, respectively, were EpCAM+, suggesting that most of Gsdmd-expressing cells were macrophages whereas the remaining were epithelial cells (Figure 5C, Figure 5—figure supplement 2A). Interestingly, nearly all (i.e., 98.5% and 95.5% in WT and Gng13-cKO, respectively) of Gsdme+ cells were F4/80+, indicating that almost all Gsdme-expressing cells were macrophages in the two genotypes at 25 dpi (Figure 5—figure supplement 2B). Finally, the density of cells expressing IL-1β and caspase 3 were 153 cells/mm2 and 267 cells/mm2, respectively, in Gng13-cKO mice at 25 dpi, and significantly greater than the corresponding density of 63 and 66 cells/ mm2 in WT mice (Figure 5, D and E; Figure 5— figure supplement 1, C and D), indicating more programmed cell death in the H1N1-infected Gng13-cKO lung. Similar to that of Gsdme, the majority, i.e., 93% and 91%, of IL-1β-expressing cells in WT and Gng13-cKO mice, respectively, were F4/80+ macrophages (Figure 5—figure supplement 2C).

Immunohistochemical analyses of gasdermin D and IL-1β expressing cells.

(A, B) Immunostaining of H1N1-infected lung sections Figure 5 continued on next page at 25 dpi with an antibody to gasdermin D (Gsdmd) indicates significantly more gasdermin D-expressing cells in the Gng13-cKO sections than in WT. Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. (C) Co-immunostaining shows that nearly the same percentage, i.e., 98.5% and 95.5% of Gsdmd+ cells in WT and Gng13-cKO, respectively, were F4/80+. (D) Immunostaining of H1N1-infected lung sections at 25 dpi with an antibody to IL-1β indicates significantly more IL-1β-expressing cells in the Gng13-cKO sections than in WT. Data are presented as means ± SD (n=4), and unpaired two-tailed t tests were performed. ***P<0.001, ****P<0.0001.

Figure 5—figure supplement 1. Expression of gasdermin E and caspase 3 in the H1N1-injured lungs.

Figure 5—figure supplement 2. Partial colocalization of gasdermin D, E and IL-1β to epithelial cells and macrophages.

To obtain a snapshot of cells indeed undergoing pyroptosis, we intranasally administrated Sytox to label dead cells (Xi et al. 2021). Histochemical analysis showed that Gng13-cKO mice at 20 and 25 dpi had 211 and 96 Sytox+ cells/mm2, respectively, significantly more than the corresponding 96 and 62 cells/mm2 in WT (Figure 6, A and B). Double staining of Sytox with an anti-CD64 antibody indicated that more than half of CD64+ macrophages, i.e., 61.5% and 57% in WT and Gng13-cKO mice at 25 dpi, respectively, were Sytox+, undergoing pyroptosis (Figure 6C). Together, these data indicate that the injured lungs exhibited immune cell infiltration, expressed pyroptotic genes, and cell death of macrophages and some epithelial cells; and from 20 to 25 dpi, the injured animals showed some inflammation resolution; but at both time points, the Gng13-cKO mice showed stronger inflammatory responses and slower inflammation resolution than WT.

Gng13 conditional nullification leads to more pyroptotic cells in the injured lungs.

(A, B) Histological analysis of lung tissue sections following Sytox administration reveals that significantly more Sytox+ cells were found in the mutant lung sections than the corresponding WT tissues at both 20 and 25 dpi while in both types of mice the numbers of the Sytox+ cells at 25 dpi seemed to be fewer than that at 20 dpi. (C) Immunostaining with an antibody to the macrophage marker CD64 indicates that most of Sytox+ cells were macrophages in both WT and Gng13-cKO mice (circled by dashed lines in the bright field images on the far right). Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. *P<0.05, **P<0.01.

Gγ13 disruption leads to severer leakage of the lung epithelia

Increased pyroptosis in the mutant lungs prompted us to assess the integrity of lung epithelia. We collected total bronchoalveolar lavage fluid (BALF) from both WT and Gng13-cKO mice at 0, 20 and 25 dpi, and the average numbers of cells per mouse present in the BALF from WT and Gng13-cKO mice at 0 dpi were determined to be 0.23×106 and 0.37×106 cells, respectively, and no significant difference was found between these BALF samples. However, the cell numbers were significantly increased to 1.22 x106 and 5.19 x106 in WT and Gng13-cKO BALF at 20 dpi, and 0.93 x106 and 2.63 x106 in WT and Gng13-cKO BALF at 25 dpi, respectively, and the latter was still significantly more than the former (Figure 7A). The BALF protein content also showed a similar pattern, i.e., a basal amount of protein was found in both WT and Gng13-cKO BALF at 0 dpi, that is, 0.36 and 0.33 mg/mL, respectively, and no significant difference was found between the two groups. At 20 dpi, however, both WT and Gng13-cKO BALF contained much more proteins than their corresponding samples at 0 dpi, but at 20 dpi, the Gng13-cKO BALF had more proteins than that of WT (2.01 versus 0.88 mg/mL), and the same is true at 25 dpi (1.65 versus 0.52 mg/mL) (Figure 7B). In WT, the BALF protein concentration was significantly decreased from 0.88 mg/mL at 20 dpi to 0.52 mg/mL at 25 dpi whereas the decrease of the BALF protein concentration in the Gng13-cKO was insignificant over the same period. The lactate dehydrogenase (LDH) activity and IL-1β content in the BALF also showed similar patterns. No significant difference in the basal LDH activity or IL-1β content was found between WT BALF (LDH, 33.81 U/L; IL-1β, 3.60 pg/ml) and Gng13-cKO BALF (LDH, 32.17 U/L; IL-1β, 3.35 pg/ml) at 0 dpi. The LDH activity then rose significantly to 73.86 U/L and 91.67 U/L at 20 dpi in WT and Gng13-cKO BALF, respectively, followed by a significant decrease to 53.12 U/L and 76.40 U/L, respectively, at 25 dpi. However, at both time points, the Gng13-cKO LDH activities were always greater than the corresponding WT activities (Figure 7C). WT BALF IL-1β content also significantly increased to 38.82 pg/ml at 20 dpi, then significantly decreased to 22.03 pg/ml at 25 dpi, while the Gng13-cKO BALF IL-1β levels at the two time points were 75.57 and 82.12 pg/ml, respectively, significantly higher than the corresponding WT levels. Interestingly, different from the LDH activity dynamics, the Gng13-cKO BALF IL-1β level at 25 dpi did not appear to be reduced comparing with that at 20 dpi (Figure 7D). These results indicate that while the WT injured lungs showed significant improvements in the epithelial leakage from 20 to 25 dpi, the Gng13-cKO injured lungs displayed higher levels of leakage and a slower repair with a significant reduction only in the BALF LDH level, but not in the total BALF protein or IL-1β level, from 20 to 25 dpi.

Bronchoalveolar lavage fluid (BALF) assays.

(A) Average numbers of cells per mouse found in both WT and Gng13-cKO BALFs at both 20 and 25 dpi were significantly more than those at 0 dpi while those of the Gng13-cKO BALF were even more than those of WT at both 20 and 25 dpi. Data are presented as means ± SD (n=6), and unpaired two-tailed student t tests were performed. (B) The protein contents of both WT BALF and Gng13-cKO BALF at 20 dpi and 25 dpi were significantly more than their corresponding ones at 0 dpi, but the Gng13-cKO BALF had significantly more proteins than WT at both 20 and 25 dpi. And in WT, the protein content was significantly reduced from 20 dpi to 25 dpi, but the reduction in the Gng13-cKO was insignificant. Data are presented as means ± SD (n=5), and unpaired two-tailed student t tests were performed. (C) Lactate dehydrogenase (LDH) activity assays showed that at both 20 and 25 dpi, LDH activities in both WT and Gng13-cKO BALF were significantly higher than their corresponding ones at 0 dpi, but the Gng13-cKO BALF at both 20 and 25 dpi showed higher activities than the corresponding WT BALF. And both WT and Gng13-cKO BALF displayed significant LDH activity reduction from 20 to 25 dpi. Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. (D) IL-1β ELISA assays showed a similar patter to that of LDH, except that the reduction in the IL-1β concentration from 20 to 25 dpi was insignificant in the Gng13-cKO BALF. Data are presented as means ± SD (n=4), and unpaired two-tailed student t tests were performed. *P<0.05, **P<0.01, ***P<0.001 ****P<0.0001.

Gng13 conditional nullification exacerbates H1N1-induced fibrosis

To assay how the Gng13 mutation affects the recovery of the severely injured lungs, we first conducted qRT-PCR to determine the expression levels of fibrotic genes: Col1a1, Fn1 and Timp1. The results showed that at 0 dpi, there was no significant difference in the expression levels of these three genes between WT and Gng13-cKO lungs; and at 25 dpi, however, the expression levels of these genes were significantly increased in both genotypes comparing with their corresponding levels at 0dpi. Comparison between the two genotypes revealed that at 25 dpi, the expression levels of all these three genes were significantly higher in the Gng13-cKO lungs than in WT. At 50 dpi, the expression levels of these genes showed a downward trend in both WT and Gng13-cKO lungs, consistent with the phase of return to homeostasis of the remodeled lungs. More specifically, the Col1a1 expression was reduced to a similar level in WT and Gng13-cKO, but still significantly higher than those at 0 dpi. For Fn1 and Timp1, their expression in WT at 50 dpi was much reduced, to a level similar to those at 0 dpi; and their expression in Gng13-cKO at 50 dpi was reduced comparing with that at 25 dpi, but still higher than that at 0 dpi or in WT at 50 dpi (Figure 8A).

More severe fibrosis in the H1N1-infected Gng13-cKO lungs.

(A) qRT-PCR analysis of Col1a1, Fn1 and Timp1 expression. While the expression levels of these three genes at 20 dpi in both WT and Gng13-cKO injured lungs were significantly greater than their corresponding ones at 0 dpi, the expression in Gng13-cKO at 25 dpi was much stronger than WT at 25 dpi. At 50 dpi, the expression levels of all these three genes were reduced comparing with those at 25 dpi in both WT and Gng13-cKO mice, but those of Fn1 and Timp1 in Gng13-cKO mice at 50 dpi were still significantly higher than WT at 50 dpi. Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. (B, C) Masson’s trichrome staining showed that at 25 dpi, the fibrosis percentages in WT and Gng13-cKO were significantly more than their corresponding ones at 0 dpi, which were reduced at 50 dpi. Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001.

Masson’s trichrome staining was used to determine the extent of fibrosis. The results showed that 3.2% and 3.1% of fibrosis areas were found in WT and Gng13-cKO lungs at 0 dpi, respectively, and no significant difference was found between these two genotypes. At 25 dpi, the percentages of fibrosis areas in WT (12.3%) and Gng13-cKO (26.3%) were significantly more than their corresponding ones at 0 dpi. And between the two genotypes, the fibrosis in Gng13-cKO was even more than that in WT at 25 dpi. At 50 dpi, the fibrosis in WT (7.9%) and Gng13-cKO (15%) appeared to be reduced comparing with those at 25 dpi, but they were still more than those at 0 dpi; and furthermore, fibrosis in Gng13-cKO at 50 dpi seemed to be more than that of WT at 50 dpi (Figure 8 B and C). These results indicate that Gng13-cKO mice recovered not as fast as WT mice and displayed a higher level of fibrosis.

Discussion

Multiple subtypes of tuft cells have been found in various organs, but their distinct functions are not well understood. In this study, we used a conditional gene knockout mouse strain to nullify Gγ13 expression in the ChAT-expressing tuft cells, and found that the number of ectopic tuft cells was reduced in the H1N1 virus-injured lung areas, and all of these remaining tuft cells were Gγ13-negative. Furthermore, the Gng13-cKO mice displayed larger injured areas in the infected lungs, along with heavier immune cell infiltration, severer and prolonged inflammation, increased pyroptosis and cell death, exacerbated lung epithelial leakage and aggravated fibrosis, increased body-weight loss and heightened fatality. These results imply that the subsets of Gγ13-positive versus negative ectopic tuft cells may play opposite but perhaps balanced roles in the inflammation resolution, tissue repair and functional recovery, and the absence of Gγ13-positive tuft cells engenders hyper inflammation, consequently leading to aggravated outcomes.

Many signaling proteins are shared by tuft cells and taste receptor cells, including the G protein-coupled Tas1r and Tas2r receptors (Zhao et al. 2003; Max et al. 2001; Matsunami et al. 2000; Chandrashekar et al. 2000; Adler et al. 2000), the heterotrimeric G proteins subunits Gα-gustducin, Gα14, Gβ1, Gβ3 and Gγ13 (Huang et al. 1999; Wong et al. 1996; Tizzano et al. 2008), phospholipase Cβ2 (Plcβ2) and Trpm5. Tuft cells can be deemed as solitary chemosensory cells that are dispersed over many organs of the body, monitoring internal environments, whereas taste receptor cells are aggregated into specialized end organs of taste, taste buds, located in the oral cavity of mammals, evaluating food’s nutritional value or potential toxicity before ingestion (Lindemann 1996;O’Leary et al. 2019). However, there are some important differences between taste buds and tuft cells. For example, distinct subsets of taste bud cells express Tas1rs and Tas2rs, which are employed to detect high-energy carbohydrates, nutritious amino acids, or potentially harmful and toxic substances, respectively, conveying the chemical information via different gustatory nerve fibers to the central nervous system, generating hedonic sweet and umami taste perception or unpleasant bitter taste (Doyle et al. 2023). In contrast, individual tuft cells can express both Tas1rs and Tas2rs, and the latter are indeed also involved in the detection of infectious microbes and invading parasites that pose risks to the host’s health whereas the function of the former is less clear, and seems to regulate Tas2rs’ signaling or tuft cell homeostasis (Howitt et al. 2020). Results of our present study indicate that Tas2rs expressed by the ectopic tuft cells are also functional, and could be involved in monitoring changes in the extracellular milieu of the injured areas during the switch from the inflammatory response to inflammation resolution, regulating disease progression.

Our results show that ectopic tuft cells first appear around day 12 to 14 post influenza infection and then the tuft cell number increases over the next two months, which is consistent with a previous report (Rane et al. 2019). However, our study further indicates that the conditional knockout of Gng13 in the ChAT-expressing tuft cells did not affect the generation of dysplastic tuft cells up to 14 days post H1N1 infection, but did decrease ectopic tuft cells’ capability to expand at 20 dpi or later, suggesting that different molecular mechanisms are employed for the initial generation versus subsequent expansion of these tuft cells: the former does not require the contribution from Gγ13-mediated signal transduction whereas the latter does require positive feedback signals released by Gγ13-mediated signaling pathways. It is still unclear how chemically or virally induced major injuries generate these dysplastic tuft cells in the distal lung tissues. Severe injuries can activate quiescent lineage-negative stem/progenitor cells in the pulmonary epithelia via a Notch signaling pathway to proliferate and differentiate into different types of pulmonary epithelial cells to repair and remodel the wounded tissues (Rane et al. 2019; Vaughan et al. 2015). But the exact signals that activate the quiescent stem cells are yet to be identified, although these signals can presumably originate from the apoptotic cells or activated immune cells following the chemical injuries or viral infection. Once activated, the stem/progenitor cells express the transcription factor Pou2f3 at a cell-type specification stage, which has been shown to be critical to the formation of both tuft cells throughout the body and Tas1r- and Tas2r-expressing receptor cells in taste buds, as nullification of the Pou2f3 gene ablates all these taste receptor cells and tuft cells (Gerbe et al. 2016; Matsumoto et al. 2011; Ualiyeva et al. 2020; Yamashita et al. 2017). Pou2f3 may interact with other accessory proteins such as POU domain class 2 associating factors 2 and 3, i.e., Pou2af2 and Pou2af3, respectively, to fine-tune gene expression patterns and direct these cells to differentiate into subtypes of tuft cells or taste receptor cells (Wu et al. 2022). On the other hand, a negative regulatory mechanism has been reported, suggesting that the activating transcription factor 5 (ATF5) appears to inhibit the generation of intestinal tuft cells since more tuft cells were found in the ATF5-deficient ileum (Nakano et al. 2023). Presumably, some intrinsic and extrinsic factors from the severely damaged lung tissues may regulate the activation of Notch, Pou2f3 and ATF5, and together with additional transcription factors such as Gfi1b and Spdef, control the generation of subtypes of ectopic tuft cells in the injured lung.

The proliferation of intestinal tuft cells is dependent on Trpm5 (Howitt et al. 2016; Lei et al. 2018; Luo et al. 2019), but the tracheal tuft cells or sweet, umami and bitter taste receptor cells in taste buds do not require Trpm5 for their expansion (Damak et al. 2006; Zhang et al. 2003; Finger et al. 2005; Hollenhorst et al. 2022). We and others have found that Trpm5 is not required for the proliferation of the dysplastic pulmonary tuft cells (Barr et al. 2022; Huang et al. 2022), but conditional nullification of Gng13 ablated part of these ectopic tuft cells and reduced the total number of these cells, indicating that different tuft cells utilize different mechanisms to control their proliferation. In the small intestine, the hyperplasia of tuft cells as well as goblet cells is regulated by the Trpm5-IL25-ILC2-IL13 molecular circuit. Pathogens or microbial metabolites activate the succinate receptor Sucnr1, Tas2rs, Vmn2r26 and prostaglandin receptors and possibly other GPCRs, which then stimulate the heterotrimeric G proteins composed of Gα-gustducin, Gβ1, Gβ3, and Gγ13. Single-cell RNAseq data have shown that many of these GPCRs and G protein subunits are co-expressed in the intestinal as well as ectopic pulmonary tuft cells although additional G protein subunits could also be involved in the signal transduction (Barr et al., 2022). Upon activation, the heterotrimeric Gαβγ dissociates into the Gα and Gβγ13 moieties, each of which can act on their own effectors. While Gα-gustducin’s effectors are still to be determined, Gβγ13 can act on the phospholipase Plcβ2 or other effectors such as adenylate cyclases, phospholipases, receptor kinases, voltage-gated ion channels and inward rectifying ion channels (Campbell and Smrcka 2018; Kankanamge et al. 2022). Among them, Plcβ2 can generate the second messengers diacylglycerol (DAG) and inositol trisphosphate (IP3) (Huang et al. 1999). DAG can activate a number of protein kinases, but their eventual physiological effect is yet to be elucidated. On the other hand, IP3 can bind to and open its channel receptor IP3R on the endoplasmic reticulum (ER), releasing the calcium ions from the ER and consequently increasing the cytosolic calcium concentration, leading to the opening of Trpm5, a calcium-activated non-selective monovalent cation channel, depolarizing the membrane potential and releasing certain output signals (Zhang et al. 2007; Liman 2014). The release of the cytokine IL-25, acetylcholine (ACh), ATP, calcitonin gene-related peptide (CGRP) and component C3 from tuft cells and taste bud cells is dependent on Trpm5, but that of some antimicrobial substances and eicosanoids does not require Trpm5 (Lee et al. 2014), indicating that Trpm5 regulates only some output signals’ release. Since the expansion of ectopic tuft cells is independent of Trpm5, we believe that Trpm5-dependent output signaling molecules may not be involved in the regulation of these tuft cells’ proliferation. On the other hand, the elimination of Gγ13-expressing pulmonary tuft cells in the Gng13-cKO mice suggests that the effectors activated by the Gβγ13 moiety are involved in the regulation of output signals that promote the proliferation of Gγ13-expressing tuft cells. The differential impact of Trpm5 knockout versus Gng13 abolishment on tuft cell proliferation is probably attributed to the fact that Gγ13 is upstream of Trpm5 in the signaling cascades, thus regulating more signaling pathways, and when mutated, having a more profound effect than Trpm5. Further investigation is needed to determine what signals the Gγ13-positive and -negative ectopic pulmonary tuft cells are needed to expand, and with this knowledge, one could intentionally manipulate increasing the clinically beneficial subtypes while decreasing the harmful ones of dysplastic tuft cells.

Previous reports indicate that elimination of all ectopic tuft cells by knocking out the Pou2f3 gene did not significantly affect the pathogenesis and outcome of H1N1-induced severe injury (Barr et al. 2022;Huang et al. 2022). In our study, elimination of Gγ13-expressing ectopic pulmonary tuft cells rendered severer symptoms following H1N1 infection, ranging from more immune cell infiltration, increased pyroptotic gene expression and cell death, larger injury areas in the lungs, more BALF protein contents and immune cells, and larger areas of fibrosis, to a slower recovery process and higher fatality rate. The striking discrepancy may be caused by the imbalance of tuft cells’ diverse functions in the Gng13-cKO lung. Since the ectopic tuft cells are absent until 12-14 dpi, these cells may not be involved in the initial inflammatory responses to viral infection or the subsequent viral clearance. When the pathogeneses enter the later stages, the Gng13-cKO lungs have a much reduced number of tuft cells with only a few of Gγ13-negative tuft cells remaining, leading to much severer symptoms, indicating that the missing tuft cells may play a major role in the transition from the inflammatory responses to inflammation resolution and subsequent tissue repair processes while the remaining tuft cells are unable to effectively contribute to these processes, or even worse, stimulate inflammation.

Indeed, some tuft cells have been shown to produce pro-inflammatory cytokines and lipid mediators such as IL-25, prostaglandins and leukotrienes that activate innate and adaptive immune responses (Kotas et al. 2022; Ualiyeva et al. 2021; Hollenhorst and Krasteva-Christ 2023), and single cell RNAseq analyses indicate that pulmonary ectopic tuft cells express genes that encode enzymes producing specialized pro-resolving mediators (SPMs), for example, Alox5, Alox12e, Alox15, Ptgs2, encoding arachidonate 5-lipoxygenase (5-LOX), arachidonate 12-lipoxygenase (12-LOX), arachidonate 15-lipoxygenase (15-LOX), prostaglandin-endoperoxide synthase II (also known as cyclooxygenase II, COX-2), respectively, as well as the genes encoding various isoforms of cytochrome P450: Cyp51, Cyp2f2, Cyp7b1 and Cyp2j6 (Barr et al. 2022). These enzymes can utilize membrane lipids such as docosahexaenoid acid, eicosapentaenoid acid and arachidonic acid to biosynthesize SPMs, including lipoxins, resolvins, protectins, and maresins, which can limit further infiltration of neutrophils, recruit and stimulate macrophage phagocytosis to remove dead cells and cell debris, suppress the production of pro-inflammatory mediators and cytokines, and increase the production of anti-inflammatory mediators and cytokines, reduce fibrosis and accelerate tissue repair and remodeling (Serhan et al. 2020; Panigrahy et al. 2021; Serhan et al. 2018). So far, the SPMs maresin 1 and resolvin D1 have been shown to inhibit virus- and bacterium-induced inflammation in the lung, respectively (Krishnamoorthy et al. 2023; Codagnone et al. 2018). Together, our data from this study imply that nullification of Gng13 not only significantly reduces the number of ectopic tuft cells but also disrupts the switch from producing proinflammatory lipid mediators to producing anti-inflammatory SPMs. Given the heterogeneity of tuft cells (Xiong et al. 2022; Bankova et al. 2018; Haber et al. 2017), we hypothesize that the remaining dysplastic tuft cells may continue to produce proinflammatory cytokines and lipid mediators whereas the anti-inflammatory agents-producing tuft cells are absent in the Gng13-cKO mice, consequently leading to severer outcomes.

In conclusion, to our knowledge, it is the first time to pustulate that different subtype of tuft cells may play different subtype-specific roles at the phases of infection, response, resolution, and recovery upon H1N1 infection or other severe injuries. Our data reveal that Gγ13-expressing ectopic tuft cells promote the inflammation resolution and recovery while other dysplastic lung tuft cells contribute to extended inflammation and perhaps immune storms as well.

Methods and Materials

Experimental design

This study was designed to investigate the possible roles of the taste signaling proteins, the heterotrimeric G protein subunit Gγ13 and the transient receptor potential ion channel Trpm5, in the host response to H1N1 influenza virus infection as well as in the subsequent inflammation resolution, tissue remodeling and recovery using transgenic and conditional gene knockout mouse models with wild-type (WT) mice as control. Animals were intranasally inoculated with H1N1 viruses at a sublethal dosage, their bodyweights were measured, mortality rates determined; the lung tissues and bronchoalveolar lavage fluid were molecularly, biochemically, immunologically, and statistically analyzed.

Animals

C57BL/6 mice were purchased from the Shanghai SLAC Laboratory animal company. ChAT-IRES-Cre, Trpm5-knockout / lac Z knockin (i.e., Trpm5−∕−), and Ai9 (Jax stock numbers 006410, 005848 and 7909) were obtained from the Jackson laboratory. The Gng13flox∕flox mice were generated previously (Li et al. 2013) and bred with the ChAT-IRES-Cre mice to conditionally knock out the Gng13 gene in the choline acetyl-transferase (ChAT)-expressing cells of the Gng13-cKO mutant mice. Mice were bred to more than 98% of C57BL/6 genetic background and maintained under a 12-hour light/dark cycle with access to water and food ad libitum in the Laboratory Animal Center of Zhejiang University. Progeny was genotyped by PCR, and both male and female mice at 6 to 8 weeks old were used in the experiments. Studies involving animals were approved by the Zhejiang University Institutional Animal Care and Use Committee, and performed following the NIH “Guidelines for the Care and Use of Laboratory Animals”.

Animal models of virally or chemically induced lung injury

Mice were anesthetized with 5% chloral hydrate (Sangon Biotech, A600288), and then intranasally administrated with 120 pfu of influenza virus A/Puerto Rico/8/1934 in 25 µl PBS at 0 day post infection (dpi). Control mice were intranasally instilled with 25 µl PBS only. The mice were checked every day and weighted every two days and euthanized for analysis at specific time points.

For chemically injured mouse models, mice were intraperitoneally administrated with a single dose of bleomycin (Selleck, s1214) at 6 mg/kg bodyweight, or first anesthetized and then injected into trachea with the house dust mite extract (Biolead, XPb70D3A2.5) at 2 mg/kg bodyweight, one dose per day for three consecutive days, or injected through trachea a single dose of lipopolysaccharide (LPS) (Sigma, L2630) at 1 mg/kg bodyweight. The lung tissues were dissected out for studies 14 days after the first injection.

Fluorescence-activated cell sorting

The heterogeneous ChAT-Cre: Ai9 mice at 25 dpi were anesthetized and transcardially perfused with cold PBS; their lungs were then isolated and dissociated into single cells as described previously (Zhao et al. 2020; Barkauskas et al. 2013). Briefly, the lung tissues were minced into pieces of 1 mm3 by a pair of ophthalmic scissors, which were then digested with 1 mg/mL dispase II (Sigma, D4693), 3 mg/mL collagenase I (Sigma, C0130) and 0.5 mg/mL DNase I (Sigma, DN25) in DMEM/F12 medium (Thermo Fisher Scientific, C11330500BT) with 1% P/S for 50 min at 37 C on a shaker. The cell suspension was then filtered by a 70-µm cell strainer (Biosharp, BS-70-XBS) and treated by red blood cell lysis buffer (Beyotime, C3702) for 3 min at RT. The cell suspension was then again filtered by a 40-µm cell strainer (Biosharp, BS-40-XBS), and the cells were collected and resuspended in the collection buffer (DPBS plus 0.04% BSA) (Figure 2—figure supplement 2). Cell sorting and data analysis were performed on a Beckman Coulter FACS flow cytometer.

Immunohistochemistry

Lungs were harvested and processed as described previously (Zhao et al. 2020). Briefly, freshly dissected mouse lungs were fixed with 4% paraformaldehyde (PFA) at RT for 2.5 h on a shaker, followed by cryoprotection in 30% sucrose in PBS at 4 C overnight. The lungs were incubated in 30% sucrose in PBS pre-mixed with an equal volume of optimal cutting temperature compound (OCT, Tissue-Tek) at RT for 2 h on a shaker, and then embedded in 100% OCT and placed in a -80 C freezer for 2 h, which were then sliced into 12 µm-thick sections on a cryostat (Thermo scientific, HM525). The tissue sections were blocked in the blocking buffer (PBS plus 3% BSA, 0.3% Triton X-100, 2% donkey serum and 0.1% sodium azide) for 1 h at RT, followed by incubation at 4C overnight with primary antibodies diluted in the blocking buffer (rabbit anti-Dclk1 1:1000, Abcam, ab31704; rabbit anti-Dclk1-488 1:1000, Abcam, ab202754; rabbit anti-Gnat3 1:500, Santa Cruz Biotechnology, sc-395; rabbit anti-Plcβ2 1:500, Santa Cruz Biotechnology, sc-206; rabbit anti-CD45 1:500, Abcam, ab10558; rabbit anti-CD64 1:200, Thermo Fisher Scientific, MA5-29706; rabbit anti-Gsdmd 1:400, Affinity, AF4012; rabbit anti-Gsdme 1:500, Proteintech, 13075-1-AP; rabbit anti-F4/80-FITC 1:500, Thermo Fisher Scientific, 11-4801-81; rabbit anti-EpCAM 1:500, Thermo Fisher Scientific, 11-5791-82; rabbit anti-krt5 1:500, BioLegend, 905504). The slides were washed and incubated with the secondary antibodies: donkey anti-rabbit IgG H&L (Alexa Fluor 488) 1:500, Abcam, ab150073; donkey anti-rabbit IgG H&L (Alexa Fluor 568) 1:500, Abcam, ab175470, for 1.5 h at RT. Fluorescent images were captured using an FV3000 laser scanning confocal microscope (Olympus) and VS200 slide scanner (Olympus).

Calcium imaging

TdTomato+ cells were isolated and FACS-sorted as described above from the H1N1-infected ChAT-Cre: Ai9 mice. These cells were seeded on laminin-coated chambers and cultured in the imaging buffer (1X HBSS, 10 mM HEPES, 1mM sodium pyruvate) for 30 min. Cells were loaded with 5 µM fluorescent Ca2+ indicator Fluo-4 AM (Dojindo, F312) and Pluronic F-127 for 40 min at 37 C, 5% CO2, following the previously reported protocol (Luo et al. 2019). The cells were gently washed twice with HBSS to remove any excess Fluo-4 AM and Pluronic F-127. Denatonium benzoate (D.B. Sigma, D5765) and quinine (MedChemExpress, HY-D0143) were used as bitter substances. Allyl isothiocyanate (AITC, Sigma, 36682), gallein (APExBIO, B7271) and U73122 (Sigma, U6756) were used as inhibitors for bitter taste receptor, G protein βγ subunit and Plcβ2 activity, respectively. The cells were first stimulated by the bitter compounds of 1 mM D.B. or 100 µM quinine, then blocked by incubating with 3 mM AITC, 100 µM gallein or 10 µM U73122 for 10 min before the inhibitors were removed; and D.B. or quinine applied again to assess the cells’ calcium responsiveness. Images were captured with an excitation wavelength of 494 nm and emission wavelength of 516 nm using an IX83 total internal reflection fluorescence microscope (Olympus).

Bronchoalveolar lavage fluid (BALF) assays

BALF were prepared as previously described (Akbari et al. 2003). Briefly, after the trachea was exposed, 0.5 mL ice-cold PBS was gently instilled into the lung and then all the fluid was collected, which was repeated 3 times. A total of 1.5 mL BALF was collected, which was centrifuged at 300x g for 10 min at 4 C. The supernatants were used to determine lactate dehydrogenase activity and IL-1β level using LDH-GloTM Cytotoxicity Assay (Promega, J2380) and Mouse IL-1 beta DuoSet ELISA (R&D Systems, DY401) following the manufacturer’s instruction. The cells in the pellet were resuspended, a portion of which was used to count cells whereas the rest was used to extract proteins using RIPA buffer (Sangon Biotech, C500005); and the protein concentrations were determined by a BCA protein assay kit (Sangon Biotech, C503021) following the manufacturer’s instruction.

Lung wet/dry weight (W/D) ratio

Mice were sacrificed with a lethal dose of 5% chloral hydrate. The lungs were isolated and cleaned to remove any adherent blood. The wet weight of lungs was immediately determined, followed by drying in a heated stove at 65 C for 48 h. The W/D ratio was calculated by dividing the wet lung weight by its corresponding dry weight.

Quantitative reverse transcription-PCR

Lungs from H1N1-infected and -uninfected mice were dissected out and washed 3 times with ice-cold PBS. The damaged lung areas of the infected mice and the corresponding lung areas of the uninfected control mice were excised to isolate total RNAs using TaKaRa MiniBEST Universal RNA Extraction Kit (TaKaRa, 9767), which were reverse transcribed into cDNAs using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, K1622). qPCR reactions were set up using iQTM SYBR Green Supermix (Bio-rad, 1708884), and run on the CFX ConnectTM Real-Time System. Relative expressions were calculated and normalized to the expression of an internal control gene Actb. The primer pairs used in these qPCR reactions are listed in Table 1.

Western blotting

Lung lobes from H1N1-infected and -uninfected control mice were harvested. The damaged lung tissues from the infected mice and the tissues from the corresponding areas of the control lungs were collected to extract proteins using RIPA buffer (Sangon Biotech, C500005), which were centrifuged at 14,000×g for 25 min to remove any debris. The protein concentrations in the supernatants were determined using a BCA protein assay kit (Sangon Biotech, C503021), and an equal amount of protein, i.e., 20 or 40 µg from each sample was loaded onto a 10-12% SDS-PAGE gel for protein separation. The proteins on the gel were then transferred to PVDF membranes (Millipore, IPVH00010), which were blocked in 5% BSA for 1 h at RT. The membranes then were incubated at 4C overnight with the primary antibodies diluted in the blocking buffer (rabbit anti-NLRP3 1:1000, Abcam, ab263899; rabbit anti-Gsdmd 1:1000, Abcam, ab219800; mouse anti-caspase-1 1:1000, Santa Cruz Biotechnology, sc-56036; goat anti-IL-1β 1:800, R&D Systems, AF-401-SP; mouse anti-β-actin 1:1000, Beyotime, AA128-1), followed by incubation with the HRP-conjugated secondary antibodies (goat anti-rabbit 1:2000, Beyotime, A0208; donkey anti-goat 1:2000, Beyotime, A0181; goat anti-mouse 1:2000, Beyotime, A0216) for 1.5 h at RT. Protein bands on the membranes were visualized using a chemiluminescent imaging system (Tanon, 5200) and quantified using the ImageJ software.

Collagen deposition assay

Lung tissue sections were analyzed with hematoxylin and eosin (H&E) staining and Masson’s trichrome staining for the assessment of the proportion of collagen content (Jia et al. 2019). Modified Masson’s Trichrome Stain Kit (absin, abs9348) was used according to the manufacturer’s instructions. The sections then were photographed using a multizoom AZ100 microscope (Nikon) and the Masson’s trichrome-stained collagen areas were determined using the ImageJ software, divided by the total area of the section. For each treatment, 3 or more mice were used.

Sytox green staining

Each mouse was intranasally administrated with 1 mL 1× Sytox Green (Beyotime, C1070S) (Halverson et al. 2015). Ten minutes later, the lungs were harvested and cut into 12 µm-thick sections as previously described (Zhao et al., 2020). Green fluorescence (Ex. 490; Em. 535) was measured using FV3000 laser scanning confocal microscope (Olympus).

Statistical analysis

Experimental data were obtained from three or more biological replicates, presented as mean ± SD using Graph Prism 9 software. Unpaired two-tailed Student’s t-tests were performed and P values <0.05 were considered statistically significant. Animal survival rates were determined using the Kaplan-Meier method and data were analyzed with the log-rank test.

Acknowledgements

We thank all Huang lab members for their helpful discussion and support.

Additional imformations

Funding

This research was in-part supported by

Starry Night Science Fund of Zhejiang University Shanghai Institute for Advanced Study Grant SN-ZJU-SIAS-005 (to LHuang) and SN-ZJU-SIAS-003 (to RZ)

National Key Research and Development Program of China Grant 2021YFF1200803 (to LHuang), 2021YFF1200404 and 2021YFA1201200 (to RZ)

National Natural Science Foundation of China U1967217 (to RZ)

National Center of Technology Innovation for Biopharmaceuticals NCTIB2022HS02010 (to RZ) National Independent Innovation Demonstration Zone Shanghai Zhangjiang Major Projects ZJZX2020014 (to RZ)

Author contributions

Conceptualization: LHuang, RZ, YHL

Methodology: LHuang, YHL, YSY, YBX, HLei, YY, SSZ

Investigation: YHL, YSY, YBX, HLei, SSZ, JQ

Visualization: YHL, YSY, YBX, JQ

Supervision: LHuang, RZ

Writing—original draft: YHL, YSY, LHuang, RZ

Writing—review & editing: LHuang, RZ

Competing interests

The authors declare that they have no competing interests.

Data and materials availability

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Additional files

This PDF file includes: Table 1, Figure 1—figure supplement 1 to Figure 5—figure supplement 2

Sequences of primers used for qPCR

Chemically induced severe injury also engenders ectopic formation of lung tuft cells.

Top three rows: Immunostaining with an antibody to Dclk1 (green) on the tissue sections of the injured lungs induced by bleomycin (BLM), extract of house dust mite (HDM) or lipopolysaccharides (LPS), respectively. Bottom row: double-immunostaining an LPS-injured lung section with antibodies to Dclk1 and Gαgust (Gα-gustducin). Scale bar: 50 µm.

Expression of taste2 receptors in tuft cell clusters.

Reanalysis of the single cell RNA sequencing data GSE197163 (Barr et al. 2022) using the Uniform manifold approximation and projection (UMAP) method. Among 613 Trpm5-GFP+ tuft cells, Tas2r108 was expressed in the greatest number of cells, i.e., 67 cells, followed by Tas2r105, Tas2R138, Tas2r137, Tas2r118 and Tas2r102, which were detected in 11,10, 10, 5 and 4 cells, respectively. The clustering analysis results also indicate that different cells may express a different subset of Tas2rs.

FACs sorting of ectopic tuft cells.

ChAT-tdTomato reporter expression in live lung epithelial cells post-influenza via FACs.

Statistical analysis of the injured tissue volumes on the lungs of WT, Trpm5−∕− and Gng13-cKO mice at 25 dpi.

(A)Sections of WT, Trpm5−∕− and Gng13-cKO mice at 25 dpi. Areas inside the dotted lines were determined as the severely damaged area. Scale bar: 0.5 mm. (B)Statistical analysis of the injured tissue volumes on the lungs of WT, Trpm5−∕− and Gng13-cKO mice at 25 dpi. The thickness of each section was 12 µm, and the area of lung injury in each section was determined using imageJ. The volume of the lung injury in each section was the area of the lung injury area multiplied by the section thickness. Each section of lung tissue was collected, and the sum of the lung injury volumes of all sections was the final volume for the whole lung. Data are presented as means ± SD (n=3), and unpaired two-tailed student t tests were performed. **P<0.01.

Ectopic generation of tuft cells in WT, Trpm5−∕− and Gng13-cKO mice over days post H1N1 infection.

(A, B, C, D, E) Immunostaining and quantitative analyses of the abundance of Dclk1+ tuft cells in the Krt5+ regions of the injured lung tissues at 14 and 25 dpi, respectively. (F) Comparative analysis of Dclk1+ tuft cell numbers in the Krt5+ regions of WT versus Gng13-cKO over days post infection. (G,H) Double immunostaining of lung tissue sections of 25 dpi with antibodies to Dclk1 and Gγ13 shows a subset of Dclk1+ tuft cells expressing Gγ13 in WT but not in Gng13-cKO mice. (I) UMAP analyses of the single-cell RNAseq dataset GSE197163 (Barr et al. 2022) indicate that out of 613 Trpm5-GFP+ tuft cells, 350, i.e., 57% of them, expressed Gγ13. (J) UMAP analyses of the single-cell RNAseq dataset GSE197163 (Barr et al. 2022) indicate that out of 350 Gγ13+ tuft cells, 277, i.e., 79% of them expressed Alox5.

qRT-PCR analysis of gasdermin gene expression in WT and Gng13-cKO lungs.

Expression of gasdermin E and caspase 3 in the H1N1-injured lungs.

Partial colocalization of gasdermin D, E and IL-1β to epithelial cells and macrophages.