Intrinsic differences between authentic and cryptic 5' splice sites

Progerin mRNA expression in non-HGPS patients is correlated with widespread shifts in transcript isoforms

Hutchinson-Gilford Progeria Syndrome (HGPS) is a premature aging disease caused primarily by a C1824T mutation in LMNA. This mutation activates a cryptic splice donor site, producing a lamin variant called progerin. Interestingly, progerin has also been detected in cells and tissues of non-HGPS patients. Here, we investigated progerin expression using publicly available RNA-seq data from non-HGPS patients in the GTEx project. We found that progerin expression is present across all tissue types in non-HGPS patients and correlated with telomere shortening in the skin. Transcriptome-wide correlation analyses suggest that the level of progerin expression is correlated with switches in gene isoform expression patterns. Differential expression analyses show that progerin expression is correlated with significant changes in genes involved in splicing regulation and mitochondrial function. Interestingly, 5' splice sites whose use is correlated with progerin expression have significantly altered frequencies of consensus trinucleotides within the core 5' splice site. Furthermore, introns whose alternative splicing correlates with progerin have reduced GC content. Our study suggests that progerin expression in non-HGPS patients is part of a global shift in splicing patterns.

U-rich elements drive pervasive cryptic splicing in 3' UTR massively parallel reporter assays

Non-coding RNA sequences play essential roles in orchestrating gene expression. However, the sequence codes and mechanisms underpinning post-transcriptional regulation remain incompletely understood. Here, we revisit the finding from a prior massively parallel reporter assay (MPRA) that AU-rich (U-rich) elements in 3' untranslated regions (3' UTRs) can drive upregulation or downregulation of mRNA expression depending on 3' UTR context. We unexpectedly discover that this variable regulation arises from widespread cryptic splicing, predominately from an unannotated splice donor in the coding sequence of GFP to diverse acceptor sites in reporter 3' UTRs. Splicing is activated by U-rich sequences, which function as potent position-dependent regulators of 5' and 3' splice site choice and overall splicing efficiency. Splicing has diverse impacts on reporter expression, causing both increases and decreases in reporter expression via multiple mechanisms. We further provide evidence that cryptic splicing impacts between 10 to 50% of measurements made by other published 3' UTR MPRAs. Overall, our work emphasizes U-rich sequences as principal drivers of splicing and provides strategies to minimize cryptic splicing artifacts in reporter assays.

Pleiotropic brain function of whirlin identified by a novel mutation

Despite some evidence indicating diverse roles of whirlin in neurons, the functional corollary of whirlin gene function and behavior has not been investigated or broadly characterized. A single nucleotide variant was identified from our recessive ENU-mutagenesis screen at a donor-splice site in whirlin, a protein critical for proper sensorineural hearing function. The mutation (head-bob, hb) led to partial intron-retention causing a frameshift and introducing a premature termination codon. Mutant mice had a head-bobbing phenotype and significant hyperactivity across several phenotyping tests. Lack of complementation of head-bob with whirler mutant mice confirmed the head-bob mutation as functionally distinct with compound mutants having a mild-moderate hearing defect. Utilizing transgenics, we demonstrate rescue of the hyperactive phenotype and combined with the expression profiling data conclude whirlin plays an essential role in activity-related behaviors. These results highlight a pleiotropic role of whirlin within the brain and implicate alternative, central mediated pathways in its function.

Tumor Antigens beyond the Human Exome

With the advent of immunotherapeutics, a new era in the combat against cancer has begun. Particularly promising are neo-epitope-targeted therapies as the expression of neo-antigens is tumor-specific. In turn, this allows the selective targeting and killing of cancer cells whilst healthy cells remain largely unaffected. So far, many advances have been made in the development of treatment options which are tailored to the individual neo-epitope repertoire. The next big step is the achievement of efficacious "off-the-shelf" immunotherapies. For this, shared neo-epitopes propose an optimal target. Given the tremendous potential, a thorough understanding of the underlying mechanisms which lead to the formation of neo-antigens is of fundamental importance. Here, we review the various processes which result in the formation of neo-epitopes. Broadly, the origin of neo-epitopes can be categorized into three groups: canonical, noncanonical, and viral neo-epitopes. For the canonical neo-antigens that arise in direct consequence of somatic mutations, we summarize past and recent findings. Beyond that, our main focus is put on the discussion of noncanonical and viral neo-epitopes as we believe that targeting those provides an encouraging perspective to shape the future of cancer immunotherapeutics.

Characterization of NEB pathogenic variants in patients reveals novel nemaline myopathy disease mechanisms and omecamtiv mecarbil force effects

Nebulin, a critical protein of the skeletal muscle thin filament, plays important roles in physiological processes such as regulating thin filament length (TFL), cross-bridge cycling, and myofibril alignment. Pathogenic variants in the nebulin gene (NEB) cause NEB-based nemaline myopathy (NEM2), a genetically heterogeneous disorder characterized by hypotonia and muscle weakness, currently lacking curative therapies. In this study, we examined a cohort of ten NEM2 patients, each with unique pathogenic variants, aiming to understand their impact on mRNA, protein, and functional levels. Results show that pathogenic truncation variants affect NEB mRNA stability and lead to nonsense-mediated decay of the mutated transcript. Moreover, a high incidence of cryptic splice site activation was found in patients with pathogenic splicing variants that are expected to disrupt the actin-binding sites of nebulin. Determination of protein levels revealed patients with either relatively normal or markedly reduced nebulin. We observed a positive relation between the reduction in nebulin and a reduction in TFL, or reduction in tension (both maximal and submaximal tension). Interestingly, our study revealed a pathogenic duplication variant in nebulin that resulted in a four-copy gain in the triplicate region of NEB and a much larger nebulin protein and longer TFL. Additionally, we investigated the effect of Omecamtiv mecarbil (OM), a small-molecule activator of cardiac myosin, on force production of type 1 muscle fibers of NEM2 patients. OM treatment substantially increased submaximal tension across all NEM2 patients ranging from 87 to 318%, with the largest effects in patients with the lowest level of nebulin. In summary, this study indicates that post-transcriptional or post-translational mechanisms regulate nebulin expression. Moreover, we propose that the pathomechanism of NEM2 involves not only shortened but also elongated thin filaments, along with the disruption of actin-binding sites resulting from pathogenic splicing variants. Significantly, our findings highlight the potential of OM treatment to improve skeletal muscle function in NEM2 patients, especially those with large reductions in nebulin levels.

Possibilities and limitations of antisense oligonucleotide therapies for the treatment of monogenic disorders

Antisense oligonucleotides (ASOs) are incredibly versatile molecules that can be designed to specifically target and modify RNA transcripts to slow down or halt rare genetic disease progression. They offer the potential to target groups of patients or can be tailored for individual cases. Nonetheless, not all genetic variants and disorders are amenable to ASO-based treatments, and hence, it is important to consider several factors before embarking on the drug development journey. Here, we discuss which genetic disorders have the potential to benefit from a specific type of ASO approach, based on the pathophysiology of the disease and pathogenic variant type, as well as those disorders that might not be suitable for ASO therapies. We further explore additional aspects, such as the target tissues, intervention time points, and potential clinical benefits, which need to be considered before developing a compound. Overall, we provide an overview of the current potentials and limitations of ASO-based therapeutics for the treatment of monogenic disorders.

