Order of intron removal influences multiple splice outcomes, including a two-exon skip, in a COL5A1 acceptor-site mutation that results in abnormal pro-alpha1(V) N-propeptides and Ehlers-Danlos syndrome type I

Identification of a novel splice variant in SEC23B gene in a patient with concomitant presence of congenital dyserythropoietic anemia II and Gilbert's syndrome

BACKGROUND

Congenital dyserythropoietic anemia Ⅱ (CDA Ⅱ) is a rare inherited disorder of defective erythropoiesis caused by SEC23B gene mutation. CDA Ⅱ is often misdiagnosed as a more common type of clinically related anemia, or it remains undiagnosed due to phenotypic variability caused by the coexistence of inherited liver diseases, including Gilbert's syndrome (GS) and hereditary hemochromatosis.

METHODS

We describe the case of a boy with genetically undetermined severe hemolytic anemia, hepatosplenomegaly, and gallstones whose diagnosis was achieved by targeted next generation sequencing.

RESULTS

Molecular analysis revealed a maternally inherited novel intronic variant and a paternally inherited missense variant, c.[994-3C > T];[1831C > T] in the SEC23B gene, confirming diagnosis of CDA Ⅱ. cDNA analysis verified that the splice acceptor site variant results in two mutant transcripts, one with an exon 9 skip and one in which exons 9 and 10 are deleted. SEC23B mRNA levels in the patient were lower than those in healthy controls. The patient was also homozygous for the UGT1A1*6 allele, consistent with GS.

CONCLUSION

Identification of the novel splice variant in this study further expands the spectrum of known SEC23B gene mutations. Molecular genetic approaches can lead to accurate diagnosis and management of CDA Ⅱ patients, particularly for those with GS coexisting.

Single-Base Substitution Causing Dual-Exon Skipping Event in PKD2 Gene: Unusual Molecular Finding from Exome Sequencing in a Patient with Autosomal Dominant Polycystic Kidney Disease

Background: Pathogenic variants in the Polycystic Kidney Disease 2 (PKD2) gene are associated with Autosomal Dominant Polycystic Kidney Disease (ADPKD) in approximately 30% of cases. In recent years, the high-throughput sequencing techniques have significantly increased the number of variants identified in affected patients. Here, we described the peculiar effect of a PKD2 splicing variant, the c.1717-2A>G, identified in an Italian male patient with ADPKD. This variant led to the unusual and rare skipping of two consecutive exons, causing a large in-frame deletion. Methods: The genetic evaluation of the patient was performed using the Next-Generation Sequencing (NGS) assay Clinical Exome Solution® (SOPHiA Genetics). Bioinformatics analysis was performed using the SOPHiA DDM platform (SOPHiA Genetics). Prediction of pathogenicity was carried out by integrating several in silico tools. RNA evaluation was performed to test the effect of the variant on the PKD2 splicing using a Reverse-Transcription PCR coupled with cDNA sequencing. Results: NGS revealed the presence of the PKD2 c.1717-2A>G variant that lies in the canonical splice site of intron 7. This rare variant was predicted to have a significant impact on the splicing, proved by the RNA-based analysis. We identified the presence of a transcript characterised by the simultaneous skipping of exons 8 and 9, with a retained reading frame and the merging of exons 7-10. Conclusions: We described for the first time a dual-exon skip event related to the presence of a single-base substitution in the PKD2 gene in an ADPKD-affected patient. We assumed that the molecular basis of such a rare mechanism lies in the specific order of intron removal. The finding represents novel evidence of an alternative and unusual splicing mechanism in the PKD2 gene, adding insights to the pathogenesis of the ADPKD.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays

BACKGROUND

Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes.

RESULTS

We benchmark eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compare experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, is lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieve the best overall performance at distinguishing disruptive and neutral variants, and controlling for overall call rate genome-wide, SpliceAI and Pangolin have superior sensitivity. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues.

CONCLUSION

SpliceAI and Pangolin show the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

Case report: Genetic analysis of a novel intronic inversion variant in the SPTB gene associated with hereditary spherocytosis

Background: Hereditary spherocytosis (HS) is a congenital haemolytic anaemia attributed to dysregulation or abnormal quantities of erythrocyte membrane proteins. Currently, the most common erythrocytic gene, spectrin β (SPTB), variants are located in exons and give rise to mRNA defects. However, the genetic characteristics and pathogenic mechanisms of SPTB intronic variants are not completely understood. This study aimed to analyse a rare intronic inversion variant in the SPTB gene associated with HS, and explore the impact of the variant on SPTB mRNA splicing. Method: The clinical manifestations of the patient were summarised and analysed for spherocytosis phenotype diagnosis. The pathogenic variant was identified in the proband using targeted next-generation and Sanger sequencing. RNA sequencing was performed to analyse whether SPTB gene splicing and expression were affected. Results: Targeted next-generation sequencing identified a novel disease-associated intronic inversion variant of the SPTB gene in the proband. The inversion variant was located between intron 19 and 20, and contained the entire exon 20 and partial sequences of adjacent introns. Sanger sequencing confirmed that the intronic inversion variant only appeared in the genome of the proband, not in his parents. RNA sequencing revealed that the variant could result in the skipping of exon 20 and reduced expression of SPTB mRNA. Conclusion: This study identifies a rare intronic inversion variant in the SPTB gene associated with hereditary spherocytosis. The pathogenic variant can lead to exon 20 skipping and decreased SPTB gene expression. This finding has not been previously reported in any literature. This study can expand the intronic variant spectrum of the SPTB gene, deepen our understanding of HS pathogenesis, and contribute to the genetic diagnosis and clinical management of patients.

Precision medicine using whole genome sequencing in a cat identifies a novel COL5A1 variant for classical Ehlers-Danlos syndrome

BACKGROUND

Ehlers-Danlos syndromes (EDS) are a heterogeneous group of heritable connective tissue disorders occurring in both human and veterinary patients. The genetics of these disorders are poorly described in small animal patients.

HYPOTHESIS/OBJECTIVES

Define the clinical manifestations and genetic cause of a suspected form of EDS in a cat.

ANIMALS

A 14-week-old male domestic medium hair cat was presented with skin hyperextensibility and fragility. The classic tragic facial expression was observed as well as chronic pruritus and mild hyperesthesia.

METHODS

Blood samples and a skin biopsy sample were collected from the affected cat. Clinical examinations, histology, electron microscopy and whole genome sequencing were conducted to characterize the clinical presentation and identify possible pathogenic DNA variants to support a diagnosis. Criteria defining variant pathogenicity were examined including human disease variant databases.

