Delineating the early transcriptional specification of the mammalian trachea and esophagus

  1. Akela Kuwahara
  2. Ace E Lewis
  3. Coohleen Coombes
  4. Fang-Shiuan Leung
  5. Michelle Percharde
  6. Jeffrey O Bush  Is a corresponding author
  1. Program in Craniofacial Biology, University of California San Francisco, United States
  2. Department of Cell and Tissue Biology, University of California San Francisco, United States
  3. Institute for Human Genetics, University of California San Francisco, United States
  4. Eli and Edythe Broad Center of Regeneration Medicine and Stem Cell Research, University of California San Francisco, United States
  5. Developmental and Stem Cell Biology Graduate Program, University of California San Francisco, United States
  6. Department of Biology, San Francisco State University, United States
  7. MRC London Institute of Medical Sciences (LMS), United Kingdom
  8. Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, United Kingdom
7 figures, 1 table and 1 additional file

Figures

Figure 1 with 5 supplements
Single-cell transcriptomics of the developing foregut epithelium identifies distinct dorsoventral populations.

Dissected, FACS-purified foregut epithelial cells were subjected to droplet-based single-cell RNA sequencing at E10.5 and E11.5. (a) UMAP representation of 6,407 cells identified at E10.5 and (b) …

Figure 1—source data 1

Top 20 markers of foregut cell populations identified by scRNA-seq of E10.5 and E11.5 foreguts.

https://cdn.elifesciences.org/articles/55526/elife-55526-fig1-data1-v1.xlsx
Figure 1—figure supplement 1
Quality control metrics and characterization of single-cell RNA sequencing experiments.

(a) Representative gating strategy for FACS of mouse embryonic foreguts. Top row shows gating for single cells, and bottom row shows gating for sytox-negative live cells (left) and epCAM-positive …

Figure 1—figure supplement 2
Additional cell type-specific markers of dorsoventral populations identified by scRNA-seq at E10.5.

(a-m) RNA localization of additional marker genes Wnt7b, Pcdh10, Ly6h, and Klf5 identified from differential expression analysis of cell populations in scRNA-seq data at E10.5. First column shows …

Figure 1—figure supplement 3
Cell type-specific markers of trachea, esophagus, and lung populations identified by scRNA-seq at E11.5.

(a) RNA localization of marker genes identified from differential expression analysis of cell populations in scRNA-seq data at E11.5 with immunofluorescent staining of NKX2-1. First three columns …

Figure 1—figure supplement 4
Newly identified esophageal marker genes are more dorsally restricted than SOX2.

RNA localization of esophageal markers genes Klf5, Krt19, and Dcn and protein localization of LRIG1 (green) with immunofluorescent staining of NKX2-1 (magenta) and SOX2 (cyan) in the unseparated …

Figure 1—figure supplement 5
Validation of pharynx and UBB cluster identities.

(a) RNA localization of pharynx marker gene Foxe1 identified from differential expression analysis of cell populations in scRNA-seq data at E10.5 and E11.5. Top row: RNA-localization of Foxe1

Figure 2 with 1 supplement
Transcriptomic analysis of Nkx2-1-/- foreguts reveals NKX2-1 regulated tracheal and esophageal transcriptional programs.

(a) Dissected, FACS-purified epithelium from E11.5 Nkx2-1-/- and WT foreguts was sequenced and analyzed for differential gene expression. 109 genes were increased and 148 genes were decreased in Nkx2…

Figure 2—source data 1

Differentially expressed genes between E11.5 Nkx2-1-/- and WT foregut epithelium identified with bulk RNA-seq.

https://cdn.elifesciences.org/articles/55526/elife-55526-fig2-data1-v1.xlsx
Figure 2—figure supplement 1
NKX2-1-dependent gene expression in E11.5 scRNA-seq.

(a-b) Combined expression of all genes that increase (a) or decrease (b) in Nkx2-1-/- foreguts as a percentage of all reads in E11.5 scRNA-seq data, projected onto the E11.5 UMAP.

Figure 3 with 1 supplement
Identification of an NKX2-1-independent tracheal and esophageal program.

(a-l) RNA localization of NKX2-1-independent tracheal (Irx2, Ly6h) and esophageal (Dcn, Ackr3) makers identified by scRNA-seq and Nkx2-1-/- mutant RNA-seq data at E11.5. First column shows …

Figure 3—source data 1

Differentially expressed genes between E11.5 WT trachea and esophagus epithelium identified with bulk RNA-seq.

