Abstract

Cystic fibrosis (CF) is a genetic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) that controls chloride current. A number of different CFTR transgenic mouse lines have been developed and subjected to both acute and chronic infection models. However, prior studies showed no substantial differences in bacterial clearance between CF and non-CF mice after single inoculations. Here, using F508del transgenic CF mice, we examined the role of repeated acute Pseudomonas aeruginosa (PA) infection, with the second inoculation 7 days after the first. We found that CF mice were more susceptible to PA infection than non-CF mice following the second inoculation, with non-CF mice showing better neutrophil recruitment and effector functions. We further investigated the characteristics of lung immune cells using single-cell RNA sequencing, finding that non-CF lung neutrophils had more prominent upregulation of adhesion molecules including intercellular adhesion molecule-1 (ICAM-1) compared to CF lung neutrophils. Although people with CF are often colonized with bacteria and have high numbers of neutrophils in the airways during chronic infection, these data suggest that CF neutrophils have deficient effector functions in the setting of repeated acute infection.

1. Introduction

Cystic fibrosis (CF) is a genetic disorder caused by mutations in the gene encoding the cystic fibrosis transmembrane conductance regulator (CFTR) epithelial chloride channel [1]. Progressive pulmonary disease due to chronic infection and inflammation is the major cause of morbidity and mortality among people with CF [2]. Despite the marked improvement in lung function and overall health of people with CF in the era of CFTR modulator therapies, chronic lung infection persists, albeit with somewhat lower bacterial loads [3]. Transgenic CF KO mice were first developed to study CFTR function and CF, but they suffered from high mortality [4]. Other groups developed the so-called “gut-corrected” CF mouse strain which has a stop codon in the murine cftr gene (S489X) but also expresses human CFTR in the gut epithelium due to transgenic introduction of human Cftr under the control of the fatty acid binding protein promoter [5]. In other approaches, the F508 deletion mutation (F508del), which is the most commonly seen in people with CF, was introduced into the murine cftr gene with a lower mortality [6]. Pseudomonas aeruginosa (PA) is a major microbe that infects CF patients. Lipopolysaccharide- (LPS-) rough mucoid PA strains are seen in chronic CF infection, but the initial colonizing PA strains in people with CF are LPS-smooth and nonmucoid [7]. LPS-rough strains lack the serogroup-determining O antigen polysaccharide (also known as O side chain) attached to the LPS core. Acute infection murine models induced by a one-time infection of PA have been reported, but there was no difference in lung bacterial colony forming units (CFUs) between CF and non-CF mice in most studies [8, 9], although some did report increased PA loads in the lungs of CF mice after single inoculations with less internalization of PA in epithelial cells in the CF mice when LPS-smooth strains of PA were inoculated. In the case of LPS-rough mucoid strains, no difference was observed [10].

Chronic infection models have been reported as models to show CF susceptibility to PA infection. Coleman et al. provided PA-containing drinking water to the gut-corrected CF mice for one week and found more oropharyngeal colonization in CF mice [11]. Paroni et al. injected PA-embedded agar bead intratracheally [12]. In this case, beads may serve as a platform for biofilm type of infection. However, there was no difference in CFU in the lungs between CF and non-CF strains.

Here, we hypothesized that repeated introduction of the same colonizing bacterial strain into the lower airways via aspiration likely leads to the transition to defective bacterial clearance on repeated aspiration. Thus, we examined if repeated instillation of PA would be associated with less bacterial clearance and/or altered host response in CF mice.

2. Methods

2.1. Mouse

F508del mice (Cftrtm1/Eur) in the FVB background were provided as breeding pairs by Dr. Bob Scholte (Erasmus University Medical Center, Rotterdam, The Netherlands) and housed under specific pathogen-free condition, with 12-hour light and draft cycles [6, 13, 14]. Mice were given acidified water (sodium acetate, pH 4) to prevent acquisition of PA from the drinking water. All animal protocols were approved by the Institutional Animal Care and Use Committee (IACUC) at Boston Children’s Hospital.

