Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Transcriptional Program of Ciliated Epithelial Cells Reveals New Cilium and Centrosome Components and Links to Human Disease

  • Ramona A. Hoh,

    Affiliation Department of Biology, Stanford University, Stanford, California, United States of America

    ⨯
  • Timothy R. Stowe,

    Affiliation Department of Biology, Stanford University, Stanford, California, United States of America

    ⨯
  • Erin Turk,

    Affiliation Department of Biology, Stanford University, Stanford, California, United States of America

    ⨯
  • Tim Stearns

    [email protected]

    Affiliations Department of Biology, Stanford University, Stanford, California, United States of America, Department of Genetics, Stanford University School of Medicine, Stanford, California, United States of America

    ⨯

Abstract

Defects in the centrosome and cilium are associated with a set of human diseases having diverse phenotypes. To further characterize the components that define the function of these organelles we determined the transcriptional profile of multiciliated tracheal epithelial cells. Cultures of mouse tracheal epithelial cells undergoing differentiation in vitro were derived from mice expressing GFP from the ciliated-cell specific FOXJ1 promoter (FOXJ1:GFP). The transcriptional profile of ciliating GFP+ cells from these cultures was defined at an early and a late time point during differentiation and was refined by subtraction of the profile of the non-ciliated GFP- cells. We identified 649 genes upregulated early, when most cells were forming basal bodies, and 73 genes genes upregulated late, when most cells were fully ciliated. Most, but not all, of known centrosome proteins are transcriptionally upregulated early, particularly Plk4, a master regulator of centriole formation. We found that three genes associated with human disease states, Mdm1, Mlf1, and Dyx1c1, are upregulated during ciliogenesis and localize to centrioles and cilia. This transcriptome for mammalian multiciliated epithelial cells identifies new candidate centrosome and cilia proteins, highlights similarities between components of motile and primary cilia, and identifies new links between cilia proteins and human disease.

Citation: Hoh RA, Stowe TR, Turk E, Stearns T (2012) Transcriptional Program of Ciliated Epithelial Cells Reveals New Cilium and Centrosome Components and Links to Human Disease. PLoS ONE 7(12): e52166. https://doi.org/10.1371/journal.pone.0052166

Editor: Panayiotis V. Benos, University of Pittsburgh, United States of America

Received: November 1, 2012; Accepted: November 15, 2012; Published: December 31, 2012

Copyright: © 2012 Hoh et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported by an NSF graduate research fellowship to T.R.S., the Gabilan Stanford Graduate Fellowship to E.T., and United States National Institutes of Health grant GM52022 to T.S. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Cilia are highly conserved microtubule-based structures that are found in all branches of the eukaryotic tree [1]. The ciliary axoneme is nucleated and anchored at the cell membrane by a modified centriole called the basal body, and is surrounded by a specialized membrane that is continuous with the plasma membrane. There are several types of cilia, differing in the specifics of their structure and function. The primary cilium, present in one copy on many cell types in the body, has a 9+0 axonemal structure and is important for sensing chemical and mechanical stimuli. Important developmental signaling pathways such as hedgehog and non-canonical Wnt are transduced through primary cilia [2], [3]. Motile cilia usually have an additional central pair of microtubules in the axoneme (9+2), as well as accessory structures including dynein arms and radial spokes associated with ciliary beating [4]. In addition to generating force, motile cilia also have sensory functions [5], [6], [7], [8].

Centrioles are a conserved feature of all organisms that have ciliated cells, and one of the main functions of centrioles is to serve as basal bodies [9], [10], [11], [12], [13]. In most cycling cells, centrioles lie within the larger centrosome, an assemblage of many proteins that functions as a microtubule organizer and a signaling center. Centrioles duplicate once per cell cycle and segregate at mitosis, maintaining the capacity for cilium formation in most cells. Although a number of human diseases have now been linked to defects in centriolar or centrosomal proteins, cell division can often proceed normally in the absence of centrioles, and adult flies generated without centrioles are viable but have defects associated with an inability to make cilia [9], [14]. We consider that it is most useful to consider the centrioles, centrosome and cilium as a single functional complex in animal cells, and will refer to it as the centrosome/cilium when appropriate.

Genetic defects in the centrosome/cilium are associated with a range of phenotypes in humans, consistent with their broad set of functions. The most clear-cut defects are those associated with ciliary motility. Motile cilia are found both in single copy, in sperm and in specialized epithelial cells in the embryonic node, and in hundreds of copies, on the surface of ciliated epithelial cells lining the respiratory airways, the ventricles of the brain, and the oviduct. Defects in motile cilia result in failure of embryonic turning, respiratory failure, infertility, hydrocephalus and randomized left-right asymmetry [15]. More recently it has been found that there is a set of human disease phenotypes that are caused by defects in the centrosome/cilium, but which are distinct from those strictly related to motile cilia [16]. These phenotypes include retinal degeneration, obesity, sterility, polydactyly, polycystic kidney disease and mental retardation. The diversity of pathologies reflect the important roles that cilia play in a multitude of tissues in the body [15], and in many cases the nature of the connection between centrosome/cilium and phenotype is not understood. Collectively, diseases caused by an underlying defect in cilia are termed ciliopathies.

Previously, we examined the process of basal body and cilium formation in multiciliated tracheal epithelial cells (MTECs) [17], [18]. Multiciliated cells are unique in that they produce 200–300 cilia and centrioles during differentiation, whereas most G1 phase cycling cells have two centrioles, and either lack a cilium or have a single primary cilium. Most of the centrioles in multiciliated cells are generated by a little-studied process in which multiple centrioles grow simultaneously from the surface of deuterosomes, structures of unknown origin [19], [20]. This duplication process is, at least outwardly, different from that in cycling cells, in which single new centrioles grow from the sides of the two existing centrioles [13]. Despite the differences in function and manner of formation, basal bodies and cycling cell centrioles share many of the same components [17]. Moreover, depletion of SAS-6 or CEP120, which are required for centriole formation in human and mouse cells [18], [21], blocked basal body generation in multiciliated cells [17], [18].

To further understand centrosome/cilium structure and function, we here exploit the multiciliated epithelial cell system by determining the transcriptional profile of ciliating cells during differentiation. This approach combines the in vitro culture model that we previously described [17] with a transgenic mouse that expresses GFP from the ciliated cell-specific promoter of the human forkhead-box transcription factor FOXJ1 [22]. The mouse tracheal epithelium is a pseudostratified epithelium that consists of a complex mixture of cells including basal cells, goblet cells and multiciliated cells. In culture, ciliogenesis is initiated by the establishment of an air-liquid interface (ALI) [23] and the process proceeds in discrete stages including the formation of basal body precursors in the cytoplasm, the migration of nascent basal bodies to the apical surface of the cells, docking at the plasma membrane and axoneme extension [17]. Importantly, ciliogenesis in culture is semi-synchronous, such that most FOXJ1-expressing cells at four days after ALI (ALI+4) are in the process of basal body formation, whereas by ALI+12, most are fully ciliated [17].

We find that the set of upregulated genes in multiciliated epithelial cells reveals the similarities between primary and motile cilia, and suggests how this cell type is uniquely able to generate hundreds of centrioles and cilia. In addition, we identify new components of these structures, including previously uncharacterized proteins and proteins associated with human disease. These results establish new links between the centrosome, cilium and genetic diseases with poorly understood molecular etiology.

Materials and Methods

Animals and Animal Care

MTECs were derived from wild-type C3H X C57Bl/6J F1 hybrid or FOXJ1:GFP transgenic mice (a gift from L. Ostrowski, University of North Carolina at Chapel Hill, Chapel Hill, NC), which were generated on C3H X C57B1/6J F1 hybrid background [22]. Heterozygous FOXJ1:GFP mice were obtained as described [17]. PCR genotyping was performed using a forward primer specific to the FOXJ1 promoter region (5′-GCAGGCACCACATACTTATTCGGAGG-3′) and a reverse primer specific to GFP (5′-CGTCCTTGAAGAAGATGGTGCG-3′). Mice were sacrificed by CO2 anesthesia and tracheas were surgically removed. All work was approved by the Stanford University Administrative Panel for Laboratory Animal Care (SUAPLAC Protocol-11659) and was performed in accordance with SUAPLAC guidelines.

MTEC Culture and Sorting

MTECs for the microarray experiments were derived from tracheas removed from male and female mice between 6–20 weeks of age from multiple litters. For microarray experiments, MTECs were pooled and cultured on 6-well plates containing Transwell-Clear permeable membrane supports (Corning) for the microarray experiments. Cells dissociated from the tracheas of approximately 12–15 mice were seeded per 6-well plate. Culture and differentiation of MTECs was performed as previously described [17]. For each biological replicate, fourteen to sixteen wells were pooled and sorted for ALI+4; four to six wells were pooled for ALI+12. Technical replicates were performed using RNA from the same sample on a separate array. To prepare MTECs for FACS, cells were first removed from the filter by incubating at 37°C for 30 minutes in a 1∶1 mixture of Cell Dissociation Solution (Sigma-Aldrich) and 0.5% Trypsin/EDTA (Invitrogen). After washing 2X in PBS, cells were resuspended in PBS at 106 cells/ml and passed through a 100 µm nylon mesh cell strainer (BD Falcon) to remove clumped cells. GFP+ and GFP- (from ALI+12) populations were collected in serum. GFP+ MTECs derived from FOXJ1:GFP mice had at least 10-fold higher FITC intensity than the mean FITC value of MTECs derived from wild-type FOXJ1−/−mice. Sorting was performed at the Stanford Shared FACS Facility on a Vantage Vanford sorter. Data were acquired using CellQuest software.