Characterization of NEB mutations in patients reveals novel nemaline myopathy disease mechanisms and omecamtiv mecarbil force effects

Nebulin, a critical protein of the skeletal muscle thin filament, plays important roles in physiological processes such as regulating thin filament length (TFL), cross-bridge cycling, and myofibril alignment. Mutations in the nebulin gene ( NEB ) cause NEB-based nemaline myopathy (NEM2), a genetically heterogeneous disorder characterized by hypotonia and muscle weakness, currently lacking therapies targeting the underlying pathological mechanisms. In this study, we examined a cohort of ten NEM2 patients, each with unique mutations, aiming to understand their impact on mRNA, protein, and functional levels. Results show that truncation mutations affect NEB mRNA stability and lead to nonsense-mediated decay of the mutated transcript. Moreover, a high incidence of cryptic splice site activation was found in patients with splicing mutations which is expected to disrupt the actin-binding sites of nebulin. Determination of protein levels revealed patients with relatively normal nebulin levels and others with markedly reduced nebulin. We observed a positive relation between the reduction in nebulin and a reduction in TFL, and a positive relation between the reduction in nebulin level and the reduction in tension (both maximal and submaximal tension). Interestingly, our study revealed a duplication mutation in nebulin that resulted in a larger nebulin protein and longer TFL. Additionally, we investigated the effect of Omecamtiv mecarbil (OM), a small-molecule activator of cardiac myosin, on force production of type I muscle fibers of NEM2 patients. OM treatment substantially increased submaximal tension across all NEM2 patients ranging from 87-318%, with the largest effects in patients with the lowest level of nebulin. In summary, this study indicates that post-transcriptional or post-translational mechanisms regulate nebulin expression. Moreover, we propose that the pathomechanism of NEM2 involves not only shortened but also elongated thin filaments, along with the disruption of actin-binding sites resulting from splicing mutations. Significantly, our findings highlight the potential of OM treatment to improve skeletal muscle function in NEM2 patients, especially those with large reductions in nebulin levels.

Hypoparathyroidism, deafness and renal dysplasia syndrome caused by a GATA3 splice site mutation leading to the activation of a cryptic splice site

The HDR syndrome is a rare autosomal dominant disorder characterised by Hypoparathyroidism, Deafness, and Renal dysplasia, and is caused by inactivating heterozygous germline mutations in the GATA3 gene. We report an 11-year-old girl with HDR syndrome caused by a heterozygous mutation located at the splice acceptor site of exon 5 of the GATA3 gene (NM_001002295.2: c.925-1G>T). Functional studies using a minigene assay showed that this splice site mutation abolished the normal splicing of the GATA3 pre-mRNA and led to the use of a cryptic splice acceptor site, resulting in the loss of the first seven nucleotides (TCTGCAG) of exon 5 in the GATA3 mRNA. These findings increase the understanding of the mechanisms by which GATA3 splicing mutations can cause HDR syndrome.

Identification and interpretation of TET2 noncanonical splicing site intronic variants in myeloid neoplasm patients

Background: DNA hypermethylation and instability due to inactivation mutations in Ten-eleven translocation 2 (TET2) is a key biomarker of hematological malignancies. This study aims at characterizing two intronic noncanonical splice-site variants, c.3954+5_3954+8delGTTT and c.3954+5G>A. Methods: We used in silico prediction tools, reverse transcription (RT)-PCR, and Sanger sequencing on blood/bone marrow-derived RNA specimens to determine the aberrant splicing. Results: In silico prediction of both variants exhibited reduced splicing strength at the TET2 intron 7 splicing donor site. RT-PCR and Sanger sequencing identified a 62-bp deletion at the exon 7, producing a frameshift mutation, p.Cys1298*. Conclusion: This study provides functional evidence for two intronic TET2 variants that cause alternative splicing and frameshift mutation.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Splicing accuracy varies across human introns, tissues and age

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigated mis-splicing using RNA-sequencing data from ~14K control samples and 42 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that mis-splicing occurs at different rates across introns and tissues and that these splicing inaccuracies are primarily affected by the abundance of core components of the spliceosome assembly and its regulators. Using publicly available data on short-hairpin RNA-knockdowns of numerous spliceosomal components and related regulators, we found support for the importance of RNA-binding proteins in mis-splicing. We also demonstrated that age is positively correlated with mis-splicing, and it affects genes implicated in neurodegenerative diseases. This in-depth characterisation of mis-splicing can have important implications for our understanding of the role of splicing inaccuracies in human disease and the interpretation of long-read RNA-sequencing data.

Isoform changes of action potential regulators in the ventricles of arrhythmogenic phospholamban-R14del humanized mouse hearts

Arrhythmogenic cardiomyopathy (ACM) is characterized by life-threatening ventricular arrhythmias and sudden cardiac death and affects hundreds of thousands of patients worldwide. The deletion of Arginine 14 (p.R14del) in the phospholamban (PLN) gene has been implicated in the pathogenesis of ACM. PLN is a key regulator of sarcoplasmic reticulum (SR) Ca2+ cycling and cardiac contractility. Despite global gene and protein expression studies, the molecular mechanisms of PLN-R14del ACM pathogenesis remain unclear. Using a humanized PLN-R14del mouse model and human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs), we investigated the transcriptome-wide mRNA splicing changes associated with the R14del mutation. We identified >200 significant alternative splicing (AS) events and distinct AS profiles were observed in the right (RV) and left (LV) ventricles in PLN-R14del compared to WT mouse hearts. Enrichment analysis of the AS events showed that the most affected biological process was associated with "cardiac cell action potential", specifically in the RV. We found that splicing of 2 key genes, Trpm4 and Camk2d, which encode proteins regulating calcium homeostasis in the heart, were altered in PLN-R14del mouse hearts and human iPSC-CMs. Bioinformatical analysis pointed to the tissue-specific splicing factors Srrm4 and Nova1 as likely upstream regulators of the observed splicing changes in the PLN-R14del cardiomyocytes. Our findings suggest that aberrant splicing may affect Ca2+-homeostasis in the heart, contributing to the increased risk of arrythmogenesis in PLN-R14del ACM.

A genetic screen in Arabidopsis reveals the identical roles for RBP45d and PRP39a in 5' cryptic splice site selection