2.2. Intranasal Instillation of Pseudomonas aeruginosa and Bacterial Load Analysis in the Lungs

All the experimental procedures complied with Animal Research Reporting of In Vivo Experiment (ARRIVE) guidelines regarding the use of animals in research [15] and the National Institutes of Health guideline for the care and use of laboratory animals and were approved by Boston Children’s Hospital animal care and use committee. Pseudomonas aeruginosa strain PAK (an LPS-smooth nonmucoid strain) [16] was administered intranasally (dose of  CFU in 20 μL) to mice anesthetized by intraperitoneal injection of ketamine/xylazine. Mice (both males and females) were infected at 8 weeks of age. At six hours after the instillation, mice were euthanized, and lungs were explanted and homogenized (PAx1 group). This time point was used to evaluate the role of innate immune responses as we previously performed [17]. Lung homogenates were serially diluted and plated on tryptic soy (TS) agar and incubated overnight for colony forming unit (CFU) counting. A group of mice received repeated intranasal instillation a week later. Similarly, these mice were euthanized at six hours after the second instillation (PAx2 group), and lungs were subjected to bacterial load analysis as above.

2.3. Flow Cytometry Analysis

Leukocyte counts in the lungs were analyzed using flow cytometry at time points described in the corresponding figure legend. Lungs were explanted following intravascular perfusion of PBS. Whole lungs were subjected to lung cell isolation. Briefly, lungs were digested with collagenase D and DNase I, and red blood cells were subjected to lysis using low and high osmotic NaCl. Then, cells were stained with Ly6G (1A8) for neutrophils, F4/18 (BM8) for macrophages, NK1.1 (PK136) for NK cells, CD4 (RM4-5) for CD4 cells, CD8 (53-6.7) for CD8 cells, and B220 (RA3-6B2) for B cells. Antibodies were purchased from BioLegend (San Diego, CA). BD Accuri C6 (BD Biosciences, Franklin Lakes, NJ) was used for flow cytometry analysis.

2.4. Histology

Lungs were fixed using 4% paraformaldehyde for histological analysis and subjected to hematoxylin and eosin (H&E) staining.

2.5. Reverse Transcription and Real-Time Quantitative Polymerase Chain Reaction

Following harvest at time points described in the corresponding figure legend, lungs were snap frozen and stored at -80°C until use. Total RNA was extracted using TRIzol reagent (Life Technologies) per the company’s protocol. 1 μg of RNA was subjected to reverse transcription using SuperScript III RNase reverse transcriptase (Life Technologies, San Diego, CA) to create cDNA. Real-time quantitative polymerase chain reaction (PCR) was performed using SYBR green master mix (Applied Biosystems, San Francisco, CA). We measured tumor necrosis factor- (TNF-) α, interleukin- (IL-) 1β, IL-6, IL-10, IL-12p35, IL-12p40, interferon- (IFN-) γ, and keratinocyte chemoattractant (KC). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a housekeeping gene. The sequences of primers are listed in Supplementary Table 1.

2.6. Measurement of Lung TNF-α and IL-6 Levels

Lungs were subjected to homogenization in lysis buffer. Then, supernatants were collected for TNF-α and IL-6 level measurement using ELISA kit per company protocol (R&D Systems, Minneapolis, MN).

2.7. Neutrophil Phagocytic Function Analysis

Neutrophil phagocytosis was examined using Phagotest kit (Glycotope Biotechnology, Heidelberg, Germany). Lung cells were incubated in complete RPMI 1640 on ice for 10 min, followed by the addition of FITC-labeled E. coli. Then, neutrophil suspension was kept on ice as cold control or was placed into 37°C water bath for 20 min. At the end of incubation, cells were transferred back on ice, quenched with ethylene blue, and washed with PBS. Neutrophils were probed with Ly6G. Cells were suspended in PBS/1% PFA and measured by flow cytometry using BD Accuri C6 (BD Biosciences).