RNA Extraction

Total RNA was extracted using the RNeasy Kit (Qiagen) according to manufacturer’s protocols. Briefly, FACS sorted samples in serum were diluted 5∶1 in PBS and centrifuged for 5 minutes at 1200×g. The pellet was resuspended in 250 µl of Buffer RT and spun through a QiaShredder Column (Qiagen) to disrupt cells before proceeding with the rest of the protocol. RNA was eluted in 30 µl of RNase free water. Total RNA concentration was measured by NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, DE) and stored at –80°C until amplification.

Probe Preparation

Amplification and coupling of RNA to Cy3/Cy5 was performed with amino allyl MessageAmp II aRNA amplification kit (Ambion) according to the manufacturer’s protocols. Briefly, total RNA from MTECs or Universal Mouse Reference RNA (Stratagene) were reverse transcribed using the supplied T7 oligo-dT primer to produce first-strand cDNA. After second strand synthesis, cDNA was subjected to a single round of in vitro amplification with amino allyl-modified UTP. The modified RNA was purified, dried and coupled to either Cy3 or Cy5 reactive dyes. Dye incorporation and final probe concentration was checked by NanoDrop-1000 and 2–10 µg labeled RNA was fractionated using RNA Fragmentation Reagent (Ambion) before proceeding to hybridization.

Microarray Hybridization and Image Extraction

Hybridization was done on MEEBO oligonucleotide arrays printed at the Stanford Functional Genomics Facility (www.microarray.org). MEEBO slides were post-processed the same day of hybridization. Fluorescently labeled amplified RNA samples from MTECs were hybridized to the array with an amplified reference RNA sample (Stratagene) for 16 hours at 65°C. Arrays were scanned on an Agilent G2565AA scanner immediately after washing. Features were extracted using GenePix6.0 (Axon Instruments). Detailed microarray protocols are available at http://brownlab.stanford.edu/protocols.html.

MIAME

Data were submitted to Gene Expression Omnibus (GEO) according to MIAME (Minimum Information About a Microarray Experiment) guidelines and can be accessed at GSE42500.

Data Analysis

Linear Models for MicroArray (LIMMA) LIMMA software [24], an R-based program the employs Empirical Bayes for the analysis of microarray data, was used to identify differentially expressed genes and for data normalization. LIMMA software is available from the BioConductor website http://www.bioconductor.org/packages/release/bioc/html/limma.html. Genes were clustered by Pearson Correlation using Cluster 3.0 and visualized using TreeView as implemented by SMD.

