Minor class splicing shapes the zebrafish transcriptome during development

Taxonomy of introns and the evolution of minor introns

Classification of introns, which is crucial to understanding their evolution and splicing, has historically been binary and has resulted in the naming of major and minor introns that are spliced by their namesake spliceosome. However, a broad range of intron consensus sequences exist, leading us to here reclassify introns as minor, minor-like, hybrid, major-like, major and non-canonical introns in 263 species across six eukaryotic supergroups. Through intron orthology analysis, we discovered that minor-like introns are a transitory node for intron conversion across evolution. Despite close resemblance of their consensus sequences to minor introns, these introns possess an AG dinucleotide at the -1 and -2 position of the 5' splice site, a salient feature of major introns. Through combined analysis of CoLa-seq, CLIP-seq for major and minor spliceosome components, and RNAseq from samples in which the minor spliceosome is inhibited we found that minor-like introns are also an intermediate class from a splicing mechanism perspective. Importantly, this analysis has provided insight into the sequence elements that have evolved to make minor-like introns amenable to recognition by both minor and major spliceosome components. We hope that this revised intron classification provides a new framework to study intron evolution and splicing.

Role of U11/U12 minor spliceosome gene ZCRB1 in Ciliogenesis and WNT Signaling

Despite the fact that 0.5% of human introns are processed by the U11/U12 minor spliceosome, the latter influences gene expression across multiple cellular processes. The ZCRB1 protein is a recently described core component of the U12 mono-snRNP minor spliceosome, but its functional significance to minor splicing, gene regulation, and biological signaling cascades is poorly understood. Using CRISPR-Cas9 and siRNA targeted knockout and knockdown strategies, we show that human cell lines with a partial reduction in ZCRB1 expression exhibit significant dysregulation of the splicing and expression of U12-type genes, primarily due to dysregulation of U12 mono-snRNA. RNA-Seq and targeted analyses of minor intron-containing genes indicate a downregulation in the expression of genes involved in ciliogenesis, and consequentially an upregulation in WNT signaling. Additionally, zcrb1 CRISPR-Cas12a knockdown in zebrafish embryos led to gross developmental and body axis abnormalities, disrupted ciliogenesis, and upregulated WNT signaling, complementing our human cell studies. This work highlights a conserved and essential biological role of the minor spliceosome in general, and the ZCRB1 protein specifically in cellular and developmental processes across species, shedding light on the multifaceted relationship between splicing regulation, ciliogenesis, and WNT signaling.

Minor Spliceosomal 65K/RNPC3 Interacts with ANKRD11 and Mediates HDAC3-Regulated Histone Deacetylation and Transcription

RNA splicing is crucial in the multilayer regulatory networks for gene expression, making functional interactions with DNA- and other RNA-processing machineries in the nucleus. However, these established couplings are all major spliceosome-related; whether the minor spliceosome is involved remains unclear. Here, through affinity purification using Drosophila lysates, an interaction is identified between the minor spliceosomal 65K/RNPC3 and ANKRD11, a cofactor of histone deacetylase 3 (HDAC3). Using a CRISPR/Cas9 system, Deletion strains are constructed and found that both Dm65KΔ/Δ and Dmankrd11Δ/Δ mutants have reduced histone deacetylation at Lys9 of histone H3 (H3K9) and Lys5 of histone H4 (H4K5) in their heads, exhibiting various neural-related defects. The 65K-ANKRD11 interaction is also conserved in human cells, and the HsANKRD11 middle-uncharacterized domain mediates Hs65K association with HDAC3. Cleavage under targets and tagmentation (CUT&Tag) assays revealed that HsANKRD11 is a bridging factor, which facilitates the synergistic common chromatin-binding of HDAC3 and Hs65K. Knockdown (KD) of HsANKRD11 simultaneously decreased their common binding, resulting in reduced deacetylation of nearby H3K9. Ultimately, this study demonstrates that expression changes of many genes caused by HsANKRD11-KD are due to the decreased common chromatin-binding of HDAC3 and Hs65K and subsequently reduced deacetylation of H3K9, illustrating a novel and conserved coupling mechanism that links the histone deacetylation with minor spliceosome for the regulation of gene expression.

Congenital hydrocephalus: new Mendelian mutations and evidence for oligogenic inheritance

BACKGROUND

Congenital hydrocephalus is characterized by ventriculomegaly, defined as a dilatation of cerebral ventricles, and thought to be due to impaired cerebrospinal fluid (CSF) homeostasis. Primary congenital hydrocephalus is a subset of cases with prenatal onset and absence of another primary cause, e.g., brain hemorrhage. Published series report a Mendelian cause in only a minority of cases. In this study, we analyzed exome data of PCH patients in search of novel causal genes and addressed the possibility of an underlying oligogenic mode of inheritance for PCH.

MATERIALS AND METHODS

We sequenced the exome in 28 unrelated probands with PCH, 12 of whom from families with at least two affected siblings and 9 of whom consanguineous, thereby increasing the contribution of genetic causes. Patient exome data were first analyzed for rare (MAF < 0.005) transmitted or de novo variants. Population stratification of unrelated PCH patients and controls was determined by principle component analysis, and outliers identified using Mahalanobis distance 5% as cutoff. Patient and control exome data for genes biologically related to cilia (SYScilia database) were analyzed by mutation burden test.

RESULTS

In 18% of probands, we identify a causal (pathogenic or likely pathogenic) variant of a known hydrocephalus gene, including genes for postnatal, syndromic hydrocephalus, not previously reported in isolated PCH. In a further 11%, we identify mutations in novel candidate genes. Through mutation burden tests, we demonstrate a significant burden of genetic variants in genes coding for proteins of the primary cilium in PCH patients compared to controls.

CONCLUSION

Our study confirms the low contribution of Mendelian mutations in PCH and reports PCH as a phenotypic presentation of some known genes known for syndromic, postnatal hydrocephalus. Furthermore, this study identifies novel Mendelian candidate genes, and provides evidence for oligogenic inheritance implicating primary cilia in PCH.

ahctf1 and kras mutations combine to amplify oncogenic stress and restrict liver overgrowth in a zebrafish model of hepatocellular carcinoma

The nucleoporin (NUP) ELYS, encoded by AHCTF1, is a large multifunctional protein with essential roles in nuclear pore assembly and mitosis. Using both larval and adult zebrafish models of hepatocellular carcinoma (HCC), in which the expression of an inducible mutant kras transgene (kras^G12V) drives hepatocyte-specific hyperplasia and liver enlargement, we show that reducing ahctf1 gene dosage by 50% markedly decreases liver volume, while non-hyperplastic tissues are unaffected. We demonstrate that in the context of cancer, ahctf1 heterozygosity impairs nuclear pore formation, mitotic spindle assembly, and chromosome segregation, leading to DNA damage and activation of a Tp53-dependent transcriptional programme that induces cell death and cell cycle arrest. Heterozygous expression of both ahctf1 and ranbp2 (encoding a second nucleoporin), or treatment of heterozygous ahctf1 larvae with the nucleocytoplasmic transport inhibitor, Selinexor, completely blocks kras^G12V-driven hepatocyte hyperplasia. Gene expression analysis of patient samples in the liver hepatocellular carcinoma (LIHC) dataset in The Cancer Genome Atlas shows that high expression of one or more of the transcripts encoding the 10 components of the NUP107-160 subcomplex, which includes AHCTF1, is positively correlated with worse overall survival. These results provide a strong and feasible rationale for the development of novel cancer therapeutics that target ELYS function and suggest potential avenues for effective combinatorial treatments.

Roles of minor spliceosome in intron recognition and the convergence with the better understood major spliceosome