https://cdn.elifesciences.org/articles/55526/elife-55526-fig3-data1-v1.xlsx
Figure 3—figure supplement 1
Nkx2.5-cre recombination, additional validation of NKX2-1-independent genes.

(a) Lineage trace of Nkx2.5-cre recombination pattern in E10.5 (top) and E11.5 (bottom) foreguts using the ROSA26mTmG reporter (red/green) and NKX2-1 immunofluorescent staining (magenta). (b) RNA …

NKX2-1 binds directly to tracheoesophageal genes in the developing foregut.

(a) ChIP-seq for NKX2-1 in E11.5 foreguts identified 15,861 NKX2-1-bound genomic regions (peaks) shared between replicates (FDR < 0.00001). n = 2 biological replicates with 175 pooled …

NKX2-1 regulates target genes in a SOX2-dependent and independent manner.

(a-c) Schematic of Nkx2-1; Sox2 compound mutant analysis phenotypes and resulting NKX2-1 and SOX2 expression patterns in E11.5 control embryos (a), Nkx2-1-TrKO mutants (b), Nkx2-1-TrKO; Sox2-TrKO …

Figure 6 with 1 supplement
NKX2-1 regulates mesenchymal specification via Wnt7b and Shh-mediated epithelial-mesenchymal crosstalk.

(a) Alcian blue staining of tracheal cartilage in E18.5 control and Nkx2-1-/- embryos. 3/3 embryos examined exhibited similar cartilage phenotype. (b) ChIP-seq of direct NKX2-1 binding near Wnt7b

Figure 6—figure supplement 1
NKX2-1 regulation of epithelial-mesenchymal signaling.

(a) RNA localization of Wnt7b (red) and Axin2 (yellow), and (b) Shh (red) and Ptch1 (yellow) in WT and Nkx2-1-/- E11.5 foreguts. (c) Schematic of rostral-caudal position of images in (d-g). (d-g) …

Proposed model of NKX2-1-dependent and -independent tracheoesophageal specification.

(a) Schematic of NKX2-1 and SOX2 regulation of gene expression in WT trachea and esophagus, and epithelial-mesenchymal signaling downstream of NKX2-1. NKX2-1 negatively regulates esophageal genes Shh…