Functional annotation of differentially expressed genes was performed using bioinformatics tools found online at DAVID (Database for Annotation, Visualization and Integrated Discovery, version 6.7, from the National Institute of Allergy and Infectious Diseases (NIAID), NIH (http://david.abcc.ncifcrf.gov/)) [25]. Genes that passed a fold-change cutoff of 2-fold or greater in addition to an adjusted P value <0.05 in the subtracted MTEC transcriptome from ALI+4 and ALI+12, or in the GFP- non-ciliated cells, were uploaded to DAVID, and functional groups that are overrepresented in these sets of genes relative to the human genome were identified.

Plasmids

cDNA encoding TTC12 (Accession BC067297) was obtained from Open Biosystems. cDNA encoding MDM1 was PCR amplified from a cDNA library generated from MTEC cultures harvested at ALI+10. The MTEC cDNA library was generated using Superscript II RT (Invitrogen) according to the manufacturer’s protocol. cDNAs were cloned in frame to pEGFP-N1 (Invitrogen) to generate C-terminal GFP fusion expression constructs, or to the lentiviral transfer vector pRRL.sin-18.PPT.PGK.IRES.GFP.pre [17] using either the Age I site or the Age I site in combination with Bam HI. All PCR generated clones were sequenced to ensure in-frame fusion to GFP. The rat DYX1C1-GFP fusion in the pCAGGS plasmid was a gift from Dr. J.J. LoTurco. The cDNA library was generated using Superscript II RT (Invitrogen) according to the manufacturer’s protocols.

Cell Culture and Transfection

Mouse embryonic fibroblast NIH/3T3 cells (ATCC Catalog CRL-1658) and hTERT-immortalized human retinal pigment epithelium hTERT-RPE1 (ATCC Catalog CRL-4000) cells were cultured in DMEM with 10% FBS (Invitrogen). For transfection, cells were plated onto poly-L-lysine coated coverslips at 80% confluency the night before transfection. The day of transfection, cells were rinsed 2× with PBS and incubated with media containing plasmid plus Lipofectamine 2000 (Invitrogen) for 6 h before replacing with fresh media. Cells were imaged 24–48 h later.

Indirect Immunofluorescence

NIH/3T3 cells and MTECs were prepared for indirect IF by washing twice with PBS and fixing in 4% PFA at room temperature for 10 min. or methanol at −20°C for 10 min. After fixation, cells were rinsed with PBS and blocked in PBS-BT (PBS +3% BSA +0.1% Triton X-100) for 30 min. at room temperature. MTEC filters were cut at this point from the solid supports for staining. MTECs were incubated with primary antibody for 1 h at room temperature, while tissue culture cells were incubated for 30 min. at room temperature. Incubations with secondary antibodies were for 30 min. at room temperature for all cells. To visualize DNA, cells were incubated with 4′,6-diamidino-2-phenylindole (DAPI, Molecular Probes) for 5 min. at room temperature using a 1 µg/ml dilution for NIH/3T3 cells and a 10 µg/ml dilution for MTECs. In between incubations, cells were washed 3× with PBS-BT for 5 minutes each. The primary antibodies used were mouse anti-acetylated alpha-tubulin (1∶2000, Sigma), mouse anti-gamma-tubulin (1∶500, Sigma), mouse anti-MLF1 (1∶400, Abnova) and rabbit anti-polaris (1∶500, gift from Bradley Yoder). The secondary antibodies used were Alexa 488 and 594-conjugated goat anti-mouse and goat anti-rabbit (1∶200, Invitrogen). All antibodies were diluted in PBS-BT. Cells were then observed under a Zeiss Axiovert 200M microscope (Zeiss, Thornwood, NJ) using a 100× oil-immersion objective.

Results

Transcriptional Profile of Ciliogenesis in Tracheal Epithelial Cells

We used microarrays to identify genes that are differentially expressed in mouse tracheal epithelial cells (MTECs) differentiated in vitro. The cells were derived from a mouse strain expressing GFP under control of the FOXJ1 ciliated cell-specific promoter [22]. This promoter is activated early in the differentiation process [17] (Figure 1A), allowing the use of fluorescence activated cell sorting (FACS) to sort dissociated cells from the MTEC cultures into ciliating and non-ciliating populations (Figure 1B). MTECs were harvested at two time points after transition to an air-liquid interface culture (ALI), which induces differentiation of ciliated cells in the columnar epithelium [26]. The two time points chosen were based on previous analysis of landmark events in these cultures [17]; four days after ALI (ALI+4) to enrich for genes involved in initial steps of centriole duplication, and ALI+12 to enrich for genes expressed when cilia are mature (Figure 1A). At ALI+12 the majority of cells expressing GFP are fully ciliated, with a small fraction of cells at an earlier stage of ciliogenesis [17]. Importantly, the in vitro-differentiated epithelium accurately recapitulates the complexity of the in vivo tissue with respect to gene expression [27], with several cell types in addition to the ciliated epithelial cells of interest to this study. To separate ciliated from non-ciliated cells in our analysis, single-cell suspensions of MTECs were sorted into GFP+ and GFP- populations (Figure 1B) and enrichment for ciliated cells in the GFP+ population was demonstrated by immunofluorescence labeling of cilia (Figure 1C). RNA from GFP+ and GFP- samples was then hybridized to mouse predicted-exon oligo microarrays and compared against a universal mouse reference RNA pool. This yielded the transcriptional profile of GFP+ MTECs (ciliating or ciliated) at ALI+4 and ALI+12 and of non-ciliated, GFP- MTECs.

thumbnail
Download:
Figure 1. Microarray analysis of ciliating mouse tracheal epithelial cells. (

A) Immunofluorescence staining showing representative images of GFP expression and ciliogenesis in FOXJ1:GFP MTEC cultures before establishment of the air-liquid interface (pre-ALI), at 4 days after ALI (ALI+4) and at 12 days after ALI (ALI+12). MTEC cultures were stained with antibodies against GFP, against glutamylated tubulin to mark cilia and basal bodies, and with DAPI to mark nuclei. GFP is expressed from the FOXJ1 promoter early during ciliogenesis, but is not expressed pre-ALI. At ALI+4, most FOXJ1:GFP+ cells are undergoing centriole formation but have not formed cilia. At ALI+12, most FOXJ1:GFP+ cells are forming cilia or have completed ciliogenesis. Scale bars, 5 µm. (B) FACS analysis of cells dissociated from a trachea from a wild type mouse (left panel) and a trachea from a FOXJ1:GFP-expressing mouse (right panel). The red and blue rectangles are representative gates used to sort GFP+ and GFP- populations, respectively. (C) Sorted cell populations were stained with DAPI (blue) and combined acetylated alpha-tubulin and gamma-tubulin antibodies (green) to detect cilia and basal bodies. After sorting, 90–95% of cells in the GFP+ population at ALI+12 had observable cilia and or amplified basal bodies, whereas <1% of cells in the GFP- population stained positive for these markers, although some had a single primary cilium (arrowhead).

https://doi.org/10.1371/journal.pone.0052166.g001

The experimental design is shown in Figure 2A. The GFP+ MTECs represent a unique cell type that would be expected to have several transcriptional signatures, including those common to lung epithelial cells in general. To identify the signature specifically associated with ciliogenesis, ALI+4 (Table S1) and ALI+12 (Table S2) transcriptomes were generated by subtracting the non-ciliated, GFP- cell transcriptome (Table S3) from the transcriptomes of GFP+ ciliated cells at both timepoints. This resulted in a set of 649 unique upregulated genes at ALI+4 and 73 unique upregulated genes at ALI+12, where significant differential expression was defined as a fold-change cutoff of 2-fold or greater in addition to an adjusted P value <0.05 (Figure 2B). Importantly, subtraction decreased both the number of differentially expressed genes at ALI+4 and ALI+12 (Figure 2B) as well as the overlap between genes that are differentially expressed in both ciliated and non-ciliated cells relative to the reference (Figure 2C).

thumbnail
Download:
Figure 2. Subtractive analysis of MTEC microarray data.

(A) Experimental flowchart detailing microarray analysis. (B) The number of differentially expressed genes in ciliated MTECs at adj. p-values <0.05, <0.01 and <0.001 before subtraction (non-subtracted) and after subtraction. The non-subtracted gene list is the transcriptional profile of GFP+ cells compared to a universal mouse reference RNA. The subtracted gene list was derived analytically by comparing the transcriptional profile of GFP+ cells to that of GFP- cells from MTEC cultures. The number of genes that show 2-fold or greater differential expression and have adjusted p-values of 0.05, 0.01 and 0.001 are shown before and after subtractive analysis. After subtraction, 649 genes were significantly upregulated at ALI+4 with respect to GFP- cells, and 73 genes at ALI+12. 143 genes were downregulated at ALI+4 and none at ALI+12 (adj. p<0.05). (C) The overlap of significantly upregulated and downregulated genes (adj. p<0.05) from ALI+4, ALI+12 and GFP- pools decreases after subtraction. The GFP- group in both panels represents genes that are differentially expressed in non-ciliated cells relative to the universal reference RNA. M = log2(fold change).

https://doi.org/10.1371/journal.pone.0052166.g002

Functional Annotation of MTEC Transcriptome

We next assessed the effects of subtraction on the functional annotations recovered from these gene sets. Bioinformatics tools provided online at DAVID (Database for Annotation, Visualization and Integrated Discovery, http://david.abcc.ncifcrf.gov/) were used to mine known and predicted functional annotation information from databases including Gene Ontology (GO), SMART, INTERPRO and KEGG Pathway, for the subsets of genes expressed 5-fold or greater in FOXJ1:GFP- cells and FOXJ1:GFP+ cells at ALI+4 and ALI+12 (Table S4) [25]. As expected, FOXJ1:GFP- MTECs were enriched in genes associated with epithelial cell functions, including epithelial cell differentiation and cell adhesion, as well as cellular compartments common to all epithelial cells, such as the apical junction complex, desmosomes and tight junctions. The functional clusters most enriched in FOXJ1:GFP+ cells at ALI+12 were also associated with epithelial cell functions and structures. Additionally, the FOXJ1:GFP+ ALI+12 subset was enriched in components of the glutathione S-transferase pathway, which may be upregulated in these cells to aid in detoxification during the approximately 17–20 days of in vitro culture that the MTEC primary culture and differentiation protocol requires. However, after subtraction of the FOXJ1:GFP+ ALI+12 transcriptomic data, a strong enrichment in groups characteristic of cilia and microtubules was revealed, indicating that this analytical method successfully removes genes that are common to epithelial cells, leaving genes that are specific to ciliated cells. A similar strong enrichment in cilia and microtubule-related functional groups was observed in FOXJ1:GFP+ ALI+4 cells. The enrichment was apparent without subtraction, although subtraction increased the number of genes found in these clusters as well as the statistical significance associated with enrichment of these clusters. Based on these results, we focused our subsequent work on the subtracted dataset, and references to differentially expressed genes at ALI+4 and ALI+12 pertain to this subtracted data. The full dataset with subtracted and non-subtracted data can be accessed online at GEO (GEO Accession number GSE42500).

Comparison to Related Datasets

Several approaches have been taken to identify the protein components of centrosomes and cilia. These include comparative genomics and proteomics approaches, transcriptional assays, and proteomics analyses of motile and non-motile cilia and basal bodies and centrosomes in a variety of ciliated organisms including Trypanosoma, Chlamydomonas, Caenorhabditis, Drosophila, Dictyostelium, mice and humans [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], [45], [46]. We assessed the quality of our dataset through comparison to these other published datasets.

First, we assessed the overlap of highly upregulated genes from our ciliated MTEC dataset with two existing cilia proteome datasets, which are compendia of ten proteomic, genomic and comparative genomic studies of cilia components [47], [48]. These include a transcriptomic analysis of Chlamydomonas during flagellar regeneration [29], a process similar to growth of motile cilia in MTEC cultures, as well as a proteomic analysis of ciliary axonemes isolated from human bronchial epithelial cells [36]. In the ALI+4 dataset, 127 of 288 genes that were upregulated 10-fold or more (44%) were previously identified in one or more studies in the cilia proteome databases, and in the ALI+12 dataset, 67/123 (54%) of such genes were previously identified (Table S5). Thus, almost half of the genes that were upregulated 10-fold or more in ciliated MTECs at ALI+4 and ALI+12 were not found in any of the previous studies included in the cilia proteome databases. Second, we tested whether our dataset was enriched for targets of the FOXJ1 transcription factor, as predicted by our use of FOXJ1:GFP to isolate ciliating cells. Among 86 putative FOXJ1 targets identified in zebrafish and Xenopus [49], [50], the mouse orthologs of 42 of these were identified as upregulated in the ALI+4 dataset, and 31 of these were upregulated in the ALI+12 dataset (Table S5). Lastly, we interrogated our dataset using a test set of 99 mouse genes predicted to be linked to cilium formation by their expression in tissues that contain numerous ciliated cells such as olfactory epithelium, testes, vomeronasal organ, trachea, and lung [46]. 39 out of the 99 genes identified by tissue expression patterns were significantly upregulated (p<0.05) at ALI+4, and 16 out of 99 genes were significantly upregulated at ALI+12 (Table S5). Interestingly, when we performed the reverse analysis, we found that only 6% of upregulated genes from ALI+4 were previously identified in this study, whereas 21% of upregulated genes from ALI+12 were previously identified. This likely reflects the fully differentiated state of the majority of ciliated cells found in the tissues assayed in previous studies, and suggests that our dataset may hold clues to the early stages of ciliogenesis that have not been revealed by other studies.

In Situ Localization of Upregulated Proteins in the Human Protein Atlas

We next assessed which of the proteins encoded by genes that were upregulated in our ciliated MTEC datasets can be detected at the cilia or the basal bodies of multiciliated cells. To do so, we scanned immunohistochemistry data from the Human Protein Atlas (HPA), an online compendium of tissue microarray data from normal tissues, pathology specimens and cell lines [51]. The HPA dataset consists of tissue sections stained with antibodies raised against human proteins. HPA data were recently used to validate results from a mass spectrometry-based proteomics screen for candidate centrosome proteins [43], as well as a second in silico identification of cilia proteins based on tissue RNA expression [52]. We used this tool to evaluate MTEC candidate genes.

The HPA includes tissue sections from bronchus, nasopharynx and the fallopian tube, which are tissues containing multiciliated epithelia. We found that we could discern four patterns of antibody staining in these sections relevant to the structure and function of ciliated epithelial cells: 1) cilia, visible as a tuft on the apical surface; 2) basal body layer, visible as a line at the junction of the ciliary tuft and the cell body, 3) the nucleus of ciliated cells, 4) whole cell staining of ciliated cells, relative to adjacent non-ciliated cells. We analyzed 8,038 HPA images representing 923 proteins whose orthologous mouse proteins are encoded by genes that were upregulated 4-fold or more in ciliated MTEC (Table S6). Of the 923 proteins represented, 504 (55%) revealed staining of ciliated cells: 195 (21%) localized to basal bodies, 163 (17%) localized to cilia, 22 (2%) localized to ciliated cell nuclei and 143 (15%) localized to other regions of ciliated cells in at least one of the three tissues examined. These data confirm by an independent method that a high percentage of genes detected as differentially expressed in ciliated MTECs by our subtractive transcriptional analysis localize to structures specifically in the ciliated epithelial cells of multiple human tissues.