Catalyzed by spliceosomes in the nucleus, RNA splicing removes intronic sequences from precursor RNAs in eukaryotes to generate mature RNA, which also significantly increases proteome complexity and fine-tunes gene expression. Most metazoans have two coexisting spliceosomes; the major spliceosome, which removes >99.5% of introns, and the minor spliceosome, which removes far fewer introns (only 770 at present have been predicted in the human genome). Both spliceosomes are large and dynamic machineries, each consisting of five small nuclear RNAs (snRNAs) and more than 100 proteins. However, the dynamic assembly, catalysis, and protein composition of the minor spliceosome are still poorly understood. With different splicing signals, minor introns are rare and usually distributed alone and flanked by major introns in genes, raising questions of how they are recognized by the minor spliceosome and how their processing deals with the splicing of neighboring major introns. Due to large numbers of introns and close similarities between the two machinery, cooperative, and competitive recognition by the two spliceosomes has been investigated. Functionally, many minor-intron-containing genes are evolutionarily conserved and essential. Mutations in the minor spliceosome exhibit a variety of developmental defects in plants and animals and are linked to numerous human diseases. Here, we review recent progress in the understanding of minor splicing, compare currently known components of the two spliceosomes, survey minor introns in a wide range of organisms, discuss cooperation and competition of the two spliceosomes in splicing of minor-intron-containing genes, and contributions of minor splicing mutations in development and diseases. This article is categorized under: RNA Processing > Processing of Small RNAs RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Addition of an affected family member to a previously ascertained autosomal recessive nonsyndromic hearing loss pedigree and systematic phenotype-genotype analysis of splice-site variants in MYO15A

Pathogenic variants in MYO15A are known to cause autosomal recessive nonsyndromic hearing loss (ARNSHL), DFNB3. We have previously reported on one ARNSHL family including two affected siblings and identified MYO15A c.5964+3G > A and c.8375 T > C (p.Val2792Ala) as the possible deafness-causing variants. Eight year follow up identified one new affected individual in this family, who also showed congenital, severe to profound sensorineural hearing loss. By whole exome sequencing, we identified a new splice-site variant c.5531+1G > C (maternal allele), in a compound heterozygote with previously identified missense variant c.8375 T > C (p.Val2792Ala) (paternal allele) in MYO15A as the disease-causing variants. The new affected individual underwent unilateral cochlear implantation at the age of 1 year, and 5 year follow-up showed satisfactory speech and language outcomes. Our results further indicate that MYO15A-associated hearing loss is good candidates for cochlear implantation, which is in accordance with previous report. In light of our findings and review of the literatures, 58 splice-site variants in MYO15A are correlated with a severe deafness phenotype, composed of 46 canonical splice-site variants and 12 non-canonical splice-site variants.

Emerging Roles of RNA-Binding Proteins in Inner Ear Hair Cell Development and Regeneration

RNA-binding proteins (RBPs) regulate gene expression at the post-transcriptional level. They play major roles in the tissue- and stage-specific expression of protein isoforms as well as in the maintenance of protein homeostasis. The inner ear is a bi-functional organ, with the cochlea and the vestibular system required for hearing and for maintaining balance, respectively. It is relatively well documented that transcription factors and signaling pathways are critically involved in the formation of inner ear structures and in the development of hair cells. Accumulating evidence highlights emerging functions of RBPs in the post-transcriptional regulation of inner ear development and hair cell function. Importantly, mutations of splicing factors of the RBP family and defective alternative splicing, which result in inappropriate expression of protein isoforms, lead to deafness in both animal models and humans. Because RBPs are critical regulators of cell proliferation and differentiation, they present the potential to promote hair cell regeneration following noise- or ototoxin-induced damage through mitotic and non-mitotic mechanisms. Therefore, deciphering RBP-regulated events during inner ear development and hair cell regeneration can help define therapeutic strategies for treatment of hearing loss. In this review, we outline our evolving understanding of the implications of RBPs in hair cell formation and hearing disease with the aim of promoting future research in this field.

The emerging significance of splicing in vertebrate development

Splicing is a crucial regulatory node of gene expression that has been leveraged to expand the proteome from a limited number of genes. Indeed, the vast increase in intron number that accompanied vertebrate emergence might have aided the evolution of developmental and organismal complexity. Here, we review how animal models for core spliceosome components have provided insights into the role of splicing in vertebrate development, with a specific focus on neuronal, neural crest and skeletal development. To this end, we also discuss relevant spliceosomopathies, which are developmental disorders linked to mutations in spliceosome subunits. Finally, we discuss potential mechanisms that could underlie the tissue-specific phenotypes often observed upon spliceosome inhibition and identify gaps in our knowledge that, we hope, will inspire further research.

Genetic analysis of human RNA binding motif protein 48 (RBM48) reveals an essential role in U12-type intron splicing

U12-type or minor introns are found in most multicellular eukaryotes and constitute ∼0.5% of all introns in species with a minor spliceosome. Although the biological significance for the evolutionary conservation of U12-type introns is debated, mutations disrupting U12 splicing cause developmental defects in both plants and animals. In human hematopoietic stem cells, U12 splicing defects disrupt proper differentiation of myeloid lineages and are associated with myelodysplastic syndrome, predisposing individuals to acute myeloid leukemia. Mutants in the maize ortholog of RNA binding motif protein 48 (RBM48) have aberrant U12-type intron splicing. Human RBM48 was recently purified biochemically as part of the minor spliceosome and shown to recognize the 5' end of the U6atac snRNA. In this report, we use CRISPR/Cas9-mediated ablation of RBM48 in human K-562 cells to show the genetic function of RBM48. RNA-seq analysis comparing wild-type and mutant K-562 genotypes found that 48% of minor intron-containing genes have significant U12-type intron retention in RBM48 mutants. Comparing these results to maize rbm48 mutants defined a subset of minor intron-containing genes disrupted in both species. Mutations in the majority of these orthologous minor intron-containing genes have been reported to cause developmental defects in both plants and animals. Our results provide genetic evidence that the primary defect of human RBM48 mutants is aberrant U12-type intron splicing, while a comparison of human and maize RNA-seq data identifies candidate genes likely to mediate mutant phenotypes of U12-type splicing defects.

Zrsr2 Is Essential for the Embryonic Development and Splicing of Minor Introns in RNA and Protein Processing Genes in Zebrafish

ZRSR2 (zinc finger CCCH-type, RNA binding motif and serine/arginine rich 2) is an essential splicing factor involved in 3' splice-site recognition as a component of both the major and minor spliceosomes that mediate the splicing of U2-type (major) and U12-type (minor) introns, respectively. Studies of ZRSR2-depleted cell lines and ZRSR2-mutated patient samples revealed its essential role in the U12-dependent minor spliceosome. However, the role of ZRSR2 during embryonic development is not clear, as its function is compensated for by Zrsr1 in mice. Here, we utilized the zebrafish model to investigate the role of zrsr2 during embryonic development. Using CRISPR/Cas9 technology, we generated a zrsr2-knockout zebrafish line, termed zrsr2hg129/hg129 (p.Trp167Argfs*9) and examined embryo development in the homozygous mutant embryos. zrsr2hg129/hg129 embryos displayed multiple developmental defects starting at 4 days post fertilization (dpf) and died after 8 dpf, suggesting that proper Zrsr2 function is required during embryonic development. The global transcriptome analysis of 3 dpf zrsr2hg129/hg129 embryos revealed that the loss of Zrsr2 results in the downregulation of essential metabolic pathways and the aberrant retention of minor introns in about one-third of all minor intron-containing genes in zebrafish. Overall, our study has demonstrated that the role of Zrsr2 as a component of the minor spliceosome is conserved and critical for proper embryonic development in zebrafish.

Zrsr2 and functional U12-dependent spliceosome are necessary for follicular development

ZRSR2 is a splicing factor involved in recognition of 3'-intron splice sites that is frequently mutated in myeloid malignancies and several tumors; however, the role of mutations of Zrsr2 in other tissues has not been analyzed. To explore the biological role of ZRSR2, we generated three Zrsr2 mutant mouse lines. All Zrsr2 mutant lines exhibited blood cell anomalies, and in two lines, oogenesis was blocked at the secondary follicle stage. RNA-seq of Zrsr2 mu secondary follicles showed aberrations in gene expression and showed altered alternative splicing (AS) events involving enrichment of U12-type intron retention (IR), supporting the functional Zrsr2 action in minor spliceosomes. IR events were preferentially associated with centriole replication, protein phosphorylation, and DNA damage checkpoint. Notably, we found alterations in AS events of 50 meiotic genes. These results indicate that ZRSR2 mutations alter splicing mainly in U12-type introns, which may affect peripheral blood cells, and impede oogenesis and female fertility.

