Introduction
Stem cells are unique in that they both self-renew and give rise to a variety of differentiated cell types. Embryonic stem cells (ESCs) are isolated from the inner cell mass of mammalian blastocysts and propagated in culturing systems (
Evans and Kaufman, 1981;
Martin, 1981). ESCs can be induced to differentiate into distinct cell types and therefore hold the promise for therapeutic applications (
Smith, 2001). More excitingly, in both mouse and human, differentiated somatic cells can be reprogrammed to become ESC-like iPSCs (
Takahashi and Yamanaka, 2006;
Takahashi et al., 2007;
Yu et al., 2007;
Park et al., 2008). Moreover, mouse male germ cells have been shown to retain the potential to become ESC-like cells
in vitro and
in vivo (
Kanatsu-Shinohara et al., 2004;
Guan et al., 2006). These exciting discoveries, which involve the reversal of lineage commitment for both somatic and germ cells, open an entirely new avenue for applying stem cells in regenerative medicine.
Another type of stem cell is the naturally existing adult stem cells, which are required for tissue homeostasis and regeneration in adulthood (
Morrison and Spradling, 2008;
Losick et al., 2011). The adult stem cell normally divides asymmetrically, either as a single cell or as a group of cells, to self-renew and give rise to daughter cells that initiate differentiation (
Yamashita et al., 2010). Before turning on the terminal differentiation program, progenitor cells derived from stem cells normally undergo proliferation to expand their population for high-throughput yield of terminally differentiated cell types. Such activities of stem cells need to be tightly controlled. Imbalances that arise between stem cell self-renewal versus differentiation, or between progenitor cell proliferation versus differentiation, are common causes of many human diseases, including cancers, tissue dystrophy and infertility (
Morrison and Kimble, 2006;
Rossi et al., 2008).
To effectively use stem cells in regenerative medicine, the molecular mechanisms underlying stem cell maintenance, proper differentiation, and reprogramming need to be thoroughly understood. Particularly, the stem cell transcriptome needs to be resolved at the levels of both transcript abundance and isoform structure. In addition, transcriptional network and chromatin landscape need to be delineated to understand the regulatory circuitries that maintain stem cell pluripotency and prevent precocious differentiation. Indeed, recent years have witnessed substantial progress in these areas (
Chen, 2008).
Recently, new technologies, including exon tiling arrays and high-throughput mRNA sequencing (RNA-seq), have revolutionized our understanding of transcription and alternative splicing (AS) in various stem cell lineages, including ESCs and adult stem cells in different organisms. In particular, several new studies have revealed that differential gene expression in stem cell lineages involves mRNA isoform switching from stem cells to differentiated cells, ranging from ESCs to various adult stem cells, and across the animal kingdom from
Drosophila to humans (
Gan et al., 2010;
Pritsker et al., 2005;
Yeo et al., 2007;
Salomonis et al., 2009;
Gabut et al., 2011;
Salomonis et al., 2010;
Wu et al., 2010). This switching either leads to changes in particular coding sequences, which may affect protein structure and/or subcellular localization, or it results in changes at the non-coding sequences, which may alter post-transcriptional regulation, such as microRNA-mediated regulatory mechanisms or nonsense-mediated mRNA decay (NMD) (
Lareau et al., 2007;
Ni et al., 2007;
Barash et al., 2010). In this review, we discuss these recent results and hypothesize that mRNA isoform switching may serve as a common molecular mechanism that acts cooperatively with other epigenetic mechanisms, such as histone modifications, in regulating stem cell function.
Alternative splicing in the regulation of ESC pluripotency and lineage-specific differentiation
Molecular mechanisms underlying the pluripotency of ESCs have been extensively investigated, but most studies have focused on dissecting the key transcription factor-regulated transcriptional network (
Boyer et al., 2005;
Boyer et al., 2006a;
Loh et al., 2006;
Zhou et al., 2007) and chromatin structure (
Azuara et al., 2006;
Bernstein et al., 2006;
Boyer et al., 2006b;
Lee et al., 2006b;
Guenther et al., 2007;
Mikkelsen et al., 2007;
Stock et al., 2007;
Kim et al., 2008). However, recent studies using new technologies, including RNA-seq, have shown the importance of ESC-specific mRNA isoforms (
Kunarso et al., 2008;
Das et al., 2011), as well as changes of isoforms during ESC differentiation (
Pritsker et al., 2005;
Yeo et al., 2007;
Salomonis et al., 2009;
Salomonis et al., 2010;
Wu et al., 2010;
Gabut et al., 2011).
Interestingly, several studies have demonstrated high diversity of isoforms in stem cells, but low diversity of isoforms in differentiated cells, a phenomenon called “isoform specialization” (
Wu et al., 2010). To account for this, it is speculated that stem cells may require a high diversity of isoforms to maintain their identities, while reduction into specific isoforms in differentiating cells may ensure proper differentiation. The diverse isoforms co-expressed in stem cells have non-redundant functions. For example, a pluripotency transcription factor, Sal-like protein 4 (Sall4), has two isoforms,
Sall4a (long isoform) and
Sall4b (short isoform), in mouse ESCs (
Rao et al., 2010). Although both Sall4a and Sall4b interact with one of the core transcription factors Nanog in ESCs (
Wu et al., 2006), they have distinct roles. Specifically, Sall4a represses expression of differentiation genes, while Sall4b essentially maintains high expression of pluripotency genes. This phenomenon whereby a critical transcription factor sustains multiple isoforms in ESCs also applies to the Wnt-responsive transcription factor TCF3. Here, again, both
TCF3(l) (long isoform) and
TCF3(s) (short isoform) are co-expressed in ESCs and regulate expression of different sets of target genes (
Salomonis et al., 2010). Other examples include critical cell cycle regulators, such as the mitogen-activated protein kinase (MAPK) pathway components MAP4K2 and MNK1 kinases, which also have different isoforms in stem cells (
Pritsker et al., 2005). In these examples, it is possible that coordinated activities of proteins encoded by different isoforms are required for maintaining ESC pluripotency.
In other cases, changes in alternative splicing during ESC differentiation are accompanied by isoform switching, in which stem cells and differentiated cells have distinct isoforms. Isoform switching is probably required for lineage specification during ESC differentiation. An earlier study using the EST (expressed sequence tag) library showed that genes required for lineage-specific differentiation undergo AS at a higher frequency than housekeeping genes (
Pritsker et al., 2005). Recent RNA-seq data obtained from distinct human organs have also revealed abundant isoform switching for tissue-specific genes (
Wang et al., 2008).
Isoform switching during ESC differentiation has been shown to mainly affect coding sequences, resulting in protein domain truncation, disruption, exchange, or modification (
Salomonis et al., 2010). Interestingly, many of the genes that show this isoform switching phenomenon encode critical regulators. For example, the master transcription factor OCT4 that regulates ESC self-renewal has an ESC-specific isoform (
OCT4A). Upon differentiation,
OCT4A expression is greatly reduced, while two other isoforms,
OCT4B and
OCT4B1, take over (
Atlasi et al., 2008). It has been shown that OCT4B has distinct subcellular localization and DNA binding affinity different from OCT4A (
Lee et al., 2006a), suggesting that OCT4B and OCT4A regulate different target genes in differentiating cells and stem cells, respectively. A recent study reveals that the Forkhead box (FOX) transcription factor FOXP1 has an ESC-specific isoform, which is different from the more broadly expressed isoform (
Gabut et al., 2011). The particular switching of an exon in
FOXP1 transcript results in amino acid changes within the critical forkhead DNA binding domain, and this leads to different choices of target genes. In this case, the ESC-specific FOXP1 (FOXP1-ES) acts as a transcriptional activator for pluripotency genes, such as OCT4 and NANOG, while it acts as a transcriptional repressor for differentiation genes. Interestingly, this dual role of FOXP1-ES is also required for efficient reprogramming of somatic cells to iPSCs, suggesting a critical role for this isoform switching in maintaining, as well as restoring, ESC function. In other cases, switched mRNA isoforms even encode proteins with antagonistic functions. For example, fibroblast growth factor 4 (FGF4) acts as an important autocrine factor to maintain human ESC pluripotency, whereas its truncated isoform encoded by
FGF4si antagonizes FGF4 function and promotes differentiation (
Mayshar et al., 2008). Genome-wide analysis reveals different AS patterns in different human tissues and cell lines (
Sultan et al., 2008;
Wang et al., 2008). Further functional studies in model organisms are needed to specifically knockdown or overexpress particular isoforms to determine the biologic significance of these changes.
Isoform switching can also lead to changes at regulatory regions, such as 3′UTRs (Untranslated Regions). Different 3′UTRs may subsequently allow post-transcriptional regulation via distinct microRNAs. For instance, the
Serca2 gene encodes a calcium pump essential for cardiac muscle cell function. Its protein expression is repressed by microRNAs in ESCs. However, during mouse ESC cardiac differentiation, a specific
Serca2 isoform with shorter 3′UTR is produced through exon exclusion, which allows it to escape microRNA binding and repression (
Salomonis et al., 2010). In addition to microRNAs, NMD adds another mechanism to regulate transcript abundance in stem cells and differentiated cells. For example, a particular set of premature termination codon-introducing exons are excluded in stem cells or embryonic tissues in order to maintain gene expression at a high level, whereas they are included in differentiated tissues to activate NMD and reduce expression level for the corresponding genes (
Barash et al., 2010).
Alternative splicing in the regulation of adult stem cell maintenance and differentiation
Adult stem cells are multi-potent or uni-potent cells that exist under physiologic conditions and function to maintain tissue homeostasis by replenishing short-lived cells, such as skin, muscle, intestinal, and blood cells, as well as sperm (
Morrison and Spradling, 2008). Understanding the molecular characteristics of adult stem cells will greatly facilitate targeted differentiation of ESCs to lineage-committed adult stem cells, or the reprogramming of differentiated cells to adult stem cells. This could be followed by transplantation of adult stem cells to their physiologic environments, thus serving as basis for curing genetic diseases and age-dependent tissue degeneration. Compared to ESCs, very little is known about how alternative splicing regulates adult stem cell maintenance and differentiation, partly because of their extremely small numbers
in vivo, as well as technical difficulties associated with their isolation and purification (
Eun et al., 2010).
The
Drosophila germline stem cells (GSCs) provide paradigmatic model systems to study adult stem cell identity and activity (
Fuller and Spradling, 2007;
Losick et al., 2011). During
Drosophila male GSC differentiation, AS events are significantly more abundant in undifferentiated cell-enriched
bag-of-marbles (
bam) mutant testes than in
wild-type testes which mainly contain differentiating cells (
Gan et al., 2010). The infrequency of splicing events in differentiating germ cells is partly achieved by turning on a large set of single-isoform or intronless genes, which probably serves as a mechanism to ensure robust expression of nearly a thousand genes for terminal differentiation (
Gan et al., 2010). However, isoform switching of multi-isoform genes is also detected at the transition from mitotic undifferentiated germ cells to meiotic differentiating germ cells [(
Gan et al., 2010) and unpublished results]. Further analyses will be performed to address their functions in regulating spermatogenesis.
Hematopoiesis and myogenesis are two other well-characterized cellular differentiation processes that rely on the activities of adult stem cells. Genome-wide studies of the mammalian hematopoietic stem cell lineage reveal that lineage-specific genes and signaling pathway components have a higher propensity to undergo splicing than housekeeping genes, providing great opportunities for functional studies (
Lemischka and Pritsker, 2006;
Pritsker et al., 2005). For example, the transcription factor-encoding
Ikaros (
Ik) gene undergoes dynamic AS to produce various isoforms that restrict hematopoietic stem cells to differentiate along the lymphoid pathway. Distinct
Ik isoforms differ in their protein-coding sequences, resulting in distinct subcellular localizations of Ik proteins responsible for their different transcription activation abilities (
Molnár and Georgopoulos, 1994).
Using splicing-sensitive microarray (
Bland et al., 2010) and RNA-seq (
Trapnell et al., 2010) techniques, recent studies have also demonstrated robust splicing switching of several hundred genes during mammalian myoblast differentiation. Subsequent studies reveal that the majority of these genes switch isoforms in coordination with the myogenic differentiation program and might be temporally controlled. Moreover, the splicing switching during myogenesis occurs even before expression of myogenic transcription factors, suggesting that isoform switching may poise myoblast progenitor cells for myogenic differentiation (
Bland et al., 2010). Finally, isoform changes for many genes, including the crucial Myc transcription factor, are so dynamic that they cannot be ascribed to a simple “pattern switching” (
Trapnell et al., 2010). Therefore, more detailed analyses of individual genes are required to thoroughly understand how AS regulates muscle stem cell differentiation.
Regulation of isoform switching during stem cell differentiation by splicing factors
The regulation of AS pattern switching in stem cell lineages at the molecular level remains an open question. However, we do know that AS is regulated by RNA-protein complexes called spliceosomes, consisting of five small nuclear RNAs (U1, U2, U4, U5, and U6 snRNAs) and hundreds of protein components (
Jurica and Moore, 2003;
Wahl et al., 2009). The
trans-acting RNA binding proteins (RBPs) act by directly recognizing and binding to splicing regulatory elements in pre-mRNAs, either as splicing enhancers or silencers, to determine inclusion or exclusion of a particular exon (
Wang and Burge, 2008). Therefore, the ultimate choice of AS sites is defined by interactions between
trans-acting splicing factors and
cis-acting splicing elements at pre-mRNAs.
Furthermore, several splicing factors have been demonstrated to regulate AS events in stem cell maintenance and differentiation. For example, using CLIP-seq (cross-linking and immunoprecipitation followed by high-throughput sequencing) technology, the FOX2 splicing factor has been shown to bind to sequences that encode many splicing regulators in human ESCs, demonstrating a potential upstream role of FOX2 in regulating AS in ESCs (
Yeo et al., 2009). Consistent with its important roles in ESCs, knockdown of FOX2 leads to ESC cell death (
Yeo et al., 2009). The conserved FOX1/2 binding motif has also been found near alternatively spliced sites for genes expressed in neural progenitor cells (
Yeo et al., 2007) and various human tissues (
Wang et al., 2008;
Zhang et al., 2008), suggesting that FOX1/2 might have a broad role in regulating tissue-specific AS events. Another splicing factor, polypyrimidine tract binding protein (PTB), is important for ESC proliferation and early embryonic development (
Shibayama et al., 2009), and it also regulates both neuronal and muscle cell differentiation. In the central nervous system, PTB is restrictively expressed in neuronal precursor cells, whereas its paralogue nPTB is specifically expressed in post-mitotic neurons. Since PTB and nPTB regulate distinct AS events, they are major regulators of isoform switching during neural differentiation (
Boutz et al., 2007b). In muscle cells, PTB and nPTB are highly expressed in precursor myoblasts where they repress inclusion of certain exons of target genes, whereas in differentiated myotubes, expression of PTB and nPTB is repressed, leading to inclusion of the particular exons for target genes (
Boutz et al., 2007a). Another example is the 68 kDa Src substrate associated during mitosis (SAM68) RBP, which also regulates AS events critical for neural differentiation (
Chawla et al., 2009) and mesenchymal cell differentiation (
Richard et al., 2005).
Intriguingly, the splicing factors themselves might be subject to dynamic regulation as well, adding another layer of complexity. For example, in the
Drosophila male germline lineage, splicing factors are highly enriched in undifferentiated cells, but they are downregulated in differentiating cells, suggesting that their expression levels are controlled by developmental programming (
Gan et al., 2010). In addition, during both neural and cardiac differentiation of ESCs, genes regulating RNA splicing themselves undergo robust AS (
Salomonis et al., 2009), for example, PTB and nPTB which have mutually exclusive expression in the nervous system. This expression pattern is partly regulated by AS. That is, in precursor cells, PTB induces splicing of nPTB into an isoform that subsequently undergoes NMD (
Boutz et al., 2007b). In addition, during myoblast differentiation, the level of splicing regulators also changes dynamically (
Bland et al., 2010). For instance, nPTB is downregulated during muscle differentiation by a muscle-specific microRNA, miR-133, providing a great example illustrating how these two post-transcriptional mechanisms (i.e., AS and microRNA) coordinate in regulating cellular differentiation (
Bland and Cooper, 2007;
Boutz et al., 2007a). These observations, together with the identification of many splicing factor binding sites at regions that are differentially spliced, demonstrate that AS events are tightly controlled during cellular differentiation and that the expression or AS of splicing factors may contribute to this.
Regulation of isoform switching during stem cell differentiation by epigenetic mechanisms: splicing code vs. epigenetic code
Increasing evidence demonstrates that splicing occurs concurrently with mRNA transcription (
Allemand et al., 2008). This feature allows regulators of transcription, such as chromatin factors, to control splicing. Recent studies demonstrate two major epigenetic mechanisms in regulating AS: histone modifications and control of RNA Pol II processivity. Because histone modifications also play a profound role in regulating Pol II processivity, these two mechanisms may exhibit crosstalk.
In stem cells, modified histones could determine choice of splicing sites by directly recruiting chromatin-associated proteins and splicing factors. One example comes from studying chromodomain-helicase-DNA binding protein-1 (CHD1), which is required for maintaining the chromatin landscape and pluripotency of mouse ESCs (
Gaspar-Maia et al., 2009). CHD1 recognizes the H3K4me3 histone modification. Interestingly, CHD1 physically interacts with spliceosome components. Furthermore, depletion of CHD1 results in malfunction of spliceosomes, which leads to failure of AS at a set of target genes (
Sims et al., 2007). Another example comes from studying human mesenchymal stem cells. In this case, the alternatively spliced third exon of the
fibroblast growth factor receptor 2 gene (
FGFR2) is labeled by H3K36me3, a histone modification for transcriptional elongation (
Barski et al., 2007). The H3K36me3 mark is subsequently recognized by the histone tail binding MORF-related protein 15, followed by the recruitment of PTB to a splicing silencer element, which regulates the exclusion of this exon in mesenchymal stem cells (
Luco et al., 2010). Since splicing is coupled with transcription, the particular histone modifications that act with splicing machinery are usually associated with transcriptional activation, such as the H3K4me3 and the H3K36me3 histone modifications. We anticipate that future studies will address how epigenetic mechanisms contribute to specific splicing codes for stem cell pluripotency, as well as their proper differentiation.
The processivity of Pol II determines the rate of transcription and may also affect choices of AS, as has been reported in stressed cells (
Muñoz et al., 2009), neurons (
Schor et al., 2009), and small interfering RNA (siRNA)-treated cells (
Alló et al., 2009). It has been shown that Pol II is more abundant at the exonic region than at introns, suggesting that slow movement of Pol II along exons may allow precise determination of exon-intron boundaries (
Schwartz et al., 2009). On the other hand, nucleosomes are found to have a higher occupancy at exons than at introns in various organisms [reviewed by (
Schwartz and Ast, 2010)]. Because high nucleosome density limits Pol II elongation rate, differential occupancy of nucleosome along the genome of different cells may lead to their distinct splicing pattern. Intriguingly, recent studies reveal that a DNA binding protein, CTCF, promotes inclusion of a weak exon by pausing Pol II at specific splicing sites. On the other hand, such a pause of Pol II is prevented by DNA methylation through inhibition of CTCF binding, providing strong evidence that epigenetic mechanisms, such as DNA methylation, regulate AS through a DNA binding protein (
Shukla et al., 2011). Based on these reports, it is likely that epigenetic mechanisms not only determine gene expression level but also regulate splicing pattern. However, since most of the data have been obtained in cultured cells, more research will be needed to understand how epigenetic mechanisms regulate stem cell-specific splicing events in various organisms [Fig. 1 and reviews (
Eun et al., 2010;
Luco et al., 2011;
Luco and Misteli, 2011)].
Differential splicing and information flow in stem cell lineages
In the previous sections, we reviewed the involvement of AS in stem cell maintenance and differentiation. In this section, we will attempt to unify these different events using information and communication theory metaphors.
Recently, some basic concepts from information theory have been used to characterize the complexity and dynamics of alternative splicing. For example, it was reported that isoform entropy was significantly higher in cancer cells than in normal tissues (
Ritchie et al., 2008). We also observed that isoform entropy decreases upon differentiation of male GSCs in
Drosophila (
Gan et al., 2010). A related observation in ESC differentiation was made in (
Wu et al., 2010). To describe the generic isoform switching that should involve at least two isoforms, a recent paper uses “Jensen-Shannon distance” as an information-theoretic approach to characterize changes in information content of transcriptomes during cellular differentiation (
Trapnell et al., 2010).
The signals received by stem cells from their niches allow them to maintain their identity and activity. These signals comprise different types of interactions such as direct cell-cell contacts and cell-extracellular matrix contacts. Shannon entropy is a measure of uncertainty about a system (
Cover and Thomas, 1991). We speculate that decrease in isoform entropy is related to information inflow that guides the stem cell differentiation process. The pre-mRNA splicing process is tightly coupled to the cellular information processing systems (CIPS). Therefore, the functional meaning of AS in stem cell maintenance and differentiation can be precisely understood only in the context of the CIPS.
The pre-mRNA splicing process in the context of the CIPS is schematically depicted in Fig. 1A. The extracellular signals activate cell-surface receptors, which, in turn, relay the information down the signal transduction pathways into the stem cell. The cell-surface receptors and the signal transduction pathways are collectively denoted as “sensor.” The sensor also receives input from intracellular processes; thus, the output of the sensor also depends on the intrinsic state of the cell. The signal transduction pathways impinge upon transcriptional, post-transcriptional, and post-translational gene-regulatory networks, which are collectively denoted as “controller.” This is depicted as the output of the sensor feeding into the controller. Transcription factors or chromatin modifiers that bind to the regulatory elements, such as the enhancer and promoter regions of genes, or modify chromatin are themselves subject to regulation by the controller (shown as a thick red line with an arrow toward transcription factors or chromatin modifiers). Based on the co-transcriptional nature of alternative splicing of pre-mRNAs, some of the histone modifications involved in transcription initiation and elongation may influence pre-mRNA splicing process. The transmission of information from the chromatin template to pre-mRNA of a gene is depicted as arrows extending from “gene” to “pre-mRNA.” The binding of RBPs, such as FOX2 and PTB discussed above, to splicing regulatory sequences, such as splicing enhancers or silencers, modulate pre-mRNA splicing. RBPs are themselves subject to regulation by the controller.
The output of pre-mRNA splicing is a set of spliced mRNA isoforms (denoted as 1 to 7 in Fig. 1A). The relative levels of these isoforms are determined by upstream regulatory processes. In other words, the information inflow from upstream processes shapes the information content of mRNA isoforms. The spliced isoforms then pass through various “filters.” NMD can be interpreted as a filter (e.g., isoform 5 does not pass the NMD filter and gets degraded in Fig. 1A). Translational inhibition or mRNA degradation by microRNAs can be treated as another filter (e.g., isoforms 3 and 6 do not pass the microRNA filter and get degraded and/or are translationally inhibited in Fig. 1A). The Serca2 gene in mouse ESC, as discussed above, is an example of a microRNA filter. The microRNAs are themselves subject to regulation by the controller. The products of mRNA isoforms that pass all post-transcriptional and post-translational filters go on to produce functional proteins that act in various cellular processes. Examples of such cellular processes include maintenance of pluripotency and repression of differentiation gene expression in stem cells. The outputs of these cellular processes then have feedback to the CIPS.
The schematic framework shown in Fig. 1A may accommodate many scenarios regarding how AS contributes to the self-renewal and pluripotency of stem cells, as well as their cellular differentiation. All examples discussed in this review can be understood using this framework. Some additional possible scenarios are shown in Fig. 1B.
The information/communication theory metaphor of CIPS described here is qualitative and does not necessarily reflect all events at the molecular level. In order for such a framework to be useful, a quantitative theory of information flow dynamics in CIPS is needed. We predict that as more high-throughput measurements of the dynamic activity of components in CIPS become available in cellular differentiation, the information-theoretic methods will be increasingly used to characterize information flow in biologic networks.
Conclusions and perspectives
In summary, new high-throughput technology allows digital inventory of RNA abundance and isoforms, which provides new opportunities to study how changing AS patterns accompany or regulate stem cell differentiation processes. Although most of the current work cannot determine whether the observed pattern switching correlates with, or contributes to, stem cell maintenance and differentiation (
Nelles and Yeo, 2010), molecular genetics assays using animal models or cultured stem cells can test the biologic consequence of disrupting such patterns in stem cell systems,
in vivo or
in vitro. Furthermore, recent research on human iPSCs has paved the way for regenerative medicine. It will be interesting to determine whether the stem cell-specific splicing code discussed in this review also exists in iPSCs. If it does, we will want to know what mechanism underlies isoform switching during normal stem cell differentiation such that it then reverses itself during the reprogramming process to become iPSCs. We anticipate more exciting work will emerge to understand the biologic significance of mRNA isoform switching during stem cell differentiation, thus facilitating the application of stem cells in regenerative medicine.
Higher Education Press and Springer-Verlag Berlin Heidelberg