Cryptic splice sites in eukaryotic genome are generally dormant unless activated by mutation of authentic splice sites or related splicing factors. How cryptic splice sites are used remains unclear in plants. Here, we identified two cryptic splicing regulators, RBP45d and PRP39a that are homologs of yeast U1 auxiliary protein Nam8 and Prp39, respectively, via genetic screening for suppressors of the virescent sot5 mutant, which results from a point mutation at the 5' splice site (5' ss) of SOT5 intron 7. Loss-of-function mutations in RBP45d and PRP39a significantly increase the level of a cryptically spliced variant that encodes a mutated but functional sot5 protein, rescuing sot5 to the WT phenotype. We furtherly demonstrated that RBP45d and PRP39a interact with each other and also with the U1C, a core subunit of U1 snRNP. We found that RBP45d directly binds to the uridine (U)-rich RNA sequence downstream the 5' ss of SOT5 intron 7. However, other RBP45/47 members do not function redundantly with RBP45d, at least in regulation of cryptic splicing. Taken together, RBP45d promotes U1 snRNP to recognize the specific 5' ss via binding to intronic U-rich elements in plants.

Curating the gnomAD database: Report of novel variants in the thyroid peroxidase gene using in silico bioinformatics algorithms and a literature review

Thyroid peroxidase (TPO) is a membrane-bound glycoprotein located at the apical side of the thyroid follicular cells that catalyzes both iodination and coupling of iodotyrosine residues within the thyroglobulin molecule, leading to the synthesis of thyroid hormone. Variants in TPO cause congenital hypothyroidism (CH) by iodide organification defect and are commonly inherited in an autosomal recessive fashion. In the present work, we report a detailed population analysis and bioinformatic prediction of the TPO variants indexed in the Genome Aggregation Database (gnomAD) v2.1.1. The proportion of missense cysteine variants and nonsense, frameshift, and splice acceptor/donor variants were analyzed in each ethnic group (European (Non-Finnish), European (Finnish), African/African Americans, Latino/Admixed American, East Asian, South Asian, Ashkenazi Jewish, Other). The results showed a clear predominance of frameshift variants in the East Asian (82%) and European (Finnish) (75%) population, whereas the splice site variants predominate in African/African Americans (99.46%), Other (96%), Latino/Admixed American (94%), South Asian (86%), European (Non-Finnish) (56%) and Ashkenazi Jewish (56%) populations. The analysis of the distribution of the variants indexed in gnomAD v2.1.1 database revealed that most missense variants identified in the An peroxidase domain map in exon 8, followed by exons 11, 7 and 9, and finally in descending order by exons 10, 6, 12 and 5. In total, 183 novel TPO variants were described (13 missense cysteine's variants, 158 missense variants involving the An peroxidase domain and 12 splicing acceptor or donor sites variants) which were not reported in the literature and that would have deleterious effects on prediction programs. In the gnomAD v2.1.1 population, the estimated prevalence of heterozygous carriers of the potentially damaging variants was 1:77. In conclusion, we provide an updated and curated reference source of new TPO variants for application in clinical diagnosis and genetic counseling. Also, this work contributes to elucidating the molecular basis of CH associated with TPO defects.

Cooperative engagement and subsequent selective displacement of SR proteins define the pre-mRNA 3D structural scaffold for early spliceosome assembly

We recently reported that serine–arginine-rich (SR) protein-mediated pre-mRNA structural remodeling generates a pre-mRNA 3D structural scaffold that is stably recognized by the early spliceosomal components. However, the intermediate steps between the free pre-mRNA and the assembled early spliceosome are not yet characterized. By probing the early spliceosomal complexes in vitro and RNA-protein interactions in vivo, we show that the SR proteins bind the pre-mRNAs cooperatively generating a substrate that recruits U1 snRNP and U2AF65 in a splice signal-independent manner. Excess U1 snRNP selectively displaces some of the SR protein molecules from the pre-mRNA generating the substrate for splice signal-specific, sequential recognition by U1 snRNP, U2AF65 and U2AF35. Our work thus identifies a novel function of U1 snRNP in mammalian splicing substrate definition, explains the need for excess U1 snRNP compared to other U snRNPs in vivo, demonstrates how excess SR proteins could inhibit splicing, and provides a conceptual basis to examine if this mechanism of splicing substrate definition is employed by other splicing regulatory proteins.

Exon junction complex-associated multi-adapter RNPS1 nucleates splicing regulatory complexes to maintain transcriptome surveillance

The exon junction complex (EJC) is an RNA-binding multi-protein complex with critical functions in post-transcriptional gene regulation. It is deposited on the mRNA during splicing and regulates diverse processes including pre-mRNA splicing and nonsense-mediated mRNA decay (NMD) via various interacting proteins. The peripheral EJC-binding protein RNPS1 was reported to serve two insufficiently characterized functions: suppressing mis-splicing of cryptic splice sites and activating NMD in the cytoplasm. The analysis of transcriptome-wide effects of EJC and RNPS1 knockdowns in different human cell lines supports the conclusion that RNPS1 can moderately influence NMD activity, but is not a globally essential NMD factor. However, numerous aberrant splicing events strongly suggest that the main function of RNPS1 is splicing regulation. Rescue analyses revealed that the RRM and C-terminal domain of RNPS1 both contribute partially to regulate RNPS1-dependent splicing events. We defined the RNPS1 core interactome using complementary immunoprecipitations and proximity labeling, which identified interactions with splicing-regulatory factors that are dependent on the C-terminus or the RRM domain of RNPS1. Thus, RNPS1 emerges as a multifunctional splicing regulator that promotes correct and efficient splicing of different vulnerable splicing events via the formation of diverse splicing-promoting complexes.

Congenital nephrotic syndrome in a Hispanic Guatemalan newborn associated with a NPHS1 variant: A case report

Congenital nephrotic syndrome (CNS) is an autosomal recessive disorder usually detected in the first 3 months of life when the syndromes effects manifest, including edema and a failure to gain weight. A baby boy was admitted to the Neonatal Intensive Care Unit for prematurity (35 weeks) with unremarkable maternal prenatal laboratory tests. The patient had persistent systemic hypertension, hypoproteinemia, hypoalbuminemia and nephrotic range proteinuria. CNS was diagnosed, and genetic testing showed a homozygous variant, c.3024A>G (AGA>AGG) in exon 22 of the nephrin locus. Bioinformatics analysis suggested the genetic condition was likely a result of malfunctional DNA binding sites of transcription factors FOXL1 and FOXC1.

Novel compound heterozygous mutations in the CYP4F22 gene in a patient with autosomal recessive congenital ichthyosis

Autosomal recessive congenital ichthyosis (ARCI) is a rare form of keratinization disorder of the skin, which can be caused by mutations in 14 ARCI genes. We present a rare case of ARCI that carried a novel null mutation and a novel splice site mutation in the CYP4F22 gene.

Curating the gnomAD database: Report of novel variants in the thyrogobulin gene using in silico bioinformatics algorithms

Thyroglobulin (TG) is a large glycosylated protein of 2767 amino acids, secreted by the thyrocytes into the follicular lumen. It plays an essential role in the process of thyroid hormone synthesis. TG gene variants lead to permanent congenital hypothyroidism. In the present work, we report a detailed population and bioinformatic prediction analyses of the TG variants indexed in the Genome Aggregation Database (gnomAD). The results showed a clear predominance of nonsense variants in the European (Finnish), European (Non-Finnish) and Ashkenazi Jewish ethnic groups, whereas the splice site variants predominate in South Asian and African/African-American populations. In total, 282 novel TG variants were described (47 missense involving the wild-type cysteine residues, 177 missense located in the ChEL domain and 58 splice site variants) which were not reported in the literature and that would have deleterious effects in prediction programs. In the gnomAD population, the estimated prevalence of heterozygous carriers of the potentially damaging variants was 1:320. In conclusion, we provide an updated and curated reference source for the diagnosis of thyroid disease, mainly to congenital hypothyroidism due to TG deficiency. The identification and characterization of TG variants is undoubtedly a valuable approach to study the TG structure/function relations and an important tool for clinical diagnosis and genetic counseling.

Intron exon boundary junctions in human genome have in-built unique structural and energetic signals

Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.

A novel mutation in intron 11 donor splice site, responsible of a rare genotype in thyroglobulin gene by altering the pre-mRNA splincing process. Cell expression and bioinformatic analysis

Thyroglobulin (TG) is a homodimeric glycoprotein synthesized by the thyroid gland. To date, two hundred twenty-seven variations of the TG gene have been identified in humans. Thyroid dyshormonogenesis due to TG gene mutations have an estimated incidence of approximately 1 in 100,000 newborns. The clinical spectrum ranges from euthyroid to mild or severe hypothyroidism. The purpose of the present study was to identify and characterize new variants in the TG gene. We report an Argentine patient with congenital hypothyroidism, enlarged thyroid gland and low levels of serum TG. Sequencing of DNA, expression of chimeric minigenes as well as bioinformatics analysis were performed. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a c.3001+5G > A, whereas the paternal mutation consists of p.Arg296*. Minigen analysis of the variant c.3001+5A performed in HeLa, CV1 and Hek293T cell lines, showed a total lack of transcript expression. So, in order to validate that the loss of expression was caused by such variation, site-directed mutagenesis was performed on the mutated clone, which previously had a pSPL3 vector change, to give rise to a wild-type clone c.3001+5G, endorsing that the mutation c.3001+5G > A is the cause of the total lack of expression. In conclusion, we demonstrate that the c.3001+5G > A mutation causes a rare genotype, altering the splicing of the pre-mRNA. This work contributes to elucidating the molecular bases of TG defects associated with congenital hypothyroidism and expands our knowledge in relation to the pathologic roles of the position 5 in the donor splice site.

Transposon clusters as substrates for aberrant splice-site activation

Transposed elements (TEs) have dramatically shaped evolution of the exon-intron structure and significantly contributed to morbidity, but how recent TE invasions into older TEs cooperate in generating new coding sequences is poorly understood. Employing an updated repository of new exon-intron boundaries induced by pathogenic mutations, termed DBASS, here we identify novel TE clusters that facilitated exon selection. To explore the extent to which such TE exons maintain RNA secondary structure of their progenitors, we carried out structural studies with a composite exon that was derived from a long terminal repeat (LTR78) and AluJ and was activated by a C > T mutation optimizing the 5ʹ splice site. Using a combination of SHAPE, DMS and enzymatic probing, we show that the disease-causing mutation disrupted a conserved AluJ stem that evolved from helix 3.3 (or 5b) of 7SL RNA, liberating a primordial GC 5ʹ splice site from the paired conformation for interactions with the spliceosome. The mutation also reduced flexibility of conserved residues in adjacent exon-derived loops of the central Alu hairpin, revealing a cross-talk between traditional and auxilliary splicing motifs that evolved from opposite termini of 7SL RNA and were approximated by Watson-Crick base-pairing already in organisms without spliceosomal introns. We also identify existing Alu exons activated by the same RNA rearrangement. Collectively, these results provide valuable TE exon models for studying formation and kinetics of pre-mRNA building blocks required for splice-site selection and will be useful for fine-tuning auxilliary splicing motifs and exon and intron size constraints that govern aberrant splice-site activation.

Expression vector cassette engineering for recombinant therapeutic production in mammalian cell systems

Human tissue plasminogen activator was the first recombinant therapy protein that successfully produced in Chinese hamster ovary cells in 1986 and approved for clinical use. Since then, more and more therapeutic proteins are being manufactured in mammalian cells, and the technologies for recombinant protein production in this expression system have developed rapidly, with the optimization of both upstream and downstream processes. One of the most promising strategies is expression vector cassette optimization based on the expression vector cassette. In this review paper, these approaches and developments are summarized, and the future strategy on the utilizing of expression cassettes for the production of recombinant therapeutic proteins in mammalian cells is discussed.

Comprehensive genetic diagnosis of patients with Duchenne/Becker muscular dystrophy (DMD/BMD) and pathogenicity analysis of splice site variants in the DMD gene

Duchenne muscular dystrophy (DMD) and Becker muscular dystrophy (BMD) are caused by mutations in the DMD gene. The aim of this study is to identify pathogenic DMD variants in probands and reduce the risk of recurrence of the disease in affected families. Variations in 100 unrelated DMD/BMD patients were detected by multiplex ligation-dependent probe amplification (MLPA) and next-generation sequencing (NGS). Pathogenic variants in DMD were successfully identified in all cases, and 11 of them were novel. The most common mutations were intragenic deletions (69%), with two hotspots located in the 5' end (exons 2-19) and the central of the DMD gene (exons 45-55), while point mutations were observed in 22% patients. Further, c.1149+1G>A and c.1150-2A>G were confirmed by hybrid minigene splicing assay (HMSA). This two splice site mutations would lead to two aberrant DMD isoforms which give rise to severely truncated protein. Therefore, the clinical use of MLPA, NGS, and HMSA is an effective strategy to identify variants. Importantly, eight embryos were terminated pregnancies according to prenatal diagnosis and a healthy boy was successfully delivered by preimplantation genetic diagnosis (PGD). Early and accurate genetic diagnosis is essential for prenatal diagnosis/PGD to reduce the risk of recurrence of DMD in affected families.

Succession of splicing regulatory elements determines cryptic 5΄ss functionality

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5’ end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bβ-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3’ exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3’ exon end definition.

Identification of the Underlying Androgen Receptor Defect in the Dallas Reifenstein Family

CONTEXT

The Dallas Reifenstein family - first described in 1965 - includes 14 known members with partial androgen insensitivity syndrome (PAIS). However, the underlying molecular defect was never identified.

OBJECTIVE

To identify the underlying genetic defect for PAIS in the Dallas Reifenstein family.

DESIGN

DNA was purified from scrotal skin fibroblasts, and whole exome sequencing was then performed in four affected men in the family. Additional family members - both affected and unaffected - were subsequently recruited to confirm segregation of the candidate mutations with the PAIS phenotype.

PATIENTS

The affected men have PAIS with infertility associated with azoospermia, hypospadias, and gynecomastia.

RESULTS

All four men harbored an intronic variant NC_000023.10:g.66788676A>C between exon 1 and exon 2 of the androgen receptor (AR) canonical transcript NM_000044 (complementary DNA position NM_000044: c.1616+22072A>C) predicted to cause an alternatively spliced AR transcript. Reverse transcription (RT) polymerase chain (PCR) experiments detected the predicted PCR product of the alternatively spliced AR transcript, and the mutation segregated with the PAIS phenotype in this family. The transcript includes the insertion of 185 nucleotides with a premature stop codon at chrX:66863131-66863133, likely resulting in a reduction in AR protein expression due to nonsense-mediated decay.

CONCLUSIONS

An intronic AR mutation was identified in the Dallas Reifenstein family. The findings suggest that in cases of PAIS without identifiable AR mutations in coding regions, intronic AR mutations should be considered.

A novel NHS mutation causes Nance-Horan Syndrome in a Chinese family

Background

Nance-Horan Syndrome (NHS) (OMIM: 302350) is a rare X-linked developmental disorder characterized by bilateral congenital cataracts, with occasional dental anomalies, characteristic dysmorphic features, brachymetacarpia and mental retardation. Carrier females exhibit similar manifestations that are less severe than in affected males.

Methods

Here, we report a four-generation Chinese family with multiple affected individuals presenting Nance-Horan Syndrome. Whole-exome sequencing combined with RT-PCR and Sanger sequencing was used to search for a genetic cause underlying the disease phenotype.

Results

Whole-exome sequencing identified in all affected individuals of the family a novel donor splicing site mutation (NM_198270: c.1045 + 2T > A) in intron 4 of the gene NHS, which maps to chromosome Xp22.13. The identified mutation results in an RNA processing defect causing a 416-nucleotide addition to exon 4 of the mRNA transcript, likely producing a truncated NHS protein.

Conclusions

The donor splicing site mutation NM_198270: c.1045 + 2T > A of the NHS gene is the causative mutation in this Nance-Horan Syndrome family. This research broadens the spectrum of NHS gene mutations, contributing to our understanding of the molecular genetics of NHS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-016-0360-9) contains supplementary material, which is available to authorized users.

