Regulation of the MID1 protein function is fine-tuned by a complex pattern of alternative splicing

Perspective in Alternative Splicing Coupled to Nonsense-Mediated mRNA Decay

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is a cellular post-transcriptional process that generates protein isoform diversity. Nonsense-mediated RNA decay (NMD) is an mRNA surveillance pathway that recognizes and selectively degrades transcripts containing premature translation-termination codons (PTCs), thereby preventing the production of truncated proteins. Nevertheless, NMD also fine-tunes the gene expression of physiological mRNAs encoding full-length proteins. Interestingly, around one third of all AS events results in PTC-containing transcripts that undergo NMD. Numerous studies have reported a coordinated action between AS and NMD, in order to regulate the expression of several genes, especially those coding for RNA-binding proteins (RBPs). This coupling of AS to NMD (AS-NMD) is considered a gene expression tool that controls the ratio of productive to unproductive mRNA isoforms, ultimately degrading PTC-containing non-functional mRNAs. In this review, we focus on the mechanisms underlying AS-NMD, and how this regulatory process is able to control the homeostatic expression of numerous RBPs, including splicing factors, through auto- and cross-regulatory feedback loops. Furthermore, we discuss the importance of AS-NMD in the regulation of biological processes, such as cell differentiation. Finally, we analyze interesting recent data on the relevance of AS-NMD to human health, covering its potential roles in cancer and other disorders.

The MID1 gene product in physiology and disease

MID1 is an E3 ubiquitin ligase of the Tripartite Motif (TRIM) subfamily of RING-containing proteins, hence also known as TRIM18. MID1 is a microtubule-binding protein found in complex with the catalytic subunit of PP2A (PP2Ac) and its regulatory subunit alpha 4 (α4). To date, several substrates and interactors of MID1 have been described, providing evidence for the involvement of MID1 in a plethora of essential biological processes, especially during embryonic development. Mutations in the MID1 gene are responsible of the X-linked form of Opitz syndrome (XLOS), a multiple congenital disease characterised by defects in the development of midline structures during embryogenesis. Here, we review MID1-related physiological mechanisms as well as the pathological implication of the MID1 gene in XLOS and in other clinical conditions.

RNA Degradation in Neurodegenerative Disease

Ribonucleic acid (RNA) homeostasis is dynamically modulated in response to changing physiological conditions. Tight regulation of RNA abundance through both transcription and degradation determines the amount, timing, and location of protein translation. This balance is of particular importance in neurons, which are among the most metabolically active and morphologically complex cells in the body. As a result, any disruptions in RNA degradation can have dramatic consequences for neuronal health. In this chapter, we will first discuss mechanisms of RNA stabilization and decay. We will then explore how the disruption of these pathways can lead to neurodegenerative disease.

Identification of a classic nuclear localization signal at the N terminus that regulates the subcellular localization of Rbfox2 isoforms during differentiation of NMuMG and P19 cells

Nuclear localization of the alternative splicing factor Rbfox2 is achieved by a C-terminal nuclear localization signal (NLS) which can be excluded from some Rbfox2 isoforms by alternative splicing. While this predicts nuclear and cytoplasmic localization, Rbfox2 is exclusively nuclear in some cell types. Here, we identify a second NLS in the N terminus of Rbfox2 isoform 1A that is not included in Rbfox2 isoform 1F. Rbfox2 1A isoforms lacking the C-terminal NLS are nuclear, whereas equivalent 1F isoforms are cytoplasmic. A shift in Rbfox2 expression toward cytoplasmic 1F isoforms occurs during epithelial to mesenchymal transition (EMT) and could be important in regulating the activity and function of Rbfox2.

Differential Transcriptome Networks between IDO1-Knockout and Wild-Type Mice in Brain Microglia and Macrophages

