Co-regulation of alternative splicing by diverse splicing factors in Caenorhabditis elegans

The SR protein RSP-2 influences expression of the truncated insulin receptor DAF-2B in Caenorhabditis elegans

The alternatively spliced daf-2b transcript in Caenorhabditis elegans encodes a truncated isoform of the nematode insulin receptor that retains the extracellular ligand binding domain but lacks the intracellular signaling domain and is therefore unable to transduce a signal. To identify factors that influence expression of daf-2b, we performed a targeted RNA interference screen of rsp genes, which encode splicing factors from the serine/arginine protein family. Loss of rsp-2 significantly increased the expression of a fluorescent daf-2b splicing reporter, as well as increasing expression of endogenous daf-2b transcripts. Correspondingly, rsp-2 mutants exhibited similar phenotypes to those previously observed with DAF-2B overexpression, namely suppression of pheromone-induced dauer formation, enhancement of dauer entry in insulin signaling mutants, inhibition of dauer recovery, and increased lifespan. However, the epistatic relationship between rsp-2 and daf-2b varied according to the experimental context. Increased dauer entry and delayed dauer exit of rsp-2 mutants in an insulin signaling mutant background were partially dependent on daf-2b. Conversely, suppression of pheromone-induced dauer formation and increased lifespan in rsp-2 mutants were independent of daf-2b. These data demonstrate that C. elegans RSP-2, an ortholog of human splicing factor protein SRSF5/SRp40, is involved in regulating the expression of the truncated DAF-2B isoform. However, we also find that RSP-2 can influence dauer formation and lifespan independently of DAF-2B.

Systematic analysis of alternative splicing in time course data using Spycone

MOTIVATION

During disease progression or organism development, alternative splicing may lead to isoform switches that demonstrate similar temporal patterns and reflect the alternative splicing co-regulation of such genes. Tools for dynamic process analysis usually neglect alternative splicing.

RESULTS

Here, we propose Spycone, a splicing-aware framework for time course data analysis. Spycone exploits a novel IS detection algorithm and offers downstream analysis such as network and gene set enrichment. We demonstrate the performance of Spycone using simulated and real-world data of SARS-CoV-2 infection.

AVAILABILITY AND IMPLEMENTATION

The Spycone package is available as a PyPI package. The source code of Spycone is available under the GPLv3 license at https://github.com/yollct/spycone and the documentation at https://spycone.readthedocs.io/en/latest/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

The combinatorial control of alternative splicing in C. elegans

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites—for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes—they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur.

Alternative splicing (AS) is a highly regulated process that is crucial for normal development. It requires the core splicing machinery, but the specific choice of splice site during AS is controlled by splicing factors (SFs) such as ELAV or RBFOX proteins that bind to specific sequences in pre-mRNAs to regulate usage of different splice sites. AS varies across the C. elegans life cycle and here we study how diverse SFs combine to regulate AS during C. elegans development. We selected 4 RRM-containing SFs that are all well studied and that have well-characterised loss-of-function genetic alleles. We find that these SFs regulate many of the same targets, and that combinatorial interactions between these SFs affect both individual splicing events and organism-level phenotypes including specific effects on the neuromuscular system. We further show that SFs combine to regulate splicing of an individual pre-mRNA in two distinct modes—either by binding in close proximity or by binding in a defined order on the pre-mRNA. Finally, we find that the genes whose splicing are most likely to be regulated by multiple SFs are genes that are required for the proper function of the neuromuscular system. These genes are also most likely to have changing AS patterns across development, suggesting that their splicing regulation is highly complex and developmentally regulated. Taken together, our data show that the precise splice variant expressed at any point in development is often the outcome of regulation by multiple SFs.

Quantitative RNA-seq meta-analysis of alternative exon usage in C. elegans