2.8. Neutrophil Reactive Oxygen Species Formation Analysis

Lung cells were cultured in complete RPMI 1640 at 37°C for 30 min. Dihydrorhodamine 123 (1 μM; Sigma-Aldrich, St. Louis, MO) was added for 5 min at 37°C to lung cells. After washing with RPMI 1640, phorbol 12-myristate 13-acetate (PMA, 100 nM) was added, and the cells was incubated for additional 15 min at 37°C. After one wash, the cells were resuspended in cold PBS with 1% FBS. Neutrophils were probed with Ly6G. Reactive oxygen species (ROS) formation was detected on BD Accuri C6 (BD Biosciences).

2.9. Single-Cell RNA Sequencing
2.9.1. 3 RNA Library Preparation and Sequencing

We performed as we previously described [18]. In brief, following fast thaw in 37°C water bath, cell viability was checked using Thermo Fisher Countess II FL with NucBlue staining. Following cell viability counting, all samples below 70% viability were subjected to dead cell removal via Annexin V bead conjugation to improve quality. Single-cell RNA libraries were generated using the Chromium Single-Cell 3 kit (10x Genomics, Pleasanton, CA). The cells were counted and loaded onto the Chromium Controller. Loading was performed to target capture of ~3000 gel-in emulsions (GEMs) per sample for downstream analysis, and samples were processed through the Chromium Controller following the standard manufacturer’s specifications.

The sequencing libraries were evaluated for quality on the Agilent TapeStation (Agilent Technologies, Palo Alto, CA), and quantified using Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA). Pooled libraries were quantified using qPCR (Applied Biosystems) prior to loading onto an Illumina sequencing platform. The samples were sequenced at a configuration compatible with the recommended guidelines as outlined by 10x Genomics.

Raw sequence data (.bcl files) were converted into fastq files and demultiplexed using the 10x Genomics’ cellranger mkfastq command. Subsequent UMI and cell barcode deconvolution along with mapping to respective genome were performed using 10x Genomics’ Cell Ranger software package to generate the final digital gene expression matrices and cloupe files.

2.10. Sequenced Data Analysis

Cloupe files were merged using 10x cloud analysis system. Cell populations were annotated using canonical markers S100A8 and S100A9 (neutrophils); Ear2 and Mrc1 (macrophages); CD3d (T cells); CD79a, CD79b, Pax5, and Ms4a1 (B cells); and F3, Lamb3, and Vim (alveolar type 1/2 cells; AT cells). Differentially expressed genes (DEGs) were defined by and . Up- and downregulated DEGs were identified and subjected to KEGG pathway analysis using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) [19, 20].

2.11. Statistical Analysis

Statistical analysis was performed using PRISM9 software (GraphPad, La Jolla, CA). Detailed statistical analyses are described in the corresponding figure legends. Statistical significance was defined as .

3. Results

3.1. CF Mice Showed Susceptibility to Acute Pseudomonas aeruginosa Infection following Repeated Instillation

We first examined bacterial loads in the lungs at 6 hours after the first intranasal instillation. However, we did not see any significant difference in PA loads between CF and non-CF lungs (Figure 1(a)). However, at 6 hours after repeated instillation of PA, CF mice had higher bacterial loads in the lungs (Figure 1(b)), suggesting that CF mice were more susceptible to repeated PA infection.

3.2. CF Mice Showed Less Robust Neutrophil Recruitment after the Second PA Instillation Compared to Non-CF Mice

To understand the underlying mechanism of how CF mice were more susceptible to repeated PA infection, we examined leukocyte subset numbers in the lungs. We did not find any significant difference in the number of macrophages, NK cells, CD4 T cells, CD8 T cells, or B cells between CF and non-CF mice at various time points (Figures 2(b)2(f)). However, neutrophil numbers were significantly higher in non-CF mice after the second instillation (Figure 2(a)). This finding was compatible with our finding from the histological analysis (Figures 3(a)3(f)). H&E-stained lung specimens after the second instillation showed more neutrophil infiltration in the lungs of non-CF mice (Figure 3(f)) compared to CF lung.