Microglia in the brain and macrophages in peripheral organs are cell types responsible for immune response to challenges. Indoleamine 2,3-dioxygenase 1 (IDO1) is an immunomodulatory enzyme of the tryptophan pathway that is expressed in the brain. The higher activity of IDO1 in response to immune challenge has been implicated in behavioral disorders. The impact of IDO1 depletion on the microglia transcriptome has not been studied. An investigation of the transcript networks in the brain microglia from IDO1-knockout (IDO1-KO) mice was undertaken, relative to peripheral macrophages and to wild-type (WT) mice under unchallenged conditions. Over 105 transcript isoforms were differentially expressed between WT and IDO1-KO within cell type. Within microglia, Saa3 and Irg1 were over-expressed in IDO1-KO relative to WT. Within macrophages, Csf3 and Sele were over-expressed in IDO1-KO relative to WT. Among the genes differentially expressed between strains, enriched biological processes included ion homeostasis and ensheathment of neurons within microglia, and cytokine and chemokine expression within macrophages. Over 11,110 transcript isoforms were differentially expressed between microglia and macrophages and of these, over 10,800 transcripts overlapped between strains. Enriched biological processes among the genes over- and under-expressed in microglia relative to macrophages included cell adhesion and apoptosis, respectively. Detected only in microglia or macrophages were 421 and 43 transcript isoforms, respectively. Alternative splicing between cell types based on differential transcript isoform abundance was detected in 210 genes including Phf11d, H2afy, and Abr. Across strains, networks depicted a predominance of genes under-expressed in microglia relative to macrophages that may be a precursor for the different response of both cell types to challenges. The detected transcriptome differences enhance the understanding of the role of IDO1 in the microglia transcriptome under unchallenged conditions.

The RNA-binding protein Rbfox2: an essential regulator of EMT-driven alternative splicing and a mediator of cellular invasion

The epithelial-mesenchymal transition (EMT), a prerequisite for cancer progression and metastasis formation, is regulated not only at the transcriptional but also at the post-transcriptional level, including at the level of alternative pre-mRNA splicing. Several recent studies have highlighted the involvement of splicing factors, including epithelial splicing regulatory proteins (Esrps) and RNA-binding Fox protein 2 (Rbfox2), in this process. Esrps regulate epithelial-specific splicing, and their expression is downregulated during EMT. By contrast, the role of Rbfox2 is controversial because Rbfox2 regulates epithelial as well as mesenchymal splicing events. Here, we have used several established cell culture models to investigate the functions of Rbfox2 during EMT. We demonstrate that induction of an EMT upregulates the expression of Rbfox2, which correlates with an increase in Rbfox2-regulated splicing events in the cortactin (Cttn), Pard3 and dynamin 2 (Dnm2) transcripts. At the same time, however, the epithelial-specific ability to splice the Enah, Slk and Tsc2 transcripts is either reduced or lost completely by Rbfox2, which might be due, in part, to downregulation of the expression of the Esrps cooperative factors. Depletion of Rbfox2 during EMT did not prevent the activation of transforming growth factor-β signaling, the upregulation of mesenchymal markers or changes in cell morphology toward a mesenchymal phenotype. In addition, this depletion did not influence cell migration. However, depletion of Rbfox2 in cells that have completed an EMT significantly reduced their invasive potential. Taken together, our results suggest that during an EMT, Rbfox2-regulated splicing shifts from epithelial-to mesenchymal-specific events, leading to a higher degree of tissue invasiveness.

A pronounced evolutionary shift of the pseudoautosomal region boundary in house mice

The pseudoautosomal region (PAR) is essential for the accurate pairing and segregation of the X and Y chromosomes during meiosis. Despite its functional significance, the PAR shows substantial evolutionary divergence in structure and sequence between mammalian species. An instructive example of PAR evolution is the house mouse Mus musculus domesticus (represented by the C57BL/6J strain), which has the smallest PAR among those that have been mapped. In C57BL/6J, the PAR boundary is located just ~700 kb from the distal end of the X chromosome, whereas the boundary is found at a more proximal position in Mus spretus, a species that diverged from house mice 2-4 million years ago. In this study we used a combination of genetic and physical mapping to document a pronounced shift in the PAR boundary in a second house mouse subspecies, Mus musculus castaneus (represented by the CAST/EiJ strain), ~430 kb proximal of the M. m. domesticus boundary. We demonstrate molecular evolutionary consequences of this shift, including a marked lineage-specific increase in sequence divergence within Mid1, a gene that resides entirely within the M. m. castaneus PAR but straddles the boundary in other subspecies. Our results extend observations of structural divergence in the PAR to closely related subspecies, pointing to major evolutionary changes in this functionally important genomic region over a short time period.

Reduced Mid1 Expression and Delayed Neuromotor Development in daDREAM Transgenic Mice