Almost 20 years after the completion of the C. elegans genome sequence, gene structure annotation is still an ongoing process with new evidence for gene variants still being regularly uncovered by additional in-depth transcriptome studies. While alternative splice forms can allow a single gene to encode several functional isoforms, the question of how much spurious splicing is tolerated is still heavily debated. Here we gathered a compendium of 1682 publicly available C. elegans RNA-seq data sets to increase the dynamic range of detection of RNA isoforms, and obtained robust measurements of the relative abundance of each splicing event. While most of the splicing reads come from reproducibly detected splicing events, a large fraction of purported junctions is only supported by a very low number of reads. We devised an automated curation method that takes into account the expression level of each gene to discriminate robust splicing events from potential biological noise. We found that rarely used splice sites disproportionately come from highly expressed genes and are significantly less conserved in other nematode genomes than splice sites with a higher usage frequency. Our increased detection power confirmed trans-splicing for at least 84% of C. elegans protein coding genes. The genes for which trans-splicing was not observed are overwhelmingly low expression genes, suggesting that the mechanism is pervasive but not fully captured by organism-wide RNA-seq. We generated annotated gene models including quantitative exon usage information for the entire C. elegans genome. This allows users to visualize at a glance the relative expression of each isoform for their gene of interest.

CRISPR-mediated genetic interaction profiling identifies RNA binding proteins controlling metazoan fitness

Genetic interaction screens have aided our understanding of complex genetic traits, diseases, and biological pathways. However, approaches for synthetic genetic analysis with null-alleles in metazoans have not been feasible. Here, we present a CRISPR/Cas9-based Synthetic Genetic Interaction (CRISPR-SGI) approach enabling systematic double-mutant generation. Applying this technique in Caenorhabditis elegans, we comprehensively screened interactions within a set of 14 conserved RNA binding protein genes, generating all possible single and double mutants. Many double mutants displayed fitness defects, revealing synthetic interactions. For one interaction between the MBNL1/2 ortholog mbl-1 and the ELAVL ortholog exc-7, double mutants displayed a severely shortened lifespan. Both genes are required for regulating hundreds of transcripts and isoforms, and both may play a critical role in lifespan extension through insulin signaling. Thus, CRISPR-SGI reveals a rich genetic interaction landscape between RNA binding proteins in maintaining organismal health, and will serve as a paradigm applicable to other biological questions.

Splicing factor SUP-12 and the molecular complexity of apparent cooperativity

The splicing factor SUP-12 from C. elegans, in combination with either ASD-1 or FOX-1 from the Fox-1 (RBFOX) family, is required for generating a muscle-specific isoform of the fibroblast growth factor receptor EGL-15. Biophysical techniques have revealed the sequence preference for the RNA Recognition Motif (RRM) domain from SUP-12 as well as the structural details of the RNA-bound complex. Detailed genetics have identified a requisite need for the presence of both SUP-12 and ASD-1/FOX-1 to regulate the alternative splicing event, prompting speculation of a cooperative mechanism between these proteins on binding RNA. In contrast, the interplay between SUP-12 and ASD-1 suggests that although the RRM domains from each protein are in direct contact on the egl-15 pre-mRNA, there is no simple contribution of binding cooperativity. Evidence for an independent binding mechanism by SUP-12 and ASD-1 will be discussed, including a model in which both positive and negative contributions are balanced during complex assembly. The ability to monitor tissue-specific alternative splicing in live nematodes will continue to provide a powerful method to test in vivo mechanistic models derived from atomic-level investigation.

Unusual regulation of splicing of the cholinergic locus in Caenorhabditis elegans

The essential neurotransmitter acetylcholine functions throughout the animal kingdom. In Caenorhabditis elegans, the acetylcholine biosynthetic enzyme [choline acetyltransferase (ChAT)] and vesicular transporter [vesicular acetylcholine transporter (VAChT)] are encoded by the cha-1 and unc-17 genes, respectively. These two genes compose a single complex locus in which the unc-17 gene is nested within the first intron of cha-1, and the two gene products arise from a common pre-messenger RNA (pre-mRNA) by alternative splicing. This genomic organization, known as the cholinergic gene locus (CGL), is conserved throughout the animal kingdom, suggesting that the structure is important for the regulation and function of these genes. However, very little is known about CGL regulation in any species. We now report the identification of an unusual type of splicing regulation in the CGL of C. elegans, mediated by two pairs of complementary sequence elements within the locus. We show that both pairs of elements are required for efficient splicing to the distal acceptor, and we also demonstrate that proper distal splicing depends more on sequence complementarity within each pair of elements than on the sequences themselves. We propose that these sequence elements are able to form stem-loop structures in the pre-mRNA; such structures would favor specific splicing alternatives and thus regulate CGL splicing. We have identified complementary elements at comparable locations in the genomes of representative species of other animal phyla; we suggest that this unusual regulatory mechanism may be a general feature of CGLs.

