Geometry of gene regulatory dynamics.

Embryonic development leads to the reproducible and ordered appearance of complexity from egg to adult. The successive differentiation of different cell types that elaborate this complexity results from the activity of gene networks and was likened by Waddington to a flow through a landscape in which valleys represent alternative fates. Geometric methods allow the formal representation of such landscapes and codify the types of behaviors that result from systems of differential equations. Results from Smale and coworkers imply that systems encompassing gene network models can be represented as potential gradients with a Riemann metric, justifying the Waddington metaphor. Here, we extend this representation to include parameter dependence and enumerate all three-way cellular decisions realizable by tuning at most two parameters, which can be generalized to include spatial coordinates in a tissue. All diagrams of cell states vs. model parameters are thereby enumerated. We unify a number of standard models for spatial pattern formation by expressing them in potential form (i.e., as topographic elevation). Turing systems appear nonpotential, yet in suitable variables the dynamics are low dimensional and potential. A time-independent embedding recovers the original variables. Lateral inhibition is described by a saddle point with many unstable directions. A model for the patterning of the Drosophila eye appears as relaxation in a bistable potential. Geometric reasoning provides intuitive dynamic models for development that are well adapted to fit time-lapse data.

Genome-wide analyses of light-regulated genes in Aspergillus nidulans reveal a complex interplay between different photoreceptors and novel photoreceptor functions.

Fungi sense light of different wavelengths using blue-, green-, and red-light photoreceptors. Blue light sensing requires the "white-collar" proteins with flavin as chromophore, and red light is sensed through phytochrome. Here we analyzed genome-wide gene expression changes caused by short-term, low-light intensity illumination with blue-, red- or far-red light in Aspergillus nidulans and found that more than 1100 genes were differentially regulated. The largest number of up- and downregulated genes depended on the phytochrome FphA and the attached HOG pathway. FphA and the white-collar orthologue LreA fulfill activating but also repressing functions under all light conditions and both appear to have roles in the dark. Additionally, we found about 100 genes, which are red-light induced in the absence of phytochrome, suggesting alternative red-light sensing systems. We also found blue-light induced genes in the absence of the blue-light receptor LreA. We present evidence that cryptochrome may be part of this regulatory cue, but that phytochrome is essential for the response. In addition to in vivo data showing that FphA is involved in blue-light sensing, we performed spectroscopy of purified phytochrome and show that it responds indeed to blue light.

LuxT controls specific quorum-sensing-regulated behaviors in Vibrionaceae spp. via repression of qrr1, encoding a small regulatory RNA.

Quorum sensing (QS) is a process of chemical communication bacteria use to transition between individual and collective behaviors. QS depends on the production, release, and synchronous response to signaling molecules called autoinducers (AIs). The marine bacterium Vibrio harveyi monitors AIs using a signal transduction pathway that relies on five small regulatory RNAs (called Qrr1-5) that post-transcriptionally control target genes. Curiously, the small RNAs largely function redundantly making it difficult to understand the necessity for five of them. Here, we identify LuxT as a transcriptional repressor of qrr1. LuxT does not regulate qrr2-5, demonstrating that qrr genes can be independently controlled to drive unique downstream QS gene expression patterns. LuxT reinforces its control over the same genes it regulates indirectly via repression of qrr1, through a second transcriptional control mechanism. Genes dually regulated by LuxT specify public goods including an aerolysin-type pore-forming toxin. Phylogenetic analyses reveal that LuxT is conserved among Vibrionaceae and sequence comparisons predict that LuxT represses qrr1 in additional species. The present findings reveal that the QS regulatory RNAs can carry out both shared and unique functions to endow bacteria with plasticity in their output behaviors.

Multilevel prioritization of gene regulators associated with consensus molecular subtypes of colorectal cancer.