RESULTS

Histology showed sparse, disorganized collagen and an increase in cutaneous mast cells. Electron microscopy identified ultrastructural defects commonly seen in collagen type V alpha 1 chain (COL5A1) variants including flower-like collagen fibrils in cross-section. Whole genome sequencing and comparison with 413 cats in the 99 Lives Cat Genome Sequencing Consortium database identified a novel splice acceptor site variant at exon 4 in COL5A1 (c.501-2A>C).

CONCLUSIONS AND CLINICAL IMPORTANCE

Our report broadens the current understanding of EDS in veterinary patients and supports the use of precision medicine techniques in clinical veterinary practice. The classification of variants for pathogenicity should be considered in companion animals.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays

Background

Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes.

Results

We benchmarked eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compared experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, was lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieved the best overall performance at distinguishing disruptive and neutral variants. Controlling for overall call rate genome-wide, SpliceAI and Pangolin also showed superior overall sensitivity for identifying SDVs. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues.

Conclusion

SpliceAI and Pangolin showed the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

Clinical features and molecular characterization of Chinese patients with FKBP10 variants

BACKGROUND

Osteogenesis imperfecta (OI) is a group of rare skeletal dysplasia. Long bone deformity and scoliosis are often associated with progressively deforming types of OI. FKBP65 (encoded by FKBP10, OMIM *607063) plays a crucial role in the processing of type I procollagen. Autosomal recessive variants in FKBP10 result in type XI osteogenesis imperfecta.

METHODS

Patients diagnosed with OI were recruited for a genetic test. RT-PCR and Sanger sequencing were applied to confirm the splicing defect in FKBP10 mRNA with the splice-site variant. The bone structure was characterized by Goldner's trichrome staining. Bioinformatic analyses of bulk RNA sequencing data were performed to examine the effect of the FKBP10 variant on gene expression.

RESULTS

Here we reported three children from a consanguineous family harboured a homozygous splice-site variant (c.918-3C > G) in FKBP10 intron and developed long bone deformity and early onset of scoliosis. We also observed frequent long bone fractures and spinal deformity in another 3 OI patients with different FKBP10 variants. The homozygous splicing variant identified in the fifth intron of FKBP10 (c.918-3C > G) led to abnormal RNA processing and loss of FKBP65 protein and consequently resulted in aberrant collagen alignment and porous bone morphology. Analysis of transcriptomic data indicated that genes involved in protein processing and osteoblast differentiation were significantly affected in the patient-derived osteoblasts.

CONCLUSION

Our study characterized the clinical features of OI patients with FKBP10 variants and revealed the pathogenesis of the c.918-3C > G variant. The molecular analyses helped to gain insight into the deleterious effects of FKBP10 variants on collagen processing and osteoblast differentiation.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

Identification of two aberrant transcripts by RNA sequencing for a novel variant c.3354 + 5 G > A of MED12 in a Chinese girl with non-syndromic intellectual disability

BACKGROUND

Missense variants in MED12 are associated with MED12-related disorders. We aimed to clarify the molecular level changes and underlying pathogenic mechanism of a female patient in our study.

METHODS

We reported a Chinese girl with clinical characteristics similar to MED12-related disorders. Trio whole exome sequencing (WES) was performed to identify related pathogenic variant(s) and RNA sequencing (RNA-seq) was subsequently applied to evaluate the effect of identified variant(s) on mRNA splicing. Moreover, X-chromosome inactivation (XCI) assay based on AR and RP2 was performed to reveal the XCI pattern of the female patient.

RESULTS

The proband manifested mainly as mental retardation and language impairment. Trio WES revealed a novel heterozygous variant c.3354 + 5 G > A in intron 23 of MED12. RNA-seq identified two aberrant transcripts. XCI assay on AR revealed a homozygous result, while XCI based on RP2 showed random pattern in peripheral blood.

CONCLUSION

In conclusion, we identified a novel variant c.3354 + 5 G > A by WES combined with RNA-seq, which extends the spectrum of MED12 variants and provide a basis for further genetic counseling. According to the result of two aberrant transcripts by RNA-seq, we speculate that our patient's milder clinical feature may be the consequence of multiple different transcripts.

Donor Splice Site Variant in SLC9A6 Causes Christianson Syndrome in a Lithuanian Family: A Case Report

Background and Objectives: The pathogenic variants of SLC9A6 are a known cause of a rare, X-linked neurological disorder called Christianson syndrome (CS). The main characteristics of CS are developmental delay, intellectual disability, and neurological findings. This study investigated the genetic basis and explored the molecular changes that led to CS in two male siblings presenting with intellectual disability, epilepsy, behavioural problems, gastrointestinal dysfunction, poor height, and weight gain. Materials and Methods: Next-generation sequencing of a tetrad was applied to identify the DNA changes and Sanger sequencing of proband’s cDNA was used to evaluate the impact of a splice site variant on mRNA structure. Bioinformatical tools were used to investigate SLC9A6 protein structure changes. Results: Sequencing and bioinformatical analysis revealed a novel donor splice site variant (NC_000023.11(NM_001042537.1):c.899 + 1G > A) that leads to a frameshift and a premature stop codon. Protein structure modelling showed that the truncated protein is unlikely to form any functionally relevant SLC9A6 dimers. Conclusions: Molecular and bioinformatical analysis revealed the impact of a novel donor splice site variant in the SLC9A6 gene that leads to truncated and functionally disrupted protein causing the phenotype of CS in the affected individuals.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

The upstream 5' splice site remains associated to the transcription machinery during intron synthesis

In the earliest step of spliceosome assembly, the two splice sites flanking an intron are brought into proximity by U1 snRNP and U2AF along with other proteins. The mechanism that facilitates this intron looping is poorly understood. Using a CRISPR interference-based approach to halt RNA polymerase II transcription in the middle of introns in human cells, we discovered that the nascent 5′ splice site base pairs with a U1 snRNA that is tethered to RNA polymerase II during intron synthesis. This association functionally corresponds with splicing outcome, involves bona fide 5′ splice sites and cryptic intronic sites, and occurs transcriptome-wide. Overall, our findings reveal that the upstream 5′ splice sites remain attached to the transcriptional machinery during intron synthesis and are thus brought into proximity of the 3′ splice sites; potentially mediating the rapid splicing of long introns.

We know that most splicing reactions take place co-transcriptionally, but how the transcription machinery facilitate splicing of introns is unknown. Here the authors show that the 5′ splice site remains associated with the transcription machinery during intron synthesis through U1 snRNP, providing a basis for the rapid splicing reaction of introns.