Backbone-independent nucleic acid binding by splicing factor SUP-12 reveals key aspects of molecular recognition

Cellular differentiation is frequently accompanied by alternative splicing, enabled by the expression of tissue-specific factors which bind to pre-mRNAs and regulate exon choice. During Caenorhabditis elegans development, muscle-specific expression of the splicing factor SUP-12, together with a member of the Fox-1 family of splicing proteins, generates a functionally distinct isoform of the fibroblast growth factor receptor EGL-15. Using a combination of NMR spectroscopy and isothermal titration calorimetry, we determined the mode of nucleic acid binding by the RNA recognition motif domain of SUP-12. The calculated structures provide the first atomic details of RNA and DNA binding by the family of proteins that include SUP-12, RBM24, RBM38/RNPC1, SEB-4 and XSeb4R. This information was further used to design strategic mutations to probe the interaction with ASD-1 and to quantitatively perturb splicing in vivo.

A pair of RNA-binding proteins controls networks of splicing events contributing to specialization of neural cell types

Alternative splicing is important for the development and function of the nervous system, but little is known about the differences in alternative splicing between distinct types of neurons. Furthermore, the factors that control cell-type-specific splicing and the physiological roles of these alternative isoforms are unclear. By monitoring alternative splicing at single-cell resolution in Caenorhabditis elegans, we demonstrate that splicing patterns in different neurons are often distinct and highly regulated. We identify two conserved RNA-binding proteins, UNC-75/CELF and EXC-7/Hu/ELAV, which regulate overlapping networks of splicing events in GABAergic and cholinergic neurons. We use the UNC-75 exon network to discover regulators of synaptic transmission and to identify unique roles for isoforms of UNC-64/Syntaxin, a protein required for synaptic vesicle fusion. Our results indicate that combinatorial regulation of alternative splicing in distinct neurons provides a mechanism to specialize metazoan nervous systems.

The multiplicity of alternative splicing decisions in Caenorhabditis elegans is linked to specific intronic regulatory motifs and minisatellites

Background

Alternative splicing diversifies the pool of messenger RNA molecules encoded by individual genes. This diversity is particularly high when multiple splicing decisions cause a combinatorial arrangement of several alternate exons. We know very little on how the multiple decisions occurring during the maturation of single transcripts are coordinated and whether specific sequence elements might be involved.

Results

Here, the Caenorhabditis elegans genome was surveyed in order to identify sequence elements that might play a specific role in the regulation of multiple splicing decisions. The introns flanking alternate exons in transcripts whose maturation involves multiple alternative splicing decisions were compared to those whose maturation involves a single decision. Fifty-eight penta-, hexa-, and hepta-meric elements, clustered in 17 groups, were significantly over-represented in genes subject to multiple alternative splicing decisions. Most of these motifs relate to known splicing regulatory elements and appear to be well conserved in the related species Caenorhabditis briggsae. The usage of specific motifs is not linked to the gene product function, but rather depends on the gene structure, since it is influenced by the distance separating the multiple splicing decision sites. Two of these motifs are part of the CeRep25B minisatellite, which is also over-represented at the vicinity of alternative splicing regions. Most of the remaining motifs are not part of repeated sequence elements, but tend to occur in specific heterologous pairs in genes subject to multiple alternative splicing decisions.

Conclusions

The existence of specific intronic sequence elements linked to multiple alternative splicing decisions is intriguing and suggests that these elements might have some specialized regulatory role during splicing.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-364) contains supplementary material, which is available to authorized users.

Alternative Splicing Regulation of Cancer-Related Pathways in Caenorhabditis elegans: An In Vivo Model System with a Powerful Reverse Genetics Toolbox

Alternative splicing allows for the generation of protein diversity and fine-tunes gene expression. Several model systems have been used for the in vivo study of alternative splicing. Here we review the use of the nematode Caenorhabditis elegans to study splicing regulation in vivo. Recent studies have shown that close to 25% of genes in the worm genome undergo alternative splicing. A big proportion of these events are functional, conserved, and under strict regulation either across development or other conditions. Several techniques like genome-wide RNAi screens and bichromatic reporters are available for the study of alternative splicing in worms. In this review, we focus, first, on the main studies that have been performed to dissect alternative splicing in this system and later on examples from genes that have human homologs that are implicated in cancer. The significant advancement towards understanding the regulation of alternative splicing and cancer that the C. elegans system has offered is discussed.