Downstream regulatory element antagonist modulator (DREAM) is a Ca²⁺-binding protein that binds DNA and represses transcription in a Ca²⁺-dependent manner. Previous work has shown a role for DREAM in cerebellar function regulating the expression of the sodium/calcium exchanger 3 (NCX3) in cerebellar granular neurons to control Ca²⁺ homeostasis and survival of these neurons. To achieve a global view of the genes regulated by DREAM in the cerebellum, we performed a genome-wide analysis in transgenic cerebellum expressing a Ca²⁺-insensitive/CREB-independent dominant active mutant DREAM (daDREAM). Here we show that DREAM regulates the expression of the midline 1 (Mid1) gene early after birth. As a consequence, daDREAM mice exhibit a significant shortening of the rostro-caudal axis of the cerebellum and a delay in neuromotor development early after birth. Our results indicate a role for DREAM in cerebellar function.

A five-gene signature as a potential predictor of metastasis and survival in colorectal cancer

To understand the molecular mechanisms of metastasis and prognosis of colorectal cancer (CRC), we isolated single cell-derived progenies (SCPs) from SW480 cells in vitro and compared their metastatic potential in an orthotopic CRC tumour model in vivo. Two groups of SCPs with the capability of high and low metastasis, respectively, were obtained. By analysing the gene expression profiles of high (SCP51), low (SCP58) metastatic SCPs, and their parental cell line (SW480/EGFP), we demonstrated that 143 genes were differentially expressed either between SCP51 and SCP58 or between SCP58 and SW480/EGFP. Gene-annotation enrichment analysis of DAVID revealed 80 genes in the top ten clusters of the analysis (gene enrichment score > 1). Of the 80-gene set, 32 genes are potentially involved in metastasis, as revealed by Geneclip. Five putative metastatic genes (LYN, SDCBP, MAP4K4, DKK1, and MID1) were selected for further validations. Immunohistochemical analysis in a cohort of 181 CRC clinical samples showed that the individual expression of LYN, MAP4K4, and MID1, as well as the five-gene signature, was closely correlated with lymph node metastasis in CRC patients. More importantly, the individual expression of LYN, MAP4K4, SDCBP, and MID1, as well as the five-gene signature, was significantly correlated with overall survival in CRC patients. Thus, our five-gene signature may be able to predict metastasis and survival of CRC in the clinic, and opens new perspectives on the biology of CRC.

Point mutations in GLI3 lead to misregulation of its subcellular localization

Background

Mutations in the transcription factor GLI3, a downstream target of Sonic Hedgehog (SHH) signaling, are responsible for the development of malformation syndromes such as Greig-cephalopolysyndactyly-syndrome (GCPS), or Pallister-Hall-syndrome (PHS). Mutations that lead to loss of function of the protein and to haploinsufficiency cause GCPS, while truncating mutations that result in constitutive repressor function of GLI3 lead to PHS. As an exception, some point mutations in the C-terminal part of GLI3 observed in GCPS patients have so far not been linked to loss of function. We have shown recently that protein phosphatase 2A (PP2A) regulates the nuclear localization and transcriptional activity a of GLI3 function.

Principal Findings

We have shown recently that protein phosphatase 2A (PP2A) and the ubiquitin ligase MID1 regulate the nuclear localization and transcriptional activity of GLI3. Here we show mapping of the functional interaction between the MID1-α4-PP2A complex and GLI3 to a region between amino acid 568-1100 of GLI3. Furthermore we demonstrate that GCPS-associated point mutations, that are located in that region, lead to misregulation of the nuclear GLI3-localization and transcriptional activity. GLI3 phosphorylation itself however appears independent of its localization and remains untouched by either of the point mutations and by PP2A-activity, which suggests involvement of an as yet unknown GLI3 interaction partner, the phosphorylation status of which is regulated by PP2A activity, in the control of GLI3 subcellular localization and activity.

Conclusions

The present findings provide an explanation for the pathogenesis of GCPS in patients carrying C-terminal point mutations, and close the gap in our understanding of how GLI3-genotypes give rise to particular phenotypes. Furthermore, they provide a molecular explanation for the phenotypic overlap between Opitz syndrome patients with dysregulated PP2A-activity and syndromes caused by GLI3-mutations.

Stochastic noise in splicing machinery

The number of known alternative human isoforms has been increasing steadily with the amount of available transcription data. To date, over 100 000 isoforms have been detected in EST libraries, and at least 75% of human genes have at least one alternative isoform. In this paper, we propose that most alternative splicing events are the result of noise in the splicing process. We show that the number of isoforms and their abundance can be predicted by a simple stochastic noise model that takes into account two factors: the number of introns in a gene and the expression level of a gene. The results strongly support the hypothesis that most alternative splicing is a consequence of stochastic noise in the splicing machinery, and has no functional significance. The results are also consistent with error rates tuned to ensure that an adequate level of functional product is produced and to reduce the toxic effect of accumulation of misfolding proteins. Based on simulation of sampling of virtual cDNA libraries, we estimate that error rates range from 1 to 10% depending on the number of introns and the expression level of a gene.