3.3. CF Lungs Showed More Proinflammatory Profiles after the Second Instillation

To understand inflammatory status, we examined mRNA expression of a group of cytokines. We found that TNF-α and IL-6 mRNA levels were significantly elevated in CF lungs after the second instillation (Figure 4). This may be due to higher bacterial loads in CF lungs following the second instillation. We also determined the protein levels of TNF-α and IL-6 using lung lysates. In line with the mRNA findings, TNF-α and IL-6 levels were higher in the lungs of CF mice following the second instillation (Suppl. Figure 1).

3.4. Non-CF Lung Neutrophils Showed Better Phagocytic Function following Repeated Instillation

There were more neutrophils in the lungs of non-CF mice after the second instillation. We examined lung neutrophil phagocytic function and ROS production. We did not see any differences in ROS production and phagocytic function after the first instillation between CF and non-CF mice (Figure 5(a)). However, we found that phagocytic function was significantly enhanced in non-CF lung neutrophils after the second instillation compared to their counterpart (Figure 5(b)). Higher lung neutrophil counts with better phagocytic functions in non-CF mice are in line with our bacterial load results (Figure 1(b)). We also examined neutrophil effector functions in naïve mice. We found that ROS was enhanced in non-CF naïve neutrophils compared to CF naïve neutrophils, but not phagocytosis (Figure 5(c)). Difference in neutrophil effector function profiles among naïve, post one-time instillation, and two-time instillation might have indicated that neutrophil phenotypical changes occurred over time following PA infection. When we examined bacterial loads in naïve blood, we did not see statistical difference although there was a trend that the average CFU was higher in CF mice (Figure 5(d)). This result is in line with the bacterial load result in the lungs following one-time instillation (Figure 1(a)).

3.5. Single-Cell RNA Sequencing Results Show the Upregulation of Adhesion Molecules in Non-CF Neutrophils following Repeated Instillation

To obtain further insights into the function of lung cells, we performed single-cell RNA sequencing analysis of mouse lung cells at baseline as well as after the second instillation.

We showed cell annotations in Figure 6(a). Neutrophil transcriptomic heatmap is shown in Figure 6(b), suggesting that the post-PAx2 instillation non-CF neutrophil transcriptomic profile was distinct. We determined up- and downregulated differentially expressed genes (DEGs) (Figure 6(c)). We identified most DEGs in non-CF neutrophils post-PAx2 instillation. We also performed ontology analysis. Chemotaxis and apoptosis were annotated from both up- and downregulated DEGs in CF and non-CF mice. In contrast, cell adhesion ontology was upregulated only in non-CF lung neutrophils following repeated instillation and downregulated in CF lung neutrophils (Figures 6(d) and 6(e)). Adhesion molecule expression profile was shown in the heatmap, demonstrating that they were highly upregulated only in non-CF neutrophils post-PAx2 instillation (Figure 2(f)). For example, ICAM-1 is highly involved in neutrophil effector functions. Since intercellular adhesion molecule-1 (ICAM-1) expressing neutrophils exhibit enhanced phagocytosis [21], this result is in line with the increased phagocytosis we observed in non-CF lung neutrophils. LPS on PA primarily targets Toll-like receptor 4 (TLR4), which would induce ICAM-1 expression [21, 22]. As shown, TLR4 expression was higher in non-CF lung neutrophils following repeated PA exposure (Figure 6(g)).