Consensus molecular subtypes (CMSs) are emerging as critical factor for prognosis and treatment of colorectal cancer. Gene regulators, including chromatin regulator, RNA-binding protein and transcriptional factor, are critical modulators of cancer hallmark, yet little is known regarding the underlying functional mechanism in CMSs. Herein, we identified a core set of 235 functional gene regulators (FGRs) by integrating genome, epigenome, transcriptome and interactome of CMSs. FGRs exhibited significant multi-omics alterations and impacts on cell lines growth, as well as significantly enriched cancer driver genes and pathways. Moreover, common FGRs played different roles in the context of CMSs. In accordance with the immune characteristics of CMSs, we found that the anti-tumor immune pathways were mainly activated by FGRs (e.g. STAT1 and CREBBP) in CMS1, while inhibited by FGRs in CMS2-4. FGRs mediated aberrant expression of ligands, which bind to receptor on immune cells, and modulated tumor immune microenvironment of subtypes. Intriguingly, systematic exploration of datasets using genomic and transcriptome co-similarity reveals the coordinated manner in FGRs act in CMSs to orchestrate their pathways and patients' prognosis. Expression signatures of the FGRs revealed an optimized CMS classifier, which demonstrated 88% concordance with the gold-standard classifier, but avoiding the influence of sample composition. Overall, our integrative analysis identified FGRs to regulate core tumorigenic processes/pathways across CMSs.

Controlling gene expression timing through gene regulatory architecture.

Gene networks typically involve the regulatory control of multiple genes with related function. This connectivity enables correlated control of the levels and timing of gene expression. Here we study how gene expression timing in the single-input module motif can be encoded in the regulatory DNA of a gene. Using stochastic simulations, we examine the role of binding affinity, TF regulatory function and network size in controlling the mean first-passage time to reach a fixed fraction of steady-state expression for both an auto-regulated TF gene and a target gene. We also examine how the variability in first-passage time depends on these factors. We find that both network size and binding affinity can dramatically speed up or slow down the response time of network genes, in some cases predicting more than a 100-fold change compared to that for a constitutive gene. Furthermore, these factors can also significantly impact the fidelity of this response. Importantly, these effects do not occur at "extremes" of network size or binding affinity, but rather in an intermediate window of either quantity.

A non-coding genetic variant associated with abdominal aortic aneurysm alters ERG gene regulation.

Abdominal aortic aneurysm (AAA) is a major cause of sudden death in the elderly. While AAA has some overlapping genetic and environmental risk factors with atherosclerosis, there are substantial differences, and AAA-specific medication is lacking. A recent meta-analysis of genome-wide association studies has identified four novel single-nucleotide polymorphisms (SNPs) specifically associated with AAA. Here, we investigated the gene regulatory function for one of four non-coding SNPs associated with AAA, rs2836411, which is located in an intron of the ERG gene. Rs2836411 resides within a >70 kb super-enhancer that has high levels of H3K27ac and H3K4me1 in vascular endothelial and haematopoietic cell types. Enhancer luciferase assays in cell lines showed that the risk allele significantly alters enhancer activity. The risk allele also correlates with reduced ERG expression in aortic and other vascular tissues. To identify whether rs2836411 directly contacts the promoters of ERG and/or of genes further away, we performed allele-specific circular chromosome conformation capture sequencing. In vascular endothelial cells, which express ERG, the SNP region interacts highly within the super-enhancer, while in vascular smooth muscle cells, which do not express ERG, the interactions are distributed across a wider region that includes neighbouring genes. Furthermore, the risk allele has fewer interactions within the super-enhancer compared to the protective allele. In conclusion, our results indicate that rs2836411 likely affects ERG expression by altering enhancer activity and changing local chromatin interactions. ERG is involved in vascular development, angiogenesis, and inflammation in atherosclerosis; therefore mechanistically, rs2836411 could contribute to AAA by modulating ERG levels.

Mrr1 regulation of methylglyoxal catabolism and methylglyoxal-induced fluconazole resistance in Candida lusitaniae.

Transcription factor Mrr1, best known for its regulation of Candida azole resistance genes such as MDR1, regulates other genes that are poorly characterized. Among the other Mrr1-regulated genes are putative methylglyoxal reductases. Methylglyoxal (MG) is a toxic metabolite that is elevated in diabetes, uremia, and sepsis, which are diseases that increase the risk for candidiasis, and MG serves as a regulatory signal in diverse organisms. Our studies in Clavispora lusitaniae, also known as Candida lusitaniae, showed that Mrr1 regulates expression of two paralogous MG reductases, MGD1 and MGD2, and that both participate in MG resistance and MG catabolism. Exogenous MG increased Mrr1-dependent expression of MGD1 and MGD2 as well as expression of MDR1, which encodes an efflux pump that exports fluconazole. MG improved growth in the presence of fluconazole and this was largely Mrr1-dependent with contributions from a secondary transcription factor, Cap1. Increased fluconazole resistance was also observed in mutants lacking Glo1, a Mrr1-independent MG catabolic enzyme. Isolates from other Candida species displayed heterogeneity in MG resistance and MG stimulation of azole resistance. We propose endogenous and host-derived MG can induce MDR1 and other Mrr1-regulated genes causing increased drug resistance, which may contribute to some instances of fungal treatment failure.