Structural implication of splicing stochastics

Even though nearly every human gene has at least one alternative splice form, very little is so far known about the structure and function of resulting protein products. It is becoming increasingly clear that a significant fraction of all isoforms are products of noisy selection of splice sites and thus contribute little to actual functional diversity, and may potentially be deleterious. In this study, we examine the impact of alternative splicing on protein sequence and structure in three datasets: alternative splicing events conserved across multiple species, alternative splicing events in genes that are strongly linked to disease and all observed alternative splicing events. We find that the vast majority of all alternative isoforms result in unstable protein conformations. In contrast to that, the small subset of isoforms conserved across species tends to maintain protein structural integrity to a greater extent. Alternative splicing in disease-associated genes produces unstable structures just as frequently as all other genes, indicating that selection to reduce the effects of alternative splicing on this set is not especially pronounced. Overall, the properties of alternative spliced proteins are consistent with the outcome of noisy selection of splice sites by splicing machinery.

MID1 mutations in patients with X-linked Opitz G/BBB syndrome

Mutations in the MID1 gene are responsible for the X-linked form of Opitz G/BBB syndrome (OS), a disorder that affects the development of midline structures. OS is characterized by hypertelorism, hypospadias, laryngo-tracheo-esophageal (LTE) abnormalities, and additional midline defects. Cardiac, anal, and neurological defects are also present. The expressivity of OS is highly variable, even within the same family. We reviewed all the MID1 mutations reported so far, in both familial and sporadic cases. The mutations are scattered along the entire length of the gene and consist of missense and nonsense mutations, insertions and deletions, either in-frame or causing frameshifts, and deletions of either single exons or the entire MID1 coding region. The variety of described mutations and the lack of a strict genotype-phenotype correlation confirm the previous suggestion of the OS phenotype being caused by a loss-of-function mechanism. However, although a specific mutation cannot entirely account for the observed phenotype, we observed preferential association between some types of mutation and specific clinical manifestations, e.g., brain anatomical defects and truncating mutations. This may suggest that the pathogenetic mechanism underlying the OS phenotype is more complex and may vary among the affected organs.

Alternative polyadenylation signals and promoters act in concert to control tissue-specific expression of the Opitz Syndrome gene MID1

BACKGROUND

Mutations in the X-linked MID1 gene are responsible for Opitz G/BBB syndrome, a malformation disorder of developing midline structures. Previous Northern blot analyses revealed the existence of at least three MID1 transcripts of differing lengths.

RESULTS

Here we show that alternative polyadenylation generates the size differences observed in the Northern blot analyses. Analysis of EST data together with additional Northern blot analyses proved tissue-specific usage of the alternative polyadenylation sites. Bioinformatic characterization of the different 3'UTRs of MID1 revealed numerous RNA-protein interaction motifs, several of which turned out to be conserved between different species. Furthermore, our data suggest that mRNA termination at different polyadenylation sites is predetermined by the choice of alternative 5'UTRs and promoters of the MID1 gene, a mechanism that efficiently allows synergistic function of 5' and 3'UTRs.

CONCLUSION

MID1 expression is tightly regulated through concerted action of alternative promoters and alternative polyadenylation signals both during embryonic development and in the adult.

The coupling of alternative splicing and nonsense-mediated mRNA decay

Most human genes exhibit alternative splicing, but not all alternatively spliced transcripts produce functional proteins. Computational and experimental results indicate that a substantial fraction of alternative splicing events in humans result in mRNA isoforms that harbor a premature termination codon (PTC). These transcripts are predicted to be degraded by the nonsense-mediated mRNA decay (NMD) pathway. One explanation for the abundance of PTC-containing isoforms is that they represent splicing errors that are identified and degraded by the NMD pathway. Another potential explanation for this startling observation is that cells may link alternative splicing and NMD to regulate the abundance of mRNA transcripts. This mechanism, which we call "Regulated Unproductive Splicing and Translation" (RUST), has been experimentally shown to regulate expression of a wide variety of genes in many organisms from yeast to human. It is frequently employed for autoregulation of proteins that affect the splicing process itself. Thus, alternative splicing and NMD act together to play an important role in regulating gene expression.

A structure-function study of MID1 mutations associated with a mild Opitz phenotype

