Splicing in a single neuron is coordinately controlled by RNA binding proteins and transcription factors

  1. Morgan Thompson
  2. Ryan Bixby
  3. Robert Dalton
  4. Alexa Vandenburg
  5. John A Calarco
  6. Adam D Norris  Is a corresponding author
  1. Southern Methodist University, United States
  2. University of Toronto, Canada

Abstract

Single-cell transcriptomes are established by transcription factors (TFs), which determine a cell's gene-expression complement. Post-transcriptional regulation of single-cell transcriptomes, and the RNA binding proteins (RBPs) responsible, are more technically challenging to determine, and combinatorial TF-RBP coordination of single-cell transcriptomes remains unexplored. We used fluorescent reporters to visualize alternative splicing in single Caenorhabditis elegans neurons, identifying complex splicing patterns in the neuronal kinase sad-1. Most neurons express both isoforms, but the ALM mechanosensory neuron expresses only the exon-included isoform, while its developmental sister cell the BDU neuron expresses only the exon-skipped isoform. A cascade of three cell-specific TFs and two RBPs are combinatorially required for sad-1 exon inclusion. Mechanistically, TFs combinatorially ensure expression of RBPs, which interact with sad-1 pre-mRNA. Thus a combinatorial TF-RBP code controls single-neuron sad-1 splicing. Additionally, we find ‘phenotypic convergence,’ previously observed for TFs, also applies to RBPs: different RBP combinations generate similar splicing outcomes in different neurons.

https://doi.org/10.7554/eLife.46726.001

eLife digest

All the cells in the human nervous system contain the same genetic information, and yet there are many kinds of neurons, each with different features and roles in the body. Proteins known as transcription factors help to establish this diversity by switching on different genes in different types of cells.

A mechanism known as RNA splicing, which is regulated by RNA binding proteins, can also provide another layer of regulation. When a gene is switched on, a faithful copy of its sequence is produced in the form of an RNA molecule, which will then be ‘read’ to create a protein. However, the RNA molecules may first be processed to create templates that can differ between cell types: this means that a single gene can code for slightly different proteins, some of them specific to a given cell type. Yet, very little is known about how RNA splicing can generate more diversity in the nervous system.

To investigate, Thompson et al. developed a fluorescent reporter system that helped them track how the RNA of a gene called sad-1 is spliced in individual neurons of the worm Caenorhabditis elegans. This showed that sad-1 was turned on in all neurons, but the particular spliced versions varied widely between different types of nerve cells.

Additional experiments combined old school and cutting-edge genetics technics such as CRISPR/Cas9 to identify the proteins that control the splicing of sad-1 in different kinds of neurons. Despite not directly participating in RNA splicing, a number of transcription factors were shown to be involved. These molecular switches were turning on genes that code for RNA binding proteins differently between types of neurons, which in turn led sad-1 to be spliced according to neuron-specific patterns.

The findings by Thompson et al. could provide some insight into how mammals can establish many types of neurons; however, a technical hurdle stands in the way of this line of research, as it is still difficult to detect splicing in single neurons in these species.

https://doi.org/10.7554/eLife.46726.002

Introduction

The complement of genes expressed in an individual cell type controls its identity, development, and function. While transcriptional regulation is a major component of gene expression, post-transcriptional regulation can further shape cellular attributes by, for example, determining which gene isoforms are expressed in a cell. Much recent work has gone into cataloging gene expression networks in single cells, particularly those of specific neuronal types (Tanay and Regev, 2017; Zeng and Sanes, 2017). Molecular studies have also identified mechanisms by which transcription factors (TFs) shape gene expression networks in single neurons. Due to technical limitations, less is known about post-transcriptional regulation at the level of single neurons, or about the RNA binding proteins (RBPs) mediating post-transcriptional regulation (Gracida et al., 2016).

It is also unknown to what extent transcriptional and post-transcriptional gene regulatory networks are coordinated in single cells. A number of studies have identified individual RBPs that affect the splicing of a TF, thus altering the activity or specificity of that TF (Calarco et al., 2009; Han et al., 2013; Linares et al., 2015; Raj et al., 2011). These results suggest that there may be extensive cross-talk between transcriptional and post-transcriptional regulatory layers.

The nematode Caenorhabditis elegans has been used extensively as a model to reveal underlying principles by which TFs shape the transcriptomes of individual neurons. The worm’s invariant cell lineage, coupled with genetic tools and a transparent body, enables systematic in vivo analysis of gene expression in single neurons, and identification of TFs responsible for cell-specific gene expression. This type of analysis has revealed a number of gene regulatory principles, including the concept of a ‘combinatorial code’ of TFs which can be re-used in different neuron types, with particular combinations of TFs determining specific cell fates (Gendrel et al., 2016; Gordon and Hobert, 2015; Pereira et al., 2015). Another example is the concept of ‘phenotypic convergence’ by which various neurons express similar gene networks but the TFs driving the networks are different for each neuron type (Gendrel et al., 2016; Pereira et al., 2015). These principles appear to apply to the nervous systems of other organisms as well (Konstantinides et al., 2018). However, it remains unknown whether similar mechanistic principles apply to post-transcriptional regulation by RBPs in the nervous system.

Here we use single-cell in vivo fluorescent splicing reporters to investigate the cell-specific splicing of sad-1, a conserved neuronal kinase. The C. elegans sad-1 gene encodes two isoforms that differ in their ability to interact with the F-actin binding protein NAB-1/Neurabin (Hung et al., 2007), and have different roles in synapse formation and development (Kim et al., 2010). We find that sad-1 undergoes unique splicing patterns in various neuron types, and that developmentally-related cell types (the ALM touch-sensing neuron and the BDU neuron) exhibit opposing patterns of splicing (exon inclusion vs. exon skipping). A combination of unbiased genetic screens and candidate targeted mutations identified a cascade of three cell-fate determining TFs and two neuronal RBPs required for proper splicing of sad-1 in ALM neurons. Mechanistic dissection revealed that the three TFs function to establish cell-specific expression of the two RBPs in the ALM neuron, and that the two RBPs in turn directly bind to sad-1 intronic regions to mediate exon inclusion in the ALM neuron. Finally, we find that in other neuron types, similar principles apply but with different combinations of TFs and RBPs mediating sad-1 exon inclusion. These results indicate that neuronal RBPs, like TFs, are employed in a combinatorial code to shape neuron-specific splicing patterns, and demonstrate phenotypic convergence by which different RBPs mediate similar splicing outcomes in various neurons.

Results

Alternative splicing of the neuronal kinase sad-1 in specific cell types

To identify alternative splicing regulation in individual neuronal cell types, we created two-color splicing reporters that provide a fluorescent readout of splicing regulation in vivo in single cells (Kuroyanagi et al., 2006; Orengo et al., 2006). A minigene representing an alternative splicing event of interest is cloned upstream of a dual GFP/RFP cassette (Figure 1A–B). The GFP and RFP coding sequences reside in alternative reading frames. The alternative exon is engineered to shift the reading frame by +1 nucleotide such that splicing of the alternative exon determines the reading frame, and therefore the translation of GFP versus RFP. Application of two-color fluorescent reporters to transparent organisms such as C. elegans enables in vivo imaging of alternative splicing in individual cells. We have created reporters for splicing events in a number of neuronal genes, and uncovered a rich variety of splicing patterns in single neurons (Norris et al., 2014).

Figure 1 with 1 supplement see all
sad-1 is alternatively spliced in single neurons.

(A) The sad-1 gene. Alternative cassette exon in blue. (B) Two-color splicing reporter schematic for sad-1 cassette exon. The cassette exon encodes a + 1 nt frameshift so that when skipped, GFP is produced with an in frame stop codon. When skipped, GFP is read out of frame without stop codons, followed by in-frame translation of RFP. (C) Whole worm fluorescent micrograph demonstrating both exon inclusion (RFP) and skipping (GFP) in many neurons, while certain neurons express only the included (ALM) of skipped (BDU) isoforms. (D–E) Higher magnification focusing on ALM and BDU neurons. (F) BDU and ALM are both paired neurons present on the left and right side of the worm. Each BDU neuron is a sister cell to an ALM neuron, derived from the same neuroblast. Scale bar represents 10 µm.