Liver Failure of Wilson's Disease With Manifestations Similar to Porphyria and Uncommon ATP7B Gene Mutation: A Case Report and Literature Review

Background: Wilson's disease (WD) is a rare condition; its diagnosis is challenging owing to a wide spectrum of ATP7B genotypes and variable clinical phenotypes, along with environmental factors. Few cases of WD with presentation of skin lesions and acute neurovisceral symptoms have been reported in the literature. To our knowledge, this is the first reported case of WD with an uncommon ATP7B gene mutation and rare symptoms of photosensitivity, sensation abnormality, and skin eruption occurring in a 19-year-old woman. Case presentation: We report the case of a 19-year-old woman with WD presenting with liver failure, skin manifestations, and acute neurovisceral symptoms.The rare mutation in intron 1 of ATP7B (c.51+2T > G) was further confirmed by gene sequencing. The patients' symptoms improved after administration of penicillamine and zinc therapy combined with plasma exchange. She received long-term penicillamine treatment, and her liver function was within the normal range at 1 year after discharge. However, she underwent liver transplantation at 1.5 years after discharge. Conclusions: We present a case of WD with a novel ATP7B gene mutation that may serve as a reference to generalists and specialists in hepatology or neurology of the rare clinical characteristics of WD, to prevent misdiagnosis and aid in the early diagnosis and treatment of the condition.

Clinical Characteristics of POC1B-Associated Retinopathy and Assignment of Pathogenicity to Novel Deep Intronic and Non-Canonical Splice Site Variants

Mutations in POC1B are a rare cause of inherited retinal degeneration. In this study, we present a thorough phenotypic and genotypic characterization of three individuals harboring putatively pathogenic variants in the POC1B gene. All patients displayed a similar, slowly progressive retinopathy (cone dystrophy or cone-rod dystrophy) with normal funduscopy but disrupted outer retinal layers on optical coherence tomography and variable age of onset. Other symptoms were decreased visual acuity and photophobia. Whole genome sequencing revealed a novel homozygous frameshift variant in one patient. Another patient was shown to harbor a novel deep intronic variant in compound heterozygous state with a previously reported canonical splice site variant. The third patient showed a novel nonsense variant and a novel non-canonical splice site variant. We aimed to validate the effect of the deep intronic variant and the non-canonical splice site variant by means of in vitro splice assays. In addition, direct RNA analysis was performed in one patient. Splicing analysis revealed that the non-canonical splice site variant c.561-3T>C leads to exon skipping while the novel deep intronic variant c.1033-327T>A causes pseudoexon activation. Our data expand the genetic landscape of POC1B mutations and confirm the benefit of genome sequencing in combination with downstream functional validation using minigene assays for the analysis of putative splice variants. In addition, we provide clinical multimodal phenotyping of the affected individuals.

The Exon Junction Complex and intron removal prevent re-splicing of mRNA

Accurate splice site selection is critical for fruitful gene expression. Recently, the mammalian EJC was shown to repress competing, cryptic, splice sites (SS). However, the evolutionary generality of this remains unclear. Here, we demonstrate the Drosophila EJC suppresses hundreds of functional cryptic SS, even though most bear weak splicing motifs and are seemingly incompetent. Mechanistically, the EJC directly conceals cryptic splicing elements by virtue of its position-specific recruitment, preventing aberrant SS definition. Unexpectedly, we discover the EJC inhibits scores of regenerated 5' and 3' recursive SS on segments that have already undergone splicing, and that loss of EJC regulation triggers faulty resplicing of mRNA. An important corollary is that certain intronless cDNA constructs yield unanticipated, truncated transcripts generated by resplicing. We conclude the EJC has conserved roles to defend transcriptome fidelity by (1) repressing illegitimate splice sites on pre-mRNAs, and (2) preventing inadvertent activation of such sites on spliced segments.

Clinical and genetic analysis of classical Ehlers-Danlos syndrome patient caused by synonymous mutation in COL5A2

Background

Classical Ehlers‐Danlos syndrome (cEDS) is a heterogeneous connective tissue disorder that mainly results from the germline mutation of COL5A1 and COL5A2. The majority of the COL5A2 mutations reported to date represent structural mutations, including missense or in‐frame exon‐skipping splice mutations. The only reported synonymous mutation was expected to affect on splicing of exon 29 by prediction programs which should be further confirmed.

Methods

Whole exome sequencing was performed to identify the genetic variants of a Chinese boy who was characterized by skin hyperextensibility, abnormal scarring, hypermobile joints and scoliosis. Sanger sequencing was used to validate the variants in his parents. Reverse transcription polymerase chain reaction (RT‐PCR) was performed to analyze the functional effects of the variant.

Results

A de novo heterozygous synonymous variant (NM_000393.5:c.1977 G>A) of COL5A2 gene was identified in the patient. The results of RT‐PCR revealed that the synonymous variant led to skipping of exon 29 in the RNA transcript.

Conclusions

Our study supplies further supporting evidence that the synonymous COL5A2 mutation c.1977 G>A can cause skipping of exon 29 in the RNA transcript, thus resulting in the production of mutant α2(V)‐chains and clinical phenotype of cEDS. This result highlights the need to include splicing‐altering synonymous mutations into the screening for cEDS.

Revealing nascent RNA processing dynamics with nano-COP

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

Hybrid Minigene Assay: An Efficient Tool to Characterize mRNA Splicing Profiles of NF1 Variants

Neurofibromatosis type 1 (NF1) is caused by heterozygous loss of function mutations in the NF1 gene. Although patients are diagnosed according to clinical criteria and few genotype-phenotype correlations are known, molecular analysis remains important. NF1 displays allelic heterogeneity, with a high proportion of variants affecting splicing, including deep intronic alleles and changes outside the canonical splice sites, making validation problematic. Next Generation Sequencing (NGS) technologies integrated with multiplex ligation-dependent probe amplification (MLPA) have largely overcome RNA-based techniques but do not detect splicing defects. A rapid minigene-based system was set up to test the effects of NF1 variants on splicing. We investigated 29 intronic and exonic NF1 variants identified in patients during the diagnostic process. The minigene assay showed the coexistence of multiple mechanisms of splicing alterations for seven variants. A leaky effect on splicing was documented in one de novo substitution detected in a sporadic patient with a specific phenotype without neurofibromas. Our splicing assay proved to be a reliable and fast method to validate novel NF1 variants potentially affecting splicing and to detect hypomorphic effects that might have phenotypic consequences, avoiding the requirement of patient's RNA.