Alternative RNA splicing in tumour heterogeneity, plasticity and therapy

Alternative splicing is a process by which a single gene is able to encode multiple different protein isoforms. It is regulated by the inclusion or exclusion of introns and exons that are joined in different patterns prior to protein translation, thus enabling transcriptomic and proteomic diversity. It is now widely accepted that alternative splicing is dysregulated across nearly all cancer types. This widespread dysregulation means that nearly all cellular processes are affected - these include processes synonymous with the hallmarks of cancer - evasion of apoptosis, tissue invasion and metastasis, altered cellular metabolism, genome instability and drug resistance. Emerging evidence indicates that the dysregulation of alternative splicing also promotes a permissive environment for increased tumour heterogeneity and cellular plasticity. These are fundamental regulators of a patient's response to therapy. In this Review, we introduce the mechanisms of alternative splicing and the role of aberrant splicing in cancer, with particular focus on newfound evidence of alternative splicing promoting tumour heterogeneity, cellular plasticity and altered metabolism. We discuss recent in vivo models generated to study alternative splicing and the importance of these for understanding complex tumourigenic processes. Finally, we review the effects of alternative splicing on immune evasion, cell death and genome instability, and how targeting these might enhance therapeutic efficacy.

Minor Intron Splicing from Basic Science to Disease

Pre-mRNA splicing is an essential step in gene expression and is catalyzed by two machineries in eukaryotes: the major (U2 type) and minor (U12 type) spliceosomes. While the majority of introns in humans are U2 type, less than 0.4% are U12 type, also known as minor introns (mi-INTs), and require a specialized spliceosome composed of U11, U12, U4atac, U5, and U6atac snRNPs. The high evolutionary conservation and apparent splicing inefficiency of U12 introns have set them apart from their major counterparts and led to speculations on the purpose for their existence. However, recent studies challenged the simple concept of mi-INTs splicing inefficiency due to low abundance of their spliceosome and confirmed their regulatory role in alternative splicing, significantly impacting the expression of their host genes. Additionally, a growing list of minor spliceosome-associated diseases with tissue-specific pathologies affirmed the importance of minor splicing as a key regulatory pathway, which when deregulated could lead to tissue-specific pathologies due to specific alterations in the expression of some minor-intron-containing genes. Consequently, uncovering how mi-INTs splicing is regulated in a tissue-specific manner would allow for better understanding of disease pathogenesis and pave the way for novel therapies, which we highlight in this review.

Multi-Dimensional Transcriptome Analysis Reveals Modulation of Cholesterol Metabolism as Highly Integrated Response to Brain Injury

Zebrafish is an attractive model to investigate regeneration of the nervous system. Despite major progress in our understanding of the underlying processes, the transcriptomic changes are largely unknown. We carried out a computational analysis of the transcriptome of the regenerating telencephalon integrating changes in the expression of mRNAs, their splice variants and investigated the putative role of regulatory RNAs in the modulation of these transcriptional changes. Profound changes in the expression of genes and their splice variants engaged in many distinct processes were observed. Differential transcription and splicing are important processes in response to injury of the telencephalon. As exemplified by the coordinated regulation of the cholesterol synthesizing enzymes and transporters, the genome responded to injury of the telencephalon in a multi-tiered manner with distinct and interwoven changes in expression of enzymes, transporters and their regulatory molecules. This coordinated genomic response involved a decrease of the mRNA of the key transcription factor SREBF2, induction of microRNAs (miR-182, miR-155, miR-146, miR-31) targeting cholesterol genes, shifts in abundance of splice variants as well as regulation of long non-coding RNAs. Cholesterol metabolism appears to be switched from synthesis to relocation of cholesterol. Based on our in silico analyses, this switch involves complementary and synergistic inputs by different regulatory principles. Our studies suggest that adaptation of cholesterol metabolism is a key process involved in regeneration of the injured zebrafish brain.

Disruption of exon-bridging interactions between the minor and major spliceosomes results in alternative splicing around minor introns

Vertebrate genomes contain major (>99.5%) and minor (<0.5%) introns that are spliced by the major and minor spliceosomes, respectively. Major intron splicing follows the exon-definition model, whereby major spliceosome components first assemble across exons. However, since most genes with minor introns predominately consist of major introns, formation of exon-definition complexes in these genes would require interaction between the major and minor spliceosomes. Here, we report that minor spliceosome protein U11-59K binds to the major spliceosome U2AF complex, thereby supporting a model in which the minor spliceosome interacts with the major spliceosome across an exon to regulate the splicing of minor introns. Inhibition of minor spliceosome snRNAs and U11-59K disrupted exon-bridging interactions, leading to exon skipping by the major spliceosome. The resulting aberrant isoforms contained a premature stop codon, yet were not subjected to nonsense-mediated decay, but rather bound to polysomes. Importantly, we detected elevated levels of these alternatively spliced transcripts in individuals with minor spliceosome-related diseases such as Roifman syndrome, Lowry–Wood syndrome and early-onset cerebellar ataxia. In all, we report that the minor spliceosome informs splicing by the major spliceosome through exon-definition interactions and show that minor spliceosome inhibition results in aberrant alternative splicing in disease.

Minor intron retention drives clonal hematopoietic disorders and diverse cancer predisposition

Most eukaryotes harbor two distinct pre-mRNA splicing machineries: the major spliceosome, which removes >99% of introns, and the minor spliceosome, which removes rare, evolutionarily conserved introns. Although hypothesized to serve important regulatory functions, physiologic roles of the minor spliceosome are not well understood. For example, the minor spliceosome component ZRSR2 is subject to recurrent, leukemia-associated mutations, yet functional connections among minor introns, hematopoiesis and cancers are unclear. Here, we identify that impaired minor intron excision via ZRSR2 loss enhances hematopoietic stem cell self-renewal. CRISPR screens mimicking nonsense-mediated decay of minor intron-containing mRNA species converged on LZTR1, a regulator of RAS-related GTPases. LZTR1 minor intron retention was also discovered in the RASopathy Noonan syndrome, due to intronic mutations disrupting splicing and diverse solid tumors. These data uncover minor intron recognition as a regulator of hematopoiesis, noncoding mutations within minor introns as potential cancer drivers and links among ZRSR2 mutations, LZTR1 regulation and leukemias.

Splicing factor DHX15 affects tp53 and mdm2 expression via alternate splicing and promoter usage

DHX15, a DEAH box containing RNA helicase, is a splicing factor required for the last step of splicing. Recent studies identified a recurrent mutational hotspot, R222G, in DHX15 in ∼ 6% of acute myeloid leukemia (AML) patients that carry the fusion protein RUNX1-RUNX1T1 produced by t (8;21) (q22;q22). Studies using yeast mutants showed that substitution of G for the residue equivalent to R222 leads to loss of its helicase function, suggesting that it is a loss-of-function mutation. To elucidate the role of DHX15 during development, we established the first vertebrate knockout model with CRISPR/Cas9 in zebrafish. Our data showed that dhx15 expression is enriched in the brain, eyes, pectoral fin primordia, liver and intestinal bulb during embryonic development. Dhx15 deficiency leads to pleiotropic morphological phenotypes in homozygous mutant embryos starting at 3 days post fertilization (dpf) that result in lethality by 7 dpf, revealing an essential role during embryonic development. RNA-seq analysis suggested important roles of Dhx15 in chromatin and nucleosome assembly and regulation of the Mdm2-p53 pathway. Interestingly, exons corresponding to the alternate transcriptional start sites for tp53 and mdm2 were preferentially expressed in the mutant embryos, leading to significant upregulation of their alternate isoforms, Δ113p53 (orthologous to Δ133p53 isoform in human) and mdm2-P2 (isoform using distal promoter P2), respectively. We speculate that these alterations in the Mdm2-p53 pathway contribute to the development of AML in patients with t(8;21) and somatically mutated DHX15.