CELF family RNA-binding protein UNC-75 regulates two sets of mutually exclusive exons of the unc-32 gene in neuron-specific manners in Caenorhabditis elegans

An enormous number of alternative pre–mRNA splicing patterns in multicellular organisms are coordinately defined by a limited number of regulatory proteins and cis elements. Mutually exclusive alternative splicing should be strictly regulated and is a challenging model for elucidating regulation mechanisms. Here we provide models of the regulation of two sets of mutually exclusive exons, 4a–4c and 7a–7b, of the Caenorhabditis elegans uncoordinated (unc)-32 gene, encoding the a subunit of V₀ complex of vacuolar-type H⁺-ATPases. We visualize selection patterns of exon 4 and exon 7 in vivo by utilizing a trio and a pair of symmetric fluorescence splicing reporter minigenes, respectively, to demonstrate that they are regulated in tissue-specific manners. Genetic analyses reveal that RBFOX family RNA–binding proteins ASD-1 and FOX-1 and a UGCAUG stretch in intron 7b are involved in the neuron-specific selection of exon 7a. Through further forward genetic screening, we identify UNC-75, a neuron-specific CELF family RNA–binding protein of unknown function, as an essential regulator for the exon 7a selection. Electrophoretic mobility shift assays specify a short fragment in intron 7a as the recognition site for UNC-75 and demonstrate that UNC-75 specifically binds via its three RNA recognition motifs to the element including a UUGUUGUGUUGU stretch. The UUGUUGUGUUGU stretch in the reporter minigenes is actually required for the selection of exon 7a in the nervous system. We compare the amounts of partially spliced RNAs in the wild-type and unc-75 mutant backgrounds and raise a model for the mutually exclusive selection of unc-32 exon 7 by the RBFOX family and UNC-75. The neuron-specific selection of unc-32 exon 4b is also regulated by UNC-75 and the unc-75 mutation suppresses the Unc phenotype of the exon-4b-specific allele of unc-32 mutants. Taken together, UNC-75 is the neuron-specific splicing factor and regulates both sets of the mutually exclusive exons of the unc-32 gene.

Tissue-specific and mutually exclusive alternative pre–mRNA splicing is a challenging model for elucidating regulation mechanisms. We previously demonstrated that evolutionarily conserved RBFOX family RNA–binding proteins ASD-1 and FOX-1 and a muscle-specific RNA–binding protein SUP-12 cooperatively direct muscle-specific selection of exon 5B of the C. elegans egl-15 gene. Here we demonstrate that two sets of mutually exclusive exons, 4a–4c and 7a–7b, of the unc-32 gene are regulated in tissue-specific manners and that ASD-1 and FOX-1, expressed in a variety of tissues, can regulate the neuron-specific selection of unc-32 exon 7a in combination with the neuron-specific CELF family RNA–binding protein UNC-75. We determine the cis-elements for the RBFOX family and UNC-75, which separately reside in intron 7b and intron 7a, respectively. By analyzing the partially spliced RNA species, we propose the orders of intron removal and the sites of action for the RBFOX family and UNC-75 in the mutually exclusive selection of exon 7a and exon 7b. We also demonstrate that UNC-75 regulates the neuron-specific selection of exon 4b and propose the models of the mutually exclusive selection of exons 4a, 4b, and 4c. These studies thus provide novel modes of regulation for tissue-specific and mutually exclusive alternative splicing in vivo.

Position-dependent and neuron-specific splicing regulation by the CELF family RNA-binding protein UNC-75 in Caenorhabditis elegans

A large fraction of protein-coding genes in metazoans undergo alternative pre-mRNA splicing in tissue- or cell-type-specific manners. Recent genome-wide approaches have identified many putative-binding sites for some of tissue-specific trans-acting splicing regulators. However, the mechanisms of splicing regulation in vivo remain largely unknown. To elucidate the modes of splicing regulation by the neuron-specific CELF family RNA-binding protein UNC-75 in Caenorhabditis elegans, we performed deep sequencing of poly(A)⁺ RNAs from the unc-75(+)- and unc-75-mutant worms and identified more than 20 cassette and mutually exclusive exons repressed or activated by UNC-75. Motif searches revealed that (G/U)UGUUGUG stretches are enriched in the upstream and downstream introns of the UNC-75-repressed and -activated exons, respectively. Recombinant UNC-75 protein specifically binds to RNA fragments carrying the (G/U)UGUUGUG stretches in vitro. Bi-chromatic fluorescence alternative splicing reporters revealed that the UNC-75-target exons are regulated in tissue-specific and (G/U)UGUUGUG element-dependent manners in vivo. The unc-75 mutation affected the splicing reporter expression specifically in the nervous system. These results indicate that UNC-75 regulates alternative splicing of its target exons in neuron-specific and position-dependent manners through the (G/U)UGUUGUG elements in C. elegans. This study thus reveals the repertoire of target events for the CELF family in the living organism.