The regulatory gene wblA is a target of the orphan response regulator OrrA in Streptomyces coelicolor.

Our previous study using transposon mutagenesis indicated that disruption of the putative response regulator gene orrA impacted antibiotic production in Streptomyces coelicolor. In this study, the role of OrrA was further characterized by comparing the phenotypes and transcriptomic profiles of the wild-type S. coelicolor strain M145 and ΔorrA, a strain with an inactivated orrA gene. Chromatin immunoprecipitation using a strain expressing OrrA fused with FLAG showed that OrrA binds the promoter of wblA, whose expression was downregulated in ΔorrA. The interaction of OrrA with the wblA promoter was further validated by a pull-down assay. Similar to ΔorrA, the deletion mutant of wblA (ΔwblA) was defective in development, and developmental genes were expressed at similar levels in ΔorrA and ΔwblA. Although both OrrA and WblA downregulated actinorhodin and undecylprodigiosin, their roles in regulation of the calcium-dependent antibiotic and yellow-pigmented type I polyketide differed. sco1375, a gene of unknown function, was identified as another OrrA target, and overexpression of either sco1375 or wblA in ΔorrA partially restored the wild-type phenotype, indicating that these genes mediate some of the effects of OrrA. This study revealed targets of OrrA and provided more insights into the role of the orphan response regulator OrrA in Streptomyces.

Maternal and zygotic gene regulatory effects of endogenous RNAi pathways.

Endogenous small RNAs (sRNAs) and Argonaute proteins are ubiquitous regulators of gene expression in germline and somatic tissues. sRNA-Argonaute complexes are often expressed in gametes and are consequently inherited by the next generation upon fertilization. In Caenorhabditis elegans, 26G-RNAs are primary endogenous sRNAs that trigger the expression of downstream secondary sRNAs. Two subpopulations of 26G-RNAs exist, each of which displaying strongly compartmentalized expression: one is expressed in the spermatogenic gonad and associates with the Argonautes ALG-3/4; plus another expressed in oocytes and in embryos, which associates with the Argonaute ERGO-1. The determinants and dynamics of gene silencing elicited by 26G-RNAs are largely unknown. Here, we provide diverse new insights into these endogenous sRNA pathways of C. elegans. Using genetics and deep sequencing, we dissect a maternal effect of the ERGO-1 branch of the 26G-RNA pathway. We find that maternal primary sRNAs can trigger the production of zygotic secondary sRNAs that are able to silence targets, even in the absence of zygotic primary triggers. Thus, the interaction of maternal and zygotic sRNA populations, assures target gene silencing throughout animal development. Furthermore, we explore other facets of 26G-RNA biology related to the ALG-3/4 branch. We find that sRNA abundance, sRNA pattern of origin and the 3' UTR length of target transcripts are predictors of the regulatory outcome by the Argonautes ALG-3/4. Lastly, we provide evidence suggesting that ALG-3 and ALG-4 regulate their own mRNAs in a negative feedback loop. Altogether, we provide several new regulatory insights on the dynamics, target regulation and self-regulation of the endogenous RNAi pathways of C. elegans.

Identification of the Regulatory Genes of UV-B-Induced Anthocyanin Biosynthesis in Pepper Fruit.

Fruit peels of certain pepper (Capsicum annum L.) varieties accumulate a large amount of anthocyanins and exhibit purple color under medium-wave ultraviolet (UV-B) conditions, which severely impacts the commodity value of peppers. However, the regulatory mechanism of the above process has not been well studied so far. To explore which key genes are involved in this regulatory mechanism, pepper variety 19Q6100, the fruit peels of which turn purple under UV-B conditions, was investigated in this study. Transcription factors with expression levels significantly impacted by UV-B were identified by RNA-seq. Those genes may be involved in the regulation of UV-B-induced anthocyanin biosynthesis. Yeast one-hybrid results revealed that seven transcription factors, CabHLH143, CaMYB113, CabHLH137, CaMYBG, CaWRKY41, CaWRKY44 and CaWRKY53 directly bound to the putative promotor regions of the structural genes in the anthocyanin biosynthesis pathway. CaMYB113 was found to interact with CabHLH143 and CaHY5 by yeast two-hybrid assay, and those three genes may participate collaboratively in UV-B-induced anthocyanin biosynthesis in pepper fruit. Virus-induced gene silencing (VIGS) indicated that fruit peels of CaMYB113-silenced plants were unable to turn purple under UV-B conditions. These findings could deepen our understanding of UV-B-induced anthocyanin biosynthesis in pepper.