Minor snRNA gene delivery improves the loss of proprioceptive synapses on SMA motor neurons

Spinal muscular atrophy (SMA) is an inherited neuromuscular disorder caused by reduced expression of the survival motor neuron (SMN) protein. SMN has key functions in multiple RNA pathways, including the biogenesis of small nuclear ribonucleoproteins that are essential components of both major (U2-dependent) and minor (U12-dependent) spliceosomes. Here we investigated the specific contribution of U12 splicing dysfunction to SMA pathology through selective restoration of this RNA pathway in mouse models of varying phenotypic severity. We show that virus-mediated delivery of minor snRNA genes specifically improves select U12 splicing defects induced by SMN deficiency in cultured mammalian cells, as well as in the spinal cord and dorsal root ganglia of SMA mice without increasing SMN expression. This approach resulted in a moderate amelioration of several parameters of the disease phenotype in SMA mice, including survival, weight gain, and motor function. Importantly, minor snRNA gene delivery improved aberrant splicing of the U12 intron-containing gene Stasimon and rescued the severe loss of proprioceptive sensory synapses on SMA motor neurons, which are early signatures of motor circuit dysfunction in mouse models. Taken together, these findings establish the direct contribution of U12 splicing dysfunction to synaptic deafferentation and motor circuit pathology in SMA.

Minor Splicing Factors Zrsr1 and Zrsr2 Are Essential for Early Embryo Development and 2-Cell-Like Conversion

Minor splicing plays an important role in vertebrate development. Zrsr1 and Zrsr2 paralog genes have essential roles in alternative splicing, mainly participating in the recognition of minor (U12) introns. To further explore their roles during early embryo development, we produced Zrsr1mu and Zrsr2mu mutant mice, containing truncating mutations within the second zinc finger domain. Both homozygous mutant mice were viable with a normal lifespan. When we crossed a homozygous Zrsr2mu/mu female with Zrsr1mu/mu male, the double heterozygotes were non-viable, giving rise to embryos that stopped developing mainly between the 2- and 4-cell stages, just after zygotic gene activation. RNA-seq analysis of Zrsr1/2mu 2-cell embryos showed altered gene and isoform expression of thousands of genes enriched in gene ontology terms and biological pathways related to ribosome, RNA transport, spliceosome, and essential zygotic gene activation steps. Alternative splicing was analyzed, showing a significant increase in intron retention in both U2 and U12 intron-containing genes related to cell cycle and mitotic nuclear division. Remarkably, both Zrsr1 and Zrsr2 were required for the conversion of mouse-induced pluripotent stem cells into 2C-like cells. According to our results, Zrsr1 or Zrsr2 are necessary for ZGA and both are indispensable for the conversion of induced pluripotent stem cells into 2C-like cells.

The Molecular Basis of Congenital Hypopituitarism and Related Disorders

CONTEXT

Congenital hypopituitarism (CH) is characterized by the presence of deficiencies in one or more of the 6 anterior pituitary (AP) hormones secreted from the 5 different specialized cell types of the AP. During human embryogenesis, hypothalamo-pituitary (HP) development is controlled by a complex spatio-temporal genetic cascade of transcription factors and signaling molecules within the hypothalamus and Rathke's pouch, the primordium of the AP.

EVIDENCE ACQUISITION

This mini-review discusses the genes and pathways involved in HP development and how mutations of these give rise to CH. This may present in the neonatal period or later on in childhood and may be associated with craniofacial midline structural abnormalities such as cleft lip/palate, visual impairment due to eye abnormalities such as optic nerve hypoplasia (ONH) and microphthalmia or anophthalmia, or midline forebrain neuroradiological defects including agenesis of the septum pellucidum or corpus callosum or the more severe holoprosencephaly.

EVIDENCE SYNTHESIS

Mutations give rise to an array of highly variable disorders ranging in severity. There are many known causative genes in HP developmental pathways that are routinely screened in CH patients; however, over the last 5 years this list has rapidly increased due to the identification of variants in new genes and pathways of interest by next-generation sequencing.

CONCLUSION

The majority of patients with these disorders do not have an identified molecular basis, often making management challenging. This mini-review aims to guide clinicians in making a genetic diagnosis based on patient phenotype, which in turn may impact on clinical management.

Spliceosomopathies and neurocristopathies: Two sides of the same coin?

Mutations in core components of the spliceosome are responsible for a group of syndromes collectively known as spliceosomopathies. Patients exhibit microcephaly, micrognathia, malar hypoplasia, external ear anomalies, eye anomalies, psychomotor delay, intellectual disability, limb, and heart defects. Craniofacial malformations in these patients are predominantly found in neural crest cells-derived structures of the face and head. Mutations in eight genes SNRPB, RNU4ATAC, SF3B4, PUF60, EFTUD2, TXNL4, EIF4A3, and CWC27 are associated with craniofacial spliceosomopathies. In this review, we provide a brief description of the normal development of the head and the face and an overview of mutations identified in genes associated with craniofacial spliceosomopathies. We also describe a model to explain how and when these mutations are most likely to impact neural crest cells. We speculate that mutations in a subset of core splicing factors lead to disrupted splicing in neural crest cells because these cells have increased sensitivity to inefficient splicing. Hence, disruption in splicing likely activates a cellular stress response that includes increased skipping of regulatory exons in genes such as MDM2 and MDM4, key regulators of P53. This would result in P53-associated death of neural crest cells and consequently craniofacial malformations associated with spliceosomopathies.

An Integrated Model of Minor Intron Emergence and Conservation

Minor introns constitute <0.5% of the introns in the human genome and have remained an enigma since their discovery. These introns are removed by a distinct splicing complex, the minor spliceosome. Both are ancient, tracing back to the last eukaryotic common ancestor (LECA), which is reflected by minor intron enrichment in specific gene families, such as the mitogen activated-protein kinase kinases, voltage-gated sodium and calcium ion channels, and E2F transcription factors. Most minor introns occur as single introns in genes with predominantly major introns. Due to this organization, minor intron-containing gene (MIG) expression requires the coordinated action of two spliceosomes, which increases the probability of missplicing. Thus, one would expect loss of minor introns via purifying selection. This has resulted in complete minor intron loss in at least nine eukaryotic lineages. However, minor introns are highly conserved in land plants and metazoans, where their importance is underscored by embryonic lethality when the minor spliceosome is inactivated. Conditional inactivation of the minor spliceosome has shown that rapidly dividing progenitor cells are highly sensitive to minor spliceosome loss. Indeed, we found that MIGs were significantly enriched in a screen for genes essential for survival in 341 cycling cell lines. Here, we propose that minor introns inserted randomly into genes in LECA or earlier and were subsequently conserved in genes crucial for cycling cell survival. We hypothesize that the essentiality of MIGs allowed minor introns to endure through the unicellularity of early eukaryotic evolution. Moreover, we identified 59 MIGs that emerged after LECA, and that many of these are essential for cycling cell survival, reinforcing our essentiality model for MIG conservation. This suggests that minor intron emergence is dynamic across eukaryotic evolution, and that minor introns should not be viewed as molecular fossils. We also posit that minor intron splicing was co-opted in multicellular evolution as a regulatory switch for en masse control of MIG expression and the biological processes they regulate. Specifically, this mode of regulation could control cell proliferation and thus body size, an idea supported by domestication syndrome, wherein MIGs are enriched in common candidate animal domestication genes.

Impaired Spermatogenesis, Muscle, and Erythrocyte Function in U12 Intron Splicing-Defective Zrsr1 Mutant Mice