Functional and Physical Association of Proteins Upregulated during Ciliogenesis

Other tools we used to evaluate the ciliated MTEC transcriptomic data were Mitocheck (www.Mitocheck.org), the Human Protein Reference Database (HPRD) [53] and the biological general repository for interaction datasets (BioGrid) repository [54]. Previously uncharacterized proteins that were upregulated during ciliogenesis were organized by known physical and genetic interactions (File S1). Mapping in this manner revealed nodes of genes with defined functional or structural links to cilia or centrosomes, such as cell cycle regulators (Aurka, Brca1, Ccna1, Figure 3A), hedgehog signaling components (Stk36, Gli2, Figure 3B), and the PCM1 complex (Pcm1, Bbs4, Hap1, Figure 3C) [55], [56], [57], [58], [59], [60]. Based on these results, we identified the uncharacterized gene Klhdc9 to be a potential candidate for the cell cycle control of ciliogenesis, and Mast2 to participate in hedgehog signaling. Additionally, the cluster containing the PCM1 complex contains AZI1, also known as CEP131 in humans and Dilatory (DILA) in flies, which are proteins that are involved in ciliogenesis [61]–[62]. AZI1, like PCM1 and BBS4 [63], [64], has recently been confirmed to be a component of centriolar satellites [65]. This example illustrates how combining the transcriptional data on MTEC ciliogenesis with publicly available and curated protein interaction data can help to identify bona fide components of structural and functional complexes.

thumbnail
Download:
Figure 3. Functional mapping of candidate centrosome/cilia proteins and identification of novel centrosome components.

Candidate genes were organized by known physical and genetic interactions mined from Mitocheck, the Human Protein Reference Database and BioGrid. Clusters containing genes associated with cell cycle control (A), hedgehog signaling (B) and the PCM1 complex (C) are shown.

https://doi.org/10.1371/journal.pone.0052166.g003

A Common Set of Genes Associated with Centriole Duplication and Basal Body Generation

The basal bodies of ciliated epithelial cells are morphologically similar to the centrioles in cycling cells, but the degree of that similarity with respect to composition and mechanism of formation is unclear. To explore this relationship, we first asked whether the genes encoding a core set of proteins known to be required for centriole duplication in cycling cells are upregulated during ciliogenesis. This set includes PLK4, SAS-6, CENPJ/SAS-4, CEP152, CEP120, CEP135, CP110, centrobin, delta-tubulin and epsilon-tubulin, all of which are physically associated with centrioles [18], [66], [67], [68]. All of the genes encoding these proteins were upregulated during MTEC ciliogenesis (Figure 4A), suggesting that centriole formation in cycling cells and basal body duplication in ciliated epithelial cells share a common mechanism. As predicted, the degree of the upregulation for core centriolegenesis genes was greater at ALI+4 when the majority of ciliating cells are generating basal bodies than at ALI+12, when centriolegenesis is complete in most cells [69]. Most striking among these results was that Plk4, a master regulator of centriole formation ([12], [70], [71]), was upregulated more than 20-fold at the ALI+4 timepoint (Figure 3A).

thumbnail
Download:
Figure 4. Genes encoding centrosome proteins are differentially regulated during basal body formation.

(A) Upregulation of regulatory and structural proteins required for centriole duplication in cycling cells and (B) heat map of putative and known centrosomal components. Genes are shown in rows and replicate arrays are shown in columns. Three biological replicates were performed for ALI+4, ALI+12 and GFP-. In addition, two technical replicates were performed for two of the ALI+4 biological replicates (columns 1,2 and 4,5) for a total of five columns. Data were zero-transformed against non-ciliated (GFP-) cells. The scale indicates the fold change in expression (log2). Grey boxes represent gene spots that failed to pass quality control filters for the indicated array.

https://doi.org/10.1371/journal.pone.0052166.g004

In addition to genes involved in centriole duplication, many other genes encoding known centrosomal proteins were also upregulated in the ciliated MTEC dataset (Figure 4B). In some cases the upregulated proteins had only previously been reported to play a role in centrosome functions in cycling cells seemingly unrelated to basal body function. For example, CEP63 and CEP72 have been shown to function in formation of bipolar spindles in mitotic cells [72], [73]. We speculate that the upregulation of such proteins during ciliogenesis in ciliated MTEC reflects a cilium-related function in the interphase centrosome of cycling cells, in which the mother centriole acts as a basal body for the primary cilium. Some centrosomal components identified in cycling cells were not induced during MTEC ciliogenesis. These include CEP170, CEP192, CEP250/C-NAP1, CEP27/HAUS2 and CEP68, all of which are thought to have functions associated with centrosomes, and not basal bodies, including microtubule and spindle organization, recruitment of pericentriolar material and centrosome cohesion [74], [75], [76], [77], [78], [79]. Similarly, the gamma-tubulin ring complex mediates microtubule nucleation at the centrosome [80], [81], and most of the members of this complex were not significantly upregulated relative to non-ciliated cells, consistent with the notion that a limited subset of centrosome functions is associated with the centrioles serving as basal bodies in ciliated MTECs.

In the interest of identifying new components of the centrosome in our dataset, we reasoned that such genes might have the pattern of high upregulation at ALI+4, during basal body formation, and lower expression at ALI+12, after ciliogenesis is complete. TTC12 (tetratricopeptide repeat domain 12) is one such protein that fits this criterion and was previously identified by comparative genomics of ciliated and non-ciliated organisms [30]. Ttc12 was upregulated 2.3-fold at ALI+4 and decreased to 1.5-fold at ALI+12 (Tables S1, S2). To test whether TTC12 is centrosome protein, TTC12-GFP was expressed in NIH/3T3 cells. We found that TTC12 localizes to the centrosome as defined by gamma-tubulin labeling (Figure 5), demonstrating that the transcriptional profile of ciliated MTEC can be used to identify novel centrosome components.

thumbnail
Download:
Figure 5. Identification of a novel centrosome component.

NIH/3T3 cells were transfected with a plasmid encoding TTC12-GFP. TTC12-GFP localizes to the centrosome as shown by overlap with gamma-tubulin staining. Green, TTC12-GFP. Red, gamma-tubulin.

https://doi.org/10.1371/journal.pone.0052166.g005

Identification of Novel Links to Human Disease

Historically, the human diseases associated with cilia have been considered as two separate groups: those associated with defects in the beating of motile cilia, and those associated with non-motile cilia, i.e. the primary cilium. It is unclear whether the proteins associated with the signaling functions of non-motile cilia are unique to those cilia; however, the motile flagella of Chlamydomonas and its components clearly have signaling functions [82], [83], [84], and motile cilia in mammalian ciliated epithelia are competent for specific signaling events [85], [86]. We examined the transcriptional data for genes with known links to diseases of motile and non-motile cilia. As expected from the source of cells in our experiments, many of the genes associated with motile cilium diseases were upregulated in the ciliated MTEC dataset (Figure 6 top panel). Notably, many genes associated with non-motile cilium diseases were also upregulated, including Rpgrip1l/Nphp8, Alms1, Nphp1 and Nphp4, Mks8, Ofd1 and Lca5 (Figure 6 middle panel). This suggests that these genes encode components common to nonmotile and motile cilia, and that mutations in them might also display a phenotype in motile cilia. This is consistent with the results for the Bardet-Biedl syndrome (BBS) genes, most of which are similarly upregulated in the ciliated MTEC dataset, and have been linked to defects in both motile and non-motile cilia [7] (Figure 6 bottom panel).

thumbnail
Download:
Figure 6. Differential regulation of disease genes associated with motile and nonmotile ciliopathies.

Heat maps showing the expression of genes associated with motile (top) and nonmotile (middle) ciliopathies, as well as genes linked to Bardet-Biedl syndrome (bottom), a model ciliopathy.

https://doi.org/10.1371/journal.pone.0052166.g006

We examined three cases in which a gene upregulated in ciliated MTEC had been previously associated with a mammalian disease phenotype, but not yet associated with motile cilia, primary cilia, or the centrosome.

Mdm1.