Chromatin regulatory genes differentially interact in networks to facilitate distinct GAL1 activity and noise profiles.

Controlling chromatin state constitutes a major regulatory step in gene expression regulation across eukaryotes. While global cellular features or processes are naturally impacted by chromatin state alterations, little is known about how chromatin regulatory genes interact in networks to dictate downstream phenotypes. Using the activity of the canonical galactose network in yeast as a model, here, we measured the impact of the disruption of key chromatin regulatory genes on downstream gene expression, genetic noise and fitness. Using Trichostatin A and nicotinamide, we characterized how drug-based modulation of global histone deacetylase activity affected these phenotypes. Performing epistasis analysis, we discovered phenotype-specific genetic interaction networks of chromatin regulators. Our work provides comprehensive insights into how the galactose network activity is affected by protein interaction networks formed by chromatin regulators.

dCas9-targeted locus-specific protein isolation method identifies histone gene regulators.

Eukaryotic gene regulation is a complex process, often coordinated by the action of tens to hundreds of proteins. Although previous biochemical studies have identified many components of the basal machinery and various ancillary factors involved in gene regulation, numerous gene-specific regulators remain undiscovered. To comprehensively survey the proteome directing gene expression at a specific genomic locus of interest, we developed an in vitro nuclease-deficient Cas9 (dCas9)-targeted chromatin-based purification strategy, called "CLASP" (Cas9 locus-associated proteome), to identify and functionally test associated gene-regulatory factors. Our CLASP method, coupled to mass spectrometry and functional screens, can be efficiently adapted for isolating associated regulatory factors in an unbiased manner targeting multiple genomic loci across different cell types. Here, we applied our method to isolate the Drosophila melanogaster histone cluster in S2 cells to identify several factors including Vig and Vig2, two proteins that bind and regulate core histone H2A and H3 mRNA via interaction with their 3' UTRs.

Structural basis for the regulation of ß-glucuronidase expression by human gut Enterobacteriaceae.

The gut microbiota harbor diverse ß-glucuronidase (GUS) enzymes that liberate glucuronic acid (GlcA) sugars from small-molecule conjugates and complex carbohydrates. However, only the Enterobacteriaceae family of human gut-associated Proteobacteria maintain a GUS operon under the transcriptional control of a glucuronide repressor, GusR. Despite its potential importance in Escherichia, Salmonella, Klebsiella, Shigella, and Yersinia opportunistic pathogens, the structure of GusR has not been examined. Here, we explore the molecular basis for GusR-mediated regulation of GUS expression in response to small-molecule glucuronides. Presented are 2.1-Å-resolution crystal structures of GusRs from Escherichia coli and Salmonella enterica in complexes with a glucuronide ligand. The GusR-specific DNA operator site in the regulatory region of the E. coli GUS operon is identified, and structure-guided GusR mutants pinpoint the residues essential for DNA binding and glucuronide recognition. Interestingly, the endobiotic estradiol-17-glucuronide and the xenobiotic indomethacin-acyl-glucuronide are found to exhibit markedly differential binding to these GusR orthologs. Using structure-guided mutations, we are able to transfer E. coli GusR's preferential DNA and glucuronide binding affinity to S. enterica GusR. Structures of putative GusR orthologs from GUS-encoding Firmicutes species also reveal functionally unique features of the Enterobacteriaceae GusRs. Finally, dominant-negative GusR variants are validated in cell-based studies. These data provide a molecular framework toward understanding the control of glucuronide utilization by opportunistic pathogens in the human gut.

Response of Pseudomonas aeruginosa to the Innate Immune System-Derived Oxidants Hypochlorous Acid and Hypothiocyanous Acid.