The U2AF35-like ZRSR1 has been implicated in the recognition of 3' splice site during spliceosome assembly, but ZRSR1 knockout mice do not show abnormal phenotypes. To analyze ZRSR1 function and its precise role in RNA splicing, we generated ZRSR1 mutant mice containing truncating mutations within its RNA-recognition motif. Homozygous mutant mice exhibited severe defects in erythrocytes, muscle stretch, and spermatogenesis, along with germ cell sloughing and apoptosis, ultimately leading to azoospermia and male sterility. Testis RNA sequencing (RNA-seq) analyses revealed increased intron retention of both U2- and U12-type introns, including U12-type intron events in genes with key functions in spermatogenesis and spermatid development. Affected U2 introns were commonly found flanking U12 introns, suggesting functional cross-talk between the two spliceosomes. The splicing and tissue defects observed in mutant mice attributed to ZRSR1 loss of function suggest a physiological role for this factor in U12 intron splicing.

New insights into minor splicing-a transcriptomic analysis of cells derived from TALS patients

Minor intron splicing plays a central role in human embryonic development and survival. Indeed, biallelic mutations in RNU4ATAC, transcribed into the minor spliceosomal U4atac snRNA, are responsible for three rare autosomal recessive multimalformation disorders named Taybi-Linder (TALS/MOPD1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes, which associate numerous overlapping signs of varying severity. Although RNA-seq experiments have been conducted on a few RFMN patient cells, none have been performed in TALS, and more generally no in-depth transcriptomic analysis of the ∼700 human genes containing a minor (U12-type) intron had been published as yet. We thus sequenced RNA from cells derived from five skin, three amniotic fluid, and one blood biosamples obtained from seven unrelated TALS cases and from age- and sex-matched controls. This allowed us to describe for the first time the mRNA expression and splicing profile of genes containing U12-type introns, in the context of a functional minor spliceosome. Concerning RNU4ATAC-mutated patients, we show that as expected, they display distinct U12-type intron splicing profiles compared to controls, but that rather unexpectedly mRNA expression levels are mostly unchanged. Furthermore, although U12-type intron missplicing concerns most of the expressed U12 genes, the level of U12-type intron retention is surprisingly low in fibroblasts and amniocytes, and much more pronounced in blood cells. Interestingly, we found several occurrences of introns that can be spliced using either U2, U12, or a combination of both types of splice site consensus sequences, with a shift towards splicing using preferentially U2 sites in TALS patients' cells compared to controls.

Minor intron splicing revisited: identification of new minor intron-containing genes and tissue-dependent retention and alternative splicing of minor introns

BACKGROUND

Mutations in minor spliceosome components such as U12 snRNA (cerebellar ataxia) and U4atac snRNA (microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1)) result in tissue-specific symptoms. Given that the minor spliceosome is ubiquitously expressed, we hypothesized that these restricted phenotypes might be caused by the tissue-specific regulation of the minor spliceosome targets, i.e. minor intron-containing genes (MIGs). The current model of inefficient splicing is thought to apply to the regulation of the ~ 500 MIGs identified in the U12DB. However this database was created more than 10 years ago. Therefore, we first wanted to revisit the classification of minor introns in light of the most recent reference genome. We then sought to address specificity of MIG expression, minor intron retention, and alternative splicing (AS) across mouse and human tissues.

RESULTS

We employed position-weight matrices to obtain a comprehensive updated list of minor introns, consisting of 722 mouse and 770 human minor introns. These can be found in the Minor Intron DataBase (MIDB). Besides identification of 99% of the minor introns found in the U12DB, we also discovered ~ 150 new MIGs. We then analyzed the RNAseq data from eleven different mouse tissues, which revealed tissue-specific MIG expression and minor intron retention. Additionally, many minor introns were efficiently spliced compared to their flanking major introns. Finally, we identified several novel AS events across minor introns in both mouse and human, which were also tissue-dependent. Bioinformatics analysis revealed that several of the AS events could result in the production of novel tissue-specific proteins. Moreover, like the major introns, we found that these AS events were more prevalent in long minor introns, while retention was favoured in shorter introns.

CONCLUSION

Here we show that minor intron splicing and AS across minor introns is a highly organised process that might be regulated in coordination with the major spliceosome in a tissue-specific manner. We have provided a framework to further study the impact of the minor spliceosome and the regulation of MIG expression. These findings may shed light on the mechanism underlying tissue-specific phenotypes in diseases associated with minor spliceosome inactivation. MIDB can be accessed at https://midb.pnb.uconn.edu .

Sex-Dimorphic Behavioral Alterations and Altered Neurogenesis in U12 Intron Splicing-Defective Zrsr1 Mutant Mice

Mutant mice with respect to the splicing factor Zrsr1 present altered spermatogenesis and infertility. To investigate whether Zrsr1 is involved in the homeostatic control that the hypothalamus exerts over reproductive functions, we first analyzed both differential gene and isoform expression and alternative splicing alterations in Zrsr1 mutant (Zrsr1^mu) hypothalamus; second, we analyzed the spontaneous and social behavior of Zrsr1^mu mice; and third, we analyzed adult cell proliferation and survival in the Zrsr1^mu hypothalamus. The Zrsr1^mu hypothalamus showed altered expression of genes and isoforms related to the glutathione metabolic process, synaptonemal complex assembly, mRNA transport, and altered splicing events involving the enrichment of U12-type intron retention (IR). Furthermore, increased IR in U12-containing genes related with the prolactin, progesterone, and gonadotropin-releasing hormone (GnRH) reproductive signaling pathway was observed. This was associated with a hyperactive phenotype in both males and females, with an anxious phenotype in females, and with increased social interaction in males, instead of the classical aggressive behavior. In addition, Zrsr1^mu females but not males exhibited reduced cell proliferation in both the hypothalamus and the subventricular zone. Overall, these results suggest that Zrsr1 expression and function are relevant to organization of the hypothalamic cell network controlling behavior.

Differential isoform expression and alternative splicing in sex determination in mice

Background

Alternative splicing (AS) may play an important role in gonadal sex determination (GSD) in mammals. The present study was designed to identify differentially expressed isoforms and AS modifications accompanying GSD in mice.

Results

Using deep RNA-sequencing, we performed a transcriptional analysis of XX and XY gonads during sex determination on embryonic days 11 (E11) and 12 (E12). Analysis of differentially expressed genes (DEG) identified hundreds of genes related to GSD and early sex differentiation that may represent good candidates for sex reversal. Expression at time point E11 in males was significantly enriched in RNA splicing and mRNA processing Gene Ontology terms. Differentially expressed isoform analysis identified hundreds of specific isoforms related to GSD, many of which showed no differences in the DEG analysis. Hundreds of AS events were identified as modified at E11 and E12. Female E11 gonads featured sex-biased upregulation of intron retention (in genes related to regulation of transcription, protein phosphorylation, protein transport and mRNA splicing) and exon skipping (in genes related to chromatin repression) suggesting AS as a post-transcription mechanism that controls sex determination of the bipotential fetal gonad.

Conclusion

Our data suggests an important role of splicing regulatory mechanisms for sex determination in mice.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5572-x) contains supplementary material, which is available to authorized users.

RNA Binding Motif Protein 48 Is Required for U12 Splicing and Maize Endosperm Differentiation

The last eukaryotic common ancestor had two classes of introns that are still found in most eukaryotic lineages. Common U2-type and rare U12-type introns are spliced by the major and minor spliceosomes, respectively. Relatively few splicing factors have been shown to be specific to the minor spliceosome. We found that the maize (Zea mays) RNA binding motif protein 48 (RBM48) is a U12 splicing factor that functions to promote cell differentiation and repress cell proliferation. RBM48 is coselected with the U12 splicing factor, zinc finger CCCH-type, RNA binding motif, and Ser/Arg rich 2/Rough endosperm 3 (RGH3). Protein-protein interactions between RBM48, RGH3, and U2 Auxiliary Factor (U2AF) subunits suggest major and minor spliceosome factors required for intron recognition form complexes with RBM48. Human RBM48 interacts with armadillo repeat containing 7 (ARMC7). Maize RBM48 and ARMC7 have a conserved protein-protein interaction. These data predict that RBM48 is likely to function in U12 splicing throughout eukaryotes and that U12 splicing promotes endosperm cell differentiation in maize.

RNA splicing analysis in genomic medicine