The connecting cilium of the photoreceptor outer segment is a specialized type of cilium that shares many features and components with motile cilia and primary cilia [87], and many genes liked to retinopathy encode proteins that localize to the connecting cilium [88]. A positional cloning approach identified a nonsense mutation in Mdm1 in the arrd2 mouse model for age-related retinal degeneration [89], indicating the Mdm1 gene product may be important for retinal function. Mdm1, which is >60-fold upregulated in ciliated MTECs (Tables S1, S2), has also been identified in proteomic analysis of human centrosomes [38] and mouse photoreceptor sensory cilia [90] and is conserved in organisms with ciliated epithelial cells, suggesting that it is a component of cilia and or basal bodies. To test whether MDM1 is a centrosome/cilia component we expressed a MDM1-GFP fusion protein in NIH/3T3 cells. MDM1-GFP localized to centrosomes, as determined by colocalization with gamma-tubulin (Figure 7A) and to the primary cilium, as determined by colocalization with acetylated alpha-tubulin (Figure 7B). Based on these results we suggest that the age-related retinal degeneration phenotype observed in the aard2 mouse results from a defect in connecting cilium function.

thumbnail
Download:
Figure 7. Identification of new links between centrosomes, cilia and human disease.

(A, B) NIH/3T3 cells were transfected with a plasmid encoding MDM1-GFP. MDM1-GFP localizes to centrosomes (A) and the primary cilium (B). NIH/3T3 cells and MTECs were stained for MLF1 with an anti-MLF1 antibody. MLF1 localizes to the primary cilium in NIH/3T3 cells (C) and motile cilia in MTECs (D). Insets show higher magnification of merged images, slightly offset as indicated by circles in upper left. Scale bars, 5 µm. Green, GFP (A, B); MLF1 (C, D). Red, gamma-tubulin (A); acetylated tubulin (B); polaris (C, D).

https://doi.org/10.1371/journal.pone.0052166.g007

Mlf1.

Mlf1 (myeloid leukemia factor 1) is associated with myeloid leukemia; Mlf1 is fused to nucleophosmin (Npm1) in 1% of acute myeloid leukemia (AML) cases [91] and is overexpressed in 25% of cases of AML associated with myelodysplastic syndrome [92]. Based on observed interactions with p53, it has been suggested that defects or dysregulation of Mlf1 might perturb cell cycle regulation and contribute to malignant transformation [93]. We identified Mlf1 as being >30-fold upregulated in ciliated MTECs (Tables S1, S2). In addition, MLF1 has been identified in proteomic analyses of the centrosome and connecting photoreceptor cilium [36], [90]. To test the localization of MLF1, we stained NIH/3T3 cells for MLF1 using a commercially available antibody (H00004291-B01, Abnova). Faint MLF1 staining was observed in primary cilia, as determined by colocalization with the intraflagellar transport component polaris (Figure 7C). MLF1 also labeled the length of motile cilia in MTECs (Figure 7D) and in some cells appeared to be concentrated at the apical tips (not shown). These results indicate that MLF1, a putative regulator of p53 and the cell cycle, is a component of both primary and motile cilia. Interestingly, the closely related protein MLF2 was identified in ciliary proteomic studies [36], [90], and tissue labeling in the HPA indicates that MLF2 localizes to cilia in ciliated epithelial cells (Table S6). We also note also that the putative Drosophila ortholog of MLF1/2 interacts with Su(fu), a component of the hedgehog signaling pathway [94] that localizes to the primary cilium in mammals [59].

Dyx1c1.

Neurocognitive impairments have been observed in individuals with ciliopathies, but the mechanistic links between the neurodevelopmental and neurocognitive diseases and cilium function are poorly understood [88]. Dyslexia is one such neurocognitive impairment that has recently been linked to primary cilia through the localization of the product of the candidate dyslexia susceptibility gene Dcdc2 to primary cilia in primary rat neurons, and the effect of its overexpresion on cilia length [95]. Several genetic loci in addition to DCDC2 have been linked to dyslexia susceptibility [96]. We noted that nine of ten candidate dyslexia genes present in three of the dyslexia-associated genomic regions [97], [98], [99] were differentially regulated in ciliated MTECs (Figure 8A). One of these genes, Dyx1c1 (Dyslexia Susceptibility Candidate 1), is linked both genetically and mechanistically to dyslexia [99], [100]. Dyx1c1 was upregulated more than 16-fold at ALI+4 in ciliated MTECs. It also was previously identified in a comparative genomics screen for genes that are found only in organisms having cilia [30], in an in silico study that used public microarray data to identify potential cilia components [45], and identified as interacting with the cytoskeleton in an assessment of molecular networks [101]. Based on these results, we examined the localization of GFP-tagged versions of rat DYX1C1 [102] that we expressed in NIH/3T3 cells. DYX1C1-GFP colocalized with gamma-tubulin at the centrosome (Figure 8B). Although the centrosome localization was most prominent, DYX1C1-GFP also localized to the primary cilium in some cells (Figure 8C). These results suggest that the neuronal migration defect of DYX1C1 depletion may reflect a function of this protein in the centrosome or cilium, adding to the growing list of such proteins that have been linked to diseases of the brain [103].

thumbnail
Download:
Figure 8. A dyslexia candidate disease gene is transcriptionally upregulated and localizes to the centrosome.

(A) Heat maps depicting the expression of dyslexia candidate genes in ciliated MTECs. Dyx1c1 is among the most highly upregulated genes, and Robo1 is strongly downregulated. (B-D) A DYX1C1-GFP fusion protein was expressed in NIH/3T3 cells. DYX1C1-GFP localizes to the centrosome (B) and in some cells the primary cilium (C). Insets show positively staining regions at higher magnification. Scale bars, 5 um. Green, GFP (B, C). Red, gamma-tubulin (B); acetylated tubulin (C).

https://doi.org/10.1371/journal.pone.0052166.g008

Discussion

We determined the transcriptional profile of ciliated tracheal epithelial cells undergoing differentiation in vitro. Previous studies focused on ciliated cells or single-cell organisms, including analyses of human bronchial epithelia [44], brain ependyma [42] and of Chlamydomonas flagellar regeneration [29] have yielded useful information, however, the combination of the unique properties of the ciliated MTECs with the cell-type specificity achieved by cell sorting make this a particularly powerful dataset. The upregulated genes comprise a dataset that both reveals aspects of function in this unique ciliated cell type, and identifies new proteins of the centriole, centrosome and cilium that include disease-associated proteins that were not previously known to be associated with centrosome and cilium defects. Here we consider the implications of these results for centrosome and cilium biology and human disease.

Ciliated epithelial cells produce hundreds of centrioles and motile cilia during the course of differentiation, whereas cycling cells produce only two centrioles and one non-motile primary cilium per cell cycle. Centriole formation in ciliated epithelial cells has been shown by electron microscopy to be morphologically distinct from that in cycling cells, with many centrioles forming simultaneously in close proximity to deuterosomes, structures unique to differentiating ciliated epithelial cells [19]. Although several previous studies have shown that particular centriole and centrosome components are upregulated during ciliated epithelial cell differentiation, the relationship of the critical regulators of centriole formation in cycling cells to the amplification process in ciliated epithelial cells was not clear. Our results show that the core set of proteins required for centriole duplication in cycling cells is upregulated during basal body formation in ciliated MTECs. Most importantly, this includes Cyclin A (Ccna1), which is required for centriole formation [104], PLK4 and SAS-6, two proteins known to be capable of inducing the formation of multiple centrioles when overexpressed [71], [105], and CEP152/ASL, which is required for normal centriole duplication and PLK4-induced overduplication in cultured human and fly cells [106], [107]. Consistent with the notion that upregulation of genes reflects a role in the unique events of ciliated MTEC differentiation, we have previously shown that SAS-6 is required for centriole formation in ciliated MTECs [17]. We suggest that basal body formation in multiciliated cells likely occurs by the same mechanism as in cycling cells, but that the tight control over centriole number in cycling cells is bypassed by strong upregulation of PLK4, allowing the formation of many centrioles at once. However, it remains possible that some unique aspects of basal body formation, such as the presence of deuterosomes, might involve gene products specifically expressed in multiciliated cells.

We found that expression of many genes associated with centriole formation was highest at ALI+4 when centrioles, or basal bodies, were forming and declined at ALI+12 when centriole formation was largely completed. We used this as a screening criterion for the identification of new centrosome components amongst the previously uncharacterized genes in the dataset. In this first-pass analysis, TTC12 was identified to be a new centrosome component. In addition, we used data from the Human Protein Atlas to assess the localization of proteins from 923 upregulated genes, and found that approximately 45% localized specifically to multiciliated epithelial cells in the tissues examined, and approximately 38% localized to the basal body layer, to cilia, or to both. Since many of these proteins were previously uncharacterized these results greatly expand the number of potential components of the centrioles and cilium. Although some of these will be specific to motile cilia or basal bodies in multiciliated cells, it is likely that many are general components of centrioles and cilia, based on the high representation of known centriolar and ciliary proteins in the dataset, and on the high frequency with which genes identified here have been identified in at least one of the other proteomic and genomic datasets associated with the centrosome and/or primary cilium. We also drew on protein interaction information from publicly available databases to map out known interactions between proteins that are upregulated during ciliogenesis, most of which have unknown function. The value of such an analysis was demonstrated by correctly identifying an interaction between centriolar satellite protein AZI1/CEP131 and the PCM1 complex, all of which are highly expressed during differentiation of multiciliated cells.