During embryogenesis, esrp1 expression is restricted to a subset of epithelial cells and is associated with splicing of a number of developmentally important genes

BACKGROUND

Development of a mature organism from a single cell requires a series of important morphological changes, which is in part regulated by alternative splicing. In this article, we report the expression of Esrp1 during early mouse embryogenesis, a splicing factor implicated in epithelial to mesenchymal transitions.

RESULTS

By qRT-PCR, we find higher expression of Esrp1 and Esrp2 in placenta compared to the embryos. We also find a correlation between the expression of Esrp1 and alternative splicing of several known target exons. Using in situ RNA hybridization we show that while Esrp1 expression is ubiquitous in embryonic day (E)6.5 mouse embryos, expression becomes restricted to the chorion and definitive endoderm starting at E7.5. Esrp1 expression was consistently restricted to a subset of epithelial cell types in developing embryos from E9.5 to E13.5.

CONCLUSIONS

Our results suggest that Esrp1 could play an important role in the morphological changes underlying embryogenesis of the placenta and embryo.

Muscle-specific splicing factors ASD-2 and SUP-12 cooperatively switch alternative pre-mRNA processing patterns of the ADF/cofilin gene in Caenorhabditis elegans

Pre-mRNAs are often processed in complex patterns in tissue-specific manners to produce a variety of protein isoforms from single genes. However, mechanisms orchestrating the processing of the entire transcript are not well understood. Muscle-specific alternative pre-mRNA processing of the unc-60 gene in Caenorhabditis elegans, encoding two tissue-specific isoforms of ADF/cofilin with distinct biochemical properties in regulating actin organization, provides an excellent in vivo model of complex and tissue-specific pre-mRNA processing; it consists of a single first exon and two separate series of downstream exons. Here we visualize the complex muscle-specific processing pattern of the unc-60 pre-mRNA with asymmetric fluorescence reporter minigenes. By disrupting juxtaposed CUAAC repeats and UGUGUG stretch in intron 1A, we demonstrate that these elements are required for retaining intron 1A, as well as for switching the processing patterns of the entire pre-mRNA from non-muscle-type to muscle-type. Mutations in genes encoding muscle-specific RNA-binding proteins ASD-2 and SUP-12 turned the colour of the unc-60 reporter worms. ASD-2 and SUP-12 proteins specifically and cooperatively bind to CUAAC repeats and UGUGUG stretch in intron 1A, respectively, to form a ternary complex in vitro. Immunohistochemical staining and RT-PCR analyses demonstrate that ASD-2 and SUP-12 are also required for switching the processing patterns of the endogenous unc-60 pre-mRNA from UNC-60A to UNC-60B in muscles. Furthermore, systematic analyses of partially spliced RNAs reveal the actual orders of intron removal for distinct mRNA isoforms. Taken together, our results demonstrate that muscle-specific splicing factors ASD-2 and SUP-12 cooperatively promote muscle-specific processing of the unc-60 gene, and provide insight into the mechanisms of complex pre-mRNA processing; combinatorial regulation of a single splice site by two tissue-specific splicing regulators determines the binary fate of the entire transcript.

A role for RNA post-transcriptional regulation in satellite cell activation

Background

Satellite cells are resident skeletal muscle stem cells responsible for muscle maintenance and repair. In resting muscle, satellite cells are maintained in a quiescent state. Satellite cell activation induces the myogenic commitment factor, MyoD, and cell cycle entry to facilitate transition to a population of proliferating myoblasts that eventually exit the cycle and regenerate muscle tissue. The molecular mechanism involved in the transition of a quiescent satellite cell to a transit-amplifying myoblast is poorly understood.