https://doi.org/10.7554/eLife.46726.003

One intriguing example of neuron-specific alternative splicing is in the conserved neuronal kinase sad-1, which plays important roles in neuronal development in both worms and mice (Kim et al., 2008; Kishi et al., 2005). In C. elegans, sad-1 is encoded by seventeen exons, and the fifteenth exon is an alternative cassette-type exon (Figure 1A). Alternative splicing of this exon changes the coding sequence and length of the sad-1 C-terminus (Kim et al., 2010). This presents an interesting parallel with mice and human genomes, which encode two separate genes homologous to sad-1 (SAD-A and SAD-B) that are nearly identical except for their C-terminal coding sequence and length.

A two-color splicing reporter for sad-1 in C. elegans revealed that many neurons express both the skipped and included isoforms (Figure 1C, Figure 1—figure supplement 1). For example, motor neurons in the ventral nerve cord express both isoforms of sad-1 (Figure 1C). On the other hand, the ALM touch-sensing neuron expresses only the included isoform, while the BDU neuron, which is the sister cell to the ALM neuron, expresses only the skipped isoform (Figure 1C–F). While different neurons exhibit differences in sad-1 splicing, the splicing pattern in a given neuron is reproducible and invariant from one animal to the next, suggesting that sad-1 splicing in various neurons is under strict regulatory control. These results led us to ask how ALM and BDU neurons, which are developmentally related (Figure 1F) and share a number of anatomical and gene-expression features, specify opposite splicing regimes.

Forward genetic screen identifies a trio of fate-determining TFs affecting sad-1 splicing in the ALM neuron

To identify regulators of sad-1 splicing in the ALM touch neuron, we performed an unbiased forward genetic screen. Parental worms harboring the sad-1 splicing reporter were mutagenized with EMS. We then screened for F2 animals (potential homozygotes) with aberrant expression of the skipped (GFP) isoform in the ALM neuron (Figure 2A). This screen identified three distinct loci that transform the splicing pattern from the ALM neuron pattern (full exon inclusion) to resemble the pattern in their BDU sister cells (full exon skipping).

Figure 2 with 1 supplement see all
Genetic screen identifies neuronal TFs affecting sad-1 splicing in the ALM neuron.

(A) Schematic of forward genetic screen to identify regulators of sad-1 splicing in the ALM touch neuron. (B–F) ALM neurons (dashed boxes) shift from complete inclusion (RFP) to skipping (GFP) in unc-86(e1416), mec-3(e1338), or alr-1(oy42) TF mutants. Splicing phenotypes fully penetrant (n = 50 animals) (G) Previously-identified roles of the three TFs in a transcriptional cascade to control touch neuron gene expression. Scale bar represents 10 µm.

https://doi.org/10.7554/eLife.46726.005

Whole-genome resequencing of the mutant strains identified loss-of-function mutations in three conserved TFs: unc-86, mec-3, and alr-1 (Figure 2B–F, Figure 2—figure supplement 1). All three genes have previously been identified as key regulators of touch-neuron cell fate (Gordon and Hobert, 2015; Topalidou et al., 2011). The three TFs function in a transcriptional cascade ensuring cell-specific expression of mec-3 in touch neurons, which then results in expression of a battery of touch-neuron specific genes (Figure 2G). Loss of the TF mec-3 results in touch neurons (ALMs) adopting certain gene-expression characteristics of their sister cells (BDUs) (Gordon and Hobert, 2015), mirroring our observation that loss of mec-3 transforms sad-1 splicing from an ALM (exon 15 included) to a BDU (exon 15 skipped) pattern.

Previous work demonstrates that the MEC-3 TF is expressed only in touch neurons, while UNC-86 and ALR-1 are expressed in various neuron types (Topalidou et al., 2011). However, we find that unc-86 and alr-1 mutants affect sad-1 splicing only in the touch neurons (Figure 2D–F). This is in accordance with previous work indicating that a major function of unc-86 and alr-1 in touch neurons is to combinatorially ensure appropriate expression of mec-3, and that all three TFs are needed for proper differentiation of touch neurons (Topalidou et al., 2011). We therefore conclude that the combinatorial activity of all three TFs is required for proper sad-1 splicing in the ALM neuron.

A pair of RNA binding proteins regulates sad-1 splicing in the ALM neuron

We were surprised to identify TFs, but not RBPs, in our forward genetic screen for regulators of sad-1 alternative splicing. We hypothesized that multiple RBPs might co-regulate sad-1 alternative splicing in the ALM neuron and therefore mutations in individual RBPs might result in mild splicing defects. We therefore examined the sequence surrounding the sad-1 alternative exon for conserved cis-elements corresponding to known in vitro RBP sequence preferences (Ray et al., 2013). We identified three candidate elements: one corresponding to the mbl-1/Mbnl1 consensus binding motif, and two corresponding to the mec-8/RBMS motif (Figure 3A–C).

Figure 3 with 2 supplements see all
Two neuronal RBPs combinatorially control sad-1 splicing in ALM neurons.

(A) Conservation scores in the introns surrounding sad-1 exon 15, basewise phyloP26way comparison of 26 nematode species (Hubisz et al., 2011). Numbers 1–3 indicate consensus binding motifs for mbl-1 and mec-8 displayed in B-C. (B–C) cis-elements matching consensus binding motifs for mbl-1 and mec-8. (D–F) mec-8 and mbl-1 mutants both cause a partial loss of sad-1 exon inclusion. (G) mec-8; mbl-1 double mutants cause complete loss of exon inclusion, phenocopying the TF mutants. Splicing phenotypes fully penetrant (n = 50 animals) Scale bar represents 10 µm.

https://doi.org/10.7554/eLife.46726.007

To test whether these RBPs affect sad-1 alternative splicing, we created deletions for each gene with CRISPR/Cas9 (Norris et al., 2017). Both mec-8 and mbl-1 mutants result in aberrant sad-1 splicing in the ALM neuron, displaying partial skipping and partial inclusion (Figure 3D–F, Figure 3—figure supplement 1). As in the case of the TF mutants, mec-8 mutants affect sad-1 splicing specifically in the ALM neurons, whereas mbl-1 mutants affect sad-1 splicing in ALM neurons as well as specific neurons in the ventral nerve cord (see Figure 6, below). To verify that the phenotypes of our CRISPR mutants were on-target effects, we crossed the sad-1 splicing reporter into existing alleles for mec-8 (e398, premature stop codon [Davies et al., 1999; Lundquist et al., 1996]) and mbl-1 (wy560, large deletion affecting multiple genes including mbl-1 [Spilker et al., 2012]). We found these alleles to affect splicing of sad-1 exactly as our CRISPR mutations (Figure 3—figure supplements 12).

Whereas TF mutants result in full skipping of the sad-1 alternative exon, RBP mutants result in only partial skipping. This provides a probable explanation for not identifying these RBPs in our genetic screen: partial exon skipping leads to dim GFP expression, which is not sufficiently bright to be noticed upon brief visual inspection. We therefore tested whether simultaneous loss of both RBPs recapitulates the full skipping of sad-1 exon 15 observed in TF mutants. We created mec-8; mbl-1 double mutants expressing the sad-1 splicing reporter. These double mutants result in complete loss of sad-1 exon inclusion in the ALM neuron, recapitulating the splicing phenotype of the single TF mutants (Figure 3G). These results led us to hypothesize that the TFs identified in our screen exert their effects on sad-1 splicing by controlling expression of both mec-8 and mbl-1.

TFs affecting sad-1 splicing are required for RBP expression in the ALM neuron