Pseudomonas aeruginosa is a significant nosocomial pathogen and is associated with lung infections in cystic fibrosis (CF). Once established, P. aeruginosa infections persist and are rarely eradicated despite host immune cells producing antimicrobial oxidants, including hypochlorous acid (HOCl) and hypothiocyanous acid (HOSCN). There is limited knowledge as to how P. aeruginosa senses, responds to, and protects itself against HOCl and HOSCN and the contribution of such responses to its success as a CF pathogen. To investigate the P. aeruginosa response to these oxidants, we screened 707 transposon mutants, with mutations in regulatory genes, for altered growth following HOCl exposure. We identified regulators of antibiotic resistance, methionine biosynthesis, catabolite repression, and PA14_07340, the homologue of the Escherichia coli HOCl-sensor RclR (30% identical), which are required for protection against HOCl. We have shown that RclR (PA14_07340) protects specifically against HOCl and HOSCN stress and responds to both oxidants by upregulating the expression of a putative peroxiredoxin, rclX (PA14_07355). Transcriptional analysis revealed that while there was specificity in the response to HOCl (231 genes upregulated) and HOSCN (105 genes upregulated), there was considerable overlap, with 74 genes upregulated by both oxidants. These included genes encoding the type 3 secretion system, sulfur and taurine transport, and the MexEF-OprN efflux pump. RclR coordinates part of the response to both oxidants, including upregulation of pyocyanin biosynthesis genes, and, in the presence of HOSCN, downregulation of chaperone genes. These data indicate that the P. aeruginosa response to HOCl and HOSCN is multifaceted, with RclR playing an essential role.IMPORTANCE The bacterial pathogen Pseudomonas aeruginosa causes devastating infections in immunocompromised hosts, including chronic lung infections in cystic fibrosis patients. To combat infection, the host's immune system produces the antimicrobial oxidants hypochlorous acid (HOCl) and hypothiocyanous acid (HOSCN). Little is known about how P. aeruginosa responds to and survives attack from these oxidants. To address this, we carried out two approaches: a mutant screen and transcriptional study. We identified the P. aeruginosa transcriptional regulator, RclR, which responds specifically to HOCl and HOSCN stress and is essential for protection against both oxidants. We uncovered a link between the P. aeruginosa transcriptional response to these oxidants and physiological processes associated with pathogenicity, including antibiotic resistance and the type 3 secretion system.

The eXceptional nature of the X chromosome.

The X chromosome is unique in the genome. In this review we discuss recent advances in our understanding of the genetics and epigenetics of the X chromosome. The X chromosome shares limited conservation with its ancestral homologue the Y chromosome and the resulting difference in X-chromosome dosage between males and females is largely compensated for by X-chromosome inactivation. The process of inactivation is initiated by the long non-coding RNA X-inactive specific transcript (XIST) and achieved through interaction with multiple synergistic silencing pathways. Identification of Xist-interacting proteins has given insight into these processes yet the cascade of events from initiation to maintenance have still to be resolved. In particular, the initiation of inactivation in humans has been challenging to study as: it occurs very early in development; most human embryonic stem cell lines already have an inactive X; and the process seems to differ from mouse. Another difference between human and mouse X inactivation is the larger number of human genes that escape silencing. In humans over 20% of X-linked genes continue to be expressed from the otherwise inactive X chromosome. We are only beginning to understand how such escape occurs but there is growing recognition that escapees contribute to sexually dimorphic traits. The unique biology and epigenetics of the X chromosome have often led to its exclusion from disease studies, yet the X constitutes 5% of the genome and is an important contributor to disease, often in a sex-specific manner.

Regulation of toxic shock syndrome toxin-1 by the accessory gene regulator in Staphylococcus aureus is mediated by the repressor of toxins.

Toxic shock syndrome toxin-1 (TSST-1) is a superantigen (SAg) produced by Staphylococcus aureus thought to be responsible for essentially all cases of menstrual-associated toxic shock syndrome (TSS). As a potent exotoxin, it is not surprising that S. aureus has evolved multiple systems to control expression of TSST-1. Although the accessory gene regulator (Agr) system is recognized to enhance TSST-1 expression, how Agr regulates TSST-1 is unclear. Using an agr-null mutant, complementation experiments demonstrated that Agr controls TSST-1 expression through the activity of the RNAIII effector molecule. RNAIII can repress translation of the repressor of toxins (Rot) regulator, and deletion of rot increased expression of TSST-1 during the exponential phase of growth. Deletion of agr did not affect rot transcription, but did result in overexpression of the Rot protein, and Rot was also shown to bind and positively regulate the rot promoter. Overexpression of Rot dramatically repressed TSST-1, and Rot bound directly to the TSST-1 promoter. Deletion of both agr and rot in S. aureus returned TSST-1 expression to wild-type levels. This work demonstrates that Agr, although widely considered to be an inducer of TSST-1, has evolved in combination with Rot, to restrict the expression of this potent SAg.

