Data from: Recapitulating apicobasal tissue polarity in extracellular matrix-incorporated airway organoids
Data files
Apr 30, 2026 version files 533.93 MB
-
Figure_1.zip
60.20 MB
-
Figure_2.zip
98.34 MB
-
Figure_3.zip
14.88 MB
-
Figure_4.zip
2.84 MB
-
Figure_5.zip
132.55 MB
-
Figure_6.zip
138.92 MB
-
README.md
21.06 KB
-
Supplementary_figures.zip
86.19 MB
Abstract
The airway epithelium is a dynamic barrier that interfaces with the external environment and internal matrix niche along its apicobasal axis. To recapitulate this tissue arrangement in an organoid format, we present the decellularized ExtraCellular Matrix-incorporated Apical-out Airway Organoid (dECM-AoAO) that integrates native matrix cues, through incorporation of human lung-derived dECM microparticles, without compromising the apical-out polarity. Incorporation of dECM effectively diversifies lineage distribution that better recapitulates native epithelial composition compared to ECM-free AoAO. Harnessing the dECM-AoAO locomotion powered by outward-facing ciliary beating, we developed an experimental and computational pipeline for swarm analysis of organoid group motility as a functional readout of ciliary function. Lastly, dECM-AoAO withstood cryopreservation and reviving with sustained viability, lineage distribution, and ciliary function, enabling future scalability and broad distribution. Together, this work establishes dECM-AoAO as a physiologically relevant and versatile model system for investigating epithelial-ECM interaction during airway homeostasis, disease pathogenesis, and injury responses.
Dataset DOI: 10.5061/dryad.5tb2rbpgq
Description of the data and file structure
This dataset contains the raw and processed image data, quantification spreadsheets, videos, and analysis code used to generate the main and supplementary figures for the dECM-AoAO Biomaterials paper. The data are organized by figure ZIP archive and internal panel folder to support reproducibility.
Note on figure order: The ZIP archives are listed below in the uploaded order: Figure 1.zip through Figure 6.zip, followed by Supplementary figures.zip. Some internal panel folders retain older manuscript panel labels, such as 2B 2C, 3A, or 4A. These folder labels are kept exactly as they appear in the dataset, so the README can be used to locate files directly.
File: Figure_1.zip
Description: This archive contains raw imaging data, analysis scripts, and quantified results for dECM microparticle labeling, fluorescence profiling, dose-dependent organoid assembly, and organoid morphometry. Internal folders are labeled as panels 2B-2I.
2B, 2C: GFP/Cy5 particle imaging
Files:
- UF-birghtfield.tif: Brightfield image of unfiltered dECM microparticles. The filename is kept as it appears in the archive.
- UF-Cy5.tif: Cy5 fluorescence channel image of unfiltered particles.
- UF-GFP.tif: GFP/autofluorescence channel image of unfiltered particles.
- particle GFP Cy5 profile.prism: Fluorescence intensity profile data for GFP and Cy5 channels.
Purpose: Visualize unfiltered dECM microparticles and compare intrinsic green autofluorescence with Cy5 labeling.
2D: GFP/Cy5 co-localization
Files:
- GFP_Cy5 particle images/: Raw GFP and Cy5 particle images from donor subfolders.
- particle Cy5 GFP coloc.cpproj: CellProfiler project used for co-localization analysis.
- Cy5-GFP coloc.prism: Processed co-localization results.
Purpose: Quantify the overlap between GFP/autofluorescence and Cy5 signals within dECM microparticles.
2E: Particle size distribution, unfiltered particles
Files:
- Unfiltered particle distribution.R: R script used to generate particle size distribution plots.
- Unfiltered particle size.xlsx: Spreadsheet containing unfiltered particle size measurements.
Purpose: Determine the baseline size distribution and morphology of unfiltered dECM microparticles.
2F: Organoid assembly with different dECM-MP doses
Files:
- TIFF images of organoids assembled with different dECM-MP concentrations in brightfield, GFP, and Cy5 channels.
- Naming examples: 10x-1000Basal-brightfield.tif, 10x-1000Basal+1ug_mL_brightfield.tif, 10x-1000Basal+10ug_mL_GFP.tif, 10x-1000Basal+100ug_mL_Cy5.tif.
- Doses represented: 0 ug/mL, 1 ug/mL, 10 ug/mL, and 100 ug/mL.
Purpose: Evaluate how dECM-MP dose affects organoid assembly, morphology, and particle incorporation.
2G, 2H: Brightfield imaging and morphometry
Files:
- 2G 2H Brightfield images/: Raw brightfield images used for organoid morphometric analysis.
- 2G 2H size+circularity.prism: Processed size and circularity data.
Purpose: Compare organoid area and circularity across dECM-MP dose conditions.
2I: E-cadherin immunostaining
Files:
- 0ug_mL dECM_AoAO-E-cad.tif
- 1ug_mL dECM_AoAO-E-cad.tif
- 10ug_mL dECM_AoAO-E-cad.tif
Purpose: Assess epithelial organization and dECM-MP incorporation in organoids generated with different particle doses.
File: Figure_2.zip
Description: This archive contains raw images, scripts, and processed data for dECM-MP filtration, particle size analysis, filtered particle count optimization, organoid morphometry, and E-cadherin intensity profiling. Internal folders are labeled as panels 3A-3K.
3A: Particle morphology before and after filtration
Files:
- donor1_before.tif: Brightfield image of donor 1 dECM particles before filtration.
- donor1_filtered.jpg: Brightfield image of donor 1 dECM particles after 40 µm filtration.
Purpose: Visually compare dECM particle morphology before and after filtration.
3B-3D: Particle size distribution and statistical comparison
Files:
- 3B particle distribution histogram.R: R script used to generate particle size histograms.
- 3C 3D particle comparison (fit log-normal distribution).R: R script used for log-normal distribution fitting and particle size comparison.
- 3C 3D.prism: GraphPad Prism file containing processed particle size analysis.
- donor1 particle size.xlsx: Particle size measurements for donor 1, including area, perimeter, and circularity.
Purpose: Quantitatively analyze the size distribution and filtration-mediated size reduction of dECM microparticles.
3E: Organoid assembly with varying filtered particle counts
Files:
- Representative brightfield images of organoids assembled with different filtered dECM-MP counts per 1,000 basal epithelial cells.
- File examples: 1000basal.tif, filtered-100+1000basal.tif, filtered-1k+1000basal.tif, filtered 10K+1000basal.tif, filtered-100K+1000Basal.tif.
Purpose: Assess the effect of filtered particle number on organoid assembly and morphology.
3F, 3G: Organoid morphometry by particle count
Files:
- Brightfield images used for analysis/: Raw brightfield images used for morphometry.
- diff counts dECM-MP formed organoid circularity+area.prism: Processed organoid area and circularity data across particle count conditions.
Purpose: Quantitatively compare organoid size and circularity across different dECM-MP loading levels.
3H: Filtration outcome
Files:
- before filtration.jpg: Representative brightfield image of particles before filtration.
- after filtration.jpg: Representative brightfield image of particles after filtration.
Purpose: Provide a visual demonstration of the filtration effect on particle size and appearance.
3I: Circularity comparison between filtered and unfiltered particle conditions
Files:
- before filtration brightfiled images/: Raw brightfield images of organoids formed with unfiltered particles. The folder name is kept as it appears in the archive.
- post-filtration brightfield images/: Raw brightfield images of organoids formed with filtered particles.
- filterd vs unfilted circularity comparison.prism: Processed circularity comparison file. The filename is kept as it appears in the archive.
Purpose: Compare organoid circularity generated from filtered versus unfiltered dECM-MPs.
3J, 3K: E-cadherin staining and longest-axis intensity profiling
Files:
- C1-filtered dECM_AoAO-E-cad-1.tif, C2-filtered dECM_AoAO-E-cad-1.tif, C3-filtered dECM_AoAO-E-cad-1.tif: E-cadherin immunofluorescence images.
- Composite_Longest_Axis.tif: Composite image used for longest-axis profiling.
- longest axis intensity.mlx: MATLAB script used for intensity profile extraction.
- Longest_Axis_Intensity_Profile.csv: Output intensity profile values.
Purpose: Assess the spatial distribution of epithelium and dECM-MPs across the organoid longest axis.
File: Figure_3.zip
Description: This archive contains raw immunofluorescence images, quantification files, CellProfiler projects, MATLAB code, and a Prism dataset for epithelial lineage characterization and ciliation analysis in AoAOs and dECM-AoAOs. Internal folders are labeled as panels 4A-4E.
4A: Immunofluorescence staining of airway epithelial markers
Files:
- AoAO images: AoAO_Ac-tub.tif, AoAO_MUC5AC.tif, AoAO_TP63.tif, AoAO_UG.tif.
- dECM-AoAO images: dECM-AoAO_Ac-tub.tif, dECM-AoAO_MUC5AC.tif, dECM-AoAO_TP63.tif, dECM-AoAO_UG.tif.
Purpose: Compare epithelial lineage composition between ECM-free AoAOs and dECM-AoAOs using markers for ciliated cells, goblet cells, basal cells, and club cells.
4B: UG/CCSP quantification
Files:
- dECM_AoAO UG quantification.xlsx: Quantified UG/CCSP-positive signal normalized to DAPI counts.
- dECM_UG.cpproj: CellProfiler project used for UG/CCSP analysis.
Purpose: Quantify club cell marker expression in AoAO and dECM-AoAO conditions.
4C: MUC5AC quantification
Files:
- Images/: Raw image slices used for MUC5AC quantification.
- dECM_MUC5.cpproj: CellProfiler project used for MUC5AC analysis.
- dECM_MUC5AC quantification.xlsx: Quantified MUC5AC-positive signal normalized to DAPI counts.
Purpose: Quantify goblet cell marker expression in AoAO and dECM-AoAO conditions.
4D: TP63 quantification
Files:
- images/: Raw image slices used for TP63 analysis.
- TP63.cpproj: CellProfiler project used for TP63 quantification.
- dECM_TP63 quantification.xlsx: Quantified TP63-positive basal cell abundance normalized to DAPI counts.
Purpose: Quantify basal cell retention in AoAO and dECM-AoAO conditions.
4E: Ciliation analysis
Files:
- images/: Raw image slices containing acetylated alpha-tubulin and DAPI channels.
- Ciliation_Analyzer.m: MATLAB script used to calculate percentage of organoid surface covered by cilia.
- dECM_AoAO ciliation.xlsx: Quantified ciliated surface coverage.
Purpose: Quantify motile cilia coverage in AoAOs and dECM-AoAOs.
Combined processed dataset
File:
- d21 dECM_AoAO differentiation.prism
Purpose: GraphPad Prism file containing processed datasets for the day 21 differentiation and lineage characterization experiments.
File: Figure_4.zip
Description: This archive contains representative immunofluorescence images and processed analysis data for IL-13 treatment of dECM-AoAOs.
4B: Representative images of control and IL-13-treated dECM-AoAOs
Files:
- Control images: control-Ac-tub.tiff, control-MUC5.tiff, control-TP63.tiff, control-UG.tiff.
- IL-13-treated images: 5ng_mL IL13-Ac-tub.tif, 5ng_mL IL13-MUC5.tif, 5ng_mL IL13-TP63.tif, 5ng_mL IL13-UG.tif.
Purpose: Compare epithelial marker expression between untreated control dECM-AoAOs and dECM-AoAOs treated with 5 ng/mL IL-13.
4C-4F: IL-13 quantification and processed analysis
Files:
- IL13 dECM-AoAO.prism: GraphPad Prism file containing processed datasets and analyses for IL-13-treated dECM-AoAOs.
Purpose: Quantify IL-13-associated changes in epithelial differentiation and ciliation-related readouts in dECM-AoAOs.
File: Figure_5.zip
Description: This archive contains representative images, raw videos, MATLAB tracking outputs, and velocity quantification files for dECM-AoAO motility assays under control conditions and after treatment with the dynein inhibitor EHNA at 1 mM.
5C: Representative images and motility videos
Files:
- images/: Representative still images of organoids from three donors under control and EHNA-treated conditions.
- Control: H61 control.png, H62 control.png, H63 control.png.
- EHNA-treated: H61 EHNA.png, H62 EHNA.png, H63 EHNA.png.
- Videos and analysis (with and without 1mM EHNA treatment)/: Donor-specific subfolders containing raw .mp4 videos and .mat MATLAB tracking files for control and EHNA-treated organoids.
- Example files: 0.0 mM.mp4, 0.0 mM.mat, 1.0 mM.mp4, 1.0 mM.mat.
Purpose: Visually and quantitatively compare cilia-driven dECM-AoAO motility before and after EHNA-mediated dynein inhibition across three donors.
5D: Velocity quantification
Files:
- velocity.xlsx: Spreadsheet containing mean velocity values in µm/s for individual organoids and group-level comparisons.
- velocity.prism: GraphPad Prism file containing velocity datasets and plotted/statistical analyses.
Purpose: Quantify dECM-AoAO motility velocity under control and EHNA-treated conditions.
File: Figure_6.zip
Description: This archive contains raw images, processed image slices, quantitative data, videos, and analysis files used to assess cryopreservation and revival of dECM-AoAOs, including viability, epithelial lineage preservation, morphology, and ciliary motility.
6B: Live/dead viability imaging
Files:
- pre-freeze live dead.tif
- 1d post-thaw live dead.tif
- 3d post-thaw live dead.tif
- 7d post-thaw live dead.tiff
Purpose: Show representative live/dead staining of dECM-AoAOs before cryopreservation and at 1, 3, and 7 days after thawing.
6C: Live/dead quantification
Files:
- Images/1d, Images/3d, Images/7d, Images/before freezing: Raw or processed image slices used for quantification.
- quantitive data/before freeze.xlsx, quantitive data/1d post-thaw.xlsx, quantitive data/d3 post-thaw.xlsx, quantitive data/d7 post thaw.xlsx: Live/dead quantification spreadsheets. The folder name is kept as it appears in the archive.
- quantitive data/live dead stain.prism: Processed live/dead viability dataset.
- live:dead intensity measurement.cpproj: CellProfiler project used for live/dead intensity measurement.
Purpose: Quantify live and dead area fractions or intensity-based viability readouts before and after cryopreservation.
6D: Post-thaw epithelial lineage marker imaging
Files:
- post-thawed dECM-AoAO-Ac-tub-7d after thawing.tif: Acetylated alpha-tubulin staining for ciliated cells.
- post-thawed dECM-AoAO-MUC5-7d after thawing.tif: MUC5AC/MUC5 staining for goblet cells.
- post-thawed dECM-AoAO-TP63-7d after thawing.tif: TP63 staining for basal cells.
- post-thawed dECM-AoAO-UG-7d after thawing.tif: UG/CCSP staining for club cells.
Purpose: Assess whether epithelial lineage markers are maintained after thawing.
6E-6H: Quantification of epithelial lineages and ciliation after thawing
Files:
- images/Ac-tub slices/, images/MUC5 slices/, images/TP63 slices/, images/UG slices/: Cropped image slices used for marker-specific quantification.
- quantitive data/Ciliation_Results.xlsx: Ciliated surface coverage quantification.
- quantitive data/MUC5AC.xlsx: MUC5AC/MUC5-positive signal normalized to DAPI counts.
- quantitive data/TP63.xlsx: TP63-positive cell counts normalized to DAPI counts.
- quantitive data/UG.xlsx: UG/CCSP-positive signal normalized to DAPI counts.
- quantitive data/post-thaw dECM_AoAO characterization.prism: GraphPad Prism file containing the processed post-thaw characterization dataset.
Purpose: Quantitatively evaluate epithelial lineage preservation and ciliation after cryopreservation.
6I: Representative brightfield images before and after thawing
Files:
- pre-thaw.png
- 1d post-thaw.png
- 3d post-thaw.png
- 7d post-thaw.png
Purpose: Compare gross organoid morphology before freezing and during post-thaw recovery.
6J: Motility analysis before and after thawing
Files:
- Videos+analysis/Pre-Thaw/: Pre-thaw raw motility video and MATLAB tracking data.
- Videos+analysis/Post-Thaw D1/, Videos+analysis/Post-Thaw D3/, Videos+analysis/Post-Thaw D7/: Post-thaw videos and MATLAB tracking outputs.
- Videos+analysis/dECM-AoAO velocity before and after thawing.xlsx: Velocity measurements in µm/s from swarm analysis.
Purpose: Quantify cilia-driven dECM-AoAO motility before freezing and at 1, 3, and 7 days after thawing.
File: Supplementary_figures.zip
Description: This archive contains raw data, processed images, analysis files, and representative images used to generate supplementary figures. The archive includes folders for Figure S2-S8 and Figure S10-S13.
Figure S2: Green autofluorescence from dECM particles in AoAOs
Files:
- 10x-1000Basal-d1-brightfield.tif: Brightfield image of ECM-free AoAOs at day 1.
- 10x-1000Basal-d1-GFP.tif: GFP channel image of ECM-free AoAOs.
- 10x-dECM-AoAO-d1-brightfield.tif: Brightfield image of dECM-AoAOs at day 1.
- 10x-dECM-AoAO-d1-GFP.tif: GFP channel image showing green autofluorescent dECM particles.
Purpose: Compare ECM-free AoAOs and dECM-incorporated AoAOs, highlighting dECM-MP-specific green autofluorescence.
Figure S3: Particle size reduction by vacuum filtration in additional donors
Files:
- images/: Brightfield images before and after filtration for donor 2 and donor 3.
- particle size+analysis/H62 size.xlsx and particle size+analysis/H63 size.xlsx: Particle size measurements for additional donors.
- particle size+analysis/particle distribution histogram.R: R script used to plot particle size distributions.
Purpose: Validate particle size reduction after filtration in additional dECM donors.
Figure S4: Dose-dependent dECM incorporation
Files:
- Analysis/diff dECM-MP AoAO GFP intensity.prism: Processed GFP intensity dataset.
- Analysis/GFP intensity diff counts dECM-AoAO.xlsx: Quantified GFP intensity values per organoid.
- images used for analysis/: Raw images used for GFP intensity quantification.
- representiative images/: Representative images for dose-dependent dECM incorporation. The folder name is kept as it appears in the archive.
Purpose: Quantify dose-dependent dECM-MP incorporation using GFP/autofluorescence intensity.
Figure S5: Effect of particle filtration in additional donors
Files:
Raw TIFF images of donor 2 and donor 3 organoids formed with filtered or unfiltered dECM-MPs at 1,000 particles per organoid.
Example filenames include 10x-donor2-1000Basal+1000-1d...tif, 10xdonor2-1000Basal+1000 unfiltered...tif, 10x-donor3-1000Basal+1000-1d...tif, and 10xdonor3-1000Basal+1000 unfiltered...tif.
Purpose: Assess whether filtration improves organoid formation and morphology across additional dECM donors.
Figure S6: dECM-AiAO morphology and epithelial marker staining
Files:
- 40x-day 7.JPG and 40x-day 28.JPG: Representative morphology images at day 7 and day 28.
- dECM-AiAO-Ac-tub.tif: Acetylated alpha-tubulin staining.
- dECM-AiAO-MUC5.tif: MUC5/MUC5AC staining.
- dECM-AiAO-TP63.tif: TP63 staining.
- dECM-AiAO-UG.tif: UG/CCSP staining.
Purpose: Characterize morphology and epithelial lineage marker expression in dECM-AiAOs.
Figure S7: FOXJ1 expression analysis
Files:
- Analysis/FOXJ1.cpproj: CellProfiler project used for FOXJ1 quantification.
- Analysis/FOXJ1 ratio.xlsx: FOXJ1-positive nuclei counts normalized to DAPI counts.
- Analysis/FOXJ.prism: Processed FOXJ1 dataset.
- images used for analysis/: Raw image slices containing FOXJ1 and DAPI channels.
- representative images/: Representative FOXJ1 immunofluorescence images.
Purpose: Quantify FOXJ1 expression and compare it with ciliation-related readouts.
Figure S8: dECM-MP localization and intermixing analysis
Files:
representative images/40x-biotin-dECM-AoAO E-cad d1.tif and representative images/40x-biotin-dECM-AoAO-E-cad-d21.tif: Representative E-cadherin and biotin-dECM-AoAO images at day 1 and day 21.
- images used for analysis/: Raw images used for localization or intermixing analysis.
- codes and output/PCC.csv: Quantified output values.
- codes and output/d1 vs d21 dECM-MP location.prism: Processed day 1 versus day 21 localization dataset.
- codes and output/intermix_folder_minimal.m and codes and output/run_intermix_Index_function.mlx: MATLAB scripts used for intermixing/localization analysis.
Purpose: Quantify dECM-MP localization or intermixing within dECM-AoAOs over time.
Figure S10: Heat-treated or inert dECM-MP control analysis
Files:
- S10 BC/10142025 H63 80C heat/: Representative images of control and 80C heat-treated particle conditions.
- S10 BC/inert particle size comparison.prism: Processed particle size comparison dataset.
- S10 D/: Representative epithelial marker images for control and 80C heat-treated conditions.
- S10 E-H/inert particles-AoAO.prism: Processed lineage and ciliation analysis dataset.
- S10 E-H/slices for analysis/: Image slices used for marker quantification.
Purpose: Evaluate whether heat-treated or inert particle conditions affect organoid morphology, epithelial differentiation, and ciliation readouts.
Figure S11: IL-13 treatment in ECM-free AoAOs
Files:
- S11 B/: Representative Ac-tub and MUC5/MUC5AC images from control and 5 ng/mL IL-13-treated AoAOs.
- S11 C-D/IL13 dECM-AoAO.prism: Processed IL-13 analysis dataset.
- S11 C-D/slices for analysis/: Raw image slices used for Ac-tub and MUC5/MUC5AC quantification.
Purpose: Compare IL-13-induced epithelial changes in ECM-free AoAOs.
Figure S12: Cryopreservation analysis of LN-AoAOs
Files:
- S12 A-B/AoAO-pre-freeze-1d.tif and S12 A-B/AoAO-1d LN-post-thaw d7-live-dead.tif: Representative live/dead images before freezing and after thawing.
- S12 A-B/slices for analysis/: Image slices used for live/dead quantification.
- S12 A-B/live dead stain.prism: Processed viability dataset.
- S12 C-D/1d-LN-AoAO pre-freeze.tif and S12 C-D/1d-LN-AoAO post-thaw 7d.tif: Representative images for pre-freeze and post-thaw comparison.
- S12 C-D/Cryo-motility Analysis.prism: Processed motility analysis dataset.
Purpose: Assess viability and motility-related outcomes of LN-AoAOs before and after cryopreservation.
Figure S13: Morphology comparison before and after thawing
Files:
- representative images/AoAO-pre-freeze.tif
- representative images/dECM-AoAO-pre-freeze.tif
- representative images/post thaw 1d-AoAO.tif
- representative images/post thaw 1d-dECM-AoAO.tif
- LN_dECM-AoAO circularity+size.prism
Purpose: Compare representative morphology and processed circularity/size measurements for AoAO and dECM-AoAO conditions before and after thawing.