Methods

Satellite cells isolated by FACS from uninjured skeletal muscle and 12 h post-muscle injury from wild type and Syndecan-4 null mice were probed using Affymetrix 430v2 gene chips and analyzed by Spotfire^tm and Ingenuity Pathway analysis to identify gene expression changes and networks associated with satellite cell activation, respectively. Additional analyses of target genes identify miRNAs exhibiting dynamic changes in expression during satellite cell activation. The function of the miRNAs was assessed using miRIDIAN hairpin inhibitors.

Results

An unbiased gene expression screen identified over 4,000 genes differentially expressed in satellite cells in vivo within 12 h following muscle damage and more than 50% of these decrease dramatically. RNA binding proteins and genes involved in post-transcriptional regulation were significantly over-represented whereas splicing factors were preferentially downregulated and mRNA stability genes preferentially upregulated. Furthermore, six computationally identified miRNAs demonstrated novel expression through muscle regeneration and in satellite cells. Three of the six miRNAs were found to regulate satellite cell fate.

Conclusions

The quiescent satellite cell is actively maintained in a state poised to activate in response to external signals. Satellite cell activation appears to be regulated by post-transcriptional gene regulation.

An integrated regulatory network reveals pervasive cross-regulation among transcription and splicing factors

Traditionally the gene expression pathway has been regarded as being comprised of independent steps, from RNA transcription to protein translation. To date there is increasing evidence of coupling between the different processes of the pathway, specifically between transcription and splicing. To study the interplay between these processes we derived a transcription-splicing integrated network. The nodes of the network included experimentally verified human proteins belonging to three groups of regulators: transcription factors, splicing factors and kinases. The nodes were wired by instances of predicted transcriptional and alternative splicing regulation. Analysis of the network indicated a pervasive cross-regulation among the nodes; specifically, splicing factors are significantly more connected by alternative splicing regulatory edges relative to the two other subgroups, while transcription factors are more extensively controlled by transcriptional regulation. Furthermore, we found that splicing factors are the most regulated of the three regulatory groups and are subject to extensive combinatorial control by alternative splicing and transcriptional regulation. Consistent with the network results, our bioinformatics analyses showed that the subgroup of kinases have the highest density of predicted phosphorylation sites. Overall, our systematic study reveals that an organizing principle in the logic of integrated networks favor the regulation of regulatory proteins by the specific regulation they conduct. Based on these results, we propose a new regulatory paradigm postulating that gene expression regulation of the master regulators in the cell is predominantly achieved by cross-regulation.

Intragenic alternative splicing coordination is essential for Caenorhabditis elegans slo-1 gene function

Alternative splicing is critical for diversifying eukaryotic proteomes, but the rules governing and coordinating splicing events among multiple alternate splice sites within individual genes are not well understood. We developed a quantitative PCR-based strategy to quantify the expression of the 12 transcripts encoded by the Caenorhabditis elegans slo-1 gene, containing three alternate splice sites. Using conditional probability-based models, we show that splicing events are coordinated across these sites. Further, we identify a point mutation in an intron adjacent to one alternate splice site that disrupts alternative splicing at all three sites. This mutation leads to aberrant synaptic transmission at the neuromuscular junction. In a genomic survey, we found that a UAAAUC element disrupted by this mutation is enriched in introns flanking alternate exons in genes with multiple alternate splice sites. These results establish that proper coordination of intragenic alternative splicing is essential for normal physiology of slo-1 in vivo and identify putative specialized cis-regulatory elements that regulate the coordination of intragenic alternative splicing.

Seeing elegance in gene regulatory networks of the worm

There has been a recent explosion in the wealth of genomic data available to C. elegans researchers, as efforts to characterize gene expression and its regulators at a molecular level have borne significant fruit. Detailed measurement of gene expression at a variety of developmental stages, and in numerous individual tissues, has dramatically increased our understanding of cell-type-specific gene expression networks. Characterization of the targets of transcription factors, chromatin-binding proteins, and miRNAs has provided genome-wide insights into the mechanisms governing gene expression. Development of new techniques have allowed this characterization to begin to shift from whole-organism studies to tissue-level, and even single-cell-level profiling, creating a first glimpse into gene regulatory circuits at the single-cell level in a living organism. Integration of these datasets has yielded novel insights into evolution, gene expression regulation, and the link between sequence and phenotype.