Calculating the most likely intron splicing orders in S. pombe, fruit fly, Arabidopsis thaliana, and humans

Background

Introns have been shown to be spliced in a defined order, and this order influences both alternative splicing regulation and splicing fidelity, but previous studies have only considered neighbouring introns. The detailed intron splicing order remains unknown.

Results

In this work, a method was developed that can calculate the intron splicing orders of all introns in each transcript. A simulation study showed that this method can accurately calculate intron splicing orders. I further applied this method to real S. pombe, fruit fly, Arabidopsis thaliana, and human sequencing datasets and found that intron splicing orders change from gene to gene and that humans contain more not in-order spliced transcripts than S. pombe, fruit fly and Arabidopsis thaliana. In addition, I reconfirmed that the first introns in humans are spliced slower than those in S. pombe, fruit fly, and Arabidopsis thaliana genome-widely. Both the calculated most likely orders and the method developed here are available on the web.

Conclusions

A novel computational method was developed to calculate the intron splicing orders and applied the method to real sequencing datasets. I obtained intron splicing orders for hundreds or thousands of genes in four organisms. I found humans contain more number of not in-order spliced transcripts.

Ovine congenital progressive muscular dystrophy (OCPMD) is a model of TNNT1 congenital myopathy

Ovine congenital progressive muscular dystrophy (OCPMD) was first described in Merino sheep flocks in Queensland and Western Australia in the 1960s and 1970s. The most prominent feature of the disease is a distinctive gait with stiffness of the hind limbs that can be seen as early as 3 weeks after birth. The disease is progressive. Histopathological examination had revealed dystrophic changes specifically in type I (slow) myofibres, while electron microscopy had demonstrated abundant nemaline bodies. Therefore, it was never certain whether the disease was a dystrophy or a congenital myopathy with dystrophic features. In this study, we performed whole genome sequencing of OCPMD sheep and identified a single base deletion at the splice donor site (+ 1) of intron 13 in the type I myofibre-specific TNNT1 gene (KT218690 c.614 + 1delG). All affected sheep were homozygous for this variant. Examination of TNNT1 splicing by RT-PCR showed intron retention and premature termination, which disrupts the highly conserved 14 amino acid C-terminus. The variant did not reduce TNNT1 protein levels or affect its localization but impaired its ability to modulate muscle contraction in response to Ca2+ levels. Identification of the causative variant in TNNT1 finally clarifies that the OCPMD sheep is in fact a large animal model of TNNT1 congenital myopathy. This model could now be used for testing molecular or gene therapies.

Exploring new genetic variants within COL5A1 intron 4-exon 5 region and TGF-β family with risk of anterior cruciate ligament ruptures

Variants within genes encoding structural and regulatory elements of ligaments have been associated with musculoskeletal soft tissue injury risk. The role of intron 4-exon 5 variants within the α1 chain of type V collagen (COL5A1) gene and genes of the transforming growth factor-β (TGF-β) family, TGFBR3 and TGFBI, was investigated on the risk of anterior cruciate ligament (ACL) ruptures. A case-control genetic association study was performed on 210 control (CON) and 249 participants with surgically diagnosed ruptures (ACL), of which 147 reported a noncontact mechanism of injury (NON). Whole-exome sequencing data were used to prioritize variants of potential functional relevance. Genotyping for COL5A1 (rs3922912 G>A, rs4841926 C>T, and rs3124299 C>T), TGFBR3 (rs1805113 G>A and rs1805117 T>C), and TGFBI (rs1442 G>C) was performed using Taqman SNP genotyping assays. Significant overrepresentation of the G allele of TGFBR3 rs1805113 was observed in CON vs ACL (P = .014) and NON groups (P = .021). Similar results were obtained in a female with the G allele (CON vs ACL: P = .029; CON vs NON: P = .016). The TGFBI rs1442 CC genotype was overrepresented in the female ACL vs CON (P = .013). Associations of inferred allele combinations were observed in line with the above results. COL5A1 intron 4-exon 5 genomic interval was not associated with the risk of ACL ruptures. Instead, this novel study is the first to use this approach to identify variants within the TGF-β signaling pathway to be implicated in the risk of ACL ruptures. A genetic susceptibility interval was identified to be explored in the context of extracellular matrix remodeling.

Splicing mutations in inherited retinal diseases

Mutations which induce aberrant transcript splicing represent a distinct class of disease-causing genetic variants in retinal disease genes. Such mutations may either weaken or erase regular splice sites or create novel splice sites which alter exon recognition. While mutations affecting the canonical GU-AG dinucleotides at the splice donor and splice acceptor site are highly predictive to cause a splicing defect, other variants in the vicinity of the canonical splice sites or those affecting additional cis-acting regulatory sequences within exons or introns are much more difficult to assess or even to recognize and require additional experimental validation. Splicing mutations are unique in that the actual outcome for the transcript (e.g. exon skipping, pseudoexon inclusion, intron retention) and the encoded protein can be quite different depending on the individual mutation. In this article, we present an overview on the current knowledge about and impact of splicing mutations in inherited retinal diseases. We introduce the most common sub-classes of splicing mutations including examples from our own work and others and discuss current strategies for the identification and validation of splicing mutations, as well as therapeutic approaches, open questions, and future perspectives in this field of research.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Genetics of Wilson disease and Wilson-like phenotype in a clinical series from eastern Spain

Wilson's disease (WD) is an autosomal recessive disorder caused by ATP7B mutations. Subjects with only one mutation may show clinical signs and individuals with biallelic changes may remain asymptomatic. We aimed to achieve a conclusive genetic diagnosis for 34 patients clinically diagnosed of WD. Genetic analysis comprised from analysis of exons to WES (whole exome sequencing), including promoter, introns, UTRs (untranslated regions), besides of study of large deletions/duplications by MLPA (multiplex ligation-dependent probe amplification). Biallelic ATP7B mutations were identified in 30 patients, so that four patients were analyzed using WES. Two affected siblings resulted to be compound heterozygous for mutations in CCDC115, which is involved in a form of congenital disorder of glycosylation. In sum, the majority of patients with a WD phenotype carry ATP7B mutations. However, if genetic diagnosis is not achieved, additional genes should be considered because other disorders may mimic WD.

Machine Learning Approaches for the Prioritization of Genomic Variants Impacting Pre-mRNA Splicing