Tables

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional
information
Strain, strain background (mouse)CD1Harlan/EnvigoCat#: 030
Genetic reagent (mouse)Nkx2-1fl/flKusakabe et al., 2006MGI: 3653645
Genetic reagent (mouse)Sox2fl/flShaham et al., 2009MGI: 4366453
Genetic reagent (mouse)Tmem163Tg(Actincre)Lewandoski et al., 1997MGI: 2176050
Genetic reagent (mouse)Nkx2.5creStanley et al., 2004MGI: 2448972
Antibodyα-NKX2-1 (rabbit polyconal)MilliporeCat#: 076011:200 for IF, 5ug for ChIP
Antibodyα-SOX2 (goat polyclonal)NeuromicsCat#: GT150981:250 for IF
Antibodyα-SOX9 (Rabbit polyclonal)Santa CruzCat#: sc-200951:250 for IF
Antibodyα-PITX1 (Rabbit polyclonal)NovusBioCat: NBP1886441:200 for IF
Antibodyα-LRIG1 (Goat polyclonal)R&D SystemsCat#: AF36881:250 for IF
Antibodyα-SMA (Rabbit polyclonal)AbcamAB56941:300 for IF
Antibodyα-SMA-cy3 (mouse monoclonal)SigmaC61981:300 for IF
Antibodyα-EpCAM (rat polyclonal)BioLegendCat#: 1182191:200 for FACS
Sequence-based reagentRNAscope probe mm-Irx2Advanced Cell DiagnosticsCat#: 519901
Sequence-based reagentRNAscope probe mm-Wnt7bAdvanced Cell DiagnosticsCat#: 401131
Sequence-based reagentRNAscope probe mm-Pcdh10Advanced Cell DiagnosticsCat#: 477781-C3
Sequence-based reagentRNAscope probe mm-Tppp3Advanced Cell DiagnosticsCat#: 586631
Sequence-based reagentRNAscope probe mm-Ly6hAdvanced Cell DiagnosticsCat#: 587811
Sequence-based reagentRNAscope probe mm-Krt19Advanced Cell DiagnosticsCat#: 402941
Sequence-based reagentRNAscope probe mm-Klf5Advanced Cell DiagnosticsCat#: 444081
Sequence-based reagentRNAscope probe mm-Foxe1Advanced Cell DiagnosticsCat#: 509641
Sequence-based reagentRNAscope probe mm-Crabp1Advanced Cell DiagnosticsCat#: 474711-C3
Sequence-based reagentRNAscope probe mm-Bmp4Advanced Cell DiagnosticsCat#: 401301-C2
Sequence-based reagentRNAscope probe mm-Nrp2Advanced Cell DiagnosticsCat#: 500661
Sequence-based reagentRNAscope probe mm-DcnAdvanced Cell DiagnosticsCat#: 413281-C3
Sequence-based reagentRNAscope probe mm-Has2Advanced Cell DiagnosticsCat#: 465171-C2
Sequence-based reagentRNAscope probe mm-Meis2Advanced Cell DiagnosticsCat#: 436371-C3
Sequence-based reagentRNAscope probe mm-Ackr3Advanced Cell DiagnosticsCat#: 482561-C2
Sequence-based reagentRNAscope probe mm-Axin2Advanced Cell DiagnosticsCat#: 400331-C3
Sequence-based reagentRNAscope probe mm-ShhAdvanced Cell DiagnosticsCat#: 314361
Sequence-based reagentRNAscope probe mm-Ptch1Advanced Cell DiagnosticsCat#: 402811-C2
Sequence-based reagentNkx2-1-geno1This paperPCR primer5'-CTA-GGG-AGG-CTA -GGA-ACT-CGG-3'
Sequence-based reagentNkx2-1-geno2This paperPCR primer5'-CCG-ACC-CAC-GTA-GAG-CC-3'
Sequence-based reagentNkx2-1-geno3This paperPCR primer5'-CTC-TTA-TCT-GGG-ATC-GCC-TGA-G-3'
Sequence-based reagentSox2-flox-geno1Steevens et al., 2019PCR primer5'-TGG-AAT-CAG-GCT-GCC-GAG-AAT-CC-3'
Sequence-based reagentSox2-flox-geno2Steevens et al., 2019PCR primer5'-TCG-TTG-TGG-CAA-CAA-GTG-CTA-AAG-C-3'
Sequence-based reagentSox2-flox-geno3Steevens et al., 2019PCR primer5'-CTG-CCA-TAG-CCA-CTC-GAG-AAG-3'
Sequence-based reagentCre-geno1Liang et al., 2005; Vauclair et al., 2005PCR primer5'-GTT-CGC-AAG-AAC-CTG-ATG-GAC-A-3'
Sequence-based reagentCre-geno2Liang et al., 2005; Vauclair et al., 2005PCR primer5'-CTA-GAG-CCT-GTT-TTG-CAC-GTT-C-3'
Commercial assay or kitMicroChIP Diapure columnsDiagenodeCat#: C03040001DNA extraction (ChIP)
Commercial assay or kitMicroplex Library Preparation Kit v2DiagenodeCat#: C05010012ChIP library prep
Commercial assay or kitRneasy Micro kitQiagenCat# 74004RNA extraction (bulk)
Commercial assay or kitSMARTer Stranded Total RNAseq kit V2TakaraCat#: 634411Bulk RNAseq library prep
Commercial assay or kitRNA 6000 Pico Bioanalyzer kitAgilentCat# 5067-1513Bioanalyzer (RNA)
Commercial assay or kitHigh Sensitivity DNA Bioanalyzer kitAgilentCat#: 5067-4626Bioanalyzer (DNA)
Commercial assay or kitQuBit dsDNA HS Assay kitInvitrogenCat#: Q32854DNA quantification
Commercial assay or kitChromium Single Cell ' Reagent Kit V210X GenomicsCat#: PN-120237, PN-120236, PN-120262Single cell RNA-seq
OtherSytox Blue nuclei acid stainThermoCat: S11348Live/dead stain for FACS, use at 1uM
OtherTSAplus FluoresceineAkoya BiosciencesNEL741001KTRNAscope, 1:1500
OtherTSAplus Cy3Akoya BiosciencesNEL744001KTRNAscope, 1:1500
OtherTSAplus Cy5Akoya BiosciencesNEL745001KTRNAscope, 1:1500
OtherDynaBeads, protein GInvitrogenCat#: 10003DChIP
OtherAMPure XP beadsBeckman CoulterCat#: A63881Library cleanup
OtherTrypLE ExpressThermoCat# 12604013Tissue dissociation

Additional files

Download links