High-throughput next-generation sequencing technologies have led to a rapid increase in the number of sequence variants identified in clinical practice via diagnostic genetic tests. Current bioinformatic analysis pipelines fail to take adequate account of the possible splicing effects of such variants, particularly where variants fall outwith canonical splice site sequences, and consequently the pathogenicity of such variants may often be missed. The regulation of splicing is highly complex and as a result, in silico prediction tools lack sufficient sensitivity and specificity for reliable use. Variants of all kinds can be linked to aberrant splicing in disease and the need for correct identification and diagnosis grows ever more crucial as novel splice-switching antisense oligonucleotide therapies start to enter clinical usage. RT-PCR provides a useful targeted assay of the splicing effects of identified variants, while minigene assays, massive parallel reporter assays and animal models can also be used for more detailed study of a particular splicing system, given enough time and resources. However, RNA-sequencing (RNA-seq) has the potential to be used as a rapid diagnostic tool in genomic medicine. By utilising data science approaches and machine learning, it may prove possible to finally understand and interpret the 'splicing code' and apply this knowledge in human disease diagnostics.

Early developmental arrest and impaired gastrointestinal homeostasis in U12-dependent splicing-defective Rnpc3-deficient mice

Splicing is an essential step in eukaryotic gene expression. While the majority of introns is excised by the U2-dependent, or major class, spliceosome, the appropriate expression of a very small subset of genes depends on U12-dependent, or minor class, splicing. The U11/U12 65K protein (hereafter 65K), encoded by RNPC3, is one of seven proteins that are unique to the U12-dependent spliceosome, and previous studies including our own have established that it plays a role in plant and vertebrate development. To pinpoint the impact of 65K loss during mammalian development and in adulthood, we generated germline and conditional Rnpc3-deficient mice. Homozygous Rnpc3^-/- embryos died prior to blastocyst implantation, whereas Rnpc3^+/- mice were born at the expected frequency, achieved sexual maturity, and exhibited a completely normal lifespan. Systemic recombination of conditional Rnpc3 alleles in adult (Rnpc3^lox/lox ) mice caused rapid weight loss, leukopenia, and degeneration of the epithelial lining of the entire gastrointestinal tract, the latter due to increased cell death and a reduction in cell proliferation. Accompanying this, we observed a loss of both 65K and the pro-proliferative phospho-ERK1/2 proteins from the stem/progenitor cells at the base of intestinal crypts. RT-PCR analysis of RNA extracted from purified preparations of intestinal epithelial cells with recombined Rnpc3^lox alleles revealed increased frequency of U12-type intron retention in all transcripts tested. Our study, using a novel conditional mouse model of Rnpc3 deficiency, establishes that U12-dependent splicing is not only important during development but is indispensable throughout life.

Minor spliceosome inactivation causes microcephaly, owing to cell cycle defects and death of self-amplifying radial glial cells

Mutation in minor spliceosome components is linked to the developmental disorder microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1). Here, we inactivated the minor spliceosome in the developing mouse cortex (pallium) by ablating Rnu11, which encodes the crucial minor spliceosome small nuclear RNA (snRNA) U11. Rnu11 conditional knockout mice were born with microcephaly, which was caused by the death of self-amplifying radial glial cells (RGCs), while intermediate progenitor cells and neurons were produced. RNA sequencing suggested that this cell death was mediated by upregulation of p53 (Trp53 - Mouse Genome Informatics) and DNA damage, which were both observed specifically in U11-null RGCs. Moreover, U11 loss caused elevated minor intron retention in genes regulating the cell cycle, which was consistent with fewer RGCs in S-phase and cytokinesis, alongside prolonged metaphase in RGCs. In all, we found that self-amplifying RGCs are the cell type most sensitive to loss of minor splicing. Together, these findings provide a potential explanation of how disruption of minor splicing might cause microcephaly in MOPD1.

Guided genetic screen to identify genes essential in the regeneration of hair cells and other tissues

Regenerative medicine holds great promise for both degenerative diseases and traumatic tissue injury which represent significant challenges to the health care system. Hearing loss, which affects hundreds of millions of people worldwide, is caused primarily by a permanent loss of the mechanosensory receptors of the inner ear known as hair cells. This failure to regenerate hair cells after loss is limited to mammals, while all other non-mammalian vertebrates tested were able to completely regenerate these mechanosensory receptors after injury. To understand the mechanism of hair cell regeneration and its association with regeneration of other tissues, we performed a guided mutagenesis screen using zebrafish lateral line hair cells as a screening platform to identify genes that are essential for hair cell regeneration, and further investigated how genes essential for hair cell regeneration were involved in the regeneration of other tissues. We created genetic mutations either by retroviral insertion or CRISPR/Cas9 approaches, and developed a high-throughput screening pipeline for analyzing hair cell development and regeneration. We screened 254 gene mutations and identified 7 genes specifically affecting hair cell regeneration. These hair cell regeneration genes fell into distinct and somewhat surprising functional categories. By examining the regeneration of caudal fin and liver, we found these hair cell regeneration genes often also affected other types of tissue regeneration. Therefore, our results demonstrate guided screening is an effective approach to discover regeneration candidates, and hair cell regeneration is associated with other tissue regeneration.

A study on zebrafish has genetically screened 254 genes and identified 7 genes implicated in the development and regeneration of hair cells and other tissues. Humans and other mammals cannot regrow hair cells—inner-ear sensory receptors that enable hearing—whereas non-mammalian vertebrates, including zebrafish, can regrow these following injury. Researchers from the United States, led by the National Institutes of Health’s Shawn Burgess, screened adult zebrafish for genes active during the regeneration of inner-ear epithelium. The researchers then produced zebrafish without these genes to study their functions. The studies tested 254 genes known to respond during regeneration, and identified seven specifically impacting regeneration. Most of these seven genes also functioned in liver and fin tissue regeneration. Understanding the mechanisms of these genes may enable future research into regenerative therapies in humans.

Mutations in the U11/U12-65K protein associated with isolated growth hormone deficiency lead to structural destabilization and impaired binding of U12 snRNA

Mutations in the components of the minor spliceosome underlie several human diseases. A subset of patients with isolated growth hormone deficiency (IGHD) harbors mutations in the RNPC3 gene, which encodes the minor spliceosome-specific U11/U12-65K protein. Although a previous study showed that IGHD patient cells have defects in U12-type intron recognition, the biochemical effects of these mutations on the 65K protein have not been characterized. Here, we show that a proline-to-threonine missense mutation (P474T) and a nonsense mutation (R502X) in the C-terminal RNA recognition motif (C-RRM) of the 65K protein impair the binding of 65K to U12 and U6atac snRNAs. We further show that the nonsense allele is targeted to the nonsense-mediated decay (NMD) pathway, but in an isoform-specific manner, with the nuclear-retained 65K long-3'UTR isoform escaping the NMD pathway. In contrast, the missense P474T mutation leads, in addition to the RNA-binding defect, to a partial defect in the folding of the C-RRM and reduced stability of the full-length protein, thus reducing the formation of U11/U12 di-snRNP complexes. We propose that both the C-RRM folding defect and NMD-mediated decrease in the levels of the U11/U12-65K protein reduce formation of the U12-type intron recognition complex and missplicing of a subset of minor introns leading to pituitary hypoplasia and a subsequent defect in growth hormone secretion.

Minor spliceosome and disease

The U12-dependent (minor) spliceosome excises a rare group of introns that are characterized by a highly conserved 5' splice site and branch point sequence. Several new congenital or somatic diseases have recently been associated with mutations in components of the minor spliceosome. A common theme in these diseases is the detection of elevated levels of transcripts containing U12-type introns, of which a subset is associated with other splicing defects. Here we review the present understanding of minor spliceosome diseases, particularly those associated with the specific components of the minor spliceosome. We also present a model for interpreting the molecular-level consequences of the different diseases.

Spliceosomal protein eftud2 mutation leads to p53-dependent apoptosis in zebrafish neural progenitors