Defects in pre-mRNA splicing are frequently a cause of Mendelian disease. Despite the advent of next-generation sequencing, allowing a deeper insight into a patient's variant landscape, the ability to characterize variants causing splicing defects has not progressed with the same speed. To address this, recent years have seen a sharp spike in the number of splice prediction tools leveraging machine learning approaches, leaving clinical geneticists with a plethora of choices for in silico analysis. In this review, some basic principles of machine learning are introduced in the context of genomics and splicing analysis. A critical comparative approach is then used to describe seven recent machine learning-based splice prediction tools, revealing highly diverse approaches and common caveats. We find that, although great progress has been made in producing specific and sensitive tools, there is still much scope for personalized approaches to prediction of variant impact on splicing. Such approaches may increase diagnostic yields and underpin improvements to patient care.

Mapping the Tooth Enamel Proteome and Amelogenin Phosphorylation Onto Mineralizing Porcine Tooth Crowns

Tooth enamel forms in an ephemeral protein matrix where changes in protein abundance, composition and posttranslational modifications are critical to achieve healthy enamel properties. Amelogenin (AMELX) with its splice variants is the most abundant enamel matrix protein, with only one known phosphorylation site at serine 16 shown in vitro to be critical for regulating mineralization. The phosphorylated form of AMELX stabilizes amorphous calcium phosphate, while crystalline hydroxyapatite forms in the presence of the unphosphorylated protein. While AMELX regulates mineral transitions over space and time, it is unknown whether and when un-phosphorylated amelogenin occurs during enamel mineralization. This study aims to reveal the spatiotemporal distribution of the cleavage products of the most abundant AMLEX splice variants including the full length P173, the shorter leucine-rich amelogenin protein (LRAP), and the exon 4-containing P190 in forming enamel, all within the context of the changing enamel matrix proteome during mineralization. We microsampled permanent pig molars, capturing known stages of enamel formation from both crown surface and inner enamel. Nano-LC-MS/MS proteomic analyses after tryptic digestion rendered more than 500 unique protein identifications in enamel, dentin, and bone. We mapped collagens, keratins, and proteolytic enzymes (CTSL, MMP2, MMP10) and determined distributions of P173, LRAP, and P190 products, the enamel proteins enamelin (ENAM) and ameloblastin (AMBN), and matrix-metalloprotease-20 (MMP20) and kallikrein-4 (KLK4). All enamel proteins and KLK4 were near-exclusive to enamel and in excellent agreement with published abundance levels. Phosphorylated P173 and LRAP products decreased in abundance from recently deposited matrix toward older enamel, mirrored by increasing abundances of testicular acid phosphatase (ACPT). Our results showed that hierarchical clustering analysis of secretory enamel links closely matching distributions of unphosphorylated P173 and LRAP products with ACPT and non-traditional amelogenesis proteins, many associated with enamel defects. We report higher protein diversity than previously published and Gene Ontology (GO)-defined protein functions related to the regulation of mineral formation in secretory enamel (e.g., casein α-S1, CSN1S1), immune response in erupted enamel (e.g., peptidoglycan recognition protein, PGRP), and phosphorylation. This study presents a novel approach to characterize and study functional relationships through spatiotemporal mapping of the ephemeral extracellular matrix proteome.

Context matters: Regulation of splice donor usage

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Systematic analysis of dark and camouflaged genes reveals disease-relevant genes hiding in plain sight

BACKGROUND

The human genome contains "dark" gene regions that cannot be adequately assembled or aligned using standard short-read sequencing technologies, preventing researchers from identifying mutations within these gene regions that may be relevant to human disease. Here, we identify regions with few mappable reads that we call dark by depth, and others that have ambiguous alignment, called camouflaged. We assess how well long-read or linked-read technologies resolve these regions.

RESULTS

Based on standard whole-genome Illumina sequencing data, we identify 36,794 dark regions in 6054 gene bodies from pathways important to human health, development, and reproduction. Of these gene bodies, 8.7% are completely dark and 35.2% are ≥ 5% dark. We identify dark regions that are present in protein-coding exons across 748 genes. Linked-read or long-read sequencing technologies from 10x Genomics, PacBio, and Oxford Nanopore Technologies reduce dark protein-coding regions to approximately 50.5%, 35.6%, and 9.6%, respectively. We present an algorithm to resolve most camouflaged regions and apply it to the Alzheimer's Disease Sequencing Project. We rescue a rare ten-nucleotide frameshift deletion in CR1, a top Alzheimer's disease gene, found in disease cases but not in controls.

CONCLUSIONS

While we could not formally assess the association of the CR1 frameshift mutation with Alzheimer's disease due to insufficient sample-size, we believe it merits investigating in a larger cohort. There remain thousands of potentially important genomic regions overlooked by short-read sequencing that are largely resolved by long-read technologies.

New variants in COL5A1 gene among Polish patients with Ehlers-Danlos syndrome: analysis of nine cases

Introduction

The Ehlers-Danlos syndrome (EDS) is a non-inflammatory, heritable connective tissue disorder divided into 13 types according to the 2017 International Classification of the Ehlers-Danlos syndromes. One of the subtypes of EDS, classical (cEDS), is characterized by joint hypermobility, skin hyperextensibility and atrophic scars, which are major criteria of cEDS.

Aim

In this study, the first in Central Eastern Europe, 44 patients were investigated. All of them were tested for COL5A1 mutations with direct DNA sequencing.

Material and methods

The study group included 44 patients of Polish origin, all of whom fulfilled criteria for the classical type of Ehlers-Danlos syndrome. Direct sequencing of the COL5A1, COL5A2 and COL1A1 c.934C>T genes was performed for all of them. Evaluation of potential pathogenicity of detected missense mutation was conducted using SIFT (Sorting Intolerant from Tolerant), PolyPhen-2, AlignGVGD (Align Grantham Variance/Grantham Difference). The effect of the splice site mutations was predicted by Human Splicing Finder and NetGene2 tools.

Results

Among all tested patients, nine mutations of COL5A1 gene were detected (8 missense mutations and 1 splice site). The alterations identified by us are new, hitherto not described in other reports. Evaluation of the mutations by in silico tools indicate their pathogenicity.

Conclusions

Our study is the first COL5A1 gene molecular investigation conducted among cEDS patients from Central Eastern Europe. Besides new COL5A1 variant findings, we gained molecular confirmation of clinical diagnosis of cEDS. In some cases, specific and adequate evaluation and classification of EDS patients based only on clinical features, may be difficult.

Splicing mutations in human genetic disorders: examples, detection, and confirmation