The X-linked form of Opitz syndrome (OS) affects midline structures and produces a characteristic, but heterogeneous, phenotype that may include severe mental retardation, hypertelorism, broad nasal bridge, widow's peak, cleft lip/cleft palate, congenital heart disease, laryngotracheal defects, and hypospadias. The MID1 gene was implicated in OS by linkage to Xp22. It encodes a 667 amino acid protein that contains a RING finger motif, two B-box zinc fingers, a coiled-coil, a fibronectin type III (FNIII) domain, and a B30.2 domain. Several mutations in MID1 are associated with severe OS. Here, we describe an intelligent male with a milder phenotype characterized by hypertelorism, broad nasal bridge, widow's peak, mild hypospadias, pectus excavatum, and a surgically corrected tracheo-esophageal fistula. He has an above average intelligence and no cleft lip/palate or heart disease. We identified a novel mutation in MID1 (P441L) which is in exon 8 and functionally associated with the FNIII domain. While OS phenotypes have been attributed to mutations in the C-terminal part of MID1, little is currently known about the structure-function relationships of MID1 mutations, and how they affect phenotype. We find from a literature review that missense mutations within the FNIII domain of MID1 are associated with a milder presentation of OS than missense mutations elsewhere in MID1. All truncating mutations (frameshift, insertions/deletions) lead to severe OS. We used homology analysis of the MID1 FNIII domain to investigate structure-function changes caused by our missense mutation. This and other missense mutations probably cause disruption of protein-protein interactions, either within MID1 or between MID1 and other proteins. We correlate these protein structure-function findings to the absence of CNS or palatal changes and conclude that the FNIII domain of the MID1 protein may be involved in midline differentiation after neural tube and palatal structures are completed.

The importance of being divisible by three in alternative splicing

Alternative splicing events that are conserved in orthologous genes in different species are commonly viewed as reliable evidence of authentic, functionally significant alternative splicing events. Several recent bioinformatic analyses have shown that conserved alternative exons possess several features that distinguish them from alternative exons that are species-specific. One of the most striking differences between conserved and species-specific alternative exons is the high percentage of exons that preserve the reading frame (exons whose length is an exact multiple of 3, termed symmetrical exons) among the conserved alternative exons. Here, we examined conserved alternative exons and found several features that differentiate between symmetrical and non-symmetrical alternative exons. We show that symmetrical alternative exons have a strong tendency not to disrupt protein domain structures, whereas the tendency of non-symmetrical alternative exons to overlap with different fractions of protein domains is similar to that of constitutive exons. Additionally, skipping isoforms of non-symmetrical alternative exons are strongly underrepresented, compared with their including isoforms, suggesting that skipping of a large fraction of non-symmetrical alternative exons produces transcripts that are degraded by the nonsense-mediated mRNA decay mechanism. Non-symmetrical alternative exons also show a tendency to reside in the 5' half of the CDS. These findings suggest that alternative splicing of symmetrical and non-symmetrical exons is governed by different selective pressures and serves different purposes.

SVC: structured visualization of evolutionary sequence conservation

We have developed a web application for the detailed analysis and visualization of evolutionary sequence conservation in complex vertebrate genes. Given a pair of orthologous genes, the protein-coding sequences are aligned. When these sequences are mapped back onto their encoding exons in the genomes, a scaffold of the conserved gene structure naturally emerges. Sequence similarity between exons and introns is analysed and embedded into the gene structure scaffold. The visualization on the SVC server provides detailed information about evolutionarily conserved features of these genes. It further allows concise representation of complex splice patterns in the context of evolutionary conservation. A particular application of our tool arises from the fact that around mRNA editing sites both exonic and intronic sequences are highly conserved. This aids in delineation of these sites. SVC is available at .

The Vertebrate Genome Annotation (Vega) database

The Vertebrate Genome Annotation (Vega) database (http://vega.sanger.ac.uk) has been designed to be a community resource for browsing manual annotation of finished sequences from a variety of vertebrate genomes. Its core database is based on an Ensembl-style schema, extended to incorporate curation-specific metadata. In collaboration with the genome sequencing centres, Vega attempts to present consistent high-quality annotation of the published human chromosome sequences. In addition, it is also possible to view various finished regions from other vertebrates, including mouse and zebrafish. Vega displays only manually annotated gene structures built using transcriptional evidence, which can be examined in the browser. Attempts have been made to standardize the annotation procedure across each vertebrate genome, which should aid comparative analysis of orthologues across the different finished regions.