Coordinated genomic control of ciliogenesis and cell movement by RFX2
Abstract
The mechanisms linking systems-level programs of gene expression to discrete cell biological processes in vivo remain poorly understood. In this study, we have defined such a program for multi-ciliated epithelial cells (MCCs), a cell type critical for proper development and homeostasis of the airway, brain and reproductive tracts. Starting from genomic analysis of the cilia-associated transcription factor Rfx2, we used bioinformatics and in vivo cell biological approaches to gain insights into the molecular basis of cilia assembly and function. Moreover, we discovered a previously un-recognized role for an Rfx factor in cell movement, finding that Rfx2 cell-autonomously controls apical surface expansion in nascent MCCs. Thus, Rfx2 coordinates multiple, distinct gene expression programs in MCCs, regulating genes that control cell movement, ciliogenesis, and cilia function. As such, the work serves as a paradigm for understanding genomic control of cell biological processes that span from early cell morphogenetic events to terminally differentiated cellular functions.
https://doi.org/10.7554/eLife.01439.001eLife digest
Cells that have hundreds of tiny hair-like structures called cilia on their surface have important roles in our airways and also in the brain and reproductive system. By beating in a coordinated manner, the cilia cause fluid to flow in a particular direction. The development of these multiciliated cells is a complex process in which genes are expressed as proteins, with this gene expression being regulated by other proteins called transcription factors.
In invertebrates the development of the cilia is controlled by transcription factors from the RFX family, which also appear to be important for development of cilia in vertebrates. However, the details of this process—in particular, the identities of the genes that are involved and how their functions are related—are not well understood in vertebrates.
Chung et al. have sought to remedy this by analyzing the network of genes whose expression is controlled by the transcription factor Rfx2 in vertebrates. The results showed that the genes controlled by Rfx2 were involved in all aspects of cilia, including several genes that are known to be mutated in diseases caused by abnormal cilia. Chung et al. also identified genes that were not previously thought to be relevant to cilia.
As multiciliated cells are developing, but before they can generate cilia, they must first migrate from the bottom of the epithelium, the layer of tissue in which they function, to the top of this layer. Chung et al. found that Rfx2 was also involved in this process.
The approach taken by Chung et al.—which involved a combination of RNA sequence analysis, examination of Rfx2 binding sites on chromosomes, computational predictions of protein interactions and in vivo cellular imaging—could be used to perform similar systems-level analyses of other developmental and biological processes.
https://doi.org/10.7554/eLife.01439.002Introduction
A major goal of biology over the last several decades has been to understand the mechanisms that control differential gene expression. While recent advances in genomic technology have dramatically empowered these studies, we still know comparatively little about the mechanisms linking systems-level programs of gene expression to discrete cell biological processes in vivo. This gap in our understanding is important because organ and tissue function are ultimately executed by the specialized behaviors of individual cells (e.g., polarized secretion in excretory organs, coordinated contraction in muscle cells).
Recent studies of cilia and ciliated epithelial cells highlight the current disconnect between genomics and cell biology: it is clear that hundreds of proteins are required for cilia assembly and function (Gherman et al., 2006; Inglis et al., 2006), and moreover, assembly of new cilia/flagella clearly requires new transcription (Thomas et al., 2010). Nonetheless, only a handful of transcription factors have been identified that control cilia structure and function, and their associated gene regulatory networks remain largely undefined, especially in vertebrates. In C. elegans, the sole RFX family transcription factor daf-19 is the central regulator of ciliogenesis, and dozens of target genes are known to effect its action (Efimenko et al., 2005; Phirke et al., 2011; Burghoorn et al., 2012). The one Rfx factor in Drosophila is likewise well characterized (Laurencon et al., 2007; Newton et al., 2012). By contrast, multiple RFX family members are essential for ciliogenesis in vertebrates, but as yet, there has been no comprehensive genome-wide survey of Rfx-dependent gene expression as it relates to ciliogenesis (Bonnafe et al., 2004; Ashique et al., 2009; El Zein et al., 2009). This gap in our knowledge of the genomics of RFX factors is made the more important because these proteins also possess cilia-independent functions about which very little is yet known, including control neuronal and pancreatic development (e.g., Senti and Swoboda, 2008; Ait-Lounis et al., 2010; Pearl et al., 2011; Benadiba et al., 2012).
One vertebrate cell type in which Rfx factors are known to play particularly important roles is the multi-ciliated epithelial cell (MCC). These cells project dozens or hundreds of motile cilia from their apical surfaces, and the polarized beating of these cilia generates fluid flow that is essential for development and homeostasis in many organ systems (Figure 1). Such cells are central to the normal homeostasis of airway, brain and reproductive tracts (Worthington and Cathcart, 1963; Yeung et al., 1991; Lyons et al., 2006; Fahy and Dickey, 2010), and defective functioning of these cells is associated with pathologies ranging from chronic infection to hydrocephalus (Afzelius, 1976). Despite these cells’ importance and long history of study (Sharpey, 1830), the transcription factors that control MCC development and function are only now being elucidated (You et al., 2004; Stubbs et al., 2008; Marcet et al., 2011; Stubbs et al., 2012; Tan et al., 2013), and among these factors is Rfx2 (Chung et al., 2012).
In this study, we sought to combine systems biology approaches with in vivo cell biological studies in order to better define the genomic control of MCC development and function. We used high-throughput sequencing of Rfx2-regulated transcripts, systematic mapping of Rfx2 chromosomal binding sites, and bioinformatic exploration of functional protein interactions to guide our mechanistic cell biology experiments. This approach provided diverse new insights into the roles played by Rfx2 in ciliogenesis and polarized ciliary beating. Surprisingly, we also find that Rfx2 plays a central, but unanticipated, role in controlling cell movement in newborn MCCs, a process about which almost nothing is currently known. Overall, Rfx2 activates a complex program of gene expression that serves to coordinate several distinct features of MCCs, including cell migration, ciliogenesis, and cilia function, thus serving as a paradigm for genomic control of cell biological processes that span from early differentiation events to terminally differentiated cell functions.
Results
Systems-level analysis of Rfx2 function in a vertebrate mucociliary epithelium
To explore the genomic control of cell behavior in MCCs, we turned to the amphibian embryo, which has emerged as a powerful and rapidly assayable in vivo model for this cell type. Indeed, several foundational studies of specification, ciliogenesis, and planar polarization in MCCs have been performed in amphibians (Deblandre et al., 1999; Stubbs et al., 2006; Park et al., 2008; Mitchell et al., 2009; Marcet et al., 2011; Sirour et al., 2011; Werner and Mitchell, 2011; Stubbs et al., 2012). Importantly, these studies have consistently prefigured results in mammals (Marcet et al., 2011; Morimoto et al., 2010; Stubbs et al., 2012; Tsao et al., 2009; Vladar et al., 2012). We therefore exploited the Xenopus system to perform parallel RNA transcript sequencing (RNA-seq) and chromatin immunoprecipitation deep sequencing (ChIP-seq) for Rfx2, which we showed previously to be essential for the normal development of cilia in MCCs (Chung et al., 2012) (Figure 2).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.01439.005
-
Figure 2—source data 2
- https://doi.org/10.7554/eLife.01439.006
A detailed description of this approach is provided in the ‘Materials and methods’ section. Briefly, RNA-seq was performed on isolates of Xenopus mucociliary epithelium grown in organotypic culture, comparing control samples to knockdown using validated Rfx2 morpholino antisense oligonucleotides (Figure 2A; Figure 2—figure supplement 1) (Chung et al., 2012). We performed two biological replicates, each containing 100–200 individual tissue isolates; analysis of the correlation found excellent reproducibility between the two replicates (Figure 2—figure supplement 2). As we are most concerned here with results that will be relevant to human health, we confined our analysis of the RNA-seq data to Xenopus genes for which there were unambiguous human orthologues to facilitate cross-organism comparison. From this list of 11,644 genes, we identified 2750 genes that were differentially expressed following Rfx2 knockdown (>twofold; FDR<0.05) (Figure 2A, Figure 2—figure supplement 2).
We next used ChIP-seq to ask which of these differentially expressed genes were direct targets of control by Rfx2. We identified 3465 genes that were significantly bound by Rfx2 and, importantly, this set of genes overlapped significantly with the set of genes differentially expressed in Rfx2 morphants (Figure 2A, Figure 2—figure supplement 3) (p≤10−6, hypergeometric test). The intersection of these two experiments defined 911 genes (Figure 2—source data 1), corresponding to the directly regulated downstream target genes of Rfx2. We tested a subset of these directly bound, differentially expressed genes by in situ hybridization and RT-PCR analysis; 100% were confirmed as Rfx2-dependent genes (Figure 2—figure supplement 4). We then focused our further studies on this set of 911 direct targets.
Rfx2 controls assembly of diverse ciliogenic machinery
An unbiased analysis of enriched Gene Ontology (GO) terms suggested a key role for Rfx2 in the control of ciliary gene expression (Figure 2—source data 2), consistent with the known roles for Rfx2 in ciliogenesis (Chung et al., 2012). Moreover, direct examination of the 911 Rfx2 target genes revealed that components of essentially all known ciliary machinery are under the control of Rfx2 in MCCs (Figure 2B). Principal among these is the intraflagellar transport (IFT) system, the core mechanism for moving cargoes into and out of the cilium. IFT involves two separable complexes, IFT-A and IFT-B (Pedersen and Rosenbaum, 2008), and we identified the IFT–B complex components ift172 and ift88 as direct Rfx2 targets, consistent with data from flies, worms, and mice (Thomas et al., 2010). Moreover, we also identified several additional Rfx2-regulated components of both IFT-B and IFT-A complexes, as well as the IFT-A adaptor tulp3 (Mukhopadhyay et al., 2010) (Figure 2B; Supplementary file 1A). Additionally, we found that Rfx2 directly controls expression of genes encoding axonemal dynein subunits, components of the transition zone, and a component of the BBSome, and many of these genes are mutated in human ciliopathies (Figure 2B; Supplementary file 1A) (Sharma et al., 2008; Reiter et al., 2012).
Significantly, our analysis identified several ciliary systems that have not previously been associated with RFX factors. For example, Rfx2 controlled expression of tubulins, enzymes involved in tubulin post-translational modification, and microtubule binding proteins, such as map7 and spef1/clamp, which localize to the proximal and distal axonemes, respectively (Brooks and Wallingford, 2012) (Figure 2B; Supplementary file 1A). Rfx2 also controlled the expression of the planar cell polarity (PCP) effector gene fritz/wdpcp, which encodes a protein governing assembly of the septin diffusion barrier at the base of cilia (Kim et al., 2010). Rfx2 also controls many centriolar genes required for ciliogenesis, including cep164 (Graser et al., 2007). Notably, wdpcp and cep164 are both implicated in human ciliopathies (Kim et al., 2010; Chaki et al., 2012). This analysis thus provides a comprehensive view of Rfx-related ciliary gene expression in vertebrates.
Genomic analysis of Rfx2 identifies a novel regulator of intraflagellar transport
These data provide a foundation for understanding Rfx-mediated control of known cilia genes. However, a particularly striking feature of the Rfx2 target genes was the large number whose relationship to cilia (or lack thereof) is not known. Because transcriptionally co-regulated genes frequently share functions (Eisen et al., 1998; Marcotte et al., 1999), we next asked if our Rfx2 target gene set might provide a jumping-off point for exploring the molecular biology of ciliogenesis, with an aim towards implicating new proteins in this important process.
To focus our search, we used the direct Rfx2 target genes as a seed set to interrogate a probabilistic human gene network (HumanNet) that captures functionally linked genes based on observations in large-scale functional genomics and proteomics datasets (Lee et al., 2011). Using ‘guilt-by-association’ in this network, we identified ttc29 as a candidate for functional interaction with the IFT machinery (Figure 3A). We tested this prediction and found not only that Ttc29-GFP localized to ciliary axonemes, but also that knockdown elicited substantial defects in ciliogenesis (Figure 3B–D). To test the HumanNet prediction more directly, we used high-speed in vivo confocal imaging of IFT particle dynamics in Xenopus MCCs (Brooks and Wallingford, 2012). Strikingly, partial knockdown of Ttc29 resulted in a significant decrease in the mean rate of anterograde IFT (Figure 3E–G; Video 1), consistent with the bioinformatic linkage of ttc29 to components of the anterograde IFT-B complex (ift88 and ttc30a; Figure 3A). Retrograde IFT movement was not significantly affected by Ttc29 knockdown (Figure 3H). Thus, by combining our genomic dataset with functional gene networks, we have not only revealed new links between Rfx2 and ciliogenesis, but we have also identified a specific function for an uncharacterized protein in the regulation of anterograde IFT.
Genomic analysis of Rfx2 targets elucidates the mechanisms governing cilia beating
To further explore the mechanisms of Rfx2 function, we returned to the list of GO terms enriched in Rfx2 target genes. Of these, the most strongly enriched biological process GO term was ‘ciliary or flagellar motility’ (Figure 2—source data 2). A direct exploration of the Rfx2 targets identified a wide range of genes with known roles in cilia beating, including genes implicated in planar polarization of cilia (frizzled, prickle) (Mitchell et al., 2009; Vladar et al., 2012), tubulin polyglutamylases governing inner arm dynein activity (ttll7, ttll9, ttll13), dynein regulatory complex components, radial spoke components, and a variety of poorly-understood cilia beating genes, such as ropn1l, mns1, pacrg, and hydin (Lechtreck et al., 2008; Wilson et al., 2010; Fiedler et al., 2011; Zhou et al., 2012) (Figure 2B; Supplementary file 1B). Given this large array of cilia-beating genes controlled by Rfx2, we turned again to HumanNet, this time with the goal of providing new insights into the molecular control of cilia beating.
Unbiased clustering predicted that three known ciliary beating genes (dnai1, ropn1l, mns1) (Horvath et al., 2005; Fiedler et al., 2011; Zhou et al., 2012) should be functionally linked to ribc2, which encodes an uncharacterized vertebrate protein with similarity to Chlamydomonas protofilament ribbon proteins (Linck and Norrander, 2003) (Figure 4A). Using a GFP-fusion, we found that Ribc2 localized to axonemes, though knockdown did not substantially affect cilia length (Figure 4B–D). Rather, high-speed in vivo imaging of axonemes and analysis of fluid flow revealed a specific defect in ciliary beating following disruption of Ribc2 function (Figure 4E–H; Videos 3 and 4). Transmission electron microscopy revealed that loss of Ribc2 did not affect the integrity of outer doublet or central pair microtubules, but did disrupt their organization within the axoneme and the apparent number of dynein arms (Figure 4I,J, Figure 4—figure supplement 1).
Little is known about protofilament ribbon proteins in any system (Linck and Norrander, 2003), but guided by connections in HumanNet (Figure 4A), we found that Ribc2 is essential for the normal axonemal localization of another protofilament ribbon protein, Tekt2 (Figure 4K–M). Moreover, our approach also led us to assign a cilia beating function to the previously uncharacterized kinase Nme5 (Figure 4—figure supplement 2) and to find that proper localization of this kinase within axonemes is Ribc2-dependent (Figure 4N).
These data highlight the power of combining genomic data with functional gene network analysis: The approach identified new roles for Ribc2 and Nme5 in cilia beating and placed these proteins in a functional hierarchy with Tekt2. By contrast, pacrg is another direct target of Rfx2 that is required for cilia motility (Wilson et al., 2010) but is not linked to ribc2 in HumanNet, and Ribc2 knockdown did not affect the normal axonemal localization of Pacrg protein (Figure 4—figure supplement 3). Thus, our combination of genomics, bioinformatics, and in vivo cell biology provide an effective paradigm for understanding the links between system-level gene expression and discrete behaviors of individual cells.
Rfx2 is essential for the insertion of nascent MCCs into the mucociliary epithelium
Our data demonstrate that Rfx2 governs ciliogenesis and cilia beating via the expression of myriad genes, including many novel genes characterized here for the first time. Interestingly however, of the 911 Rfx2 target genes, only ∼20% are present in known cilia proteomes (Figure 5A). This finding suggests that many of the Rfx2 target genes are NOT involved in ciliogenesis or cilia function, a notion that is intriguing in light of recent reports of cilia-independent roles for ‘ciliary’ RFX genes (e.g., Senti and Swoboda, 2008; Ait-Lounis et al., 2010).
Analysis of GO terms associated with Rfx2 target genes suggested a potential role for Rfx2 in cell movement and cell morphogenesis (Figure 5A, Figure 2—source data 2), which is of interest because both Xenopus MCCs and their counterparts in mammalian airways arise from a population of p63-expressing basal precursor cells (Figure 1, red cell at right) (Evans and Moller, 1991; Drysdale and Elinson, 1992; Lu et al., 2001; Rock et al., 2010). Consistent with a role for Rfx2 in the apical movement of MCCs, we found that strong knockdowns consistently resulted in MCCs being positioned well below the apical surface of the epithelium (Figure 5B,C). This effect was not secondary to a ciliogenesis defect, as milder knockdowns with lower doses of MO did not inhibit MCC emergence but did suppress ciliogenesis (Chung et al., 2012). These results suggested a previously unrecognized role for an Rfx factor in cell movement.
A role for Rfx2 in apical movement is noteworthy, because this process is a conserved feature of newly born MCCs (Evans et al., 1989; Drysdale and Elinson, 1992; Lu et al., 2001; Rock et al., 2009), and it must involve apical movement, remodeling of junctions, and assembly of an apical cell surface to which basal bodies can dock prior to ciliogenesis. The mechanisms guiding these crucial aspects of MCC biology are largely unknown, so we used transgenic drivers to direct expression of fluorescent reporters specifically in MCCs and documented these cells’ behavior using 4D confocal imaging in vivo.
We identified two broad categories of MCC behavior associated with apical insertion. MCCs first engaged in a probing of neighboring cell–cell boundaries; subsequently, this behavior ceased and MCCs smoothly expanded their apical surfaces (Figure 5D; Videos 5 and 6). Time-lapse imaging in Rfx2 morphants revealed that the initial probing behavior was qualitatively normal, while the latter phase of apical surface expansion failed completely (Figure 5E; Video 7).
Rfx2 acts cell-autonomously to control apical surface expansion in nascent MCCs
Rfx2 is expressed strongly in MCCs (Chung et al., 2012) but RT-PCR also detected expression of Rfx2 in the neighboring goblet cells (not shown). Thus, Rfx2 might act in MCCs to control their cell movement and insertion into the epithelium or it may control behaviors in neighboring superficial cells that facilitate the process—or both. To distinguish between these possibilities, we generated mosaics in which Rfx2 was disrupted specifically in MCCs or specifically in the overlying goblet cells. Knockdown of Rfx2 in the MCCs led to a robust failure of MCC apical surface expansion while knockdown in overlying superficial cells had no effect (Figure 6). Thus our data suggest that Rfx2 cell-autonomously controls the process of MCC apical surface expansion.
Rfx2 controls insertion of nascent MCCs via dystroglycan and slit2
Insertion of MCCs into the mucociliary epithelium represents a novel cell behavior, the molecular mechanisms of which remain almost entirely obscure. Given the specific role for Rfx2 in this process, we leveraged our genomic dataset to gain insights. Consistent with the observed insertion defects, we found that Rfx2 directly regulates the extracellular matrix component dystroglycan (dag1), which is a known regulator of MCC insertion (Sirour et al., 2011). Interestingly, Dag1 is known to functionally interact with the Slit/Robo signaling system to control cellular morphogenesis in diverse settings (Medioni et al., 2008; Wright et al., 2012), so it was notable that Rfx2 directly controlled transcription of both the slit2 ligand and the intracellular effector srgap2 (Supplementary file 1C). Accordingly, we found that Slit2 knockdown disrupted insertion of nascent MCCs into the mucociliary epithelium and analysis of mosaic embryos revealed that Slit2 acts cell-autonomously in MCC insertion, as was the case for Rfx2 (Figure 7). Together, these data identify rfx2, dag1, and slit2 as a preliminary molecular framework for epithelial insertion of nascent MCCs and reveal a tight genomic coordination of cell movement, ciliogenesis, and cilia motility.
Discussion
During development, many cell types must first execute a specific set of transient behaviors before finally differentiating into a final, functional form. For example, many cells must first undergo transitions from mesenchyme to epithelium or vice-versa, must engage in migration, or must radically change their shape. Then, as they terminally differentiate, most cells will assemble specialized cellular machinery in order to execute their cell-type specific functions. Understanding how disparate cell types execute specialized cell behaviors is therefore a key challenge in developmental biology and will require us not only to delineate gene regulatory networks, but also to associate particular genes in the network with particular cellular machinery. Meeting this challenge will require the application of genomic technologies specifically to cell biological questions in vivo.
In this study, we explored the role of the Rfx2 transcription factor in MCCs, an essential cell type for homeostasis of the brain, airway, and reproductive tracts. We performed a systems-level analysis of Rfx2-dependent gene expression, and we used functional gene networks to guide exploration of the resulting datasets. By coupling this approach to validation by in vivo imaging, this study provides important insights into the control of MCC gene expression (Figure 8, left) and into the molecular control of discrete MCC cell behaviors (Figure 8, right). Moreover, the work dramatically expands our understanding of the vertebrate RFX family transcription factors in ciliogenesis, in cilia function, and surprisingly, in cell movement.
Rfx2 and the genetic network of multiciliated cell development and function
In addition to Rfx2, recent work has identified keys roles for the transcription factors Mcidas and Foxj1 in the development of MCCs (You et al., 2004; Stubbs et al., 2008, 2012). A key challenge going forward will be to understand how Rfx2 fits into this larger gene regulatory network, so we compared our set of direct Rfx2 targets with the downstream transcriptome of these other factors (Figure 8, left).
Mcidas has been implicated in the earliest stages of MCC specification (Stubbs et al., 2012), and accordingly, both Rfx2 and FoxJ1 are downstream of Mcidas, while Mcidas is not present in the target gene sets of the other two factors. Interestingly, our comparison suggests that many ciliary machines require the combined action of all or a subset of these three factors. For example, ift80 is downstream of all three factors, while Ift46 and ift57 are present only in the FoxJ1 gene set and ift140, ift172, ift81, and ift88 are targets only of Rfx2. Such combined action by Rfx factors and FoxJ1 is consistent with recent reports in mammalian airway epithelium (Didon et al., 2013) and in Drosophila (Newton et al., 2012). Nonetheless, our data also suggest that Rfx2 also plays essential roles in ciliogenesis that are independent of FoxJ1 or Mcidas, as many IFT genes and several ciliopathy-associated genes (cep164, rpgripl1 and cc2d2a) are under the sole control of Rfx2 (Figure 8, left).
Moreover, we have discovered a role here for Rfx2 in apical surface assembly of nascent MCCs, and Foxj1 has no known role in this process. Accordingly, the MCC insertion genes dag1 and slit2 were present only in the Rfx2 target gene list and were not identified as downstream targets of FoxJ1 or Mcidas. Interestingly however, Rab11 is also required for MCC insertion (Kim et al., 2012), but was not present among our Rfx2 targets. Future studies will be needed to fully understand the gene regulatory network of MCCs, including how other essential transcription factors (e.g., Myb [Tan et al., 2013]), collaborate with the factors discussed here. Finally, it will now be of great interest to ask how this emerging gene regulatory network governing cilia structure and function in multiciliated cells compares with that governing the structure and function of primary, non-motile cilia in other cell types.
From Rfx2-dependent gene expression to cell biological mechanism
Our analysis revealed that Rfx2 controls cilia structure and function by modulating expression of dozens of known genes, and this dataset serves as an important complement to genomic studies of RFX factors in invertebrates (Efimenko et al., 2005; Laurencon et al., 2007; Newton et al., 2012; Phirke et al., 2011; Burghoorn et al., 2012). However, a large proportion of genes, even in model animals, remain uncharacterized, and our Rfx2 target gene list contains dozens of genes with no known function. To overcome this hurdle, we used a functional gene network to identify important functional interactions among novel genes and known genes in the Rfx2 target gene set and we used in vivo imaging to test these interactions (Figure 8, right). This approach led us to discover an essential role for Ttc29 in the control of anterograde IFT and ciliogenesis and to begin delineating a functional hierarchy for the novel proteins Ribc2 and Nme5 in ciliary beating.
Moreover, guided by our exploration of gene networks, we found that the Rfx2 targets Tekt2 and Nme5 localize in novel proximodistally-restricted patterns along the axoneme, a result that is of interest because precise proximodistal positioning of specific dynein arm proteins along motile axonemes is central to normal cilia beating. Indeed, though defective proximodistal pattern in axonemes is associated with human primary ciliary dyskinesia (PCD) (Fliegauf et al., 2005; Panizzi et al., 2012), we know essentially nothing of how these patterns are controlled. Thus, by exploiting our genomic data to guide cell biological inquiry, this study has provided a foundation for future exploration of this problem.
Cilia-independent roles for Rfx2 and apical surface assembly in nascent MCCs
The role of RFX factors in ciliogenesis has been a key focus of recent research, so it is notable that this function is not thought to be ancestral for RFX proteins (Piasecki et al., 2010). Indeed, the yeast RFX orthologue, Sak1, controls cell cycle exit (Wu and McLeod, 1995) and even Rfx factors commonly associated with ciliogenesis (e.g., Daf-19, Rfx3) perform cilia-independent functions (Senti and Swoboda, 2008; Ait-Lounis et al., 2010). Accordingly, we found that Rfx2 also controls the assembly of the apical surface of nascent MCCs as they insert into the epithelium. Rfx2 effects this process via the known regulator dag1 and through slit2, which we show here to be essential for MCC apical surface expansion. These results are of particular interest both because Dag1 and Slit/Robo signaling also collaborate during axon growth cone guidance (Andrews et al., 2006; Wright et al., 2012) and because Rfx3 mutant mice display axon guidance defects (Benadiba et al., 2012). These data may also shed light on the mechanism of Rfx action in synapse morphogenesis in C. elegans (Senti and Swoboda, 2008), as the neuron projection morphogenesis GO term was enriched in our Rfx2 target gene set, reflecting that mediators of synapse morphogenesis, such as PTPδ and Netrin3, were present in the target set (Supplementary file 1C) (Takahashi et al., 2012). Finally, both the genomic studies and our in vivo imaging of MCC insertion provide an important foundation for future work, because mammalian MCCs—like the Xenopus MCCs studied here—arise from basally-positioned precursor cells (Evans and Moller, 1991; Rock et al., 2010).
In conclusion, this study has revealed a central role for Rfx2 at the nexus of cell movement, ciliogenesis, and cilia motility (Figure 8). The work sheds new light on Rfx2 protein functions specifically and also on the transcriptional control of cell behavior and organelle biogenesis generally. Our combined approach of systems biology, computational biology, and in vivo cell biology provides a generalizable paradigm for exploiting genomic data to advance our understanding of cell biological processes.
Materials and methods
Morpholino and RNA injection
Request a detailed protocolCapped mRNA was synthesized using mMESSAGE mMACHINE (Ambion, Austin, TX). mRNA and anti-sense morpholino were injected into ventral blastomeres at the 4-cell stage to target the epidermis (Moody, 1987). Embryos were incubated until appropriate stages and were fixed in MEMFA (Davidson and Wallingford, 2005). The embryos were embedded in 2% agarose for thick (250–300 micron) sections, which were cut with a Vibratome series 1000 (Davidson and Wallingford, 2005). Morpholino sequence and the working concentration:
Rfx2 morpholino: AATTCTGCATACTGGTTTCTCCGTC, 12 ng
Ttc29 morpholino: GTGCACTCATTCTCTTCAAGTTTGC, 40 ng
Ribc2 morpholino: CGATAGGCAGATCCAGTCGGTACAT, 21 ng
Slit2 morpholino: TTCAGGTCTCTGGGAAAACAGGAAC, 10 ng
Draft genome and annotation of Xenopus laevis
Request a detailed protocolWe used draft version (version 6, JGIv6) of the Xenopus laevis genome for analyzing genomic and ChIP-seq datasets, obtained from the International Xenopus laevis genome project consortium and now available from the XenBase FTP site (ftp://ftp.xenbase.org). We selected scaffolds longer than 10,000 bases for further analysis, for a total of 8426 scaffolds used in this analysis. Because gene models for X. laevis are not yet finalized, we employed interim gene models for the analysis described here, using the released transcriptome-derived gene models (‘Oktoberfest’ version), also provided by the International Xenopus laevis genome project. Sequences and detailed descriptions of the gene model construction pipeline are available at the project website (http://www.marcottelab.org/index.php/Xenopus_Genome_Project).
Differential expression following Rfx2 knockdown as measured by RNA-seq
Request a detailed protocolTotal RNA was collected from 100 animal caps each at stage 20 from control embryos and from Rfx2 morphants. After poly-(A)-capture, we prepared sequencing libraries using the standard manufacturer’s non-strand specific Ilumina RNA-seq protocol, and sequenced paired-end 2 × 50 bp reads by Ilumina HiSeq 2000. Reads were mapped to the longest transcripts for each of the ‘Oktoberfest’ gene models using the bowtie mapper (version 0.12.7; ‘-a -v 2’ options were applied) (Langmead et al., 2009). We then estimated normalized gene expression values for each sample by calculating RPKMs (Reads Per Kilobase per Million mapped reads). It should be noted that one of the two Rfx2 MO replicates (M2) had significantly fewer reads (30 M read pairs) than the other (72 M read pairs). However, measured expression levels (RPKMs) correlated well between the replicates (Figure 2—figure supplement 2), so we retained both for subsequent analyses. Using the edgeR package (Robinson et al., 2010), we identified significant changes in RNA abundances between control and Rfx2 knockdowns, requiring greater than twofold abundance changes and a false discovery rate (FDR) less than 5%. Out of 11,644 genes tested, we identified 2750 genes significantly differentially expressed following Rfx2 knockdown.
Rfx2 chromatin-immunoprecipitation and sequencing (ChIP-seq)
Request a detailed protocolChIP-seq was performed as described previously (Kim et al., 2011) to identify direct chromosomal binding sites of Rfx2. Briefly, ChIP-seq assays were performed with stage 20 X. laevis embryos, injected with in vitro transcribed mRNA coding for either GFP-Rfx2 or GFP alone. Samples were crosslinked with 1% formaldehyde for 1 hr and the reaction stopped by adding glycine (to 125 mM) for 10 min. The embryos were rinsed with PBS and resuspended in lysis buffer with protease inhibitor cocktail (Roche). Chromatin was sonicated to an average size of 200–600 base pair using a Branson 450 Sonifier, then immunoprecipitated using protein G magnetic beads (Invitrogen) coupled to 5 μg α-GFP antibody (ab290) at 4°C overnight. Magnetic beads were washed, the bound chromatin eluted, and crosslinks reversed. ChIP DNA was extracted with phenol-chloroform and purified with a QIAquick PCR Purification Kit (Qiagen). ChIP-seq libraries were prepared according to the standard manufacturer’s Illumina sequencing protocol and sequenced by Illumina HiSeq.
Reads were mapped to the Xenopus laevis draft genome (version 6) using bowtie (‘-m 1 -n 2’ options were applied) and peaks identified using MACS (version 1.4.2) (Zhang et al., 2008) with default options. Out of 29,448 peaks identified, we selected 6646 peaks for further study that exhibited either a FDR <5% or a fold-enrichment > 20. Genes were associated with significant ChIP-seq peaks based on proximity in the draft genome, requiring genes (specifically, the longest transcript for each gene model as mapped to the draft genome sequence) to lie within 10,000 bases from an identified peak. Out of the 6646 significant ChIP-seq peaks, 3465 peaks could be assigned to nearby genes and 911 of those putative target genes also showed significantly different gene expression after Rfx2 knockdown. We focused on these differentially expressed, directly bound genes for subsequent analyses (Figure 2—source data 1).
Analysis of RFX2 targets using the HumanNet gene network
Request a detailed protocolTo better understand the molecular mechanisms of genes regulated by RFX2, we analyzed the clustering of RFX2 target genes using a human functional gene network (‘HumanNet’) (Lee et al., 2011). To increase potential coverage, we considered X. laevis genes that were either marked by RFX2 binding sites in the ChIP-seq data or that were significantly differentially expressed following RFX2. For the human orthologs of these genes (as defined by the International Xenopus laevis genome project using phylogenetic analyses of gene models during the course of annotating X. laevis genes), we extracted HumanNet gene–gene linkage information and associated confidence scores (log likelihood scores; LLS). The resulting network contained 4609 genes with 52,714 functional linkages, and served as the basis for later analyses.
Due to an ancestral genome duplication along the X. laevis lineage, many human genes have (typically) two X. laevis orthologs, generally referred to as ‘homeologs’ or ‘alloalleles’. For the purposes of calculating gene networks among RFX2 target genes, we transferred evidence from either of the alloalleles to the single orthologous human gene. For example, in the case of having one human gene orthologous to two frog genes, if one of the frog homologs was identified as a direct target of RFX2 and the other was not, we considered the human gene as a direct target of RFX2 for the purposes of reconstructing the network. To identify functional modules regulated by RFX2 in unbiased way, we clustered this network using the clusterONE algorithm (Nepusz et al., 2012) available in Cytoscape (version 2.8.3) (Shannon et al., 2003; Cline et al., 2007), considering linkage confidence scores during the clustering. All network information is available at the following URL: http://www.marcottelab.org/index.php/ChungKwon2013_RFX2
Comparison of Rfx2 targets to the cilia proteome
Request a detailed protocolFor the comparison of Rfx2 targets to ciliary proteins in Figure 5A, we used a compiled list of ciliary proteins drawn from several studies using a combination of proteomics and comparative genomics (Gherman et al., 2006). The protein set was downloaded from this website, http://v3.ciliaproteome.org/cgi-bin/protein_browser.php, then converted to EnsEMBL gene IDs using BioMart (version 63).
In situ hybridization
Request a detailed protocolIn situ hybridization was performed as described previously (Sive et al., 2000). Bright field and low magnification fluorescence images were captured on a fluorescent stereomicroscope, Leica MZ16FA.
Immunohistochemistry
Request a detailed protocolEmbryos were fixed in MEMFA for 1 hr followed by washing in PTW (PBS+0.1% Tween 20) for 30 min (3 × 10 min). Embryos were then blocked in fetal bovine serum (FBS) solution (TBS containing 10% FBS and 5% DMSO) for 1 hr at room temperature. Embryos were then incubated with the following primary antibodies at 4°C overnight: monoclonal anti-α-tubulin antibody (1:500 dilutions, clone DM1A, Sigma), mouse anti-acetylated-α-tubulin (1:500, clone 6-11B-1, Sigma), chicken anti-GFP antibody (1:500, ab13970), and rabbit anti-RFP (1:500, ab62341). After primary antibody incubation, all samples were washed with TBST (TBS+0.1% Triton X-100) for 5 hr (1 × 5 hr). Primary antibodies were detected with Alexa Fluor 488 goat anti-mouse IgG (1:500, Molecular Probes), Alexa-555 goat anti-rabbit IgG (1:500, Molecular Probes), and Alexa-488 goat anti-chicken IgG (1:500, Molecular Probes). After secondary antibody incubation, all samples were washed with TBST for 5 hr.
Embryos prepared for confocal imaging as described (Wallingford, 2010). Images were obtained using Zeiss LSM5 Pascal and Zeiss LSM700 confocal microscope. Cilia lengths were measured with Fiji software. Images used throughout this paper have been enhanced using the Unsharp Mask filter in Adobe Photoshop.
Confocal imaging of live embryos
Request a detailed protocolHigh-speed confocal imaging was performed by time-lapse collection of single optical section at a frame rate of 370fps using a Zeiss LSM 5LIVE microscope. Images were collected from living embryos expressing membrane-GFP driven by MCC-specific promoter and from embryos injected with Ribc2 morpholino.
For filming the MCC intercalation, living embryos of either control or Rfx2 morpholino-injected animals were put on a round cover glass in custom machined dishes (Kieserman et al., 2010). Embryos were gently pushed down by a small piece of cover glass. Images were collected every 3–5 min and were then processed into a time-lapse movie using Fiji software.
For IFT imaging, embryos expressing GFP-IFT20 alone or with Ttc29 MO were mounted flank down in 0.8% LMP agarose. Single confocal slices were collected at ∼2 fps using an LSM 5LIVE confocal microscope, as previously described in Brooks and Wallingford (2012).
Cell-autonomy/transplantation assays
Request a detailed protocolEmbryos were injected ventrally at 4-cell stage with mRNA encoding either membrane-GFP or membrane-RFP. Rfx2 morpholino was injected together with membrane-RFP. At stage 10, a fine hair was used to peel off the outer layer from a region of the ectoderm of a donor embryo. This outer layer peel was then transferred onto a host embryo after removing a similar patch of outer cells. To help the healing process, a small piece of glass coverslip with clay feet was used to press down embryos. Transplantations were performed in Danilchick’s Solution for Amy (DFA) + 0.1% BSA. After healing, embryos were then transferred back to 1/3 MMR (Stubbs et al., 2006).
Data availability
-
Coordinated genomic control of ciliogenesis and cell movementID GSE50593. Publicly available at GEO (http://www.ncbi.nlm.nih.gov/geo/).
-
Xenopus laevis genome (JGI version 6.0)Publicly available at XenBase (http://www.xenbase.org/), the Xenopus model organism database.
-
Xenopus laevis gene model (Oktoberfest version)Publicly available at XenBase (http://www.xenbase.org/), the Xenopus model organism database.
-
HumanNet - probabilistic functional network of humanPublicly available on the project website (http://www.functionalnet.org/humannet/).
References
-
The transcription factor RFX3 directs nodal cilium development and left-right asymmetry specificationMolecular and Cellular Biology 24:4417–4427.https://doi.org/10.1128/MCB.24.10.4417-4427.2004
-
Control of vertebrate intraflagellar transport by the planar cell polarity effector FuzThe Journal of Biology 198:37–45.https://doi.org/10.1083/jcb.201204072
-
RFX2 is broadly required for ciliogenesis during vertebrate developmentDevelopmental Biology 363:155–165.https://doi.org/10.1016/j.ydbio.2011.12.029
-
Integration of biological networks and gene expression data using CytoscapeNature Protocols 2:2366–2382.https://doi.org/10.1038/nprot.2007.324
-
BookVisualizing cell biology and tissue movements during morphogenesis in the frog embryoIn: Yuste R, Konnerth A, editors. Imaging in neuroscience and developement. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press. pp. 125–136.
-
A two-step mechanism generates the spacing pattern of the ciliated cells in the skin of Xenopus embryosDevelopment 126:4715–4728.
-
Cell migration and induction in the development of the surface ectodermal patttern of the Xenopus laevis tadpoleDevelopment Growth & Differentiation 34:51–59.https://doi.org/10.1111/j.1440-169X.1992.00051.x
-
Cluster analysis and display of genome-wide expression patternsProceedings of the National Academy of Sciences of the United States of America 95:14863–14868.https://doi.org/10.1073/pnas.95.25.14863
-
The role of basal cells in attachment of columnar cells to the basal lamina of the tracheaAmerican Journal of Respiratory Cell and Molecular Biology 1:463–469.https://doi.org/10.1165/ajrcmb/1.6.463
-
Biology of airway basal cellsExperimental Lung Research 17:513–531.https://doi.org/10.3109/01902149109062862
-
Airway mucus function and dysfunctionThe New England Journal of Medicine 363:2233–2247.https://doi.org/10.1056/NEJMra0910061
-
Loss of ASP but not ROPN1 reduces mammalian ciliary motilityCytoskeleton (Hoboken) 69:22–32.https://doi.org/10.1002/cm.20539
-
Mislocalization of DNAH5 and DNAH9 in respiratory cells from patients with primary ciliary dyskinesiaAmerican Journal of Respiratory and Critical Care Medicine 171:1343–1349.https://doi.org/10.1164/rccm.200411-1583OC
-
Cep164, a novel centriole appendage protein required for primary cilium formationThe Journal of Biology 179:321–330.https://doi.org/10.1083/jcb.200707181
-
Identification and analysis of axonemal dynein light chain 1 in primary ciliary dyskinesia patientsAmerican Journal of Respiratory Cell and Molecular Biology 33:41–47.https://doi.org/10.1165/rcmb.2004-0335OC
-
High-magnification in vivo imaging of Xenopus embryos for cell and developmental biologyCold Spring Harbor Protocols 2010:pdb.prot5427.https://doi.org/10.1101/pdb.prot5427
-
Rab11 regulates planar polarity and migratory behavior of multiciliated cells in Xenopus embryonic epidermisDevelopmental Dynamics 241:1385–1395.https://doi.org/10.1002/dvdy.23826
-
Mutations in Hydin impair ciliary motility in miceThe Journal of Biology 180:633–643.https://doi.org/10.1083/jcb.200710162
-
Xenopus p63 expression in early ectoderm and neurectodermMechanisms of Development 102:275–278.https://doi.org/10.1016/S0925-4773(01)00315-X
-
The reproductive significance of human Fallopian tube ciliaHuman Reproduction Update 12:363–372.https://doi.org/10.1093/humupd/dml012
-
Genetic control of cell morphogenesis during Drosophila melanogaster cardiac tube formationThe Journal of Biology 182:249–261.https://doi.org/10.1083/jcb.200801100
-
Fates of the blastomeres of the 32-cell-stage Xenopus embryoDevelopmental Biology 122:300–319.https://doi.org/10.1016/0012-1606(87)90296-X
-
Functional analysis of Rfx6 and mutant variants associated with neonatal diabetesDevelopmental Biology 351:135–145.https://doi.org/10.1016/j.ydbio.2010.12.043
-
Intraflagellar transport (IFT) role in ciliary assembly, resorption and signallingCurrent Topics in Developmental Biology 85:23–61.https://doi.org/10.1016/S0070-2153(08)00802-8
-
Regulatory Factor X (RFX)-mediated transcriptional rewiring of ciliary genes in animalsProceedings of the National Academy of Sciences of the United States of America 107:12969–12974.https://doi.org/10.1073/pnas.0914241107
-
Basal cells as stem cells of the mouse trachea and human airway epitheliumProceedings of the National Academy of Sciences of the United States of America 106:12771–12775.https://doi.org/10.1073/pnas.0906850106
-
Airway basal stem cells: a perspective on their roles in epithelial homeostasis and remodelingDisease Models & Mechanisms 3:545–556.https://doi.org/10.1242/dmm.006031
-
Distinct isoforms of the RFX transcription factor DAF-19 regulate ciliogenesis and maintenance of synaptic activityMolecular Biology of the Cell 19:5517–5528.https://doi.org/10.1091/mbc.E08-04-0416
-
Ciliary dysfunction in developmental abnormalities and diseasesCurrent Topics in Developmental Biology 85:371–427.https://doi.org/10.1016/S0070-2153(08)00813-2
-
On a peculiar motion excited in fluids by the surface of certain animalsEdinburgh Medical and Surgical Journal 34:113–122.
-
Dystroglycan is involved in skin morphogenesis downstream of the Notch signaling pathwayMolecular Biology of the Cell 22:2957–2969.https://doi.org/10.1091/mbc.E11-01-0074
-
BookEarly development of Xenopus laevis: a laboratory manualCold Spring Harbor, NY: Cold Spring Harbor Press.
-
Selective control of inhibitory synapse development by Slitrk3-PTPdelta trans-synaptic interaction389–398, Nature Neuroscience, 15, S1-2, 10.1038/nn.3040.
-
Transcriptional control of genes involved in ciliogenesis: a first step in making ciliaBiology of the cell/under the auspices of the European Cell Biology Organization 102:499–513.https://doi.org/10.1042/BC20100035
-
Microtubules enable the planar cell polarity of airway ciliaCurrent Biology 22:2203–2212.https://doi.org/10.1016/j.cub.2012.09.046
-
Preparation of fixed Xenopus embryos for confocal imagingCold Spring Harbor Protocols 2010:pdb.prot5426.https://doi.org/10.1101/pdb.prot5426
-
The sak1+ gene of Schizosaccharomyces pombe encodes an RFX family DNA-binding protein that positively regulates cyclic AMP-dependent protein kinase-mediated exit from the mitotic cell cycleMolecular Cellular Biology 15:1479–1488.
-
Organization of tubules in the human caput epididymidis and the ultrastructure of their epitheliaThe American Journal of Anatomy 191:261–279.https://doi.org/10.1002/aja.1001910306
-
Role of f-box factor foxj1 in differentiation of ciliated airway epithelial cellsAmerican journal of physiology. Lung cellular and molecular physiology 286:L650–L657.https://doi.org/10.1152/ajplung.00170.2003
-
Model-based analysis of ChIP-seq (MACS)Genome Biology 9:R137.https://doi.org/10.1186/gb-2008-9-9-r137
Article and author information
Author details
Funding
National Institute of General Medical Sciences
- Julie C Baker
- Edward M Marcotte
- John B Wallingford
Howard Hughes Medical Institute
- John B Wallingford
Cancer Prevention Research Institute of Texas
- Edward M Marcotte
National Science Foundation
- Edward M Marcotte
The US Army (58343-MA)
- Edward M Marcotte
The Welch Foundation (F1515)
- Edward M Marcotte
National Heart, Lung, and Blood Institute
- John B Wallingford
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank The International Xenopus laevis Genome Project Consortium, especially the Harland, Rokhsar, and Taira labs, as well as all data contributors to the consortium. We also thank Scott Hunicke-Smith and The Genome Sequencing and Analysis Facility (GSAF) at UT Austin for RNA-seq and the Texas Advanced Computing Center (TACC) for supercomputer use. JBW is an Early Career Scientist of the Howard Hughes Medical Institute.
Ethics
Animal experimentation: All animal experiments were approved by the IACUC of the University of Texas at Austin, Protocol# AUP-2012-00156.
Copyright
© 2014, Chung et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,792
- views
-
- 414
- downloads
-
- 118
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Neuroscience
Chronic high-fat feeding triggers metabolic dysfunction including obesity, insulin resistance, and diabetes. How high-fat intake first triggers these pathophysiological states remains unknown. Here, we identify an acute microglial metabolic response that rapidly translates intake of high-fat diet (HFD) to a surprisingly beneficial effect on metabolism and spatial/learning memory. High-fat intake rapidly increases palmitate levels in cerebrospinal fluid and triggers a wave of microglial metabolic activation characterized by mitochondrial membrane activation and fission as well as metabolic skewing toward aerobic glycolysis. These effects are detectable throughout the brain and can be detected within as little as 12 hr of HFD exposure. In vivo, microglial ablation and conditional DRP1 deletion show that the microglial metabolic response is necessary for the acute effects of HFD. 13C-tracing experiments reveal that in addition to processing via β-oxidation, microglia shunt a substantial fraction of palmitate toward anaplerosis and re-release of bioenergetic carbons into the extracellular milieu in the form of lactate, glutamate, succinate, and intriguingly, the neuroprotective metabolite itaconate. Together, these data identify microglia as a critical nutrient regulatory node in the brain, metabolizing away harmful fatty acids and liberating the same carbons as alternate bioenergetic and protective substrates for surrounding cells. The data identify a surprisingly beneficial effect of short-term HFD on learning and memory.
-
- Cell Biology
- Chromosomes and Gene Expression
During oncogene-induced senescence there are striking changes in the organisation of heterochromatin in the nucleus. This is accompanied by activation of a pro-inflammatory gene expression programme - the senescence associated secretory phenotype (SASP) - driven by transcription factors such as NF-κB. The relationship between heterochromatin re-organisation and the SASP has been unclear. Here we show that TPR, a protein of the nuclear pore complex basket required for heterochromatin re-organisation during senescence, is also required for the very early activation of NF-κB signalling during the stress-response phase of oncogene-induced senescence. This is prior to activation of the SASP and occurs without affecting NF-κB nuclear import. We show that TPR is required for the activation of innate immune signalling at these early stages of senescence and we link this to the formation of heterochromatin-enriched cytoplasmic chromatin fragments thought to bleb off from the nuclear periphery. We show that HMGA1 is also required for cytoplasmic chromatin fragment formation. Together these data suggest that re-organisation of heterochromatin is involved in altered structural integrity of the nuclear periphery during senescence, and that this can lead to activation of cytoplasmic nucleic acid sensing, NF-κB signalling, and activation of the SASP.