Precise pre-mRNA splicing, essential for appropriate protein translation, depends on the presence of consensus "cis" sequences that define exon-intron boundaries and regulatory sequences recognized by splicing machinery. Point mutations at these consensus sequences can cause improper exon and intron recognition and may result in the formation of an aberrant transcript of the mutated gene. The splicing mutation may occur in both introns and exons and disrupt existing splice sites or splicing regulatory sequences (intronic and exonic splicing silencers and enhancers), create new ones, or activate the cryptic ones. Usually such mutations result in errors during the splicing process and may lead to improper intron removal and thus cause alterations of the open reading frame. Recent research has underlined the abundance and importance of splicing mutations in the etiology of inherited diseases. The application of modern techniques allowed to identify synonymous and nonsynonymous variants as well as deep intronic mutations that affected pre-mRNA splicing. The bioinformatic algorithms can be applied as a tool to assess the possible effect of the identified changes. However, it should be underlined that the results of such tests are only predictive, and the exact effect of the specific mutation should be verified in functional studies. This article summarizes the current knowledge about the "splicing mutations" and methods that help to identify such changes in clinical diagnosis.

In-frame Variants in FLNA Proximal Rod 1 Domain Associate With a Predominant Cardiac Valvular Phenotype

INTRODUCTION AND OBJECTIVES

X-linked cardiac valvular dysplasia is a rare form of male-specific congenital heart defect mainly characterized by myxomatous degeneration of the atrioventricular valves with variable hemodynamic consequences. It is caused by genetic defects in FLNA-encoded filamin A, a widely expressed actin-binding protein that regulates cytoskeleton organization. Filamin A loss of function has also been associated with often concurring neurologic and connective tissue manifestations, with mutations in the first half of the Rod 1 domain apparently expressing the full cardiac phenotype. We contribute to previous genotype-phenotype correlations with a multidisciplinary approach in a newly-described family.

METHODS

Cardiologic, dysmorphologic, and genetic evaluation of available members were complemented with transcriptional and X-chromosome inactivation studies.

RESULTS

A novel FLNA mutation c.1066-3C>G cosegregated with a male-expressed, apparently isolated, cardiac phenotype with no skewed X-inactivation pattern in female carriers. This variant was shown to result in an in-frame deletion of 8 amino acid residues near the N-terminal region of the protein.

CONCLUSIONS

A nonimprinted, partial loss of function of filamin A proximal Rod 1 domain seems to be the pathogenetic mechanism of cardiac valvular dysplasia, with some cases occasionally expressing associated extracardiac manifestations.

Widespread intra-dependencies in the removal of introns from human transcripts

Research into the problem of splice site selection has followed a reductionist approach focused on how individual splice sites are recognized. Early applications of information theory uncovered an inconsistency. Human splice signals do not contain enough information to explain the observed fidelity of splicing. Here, we conclude that introns do not necessarily contain ‘missing’ information but rather may require definition from neighboring processing events. For example, there are known cases where an intronic mutation disrupts the splicing of not only the local intron but also adjacent introns. We present a genome-wide measurement of the order of splicing within human transcripts. The observed order of splicing cannot be explained by a simple kinetic model. Simulations reveal a bias toward a particular, transcript-specific order of intron removal in human genes. We validate an extreme class of intron that can only splice in a multi-intron context. Special categories of splicing such as exon circularization, first and last intron processing, alternative 5 and 3′ss usage and exon skipping are marked by distinct patterns of ordered intron removal. Excessive intronic length and silencer density tend to delay splicing. Shorter introns that contain enhancers splice early.

Vascular phenotypes in nonvascular subtypes of the Ehlers-Danlos syndrome: a systematic review

Purpose

Within the spectrum of the Ehlers-Danlos syndromes (EDS), vascular complications are usually associated with the vascular subtype of EDS. Vascular complications are also observed in other EDS subtypes, but the reports are anecdotal and the information is dispersed. To better document the nature of vascular complications among “nonvascular” EDS subtypes, we performed a systematic review.

Methods

We queried three databases for English-language studies from inception until May 2017, documenting both phenotypes and genotypes of patients with nonvascular EDS subtypes. The outcome included the number and nature of vascular complications.

Results

A total of 112 papers were included and data were collected from 467 patients, of whom 77 presented with a vascular phenotype. Severe complications included mainly hematomas (53%), frequently reported in musculocontractural and classical-like EDS; intracranial hemorrhages (18%), with a high risk in dermatosparaxis EDS; and arterial dissections (16%), frequently reported in kyphoscoliotic and classical EDS. Other, more minor, vascular complications were reported in cardiac-valvular, arthrochalasia, spondylodysplastic, and periodontal EDS.

Conclusion

Potentially life-threatening vascular complications are a rare but important finding in several nonvascular EDS subtypes, highlighting a need for more systematic documentation. This review will help familiarize clinicians with the spectrum of vascular complications in EDS and guide follow-up and management.

An ENU-induced splice site mutation of mouse Col1a1 causing recessive osteogenesis imperfecta and revealing a novel splicing rescue

GU-AG consensus sequences are used for intron recognition in the majority of cases of pre-mRNA splicing in eukaryotes. Mutations at splice junctions often cause exon skipping, short deletions, or insertions in the mature mRNA, underlying one common molecular mechanism of genetic diseases. Using N-ethyl-N-nitrosourea, a novel recessive mutation named seal was produced, associated with fragile bones and susceptibility to fractures (spine and limbs). A single nucleotide transversion (T → A) at the second position of intron 36 of the Col1a1 gene, encoding the type I collagen, α1 chain, was responsible for the phenotype. Col1a1 ^seal mRNA expression occurred at greatly reduced levels compared to the wild-type transcript, resulting in reduced and aberrant collagen fibers in tibiae of seal homozygous mice. Unexpectedly, splicing of Col1a1 ^seal mRNA followed the normal pattern despite the presence of the donor splice site mutation, likely due to the action of a putative intronic splicing enhancer present in intron 25, which appeared to function redundantly with the splice donor site of intron 36. Seal mice represent a model of human osteogenesis imperfecta, and reveal a previously unknown mechanism for splicing "rescue."

Inherited, not acquired, Gitelman syndrome in a patient with Sjögren's syndrome: importance of genetic testing to distinguish the two forms

