$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Differentiation
Murine airway epithelial cells successfully differentiated after culturing at an air-liquid interface with a differentiation medium for 28 days. The presence of ciliated and goblet cells was demonstrated by immunofluorescence assay of cilia marker acetylated α-Tubulin (green; Figure 3A) and the goblet cell marker Mucin5AC, respectively6 (red; Figure 3B).
Determination of CSE concentration for cell differentiation
Stimulate the isolated epithelial cells with different concentrations of CSE for 24 h. The results showed that the epithelial cells declined in cell viability when the concentration of CSE was higher than 6% (p <0.05, Figure 4A). In order to maintain cell activity under long-term stimulation, culture mediums containing 2%, 4%, and 6% CSE were used to stimulate differentiation of epithelial cells for 28 days. Furthermore, the results showed that with the increase in CSE concentration, there was no significant cell death in the epithelial cells (Figure 4B), but there was a gradual change in the number of intercellular transmembrane resistance (TEER) (Figure 4C) and exfoliated cells (Figure 4D). In addition, an overall reduction in differentiated cells and ciliated cells was noted when murine airway epithelial cell cultures were chronically exposed to CSE (Figure 4E-H). Therefore, the differentiation medium containing 6% CSE was chosen to dysregulate epithelial barrier function and reduce ciliated cell numbers without a decrease in cell viability while differentiating epithelial cells for 28 days.
Multi-omics analysis
Two groups of the cells stimulated with differential medium containing 0% or 6% CSE for 28 days were harvested. Total mRNA was extracted for transcriptomics analysis. After library preparation, sequencing, quality control, reads mapping to the reference genome, and quantification of gene expression, the differential expression of genes between the two groups was analyzed. The clusterProfiler R package was utilized to conduct an enrichment analysis of the Gene Ontology (GO) for genes that exhibited differential expression, with adjustments made to account for gene length discrepancies (Figure 5A). The Kyoto Encyclopedia of Genes and Genomes (KEGG) serves as a comprehensive database that facilitates the comprehension of complex biological system functionalities and roles at various levels, including cellular, organismal, and ecological, by integrating molecular-level data derived from extensive datasets such as those produced by genomic sequencing and high-throughput experimental methodologies (Figure 5B). For proteomics, total proteins were extracted from two groups of cells at 28 days after cultures with differential mediums containing 0% or 6% CSE. After protein quality test, salt removal, protein enzymatic hydrolysis, Tandem Mass Tags (TMT) labeling of peptides, separation of fractions, LC-MS analysis, protein identification, and protein quantitation, differentially expressed proteins between two groups can be defined. GO analysis was conducted using the interproscan program against the non-redundant protein database (Figure 6A), and KEGG was used to analyze the pathway (Figure 6B). For metabolomics, the culture supernatants were collected from two groups of airway epithelial cells at 28 days after cultures with differentiation medium containing 0% or 6% CSE. After HPLC-MS/MS analysis, identification, and quantification of metabolites, the differential metabolites between the two groups can be defined. The KEGG Database (Figure 7A), Human Metabolome Database (Figure 7B), and Lipidmaps Database (Figure 7C) were used to annotate all metabolites. The KEGG database was used to analyze the functions of these differential metabolites and metabolic pathways (Figure 7D). For epigenomics, ATAC-seq was performed as previously reported9. Nuclei were extracted from two groups of cells at 28 days after culture with differential medium containing 0% or 6% CSE, and the nuclei pellet was resuspended in the Tn5 transposase reaction mix. The transposition reaction was incubated at 37 °C for 30 min. After adding an adapter, library preparation, quality assessment, clustering, and sequencing, differential epigenomics between two groups can be defined. The GO enrichment analysis of differentially peak-related genes directly reflects the number distribution of differentially peak-related genes in GO items enriched in biological processes, cellular components, and molecular function (Figure 8A). The KEGG database was used to analyze the functions of these differential genes and pathways (Figure 8B).
Genes encoding tight junction proteins
Transcriptomic data showed that after long-term stimulation of CSE, there was decreased expression of numerous genes encoding tight junction proteins, including Cgn, Cldn2, Cldn3, Cldn8, Cldn10, Cldn20, Cldn23, Jam2 and Ocln (Figure 9A), in which reduced expression of Cldn3 and Ocln were further confirmed by Western-Blot and immunohistochemistry (IHC) (Figures 9B-E). These data suggest that exposure to CSE markedly reduces the expression of tight junction proteins of airway epithelial cells.
Cytokine and chemokine expression
In addition to the breakdown of tight connections, transcriptomics also showed that after long-term stimulation of CSE, the expression of some cytokines and chemokines, including Cxcl5, Csf3, Il1a, Il34, Ccl20, and Il33 increased, but expression of Ccl5 decreased (Figure 10A). Luminex test further confirmed that there were increased concentrations of CXCL-5, CSF-3, and IL-1α (Figure 10B-D) but decreased CCL-5 (Figure 10E) in the supernatant of murine airway epithelial cells. These results suggest that exposure to CSE selectively acts on the expression of cytokines and chemokines.