To examine whether the neuronal TFs alter expression of mec-8 and mbl-1 RBPs in the ALM neuron, we created reporter lines for each RBP. To this end, each RBP was C-terminally tagged in a fosmid containing large regions of surrounding genomic context (Poser et al., 2008; Spilker et al., 2012) (Figure 4A–E). Compared to traditional transgenic reporters, fosmids are more likely to contain all regulatory information needed to drive normal expression of the gene in question. This is demonstrated in the case of the mec-8 RBP. The classical mec-8::GFP promoter fusion drives expression in a number of cells, but not in the ALM neuron (Figure 4—figure supplement 1) (Spike et al., 2002). On the other hand, we detected expression of the mec-8 fosmid reporter in many of the same cells, both neuronal and non-neuronal, plus strong expression in the ALM neuron (Figure 4A–B). A similar fosmid reporter for mbl-1 likewise exhibits expression in the ALM neuron, as well as many other neurons in the nervous system (Figure 4D, Figure 4—figure supplement 1). This is in line with previous reports on mbl-1 expression (Spilker et al., 2012).

Figure 4 with 3 supplements see all
Neuronal TFs establish expression of both mec-8 and mbl-1 to mediate splicing of sad-1 in ALM neurons.

(A–B) A mec-8 translational GFP fosmid reporter reveals strong expression in ALM neuron (strong expression in 28/31 = 90% of animals inspected). (C) In a mec-3 TF mutant, mec-8 expression is absent specifically in ALM (no detectable expression in 43/50 = 86%, dim expression in 7/50 = 14% of animals inspected). (D) mbl-1 translational RFP fosmid reporter is expressed in ALM neuron (strong expression in 19/20 = 95% of animals inspected). (E) In a mec-3 mutant, mbl-1 expression is absent specifically in ALM (no detectable expression in 19/21 = 90%, dim expression in 2/21 = 10% of animals inspected). (F–G) Aberrant splicing of sad-1 in alr-1 TF mutants is partially rescued by over-expression of either mec-8 (6/6 animals examined) or mbl-1 (6/7 animals examined) RBPs (H–I). Scale bar represents 10 µm.

https://doi.org/10.7554/eLife.46726.010

We tested expression of our reporters in the context of a mec-3 mutant to determine whether expression of mec-8 and mbl-1 in ALM neurons depends on the TF cascade uncovered in our screen. The mec-3 TF is expressed only in touch neurons, and therefore we would expect mec-3 mutants to affect RBP expression only in the touch neurons. Indeed, in mec-3 mutants, expression of both mec-8 and mbl-1 RBPs are abolished in the ALM neuron, while expression in the surrounding neurons and tissues remains unchanged (Figure 4B–E). Together these results indicate that the expression of mec-8 and mbl-1 RBPs are under the control of neuron subtype-specific TFs.

To examine whether mec-8 and mbl-1 RBPs might be under direct transcriptional control by one or more of the TFs, we used existing ChIP data for ALR-1 (Niu et al., 2011), in vitro derived consensus binding motifs for UNC-86 (Weirauch et al., 2014), and a previously-defined UNC-86/MEC-3 heterodimer binding motif (Röhrig et al., 2000; Xue et al., 1993). We did not find conserved UNC-86 binding motifs or an UNC-86/MEC-3 heterodimer binding motif in the promoters for mec-8 or mbl-1, but did find ALR-1 ChIP peaks in both promoters (Figure 4—figure supplement 2). This data suggests that alr-1 may directly control transcription of mec-8 and mbl-1 RBPs.

TFs affect sad-1 splicing by controlling RBP expression in the ALM neuron

The observations that (1) mec-8; mbl-1 RBP double mutants recapitulate the phenotype of the TF mutants, and (2) the TFs are necessary for expression of both RBPs in the ALM neuron, together suggest that the splicing defects in the TF mutants are mediated by effects on expression of the two RBPs. Further support for this hypothesis arose indirectly in the course of crossing TF and RBP mutants together. We found that while TF or RBP mutant heterozygotes exhibit normal sad-1 splicing in the ALM neuron, double heterozygotes (for example alr-1/+; mbl-1/+, or mec-3/+; mec-8/+) exhibit partial exon skipping in ALM, similar to the RBP single mutants (Figure 4—figure supplement 3). Such ‘non-allelic non-complementation’ is often interpreted to mean that the two genes function in the same complex, or, more likely in this case, function in the same pathway (Yook et al., 2001). This indirect evidence further suggests that the TFs and RBPs affect sad-1 splicing as part of the same molecular pathway.

If sad-1 splicing is controlled in a linear pathway as suggested by the above series of experiments, with upstream TFs affecting RBP expression in the ALM neuron, then over-expressing an RBP in the context of a TF mutant should partially restore splicing in ALM. To test this hypothesis we created a strain over-expressing a mec-8 transgene specifically in the touch neurons (pmec-3::mec-8). When introduced into an alr-1 mutant, this transgene partially rescues the splicing of sad-1 in the ALM neuron (Figure 4F–H). Likewise, over-expression of mbl-1 in an alr-1 mutant partially rescues splicing in the ALM neuron (Figure 4I). These results further support a linear gene regulatory pathway in which neuronal fate-determining TFs control neuron-specific expression of RBPs, which then control alternative splicing of sad-1 (Figure 4J).

RBPs directly mediate sad-1 exon inclusion through interactions with surrounding introns

To test whether mec-8 and mbl-1 directly affect splicing by binding to the sad-1 pre-mRNA, we created two-color splicing reporters in which the putative mec-8 or mbl-1 cis-elements are mutated (Figure 3A and Figure 5). If the RBPs act directly by binding the cis-element, then mutation of the cis-element should affect the splicing pattern in a manner resembling the wild-type splicing reporter in the context of the RBP deletion mutant. If the RBPs act indirectly, mutating the cis-element should have no effect on the splicing pattern.

Figure 5 with 1 supplement see all
mbl-1 and mec-8 affect sad-1 splicing by direct interaction with sad-1 introns.

(A–B) Mutation of mbl-1 consensus sequence in sad-1 splicing reporter results in aberrant splicing in ALM neurons that phenocopies an mbl-1 mutant. (C–E) Mutation of either mec-8 binding motif, or both simultaneously, likewise results in aberrant sad-1 splicing in ALM neurons. ALM splicing phenotypes fully penetrant (n = 25 animals) Scale bar represents 10 µm.

https://doi.org/10.7554/eLife.46726.014

Mutation of the mbl-1 cis-element resulted in ALM neurons with altered sad-1 splicing in which the exon is partially skipped and partially included (Figure 5A–B). This recapitulates the phenotype of mbl-1 null mutations (Figure 3F), suggesting that mbl-1 exerts its effects on splicing directly through binding a conserved cis-element in the upstream intron.

We identified two consensus mec-8 binding motifs in conserved regions in the intron downstream of the cassette exon. We therefore created splicing reporters mutant for both cis-elements as well as for each element individually. The splicing reporter mutant for both elements recapitulates the splicing phenotype of mec-8 null mutants (Figure 5E). Likewise, mutating either mec-8 binding site in isolation recapitulates a mec-8 null mutation (Figure 3E and Figure 5C–D), suggesting that mec-8 binding to both cis-elements is required for appropriate sad-1 splicing.

We tested whether mutation of a putative cis-element could be rescued by over-expression of its cognate RBP, and found that cis-element mutants were not rescued by RBP over-expression (Figure 5—figure supplement 1), providing further evidence that the RBPs act directly on the sad-1 pre-mRNA. Together these results indicate that mec-8 and mbl-1 RBPs combinatorially ensure sad-1 exon inclusion in ALM neurons through direct interactions with the neighboring introns.

sad-1 splicing in other neuron types is controlled through both distinct and overlapping mechanisms

Having identified regulatory mechanisms controlling sad-1 splicing in the ALM neuron, we next wondered whether similar principles apply in other neuron types. Most neurons besides the ALM and BDU neurons express both skipped and included sad-1 isoforms. This could represent the neuronal ‘ground state’ of splicing in the absence of cell-specific splicing regulators. On the other hand, our observations that loss of both mec-8 and mbl-1 in the ALM neuron results in full exon skipping suggest that the ground state may be complete exon skipping. This hypothesis predicts that other neurons in which sad-1 is partially included express one or more RBPs mediating exon inclusion.