Gitelman syndrome (GS) is an autosomal recessive, salt-losing renal tubulopathy caused by mutations in the SLC12A3 gene; however, it can also be acquired in patients with autoimmune disease, especially in those with Sjögren’s syndrome. Differentiating between the inherited and acquired forms of GS is clinically difficult. We report a case of inherited, not acquired, GS in a patient with Sjögren’s syndrome. A 41-year-old woman, who had been diagnosed with Sjögren’s syndrome at 27-years-old, had shown chronic hypokalemia (2.5–3.5 mmol/L). Laboratory tests showed hypokalemic alkalosis, hypomagnesemia, and hypocalciuria, corresponding to GS. Although acquired GS associated with Sjögren’s syndrome was initially suspected, a genetic test identified a novel homozygous mutation of c.1336-2A > T in the SLC12A3 gene, which resulted in aberrant splicing in the SLC12A3 transcript with the exclusion of exons 11 and 12. Thus, the GS was diagnosed as not the acquired but the inherited form. In the diagnosis of GS in patients with autoimmune disease, genetic testing of SLC12A3 is essential for differentiating the two forms.

CDHR1 mutations in retinal dystrophies

We report ophthalmic and genetic findings in patients with autosomal recessive retinitis pigmentosa (RP), cone-rod dystrophy (CRD) or cone dystrophy (CD) harboring potential pathogenic variants in the CDHR1 gene. Detailed ophthalmic examination was performed in seven sporadic and six familial subjects. Mutation screening was done using a customized next generation sequencing panel targeting 105 genes implicated in inherited retinal disorders. In one family, homozygosity mapping with subsequent candidate gene analysis was performed. Stringent filtering for rare and potentially disease causing variants following a model of autosomal recessive inheritance led to the identification of eleven different CDHR1 variants in nine index cases. All variants were novel at the time of their identification. In silico analyses confirmed their pathogenic potential. Minigene assays were performed for two non-canonical splice site variants and revealed missplicing for the mutant alleles. Mutations in CDHR1 are a rare cause of retinal dystrophy. Our study further expands the mutational spectrum of this gene and the associated clinical presentation.

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.

Human βA3/A1-crystallin splicing mutation causes cataracts by activating the unfolded protein response and inducing apoptosis in differentiating lens fiber cells

βγ-Crystallins, having a uniquely stable two domain four Greek key structure, are crucial for transparency of the eye lens,. Mutations in lens crystallins have been proposed to cause cataract formation by a variety of mechanisms most of which involve destabilization of the protein fold. The underlying molecular mechanism for autosomal dominant zonular cataracts with sutural opacities in an Indian family caused by a c.215+1G>A splice mutation in the βA3/A1-crystallin gene CRYBA1 was elucidated using three transgenic mice models. This mutation causes a splice defect in which the mutant mRNA escapes nonsense mediated decay by skipping both exons 3 and 4. Skipping these exons results in an in-frame deletion of the mRNA and synthesis of an unstable p.Ile33_Ala119del mutant βA3/A1-crystallin protein. Transgenic expression of mutant βA3/A1-crystallin but not the wild type protein results in toxicity and abnormalities in the maturation and orientation of differentiating lens fibers in c.97_357del CRYBA1 transgenic mice, leading to a small spherical lens, cataract, and often lens capsule rupture. On a cellular level, the lenses accumulated p.Ile33_Ala119del βA3/A1-crystallin with resultant activation of the stress signaling pathway - unfolded protein response (UPR) and inhibition of normal protein synthesis, culminating in apoptosis. This highlights the mechanistic contrast between mild mutations that destabilize crystallins and other proteins, resulting in their being bound by the α-crystallins that buffer lens cells against damage by denatured proteins, and severely misfolded proteins that are not bound by α-crystallin but accumulate and have a direct toxic effect on lens cells, resulting in early onset cataracts.

Specific variants in WDR35 cause a distinctive form of Ellis-van Creveld syndrome by disrupting the recruitment of the EvC complex and SMO into the cilium

Most patients with Ellis-van Creveld syndrome (EvC) are identified with pathogenic changes in EVC or EVC2, however further genetic heterogeneity has been suggested. In this report we describe pathogenic splicing variants in WDR35, encoding retrograde intraflagellar transport protein 121 (IFT121), in three families with a clinical diagnosis of EvC but having a distinctive phenotype. To understand why WDR35 variants result in EvC, we analysed EVC, EVC2 and Smoothened (SMO) in IFT-A deficient cells. We found that the three proteins failed to localize to Wdr35(-/-) cilia, but not to the cilium of the IFT retrograde motor mutant Dync2h1(-/-), indicating that IFT121 is specifically required for their entry into the ciliary compartment. Furthermore expression of Wdr35 disease cDNAs in Wdr35(-/-) fibroblasts revealed that the newly identified variants lead to Hedgehog signalling defects resembling those of Evc(-/-) and Evc2(-/-) mutants. Together our data indicate that splicing variants in WDR35, and possibly in other IFT-A components, underlie a number of EvC cases by disrupting targeting of both the EvC complex and SMO to cilia.

The effect of the common c.2299delG mutation in USH2A on RNA splicing

Recessive variants in the USH2A gene are an important cause of both Usher syndrome and nonsyndromic retinitis pigmentosa. A single base-pair deletion in exon 13 (c.2299delG, p.Glu767Serfs*21) is considered the most frequent mutation of USH2A. It is predicted to generate a premature termination codon and is presumed to lead to nonsense mediated decay. However the effect of this variant on RNA has not been formally investigated. It is not uncommon for exonic sequence alterations to cause aberrant splicing and the aim of the present report is to evaluate the effect of c.2299delG on USH2A transcripts. Nasal cells represent the simplest available tissue to study splicing defects in USH2A. Nasal brushing, RNA extraction from nasal epithelial cells and reverse transcription PCR were performed in five Usher syndrome patients who were homozygous for c.2299delG, two unaffected c.2299delG heterozygotes and seven control individuals. Primers to amplify between exons 12 and 15 and exons 10 and 14 were utilised. Significant variability was observed between different RT-PCR experiments. Importantly, in controls, PCR product of the expected size were amplified on all occasions (13/13 experiments); for patients this was true in only 4/14 experiments (Fisher exact test p = 0.0002). Bioinformatics tools predict the c.2299delG change to disrupt an exonic splicing enhancer and to create an exonic splicing silencer within exon 13. Here, we report an effect of the common c.2299delG mutation on splicing of exons 12 and 13 of USH2A. Future studies are expected to provide important insights into the contribution of this effect on the phenotype.

Validation of predicted mRNA splicing mutations using high-throughput transcriptome data