Figure 1: Growth of primary murine airway epithelial cells. (A) Newly isolated cells appear round and transparent. (B) By day 2 post-isolation, expanded cells formed into small cell islands. (C) By day 4, the formation of larger cell islands. (D) By day 6, over 90% fusion with cobblestone morphology. (4x magnification; scale bars: 250 µm). Please click here to view a larger version of this figure.

Figure 2: Immunofluorescence assay of epithelial cell markers. Cytokeratin expression in primary murine airway epithelial cells, A549 (positive control), and Raw264.7 (negative control). Nuclei stained with DAPI (63x magnification; scale bars: 40 µm). Please click here to view a larger version of this figure.

Figure 3: Immunofluorescence assay for markers of ciliated and goblet cells. The markers used were acetylated (A) α-Tubulin and (B) Mucin 5AC for confirming cilia cells and goblet cells, respectively (100 magnification; scale bars: 10 µm in panel A and 50 µm in panel B). Please click here to view a larger version of this figure.

Figure 4: Impact of CSE concentrations on epithelial cells. (A) Cell viability at 24 h at varying CSE concentrations. (B) Cell viability at 28 days. (C) TEER measurement. (D) Count of exfoliated cells. (E-H) Immunofluorescence for ciliated (green) and goblet cells (purple) at 0%, 2%, 4%, and 6% CSE concentrations (20x magnification; scale bars: 50 µm). Please click here to view a larger version of this figure.

Figure 5: GO and KEGG enrichment analyses for gene expression. (A) GO term significance in enriched biological process, cellular component, and molecular function. (B) Proportion of differentially expressed genes in KEGG pathways. Please click here to view a larger version of this figure.

Figure 6: Protein expression analysis. (A) GO enrichment bar chart showing differential protein expression. (B) KEGG pathway analysis based on differentially expressed proteins. Please click here to view a larger version of this figure.

Figure 7: Metabolite KEGG pathway analysis. (A) Annotating all metabolites using the KEGG Database. (B) Annotating all metabolites using the Human Metabolome Database. (C) Annotating all metabolites using the Lipidmaps Database. (D) The distribution of differentially expressed metabolites in various KEGG pathways. Please click here to view a larger version of this figure.

Figure 8: Differentially Peak-related gene expression analysis. (A) GO enrichment bar chart showing differential protein expression. (B) KEGG pathway analysis based on differentially expressed proteins. Please click here to view a larger version of this figure.

Figure 9: Expression of tight junction protein in murine airway epithelial cells with or without CSE stimulation. (A) Heat map of tight junction protein expression in the transcriptome. (B-E) The expression levels of (B,C) Cldn3 and (D,E) Ocln in murine airway epithelial cells were detected by Western Blot and immunohistochemistry for Cldn3 and Ocln (80 magnification; scale bars: 50 µm). Please click here to view a larger version of this figure.

Figure 10: Cytokine expression analysis. (A) Transcriptome heatmap of cytokine expression. (B-E) Luminex assay for CXCL-5, CSF-3, IL-1α, and CCL-5 levels in cells with/without CSE stimulation. Please click here to view a larger version of this figure.