We also examined transcriptomic profiles of other cell types. The heatmap for macrophages is included in Figure 7(a). We also determined up- and downregulated DEGs in macrophages (Figure 7(b)). As in the case of neutrophils, the number of DEGs was highest in non-CF macrophages post-PA instillation. Non-CF macrophages showed an upregulation of chemotaxis and cell adhesion (Figure 7(c)), while they were downregulated in CF macrophages (Figure 7(d)). However, we did not see any difference in lung macrophage number between CF and non-CF mice (Figure 2(b)). Apoptosis was downregulated in CF macrophages (Figure 7(d)), which may explain macrophage counts in the lungs. Although NK cells showed limited number of DEGs (Suppl Figure 2a), non-CF post-PA NK cells showed more distinct transcriptomic profiles than others on the heatmap (Suppl Figure 2b). Given its very limited number of DEGs, we did not perform ontology analysis for NK cells. We also performed the analysis of B cells (Suppl Figure 3a-3d), T cells (Suppl Figure 4a-4d), and AT cells (Suppl Figure 5) analyses. Based on transcriptomic profiles, there was a significant difference in post-PA instillation non-CF and CF B cells and T cells (Suppl. Fig. 3-4). For example, CF lung B cells showed significant upregulation of cell adhesion (Suppl. Figure 3c). CF lung T cells showed significant upregulation of apoptosis, mitosis, cell division, and cell cycling (Suppl. Figure 4c). CF lung AT cells showed significant upregulation of apoptosis and downregulation of inflammatory responses and immunity (Suppl. Figure 5c-5d).

4. Discussion

Here, we describe that repeated acute PA exposure worsened PA bacterial loads in the lungs of CF mice compared to non-CF mice, which was associated with less neutrophil recruitment and effector functions.

Earlier cftr KO mouse strains were developed by various investigators to recapitulate CF pathophysiology. However, their use was significantly hindered by high mortality rates associated with intestinal obstruction [4], which is less prevalent in mice homozygous for the F508 deletion mutation in the cftr gene [6]. Infection models with CF mice are divided into two categories, acute and chronic infection models. Acute lung infection models have been established by a one-time infection with PA. However, the clinical course of the acute PA lung infection was, in most studies, indistinguishable between CF and non-CF mice that share the same background [810]. An artificial bead model has been used to mimic bacterial persistence in lungs, but the model is also invasive since most investigators administer the bacteria via intratracheal instillation. In addition, no difference in lung CFUs was observed between CF and non-CF mice using the model. Acute lung infection models, however, are informative to understand how CF responds to acute infection. Clinically, people with CF succumb to chronic lung infection, which is likely initially established by repeated aspiration of colonizing bacteria from the oropharynx into the lower airways. We found that CF mice were more susceptible to the second infection, which was associated with a lower number of recruited neutrophils to the lungs as well as defective phagocytic function.

The Cftr gene is highly expressed in barrier tissues such as lung epithelial cells and gut epithelial cells [23]. The loss of CFTR function leads to an altered ion flux through the airway epithelium and ablation of mucociliary clearance. This is considered one mechanism responsible for the colonization of bacteria in the CF airway. In addition to epithelial cells, CFTR expression in phagocytes has been reported. Neutrophils express CFTR in the plasma membrane and in intracellular vesicles related to phagosomes [24]. Although chloride current is important in neutrophil physiology [25], neutrophils have chloride channels other than CFTR such as the glycine receptor [26] and voltage-gated chloride channel [27]. In fact, the role of CFTR in neutrophil functions has been conflicting in the literature. For example, ROS was attenuated in human CF neutrophils in some reports [24, 28, 29], but not in others [3032]. Similarly, neutrophil phagocytosis was impaired in CF patients in one study [33], but not in another [34]. In the current study, we have examined murine lung neutrophil ROS and phagocytosis after the first and second instillation of PA. Phagocytosis was significantly impaired only after the second instillation. These changes might have been due to the inflammatory milieu. In fact, TNF-α and IL-6 levels were significantly elevated in the CF lungs after the second instillation compared to non-CF lungs. Ontology analysis showed upregulated adhesion molecules including ICAM-1 in non-CF lung neutrophils after PA instillation. We and others showed that ICAM-1 is induced in neutrophils via LPS stimulation [21, 22]. TLR4, which senses LPS, was also upregulated in non-CF lung neutrophils after PA instillation. Because CF mice had higher bacterial loads in the lungs after the second instillation, whether this was driven by more available LPS in CF mice or CFTR affected intrinsic TLR4 signaling pathway needs to be determined. This will be a subject of future investigation.