Interpretation of variants present in complete genomes or exomes reveals numerous sequence changes, only a fraction of which are likely to be pathogenic. Mutations have been traditionally inferred from allele frequencies and inheritance patterns in such data. Variants predicted to alter mRNA splicing can be validated by manual inspection of transcriptome sequencing data, however this approach is intractable for large datasets. These abnormal mRNA splicing patterns are characterized by reads demonstrating either exon skipping, cryptic splice site use, and high levels of intron inclusion, or combinations of these properties. We present, Veridical, an in silico method for the automatic validation of DNA sequencing variants that alter mRNA splicing. Veridical performs statistically valid comparisons of the normalized read counts of abnormal RNA species in mutant versus non-mutant tissues. This leverages large numbers of control samples to corroborate the consequences of predicted splicing variants in complete genomes and exomes.

Validation of predicted mRNA splicing mutations using high-throughput transcriptome data

Interpretation of variants present in complete genomes or exomes reveals numerous sequence changes, only a fraction of which are likely to be pathogenic. Mutations have been traditionally inferred from allele frequencies and inheritance patterns in such data. Variants predicted to alter mRNA splicing can be validated by manual inspection of transcriptome sequencing data, however this approach is intractable for large datasets. These abnormal mRNA splicing patterns are characterized by reads demonstrating either exon skipping, cryptic splice site use, and high levels of intron inclusion, or combinations of these properties. We present, Veridical, an in silico method for the automatic validation of DNA sequencing variants that alter mRNA splicing. Veridical performs statistically valid comparisons of the normalized read counts of abnormal RNA species in mutant versus non-mutant tissues. This leverages large numbers of control samples to corroborate the consequences of predicted splicing variants in complete genomes and exomes.

Two-exon skipping within MLPH is associated with coat color dilution in rabbits

Coat color dilution turns black coat color to blue and red color to cream and is a characteristic in many mammalian species. Matings among Netherland Dwarf, Loh, and Lionhead Dwarf rabbits over two generations gave evidence for a monogenic autosomal recessive inheritance of coat colour dilution. Histological analyses showed non-uniformly distributed, large, agglomerating melanin granules in the hair bulbs of coat color diluted rabbits. We sequenced the cDNA of MLPH in two dilute and one black rabbit for polymorphism detection. In both color diluted rabbits, skipping of exons 3 and 4 was present resulting in altered amino acids at p.QGL[37-39]QWA and a premature stop codon at p.K40*. Sequencing of genomic DNA revealed a c.111-5C>A splice acceptor mutation within the polypyrimidine tract of intron 2 within MLPH. This mutation presumably causes skipping of exons 3 and 4. In 14/15 dilute rabbits, the c.111-5C>A mutation was homozygous and in a further dilute rabbit, heterozygous and in combination with a homozygous frame shift mutation within exon 6 (c.585delG). In conclusion, our results demonstrated a colour dilution associated MLPH splice variant causing a strongly truncated protein (p.Q37QfsX4). An involvement of further MLPH-associated mutations needs further investigations.

Clinical and molecular characterization of 40 patients with classic Ehlers-Danlos syndrome: identification of 18 COL5A1 and 2 COL5A2 novel mutations

BACKGROUND

Classic Ehlers-Danlos syndrome (cEDS) is a rare autosomal dominant connective tissue disorder that is primarily characterized by skin hyperextensibility, abnormal wound healing/atrophic scars, and joint hypermobility. A recent study demonstrated that more than 90% of patients who satisfy all of these major criteria harbor a type V collagen (COLLV) defect.

METHODS

This cohort included 40 patients with cEDS who were clinically diagnosed according to the Villefranche nosology. The flowchart that was adopted for mutation detection consisted of sequencing the COL5A1 gene and, if no mutation was detected, COL5A2 analysis. In the negative patients the presence of large genomic rearrangements in COL5A1 was investigated using MLPA, and positive results were confirmed via SNP-array analysis.

RESULTS

We report the clinical and molecular characterization of 40 patients from 28 families, consisting of 14 pediatric patients and 26 adults. A family history of cEDS was present in 9 patients. The majority of the patients fulfilled all the major diagnostic criteria for cEDS; atrophic scars were absent in 2 females, skin hyperextensibility was not detected in a male and joint hypermobility was negative in 8 patients (20% of the entire cohort). Wide inter- and intra-familial phenotypic heterogeneity was observed. We identified causal mutations with a detection rate of approximately 93%. In 25/28 probands, COL5A1 or COL5A2 mutations were detected. Twenty-one mutations were in the COL5A1 gene, 18 of which were novel (2 recurrent). Of these, 16 mutations led to nonsense-mediated mRNA decay (NMD) and to COLLV haploinsufficiency and 5 mutations were structural. Two novel COL5A2 splice mutations were detected in patients with the most severe phenotypes. The known p. (Arg312Cys) mutation in the COL1A1 gene was identified in one patient with vascular-like cEDS.

CONCLUSIONS

Our findings highlight that the three major criteria for cEDS are useful and sufficient for cEDS clinical diagnosis in the large majority of the patients. The borderline patients for whom these criteria fail can be diagnosed when minor signs of connective tissue diseases and family history are present and when genetic testing reveals a defect in COLLV. Our data also confirm that COL5A1 and COL5A2 are the major, if not the only, genes involved in cEDS.

Molecular basis of two-exon skipping (exons 12 and 13) by c.1248+5g>a in OXCT1 gene: study on intermediates of OXCT1 transcripts in fibroblasts

The molecular basis of simultaneous two-exon skipping induced by a splice-site mutation has yet to be completely explained. The splice donor site mutation c.1248+5g>a (IVS13) of the OXCT1 gene resulted predominantly in skipping of exons 12 and 13 in fibroblasts from a patient (GS23) with succinyl-CoA:3-ketoacid CoA transferase (SCOT) deficiency. We compared heteronuclear RNA (hnRNA) intermediates between controls' and GS23's fibroblasts. Our strategy was to use RT-PCR of hnRNA to detect the presence or absence of spliced exon clusters in RNA intermediates (SECRIs) comprising sequential exons. Our initial hypothesis was that a SECRI comprising exons 12 and 13 was formed first followed by skipping of this SECRI in GS23 cells. However, such a pathway was revealed to be not a major one. Hence, we compared the intron removal of SCOT transcript between controls and GS23. In controls, intron 11 was the last intron to be spliced and the removal of intron 12 was also rather slow and occurred after the removal of intron 13 in a major pathway. However, the mutation in GS23 cells resulted in retention of intron 13, thus causing the retention of introns 12 and 11. This "splicing paralysis" may be solved by skipping the whole intron 11-exon 12-intron 12-exon 13-mutated intron 13, resulting in skipping of exons 12 and 13.