Molecular Outcome, Prediction, and Clinical Consequences of Splice Variants in COL1A1, Which Encodes the proα1(I) Chains of Type I Procollagen

Approximately 10%-20% of germline pathogenic variants alter mRNA splicing, with phenotypes often dependent on the stability of the mRNA produced by the mutant allele. To better understand the relationships between genotype, mRNA splicing, and phenotype, we examined clinical and molecular data from 243 probands with osteogenesis imperfecta (OI) representing 145 unique splicing variants within the type I procollagen gene, COL1A1. All individuals with IVSX-1G>A mutations had OI type I because the substitution shifted the splice acceptor site 1 nt downstream and destabilized the mRNA. OI phenotypes were not consistent for any other splice variant identified. We sequenced all cDNA species from cultured dermal fibroblasts from 40 individuals to identify splice outcome and compared those results to splice predictions from Human Splice Finder (HSF), Spliceport (SP), and Automatic Splice Site and Exon Definition Analyses (ASSEDA). Software-based splice predictions were correct in 42%, 55%, and 74% instances for HSF, SP, and ASSEDA, respectively. As molecular diagnostics move increasingly to DNA sequence analysis, the need to understand the effects of splice site variants will increase. These data demonstrate that caution must be exercised when using splice prediction software to predict splice outcome.