In the course of examining sad-1 splicing in ALM neurons, we noticed that mbl-1 mutants affect sad-1 splicing not only in ALM, but also in the excitatory cholinergic motor neurons of the ventral nerve cord (Figure 6A–D). Whereas mbl-1 mutants cause a change in sad-1 splicing from full inclusion to partial inclusion in ALM neurons, in excitatory motor neurons mbl-1 mutants shift from partial inclusion to no inclusion (Figure 6C–D). On the other hand, the inhibitory motor neurons remain unaffected in mbl-1 mutants, expressing both the included and skipped isoforms (Figure 6D, arrowheads). This is consistent with our mbl-1 gene expression reporter, which reveals expression of mbl-1 in the excitatory motor neurons, but not in the inhibitory motor neurons (Figure 6—figure supplement 1).

Figure 6 with 2 supplements see all
sad-1 splicing in motor neurons of the ventral nerve cord is controlled by mbl-1 and msi-1 RBPs.

(A–C) In wild-type worms, sad-1 is partially included in both excitatory and inhibitory motor neurons. (D) In mbl-1 mutants, exon inclusion is lost in excitatory motor neurons, but remains in inhibitory motor neurons (arrowheads). (E) msi-1 mutants lose exon inclusion in inhibitory motor neurons (arrowheads) but not in excitatory motor neurons. (F) mbl-1; msi-1 double mutants lose exon inclusion in all motor neurons in the ventral nerve cord. Splicing phenotypes in ventral nerve cord invariant (n = 15 animals) (G) Conservation scores (determined as in Figure 3A) in the introns surrounding sad-1 exon 15. Number one indicates consensus binding motifs for msi-1. (H) cis-elements matching consensus binding motifs for msi-1. Asterisk indicates anterior-posterior position of ALM neuron as anatomical reference. Splicing phenotypes fully penetrant (n = 50 animals). Scale bar represents 10 µm.

https://doi.org/10.7554/eLife.46726.016

We did not detect mec-8 expression in motor neurons of the ventral nerve cord, and mec-8 mutants had no effect on splicing of sad-1 in motor neurons (Figures 3E and 4B). It therefore seems that in neurons expressing mbl-1 such as excitatory motor neurons, the presence of mbl-1 mediates partial exon inclusion. In neurons expressing both mbl-1 and mec-8 such as ALM touch neurons, the two RBPs together mediate full inclusion.

In mbl-1 mutants, sad-1 exon inclusion is lost in excitatory neurons but remains in inhibitory motor neurons. We therefore wondered whether there was an additional RBP expressed in inhibitory motor neurons mediating sad-1 inclusion. mec-8 was ruled out because it is not expressed in inhibitory motor neurons and does not affect sad-1 splicing in the nerve cord. On the other hand, the RBP msi-1/Musashi has been reported to be expressed in inhibitory but not excitatory neurons of the nerve cord (Yoda et al., 2000), which is a mutually exclusive pattern with mbl-1. We therefore tested msi-1 as a candidate for the RBP mediating sad-1 exon inclusion in the inhibitory motor neurons. We generated a msi-1 deletion mutant, which shows loss of sad-1 inclusion specifically in the inhibitory motor neurons (Figure 6E). Furthermore, msi-1; mbl-1 double mutants result in complete loss of exon inclusion in the ventral nerve cord (Figure 6F). These results indicate that mbl-1 and msi-1 act in distinct cell types to achieve partial sad-1 exon inclusion throughout the ventral nerve cord.

We suspect that msi-1, like mbl-1 and mec-8, directly affects sad-1 splicing by binding in the intronic regions surrounding the alternative exon. in vitro experiments have identified a UAG motif (Figure 6H) (Ray et al., 2013), usually in bipartite form (e.g. UAGNNUAG) (Dominguez et al., 2018), as the consensus binding motif for msi-1. There is a conserved bipartite UAG motif in the intron downstream of the sad-1 cassette exon (Figure 6G–H), and we hypothesize that msi-1 binds there to mediate exon inclusion in inhibitory motor neurons.

Phenotypic convergence of splicing regulation in diverse neuron types

Together the results from three different neuronal cell types (ALM neuron, excitatory motor neurons, and inhibitory motor neurons) constitute an example of ‘phenotypic convergence,’ in which phenotypic similarity between cells is generated by distinct molecular mechanisms. Substantial evidence of such phenotypic convergence exists for TFs controlling neuronal properties in worms and flies (Gendrel et al., 2016; Konstantinides et al., 2018; Pereira et al., 2015). Our results now extend this principle to RBPs and their control of alternative splicing, revealing phenotypic convergence in which similar splicing patterns (i.e. sad-1 exon inclusion) are generated in various neurons by diverse RBPs acting in specific neuronal subtypes (Figure 7).

Phenotypic convergence at the level of splicing regulation.

Different RBPs act in different neuron types to carry out the common function of mediating sad-1 exon inclusion.

https://doi.org/10.7554/eLife.46726.019

To further examine this principle, we tested whether ectopic expression of an RBP in a neuron type in which it is not normally expressed would be sufficient to alter sad-1 splicing in that neuron. We expressed mec-8 in excitatory motor neurons (where normally only mbl-1 is expressed) and found that mec-8 expression is sufficient to alter sad-1 splicing patterns from partial inclusion to full inclusion specifically in the excitatory motor neurons (Figure 6—figure supplement 2A). Similarly, mbl-1 expression in inhibitory motor neurons (where normally only msi-1 is expressed) results in full exon inclusion [Figure 6—figure supplement 2A]).

Finally, we asked whether phenotypic convergence occurs simultaneously at multiple levels (TFs and RBPs) with regard to sad-1 splicing. To do so we examined mutants for the TF unc-3, which controls the fate of excitatory motor neurons in the ventral nerve cord (Kratsios et al., 2012), analogous to ALM cell fate determination by unc-86/mec-3/alr-1. In unc-3 mutants, sad-1 exon inclusion is lost in excitatory motor neurons, similar to mbl-1 RBP mutants (Figure 6D, Figure 6—figure supplement 2B). However, whereas unc-86/mec-3/alr-1 mutants exhibit completely-penetrant loss of sad-1 exon inclusion, unc-3 mutants exhibit partially-penetrant defects, ranging from moderate to complete loss of sad-1 inclusion in excitatory motor neurons. Similarly, loss of unc-3 results in partially-penetrant defects in mbl-1 expression (Figure 6—figure supplement 2C).

Together these results demonstrate that phenotypic convergence among different neuron types occurs simultaneously at multiple layers of gene regulation: different TFs (e.g. mec-3 and unc-3) specify expression of different RBP complements (e.g. mbl-1 and mec-8) which have a common function of mediating sad-1 exon inclusion.

Discussion

Neuron-specific regulation of sad-1 splicing

In this study we find that sad-1 splicing undergoes precise regulation in numerous neuronal types. Although ALM and BDU neurons are sister cells, express many of the same genes, and share a number of cell-specific TFs, they have opposing patterns of sad-1 splicing. This highlights the fact that post-transcriptional control can further diversify attributes of single cells on top of the more well-known role of transcriptional control.

Our results demonstrate that sad-1 splicing is regulated according to a combinatorial RBP code, with different splicing outcomes depending on whether a cell expresses zero, one, or two neuron-specific RBPs (Figure 7).This suggests that the ‘default’ outcome of sad-1 splicing is full skipping of the cassette exon, as observed in the BDU neuron which does not express any of the sad-1-regulating RBPs. Only cells with at least one RBP mediating exon inclusion express sad-1 included isoforms. Cells with multiple such RBPs (e.g. the ALM neuron) express only the included isoform.