Among the genes upregulated during ciliogenesis in ciliated MTEC, we have identified several associated with human disorders. The disorders include those known to be associated with centrosome/cilium dysfunction, such as retinal degeneration, but also others not previously associated, such as leukemia and dyslexia. Neurocognitive deficits in general are observed in patients with mutations affecting centrosome or cilium function, ranging in severity from microcephaly and lissencephaly, in which the brain or segments of it are greatly reduced in size or complexity, to more subtle neurological impairments such as schizophrenia and bipolar disorder [108], [109]. In this context it would not perhaps be surprising that a reading disability would have a centrosome/cilium defect as its root cause. We show here that DYX1C1 is physically associated with the centrosome and the cilium, and also note that another strong dyslexia candidate gene from human genetic studies, DCDC2, encodes a protein that is a component of the primary cilium that can impact cilia length when overexpressed [95]. The molecular basis for dyslexia is unknown, but these candidates have both been shown to be required for neuronal migration in the cortex, providing one possible mechanism for development of the disease [110].

The set of human diseases referred to as ciliopathies has received much attention recently, and their analysis is an outstanding example of the reciprocal interaction between clinical and basic research. We note that the ciliated MTEC transcriptome dataset described here might be particularly useful for extending the correlations between human phenotype and centrosome/cilium function because it is based on a mammalian gene set, and is focused on a cell type that has an extreme dependence on the centrosome/cilium for its physiological function. Here we have established connections between genes identified by human or other mammalian genetic studies, upregulation of the genes in this specialized tissue, and localization of the cognate proteins to the centrosome/cilium. Future work will determine whether the pathologies caused by mutation of these genes result from defects in centrosome/cilium function, or alternatively whether localization of the disease-associated protein to the centrosome/cilium is important for the function of that particular protein or pathway.

Supporting Information

Table S1.

Subtracted ALI+4 transcriptome.

https://doi.org/10.1371/journal.pone.0052166.s001

(XLSX)

Table S2.

Subtracted ALI+12 transcriptome.

https://doi.org/10.1371/journal.pone.0052166.s002

(XLSX)

Table S3.

GFP- transcriptome.

https://doi.org/10.1371/journal.pone.0052166.s003

(XLSX)

Table S4.

Top 5 functional annotation clusters in FOXJ1:GFP + (ALI+4 and ALI+12) and FOXJ1:GFP - MTECs.

https://doi.org/10.1371/journal.pone.0052166.s004

(XLSX)

Table S5.

Overlap of the ciliated MTEC transcriptome with other proteomic, transcriptomic and genomic studies of the cilium and centrosome.

https://doi.org/10.1371/journal.pone.0052166.s005

(XLSX)

Table S6.

Mining the Human Protein Atlas for localization of human orthologs of candidate ciliated MTEC proteins in multiciliated human tissues.

https://doi.org/10.1371/journal.pone.0052166.s006

(XLSX)

File S1.

Inferring function of candidate ciliogenesis proteins by mapping known physical and genetic interactions.

https://doi.org/10.1371/journal.pone.0052166.s007

(CYS)

Acknowledgments

We thank L. Ostrowski for the FOXJ1:GFP mouse strain, J. LoTurco and B. Yoder for antibodies, Tim Knaak for assistance with flow cytometry, B. Daigle for advice on microarray experimental design, and J. Demeter for technical assistance with microarray analysis software.

Author Contributions

Conceived and designed the experiments: RAH TRS. Performed the experiments: RAH TRS. Analyzed the data: RAH TRS ET TS. Contributed reagents/materials/analysis tools: RAH TRS ET. Wrote the paper: RAH TS.

