RNA Function
RNA Function
RNA Function
*Correspondence: [email protected]
DOI 10.1016/j.cell.2007.08.003
Of the genes expressed during this developmental signals due to cross-hybridization to mRNAs with similar
window, a surprising 71% were found to encode mRNAs sequences is highly unlikely (Figure S2).
exhibiting clear subcellular distribution patterns. The fre- Following high-throughput FISH, samples were
quency and variety of localization events suggests that vir- mounted and analyzed using epifluorescence micros-
tually all aspects of cellular function are impacted by RNA copy. For each expressed gene, representative low and
trafficking pathways. We conclude that mRNA localization high magnification images were captured at key develop-
is a major mechanism for controlling cellular architecture mental stages and incorporated within a relational data-
and function. base. The first 4.5 hr of Drosophila development, spanning
embryonic stages 1–9, was chosen for analysis, as this in-
RESULTS terval is manageable in terms of data annotation and en-
compasses major developmental landmarks such as the
Method Development, Screening, midblastula transition (MBT), gastrulation and the specifi-
and Localization Database cation of many cell types. The MBT is the period during
To circumvent the deficiencies of existing in situ hybridiza- which developmental regulation switches from control
tion protocols, considerable effort was made to develop by maternally synthesized gene products to control by zy-
a procedure with optimal subcellular resolution, sensitiv- gotically transcribed genes (Tadros et al., 2007b). Impor-
ity, consistency, throughput and economy. Typical results tantly, our FISH method enables the unambiguous distinc-
obtained with the resulting method (Lécuyer et al., 2007), tion between maternal and zygotic mRNA populations.
versus traditional alkaline phosphatase-based probe de- Maternally provided transcripts are generally cytoplasmic,
tection, are illustrated in Figures 1A and S1 in the Supple- exhibit FISH signal intensities above background in stage
mental Data available with this article online. Reassuringly, 1 embryos and decrease in intensity during later stages.
control analyses using probes with increasing sequence Zygotic mRNAs, on the other hand, are always first
divergence indicate that the occurrence of false positive detected in nascent transcript foci within subsets of
transcripts is far outnumbered by those that are posteri- The posterior group of mRNAs could be clearly subdi-
orly localized (5 versus 198; Table 1). For both of these vided into three categories: (1) transcripts that localize to
categories, the sensitivity and resolution of detection early pole plasm (Figure 2C), a specialized region of cyto-
enabled the identification of distinct subgroups of patterns plasm that directs the formation of germline pole cells; (2)
(Figure 2). For example, among the anteriorly localized transcripts that reside in the pole plasm and then localize
mRNAs, which include bcd, CycB, lok, milt and asp, only further into distinctive rings around pole cell nuclei
bcd exhibits tight anterior localization (Figure 2A), con- (Figure 2D); and (3) transcripts that only begin to localize
sistent with its extensively characterized function as the posteriorly in early stage 4 embryos (Figure 2E). The bio-
primary determinant of anterior cell fate specification logical significance of these subcategories is under-
(Ephrussi and St Johnston, 2004). The other four mRNAs scored by their specific GO term enrichments (Figure 2F
exhibit a more diffuse gradient of anterior enrichment (Fig- and Table S3). For example, the category 1 and 2 pole
ure 2B). Interestingly, three of these encode proteins plasm localized transcripts, which include mRNAs such
involved in cytoskeleton organization and microtubule- as aret, eIF5, gcl, Imp, nos, osk, orb, pAbp, pum, spir,
based processes (Figure 2F and Table S3). and Tm1, are strongly enriched for cell development,
translation regulation, pole plasm assembly and RNA Subcellular Categories: Apicobasal Patterns
localization functions (Figure 2F). In contrast, the late Besides the subembryonic localization patterns, a large
posterior group, which includes genes such as aur, collection of mRNAs, either of maternal or zygotic origin,
CG14030, CycA, grp, gwl, Rbp-1, Su(var)3-9, and ttk, were found to exhibit intricate subcellular localization
is enriched for protein kinases and negative regulators patterns. Classic examples include the subset of mRNAs
of gene expression, again consistent with previous find- that localize to the apical cytoplasm within the embryonic
ings that germ cells are transcriptionally silent (Seydoux epithelium (Figure 3). Although apical mRNAs have been
and Braun, 2006; Seydoux and Dunn, 1997; Van Doren characterized previously and considered as a homoge-
et al., 1998). Taken together, these observations sug- neous group (Davis and Ish-Horowicz, 1991; Simmonds
gest the existence of distinct early and late pathways for et al., 2001; Tepass et al., 1990), many distinctive sub-
posterior transcript localization. groups of apical transcripts could be distinguished, rang-
Notably, no maternal transcripts were identified as be- ing from broad gradients of apical enrichment to tightly
ing either dorsally or ventrally localized. Instead, differen- localized clusters or foci (Figures 3A–3E). Likewise,
tial distribution of transcripts along the dorso-ventral axis a large number of basally localized mRNAs were identi-
was always a consequence of localized zygotic tran- fied, which also fall into a number of subgroups (Figures
scription. The preponderance of transcripts localized to 3G–3I). Other patterns that vary along the apico-basal
the posterior pole of the embryo, in comparison to the axis include transcripts that are excluded from the apical
other embryonic poles, seemingly reflects special re- cytoplasm or from the entire blastoderm layer (Figures 3J
quirements for germ cell specification and the sufficiency and 3K).
of existing zygotic mechanisms to define the other coor- The functional relevance of these subgroup classifica-
dinates. tions is underlined by the GO term enrichments observed.
For example, the apical clusters category, which includes or perinuclear clouds as early as stage 3, and then evolve
mRNAs such as Ama, bib, Btk29A, crb, fra, Gp150, htl, Ptr, into polygonal mosaic networks shortly thereafter (Figures
scb, sog, and smo, is enriched for GO categories for 4A and 4C). These patterns resemble subsequent actin
plasma membrane and signaling pathway components filament distributions and dynamics and precede cell
(Figure 3M and Table S4). In contrast, the diffuse apical junction formation. Several other mRNAs encoding cell
group, which contains several pair-rule gene transcripts junction components, such as Patj (Figure 4D) and dlg1
(hairy, odd, prd, run), is enriched for functions associated (Figure 4E), localize at specific sites along the basolateral
with transcriptional regulation and pattern/axis specifica- membrane. In contrast, mira transcripts localize through-
tion. out the lateral membrane of embryonic epithelial cells
In addition to the apicobasal patterns detected in the (Figure 4F). Accordingly, this category is enriched for
embryonic epithelium, several mRNAs were observed GO terms related to cytoskeleton organization and biogen-
with asymmetric patterns in neuroblasts. This category in- esis (Table S5). These observations imply a significant
cludes transcripts such as asp, Gp150, mira, odd, pros, role for mRNA localization in the nucleation and position-
and wg, some of which have been observed previously ing of cytoskeletal networks and membrane-associated
(Broadus et al., 1998; Schuldt et al., 1998), and not sur- structures.
prisingly, show GO term enrichments for asymmetric cell
division functions (Figure 3M). We also identified mRNAs Cell Division and Nuclei-Associated Patterns
from uncharacterized genes, such CG14896, which Many of the most striking subcellular patterns observed
exhibit apical localization both in the posterior embryonic occur during nuclear or cellular division (Figure 5). These
epithelium underlying the pole cells and in neuroblasts include transcripts that localize to spindle poles, centro-
that arise later in embryogenesis (Figures 3F and 3L). somes/centrioles, astral microtubules, or along the mitotic
This example suggests the likelihood that many spindles themselves during anaphase and telophase (Fig-
localization mechanisms will operate in a variety of cell ures 5A–5H). Furthermore, several mRNAs that are zygoti-
types. cally transcribed in early stage 4 embryos concentrate
around metaphase chromosomes during mitosis and of-
Membrane-Associated Patterns ten become associated with spindle midbodies (Figure
Also remarkable are transcripts that localize to mem- 5D). Intriguingly, many of these mRNAs encode trans-
brane-associated structures prior to and following cellula- criptional regulators (dpld, nullo, odd, rib, run, stwl, Taf4).
rization (Figure 4). For example, cno and anillin mRNAs The genes in these categories show GO term enrich-
(Figures 4A and 4B) associate with the embryonic cortex ments for cell division related processes (Figure S3 and
Table S5), implying important roles for localized mRNAs in homologous and in close proximity to the tandemly re-
the establishment, function and regulation of cell division peated Stellate gene cluster on the X chromosome, which
machineries. has been implicated in the maintenance of male fertility
Interestingly, several of the earliest zygotically ex- through an RNA interference process involving the Su(Ste)
pressed mRNAs, such as those encoded by the Doc gene cluster located on the Y chromosome (Aravin et al.,
and copia transposons, exhibit intricate chromatin-asso- 2001). Finally, roX1, a noncoding RNA involved in X chro-
ciated patterns. Indeed, Doc-element RNA localizes in mosome dosage compensation (Park et al., 2002), local-
the vicinity of centromeres, either along the ‘rosettes’ izes to the basal side of blastoderm nuclei, where the
formed by the polar body chromosomes or in dividing dip- X chromosome presumably resides (Figure 5M). As has
loid nuclei (Figures 5I and 5J). In contrast, mRNA encoded been demonstrated for roX1, these assorted DNA-associ-
by the Ste12DOR gene is found within large chromatin- ated RNAs may be functioning to help organize, establish
associated foci that localize to telomeric regions during and/or maintain chromatin domains. For Doc and Ste12-
anaphase (Figures 5K and 5L). Notably, this gene is highly DOR, the chromosome-associated RNA may be
functioning in the ‘repeat-associated small interfering staining for selected transcripts, their protein products
RNA’ (rasiRNA) pathway, which acts in part to suppress or relevant markers was carried out (Figure 6). Examples
transposable element activity (Slotkin and Martienssen, are shown for mRNAs that are apically localized (Figure
2007). If so, this autoregulation would add a new dimen- 6A), cell division apparatus-associated (Figures 6B and
sion to our understanding of this process. 6C), or that reside at cell junctions (Figure 6D). In each
Finally, other nuclei-associated mRNAs were observed of these cases, mRNA localization is noted prior to the
that range from mRNAs with tight perinuclear localization appearance of protein, consistent with the view that
(Figure 5N) to those that appear to be uniformly localized transcript localization generally predetermines protein
throughout the nucleus (Figures 5O and 5P). These in- distribution at most subcellular destinations. Another in-
clude cas, CG15634, CG8552, Eip71CD, hb, Jra, kuk, triguing example is CG14438 mRNA, which appears to lo-
mfas, and scw. Interestingly, some of these are only re- calize at the level of centrioles and is found nestled within
tained within nuclei that appear to be dropping out of structures labeled with the pericentriolar marker CNN
the blastoderm layer (Figures 5P and 6E). As these nuclei (Figure 6C). These types of examples support the notion
are generally observed in pairs, and ‘nuclear fallout’ may that localized transcripts, and the proteins they encode,
be a consequence of unsuccessful nuclear divisions play central nucleation functions within the cell. These ob-
(Rothwell et al., 1998), this suggests a potential function servations further suggest that the necessary translation
of these mRNAs in nuclear migration, sorting and/or and secretory machineries are generally available at each
apoptosis. of these sites.
Interestingly, a reverse correlation exists for mRNAs
Colocalization of RNAs and Proteins that show nuclear retention, as shown for kuk and cas
In many of the cases cited above, where details are known transcripts (Figures 6E and 6F). Indeed, while Kuk protein
about protein localization or gene function, there is a strik- exhibits robust nuclear envelop localization in cortical
ing correlation between transcript localization and the nuclei where the mRNA is cytoplasmic, no protein is
patterns or functions of the encoded proteins (Table S6). observed in or around the yolk nuclei in which kuk mRNA
To further illustrate some of these relationships, double- is retained (Figure 6E). Similarly, in stage 9 neuroblasts,
Kambadur, R., Koizumi, K., Stivers, C., Nagle, J., Poole, S.J., and Rubin, G.M., Hong, L., Brokstein, P., Evans-Holm, M., Frise, E., Staple-
Odenwald, W.F. (1998). Regulation of POU genes by castor and ton, M., and Harvey, D.A. (2000). A Drosophila complementary DNA
hunchback establishes layered compartments in the Drosophila resource. Science 287, 2222–2224.
CNS. Genes Dev. 12, 246–260. Schuldt, A.J., Adams, J.H., Davidson, C.M., Micklem, D.R., Haseloff,
Karr, T.L., and Alberts, B.M. (1986). Organization of the cytoskeleton in J., St Johnston, D., and Brand, A.H. (1998). Miranda mediates asym-
early Drosophila embryos. J. Cell Biol. 102, 1494–1509. metric protein and RNA localization in the developing nervous system.
Kloc, M., Wilk, K., Vargas, D., Shirato, Y., Bilinski, S., and Etkin, L.D. Genes Dev. 12, 1847–1857.
(2005). Potential structural role of non-coding and coding RNAs in Seydoux, G., and Braun, R.E. (2006). Pathway to totipotency: lessons
the organization of the cytoskeleton at the vegetal cortex of Xenopus from germ cells. Cell 127, 891–904.
oocytes. Development 132, 3445–3457.
Seydoux, G., and Dunn, M.A. (1997). Transcriptionally repressed germ
Kloc, M., Zearfoss, N.R., and Etkin, L.D. (2002). Mechanisms of sub- cells lack a subpopulation of phosphorylated RNA polymerase II in
cellular mRNA localization. Cell 108, 533–544. early embryos of Caenorhabditis elegans and Drosophila mela-
Lambert, J.D., and Nagy, L.M. (2002). Asymmetric inheritance of cen- nogaster. Development 124, 2191–2201.
trosomally localized mRNAs during embryonic cleavages. Nature 420, Siegrist, S.E., and Doe, C.Q. (2005). Microtubule-induced Pins/Galphai
682–686. cortical polarity in Drosophila neuroblasts. Cell 123, 1323–1335.
Lawrence, J.B., and Singer, R.H. (1986). Intracellular localization of
Simmonds, A.J., dosSantos, G., Livne-Bar, I., and Krause, H.M. (2001).
messenger RNAs for cytoskeletal proteins. Cell 45, 407–415.
Apical localization of wingless transcripts is required for wingless
Lécuyer, E., Parthasarathy, N., and Krause, H.M. (2007). Fluorescent In signaling. Cell 105, 197–207.
Situ Hybridization Protocols in Drosophila Embryos and Tissues. In
Slotkin, R.K., and Martienssen, R. (2007). Transposable elements and
Methods Mol. Biol., C. Dahmann, ed. (Totowa, NJ: Humana Press,
the epigenetic regulation of the genome. Nat. Rev. Genet. 8, 272–285.
Inc.), pp. 289–302.
Li, P., Yang, X., Wasser, M., Cai, Y., and Chia, W. (1997). Inscuteable St Johnston, D. (2005). Moving messages: the intracellular localization
and Staufen mediate asymmetric localization and segregation of pros- of mRNAs. Nat. Rev. Mol. Cell Biol. 6, 363–375.
pero RNA during Drosophila neuroblast cell divisions. Cell 90, 437–447. Stapleton, M., Liao, G., Brokstein, P., Hong, L., Carninci, P., Shiraki, T.,
Long, R.M., Singer, R.H., Meng, X., Gonzalez, I., Nasmyth, K., and Jan- Hayashizaki, Y., Champe, M., Pacleb, J., Wan, K., et al. (2002). The
sen, R.P. (1997). Mating type switching in yeast controlled by asym- Drosophila gene collection: identification of putative full-length cDNAs
metric localization of ASH1 mRNA. Science 277, 383–387. for 70% of D. melanogaster genes. Genome Res. 12, 1294–1300.
Mattick, J.S. (2004). RNA regulation: a new genetics? Nat. Rev. Genet. Tadros, W., Goldman, A.L., Babak, T., Menzies, F., Vardy, L., Orr-
5, 316–323. Weaver, T., Hughes, T.R., Westwood, J.T., Smibert, C.A., and Lipshitz,
Melton, D.A. (1987). Translocation of a localized maternal mRNA to the H.D. (2007a). SMAUG is a major regulator of maternal mRNA destabi-
vegetal pole of Xenopus oocytes. Nature 328, 80–82. lization in Drosophila and its translation is activated by the PAN GU
kinase. Dev. Cell 12, 143–155.
Mingle, L.A., Okuhama, N.N., Shi, J., Singer, R.H., Condeelis, J., and
Liu, G. (2005). Localization of all seven messenger RNAs for the ac- Tadros, W., Westwood, J.T., and Lipshitz, H.D. (2007b). The Mother-
tin-polymerization nucleator Arp2/3 complex in the protrusions of to-Child Transition. Dev. Cell 12, 847–849.
fibroblasts. J. Cell Sci. 118, 2425–2433. Takizawa, P.A., Sil, A., Swedlow, J.R., Herskowitz, I., and Vale, R.D.
Mueller, W., Nutt, C.L., Ehrich, M., Riemenschneider, M.J., von Deim- (1997). Actin-dependent localization of an RNA encoding a cell-fate
ling, A., van den Boom, D., and Louis, D.N. (2007). Downregulation of determinant in yeast. Nature 389, 90–93.