In previous work we found that alternative splicing of the kinase unc-16/ JIP3 in motor neurons is likewise controlled by a pair of RNA binding proteins (Norris et al., 2014). However, unc-16 splicing and sad-1 splicing in motor neurons are regulated by distinct pairs of RBPs. Whereas sad-1 splicing in motor neurons is regulated by mbl-1 and msi-1 RBPs, unc-16 is combinatorially regulated by unc-75 and exc-7 in motor neurons (Norris et al., 2014). This suggests that even within a single neuron type, different splicing events are regulated by different complements of RBPs.

Coordinated splicing regulation across layers of gene expression

The importance of TFs controlling gene expression networks in single neurons is well established, and the importance of RBPs controlling post-transcriptional networks in single cells is gaining wider appreciation (Norris and Calarco, 2012; Norris et al., 2014; Song et al., 2017; Wamsley et al., 2018). How these two modes of regulation might interact remains understudied. Here we show that the two modes of regulation interact in a traditional linear type of pathway. A combination of cell-specific TFs establishes a transcriptional network in a single neuron type. This network includes a specific combination of neuronal RBPs, and the particular combination of RBPs in a given neuron then establish a unique post-transcriptional gene regulatory network in that neuron. Multiple layers of regulatory control can thus increase the diversity of single neuron transcriptomes and fine-tune the properties of individual neurons.

In the present study we have identified a linear pathway in which TFs influence the expression of RBPs, which then influence alternative splicing in single neurons. This adds to a substantial body of literature finding that RBPs can affect the function of specific TFs by modulating their alternative splicing (Calarco et al., 2009; Han et al., 2013; Linares et al., 2015; Raj et al., 2011). In the future it will be interesting to see whether additional regulatory logics exist between TFs and RBPs. Single-neuron TF combinations have been identified with a variety of feedback and feedforward mechanisms resulting in interesting regulatory properties (Mangan and Alon, 2003), and in principle TFs and RBPs could likewise interact in complex ways, leading to an even greater array of diversification strategies (Han et al., 2017). Together this study highlights the importance of considering neuron-specific ‘combinatorial codes’ not only from the perspective of TF combinations, but the specific complement of both TFs and RBPs shaping the transcriptome of a given neuron.

Phenotypic convergence at the level of single-neuron splicing

A theme emerging from recent studies of single-neuron transcriptomes is ‘phenotypic convergence,’ in which multiple neurons share gene expression similarities, but the regulatory mechanisms by which they do so are distinct in each neuron. For example, in worms, cholinergic neuron cell fate and core cholinergic gene expression properties are controlled by different combinations of TFs in different cholinergic neuron sub-types (Pereira et al., 2015). This is also the case for other neuron types in C. elegans (Gendrel et al., 2016). More recently, phenotypic convergence has been reported for TFs in neurons of the Drosophila optic lobe (Konstantinides et al., 2018), indicating that phenotypic convergence mediated by TFs is a widespread phenomenon.

We now extend this principle of phenotypic convergence to the regulation of splicing by RBPs as well. sad-1 exon inclusion is mediated in various neuron types, with a unique complement of RBPs responsible for exon inclusion in each specific type that we have studied (ALM neuron, inhibitory motor neurons, and excitatory motor neurons; Figure 7). This likely represents phenotypic convergence on multiple levels, as the RBPs regulating splicing are different in each neuron, and the TFs regulating RBP expression are likewise different in each neuron. Each of these levels coordinately converges upon appropriate splicing of sad-1 in each neuron type. Additional neuron types with similar sad-1 splicing patterns (see Figure 1 and Figure 1—figure supplement 1) may represent additional examples of phenotypic convergence whose underlying mechanisms remain unexplored.

Materials and methods

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional
information
Strainunc-86(csb9)This studyJAC401Norris Lab. SMU. Dallas, TX.
Strainmec-3(csb10)This studyJAC402Norris Lab. SMU. Dallas, TX.
Strainalr-1(csb11)This studyJAC403Norris Lab. SMU. Dallas, TX.
Strainunc-86(e1416)CGC, University of MinnesotaCB1416
Strainmec-3(e1338)CGC, University of MinnesotaCB1338
Strainalr-1(oy42)CGC, University of MinnesotaPY1598
Strainmec-8(e398)CGC, University of MinnesotaCB398
Strainmec-8(csb22)This studdyJAC626Norris Lab. SMU. Dallas, TX.
Strainmbl-1(csb31)This studyJAC635Norris Lab. SMU. Dallas, TX.
Strainmbl-1(wy560)CGC, University of MinnesotaJAC002
Strainmsi-1(csb24)This studyJAC628Norris Lab. SMU. Dallas, TX.
Strainmec-8(csb22); mbl-1(wy560)This studyADN342Norris Lab. SMU. Dallas, TX.
Strainmec-8(csb22); mbl-1(csb31)This studyJAC670Norris Lab. SMU. Dallas, TX.
Strainmsi-1(csb24); mbl-1(csb31)This studyADN257Norris Lab. SMU. Dallas, TX.
Strainpmec-3::mec-8This studyADN431Norris Lab. SMU. Dallas, TX.
Strainpmec-3::mbl-1This studyADN514Norris Lab. SMU. Dallas, TX.
Strainpunc-25::mbl-1This studyADN515Norris Lab. SMU. Dallas, TX.
Strainpunc-17::mec-8This studyADN505Norris Lab. SMU. Dallas, TX.
StrainΔmbl-1 cis-element sad-1 splicing reporterThis studyADN319Norris Lab. SMU. Dallas, TX.
StrainΔmec-8[1] cis-element sad-1 splicing reporterThis studyADN364Norris Lab. SMU. Dallas, TX.
StrainΔmec-8[2] cis-element sad-1 splicing reporterThis studyADN377Norris Lab. SMU. Dallas, TX.
StrainΔmec-8[both] cis-element sad-1 splicing reporterThis studyADN333Norris Lab. SMU. Dallas, TX.
Strainsad-1 splicing reporterThis studyJAC017Norris Lab. SMU. Dallas, TX.
StrainMEC-8::GFP reporter fosmidThis studyJAC447Norris Lab. SMU. Dallas, TX.
StrainMBL-1::RFP reporter fosmidThis studyJAC576Norris Lab. SMU. Dallas, TX.
Strainpmec-8::GFPCGC, University of MinnesotaBC11068
Strainunc-3(e151)CGC, University of MinnesotaCB151

Strain maintenance

Request a detailed protocol

C. elegans were maintained under standard conditions (Brenner, 1974) at 20°C on nematode growth media (NGM) plates seeded with OP50 E. coli bacteria. New transgenic worms were generated by microinjection with 15 ng/μl transgene and 15 ng/μl co-injection marker (either rgef-1, unc-17, or unc-25 promoter driving BFP).

Mutant generation and genetic screening

Request a detailed protocol

The forward mutagenesis screen was performed on animals harboring the sad-1 exon 15 splicing reporter with EMS at 47 mM for 4 hr. F1s were picked onto new plates, 10 F1s per plate. After 3–4 days of growth, F2s were screened by eye on the Zeiss Axiozoom.V16 for touch cells appearing in the GFP channel (representing aberrant exon skipping) and were then verified for a concomitant loss of RFP (representing loss of exon inclusion). Such worms were picked individually onto a new plate to verify the phenotype in the F3 generation and to establish a clonal population. After outcrossing, strains were subjected to whole-genome resequencing (Illumina, 1 × 75 bp) and potential causative mutations were identified using the CloudMAP workflow on the Galaxy web platform (Minevich et al., 2012). A total of approximately 6000 haploid genomes were screened.

Targeted mutant strains were generated using CRISPR/Cas9 as previously described (Calarco and Norris, 2018; Norris et al., 2015), such that the gene of interest is deleted and is replaced with a heterologous GFP reporter under the control of a pharyngeal promoter (pmyo-2) which does not interfere with the visualization of the sad-1 splicing reporter in the ALM, BDU or ventral nerve cord neurons. Seamless gene replacement was verified by PCR amplification and Sanger sequencing of both junction boundaries.

