XB-ART-54695
J Cell Biol
2018 Jun 04;2176:1941-1955. doi: 10.1083/jcb.201706106.
Show Gene links
Show Anatomy links
TRRAP is a central regulator of human multiciliated cell formation.
Wang Z
,
Plasschaert LW
,
Aryal S
,
Renaud NA
,
Yang Z
,
Choo-Wing R
,
Pessotti AD
,
Kirkpatrick ND
,
Cochran NR
,
Carbone W
,
Maher R
,
Lindeman A
,
Russ C
,
Reece-Hoyes J
,
McAllister G
,
Hoffman GR
,
Roma G
,
Jaffe AB
.
Abstract
The multiciliated cell (MCC) is an evolutionarily conserved cell type, which in vertebrates functions to promote directional fluid flow across epithelial tissues. In the conducting airway, MCCs are generated by basal stem/progenitor cells and act in concert with secretory cells to perform mucociliary clearance to expel pathogens from the lung. Studies in multiple systems, including Xenopus laevis epidermis, murine trachea, and zebrafish kidney, have uncovered a transcriptional network that regulates multiple steps of multiciliogenesis, ultimately leading to an MCC with hundreds of motile cilia extended from their apical surface, which beat in a coordinated fashion. Here, we used a pool-based short hairpin RNA screening approach and identified TRRAP, an essential component of multiple histone acetyltransferase complexes, as a central regulator of MCC formation. Using a combination of immunofluorescence, signaling pathway modulation, and genomic approaches, we show that (a) TRRAP acts downstream of the Notch2-mediated basal progenitor cell fate decision and upstream of Multicilin to control MCC differentiation; and (b) TRRAP binds to the promoters and regulates the expression of a network of genes involved in MCC differentiation and function, including several genes associated with human ciliopathies.
PubMed ID: 29588376
PMC ID: PMC5987713
Article link: J Cell Biol
Species referenced: Xenopus laevis
Genes referenced: ccno foxj1.2 itga6 mcc mcidas muc5b myb notch2 pcbd1 pcnt rfx3 tp63 trrap
Article Images: [+] show captions
![]() |
Fig 1 A pooled shRNA screen to identify regulators of basal cell fate. (A) Workflow of the pooled shRNA screen. Human airway basal cells were transduced with a lentiviral shRNA library 1 h after they were seeded on Transwell filters at day 0. At day 21, cells were fixed, stained, and sorted into three different groups based on expression of the ciliated cell marker FOXJ1 and the basal cell marker ITGA6. The extracted gDNA from the sorted samples was analyzed for shRNA hairpin counts by NGS. (B) Gene-centric visualization of the log2 fold change and RSA score in the ciliated (ITGA6− FOXJ1+) versus secretory (ITGA6− FOXJ1−) cells. (C) Ratio of individual shRNA hairpin counts in ciliated/secretory cells from the pooled shRNA screen. KDM8 is shown as an example of a gene whose shRNA barcode counts do not differ between ciliated and secretory cells. (D) In a primary validation, the cells from two independent donors were transduced with four shRNAs for each gene candidate and sorted with FOXJ1 and ITGA6 signals. The ratios of ciliated/secretory cells from all shRNA-treated samples were normalized to those of control cells transduced with shNT. Data were acquired from cells of two donors; mean ± SEM; ***, P < 0.001; Student’s two-tailed t test. |
![]() |
Fig 2 TRRAP is required for MCC formation. (A) Schematic workflow for validation of TRRAP. Six Transwell filters of fully differentiated ALI cultures derived from cells transduced with either shNT or TRRAP shRNA (shTRRAP) were pooled together. Two thirds were fixed and stained for ciliated (FOXJ1) and basal (ITGA6) markers and analyzed by flow cytometry. The remaining third was sorted with RFP to estimate gene knockdown by qPCR. (B) Representative FACS plots from TRRAP shRNA treatment. Human airway basal cells transduced with viruses for either shNT or TRRAP shRNAs were fixed and sorted with FOXJ1 and ITGA6 signals. (C) Ratios of ciliated/secretory cells in the TRRAP shRNA–treated samples were normalized to the ones from NT controls. (D) Quantification of TRRAP knockdown. CT values from qPCR in the TRRAP shRNA transduced cells were normalized to those from shNT controls (n = 6 from two independent donors; mean ± SEM; ***, P < 0.001). (E) The cells transduced with either shNT or shTRRAP were immunostained for markers of goblet cells (MUC5AC, MUC5B), club cells (SCGB1A1), ciliated cells (AC-TUBULIN, FOXJ1), or basal cells (p63). Representative images from ALI cultures immunostained for AC-TUBULIN/MUC5B (top) or AC-TUBULIN/SCGB1A1 (bottom) are shown. Bar, 10 µm. (F) The RFP+ cells expressing the indicated cell markers were normalized to the total RFP+ cells to calculate the percentage of goblet, club, ciliated, or basal cells after TRRAP knockdown or in shNT controls. n = 6 from two independent donors; mean ± SEM; **, P < 0.01; ***, P < 0.001; Student’s two-tailed t test. |
![]() |
Fig 3 TRRAP is enriched in ciliated cells and accumulates before FOXJ1. (A–D) During the course of differentiation at ALI, cultures were stained with antibodies to TRRAP, FOXJ1 (ciliated cells), MUC5B (goblet cells), SCGB1A1 (club cells), and p63 (basal cells). Colocalization was quantified at days 0, 3, 7, and 14 at ALI (A′′′–D′′′). At ALI day 14, TRRAP colocalizes with FOXJ1 in ciliated cells (arrows in A–A′′) but rarely accumulates in other cell types (A′–D′′). (E–G) At day 0, before culture at ALI, TRRAP has begun to accumulate in nuclei with a small subset of those being coenriched for FOXJ1 (arrows; A′′′; arrow in E–E′′). At days 3 and 7 of culture at ALI, the number of cells enriched for both TRRAP and FOXJ1 increases (A′′′; arrows in F–G′′), and the number of TRRAP-only cells decreases. Notably, cells with high TRRAP accumulation have low FOXJ1 (F–G′′, arrowheads), whereas cells with high FOXJ1 accumulation have low TRRAP (F–G′′, arrows), suggesting that their expression is temporally distributed in the same cell type. Bars, 20 µm. |
![]() |
Fig 4 TRRAP acts downstream of Notch2 signaling to regulate MCC formation. (A) Representative FACS plots of shRNA-infected ALI cultures that were treated with IgG or α-N2. Human airway basal cells were transduced with either shNT or shTRRAP-a at day 0. The cells were cultured in medium containing either 1 µg/ml IgG or α-N2 from day 7 (ALI day 0). At day 21, cells were stained for FOXJ1 and ITGA6 to label ciliated and basal cells, respectively, and analyzed by flow cytometry. (B) The ratio of ciliated/secretory cells of each sample was normalized to the control cells treated with IgG and shNT (n = 6 from two independent donors; mean ± SEM). (C) Representative images of ALI cultures with the indicated treatments, stained for AC-TUBULIN, to label ciliated cells, and either MUC5B (top row) or SCGB1A1 (bottom row) to label goblet or club cells, respectively. Bar, 10 µm. (D) Quantification of the percentage of RFP+ (shRNA-containing) ciliated (AC-TUBULIN+), goblet (MUC5B+), and club (SCGB1A1+) cells. Data were obtained from six independent cultures with two independent donors. Error bars indicate SEM. *, P < 0.05; **, P < 0.01; ***, P < 0.001; Student’s two-tailed t test. |
![]() |
Fig 5 TRRAP regulates MCC formation upstream of MULTICILIN. (A) Human airway basal cells transduced with either shNT or shTRRAP-b at day 0 were differentiated at ALI and stained for RFP (red), PERICENTRIN (green), or MUC5B (white). Arrowheads highlight RFP+MUC5B− presumptive ciliated cells. Note that knockdown of TRRAP results in an inhibition of centriole replication. Bar, 10 µm. (B) Quantification of the percentage of shRNA-expressing presumptive ciliated cells (RFP+ MUC5B−) with more than two centrioles. Each condition was performed in two independent donors, and the total number of cells quantified is indicated. (C) Representative FACS plots of day 21 ALI cultures, which had been transduced with lentivirus encoding for GFP or MULTICILIN-GFP, with either shNT or shTRRAP-b at day 0. Cultures were dissociated and then sorted based on FOXJ1 and ITGA6 signals. (D) Quantification of the ratio of ciliated/secretory cells from the indicated treatments, normalized to control cells that were cotransduced with shNT and GFP (n = 3). Data shown are represented as mean ± SEM; *, P < 0.05; **, P < 0.01; Student’s two-tailed t test. |
![]() |
Fig 6 TRRAP regulates a network of genes controlling multiple steps of multiciliogenesis. (A) Schematic of the experimental design. (B) Venn diagram of the number of TRRAP DEGs, identified by RNA-seq, and the number of promoter regions (<1 kb from the transcription start site) containing TRRAP peaks. (C) Enrichment analysis of TRRAP DTGs in the indicated datasets. The number of DTGs within each dataset is indicated. Shown are the p-values calculated by Fisher’s exact test with Bonferroni correction. Note the TRRAP DTG enrichment in sets of MCC genes (Hoh et al., 2012; van Dam et al., 2013; Nemajerova et al., 2016; Quigley and Kintner, 2017) and genes that are mutated in PCD; (Horani and Ferkol, 2016) but not in the set of basal cell genes (Hackett et al., 2011). (D) Hierarchical clustering of the 271 ciliogenesis DEGs at ALI days 4 and 14 in shNT, shTRRAP-a, and shTRRAP-b samples. Shown are the z-scores for each gene across all samples. Right: Bar plot displaying the maximum log2 expression differences between shTRRAP and shNT samples. Gray, DEGs; blue, DEGs containing a TRRAP peak in the promoter. (E) ChIP-seq tracks displaying TRRAP binding at the MCIDAS, CCNO, MYB, and RFX3 genomic loci. Shown is the fold enrichment over input. |
![]() |
Figure 1. A pooled shRNA screen to identify regulators of basal cell fate. (A) Workflow of the pooled shRNA screen. Human airway basal cells were transduced with a lentiviral shRNA library 1 h after they were seeded on Transwell filters at day 0. At day 21, cells were fixed, stained, and sorted into three different groups based on expression of the ciliated cell marker FOXJ1 and the basal cell marker ITGA6. The extracted gDNA from the sorted samples was analyzed for shRNA hairpin counts by NGS. (B) Gene-centric visualization of the log2 fold change and RSA score in the ciliated (ITGA6− FOXJ1+) versus secretory (ITGA6− FOXJ1−) cells. (C) Ratio of individual shRNA hairpin counts in ciliated/secretory cells from the pooled shRNA screen. KDM8 is shown as an example of a gene whose shRNA barcode counts do not differ between ciliated and secretory cells. (D) In a primary validation, the cells from two independent donors were transduced with four shRNAs for each gene candidate and sorted with FOXJ1 and ITGA6 signals. The ratios of ciliated/secretory cells from all shRNA-treated samples were normalized to those of control cells transduced with shNT. Data were acquired from cells of two donors; mean ± SEM; ***, P < 0.001; Student’s two-tailed t test. |
![]() |
Figure 2. TRRAP is required for MCC formation. (A) Schematic workflow for validation of TRRAP. Six Transwell filters of fully differentiated ALI cultures derived from cells transduced with either shNT or TRRAP shRNA (shTRRAP) were pooled together. Two thirds were fixed and stained for ciliated (FOXJ1) and basal (ITGA6) markers and analyzed by flow cytometry. The remaining third was sorted with RFP to estimate gene knockdown by qPCR. (B) Representative FACS plots from TRRAP shRNA treatment. Human airway basal cells transduced with viruses for either shNT or TRRAP shRNAs were fixed and sorted with FOXJ1 and ITGA6 signals. (C) Ratios of ciliated/secretory cells in the TRRAP shRNA–treated samples were normalized to the ones from NT controls. (D) Quantification of TRRAP knockdown. CT values from qPCR in the TRRAP shRNA transduced cells were normalized to those from shNT controls (n = 6 from two independent donors; mean ± SEM; ***, P < 0.001). (E) The cells transduced with either shNT or shTRRAP were immunostained for markers of goblet cells (MUC5AC, MUC5B), club cells (SCGB1A1), ciliated cells (AC-TUBULIN, FOXJ1), or basal cells (p63). Representative images from ALI cultures immunostained for AC-TUBULIN/MUC5B (top) or AC-TUBULIN/SCGB1A1 (bottom) are shown. Bar, 10 µm. (F) The RFP+ cells expressing the indicated cell markers were normalized to the total RFP+ cells to calculate the percentage of goblet, club, ciliated, or basal cells after TRRAP knockdown or in shNT controls. n = 6 from two independent donors; mean ± SEM; **, P < 0.01; ***, P < 0.001; Student’s two-tailed t test. |
![]() |
Figure 3. TRRAP is enriched in ciliated cells and accumulates before FOXJ1. (A–D) During the course of differentiation at ALI, cultures were stained with antibodies to TRRAP, FOXJ1 (ciliated cells), MUC5B (goblet cells), SCGB1A1 (club cells), and p63 (basal cells). Colocalization was quantified at days 0, 3, 7, and 14 at ALI (A′′′–D′′′). At ALI day 14, TRRAP colocalizes with FOXJ1 in ciliated cells (arrows in A–A′′) but rarely accumulates in other cell types (A′–D′′). (E–G) At day 0, before culture at ALI, TRRAP has begun to accumulate in nuclei with a small subset of those being coenriched for FOXJ1 (arrows; A′′′; arrow in E–E′′). At days 3 and 7 of culture at ALI, the number of cells enriched for both TRRAP and FOXJ1 increases (A′′′; arrows in F–G′′), and the number of TRRAP-only cells decreases. Notably, cells with high TRRAP accumulation have low FOXJ1 (F–G′′, arrowheads), whereas cells with high FOXJ1 accumulation have low TRRAP (F–G′′, arrows), suggesting that their expression is temporally distributed in the same cell type. Bars, 20 µm. |
![]() |
Figure 4. TRRAP acts downstream of Notch2 signaling to regulate MCC formation. (A) Representative FACS plots of shRNA-infected ALI cultures that were treated with IgG or α-N2. Human airway basal cells were transduced with either shNT or shTRRAP-a at day 0. The cells were cultured in medium containing either 1 µg/ml IgG or α-N2 from day 7 (ALI day 0). At day 21, cells were stained for FOXJ1 and ITGA6 to label ciliated and basal cells, respectively, and analyzed by flow cytometry. (B) The ratio of ciliated/secretory cells of each sample was normalized to the control cells treated with IgG and shNT (n = 6 from two independent donors; mean ± SEM). (C) Representative images of ALI cultures with the indicated treatments, stained for AC-TUBULIN, to label ciliated cells, and either MUC5B (top row) or SCGB1A1 (bottom row) to label goblet or club cells, respectively. Bar, 10 µm. (D) Quantification of the percentage of RFP+ (shRNA-containing) ciliated (AC-TUBULIN+), goblet (MUC5B+), and club (SCGB1A1+) cells. Data were obtained from six independent cultures with two independent donors. Error bars indicate SEM. *, P < 0.05; **, P < 0.01; ***, P < 0.001; Student’s two-tailed t test. |
![]() |
Figure 5. TRRAP regulates MCC formation upstream of MULTICILIN. (A) Human airway basal cells transduced with either shNT or shTRRAP-b at day 0 were differentiated at ALI and stained for RFP (red), PERICENTRIN (green), or MUC5B (white). Arrowheads highlight RFP+MUC5B− presumptive ciliated cells. Note that knockdown of TRRAP results in an inhibition of centriole replication. Bar, 10 µm. (B) Quantification of the percentage of shRNA-expressing presumptive ciliated cells (RFP+ MUC5B−) with more than two centrioles. Each condition was performed in two independent donors, and the total number of cells quantified is indicated. (C) Representative FACS plots of day 21 ALI cultures, which had been transduced with lentivirus encoding for GFP or MULTICILIN-GFP, with either shNT or shTRRAP-b at day 0. Cultures were dissociated and then sorted based on FOXJ1 and ITGA6 signals. (D) Quantification of the ratio of ciliated/secretory cells from the indicated treatments, normalized to control cells that were cotransduced with shNT and GFP (n = 3). Data shown are represented as mean ± SEM; *, P < 0.05; **, P < 0.01; Student’s two-tailed t test. |
![]() |
Figure 6. TRRAP regulates a network of genes controlling multiple steps of multiciliogenesis. (A) Schematic of the experimental design. (B) Venn diagram of the number of TRRAP DEGs, identified by RNA-seq, and the number of promoter regions (<1 kb from the transcription start site) containing TRRAP peaks. (C) Enrichment analysis of TRRAP DTGs in the indicated datasets. The number of DTGs within each dataset is indicated. Shown are the p-values calculated by Fisher’s exact test with Bonferroni correction. Note the TRRAP DTG enrichment in sets of MCC genes (Hoh et al., 2012; van Dam et al., 2013; Nemajerova et al., 2016; Quigley and Kintner, 2017) and genes that are mutated in PCD; (Horani and Ferkol, 2016) but not in the set of basal cell genes (Hackett et al., 2011). (D) Hierarchical clustering of the 271 ciliogenesis DEGs at ALI days 4 and 14 in shNT, shTRRAP-a, and shTRRAP-b samples. Shown are the z-scores for each gene across all samples. Right: Bar plot displaying the maximum log2 expression differences between shTRRAP and shNT samples. Gray, DEGs; blue, DEGs containing a TRRAP peak in the promoter. (E) ChIP-seq tracks displaying TRRAP binding at the MCIDAS, CCNO, MYB, and RFX3 genomic loci. Shown is the fold enrichment over input. |
References [+] :
Boon,
MCIDAS mutations result in a mucociliary clearance disorder with reduced generation of multiple motile cilia.
2016, Pubmed,
Xenbase
Boon, MCIDAS mutations result in a mucociliary clearance disorder with reduced generation of multiple motile cilia. 2016, Pubmed , Xenbase
Brooks, Multiciliated cells. 2015, Pubmed
Bustamante-Marin, Cilia and Mucociliary Clearance. 2017, Pubmed
Cunningham, Ensembl 2015. 2015, Pubmed
Danahay, Notch2 is required for inflammatory cytokine-driven goblet cell metaplasia in the lung. 2015, Pubmed
Danahay, Interleukin-13 induces a hypersecretory ion transport phenotype in human bronchial epithelial cells. 2002, Pubmed
Danielian, E2f4 is required for normal development of the airway epithelium. 2007, Pubmed
El Zein, RFX3 governs growth and beating efficiency of motile cilia in mouse and controls the expression of genes involved in human ciliopathies. 2009, Pubmed
Fahy, Airway mucus function and dysfunction. 2010, Pubmed
Fortini, Notch signaling: the core pathway and its posttranslational regulation. 2009, Pubmed
Funk, Cyclin O (Ccno) functions during deuterosome-mediated centriole amplification of multiciliated cells. 2015, Pubmed
Guseh, Notch signaling promotes airway mucous metaplasia and inhibits alveolar development. 2009, Pubmed
Hackett, The human airway epithelial basal cell transcriptome. 2011, Pubmed
Herceg, Disruption of Trrap causes early embryonic lethality and defects in cell cycle progression. 2001, Pubmed
Hoffman, Functional epigenetics approach identifies BRM/SMARCA2 as a critical synthetic lethal target in BRG1-deficient cancers. 2014, Pubmed
Hogan, Repair and regeneration of the respiratory system: complexity, plasticity, and mechanisms of lung stem cell function. 2015, Pubmed
Hoh, Transcriptional program of ciliated epithelial cells reveals new cilium and centrosome components and links to human disease. 2013, Pubmed
Horani, Primary ciliary dyskinesia and associated sensory ciliopathies. 2017, Pubmed
Hrvatin, MARIS: method for analyzing RNA following intracellular sorting. 2015, Pubmed
Humbert, E2F4 is essential for normal erythrocyte maturation and neonatal viability. 2000, Pubmed
Kamentsky, Improved structure, function and compatibility for CellProfiler: modular high-throughput image analysis software. 2011, Pubmed
Kharchenko, Design and analysis of ChIP-seq experiments for DNA-binding proteins. 2009, Pubmed
König, A probability-based approach for the analysis of large-scale RNAi screens. 2007, Pubmed
Lafkas, Therapeutic antibodies reveal Notch control of transdifferentiation in the adult lung. 2016, Pubmed
Lang, E2F transcriptional activation requires TRRAP and GCN5 cofactors. 2001, Pubmed
Leachman, ATAD2B is a phylogenetically conserved nuclear protein expressed during neuronal differentiation and tumorigenesis. 2011, Pubmed
Li, The Sequence Alignment/Map format and SAMtools. 2010, Pubmed
Li, Fast and accurate short read alignment with Burrows-Wheeler transform. 2009, Pubmed
Loizou, Histone acetyltransferase cofactor Trrap is essential for maintaining the hematopoietic stem/progenitor cell pool. 2009, Pubmed
Lunyak, Epigenetic regulation of stem cell fate. 2008, Pubmed
Ma, Multicilin drives centriole biogenesis via E2f proteins. 2014, Pubmed , Xenbase
McMahon, The novel ATM-related protein TRRAP is an essential cofactor for the c-Myc and E2F oncoproteins. 1998, Pubmed
Morimoto, Different assemblies of Notch receptors coordinate the distribution of the major bronchial Clara, ciliated and neuroendocrine cells. 2013, Pubmed
MuhChyi, Epigenetic regulation of neural stem cell fate during corticogenesis. 2014, Pubmed
Murr, Orchestration of chromatin-based processes: mind the TRRAP. 2007, Pubmed
Nemajerova, TAp73 is a central transcriptional regulator of airway multiciliogenesis. 2017, Pubmed
Quigley, Rfx2 Stabilizes Foxj1 Binding at Chromatin Loops to Enable Multiciliated Cell Gene Expression. 2017, Pubmed , Xenbase
Revenko, Chromatin loading of E2F-MLL complex by cancer-associated coregulator ANCCA via reading a specific histone mark. 2010, Pubmed
Ritchie, limma powers differential expression analyses for RNA-sequencing and microarray studies. 2015, Pubmed
Rock, Basal cells as stem cells of the mouse trachea and human airway epithelium. 2009, Pubmed
Roostaee, Epigenetics in Intestinal Epithelial Cell Renewal. 2017, Pubmed
Sawan, Histone acetyltransferase cofactor Trrap maintains self-renewal and restricts differentiation of embryonic stem cells. 2014, Pubmed
Schuierer, The exon quantification pipeline (EQP): a comprehensive approach to the quantification of gene, exon and junction expression from RNA-seq data. 2017, Pubmed
Snitow, Ezh2 represses the basal cell lineage during lung endoderm development. 2015, Pubmed
Stubbs, Multicilin promotes centriole assembly and ciliogenesis during multiciliate cell differentiation. 2012, Pubmed , Xenbase
Tan, Myb promotes centriole amplification and later steps of the multiciliogenesis program. 2013, Pubmed , Xenbase
Tapias, Trrap-dependent histone acetylation specifically regulates cell-cycle gene transcription to control neural progenitor fate decisions. 2014, Pubmed
Tsao, Notch signaling controls the balance of ciliated and secretory cell fates in developing airways. 2009, Pubmed
Wallmeier, Mutations in CCNO result in congenital mucociliary clearance disorder with reduced generation of multiple motile cilia. 2014, Pubmed , Xenbase
Wang, Development and regeneration of Sox2+ endoderm progenitors are regulated by a Hdac1/2-Bmp4/Rb1 regulatory pathway. 2013, Pubmed
Wurdak, An RNAi screen identifies TRRAP as a regulator of brain tumor-initiating cell differentiation. 2010, Pubmed
Yu, ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. 2016, Pubmed
Zhang, Model-based analysis of ChIP-Seq (MACS). 2009, Pubmed
van Dam, The SYSCILIA gold standard (SCGSv1) of known ciliary components and its applications within a systems biology consortium. 2013, Pubmed