Resource and Technical AdvancePulmonology Free access | 10.1172/jci.insight.99385
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by Brewington, J. in: JCI | PubMed | Google Scholar
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by Filbrandt, E. in: JCI | PubMed | Google Scholar
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by LaRosa, F. in: JCI | PubMed | Google Scholar
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by Moncivaiz, J. in: JCI | PubMed | Google Scholar
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by Ostmann, A. in: JCI | PubMed | Google Scholar
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by Strecker, L. in: JCI | PubMed | Google Scholar
1Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati Ohio, USA.
2Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati Ohio, USA.
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Find articles by Clancy, J. in: JCI | PubMed | Google Scholar
Published July 12, 2018 - More info
Recent advances in the management of cystic fibrosis (CF) target underlying defects in the CF transmembrane conductance regulator (CFTR) protein, but efficacy analyses remain limited to specific genotype–based subgroups. Patient-derived model systems may therefore aid in expanding access to these drugs. Brushed human nasal epithelial cells (HNEs) are an attractive tissue source, but it remains unclear how faithfully they recapitulate human bronchial epithelial cell (HBE) CFTR activity. We examined this gap using paired, brushed HNE/HBE samples from pediatric CF subjects with a wide variety of CFTR mutations cultured at the air-liquid interface. Growth and structural characteristics for the two cell types were similar, including differentiation into mature respiratory epithelia. In electrophysiologic analysis, no correlation was identified between nasal and bronchial cultures in baseline resistance or epithelial sodium channel (ENaC) activity. Conversely, robust correlation was demonstrated between nasal and bronchial cultures in both stimulated and inhibited CFTR activity. There was close correlation in modulator-induced change in CFTR activity, and CFTR activity in both cell types correlated with in vivo sweat chloride measurements. These data confirm that brushed HNE cell cultures recapitulate the functional CFTR characteristics of HBEs with fidelity and are therefore an appropriate noninvasive HBE surrogate for individualized CFTR analysis.
Cystic fibrosis (CF) is the most common fatal genetic disease in persons of Northern European descent, affecting approximately 70,000 people worldwide (1, 2). Early recognition of characteristic alterations in salt/ion transport led to discovery of the genetic cause of CF, mutations in the gene encoding the CF transmembrane conductance regulator (CFTR) protein (3, 4). CFTR is a traffic ATPase chloride and bicarbonate channel present in numerous tissues that primarily functions to regulate epithelial ion and fluid homeostasis (5). Loss-of-function mutations in CFTR lead to complex, multisystem disease, with the majority of morbidity and mortality occurring in the respiratory tract (2). In the CF airway, pathologic dehydration of airway surface liquid and loss of ion transport lead to chronic obstruction, infection, inflammation, and ultimately progressive lung function decline (6).
With more than 2,000 described mutations in the CFTR gene, it is readily apparent that different genetic variants drive a variety of protein consequences and therefore require different interventions for recovery of function (2, 7). Much of this understanding was derived from heterologous models expressing high levels of normal or mutant CFTR (e.g., Fisher rat thyroid cells) and in studies of primary human bronchial epithelial cells (HBEs) from lung explants obtained at lung transplantation (8–11). Through this work, many common CFTR mutations have been classified based on the protein defect (e.g., improper folding/trafficking, poor channel conductance, altered channel gating, reduced protein synthesis). Most CFTR mutations, however, remain poorly described, with efforts ongoing to fully characterize their protein defects (12).
Small-molecule therapeutics termed “modulators” targeting dysfunctional CFTR have recently been introduced and hold promise for disease modification (13, 14). There are currently 3 FDA-approved CFTR modulators, prescribed in a genotype-based fashion. The first, ivacaftor, is approved for a group of mutations representing <10% of the CF population that cause gating, conductance, or splicing defects, as well as mutations demonstrating drug response in laboratory studies. Ivacaftor serves as the current benchmark for modulator therapies; when tested in subjects with G551D CFTR (gating defect), ivacaftor produces ~10% improvement in lung function and slows the rate of lung function decline (13, 15, 16). The second available modulator, lumacaftor, is approved in combination with ivacaftor for patients homozygous for the most common CFTR mutation, F508del (trafficking defect). This combination produced smaller improvements in lung function (~3%) and a significant reduction in pulmonary exacerbation frequency (14). Finally, a third drug, tezacaftor, was recently approved by the FDA in conjunction with ivacaftor. This combination is approved as an alternative to lumacaftor/ivacaftor for patients homozygous for F508del CFTR, with similar clinical efficacy, as well as for subjects with at least one copy of a tezacaftor/ivacaftor-responsive mutation (17, 18).
While genotype-directed therapy represents a form of personalized medicine, it remains imperfect. For those harboring CFTR mutations for which there exists an FDA-approved modulator therapy, there remains a wide range of response, suggesting future opportunities for personalized optimization (14, 16). Some of this variability represents patient-specific differences in mutant CFTR activity, as previously reported in HBEs from different donors with common genotypes (19, 20). In addition, those with less-common or poorly described mutations are currently excluded from potentially beneficial therapies due to a lack of sufficient numbers for traditional clinical trials. Both problems may be overcome through personalized preclinical model systems, wherein patient tissues are used to generate laboratory models predicting modulator response to drive precision medicine (21).
Established model systems derived from CF patients to test CFTR modulator responsiveness include explant HBEs derived at the time of transplant and intestinal organoids isolated from rectal tissue (11, 22–26). The former are not feasible for individualized testing, and the latter require rectal biopsy and specialized techniques for cell isolation, propagation, and CFTR functional analysis. There is a need, therefore, to generate individualized, portable, and noninvasive model systems able to examine the cellular responses to current and emerging CFTR modulators.
In search of an easily obtained, noninvasive source of respiratory tissue, investigators have begun to examine brushed human nasal epithelial cells (bHNEs) (27–30). The nasopharynx harbors a ciliated, pseudostratified respiratory epithelium with many similarities to the bronchial epithelium, making this a feasible model representative of the lower airway. This principle underlies the use of nasal potential difference to detect CFTR function, a well-described CF diagnostic test (31). Early ex vivo studies in surgically obtained polyp specimens confirmed robust CFTR expression that is readily measurable (32). More recently, it has been demonstrated that bHNEs are easily and safely acquired in a nonsurgical fashion through cytologic brush or curettage in individuals ranging from infants to adults (33). Propagation of bHNEs in culture was initially limited by senescence and squamous transformation, which can be mitigated through use of recently developed conditional reprogramming culture (CRC) techniques (34). This method aids in preservation of basal cells and propagation over multiple passages, allowing a single sample to produce multiple differentiated cultures (35).
While it appears logical to pursue bHNEs as a model of CFTR function, it remains unclear how faithfully these cells recapitulate HBE characteristics, including electrophysiologic behavior, expression of respiratory epithelial markers, and response to CFTR modulation. In the current study, we examined whether bHNEs and brushed HBEs (bHBEs) obtained from the same subject exhibited similar characteristics relevant to CFTR functional modulation ex vivo.
Subject characteristics. Eleven pairs of bHNE and bHBE samples were expanded and tested from 10 subjects, all carrying a clinical diagnosis of CF with a wide range of disease severity. No subjects were related, and no subjects were taking CFTR modulator therapies at the time of sample collection. CFTR genotype was determined by chart review and prior clinical assays, which varied by patient; additional genetic sequencing was not performed, resulting in 4 unidentified alleles across the 10-subject cohort (subjects D, H, and J). Clinical characteristics of the studied cohort, including CFTR genotype, are summarized in Table 1, which also includes a key to identify each subject throughout the subsequent text and figures.
Clinical characteristics of brushed human nasal and bronchial epithelial (bHNE and bHBE) cell donors
Growth characteristics and culture success. In all bHNE (n = 102) and bHBE (n = 22) cultures at our institution over a 2-year period, there was a smaller initial cell number in bHNEs (0.36 × 106 ± 0.04 × 106 cells) compared with bHBEs (0.69 × 106 ± 0.15 × 106 cells; P = 0.003, Figure 1A). Despite this, total cell expansion between bHNE and bHBE cultures was equivalent (Figure 1B), as was the amount of time required for cells to achieve confluence for each subsequent passage (Figure 1C) up to passage 4.
Center-specific growth characteristics of brushed human nasal and bronchial epithelial cells (bHNEs and bHBEs). (A) Initial sample cell count for bHNEs (white circles, n = 102) and bHBEs (gray circles, n = 22). **P < 0.01 by 2-tailed t test; bars represent mean ± SEM. (B) Cumulative cell growth. Note that the passage 0 (P0) time point reflects the initial cell counts shown in A; error bars indicate SEM. (C) Time between each passage, with no significant difference between bHNE (white bars) and bHBE (gray bars) cultures at any passage. Boxes represent median (center line), 25th and 75th quartiles (box limits), and minimum/maximum (whiskers) values. Note that in B and C, the number of ongoing cultures in each tissue type dwindles with advancing passage (P1: 87/22; P2: 32/15; P3: 13/3, P4: 3/3 for bHNE/bHBE, respectively) due both to cell-specific (e.g., contamination) and protocol-specific (e.g., systematic banking at P2) factors. Using multiple t tests (with Holm-Šidák correction), no significant difference was noted between bHNE and bHBE total cell count (B) or time to passage (C) at any passage beyond P0.
Culture success (defined as air-liquid interface [ALI] insert cultures achieving baseline transepithelial resistance >50 Ω*cm2 and demonstrating measurable amiloride- and forskolin/IBMX-induced [forskolin/ 3-isobutyl-1-methylxanthine–induced] currents) was achieved in 59% of bHNE and 73% of bHBE samples (all CF samples; subsequent experiments on bHNE samples in our laboratory showed improved culture success rate of 67% [data not shown]). The primary reason for failure was bacterial or fungal contamination during the expansion phase of culture.
Structural characteristics. Immunofluorescence of bHNE and bHBE samples demonstrated key characteristics of mature respiratory epithelia, including a pseudostratified columnar structure. Ciliated cells (acetylated α-tubulin), mucus-secreting cells (muc5AC), adherens junctions (E-cadherin), and CFTR were all detected in both monolayer types. Example images from subject J are provided in Figure 2, and additional images from subjects B and F are provided in Supplemental Figure 1; supplemental material available online with this article; https://doi.org/10.1172/jci.insight.99385DS1
Brushed human nasal and bronchial epithelial cells (bHNEs and bHBEs) share similar characteristics of a mature respiratory epithelium. (A) Detection of adherens junctions (E-cadherin, green), mucus-producing cells (muc5AC, red), and ciliated cells (acetylated α-tubulin, purple) in both bHNE (left) and bHBE (right) samples. (B) Similar CFTR (green) detection between paired bHNE (left) and bHBE (right) samples. All samples were taken from subject J (Table 1) and are representative of other subjects (3 paired samples; see Supplemental Figure 2 for additional subject images), though CFTR staining varied with mutation. Scale bars: 25 μm
Ion transport. ALI cultures were mounted in Ussing chambers and studied under voltage-clamp conditions in asymmetric chloride solutions (see Methods). Baseline resistance was variable across samples and donors, with statistically significant differences observed within 4 paired samples (Figure 3A). This resulted in an aggregate trend toward increased resistance in bHBE samples that was not significant (Figure 3B; P = 0.97). ENaC activity, as measured by amiloride-inhibited current change, was also variable; there were statistically significant differences within pairs from 5 donors (Figure 3C). Based on group comparisons, there was no difference in ENaC activity between bHNEs and bHBEs (Figure 3D; P = 0.74). There was no significant correlation of bHNE and bHBE baseline resistance or ENaC activity (Figure 3, E and F, respectively). Bland-Altman analyses revealed a bias toward larger baseline resistance measurements (+67.3 Ω*cm2) in bHBE cultures, though the line of no difference fell within the 95% CI of this mean (Supplemental Figure 2, A and B). Conversely, there was a bias toward more negative (i.e., larger absolute value) amiloride-sensitive current measurements in bHNEs (–13.7 μA/cm2), with the line of no difference falling within the 95% CI of the mean (Supplemental Figure 2, C and D).
There is poor correlation between brushed human nasal epithelial cells (bHNEs and bHBEs) for baseline resistance and epithelial sodium channel (ENaC) activity. Baseline resistance values for bHNE (white bars) and bHBE (gray bars) samples for each subject (A) and aggregate data for all donors (B). ENaC activity (as measured by amiloride-inhibited short-circuit current [Isc]) in bHNE (white bars) and bHBE (gray bars) samples for each subject (C) and aggregate data for all donors (D). (E) Baseline resistance in bHNE (x axis) and bHBE (y axis) is plotted for each subject; linear regression was nonsignificant. (F) ENaC activity is similarly plotted, with no significant correlation. n = 3–6 inserts per sample and condition; individual replicates are noted by circles in A and C to demonstrate distribution. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. Two-way ANOVA with Tukey’s correction for multiple comparisons was used for A and C, while unpaired, 2-sided t tests were used for B and D. Individual sample labels correlate with donor codes shown in Table 1. Error bars in A, C, E, and F represent SEM. Boxes in B and D represent median (center line), 25th and 75th quartiles (box limits), and minimum/maximum (whiskers) values.
CFTR activity was analyzed in two fashions; first through stimulated currents with forskolin/IBMX and VX-770 (ivacaftor) (following ENaC inhibition with amiloride), and then inhibited currents with CFTR Inhibitor-172 (CFTR Inh172). Values of stimulated CFTR short-circuit current (Isc) ranged from almost no function (subject A bHBEs, 0.1 μA/cm2) to WT CFTR function (subject J2 bHNEs, 38.5 μA/cm2), with similar CFTR Inh172 Isc (Figure 4A). Stimulated CFTR currents were statistically different within 4 paired samples, as were inhibited CFTR currents. As a group, there was no difference in stimulated (P = 0.56) or inhibited (P = 0.11) CFTR currents between bHNE and bHBE samples (Figure 4B). There was a statistically significant correlation of bHNE and bHBE pairs for both stimulated (Figure 4C; P = 0.001, R2 =0.761) and inhibited (Figure 4D; P = 0.001, R2 = 0.743) CFTR currents by linear regression. Bland-Altman analyses revealed a bias toward larger stimulated CFTR Isc measurements (+1.5 μA/cm2) in bHBE cultures, with the line of no difference within the 95% CI of this mean (Supplemental Figure 3, A and B). Similarly, there was a bias toward more negative (i.e., larger absolute values) Inh712-sensitive current measurements in bHBEs (–2.6 μA/cm2), with the line of no difference falling within the 95% CI of the mean (Supplemental Figure 3, C and D).
Stimulated and inhibited CFTR activity correlate well between brushed human nasal and bronchial epithelial cells (bHNEs and bHBEs). (A) Both stimulated (forskolin/IBMX, white; VX-770, light gray) and inhibited (Inh172, dark gray) CFTR-dependent short-circuit current (Isc) for paired bHNE (solid bars) and bHBE (hatched bars) samples for all subjects. (Note the different y-axis scales used for subjects A–F to the left and subjects G–J to the right.) (B) Aggregate data for all donors. (C) Total stimulated CFTR current (forskolin/IBMX+VX-770) in bHNEs (x axis) and bHBEs (y axis) are plotted for each subject; linear regression with 95% CI (dotted lines) is shown. (D) Inhibited CFTR currents (Inh172) are similarly plotted; linear regression with 95% CI (dotted lines) is shown. n = 3–6 inserts per sample and condition; individual replicates are noted by circles in A (stimulated: forskolin/IBMX+VX-770) to demonstrate distribution. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. Two-way ANOVA with Tukey’s correction for multiple comparisons was used to compare stimulated (forskolin/IBMX+VX-770) and inhibited (Inh172) CFTR currents in A, while unpaired, 2-sided t test was used for the same comparison in B. Individual sample labels correlate with donor codes shown in Table 1. Error bars in A, C, and D represent SEM. Boxes in B represent median (center line), 25th and 75th quartiles (box limits), and minimum/maximum (whiskers) values.
Detection of CFTR modulation. Select inserts from all paired samples, regardless of mutation, were pretreated with 3 μM VX-809 (lumacaftor) for 48–72 hours prior to analysis. VX-809 incubation did not change baseline resistance in bHNEs (control, 290.4 vs. VX-809, 294.3 Ω*cm2; P = 0.97) or bHBEs (control, 351.2 vs. VX-809, 378.5 Ω*cm2; P = 0.81). Incubation with VX-809 produced varying changes in ENaC activity that were donor dependent. These ranged from a near-complete loss of amiloride-inhibited currents (subject J) to more than a doubling (subject G). There was no statistical correlation in absolute VX-809–induced change in ENaC activity in paired nasal and bronchial cells by linear regression (P = 0.057, R2 = 0.381; Figure 5A). Bland-Altman analyses revealed a bias toward more positive absolute VX-809–induced changes in amiloride-inhibited currents (+4.7 μA/cm2) in bHBE cultures, with the line of no difference within the 95% CI of this mean (Supplemental Figure 4, A and B). When analyzed as a relative VX-809–induced change (compared with absolute Isc values), there was no correlation in percent change in ENaC activity from baseline in paired nasal and bronchial cells by linear regression (P = 0.063, R2 = 0.368; Supplemental Figure 5A). Bland-Altman analysis demonstrated a small bias toward more positive relative VX-809–induced changes in amiloride-sensitive currents (+12.1 μA/cm2), with the line of no difference falling within the 95% CI of the mean (Supplemental Figure 6, A and B).
Absolute VX-809–induced changes in CFTR currents correlate well for same-subject brushed human nasal and bronchial epithelial (bHNE and bHBE) samples. (A) Absolute change in amiloride-inhibited currents in VX-809–pretreated bHNE (x axis) and bHBE (y axis) paired samples; linear regression was nonsignificant. Absolute VX-809–induced change in stimulated (forskolin/IBMX+VX-770) (B) and inhibited (Inh172) (C) CFTR currents for paired bHNE (x axis) and bHBE (y axis) samples, along with linear regression characteristics with 95% CI (dotted lines). n = 3-6 inserts per sample and condition. Individual sample labels correlate with donor codes shown in Table 1. Error bars represent SEM.
VX-809 incubation produced variable changes in both stimulated (forskolin/IBMX and VX-770) and inhibited (Inh172) CFTR currents that were donor- and CFTR mutation–dependent. These ranged from significant reductions (subject J) to >200% improvements in both stimulated and inhibited currents (subject A). The absolute change in CFTR currents in paired bHNEs and bHBEs induced by VX-809 exposure correlated strongly by linear regression for both stimulated (P < 0.0001, R2 = 0.909; Figure 5B) and inhibited (P < 0.0001, R2 = 0.934; Figure 5C) currents. Minimal bias toward more positive changes in bHBE samples for both stimulated and inhibited CFTR current was demonstrated upon Bland-Altman analysis (Supplemental Figure 4, C–F), with the line of no difference falling within the 95% CI for these means. When analyzed as a relative VX-809–induced change, the percent change in CFTR function from baseline correlated well within subjects for both stimulated (P = 0.0004, R2 = 0.808; Supplemental Figure 5B) and inhibited (P < 0.0001, R2 = 0.909; Supplemental Figure 5C) currents. Similar to analysis of the absolute change in CFTR function, Bland-Altman analysis demonstrated minor bias toward larger stimulated and inhibited relative CFTR function changes in bHBE samples, with the line of no difference falling within the 95% CI of these mean values (Supplemental Figure 6, C–F).
Correlation with sweat chloride. Stimulated CFTR-dependent Isc (forskolin/IBMX and VX-770, no VX-809 pretreatment) for both bHNEs and bHBEs was plotted against subject baseline (non-modulator-treated) sweat chloride (SC) and compared; for the single subject sampled twice, the Isc values were averaged. For both tissue types, there was a strong correlation with SC by nonlinear regression, with similar regression characteristics (bHNEs, P = 0.0008, R2 = 0.764; bHBEs, P < 0.0001, R2 = 0.866), as shown in Figure 6. These regression curves demonstrated relatively low Isc responses in both cell types above an SC of 90, quickly rising in an exponential fashion as SC fell into the 60s. These curves suggest a highly dynamic range of functional CFTR detection in bHNEs and bHBEs for subjects with baseline SC values between 60 and 90.
Similar, nonlinear inverse correlation between CFTR activity and sweat chloride (SC) in both brushed human nasal and bronchial epithelial cells (bHNEs and bHBEs). Stimulated CFTR short-circuit current (Isc; forskolin/IBMX+VX-770, y axis) is plotted against SC values (x axis) for both bHNE (white circles) and bHBE (gray circles) samples. Similar nonlinear (log-log) regression lines are plotted for bHNE (solid line) and bHBE (dotted line) samples; regression values are as indicated. For subject J (2 paired samples), Isc values from repeated samples were averaged. n = 1 SC measurement for each subject, 3–6 inserts per tissue sample per subject. Individual sample labels correlate with donor codes shown in Table 1. Error bars represent SEM.
Primary HBEs are an established model system to examine CF disease pathology and an important tool in CFTR modulator development. In this article, we present a number of key similarities between same-subject pediatric CF bHNEs and bHBEs grown at ALI, supporting the concept of bHNE-derived, personalized models for CFTR-focused investigation. This work differs from past studies in several ways. Previous functional characterization of nasal and bronchial epithelia focused on aggregate analysis, whereas the work reported herein evaluated paired, same-subject samples (29). Additionally, most work to date has compared bHNEs with HBEs harvested from lung explants, a substantially different method. We sought to minimize bias in our comparisons by using brushed cells from the upper and lower respiratory tract that were processed, expanded, differentiated, and studied with equivalent protocols in parallel and by the same staff. Finally, this work focused on pediatric subjects and describes results from a wide range of CF mutations, spanning near-WT CFTR to severe CFTR dysfunction.
In our ALI-based culture models, we noted a number of key similarities and differences between tissue types. Growth characteristics were equivalent, despite a larger initial sample in bHBEs (Figure 1); even this difference may be related to the smaller number of bHBE samples and related skewing of that data pool due to outliers. Significant variability was noted in both cell types in the time to confluence in passages 1 and 2 (Figure 1C); this variability was primarily related to initial sample size (with larger samples reaching confluence faster). The variance declined with advancing passage due to the laboratory practice of banking samples once growth slowed, thereby selecting for more rapidly growing samples to reach passages 3 and 4. Notably, we observed a higher contamination rate in bHNEs, though this gap has closed with advancing experience (at our center in the 6 months prior to writing of this manuscript, bHNE success was at 67%, while bHBE success remained at 73% for CF cultures). Improved contamination rates have followed centralization of providers procuring samples and technicians preparing media, as well as adoption of the described antimicrobial media use for all cultures (Table 2). Any difference in bHNE contamination rates versus those of bHBEs, however, was offset by the ease of bHNE collection, and indeed repeated collection if necessary, in the non-sedated patient.
Both tissue types expressed important structural characteristics and proteins consistent with classical respiratory epithelia (Figure 2 and Supplemental Figure 1). The pseudostratified columnar phenotype seen in bHBEs was generally recapitulated in bHNEs; previous studies have also demonstrated pseudostratified columnar epithelia in bHNE cultures (29). Despite minor differences, both bHNE and bHBE samples expressed key markers of a mature respiratory epithelium, including ciliated cells, mucus-producing cells, adherens junctions, apical versus basolateral polarization, and CFTR expression. The differences noted between cell types were not consistent (e.g., between donors), as demonstrated in comparisons of the 3 donor pairs in Figure 2 and Supplemental Figure 1. Prior work has used electron microscopy to demonstrate similar characteristics in pediatric CF bHNE cultures, confirming these results (36).
There was poor within-subject correlation between bHNE and bHBE monolayers for transepithelial resistance, which likely represents cell sampling variability or variance in monolayer maturation over time (Figure 3). While there was no difference in ENaC activity between bHNEs and bHBEs in aggregate, there was also no within-subject correlation across the samples (Figure 3). Variance in ENaC activity between bHNE and bHBE samples appeared random in Bland-Altman analyses, with a slight bias toward more negative measurements in bHNE samples (though the line of no difference consistently fell within the 95% CI of the mean). Ultimately, these data imply variability in the individual monolayers, as opposed to inherent tissue differences in ENaC activity, but further study will be required for confirmation. Nonetheless, these results do raise questions regarding the suitability of bHNEs to replicate bHBE behavior for studies of ENaC.
The primary focus of this study was on electrophysiologic assay of CFTR function, which correlated well within subjects for nasal and bronchial cells. Values of stimulated CFTR function in bHNEs were similar or higher than those previously reported for CF subjects, likely reflecting the wide range of CFTR function in this diverse study group, though technical items (such as media composition) in the culture process may also drive these differences (29). In our paired analyses, both stimulated and inhibited CFTR-dependent currents correlated well, supporting the use of bHNEs as a proxy for bHBEs (Figure 4, C and D). There was a trend toward larger stimulated and inhibited CFTR-dependent currents in bHBEs, though this did not reach statistical significance. This bias toward larger CFTR currents in bHBEs was again demonstrated in Bland-Altman analyses (Supplemental Figure 3), with an average increase in bHBE stimulated CFTR Isc of 1.5 μA/cm2, and an increase in the absolute value of Inh172-sensitive current of 2.6 μA/cm2. Both of these average differences, however, were nonsignificant, as the 95% CI of the mean bias lines contained the line of no difference.
Critically, CFTR modulation with VX-809 produced detectable changes in CFTR function in both tissues across a variety of CFTR mutations, even in the subjects with the smallest baseline currents. This includes functional improvement in 3 subjects (2 F508del homozygotes [subjects C and E] and 1 with a single F508del mutation and a second folding mutation, I507del [subject A]) and severe functional decline in both cell types from both samples obtained from a subject with near-WT CFTR function (G542X/unknown CFTR mutation, subject J). VX-809–induced changes in both stimulated and inhibited CFTR function correlated well across this diverse cohort, whether assessed as an absolute change in Isc (Figure 5, B and C) or as change relative to the baseline (Supplemental Figure 5, B and C). Similarly, Bland-Altman analyses demonstrated minimal bias in both absolute and relative VX-809–induced changes in CFTR activity, and the demonstrated CIs of this bias consistently included the line of no difference. Thus, these data confirm strong correlation, though not exact matching, of pharmacologic CFTR modulation between paired bHNE and bHBE samples.
As noted in Figure 5, VX-809 induced a negative change in CFTR-dependent Isc in several subjects, ranging from very mild reduction (e.g., subject I) to almost complete loss of CFTR function (subject J, with consistent results from 2 separate sample sets). Impairment of CFTR activity following VX-809 exposure has been previously reported, though not extensively studied (29). It is possible that this reflects a toxic effect of the drug, or off-target effects; however, the same drug (same stock and batch) yielded improvement or no change in the other subjects, making this less likely. Moreover, baseline (unstimulated, before change to an apical low-chloride solution) electrical resistance was not altered in the VX-809–exposed samples from subject J, suggesting cell toxicity was not involved in the negative CFTR effects. As VX-809 is believed to interact directly with CFTR, it is possible that this interaction can have negative consequences for some CFTR mutations. This possibility has significant implications, as efforts are underway to expand existing modulator drugs to new patient groups. Thus, further study is necessary to determine the mechanism of this non-F508del CFTR impairment with VX-809 and to assess if similar effects occur with other modulators.
Both bHNE and bHBE models were assessed for correlation to SC as an in vivo biomarker of CFTR function. Interestingly, models of both cell types fit a similar nonlinear (log-log) regression line when plotted against SC values (Figure 6). A previous study has suggested a linear association between nasal potential difference and bHNE ALI cultures, but there has not been a direct link to SC reported (29). This nonlinear relationship may simply represent the small sample size in our study, but SC has also been regarded as a nonlinear measure of CFTR function (37, 38). The key inflection point of both curves suggests a highly dynamic range of CFTR detection between SC values of approximately 60 and 90 mmol/l, which is also a critical window of pathologic importance. SC values of >90 mmol/l are typically associated with minimal CFTR function and a “severe” phenotype. In contrast, SC values that approach 60 mmol/l are often observed in patients with partial-function CFTR mutations and pancreatic sufficiency. Indeed, 60 mmol/l is a diagnostic threshold for CF, and highly effective CFTR modulator therapy (e.g., ivacaftor in patients with gating CFTR mutations) can reduce SC values below this threshold (39, 40). Thus, ex vivo detection of residual CFTR function in samples from subjects with SC values between 60 and 90 mmol/l reinforces the fidelity of bHNE and bHBE models to the in vivo condition.
There are several limitations of the present work. Subjects were enrolled in an effort to study a wide range of CFTR function, and all meet clinical criteria for a diagnosis of CF. This approach, however, yields a group that is heterogeneous and not fully representative of the larger CF population (with only 2 subjects homozygous for F508del CFTR). Moreover, many samples demonstrated very small CFTR responses, which may limit the precision of correlation analyses. The interpretation of these analyses is further complicated by inclusion of 4 unidentified CFTR mutations (2 subjects homozygous for G542X and 1 with no identified mutations). While this limits genotype/function comparisons, these subjects were included to cover a wide range of CFTR activity. Importantly, this heterogeneous population represents that in which patient-derived model systems may be of greatest benefit, as they are not routinely part of modulator trials. The small number of subjects and wide variability in CFTR function included in this study also preclude the broad adoption of bHNEs without additional validation, including studies to tie model system readouts to in vivo CFTR activity. Critically, it is unclear what (if any) model system results predict patient response to CFTR modulation (e.g., absolute vs. relative change in CFTR function, etc.). Additional study will be required to vet such models as preclinical predictors of treatment success or failure. Moreover, rigorous study of bHNE models will be necessary to determine whether the variance of the test is beyond its discriminating capacity for pharmacologic CFTR modulation and to establish the repeatability of the assay within a subject.
Our studies focused primarily on electrophysiologic end points, not other airway epithelial cell–specific outcomes (e.g., transcriptomics, proteomics). It is important to note that bHNE cultures may be an inadequate model of HBEs for other aspects of disease and study, especially for pathology with clinical variance between the upper and lower airways. Previous work has demonstrated similar epithelial structure, ciliary activity, and inflammatory profiles between cell types, while a prior report noted differences in transcriptional signatures between cells derived from the upper and lower airway (41–46). Our focus on electrophysiology allows for careful comparison of CFTR-dependent ion transport but does not assess bicarbonate transport (an important aspect of CFTR function). Slower-growing cultures may be predisposed to changes in ENaC or CFTR expression, confounding results; we offset this possibility by using only cultures of passage 1 (47). Additionally, our studies utilized CRC methodologies, which may alter the native programming of the cells. A previous report indicated that employing these methods in bHNE samples enriched the population for basal cells (35). We sought to control for independent variables induced by CRC by treating bHNEs and bHBEs identically after acquisition, recognizing that subtle differences in culture conditions may optimize performance of nasal versus bronchial monolayers.
In conclusion, we have demonstrated that bHNEs recapitulate many characteristics of bHBEs within subjects by completing careful parallel cell acquisition, expansion, maturation, and analysis. While the results were not identical, there is a strong correlation between CFTR function and quantification of pharmacologic CFTR modulation across the paired cell types. These data support the use of bHNEs as a noninvasive proxy for bHBEs for individualized CFTR study.
Subject enrollment. Subjects with documented CF undergoing clinically indicated bronchoscopy procedures were enrolled, and limited clinical data were collected from the electronic medical record, deidentified, and recorded, including baseline SC, lung function, age, sex, CFTR genotype, and respiratory microbiology; all values were obtained before modulator therapy, where applicable. No subjects were on CFTR modulator therapy at the time of sample acquisition. Exclusion criteria included nasal polyposis or airway friability, in which case the clinician deemed performing the cell collection procedure to be unsafe.
Cell procurement, processing, and expansion. bHNEs were obtained through curettage (rhinoprobe, Arlington Scientific Inc.) or cytology brushing (CYB-1, Medical Packaging Corp.) of the inferior turbinate, pooling cells from both nostrils into a single 15-ml conical filled with antibiotic media (see Table 2). bHBEs were obtained by bronchoscopically brushing the bilateral mainstem bronchi under direct endoscopic visualization, pooling cells from both bronchi into a single 15-ml conical filled with antibiotic media. All bronchoscopy procedures were performed with either a BF-XP160F or BF-MP160 pediatric bronchoscope and brushing with a BC-203D-2006 1.2-mm or BC-202D-2010 2.0-mm cytology brush, accordingly (Olympus Corp.).
Cell expansion was performed as previously described, with minor modifications (28, 34, 35). bHNEs and bHBEs were kept separate, with mirrored steps performed for each sample. Cells were rinsed off of the curettes/brushes with media, collected, and centrifuged at 360 g for 5 minutes. The resultant pellets were suspended in Accutase (Innovative Cell Technologies Inc.) and centrifuged at 360 g for 5 minutes. Cells were subsequently suspended in antibiotic media and seeded into cell culture dishes precoated with VitroCol (Advanced BioMatrix) and growing irradiated mouse embryonic fibroblast feeder cells (GlobalStem). bHNEs and bHBEs were maintained in dishes until confluent, with media changed daily; media was changed from antibiotic media to expansion media (see Table 2) after 5 days.
To passage confluent cell culture dishes, media was removed, and 0.1% trypsin (Sigma-Aldrich) was added for 5 minutes to facilitate detachment. This mixture was collected and centrifuged at 360 g for 5 minutes and the supernatant removed. The resultant pellet was then passaged to a new dish or onto permeable inserts for ALI culture. All presented functional experiments utilized cells of passage 1, though cells were carried forward when possible.
ALI culture. Passaged bHNEs and bHBEs were grown at ALI as previously described (11, 28). In brief, passaged cell pellets were suspended in differentiation media (see Table 2) and counted with a hemocytometer. This mixture was seeded onto Transwell-Clear permeable supports (0.33-cm2 filters, 0.4-μm pore size; Corning Inc.) precoated with type IV collagen (Sigma-Aldrich) at approximately 260,000 cells/cm2 (80,000 cells/insert). All cells were maintained in differentiation media, changing daily, removing apical media once confluent (approximately 3–4 days). Cells were maintained with basolateral media only, changing daily, for 4–5 weeks (depending on schedule) until testing.
Ion transport. Once mature, select inserts were pretreated with VX-809 (3 μM; Selleck Chemicals) for 48 hours prior to study. Inserts were removed from media and rinsed of any pretreatment drugs, then mounted in Ussing chambers (Physiologic Instruments) and studied as previously described (28). All studies were performed in an asymmetric chloride Ringer buffer (see Table 3 for buffer solutions), producing a basolateral-to-apical Cl– secretory gradient. Under voltage-clamp conditions, cells were treated apically with 100 μM amiloride to block ENaC and sodium transport. Forskolin (10 μM) and IBMX (100 μM) were added in both compartments to increase cAMP and stimulate CFTR. VX-770 (1 μM) was added apically to potentiate CFTR. Finally, CFTR Inhibitor-172 (10 μM) was added to the apical compartment to block CFTR currents. Isc and resistance were measured using Acquire and Analyze 2.3 software (Physiologic Instruments).
Whole-mount immunofluorescence. Inserts were washed in 1× PBS for 5 minutes 3 times, then fixed overnight in 4% paraformaldehyde (Electron Microscopy Sciences) in PBS. Cells were permeabilized with 1% Triton X-100 (Sigma-Aldrich) in PBS for 15 minutes, then blocked with 2.5% horse serum (Vector Laboratories) for 3 hours, all at room temperature. Inserts were washed in 1× PBS for 5 minutes 3 times, and immunostaining was performed with the antibodies and concentrations listed in Table 4, using overnight 4°C incubation for primary antibodies and 1-hour room temperature incubation for secondary antibodies. Samples were washed with 1× PBS for 5 minutes 3 times, and the cells and permeable membranes were cut away from the insert with a razor. Membranes were mounted onto slides with ProLong Gold Antifade with DAPI (Thermo Fisher Scientific), and slides were imaged on a Nikon A1 inverted confocal microscope.
Statistics. Growth analyses include all center-specific bHNE and bHBE cultures from CF subjects over a 2-year period (n = 102 bHNE, 22 bHBE samples), while all subsequent figures include only bHBE/bHNE culture pairs (n = 11 paired samples from 10 subjects). For Isc data, a minimum of 3 culture inserts were used to generate each data point; n varied from 3 to 6 inserts per condition. All data were imported into Prism software (GraphPad Software) for analysis. Unpaired, 2-tailed t tests (2 groups, with Holm-Šidák correction for multiple comparisons where applicable) or 2-way ANOVA (multiple groups, with Tukey’s multiple comparison test where applicable) were used to compare continuous data, including growth characteristics and Isc with an a priori α (P) value less than 0.05 to determine statistical significance. Mean estimates (±SEM) are presented for comparison of continuous data. Box-and-whisker plots represent minimum/maximum (whiskers), 25th and 75th percentile (box limits), and median (central line). For all correlation analyses, data from both samples from subject J were averaged to avoid unbalanced representation of the single subject. Correlation was determined by linear and nonlinear regression analysis as indicated; for linear regression, assumption of linearity was confirmed with multiple-runs tests. For all correlation analyses, agreement was tested through Bland-Altman analysis, plotting both the absolute difference between bHNE and bHBE measurements (bHNE – bHBE) and the relative difference (bHNE – bHBE)/(mean bHNE:bHBE) against the average of the bHNE and bHBE values (48). From each resulting data set, mean values, limits of agreement (±1.96 SD), and 95% CIs of the mean were calculated and plotted.
Study approval. All studies were approved by the Cincinnati Children’s Hospital Medical Center Institutional Review Board. All subjects/caregivers provided written assent/consent prior to inclusion.
JJB and JPC conceived and designed the study. JJB, ETF, FJL, JDM, AJO, and LMS executed the study. JJB, ETF, FJL, JDM, AJO, LMS, and JPC analyzed and interpreted data. JJB and JPC drafted the manuscript for intellectual content.
The authors wish to thank Kristina Ray for her assistance in recruitment and records management and the Cincinnati Children’s Hospital Cystic Fibrosis Center for their collaboration in this work. This work was supported by Cystic Fibrosis Foundation Therapeutics (CLANCY14XX0) and by the Cystic Fibrosis Foundation (CLANCY15R0).
Address correspondence to: John J. Brewington, 3333 Burnet Avenue/MLC 2021, Cincinnati Ohio, USA. Phone: 513.636.6771; Email: john.brewington@cchmc.org.
Conflict of interest: The authors have declared that no conflict of interest exists.
Reference information: JCI Insight. 2018;3(13): e99385. https://doi.org/10.1172/jci.insight.99385.