Microscopy

Request a detailed protocol

Images were obtained with a Zeiss Axio Imager.Z1 and processed in ImageJ.

Generation of splicing reporters

Request a detailed protocol

sad-1 minigenes were created using the following primers: Forward 5’ GATAAAACTGAAACAACTTCTGC and Reverse 5’ GGGGTTGGCGATTTGTATGAGaTAGC. Restriction sites were appended to both the forward primer (XhoI) and reverse (NotI) primers to facilitate cloning into a Gateway-compatible vector as previously described (Norris et al., 2014). The reporter was then cloned downstream of a pan-neuronal rgef-1 promoter, as endogenous sad-1 has been detected broadly throughout the nervous system (Crump et al., 2001). Mutant versions of the splicing reporter were synthesized de novo then cut with XhoI and NotI and cloned as above.

Some strains were provided by the Caenorhabditis Genome Center, which is funded by the NIH Office of Research Infrastructure Programs (P40 OD010440). Other strains were provided by the National BioResource Project (Tokyo).

Data availability

All data from this study included in manuscript and supplemental materials.

References

    1. Brenner S
    (1974)
    The genetics of Caenorhabditis elegans
    Genetics 77:71–94.
    1. Davies AG
    2. Spike CA
    3. Shaw JE
    4. Herman RK
    (1999)
    Functional overlap between the mec-8 gene and five sym genes in caenorhabditis elegans
    Genetics 153:117–134.
  1. Book
    1. Gracida X
    2. Norris AD
    3. Calarco JA
    (2016)
    Regulation of Tissue-Specific Alternative Splicing: C. elegans as a Model System
    In: Yeo G. W, editors. RNA Processing. Cham: Springer International Publishing. pp. 229–261.
    1. Lundquist EA
    2. Herman RK
    3. Rogalski TM
    4. Mullen GP
    5. Moerman DG
    6. Shaw JE
    (1996)
    The mec-8 gene of C. Elegans encodes a protein with two RNA recognition motifs and regulates alternative splicing of unc-52 transcripts
    Development 122:1601–1610.
    1. Spike CA
    2. Davies AG
    3. Shaw JE
    4. Herman RK
    (2002)
    MEC-8 regulates alternative splicing of unc-52 transcripts in C. elegans hypodermal cells
    Development 129:4999–5008.
    1. Yook KJ
    2. Proulx SR
    3. Jorgensen EM
    (2001)
    Rules of nonallelic noncomplementation at the synapse in Caenorhabditis elegans
    Genetics 158:209–220.

Decision letter

  1. Oliver Hobert
    Reviewing Editor; Howard Hughes Medical Institute, Columbia University, United States
  2. James L Manley
    Senior Editor; Columbia University, United States
  3. Hidehito Kuroyanagi
    Reviewer; Tokyo Medical and Dental University, Japan

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

Thank you for submitting your article "Splicing in a single neuron is coordinately controlled by RNA binding proteins and transcription factors" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and James Manley as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Hidehito Kuroyanagi (Reviewer #2).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Summary:

In this paper, the authors use an in vivo fluorescent splicing reporter to study cell-specific post-translational regulation. Making use of the genetic amenability of C. elegans, the authors employ forward genetic screens and candidate approaches to identify the genes required for cell-specific splicing events. Through the characterization of the complex splicing patterns of the conserved neuronal kinase sad-1, the authors elegantly reveal that neuronal identity regulator transcription factors and RNA binding proteins are required for sad-1 exon inclusion. Neuron-type specific expression of RBPs is controlled by these transcription factors, and RBPs directly interact with sad-1 introns. Moreover, the authors demonstrate that sad-1 splicing in different neuron-types is controlled by different RBPs, revealing phenotypic convergence at the level of single-neuron splicing.

The authors employed CRISPR-Cas9 mediated genome editing and a smart splicing reporter to support their conclusions. The overall conclusions on transcriptional control of RBPs and phenotypic convergence appear interesting. However, there is some lack in the depth of analysis, particularly as it pertains to the analysis of the transcription factors. Moreover, there is striking, in fact almost baffling, lack of detail in data/reagent presentation throughout the manuscript. A substantial revision of the paper will allow the authors to better support what appears to be a very interesting finding.

Essential revisions:

1) There is a concern that the evidence that the three RBPs "directly" regulate pre-mRNA splicing of the sad-1 gene via their respective cis-elements is weak. The conclusion is based on high-throughput in vitro screening for consensus binding sequences of the RBPs (Ray et al., 2013), effects of the RBP mutations on the sad-1 reporter expression (Figure 3 and Figure 6) and effects of mutations in the putative cis-elements on the reporter expression (Figure 5). Direct and specific binding of recombinant RBPs to the respective RNA sequences in vitro and/or defective responses of the mutant reporters to RBP overexpression need to be provided support the direct regulation; otherwise the possibility that MBL-1 regulating sad-1 splicing through post-transcriptional regulation of MEC-8 and/or other RBP(s), and vice versa, cannot be excluded. This point is important especially for mbl-1 and msi-1 because the manuscript will be the first experimental demonstration that these genes regulate pre-mRNA splicing in C. elegans. It is also important to keep in mind that the consensus sequences experimentally identified in Ray et al., 2013 and logos shown in Figure 3, Figure 5 and Figure 6 of the manuscript are for a human orthologue of MBL-1, human and Drosophila orthologues of MSI-1 and not those in C. elegans (this should be more explicitly stated in the figures).

2) The common standard for validating that a mutant retrieved from a screen really is the causative mutant is to either rescue or use multiple alleles. The authors do neither. They need to validate the identity of their screen mutants by examining the expression of the sad-1 splicing reporter on a canonical allele of at least mec-3.

3) It would be appreciated if the authors could extend their concept of TFs controlling splicing patterns by RBPs through testing whether mbl-1 expression in excitatory neurons – and hence sad-1 splicing – depends on unc-3, the excitatory MN "alter ego" of unc-86/mec-3 in the touch neurons. This experiment is very easy to perform and would broaden the concept of cell type-specific TFs controlling cell-type specific RBPs to other neuron types.

4) Figure 7. The authors discuss that 100% exon inclusion might require multiple RBPs. Alternatively, 100% exon inclusion might require RBP binding on both sides of the exon in question (e.g. in ALM, mbl-1 seems to bind upstream of sad-1 exon 15, while mec-8 seems to bind downstream). This hypothesis can be easily tested by mec-8 overexpression in cholinergic MN (mirroring ALM neurons), or by mbl-1 overexpression in GABAergic MN (preferable since it is a new untested scenario).

5) Generally, there is a baffling, if not troubling lack of any detail in describing and properly quantifying data and of description of reagents. This pervades the entire manuscript. Below a number of examples (in no order of importance):

a) In panels Figure 2C-F, Figure 3D-G, Figure 4B-E, Figure 4F-I, Figure 5, Figure 6C-F, Figure 3—figure supplement 2, Figure 4—figure supplement 1, Figure 4—figure supplement 2, and Figure 4—figure supplement 3 the authors do not provide any quantification to support their data (incl. the number of animals examined). It is unclear whether in all cases there is a fully penetrant phenotype. In cases of a fully penetrant phenotype (e.g. Figure 2C-F?), a simple explanation on the legend accompanied with the number of animals where the phenotype was observed will suffice. It is odd that the authors do not appear to be familiar with the standards in the field. In the case of partial exon skipping phenotypes, it is unclear to this reviewer whether all animals show always same levels of GFP and RFP expression in ALM (e.g. in the double heterozygotes). For variable phenotypes, a qualitative quantification (%dark/dim/bright red/green), or any other means of quantification deemed appropriate by the authors, is required.

b) There's no description of the splicing reporter. Is the reporter driven by its own promoter or the endogenous promoter? They refer to a previous paper (Norris, 2014) which states vectors that drove expression of each reporter pan-neuronally, or with endogenous promoter sequences