Role of the regulatory genes SEF1, VMA1 and SFU1 in riboflavin synthesis in the flavinogenic yeast Candida famata (Candida flareri).

Riboflavin or vitamin B2 is an essential dietary component for humans and animals that is the precursor of flavin coenzymes flavin mononucleotide and flavin adenine dinucleotide involved in numerous enzymatic reactions. The flavinogenic yeast Candida famata overproduces riboflavin under iron starvation; however, regulation of this process is poorly understood. Regulatory gene SEF1 encoding transcription activator has been identified. Its deletion blocks yeast ability to overproduce riboflavin under iron starvation. It was shown here that the SEF1 promoters from other flavinogenic (Candida albicans) and non-flavinogenic (Candida tropicalis) yeasts fused with the open reading frame (ORF) of SEF1 gene from C. famata are able to restore riboflavin oversynthesis in sef1Δ mutants. It is known that in the pathogenic flavinogenic yeast C. albicans, Sfu1 (GATA-type transcription factor) represses SEF1. Here, we found that deletion of SFU1 gene in wild-type C. famata leads to riboflavin oversynthesis. Moreover, it was shown that disruption of VMA1 gene (coding for vacuolar ATPase subunit A) also results in riboflavin oversynthesis in C. famata.

The chromatin-remodeling enzyme BRG1 coordinates CIITA induction through many interdependent distal enhancers.

The chromatin-remodeling enzyme BRG1 is critical for interferon-gamma (IFN-gamma)-mediated gene induction. Promoter-proximal elements are sufficient to mediate BRG1 dependency at some IFN-gamma targets. In contrast, we show here that at CIITA, which encodes the 'master regulator' of induction of major histocompatibility complex class II, distal elements conferred BRG1 dependency. At the uninduced locus, many sites formed BRG1-independent loops. One loop juxtaposed a far downstream element adjacent to a far upstream site. Notably, BRG1 was recruited to the latter site, which triggered the appearance of a histone 'mark' linked to activation. This subtle change was crucial, as subsequent IFN-gamma-induced recruitment of the transcription factors STAT1, IRF1 and p300, as well as histone modifications, accessibility and additional loops, showed BRG1 dependency. Like BRG1, each remote element was critical for the induction of CIITA expression. Thus, BRG1 regulates CIITA through many interdependent remote enhancers, not through the promoter alone.

Expression of small RNAs of Bordetella pertussis colonizing murine tracheas.

We performed RNA sequencing on Bordetella pertussis, the causative agent of whooping cough, and identified nine novel small RNAs (sRNAs) that were transcribed during the bacterial colonization of murine tracheas. Among them, four sRNAs were more strongly expressed in vivo than in vitro. Moreover, the expression of eight sRNAs was not regulated by the BvgAS two-component system, which is the master regulator for the expression of genes contributing to the bacterial infection. The present results suggest a BvgAS-independent gene regulatory system involving the sRNAs that is active during B. pertussis infection.

fus/TLS orchestrates splicing of developmental regulators during gastrulation.

Here we investigated the function of the atypical RNA-binding protein fus/TLS (fused in sarcoma/translocated in sarcoma) during early frog development. We found that fus is necessary for proper mRNA splicing of a set of developmental regulatory genes during early frog development and gastrulation. Upon fus knockdown, embryos fail to gastrulate and show mesodermal differentiation defects that we connect to intron retention in fgf8 (fibroblast growth factor 8) and fgfr2 (fgf receptor 2) transcripts. During gastrulation, the animal and marginal regions dissociate, and we show that this is caused, at least in part, by intron retention in cdh1 transcripts. We confirm the specificity of splicing defects at a genomic level using analysis of RNA sequencing (RNA-seq) and show that 3%-5% of all transcripts display intron retention throughout the pre-mRNA. By analyzing gene ontology slim annotations, we show that the affected genes are enriched for developmental regulators and therefore represent a biologically coherent set of targets for fus regulation in embryogenesis. This shows that fus is central to embryogenesis and may provide information on its function in neurodegenerative disease.