References

  1. 1. Vincensini L, Blisnick T, Bastin P (2011) 1001 model organisms to study cilia and flagella. Biol Cell 103: 109–130.
  2. 2. Berbari NF, O’Connor AK, Haycraft CJ, Yoder BK (2009) The primary cilium as a complex signaling center. Curr Biol 19: R526–535.
  3. 3. Goetz SC, Anderson KV (2010) The primary cilium: a signalling centre during vertebrate development. Nat Rev Genet 11: 331–344.
  4. 4. Fisch C, Dupuis-Williams P (2011) Ultrastructure of cilia and flagella - back to the future! Biol Cell. 103: 249–270.
  5. 5. Christensen ST, Guerra CF, Awan A, Wheatley DN, Satir P (2003) Insulin receptor-like proteins in Tetrahymena thermophila ciliary membranes. Curr Biol 13: R50–52.
  6. 6. Pan J, Snell WJ (2000) Regulated targeting of a protein kinase into an intact flagellum. An aurora/Ipl1p-like protein kinase translocates from the cell body into the flagella during gamete activation in chlamydomonas. J Biol Chem 275: 24106–24114.
  7. 7. Shah AS, Farmen SL, Moninger TO, Businga TR, Andrews MP, et al. (2008) Loss of Bardet-Biedl syndrome proteins alters the morphology and function of motile cilia in airway epithelia. Proc Natl Acad Sci U S A 105: 3380–3385.
  8. 8. Bloodgood RA (2010) Sensory reception is an attribute of both primary cilia and motile cilia. J Cell Sci 123: 505–509.
  9. 9. Basto R, Lau J, Vinogradova T, Gardiol A, Woods CG, et al. (2006) Flies without centrioles. Cell 125: 1375–1386.
  10. 10. Marshall WF (2009) Centriole evolution. Curr Opin Cell Biol 21: 14–19.
  11. 11. Uetake Y, Loncarek J, Nordberg JJ, English CN, La Terra S, et al. (2007) Cell cycle progression and de novo centriole assembly after centrosomal removal in untransformed human cells. J Cell Biol 176: 173–182.
  12. 12. Bettencourt-Dias M, Rodrigues-Martins A, Carpenter L, Riparbelli M, Lehmann L, et al. (2005) SAK/PLK4 is required for centriole duplication and flagella development. Curr Biol 15: 2199–2207.
  13. 13. Dawe HR, Farr H, Gull K (2007) Centriole/basal body morphogenesis and migration during ciliogenesis in animal cells. J Cell Sci 120: 7–15.
  14. 14. Nigg EA, Raff JW (2009) Centrioles, centrosomes, and cilia in health and disease. Cell 139: 663–678.
  15. 15. Badano JL, Mitsuma N, Beales PL, Katsanis N (2006) The ciliopathies: an emerging class of human genetic disorders. Annu Rev Genomics Hum Genet 7: 125–148.
  16. 16. Tobin JL, Beales PL (2009) The nonmotile ciliopathies. Genet Med 11: 386–402.
  17. 17. Vladar EK, Stearns T (2007) Molecular characterization of centriole assembly in ciliated epithelial cells. J Cell Biol 178: 31–42.
  18. 18. Mahjoub MR, Xie Z, Stearns T (2010) Cep120 is asymmetrically localized to the daughter centriole and is essential for centriole assembly. J Cell Biol 191: 331–346.
  19. 19. Dirksen ER (1991) Centriole and basal body formation during ciliogenesis revisited. Biol Cell 72: 31–38.
  20. 20. Hagiwara H, Ohwada N, Takata K (2004) Cell biology of normal and abnormal ciliogenesis in the ciliated epithelium. Int Rev Cytol 234: 101–141.
  21. 21. Leidel S, Delattre M, Cerutti L, Baumer K, Gonczy P (2005) SAS-6 defines a protein family required for centrosome duplication in C. elegans and in human cells. Nat Cell Biol 7: 115–125.
  22. 22. Ostrowski LE, Hutchins JR, Zakel K, O’Neal WK (2003) Targeting expression of a transgene to the airway surface epithelium using a ciliated cell-specific promoter. Mol Ther 8: 637–645.
  23. 23. You Y, Richer EJ, Huang T, Brody SL (2002) Growth and differentiation of mouse tracheal epithelial cells: selection of a proliferative population. Am J Physiol Lung Cell Mol Physiol 283: L1315–1321.
  24. 24. Wettenhall JM, Smyth GK (2004) limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics 20: 3705–3706.
  25. 25. Huang da W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57.
  26. 26. You Y, Huang T, Richer EJ, Schmidt JE, Zabner J, et al. (2004) Role of f-box factor foxj1 in differentiation of ciliated airway epithelial cells. Am J Physiol Lung Cell Mol Physiol 286: L650–657.
  27. 27. Pezzulo AA, Starner TD, Scheetz TE, Traver GL, Tilley AE, et al. (2011) The air-liquid interface and use of primary cell cultures are important to recapitulate the transcriptional profile of in vivo airway epithelia. Am J Physiol Lung Cell Mol Physiol 300: L25–31.
  28. 28. Pazour GJ, Agrin N, Leszyk J, Witman GB (2005) Proteomic analysis of a eukaryotic cilium. J Cell Biol 170: 103–113.
  29. 29. Stolc V, Samanta MP, Tongprasit W, Marshall WF (2005) Genome-wide transcriptional analysis of flagellar regeneration in Chlamydomonas reinhardtii identifies orthologs of ciliary disease genes. Proc Natl Acad Sci U S A 102: 3703–3707.
  30. 30. Li JB, Gerdes JM, Haycraft CJ, Fan Y, Teslovich TM, et al. (2004) Comparative genomics identifies a flagellar and basal body proteome that includes the BBS5 human disease gene. Cell 117: 541–552.
  31. 31. Avidor-Reiss T, Maer AM, Koundakjian E, Polyanovsky A, Keil T, et al. (2004) Decoding cilia function: defining specialized genes required for compartmentalized cilia biogenesis. Cell 117: 527–539.
  32. 32. Blacque OE, Perens EA, Boroevich KA, Inglis PN, Li C, et al. (2005) Functional genomics of the cilium, a sensory organelle. Curr Biol 15: 935–941.
  33. 33. Broadhead R, Dawe HR, Farr H, Griffiths S, Hart SR, et al. (2006) Flagellar motility is required for the viability of the bloodstream trypanosome. Nature 440: 224–227.
  34. 34. Efimenko E, Bubb K, Mak HY, Holzman T, Leroux MR, et al. (2005) Analysis of xbx genes in C. elegans. Development 132: 1923–1934.
  35. 35. Keller LC, Romijn EP, Zamora I, Yates JR, 3rd, Marshall WF (2005) Proteomic analysis of isolated chlamydomonas centrioles reveals orthologs of ciliary-disease genes. Curr Biol 15: 1090–1098.
  36. 36. Ostrowski LE, Blackburn K, Radde KM, Moyer MB, Schlatzer DM, et al. (2002) A proteomic analysis of human cilia: identification of novel components. Mol Cell Proteomics 1: 451–465.
  37. 37. Schulz I, Reinders Y, Sickmann A, Graf R (2006) An improved method for Dictyostelium centrosome isolation. Methods Mol Biol 346: 479–489.
  38. 38. Andersen JS, Wilkinson CJ, Mayor T, Mortensen P, Nigg EA, et al. (2003) Proteomic characterization of the human centrosome by protein correlation profiling. Nature 426: 570–574.
  39. 39. Chen N, Mah A, Blacque OE, Chu J, Phgora K, et al. (2006) Identification of ciliary and ciliopathy genes in Caenorhabditis elegans through comparative genomics. Genome Biol 7: R126.
  40. 40. Baron DM, Ralston KS, Kabututu ZP, Hill KL (2007) Functional genomics in Trypanosoma brucei identifies evolutionarily conserved components of motile flagella. J Cell Sci 120: 478–491.
  41. 41. Merchant SS, Prochnik SE, Vallon O, Harris EH, Karpowicz SJ, et al. (2007) The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science 318: 245–250.
  42. 42. Jacquet BV, Salinas-Mondragon R, Liang H, Therit B, Buie JD, et al. (2009) FoxJ1-dependent gene expression is required for differentiation of radial glia into ependymal cells and a subset of astrocytes in the postnatal brain. Development 136: 4021–4031.
  43. 43. Jakobsen L, Vanselow K, Skogs M, Toyoda Y, Lundberg E, et al. (2011) Novel asymmetrically localizing components of human centrosomes identified by complementary proteomics methods. EMBO J 30: 1520–1535.
  44. 44. Ross AJ, Dailey LA, Brighton LE, Devlin RB (2007) Transcriptional profiling of mucociliary differentiation in human airway epithelial cells. Am J Respir Cell Mol Biol 37: 169–185.
  45. 45. Ivliev AE, ‘t Hoen PAC, van Roon-Mom WMC, Peters DJM, Sergeeva MG (2012) Exploring the Transcriptome of Ciliated Cells Using In Silico Dissection of Human Tissues. PLOS One 7: e35618.
  46. 46. McClintock TS, Glasser CE, Bose SC, Bergman DA (2008) Tissue expression patterns identify mouse cilia genes. Physiol Genomics 32: 198–206.
  47. 47. Gherman A, Davis EE, Katsanis N (2006) The ciliary proteome database: an integrated community resource for the genetic and functional dissection of cilia. Nat Genet 38: 961–962.
  48. 48. Inglis PN, Boroevich KA, Leroux MR (2006) Piecing together a ciliome. Trends Genet 22: 491–500.
  49. 49. Stubbs JL, Oishi I, Izpisua Belmonte JC, Kintner C (2008) The forkhead protein Foxj1 specifies node-like cilia in Xenopus and zebrafish embryos. Nat Genet 40: 1454–1460.
  50. 50. Yu X, Ng CP, Habacher H, Roy S (2008) Foxj1 transcription factors are master regulators of the motile ciliogenic program. Nat Genet 40: 1445–1453.
  51. 51. Barbe L, Lundberg E, Oksvold P, Stenius A, Lewin E, et al. (2008) Toward a confocal subcellular atlas of the human proteome. Mol Cell Proteomics 7: 499–508.
  52. 52. Ivliev AE, t Hoen PA, van Roon-Mom WM, Peters DJ, Sergeeva MG (2012) Exploring the transcriptome of ciliated cells using in silico dissection of human tissues. PLOS One 7: e35618.
  53. 53. Prasad TS, Kandasamy K, Pandey A (2009) Human Protein Reference Database and Human Proteinpedia as discovery tools for systems biology. Methods Mol Biol 577: 67–79.
  54. 54. Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, et al. (2006) BioGRID: a general repository for interaction datasets. Nucleic Acids Res 34: D535–539.
  55. 55. Dutertre S, Descamps S, Prigent C (2002) On the role of aurora-A in centrosome function. Oncogene 21: 6175–6183.
  56. 56. Lacey KR, Jackson PK, Stearns T (1999) Cyclin-dependent kinase control of centrosome duplication. Proc Natl Acad Sci U S A 96: 2817–2822.
  57. 57. Deng CX (2002) Roles of BRCA1 in centrosome duplication. Oncogene 21: 6222–6227.
  58. 58. Wilson CW, Nguyen CT, Chen MH, Yang JH, Gacayan R, et al. (2009) Fused has evolved divergent roles in vertebrate Hedgehog signalling and motile ciliogenesis. Nature 459: 98–102.
  59. 59. Haycraft CJ, Banizs B, Aydin-Son Y, Zhang Q, Michaud EJ, et al. (2005) Gli2 and Gli3 localize to cilia and require the intraflagellar transport protein polaris for processing and function. PLoS Genet 1: e53.
  60. 60. Engelender S, Sharp AH, Colomer V, Tokito MK, Lanahan A, et al. (1997) Huntingtin-associated protein 1 (HAP1) interacts with the p150Glued subunit of dynactin. Hum Mol Genet 6: 2205–2212.
  61. 61. Wilkinson CJ, Carl M, Harris WA (2009) Cep70 and Cep131 contribute to ciliogenesis in zebrafish embryos. BMC Cell Biol 10: 17.
  62. 62. Ma L, Jarman AP (2011) Dilatory is a Drosophila protein related to AZI1 (CEP131) that is located at the ciliary base and required for cilium formation. J Cell Sci 124: 2622–2630.
  63. 63. Kubo A, Sasaki H, Yuba-Kubo A, Tsukita S, Shiina N (1999) Centriolar satellites: molecular characterization, ATP-dependent movement toward centrioles and possible involvement in ciliogenesis. J Cell Biol 147: 969–980.
  64. 64. Kim JC, Badano JL, Sibold S, Esmail MA, Hill J, et al. (2004) The Bardet-Biedl protein BBS4 targets cargo to the pericentriolar region and is required for microtubule anchoring and cell cycle progression. Nat Genet 36: 462–470.
  65. 65. Staples CJ, Myers KN, Beveridge RD, Patil AA, Lee AJ, et al.. (2012) The centriolar satellite protein Cep131 is important for genome stability. J Cell Sci.
    • 66. Bettencourt-Dias M, Glover DM (2009) SnapShot: centriole biogenesis. Cell 136: 188–188 e181.
      • 67. Hatch EM, Kulukian A, Holland AJ, Cleveland DW, Stearns T (2010) Cep152 interacts with Plk4 and is required for centriole duplication. J Cell Biol 191: 721–729.
      • 68. Nigg EA, Stearns T (2011) The centrosome cycle: Centriole biogenesis, duplication and inherent asymmetries. Nat Cell Biol 13: 1154–1160.
      • 69. Vladar EK, Stearns T (2007) Molecular characterization of centriole assembly in ciliated epithelial cells. J Cell Biol 178: 31–42.
      • 70. Habedanck R, Stierhof YD, Wilkinson CJ, Nigg EA (2005) The Polo kinase Plk4 functions in centriole duplication. Nat Cell Biol 7: 1140–1146.
      • 71. Kleylein-Sohn J, Westendorf J, Le Clech M, Habedanck R, Stierhof YD, et al. (2007) Plk4-induced centriole biogenesis in human cells. Dev Cell 13: 190–202.
      • 72. Smith E, Dejsuphong D, Balestrini A, Hampel M, Lenz C, et al. (2009) An ATM- and ATR-dependent checkpoint inactivates spindle assembly by targeting CEP63. Nat Cell Biol 11: 278–285.
      • 73. Oshimori N, Li X, Ohsugi M, Yamamoto T (2009) Cep72 regulates the localization of key centrosomal proteins and proper bipolar spindle formation. EMBO J 28: 2066–2076.
      • 74. Guarguaglini G, Duncan PI, Stierhof YD, Holmstrom T, Duensing S, et al. (2005) The forkhead-associated domain protein Cep170 interacts with Polo-like kinase 1 and serves as a marker for mature centrioles. Mol Biol Cell 16: 1095–1107.
      • 75. Zhu F, Lawo S, Bird A, Pinchev D, Ralph A, et al. (2008) The mammalian SPD-2 ortholog Cep192 regulates centrosome biogenesis. Curr Biol 18: 136–141.
      • 76. Mayor T, Stierhof YD, Tanaka K, Fry AM, Nigg EA (2000) The centrosomal protein C-Nap1 is required for cell cycle-regulated centrosome cohesion. J Cell Biol 151: 837–846.
      • 77. Uehara R, Nozawa RS, Tomioka A, Petry S, Vale RD, et al. (2009) The augmin complex plays a critical role in spindle microtubule generation for mitotic progression and cytokinesis in human cells. Proc Natl Acad Sci U S A 106: 6998–7003.
      • 78. Lawo S, Bashkurov M, Mullin M, Ferreria MG, Kittler R, et al. (2009) HAUS, the 8-subunit human Augmin complex, regulates centrosome and spindle integrity. Curr Biol 19: 816–826.
      • 79. Graser S, Stierhof YD, Nigg EA (2007) Cep68 and Cep215 (Cdk5rap2) are required for centrosome cohesion. J Cell Sci 120: 4321–4331.
      • 80. Moritz M, Braunfeld MB, Sedat JW, Alberts B, Agard DA (1995) Microtubule nucleation by gamma-tubulin-containing rings in the centrosome. Nature 378: 638–640.
      • 81. Zheng Y, Wong ML, Alberts B, Mitchison T (1995) Nucleation of microtubule assembly by a gamma-tubulin-containing ring complex. Nature 378: 578–583.
      • 82. Bloodgood RA, Salomonsky NL (1990) Calcium influx regulates antibody-induced glycoprotein movements within the Chlamydomonas flagellar membrane. J Cell Sci 96 (Pt 1): 27–33.
      • 83. Wang Q, Pan J, Snell WJ (2006) Intraflagellar transport particles participate directly in cilium-generated signaling in Chlamydomonas. Cell 125: 549–562.
      • 84. Yang P, Diener DR, Yang C, Kohno T, Pazour GJ, et al. (2006) Radial spoke proteins of Chlamydomonas flagella. J Cell Sci 119: 1165–1174.
      • 85. Shah AS, Ben-Shahar Y, Moninger TO, Kline JN, Welsh MJ (2009) Motile Cilia of Human Airway Epithelia Are Chemosensory. Science.
        • 86. Lorenzo IM, Liedtke W, Sanderson MJ, Valverde MA (2008) TRPV4 channel participates in receptor-operated calcium entry and ciliary beat frequency regulation in mouse airway epithelial cells. Proc Natl Acad Sci U S A 105: 12611–12616.
        • 87. Liu Q, Zhang Q, Pierce EA (2010) Photoreceptor sensory cilia and inherited retinal degeneration. Adv Exp Med Biol 664: 223–232.
        • 88. Zaghloul NA, Katsanis N (2009) Mechanistic insights into Bardet-Biedl syndrome, a model ciliopathy. J Clin Invest 119: 428–437.
        • 89. Chang B, Mandal MN, Chavali VR, Hawes NL, Khan NW, et al. (2008) Age-related retinal degeneration (arrd2) in a novel mouse model due to a nonsense mutation in the Mdm1 gene. Hum Mol Genet 17: 3929–3941.
        • 90. Liu Q, Tan G, Levenkova N, Li T, Pugh EN Jr, et al. (2007) The proteome of the mouse photoreceptor sensory cilium complex. Mol Cell Proteomics 6: 1299–1317.
        • 91. Okuda M, Horn HF, Tarapore P, Tokuyama Y, Smulian AG, et al. (2000) Nucleophosmin/B23 is a target of CDK2/cyclin E in centrosome duplication. Cell 103: 127–140.
        • 92. Matsumoto N, Yoneda-Kato N, Iguchi T, Kishimoto Y, Kyo T, et al. (2000) Elevated MLF1 expression correlates with malignant progression from myelodysplastic syndrome. Leukemia 14: 1757–1765.
        • 93. Yoneda-Kato N, Kato JY (2008) Shuttling imbalance of MLF1 results in p53 instability and increases susceptibility to oncogenic transformation. Mol Cell Biol 28: 422–434.
        • 94. Martin-Lanneree S, Lasbleiz C, Sanial M, Fouix S, Besse F, et al. (2006) Characterization of the Drosophila myeloid leukemia factor. Genes Cells 11: 1317–1335.
        • 95. Massinen S, Hokkanen ME, Matsson H, Tammimies K, Tapia-Paez I, et al. (2011) Increased expression of the dyslexia candidate gene DCDC2 affects length and signaling of primary cilia in neurons. PLOS One 6: e20580.
        • 96. Williams J, O’Donovan MC (2006) The genetics of developmental dyslexia. Eur J Hum Genet 14: 681–689.
        • 97. Paracchini S, Thomas A, Castro S, Lai C, Paramasivam M, et al. (2006) The chromosome 6p22 haplotype associated with dyslexia reduces the expression of KIAA0319, a novel gene involved in neuronal migration. Hum Mol Genet 15: 1659–1666.
        • 98. Poelmans G, Engelen JJ, Van Lent-Albrechts J, Smeets HJ, Schoenmakers E, et al. (2009) Identification of novel dyslexia candidate genes through the analysis of a chromosomal deletion. Am J Med Genet B Neuropsychiatr Genet 150B: 140–147.
        • 99. Taipale M, Kaminen N, Nopola-Hemmi J, Haltia T, Myllyluoma B, et al. (2003) A candidate gene for developmental dyslexia encodes a nuclear tetratricopeptide repeat domain protein dynamically regulated in brain. Proc Natl Acad Sci U S A 100: 11553–11558.
        • 100. Threlkeld SW, McClure MM, Bai J, Wang Y, LoTurco JJ, et al. (2007) Developmental disruptions and behavioral impairments in rats following in utero RNAi of Dyx1c1. Brain Res Bull 71: 508–514.
        • 101. Tammimies K, Vitezic M, Matsson H, Le Guyader S, Burglin TR, et al.. (2012) Molecular Networks of DYX1C1 Gene Show Connection to Neuronal Migration Genes and Cytoskeletal Proteins. Biol Psychiatry.
          • 102. Wang Y, Paramasivam M, Thomas A, Bai J, Kaminen-Ahola N, et al. (2006) DYX1C1 functions in neuronal migration in developing neocortex. Neuroscience 143: 515–522.
          • 103. Bettencourt-Dias M, Hildebrandt F, Pellman D, Woods G, Godinho SA (2011) Centrosomes and cilia in human disease. Trends Genet 27: 307–315.
          • 104. Tsou MF, Stearns T (2006) Controlling centrosome number: licenses and blocks. Curr Opin Cell Biol 18: 74–78.
          • 105. Rodrigues-Martins A, Bettencourt-Dias M, Riparbelli M, Ferreira C, Ferreira I, et al. (2007) DSAS-6 organizes a tube-like centriole precursor, and its absence suggests modularity in centriole assembly. Curr Biol 17: 1465–1472.
          • 106. Cizmecioglu O, Arnold M, Bahtz R, Settele F, Ehret L, et al. (2010) Cep152 acts as a scaffold for recruitment of Plk4 and CPAP to the centrosome. J Cell Biol 191: 731–739.
          • 107. Dzhindzhev NS, Yu QD, Weiskopf K, Tzolovsky G, Cunha-Ferreira I, et al. (2010) Asterless is a scaffold for the onset of centriole assembly. Nature 467: 714–718.
          • 108. Thornton GK, Woods CG (2009) Primary microcephaly: do all roads lead to Rome? Trends Genet 25: 501–510.
          • 109. Anitha A, Nakamura K, Yamada K, Iwayama Y, Toyota T, et al. (2009) Association studies and gene expression analyses of the DISC1-interacting molecules, pericentrin 2 (PCNT2) and DISC1-binding zinc finger protein (DBZ), with schizophrenia and with bipolar disorder. Am J Med Genet B Neuropsychiatr Genet 150B: 967–976.
          • 110. Gabel LA, Gibson CJ, Gruen JR, Loturco JJ (2009) Progress towards a cellular neurobiology of reading disability. Neurobiol Dis. 38: 173–180.