c) The authors generate a number of deletion mutants for different genes using CRISPR/Cas9 (mec-8, mbl-1, not specified for msi-1) however, they don't provide any description on the nature of the molecular lesion other than a brief reference to a previous paper. A schematic of the mutations should be provided in the respective main figures or in a supplementary figure, with a more detailed explanation on the legends or in the Materials and methods section.

d) There is no list of the strains used on this study as a supplementary file. For example, which is the allele name of the new deletions?

e) Authors need to describe how new transgenic strains were generated (DNA mix concentrations, simple or complex arrays, etc.).

f) These newly generated deletion mutants via CRISPR/Cas9 seem to include a GFP reporter construct, subsection “Mutant generation and genetic screening”: "the gene of interest is deleted and is replaced with a heterologous GFP reporter construct". However, authors use fosmid reporters to analyze mec-8 and mbl-1 expression (Figure 4). Do these new mbl-1 and mec-8 deletion alleles recapitulate mbl-1 and mec-8 expression through the GFP reporter? In that case, does this GFP expression interfere with the analysis of the sad-1 splicing reporter green fluorescence?

g) Subsection “RBPs directly mediate sad-1 exon inclusion through interactions with surrounding introns”, "Most neurons besides the ALM and BDU neurons express both skipped and included sad-1 isoforms". Attending to the merge provided in Figure 1C, several neurons in the head and tail seem to exclude exon 15. It will be very informative if authors provide a more detailed expression data of their splicing reporter on Figure 1. For example, just a gross quantification of the% of neurons with included/excluded/both in the head, tail, major ganglia, etc., will greatly help the authors to better illustrate the splicing event under study.

h) Figure 4A-E, mec-8 and mbl-1 expression. Was mec-8 fosmid reporter expression detected in other neurons? Also, authors mention that the mbl-1 fosmid reporter is expressed in "most other neurons" (subsection “TFs affecting sad-1 splicing are required for RBP expression in the ALM neuron”, Figure 4D). To support their point, authors could show whole worm images or high magnification images of specific regions (head, tail, etc.). This needs to be properly documented with an image.

https://doi.org/10.7554/eLife.46726.022

Author response

Essential revisions:

1) There is a concern that the evidence that the three RBPs "directly" regulate pre-mRNA splicing of the sad-1 gene via their respective cis-elements is weak. The conclusion is based on high-throughput in vitro screening for consensus binding sequences of the RBPs (Ray et al., 2013), effects of the RBP mutations on the sad-1 reporter expression (Figure 3 and Figure 6) and effects of mutations in the putative cis-elements on the reporter expression (Figure 5). Direct and specific binding of recombinant RBPs to the respective RNA sequences in vitro and/or defective responses of the mutant reporters to RBP overexpression need to be provided support the direct regulation; otherwise the possibility that MBL-1 regulating sad-1 splicing through post-transcriptional regulation of MEC-8 and/or other RBP(s), and vice versa, can not be excluded. This point is important especially for mbl-1 and msi-1 because the manuscript will be the first experimental demonstration that these genes regulate pre-mRNA splicing in C. elegans. It is also important to keep in mind that the consensus sequences experimentally identified in Ray et al., 2013 and logos shown in Figure 3, Figure 5 and Figure 6 of the manuscript are for a human orthologue of MBL-1, human and Drosophila orthologues of MSI-1 and not those in C. elegans (this should be more explicitly stated in the figures).

We have performed the suggested in vivoexperiments, subjecting mutant reporters to RBP overexpression, and found that the defective splicing pattern remains unchanged (subsection “RBPs directly mediate sad-1 exon inclusion through interactions with surrounding introns” and Figure 5—figure supplement 1). We also now explicitly state in the figures which species provided each recombinant RBP for generation of in vitrobinding motifs.

2) The common standard for validating that a mutant retrieved from a screen really is the causative mutant is to either rescue or use multiple alleles. The authors do neither. They need to validate the identity of their screen mutants by examining the expression of the sad-1 splicing reporter on a canonical allele of at least mec-3.

We thank the reviewers for pointing out this oversight. We have remedied the error, displaying canonical alleles for each of the TFs identified in our screen in Figure 2 (with allele designations in legend), and alleles identified from our screen in Figure 2—figure supplement 1.

3) It would be appreciated if the authors could extend their concept of TFs controlling splicing patterns by RBPs through testing whether mbl-1 expression in excitatory neurons – and hence sad-1 splicing – depends on unc-3, the excitatory MN "alter ego" of unc-86/mec-3 in the touch neurons. This experiment is very easy to perform and would broaden the concept of cell type-specific TFs controlling cell-type specific RBPs to other neuron types.

These experiments have been performed, and support the general concept of TFs controlling alternative splicing via RBP expression (Figure 6—figure supplement 2). One unique aspect we have noted is that with unc-3 mutants and excitatory motor neurons the phenotypes seem to be less penetrant than the fully-penetrant defects in unc-86/mec-3/alr-1 mutants (discussed in the Discussion section).

4) Figure 7. The authors discuss that 100% exon inclusion might require multiple RBPs. Alternatively, 100% exon inclusion might require RBP binding on both sides of the exon in question (e.g. in ALM, mbl-1 seems to bind upstream of sad-1 exon 15, while mec-8 seems to bind downstream). This hypothesis can be easily tested by mec-8 overexpression in cholinergic MN (mirroring ALM neurons), or by mbl-1 overexpression in GABAergic MN (preferable since it is a new untested scenario).

Both experiments have been performed (Figure 6—figure supplement 2, and subsection “Phenotypic convergence of splicing regulation in diverse neuron types”), and indeed, ectopic expression of either mbl-1 or mec-8 in a neuron type in which they are not normally expressed is sufficient to alter sad-1 splicing patterns (in both cases changing the neuron’s splicing pattern from partial exon inclusion to ALM-like full inclusion).

5) Generally, there is a baffling, if not troubling lack of any detail in describing and properly quantifying data and of description of reagents. This pervades the entire manuscript. Below a number of examples (in no order of importance):

a) In panels Figure 2C-F, Figure 3D-G, Figure 4B-E, Figure 4F-I, Figure 5, Figure 6C-F, Figure 3—figure supplement 2, Figure 4—figure supplement 1, Figure 4—figure supplement 2, and Figure 4—figure supplement 3 the authors do not provide any quantification to support their data (incl. the number of animals examined). It is unclear whether in all cases there is a fully penetrant phenotype. In cases of a fully penetrant phenotype (e.g. Figure 2C-F?), a simple explanation on the legend accompanied with the number of animals where the phenotype was observed will suffice. It is odd that the authors do not appear to be familiar with the standards in the field. In the case of partial exon skipping phenotypes, it is unclear to this reviewer whether all animals show always same levels of GFP and RFP expression in ALM (e.g. in the double heterozygotes). For variable phenotypes, a qualitative quantification (%dark/dim/bright red/green), or any other means of quantification deemed appropriate by the authors, is required.

We apologize for the lack of various details, including quantification, in the first draft of the manuscript. Quantification has now been added in figure legends for all figures requested, as well as a few additional figures. This quantification demonstrates that most (but not all) conditions result in fully-penetrant and invariant splicing outcomes.

b) There's no description of the splicing reporter. Is the reporter driven by its own promoter or the endogenous promoter? They refer to a previous paper (Norris, 2014) which states vectors that drove expression of each reporter pan-neuronally, or with endogenous promoter sequences

This has been added in subsection “Microscopy”.

c) The authors generate a number of deletion mutants for different genes using CRISPR/Cas9 (mec-8, mbl-1, not specified for msi-1) however, they don't provide any description on the nature of the molecular lesion other than a brief reference to a previous paper. A schematic of the mutations should be provided in the respective main figures or in a supplementary figure, with a more detailed explanation on the legends or in the Materials and methods section.

This has been added as Figure 3—figure supplement 1