Haploinsufficiency of EFTUD2 (Elongation Factor Tu GTP Binding Domain Containing 2) is linked to human mandibulofacial dysostosis, Guion-Almeida type (MFDGA), but the underlying cellular and molecular mechanisms remain to be addressed. We report here the isolation, cloning and functional analysis of the mutated eftud2 (snu114) in a novel neuronal mutant fn10a in zebrafish. This mutant displayed abnormal brain development with evident neuronal apoptosis while the development of other organs appeared less affected. Positional cloning revealed a nonsense mutation such that the mutant eftud2 mRNA encoded a truncated Eftud2 protein and was subjected to nonsense-mediated decay. Disruption of eftud2 led to increased apoptosis and mitosis of neural progenitors while it had little effect on differentiated neurons. Further RNA-seq and functional analyses revealed a transcriptome-wide RNA splicing deficiency and a large amount of intron-retaining and exon-skipping transcripts, which resulted in inadequate nonsense-mediated RNA decay and activation of the p53 pathway in fn10a mutants. Therefore, our study has established that eftud2 functions in RNA splicing during neural development and provides a suitable zebrafish model for studying the molecular pathology of the neurological disease MFDGA.

Aberrant splicing in maize rough endosperm3 reveals a conserved role for U12 splicing in eukaryotic multicellular development

RNA splicing of U12-type introns functions in human cell differentiation, but it is not known whether this class of introns has a similar role in plants. The maize ROUGH ENDOSPERM3 (RGH3) protein is orthologous to the human splicing factor, ZRSR2. ZRSR2 mutations are associated with myelodysplastic syndrome (MDS) and cause U12 splicing defects. Maize rgh3 mutants have aberrant endosperm cell differentiation and proliferation. We found that most U12-type introns are retained or misspliced in rgh3 Genes affected in rgh3 and ZRSR2 mutants identify cell cycle and protein glycosylation as common pathways disrupted. Transcripts with retained U12-type introns can be found in polysomes, suggesting that splicing efficiency can alter protein isoforms. The rgh3 mutant protein disrupts colocalization with a known ZRSR2-interacting protein, U2AF2. These results indicate conserved function for RGH3/ZRSR2 in U12 splicing and a deeply conserved role for the minor spliceosome to promote cell differentiation from stem cells to terminal fates.

Modeling intestinal disorders using zebrafish

Although the zebrafish was initially developed as a model system to study embryonic development, it has gained increasing attention as an advantageous system to investigate human diseases, including intestinal disorders. Zebrafish embryos develop rapidly, and their digestive system is fully functional and visible by 5days post fertilization. There is a large degree of homology between the intestine of zebrafish and higher vertebrate organisms in terms of its cellular composition and function as both a digestive and immune organ. Furthermore, molecular pathways regulating injury and immune responses are highly conserved. In this chapter, we provide an overview of studies addressing developmental and physiological processes relevant to human intestinal disease. These studies include those related to congenital disorders, host-microbiota interactions, inflammatory diseases, motility disorders, and intestinal cancer. We also highlight the utility of zebrafish to functionally validate candidate genes identified through mutational analyses and genome-wide association studies, and discuss methodologies to investigate the intestinal biology that are unique to zebrafish.

Evolutionarily conserved exon definition interactions with U11 snRNP mediate alternative splicing regulation on U11-48K and U11/U12-65K genes

Many splicing regulators bind to their own pre-mRNAs to induce alternative splicing that leads to formation of unstable mRNA isoforms. This provides an autoregulatory feedback mechanism that regulates the cellular homeostasis of these factors. We have described such an autoregulatory mechanism for two core protein components, U11-48K and U11/U12-65K, of the U12-dependent spliceosome. This regulatory system uses an atypical splicing enhancer element termed USSE (U11 snRNP-binding splicing enhancer), which contains two U12-type consensus 5' splice sites (5'ss). Evolutionary analysis of the USSE element from a large number of animal and plant species indicate that USSE sequence must be located 25-50 nt downstream from the target 3' splice site (3'ss). Together with functional evidence showing a loss of USSE activity when this distance is reduced and a requirement for RS-domain of U11-35K protein for 3'ss activation, our data suggests that U11 snRNP bound to USSE uses exon definition interactions for regulating alternative splicing. However, unlike standard exon definition where the 5'ss bound by U1 or U11 will be subsequently activated for splicing, the USSE element functions similarly as an exonic splicing enhancer and is involved only in upstream splice site activation but does not function as a splicing donor. Additionally, our evolutionary and functional data suggests that the function of the 5'ss duplication within the USSE elements is to allow binding of two U11/U12 di-snRNPs that stabilize each others' binding through putative mutual interactions.

Structural features important for the U12 snRNA binding and minor spliceosome assembly of Arabidopsis U11/U12-small nuclear ribonucleoproteins

Although seven proteins unique to U12 intron-specific minor spliceosomes, denoted as U11/U12-65K, -59K, -48K, -35K, -31K, -25K, and -20K, have been identified in humans and the roles of some of them have been demonstrated, the functional role of most of these proteins in plants is not understood. A recent study demonstrated that Arabidopsis U11/U12-65K is essential for U12 intron splicing and normal plant development. However, the structural features and sequence motifs important for 65 K binding to U12 snRNA and other spliceosomal proteins remain unclear. Here, we demonstrated by domain-deletion analysis that the C-terminal region of the 65 K protein bound specifically to the stem-loop III of U12 snRNA, whereas the N-terminal region of the 65 K protein was responsible for interacting with the 59 K protein. Analysis of the interactions between each snRNP protein using yeast two-hybrid analysis and in planta bimolecular fluorescence complementation and luciferase complementation imaging assays demonstrated that the core interactions among the 65 K, 59 K, and 48 K proteins were conserved between plants and animals, and multiple interactions were observed among the U11/U12-snRNP proteins. Taken together, these results reveal that U11/U12-65K is an indispensible component of the minor spliceosome complex by binding to both U11/U12-59K and U12 snRNA, and that multiple interactions among the U11/U12-snRNP proteins are necessary for minor spliceosome assembly.

The Arabidopsis homolog of human minor spliceosomal protein U11-48K plays a crucial role in U12 intron splicing and plant development

The minor U12 introns are removed from precursor mRNAs by the U12 intron-specific minor spliceosome. Among the seven ribonucleoproteins unique to the minor spliceosome, denoted as U11/U12-20K, U11/U12-25K, U11/U12-31K, U11/U12-65K, U11-35K, U11-48K, and U11-59K, the roles of only U11/U12-31K and U11/U12-65K have been demonstrated in U12 intron splicing and plant development. Here, the functional role of the Arabidopsis homolog of human U11-48K in U12 intron splicing and the development of Arabidopsis thaliana was examined using transgenic knockdown plants. The u11-48k mutants exhibited several defects in growth and development, such as severely arrested primary inflorescence stems, formation of serrated leaves, production of many rosette leaves after bolting, and delayed senescence. The splicing of most U12 introns analyzed was impaired in the u11-48k mutants. Comparative analysis of the splicing defects and phenotypes among the u11/u12-31k, u11-48k, and u11/12-65k mutants showed that the severity of abnormal development was closely correlated with the degree of impairment in U12 intron splicing. Taken together, these results provide compelling evidence that the Arabidopsis homolog of human U11-48K protein, as well as U11/U12-31K and U11/U12-65K proteins, is necessary for correct splicing of U12 introns and normal plant growth and development.

Experimental approaches to studying the nature and impact of splicing variation in zebrafish

From a fixed number of genes carried in all cells, organisms create considerable diversity in cellular phenotype through differential regulation of gene expression. One prevalent source of transcriptome diversity is alternative pre-mRNA splicing, which is manifested in many different forms. Zebrafish models of splicing dysfunction due to mutated spliceosome components provide opportunity to link biochemical analyses of spliceosome structure and function with whole organism phenotypic outcomes. Drawing from experience with two zebrafish mutants: cephalophŏnus (a prpf8 mutant, isolated for defects in granulopoiesis) and caliban (a rnpc3 mutant, isolated for defects in digestive organ development), we describe the use of glycerol gradient sedimentation and native gel electrophoresis to resolve components of aberrant splicing complexes. We also describe how RNAseq can be employed to examine relatively rare alternative splicing events including intron retention. Such experimental approaches in zebrafish can promote understanding of how splicing variation and dysfunction contribute to phenotypic diversity and disease pathogenesis.