Compound heterozygous DUOX2 gene mutations (c.2335-1G>C/c.3264_3267delCAGC) associated with congenital hypothyroidism. Characterization of complex cryptic splice sites by minigene analysis

Iodide Organification defects (IOD) represent 10% of cases of congenital hypothyroidism (CH) being the main genes affected that of TPO (thyroid peroxidase) and DUOX2 (dual oxidasa 2). From a patient with clinical and biochemical criteria suggestive with CH associated with IOD, TPO and DUOX2 genes were analyzed by means of PCR-Single Strand Conformation Polymorphism analysis and sequencing. A novel heterozygous compound to the mutations c.2335-1G>C (paternal mutation, intron 17) and c.3264_3267delCAGC (maternal mutation, exon 24) was identified in the DUOX2 gene. Ex-vivo splicing assays and subsequent RT-PCR and sequencing analyses were performed on mRNA isolated from the HeLa cells transfected with wild-type and mutant pSPL3 expression vectors. The wild-type and c.2335-1G>C mutant alleles result in the complete inclusion or exclusion of exon 18, or in the activation of an exonic cryptic 5' ss with the consequent deletion of 169 bp at the end of this exon. However, we observed only a band of the expected size in normal thyroid tissue by RT-PCR. Additionally, the c.2335-1G>C mutation activates an unusual cryptic donor splice site in intron 17, located at position -14 of the authentic intron 17/exon 18 junction site, with an insertion of the last 14 nucleotides of the intron 17 in mutant transcripts with complete and partial inclusion of exon 18. The theoretical consequences of splice site mutation, predicted with the bioinformatics NNSplice, Fsplice, SPL, SPLM and MaxEntScan programs were investigated and evaluated in relation with the experimental evidence. These analyses confirm that c.2335-1G>C mutant allele would result in the abolition of the authentic splice acceptor site. The results suggest the coexistence in our patient of four putative truncated proteins of 786, 805, 806 and 1105 amino acids, with conservation of peroxidase-like domain and loss of gp91(phox)/NOX2-like domain. In conclusion a novel heterozygous compound was identified being responsible of IOD. Cryptic splicing sites have been characterized in DUOX2 gene for the first time. The use of molecular biology techniques is a valuable tool for understanding the molecular pathophysiology of this type of thyroid defects.

Complete androgen insensitivity syndrome caused by a novel splice donor site mutation and activation of a cryptic splice donor site in the androgen receptor gene

The androgen insensitivity syndrome is an X-linked recessive genetic disorder characterized by resistance to the actions of androgens in an individual with a male karyotype. We evaluated a 34-year-old female with primary amenorrhea and a 46,XY karyotype, with normal secondary sex characteristics, absence of uterus and ovaries, intra-abdominal testis, and elevated testosterone levels. Sequence analysis of the androgen receptor (AR) gene revealed a novel splice donor site mutation in intron 4 (c.2173+2T>C). RT-PCR analysis showed that this mutation resulted in the activation of a cryptic splice donor site located in the second half of exon 4 and in the synthesis of a shorter mRNA transcript and an in-frame deletion of 41 amino acids. This novel mutation associated with a rare mechanism of abnormal splicing further expands the spectrum of mutations associated with the androgen insensitivity syndrome and may contribute to the understanding of the molecular mechanisms involved in splicing defects.

Novel mutations in EVC cause aberrant splicing in Ellis-van Creveld syndrome