d) There is no list of the strains used on this study as a supplementary file. For example, which is the allele name of the new deletions?

This has been added in the Key Resource table at the beginning of the Materials and methods section.

e) Authors need to describe how new transgenic strains were generated (DNA mix concentrations, simple or complex arrays, etc.).

We have now included details of our microinjection conditions with DNA concentrations, subsection “Strain maintenance”. For all experiments, we generated simple arrays.

f) These newly generated deletion mutants via CRISPR/Cas9 seem to include a GFP reporter construct, subsection “Mutant generation and genetic screening”: "the gene of interest is deleted and is replaced with a heterologous GFP reporter construct". However, authors use fosmid reporters to analyze mec-8 and mbl-1 expression (Figure 4). Do these new mbl-1 and mec-8 deletion alleles recapitulate mbl-1 and mec-8 expression through the GFP reporter? In that case, does this GFP expression interfere with the analysis of the sad-1 splicing reporter green fluorescence?

This has been clarified (subsection “Mutant generation and genetic screening”) to explain that the GFP in the CRISPR alleles does not report on the endogenous expression of the targeted gene, but rather uses a heterologous myo-2 pharyngeal promoter to express GFP, in effect “marking” the deletion allele. Given that the GFP is expressed only in the pharynx, this does not interfere with analysis of the sad-1 reporter in the ALM of VNC neurons.

g) Subsection “RBPs directly mediate sad-1 exon inclusion through interactions with surrounding introns”, "Most neurons besides the ALM and BDU neurons express both skipped and included sad-1 isoforms". Attending to the merge provided in Figure 1C, several neurons in the head and tail seem to exclude exon 15. It will be very informative if authors provide a more detailed expression data of their splicing reporter on Figure 1. For example, just a gross quantification of the% of neurons with included/excluded/both in the head, tail, major ganglia, etc., will greatly help the authors to better illustrate the splicing event under study.

This has been added in Figure 1—figure supplement 1.

h) Figure 4A-E, mec-8 and mbl-1 expression. Was mec-8 fosmid reporter expression detected in other neurons?

We have added text (subsection “TFs affecting sad-1 splicing are required for RBP expression in the ALM neuron”) that mec-8 expression was observed in other neurons in addition to ALM.

Also, authors mention that the mbl-1 fosmid reporter is expressed in "most other neurons" (subsection “TFs affecting sad-1 splicing are required for RBP expression in the ALM neuron”, Figure 4D). To support their point, authors could show whole worm images or high magnification images of specific regions (head, tail, etc.). This needs to be properly documented with an image.

We have added a whole worm image for reference (Figure 4—figure supplement 1).

https://doi.org/10.7554/eLife.46726.023

Article and author information

Author details

  1. Morgan Thompson

    Biological Sciences, Southern Methodist University, Dallas, United States
    Contribution
    Conceptualization, Investigation, Writing—original draft, Writing—review and editing
    Competing interests
    No competing interests declared
  2. Ryan Bixby

    Biological Sciences, Southern Methodist University, Dallas, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  3. Robert Dalton

    Biological Sciences, Southern Methodist University, Dallas, United States
    Contribution
    Investigation, Writing—review and editing
    Competing interests
    No competing interests declared
  4. Alexa Vandenburg

    Biological Sciences, Southern Methodist University, Dallas, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  5. John A Calarco

    Cell & Systems Biology, University of Toronto, Toronto, Canada
    Contribution
    Conceptualization, Writing—review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2197-7801
  6. Adam D Norris

    Biological Sciences, Southern Methodist University, Dallas, United States
    Contribution
    Conceptualization, Investigation, Writing—original draft, Writing—review and editing
    For correspondence
    adnorris@smu.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0570-7414

Funding

Canadian Institutes of Health Research

  • John Calarco

Natural Sciences and Engineering Research Council of Canada

  • John Calarco

Oak Ridge Associated Universities

  • Adam D Norris

NIGMS (R35 GM133461)

  • Adam D Norris

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

Some strains were provided by the CGC, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440). JAC is supported by funding from the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canadian Institutes of Health Research (CIHR), ADN is supported by funding from Oak Ridge Associated Universities (ORAU) and NIH R35 GM133461. Other strains were provided by the National BioResource Project (Tokyo). We thank Megan Norris for critical reeding of the manuscript. MT devised experiments, performed experiments, wrote manuscript. RB performed experiments. RD performed experiments. AV performed experiments. JC devised experiments. AN devised experiments, performed experiments, wrote manuscript.

Senior Editor

  1. James L Manley, Columbia University, United States

Reviewing Editor

  1. Oliver Hobert, Howard Hughes Medical Institute, Columbia University, United States

Reviewer

  1. Hidehito Kuroyanagi, Tokyo Medical and Dental University, Japan

Publication history

  1. Received: March 11, 2019
  2. Accepted: June 20, 2019
  3. Version of Record published: July 19, 2019 (version 1)

Copyright

© 2019, Thompson 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

  • 3,178
    Page views
  • 430
    Downloads
  • 16
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Morgan Thompson
  2. Ryan Bixby
  3. Robert Dalton
  4. Alexa Vandenburg
  5. John A Calarco
  6. Adam D Norris
(2019)
Splicing in a single neuron is coordinately controlled by RNA binding proteins and transcription factors
eLife 8:e46726.
https://doi.org/10.7554/eLife.46726

Further reading

    1. Developmental Biology
    2. Computational and Systems Biology
    Erik S Schild, Shivam Gupta ... Hendrik C Korswagen
    Research Article Updated

    Many developmental processes depend on precise temporal control of gene expression. We have previously established a theoretical framework for regulatory strategies that can govern such high temporal precision, but experimental validation of these predictions was still lacking. Here, we use the time-dependent expression of a Wnt receptor that controls neuroblast migration in Caenorhabditis elegans as a tractable system to study a robust, cell-intrinsic timing mechanism in vivo. Single-molecule mRNA quantification showed that the expression of the receptor increases non-linearly, a dynamic that is predicted to enhance timing precision over an unregulated, linear increase in timekeeper abundance. We show that this upregulation depends on transcriptional activation, providing in vivo evidence for a model in which the timing of receptor expression is regulated through an accumulating activator that triggers expression when a specific threshold is reached. This timing mechanism acts across a cell division that occurs in the neuroblast lineage and is influenced by the asymmetry of the division. Finally, we show that positive feedback of receptor expression through the canonical Wnt pathway enhances temporal precision. We conclude that robust cell-intrinsic timing can be achieved by combining regulation and feedback of the timekeeper gene.

    1. Developmental Biology
    2. Stem Cells and Regenerative Medicine
    Francisco Xavier Galdos, Carissa Lee ... Sean Wu
    Research Article

    During mammalian development, the left and right ventricles arise from early populations of cardiac progenitors known as the first and second heart fields, respectively. While these populations have been extensively studied in non-human model systems, their identification and study in vivo human tissues have been limited due to the ethical and technical limitations of accessing gastrulation stage human embryos. Human induced pluripotent stem cells (hiPSCs) present an exciting alternative for modeling early human embryogenesis due to their well-established ability to differentiate into all embryonic germ layers. Here, we describe the development of a TBX5/MYL2 lineage tracing reporter system that allows for the identification of FHF- progenitors and their descendants including left ventricular cardiomyocytes. Furthermore, using single cell RNA sequencing (scRNA-seq) with oligonucleotide-based sample multiplexing, we extensively profiled differentiating hiPSCs across 12 timepoints in two independent iPSC lines. Surprisingly, our reporter system and scRNA-seq analysis revealed a predominance of FHF differentiation using the small molecule Wnt-based 2D differentiation protocol. We compared this data with existing murine and 3D cardiac organoid scRNA-seq data and confirmed the dominance of left ventricular cardiomyocytes (>90%) in our hiPSC-derived progeny. Together, our work provides the scientific community with a powerful new genetic lineage tracing approach as well as a single cell transcriptomic atlas of hiPSCs undergoing cardiac differentiation.