Focusing the Spotlight on the Zebrafish Intestine to Illuminate Mechanisms of Colorectal Cancer

Colorectal cancer, encompassing colon and rectal cancer, arises from the epithelial lining of the large bowel. It is most prevalent in Westernised societies and is increasing in frequency as the world becomes more industrialised. Unfortunately, metastatic colorectal cancer is not cured by chemotherapy and the annual number of deaths caused by colorectal cancer, currently 700,000, is expected to rise. Our understanding of the contribution that genetic mutations make to colorectal cancer, although incomplete, is reasonably well advanced. However, it has only recently become widely appreciated that in addition to the ongoing accumulation of genetic mutations, chronic inflammation also plays a critical role in the initiation and progression of this disease. While a robust and tractable genetic model of colorectal cancer in zebrafish, suitable for pre-clinical studies, is not yet available, the identification of genes required for the rapid proliferation of zebrafish intestinal epithelial cells during development has highlighted a number of essential genes that could be targeted to disable colorectal cancer cells. Moreover, appreciation of the utility of zebrafish to study intestinal inflammation is on the rise. In particular, zebrafish provide unique opportunities to investigate the impact of genetic and environmental factors on the integrity of intestinal epithelial barrier function. With currently available tools, the interplay between epigenetic regulators, intestinal injury, microbiota composition and innate immune cell mobilisation can be analysed in exquisite detail. This provides excellent opportunities to define critical events that could potentially be targeted therapeutically. Further into the future, the use of zebrafish larvae as hosts for xenografts of human colorectal cancer tissue, while still in its infancy, holds great promise that zebrafish could one day provide a practical, preclinical personalized medicine platform for the rapid assessment of the metastatic potential and drug-sensitivity of patient-derived cancers.

Regulation of gene expression through inefficient splicing of U12-type introns

U12-type introns are a rare class of nuclear introns that are removed by a dedicated U12-dependent spliceosome and are thought to regulate the expression of their target genes owing through their slower splicing reaction. Recent genome-wide studies on the splicing of U12-type introns are now providing new insights on the biological significance of this parallel splicing machinery. The new studies cover multiple different organisms and experimental systems, including human patient cells with mutations in the components of the minor spliceosome, zebrafish with similar mutations and various experimentally manipulated human cells and Arabidopsis plants. Here, we will discuss the potential implications of these studies on the understanding of the mechanism and regulation of the minor spliceosome, as well as their medical implications.

Statistically based splicing detection reveals neural enrichment and tissue-specific induction of circular RNA during human fetal development

Background

The pervasive expression of circular RNA is a recently discovered feature of gene expression in highly diverged eukaryotes, but the functions of most circular RNAs are still unknown. Computational methods to discover and quantify circular RNA are essential. Moreover, discovering biological contexts where circular RNAs are regulated will shed light on potential functional roles they may play.

Results

We present a new algorithm that increases the sensitivity and specificity of circular RNA detection by discovering and quantifying circular and linear RNA splicing events at both annotated and un-annotated exon boundaries, including intergenic regions of the genome, with high statistical confidence. Unlike approaches that rely on read count and exon homology to determine confidence in prediction of circular RNA expression, our algorithm uses a statistical approach. Using our algorithm, we unveiled striking induction of general and tissue-specific circular RNAs, including in the heart and lung, during human fetal development. We discover regions of the human fetal brain, such as the frontal cortex, with marked enrichment for genes where circular RNA isoforms are dominant.

Conclusions

The vast majority of circular RNA production occurs at major spliceosome splice sites; however, we find the first examples of developmentally induced circular RNAs processed by the minor spliceosome, and an enriched propensity of minor spliceosome donors to splice into circular RNA at un-annotated, rather than annotated, exons. Together, these results suggest a potentially significant role for circular RNA in human development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0690-5) contains supplementary material, which is available to authorized users.

Aberrant splicing of U12-type introns is the hallmark of ZRSR2 mutant myelodysplastic syndrome

Somatic mutations in the spliceosome gene ZRSR2 — located on the X chromosome — are associated with myelodysplastic syndrome (MDS). ZRSR2 is involved in the recognition of 3΄ splice site during the early stages of spliceosome assembly; however, its precise role in RNA splicing has remained unclear. Here, we characterize ZRSR2 as an essential component of the minor spliceosome (U12-dependent) assembly. shRNA mediated knockdown of ZRSR2 leads to impaired splicing of the U12-type introns, and RNA-Sequencing of MDS bone marrow reveals that loss of ZRSR2 activity causes increased mis-splicing. These splicing defects involve retention of the U12-type introns while splicing of the U2-type introns remain mostly unaffected. ZRSR2 deficient cells also exhibit reduced proliferation potential and distinct alterations in myeloid and erythroid differentiation in vitro. These data identify a specific role for ZRSR2 in RNA splicing and highlight dysregulated splicing of U12-type introns as a characteristic feature of ZRSR2 mutations in MDS.

Minor splicing snRNAs are enriched in the developing mouse CNS and are crucial for survival of differentiating retinal neurons

In eukaryotes, gene expression requires splicing, which starts with the identification of exon-intron boundaries by the small, nuclear RNA (snRNAs) of the spliceosome, aided by associated proteins. In the mammalian genome, <1% of introns lack canonical exon-intron boundary sequences and cannot be spliced by the canonical splicing machinery. These introns are spliced by the minor spliceosome, consisting of unique snRNAs (U11, U12, U4atac, and U6atac). The importance of the minor spliceosome is underscored by the disease microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1), which is caused by mutation in U4atac. Thus, it is important to understand the expression and function of the minor spliceosome and its targets in mammalian development, for which we used the mouse as our model. Here, we report enrichment of the minor snRNAs in the developing head/central nervous system (CNS) between E9.5 and E12.5, along with enrichment of these snRNAs in differentiating retinal neurons. Moreover, dynamic expression kinetics of minor intron-containing genes (MIGs) was observed across retinal development. DAVID analysis of MIGs that were cotranscriptionally upregulated embryonically revealed enrichment for RNA metabolism and cell cycle regulation. In contrast, MIGs that were cotranscriptionally upregulated postnatally revealed enrichment for protein localization/transport, vesicle-mediated transport, and calcium transport. Finally, we used U12 morpholino to inactivate the minor spliceosome in the postnatal retina, which resulted in apoptosis of differentiating retinal neurons. Taken together, our data suggest that the minor spliceosome may have distinct functions in embryonic versus postnatal development. Importantly, we show that the minor spliceosome is crucial for the survival of terminally differentiating retinal neurons.

Global analysis of the nuclear processing of transcripts with unspliced U12-type introns by the exosome

U12-type introns are a rare class of introns in the genomes of diverse eukaryotes. In the human genome, they number over 700. A subset of these introns has been shown to be spliced at a slower rate compared to the major U2-type introns. This suggests a rate-limiting regulatory function for the minor spliceosome in the processing of transcripts containing U12-type introns. However, both the generality of slower splicing and the subsequent fate of partially processed pre-mRNAs remained unknown. Here, we present a global analysis of the nuclear retention of transcripts containing U12-type introns and provide evidence for the nuclear decay of such transcripts in human cells. Using SOLiD RNA sequencing technology, we find that, in normal cells, U12-type introns are on average 2-fold more retained than the surrounding U2-type introns. Furthermore, we find that knockdown of RRP41 and DIS3 subunits of the exosome stabilizes an overlapping set of U12-type introns. RRP41 knockdown leads to slower decay kinetics of U12-type introns and globally upregulates the retention of U12-type, but not U2-type, introns. Our results indicate that U12-type introns are spliced less efficiently and are targeted by the exosome. These characteristics support their role in the regulation of cellular mRNA levels.