Ellis-van Creveld syndrome (EvC) is a rare autosomal recessive disorder characterized by disproportionate chondrodysplasia, postaxial polydactyly, nail dystrophy, dental abnormalities and in a proportion of patients, congenital cardiac malformations. Weyers acrofacial dysostosis (Weyers) is another dominantly inherited disorder allelic to EvC syndrome but with milder phenotypes. Both disorders can result from loss-of-function mutations in either EVC or EVC2 gene, and phenotypes associated with the two gene mutations are clinically indistinguishable. We present here a clinical and molecular analysis of a Chinese family manifested specific features of EvC syndrome. Sequencing of both EVC and EVC2 identified two novel heterozygous splice site mutations c.384+5G>C in intron 3 and c.1465-1G>A in intron 10 in EVC, which were inherited from mother and father, respectively. In vitro minigene expression assay, RT-PCR and sequencing analysis demonstrated that c.384+5G>C mutation abolished normal splice site and created a new cryptic acceptor site within exon 4, whereas c.1465-1G>A mutation affected consensus splice junction site and resulted in full exon 11 skipping. These two aberrant pre-mRNA splicing processes both produced in-frame abnormal transcripts that possibly led to abolishment of important functional domains. To our knowledge, this is the first report of EVC mutations that cause EvC syndrome in Chinese population. Our data revealed that EVC splice site mutations altered splicing pattern and helped elucidate the pathogenesis of EvC syndrome.

A novel mutation in the β-spectrin gene causes the activation of a cryptic 5'-splice site and the creation of a de novo 3'-splice site

The analysis of genes involved in hereditary spherocytosis, by next-generation sequencing in two patients with clinical diagnosis of the disease, showed the presence of the c.1795+1G>A mutation in the SPTB gene. cDNA amplification then revealed the occurrence of a consequent aberrant mRNA isoform produced from the activation of a cryptic 5′-splice site and the creation of a newly 3′-splice site. The mechanisms by which these two splice sites are used as a result of the same mutation should be analyzed in depth in further studies.

Enter exitrons

Exitrons are exon-like introns located within protein-coding exons. Removal or retention of exitrons through alternative splicing increases proteome complexity and thus adds to phenotypic diversity.

Generation and Characterization of a Transgenic Mouse Carrying a Functional Human β -Globin Gene with the IVSI-6 Thalassemia Mutation

Mouse models that carry mutations causing thalassemia represent a suitable tool to test in vivo new mutation-specific therapeutic approaches. Transgenic mice carrying the β-globin IVSI-6 mutation (the most frequent in Middle-Eastern regions and recurrent in Italy and Greece) are, at present, not available. We report the production and characterization of a transgenic mouse line (TG-β-IVSI-6) carrying the IVSI-6 thalassemia point mutation within the human β-globin gene. In the TG-β-IVSI-6 mouse (a) the transgenic integration region is located in mouse chromosome 7; (b) the expression of the transgene is tissue specific; (c) as expected, normally spliced human β-globin mRNA is produced, giving rise to β-globin production and formation of a human-mouse tetrameric chimeric hemoglobin ^muα-globin₂/^huβ-globin₂ and, more importantly, (d) the aberrant β-globin-IVSI-6 RNAs are present in blood cells. The TG-β-IVSI-6 mouse reproduces the molecular features of IVSI-6 β-thalassemia and might be used as an in vivo model to characterize the effects of antisense oligodeoxynucleotides targeting the cryptic sites responsible for the generation of aberrantly spliced β-globin RNA sequences, caused by the IVSI-6 mutation. These experiments are expected to be crucial for the development of a personalized therapy for β-thalassemia.

Unmasking alternative splicing inside protein-coding exons defines exitrons and their role in proteome plasticity

Alternative splicing (AS) diversifies transcriptomes and proteomes and is widely recognized as a key mechanism for regulating gene expression. Previously, in an analysis of intron retention events in Arabidopsis, we found unusual AS events inside annotated protein-coding exons. Here, we also identify such AS events in human and use these two sets to analyse their features, regulation, functional impact, and evolutionary origin. As these events involve introns with features of both introns and protein-coding exons, we name them exitrons (exonic introns). Though exitrons were detected as a subset of retained introns, they are clearly distinguishable, and their splicing results in transcripts with different fates. About half of the 1002 Arabidopsis and 923 human exitrons have sizes of multiples of 3 nucleotides (nt). Splicing of these exitrons results in internally deleted proteins and affects protein domains, disordered regions, and various post-translational modification sites, thus broadly impacting protein function. Exitron splicing is regulated across tissues, in response to stress and in carcinogenesis. Intriguingly, annotated intronless genes can be also alternatively spliced via exitron usage. We demonstrate that at least some exitrons originate from ancestral coding exons. Based on our findings, we propose a “splicing memory” hypothesis whereby upon intron loss imprints of former exon borders defined by vestigial splicing regulatory elements could drive the evolution of exitron splicing. Altogether, our studies show that exitron splicing is a conserved strategy for increasing proteome plasticity in plants and animals, complementing the repertoire of AS events.

BAP1 missense mutation c.2054 A>T (p.E685V) completely disrupts normal splicing through creation of a novel 5' splice site in a human mesothelioma cell line

BAP1 is a tumor suppressor gene that is lost or deleted in diverse cancers, including uveal mela¬noma, malignant pleural mesothelioma (MPM), clear cell renal carcinoma, and cholangiocarcinoma. Recently, BAP1 germline mutations have been reported in families with combinations of these same cancers. A particular challenge for mutation screening is the classification of non-truncating BAP1 sequence variants because it is not known whether these subtle changes can affect the protein function sufficiently to predispose to cancer development. Here we report mRNA splicing analysis on a homozygous substitution mutation, BAP1 c. 2054 A&T (p.Glu685Val), identified in an MPM cell line derived from a mesothelioma patient. The mutation occurred at the 3rd nucleotide from the 3' end of exon 16. RT-PCR, cloning and subsequent sequencing revealed several aberrant splicing products not observed in the controls: 1) a 4 bp deletion at the end of exon 16 in all clones derived from the major splicing product. The BAP1 c. 2054 A&T mutation introduced a new 5' splice site (GU), which resulted in the deletion of 4 base pairs and presumably protein truncation; 2) a variety of alternative splicing products that led to retention of different introns: introns 14-16; introns 15-16; intron 14 and intron 16; 3) partial intron 14 and 15 retentions caused by activation of alternative 3' splice acceptor sites (AG) in the introns. Taken together, we were unable to detect any correctly spliced mRNA transcripts in this cell line. These results suggest that aberrant splicing caused by this mutation is quite efficient as it completely abolishes normal splicing through creation of a novel 5' splice site and activation of cryptic splice sites. These data support the conclusion that BAP1 c.2054 A&T (p.E685V) variant is a pathogenic mutation and contributes to MPM through disruption of normal splicing.

Novel compound heterozygous Thyroglobulin mutations c.745+1G>A/c.7036+2T>A associated with congenital goiter and hypothyroidism in a Vietnamese family. Identification of a new cryptic 5' splice site in the exon 6