One of the features of the CF lung is the accumulation of neutrophils. This is considered due to chronic colonization of bacteria. However, we did not extend our study to more chronic phases of colonization. Rather, we primarily examined how CF mice responded to repeated acute PA infection. We found that CF mice become more susceptible to repeated acute PA infection associated with less recruitment of neutrophils along with less neutrophil effector functions. Although CF lungs may harbor many neutrophils in a chronically colonized state, our finding indicates that they may not have proper neutrophil recruitment when they experience repeated acute infection. In addition to neutrophils, we also examined the phenotype of other leukocytes. CF lung macrophages showed downregulation of chemotaxis and adhesion. CF T cells showed apoptosis and enhanced cell cycle activity. CF B cells showed an increase in cell adhesion. The number of these cell types was comparable between CF and non-CF mice in the lungs, but it would deserve future investigation of functional analysis.

In conclusion, we described a repeated acute PA infection model in transgenic CF mice. Following repeated instillation of PA by the intranasal route, CF lungs were more susceptible to infection. Lower neutrophil numbers/function and less adhesion molecule upregulation are likely responsible for this phenotype, but the underlying regulatory mechanism of adhesion molecule expression in CF lung neutrophils needs further investigation.

Data Availability

Upon request, we would provide all the data available.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

The author (K.Y.) received fund no. R21AI158886.

Acknowledgments

We thank Dr. Bob Scholte (Netherlands) for providing the F508del mice. This work was in part supported by the CHMC Anesthesia Foundation (K.Y.).

Supplementary Materials

Supplementary 1. Supplementary Figure 1: TNF-α and IL-6 levels in the lungs. CF and non-CF mice were subjected to PA instillation twice. At 6 hours after the second instillation, lungs were obtained and subjected to homogenization followed by TNF-α and IL-6 measurement. Data were shown as of 4 samples per group. One-way ANOVA with the Bonferroni post hoc analysis was performed. . Supplementary Figure 2: single-cell RNA sequencing analysis of lung NK cells. (a) Heatmap of lung NK cell differentially expressed genes (DEGs). (b) The number of up- and downregulated DEGs. Supplementary Figure 3: single-cell RNA sequencing analysis of lung B cells. (a) Heatmap of lung B cell differentially expressed genes (DEGs). (b) The number of up- and downregulated DEGs. (c) Upregulated DEGs in CF and non-CF lung B cells were subjected to KEGG pathway analysis and presented as bubble enrichment maps. (d) Downregulated DEGs in CF and non-CF lung B cells were subjected to KEGG pathway analysis and presented as bubble enrichment maps. Supplementary Figure 4: single-cell RNA sequencing analysis of lung T cells. (a) Heatmap of lung T cell differentially expressed genes (DEGs). (b) The number of up- and downregulated DEGs. (c) Upregulated DEGs in CF and non-CF lung T cells were subjected to KEGG pathway analysis and presented as bubble enrichment maps. (d) Downregulated DEGs in CF and non-CF lung T cells were subjected to KEGG pathway analysis and presented as bubble enrichment maps. Supplementary Figure 5: single-cell RNA sequencing analysis of lung AT cells. (a) Heatmap of lung AT cell differentially expressed genes (DEGs). (b) The number of up- and downregulated DEGs. (c) Upregulated DEGs in CF and non-CF lung AT cells were subjected to KEGG pathway analysis and presented as bubble enrichment maps. (d) Downregulated DEGs in CF and non-CF lung AT cells were subjected to KEGG pathway analysis and presented as bubble enrichment maps.

Supplementary 2. Supplementary Table 1: primer sequences for RT-qPCR.