Several patients were identified with dyshormonogenesis caused by mutations in the thyroglobulin (TG) gene. These defects are inherited in an autosomal recessive manner and affected individuals are either homozygous or compound heterozygous for the mutations. The aim of the present study was to identify new TG mutations in a patient of Vietnamese origin affected by congenital hypothyroidism, goiter and low levels of serum TG. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a novel c.745+1G>A (g.IVS6 + 1G>A), whereas the hypothetical paternal mutation consists of a novel c.7036+2T>A (g.IVS40 + 2T>A). The father was not available for segregation analysis. Ex-vivo splicing assays and subsequent RT-PCR analyses were performed on mRNA isolated from the eukaryotic-cells transfected with normal and mutant expression vectors. Minigene analysis of the c.745+1G>A mutant showed that the exon 6 is skipped during pre-mRNA splicing or partially included by use of a cryptic 5' splice site located to 55 nucleotides upstream of the authentic exon 6/intron 6 junction site. The functional analysis of c.7036+2T>A mutation showed a complete skipping of exon 40. The theoretical consequences of splice site mutations, predicted with the bioinformatics tool NNSplice, Fsplice, SPL, SPLM and MaxEntScan programs were investigated and evaluated in relation with the experimental evidence. These analyses predicted that both mutant alleles would result in the abolition of the authentic splice donor sites. The c.745+1G>A mutation originates two putative truncated proteins of 200 and 1142 amino acids, whereas c.7036+2T>A mutation results in a putative truncated protein of 2277 amino acids. In conclusion, we show that the c.745+1G>A mutation promotes the activation of a new cryptic donor splice site in the exon 6 of the TG gene. The functional consequences of these mutations could be structural changes in the protein molecule that alter the biosynthesis of thyroid hormones.

Novel and atypical splicing mutation in a compound heterozygous UNC13D defect presenting in Familial Hemophagocytic Lymphohistiocytosis triggered by EBV infection

Familial Hemophagocytic Lymphohistiocytosis type 3 (FHL3) is a genetic disorder caused by mutations in UNC13D gene, coding the granule priming factor Munc13-4 that intervenes in NK and T cell cytotoxic function. Here we report the case of a 17-month-old girl with prolonged symptomatic EBV infectious mononucleosis and clinical symptoms of hemophagocytic syndrome. In vitro functional analysis pointed to a degranulation defect. The genetic analysis of UNC13D gene identified initially a heterozygous mutation (c.753+1G>T) in the donor splice-site that resulted in exon 9 skipping (maternal allele). Mutations in other genes were considered, but additional analysis of UNC13D cDNA revealed in the paternal allele a heterozygous transition from G to A (c.2448-13G>A) at the 3' acceptor splice-site in intron 25, generating a new acceptor splice-site that leads to a frameshift and a premature STOP codon. Allele specific amplification of the cDNA confirmed the absence of a functional mRNA from the paternal allele. This case illustrates an atypical compound heterozygous UNC13D mutation affecting the RNA splicing that generates a typical FHL3 phenotype.

An FBN1 deep intronic mutation in a familial case of Marfan syndrome: an explanation for genetically unsolved cases?

Marfan syndrome (MFS) is caused by mutations in the FBN1 (fibrillin-1) gene, but approximately 10% of MFS cases remain genetically unsolved. Here, we report a new FBN1 mutation in an MFS family that had remained negative after extensive molecular genomic DNA FBN1 testing, including denaturing high-performance liquid chromatography, Sanger sequencing, and multiplex ligation-dependent probe amplification. Linkage analysis in the family and cDNA sequencing of the proband revealed a deep intronic point mutation in intron 56 generating a new splice donor site. This mutation results in the integration of a 90-bp pseudo-exon between exons 56 and 57 containing a stop codon, causing nonsense-mediated mRNA decay. Although more than 90% of FBN1 mutations can be identified with regular molecular testing at the genomic level, deep intronic mutations will be missed and require cDNA sequencing or whole-genome sequencing.

Evolution of the human gastrokine locus and confounding factors regarding the pseudogenicity of GKN3

In a screen for genes expressed specifically in gastric mucous neck cells, we identified GKN3, the recently discovered third member of the gastrokine family. We present confirmatory mouse data and novel porcine data showing that mouse GKN3 expression is confined to mucous cells of the corpus neck and antrum base and is prominently expressed in metaplastic lesions. GKN3 was proposed originally to be expressed in some human populations and a pseudogene in others. To investigate that hypothesis, we studied human GKN3 evolution in the context of its paralogous genomic neighbors, GKN1 and GKN2. Haplotype analysis revealed that GKN3 mimics GKN2 in patterns of exonic SNP allocation, whereas GKN1 appeared to be more stringently selected. GKN3 showed signatures of both directional selection and population based selective sweeps in humans. One such selective sweep includes SNP rs10187256, originally identified as an ancestral tryptophan to premature STOP codon mutation. The derived (nonancestral) allele went to fixation in Asia. We show that another SNP, rs75578132, identified 5 bp downstream of rs10187256, exhibits a second selective sweep in almost all Europeans, some Latinos, and some Africans, possibly resulting from a reintroduction of European genes during African colonization. Finally, we identify a mutation that would destroy the splice donor site in the putative exon3-intron3 boundary, which occurs in all human genomes examined to date. Our results highlight a stomach-specific human genetic locus, which has undergone various selective sweeps across European, Asian, and African populations and thus reflects geographic and ethnic patterns in genome evolution.

Targeting RNA splicing for disease therapy

Splicing of pre-messenger RNA into mature messenger RNA is an essential step for the expression of most genes in higher eukaryotes. Defects in this process typically affect cellular function and can have pathological consequences. Many human genetic diseases are caused by mutations that cause splicing defects. Furthermore, a number of diseases are associated with splicing defects that are not attributed to overt mutations. Targeting splicing directly to correct disease-associated aberrant splicing is a logical approach to therapy. Splicing is a favorable intervention point for disease therapeutics, because it is an early step in gene expression and does not alter the genome. Significant advances have been made in the development of approaches to manipulate splicing for therapy. Splicing can be manipulated with a number of tools including antisense oligonucleotides, modified small nuclear RNAs (snRNAs), trans-splicing, and small molecule compounds, all of which have been used to increase specific alternatively spliced isoforms or to correct aberrant gene expression resulting from gene mutations that alter splicing. Here we describe clinically relevant splicing defects in disease states, the current tools used to target and alter splicing, specific mutations and diseases that are being targeted using splice-modulating approaches, and emerging therapeutics.

Interplay between DMD point mutations and splicing signals in Dystrophinopathy phenotypes

DMD nonsense and frameshift mutations lead to severe Duchenne muscular dystrophy while in-frame mutations lead to milder Becker muscular dystrophy. Exceptions are found in 10% of cases and the production of alternatively spliced transcripts is considered a key modifier of disease severity. Several exonic mutations have been shown to induce exon-skipping, while splice site mutations result in exon-skipping or activation of cryptic splice sites. However, factors determining the splicing pathway are still unclear. Point mutations provide valuable information regarding the regulation of pre-mRNA splicing and elements defining exon identity in the DMD gene. Here we provide a comprehensive analysis of 98 point mutations related to clinical phenotype and their effect on muscle mRNA and dystrophin expression. Aberrant splicing was found in 27 mutations due to alteration of splice sites or splicing regulatory elements. Bioinformatics analysis was performed to test the ability of the available algorithms to predict consequences on mRNA and to investigate the major factors that determine the splicing pathway in mutations affecting splicing signals. Our findings suggest that the splicing pathway is highly dependent on the interplay between splice site strength and density of regulatory elements.

IL-10R polymorphisms are associated with very-early-onset ulcerative colitis

BACKGROUND

Interleukin-10 (IL-10) signaling genes are attractive inflammatory bowel disease (IBD) candidate genes as IL-10 restricts intestinal inflammation, IL-10 polymorphisms have been associated with IBD in genome-wide association studies, and mutations in IL-10 and IL-10 receptor (IL-10R) genes have been reported in immunodeficient children with severe infantile-onset IBD. Our objective was to determine if IL-10R polymorphisms were associated with early-onset IBD (EO-IBD) and very-early-onset IBD (VEO-IBD).

METHODS

Candidate-gene analysis of IL10RA and IL10RB was performed after initial sequencing of an infantile onset-IBD patient identified a novel homozygous mutation. The discovery cohort included 188 EO-IBD subjects and 188 healthy subjects. Polymorphisms associated with IBD in the discovery cohort were genotyped in an independent validation cohort of 422 EO-IBD subjects and 480 healthy subjects.

RESULTS

We identified a homozygous, splice-site point mutation in IL10RA in an infantile-onset IBD patient causing a premature stop codon (P206X) and IL-10 insensitivity. IL10RA and IL10RB sequencing in the discovery cohort identified five IL10RA polymorphisms associated with ulcerative colitis (UC) and two IL10RB polymorphisms associated with Crohn's disease (CD). Of these polymorphisms, two IL10RA single nucleotide polymorphisms, rs2228054 and rs2228055, were associated with VEO-UC in the discovery cohort and replicated in an independent validation cohort (odds ratio [OR] 3.08, combined P = 2 x 10(-4); and OR 2.93, P = 6 x 10(-4), respectively).

CONCLUSIONS

We identified IL10RA polymorphisms that confer risk for developing VEO-UC. Additionally, we identified the first splice site mutation in IL10RA resulting in infantile-onset IBD. This study expands the phenotype of IL10RA polymorphisms to include both severe arthritis and VEO-UC.

The chicken frizzle feather is due to an α-keratin (KRT75) mutation that causes a defective rachis

Feathers have complex forms and are an excellent model to study the development and evolution of morphologies. Existing chicken feather mutants are especially useful for identifying genetic determinants of feather formation. This study focused on the gene F, underlying the frizzle feather trait that has a characteristic curled feather rachis and barbs in domestic chickens. Our developmental biology studies identified defects in feather medulla formation, and physical studies revealed that the frizzle feather curls in a stepwise manner. The frizzle gene is transmitted in an autosomal incomplete dominant mode. A whole-genome linkage scan of five pedigrees with 2678 SNPs revealed association of the frizzle locus with a keratin gene-enriched region within the linkage group E22C19W28_E50C23. Sequence analyses of the keratin gene cluster identified a 69 bp in-frame deletion in a conserved region of KRT75, an α-keratin gene. Retroviral-mediated expression of the mutated F cDNA in the wild-type rectrix qualitatively changed the bending of the rachis with some features of frizzle feathers including irregular kinks, severe bending near their distal ends, and substantially higher variations among samples in comparison to normal feathers. These results confirmed KRT75 as the F gene. This study demonstrates the potential of our approach for identifying genetic determinants of feather forms.

With the availability of a sequenced chicken genome, the reservoir of variant plumage genes found in domestic chickens can provide insight into the molecular mechanisms underlying the diversity of feather forms. In this paper, we identify the molecular basis of the distinctive frizzle (F) feather phenotype that is caused by a single autosomal incomplete dominant gene in which heterozygous individuals show less severe phenotypes than homozygous individuals. Feathers in frizzle chickens curve backward. We used computer-assisted analysis to establish that the rachis of the frizzle feather was irregularly kinked and more severely bent than normal. Moreover, microscopic evaluation of regenerating feathers found reduced proliferating cells that give rise to the frizzle rachis. Analysis of a pedigree of frizzle chickens showed that the phenotype is linked to two single-nucleotide polymorphisms in a cluster of keratin genes within the linkage group E22C19W28_E50C23. Sequencing of the gene cluster identified a 69-base pair in-frame deletion of the protein coding sequence of the α-keratin-75 gene. Forced expression of the mutated gene in normal chickens produced a twisted rachis. Although chicken feathers are primarily composed of beta-keratins, our findings indicate that alpha-keratins have an important role in establishing the structure of feathers.

Re-splicing of mature mRNA in cancer cells promotes activation of distant weak alternative splice sites

Transcripts of the human tumor susceptibility gene 101 (TSG101) are aberrantly spliced in many cancers. A major aberrant splicing event on the TSG101 pre-mRNA involves joining of distant alternative 5′ and 3′ splice sites within exon 2 and exon 9, respectively, resulting in the extensive elimination of the mRNA. The estimated strengths of the alternative splice sites are much lower than those of authentic splice sites. We observed that the equivalent aberrant mRNA could be generated from an intron-less TSG101 gene expressed ectopically in breast cancer cells. Remarkably, we identified a pathway-specific endogenous lariat RNA consisting solely of exonic sequences, predicted to be generated by a re-splicing between exon 2 and exon 9 on the spliced mRNA. Our results provide evidence for a two-step splicing pathway in which the initial constitutive splicing removes all 14 authentic splice sites, thereby bringing the weak alternative splice sites into close proximity. We also demonstrate that aberrant multiple-exon skipping of the fragile histidine triad (FHIT) pre-mRNA in cancer cells occurs via re-splicing of spliced FHIT mRNA. The re-splicing of mature mRNA can potentially generate mutation-independent diversity in cancer transcriptomes. Conversely, a mechanism may exist in normal cells to prevent potentially deleterious mRNA re-splicing events.

Complexities of 5'splice site definition: implications in clinical analyses

In higher eukaryotes, the 5' splice site (5'ss) is initially recognized through an RNA-RNA interaction by U1 small nuclear ribonucleoprotein (U1 snRNP). This event represents one of the key steps in initial spliceosomal assembly and many disease-associated mutations in humans often disrupt this process. Beside base pair complementarity, 5'ss recognition can also be modified by additional factors such as RNA secondary structures or the specific binding of other nuclear proteins. In this work, we have focused on investigating a few examples of changes detected within the 5'ss in patients, that would not be immediately considered "disease causing mutations". We show that the splicing outcome of very similar mutations can be very different due to variations in trans-acting factor(s) interactions and specific context influences. Using several NF1 donor sites and SELEX approaches as experimental models, we have examined the binding properties of particular sequence motifs such as GGGU found in donor sites, and how the sequence context can change their interaction with hnRNPs such as H/F and A1/A2. Our results clearly show that even minor differences in local nucleotide context can differentially affect the binding ability of these factors to the GGGU core. Finally, using a previously identified mutation in KCNH2 that resulted in intron retention we show how very similar 5'ss mutations found in patients can have a very different splicing outcome due to the neighbouring sequence context, thus highlighting the general need to approach splicing problems with suitable experimental approaches.