SplicePie: a novel analytical approach for the detection of alternative, non-sequential and recursive splicing.

Alternative splicing is a powerful mechanism present in eukaryotic cells to obtain a wide range of transcripts and protein isoforms from a relatively small number of genes. The mechanisms regulating (alternative) splicing and the paradigm of consecutive splicing have recently been challenged, especially for genes with a large number of introns. RNA-Seq, a powerful technology using deep sequencing in order to determine transcript structure and expression levels, is usually performed on mature mRNA, therefore not allowing detailed analysis of splicing progression. Sequencing pre-mRNA at different stages of splicing potentially provides insight into mRNA maturation. Although the number of tools that analyze total and cytoplasmic RNA in order to elucidate the transcriptome composition is rapidly growing, there are no tools specifically designed for the analysis of nuclear RNA (which contains mixtures of pre- and mature mRNA). We developed dedicated algorithms to investigate the splicing process. In this paper, we present a new classification of RNA-Seq reads based on three major stages of splicing: pre-, intermediate- and post-splicing. Applying this novel classification we demonstrate the possibility to analyze the order of splicing. Furthermore, we uncover the potential to investigate the multi-step nature of splicing, assessing various types of recursive splicing events. We provide the data that gives biological insight into the order of splicing, show that non-sequential splicing of certain introns is reproducible and coinciding in multiple cell lines. We validated our observations with independent experimental technologies and showed the reliability of our method. The pipeline, named SplicePie, is freely available at: https://github.com/pulyakhina/splicing_analysis_pipeline. The example data can be found at: https://barmsijs.lumc.nl/HG/irina/example_data.tar.gz.

Non-sequential and multi-step splicing of the dystrophin transcript.

The dystrophin protein encoding DMD gene is the longest human gene. The 2.2 Mb long human dystrophin transcript takes 16 hours to be transcribed and is co-transcriptionally spliced. It contains long introns (24 over 10kb long, 5 over 100kb long) and the heterogeneity in intron size makes it an ideal transcript to study different aspects of the human splicing process. Splicing is a complex process and much is unknown regarding the splicing of long introns in human genes. Here, we used ultra-deep transcript sequencing to characterize splicing of the dystrophin transcripts in 3 different human skeletal muscle cell lines, and explored the order of intron removal and multi-step splicing. Coverage and read pair analyses showed that around 40% of the introns were not always removed sequentially. Additionally, for the first time, we report that non-consecutive intron removal resulted in 3 or more joined exons which are flanked by unspliced introns and we defined these joined exons as an exon block. Lastly, computational and experimental data revealed that, for the majority of dystrophin introns, multistep splicing events are used to splice out a single intron. Overall, our data show for the first time in a human transcript, that multi-step intron removal is a general feature of mRNA splicing.

High throughput AS LNA qPCR method for the detection of a specific mutation in poliovirus vaccine strains.

Sabin Inactivated Poliovirus Vaccine (sIPV) has become one of the preferred vaccination options for the last step in the Poliovirus eradication program. Sequencing of poliovirus samples is needed during the manufacturing of poliovirus vaccines to assure the safety and immunogenicity of these vaccines. Next-generation sequencing analysis is the current costly and time-consuming gold standard for monitoring the manufacturing processes. We developed a low-cost and quick, highly sensitive, and allele-specific locked nucleic acid-probe-based reverse transcription quantitative PCR alternative that can accurately detect mutations in poliovirus vaccine samples during process development, scaling up, and release. Using the frequently in vitro occurring and viral replication-impacting VP1-E295K mutation as a showcase, we show that this technology can accurately detect E295K mutations in poliovirus 2 samples to similar levels as NGS. The qPCR technology was developed employing a synthetic dsDNA fragment-based standard curve containing mixes of E295K-WT (wildtype) and Mut (mutant) synthetic dsDNA fragments ranging from 1 × 107 copies/µL to 1 × 102 copies/µL to achieve a linear correlation with R2 > 0.999, and PCR efficiencies of 95-105 %. Individual standard concentration levels achieved accuracies of ≥92 % (average 96 %) and precisions of ≤17 % (average 3.3 %) RSD. Specificity of locked nucleic acid (LNA)-probes was confirmed in the presence and absence of co-mutations in the probe-binding region. Application of the developed assay to Sabin Poliovirus type 2 production run samples, illustrated a linear relationship with an R2 of 0.994, and an average accuracy of 97.2 % of the variant (allele)-specific AS LNA qPCR result, compared to NGS. The assay showed good sensitivity for poliovirus samples, containing E295K mutation levels between 0 % and 95 % (quantification range). In conclusion, the developed AS LNA qPCR presents a valuable low-cost, and fast tool, suitable for the process development and quality control of polio vaccines.

Antisense-Mediated Skipping of Dysferlin Exons in Control and Dysferlinopathy Patient-Derived Cells.

Dysferlinopathies encompass a spectrum of progressive muscular dystrophies caused by the lack of dysferlin due to missense mutations in the dysferlin gene or mutations causing premature truncation of protein translation. Dysferlin is a modular protein, and dysferlins lacking one or more repetitive domains have been shown to retain functionality. As such, antisense-mediated exon skipping has been proposed as a therapy for dysferlinopathy. By skipping the mutated exon, the reading frame would be maintained, while the mutation would be bypassed, thus allowing production of an internally deleted, but partially functional, dysferlin. We previously showed that dysferlin exon skipping is feasible in control cell lines. We here evaluated exon skipping and dysferlin protein restoration in patient-derived cells requiring the skipping of exon 9, 29, 30, or 34. Exon 30 skipping was possible at high efficiency, but did not result in increased dysferlin. We discovered that the alleged exon 30 mutation was in fact a polymorphism and identified a splicing mutation in intron 28 as the disease-causing mutation. While exon skipping was feasible for each of the other cell lines, no increases in dysferlin protein could be detected by western blotting.

Prognostic relevance of MLH1 and MSH2 mutations in hereditary non-polyposis colorectal cancer patients.

AIMS AND BACKGROUND: Colorectal carcinoma patients from hereditary non-polyposis colorectal cancer families are suggested to have a better prognosis than sporadic colorectal carcinoma cases. Since the majority of hereditary non-polyposis colorectal cancer-related colorectal carcinomas are characterized by microsatellite instability due to germline mutations in DNA mismatch repair genes, this is consistent with the prolonged survival observed in sporadic microsatellite instability-positive colorectal carcinoma compared to microsatellite stable cases. However, a fraction of colorectal carcinoma cases belongs to families that, despite fulfilling the clinical criteria for hereditary non-polyposis colorectal cancer, do not carry mismatch repair gene mutations. Our aim was to verify to what extent the genotypic heterogeneity influences the prognosis of hereditary non-polyposis colorectal cancer patients. METHODS: A survival analysis was performed on 526 colorectal carcinoma cases from 204 Amsterdam Criteria-positive hereditary non-polyposis colorectal cancer families. Enrolled cases were classified as MLH1-positive, MSH2-positive and mutation-negative, according to the results of genetic testing in each family. RESULTS: Five-year survival rates were 0.73 (95% CI, 0.66-0.80), 0.75 (95% CI, 0.66-0.84) and 0.62 (95% CI, 0.55-0.68) for MLH1-positive, MSH2-positive and mutation-negative groups, respectively (logrank test, P = 0.01). Hazard ratio, computed using Cox regression analysis and adjusted for age, sex, tumor site and stage, was 0.71 (95% CI, 0.51-0.98) for the mutation-positive compared to the mutation-negative group. Moreover, in the latter group, patients with microsatellite instability-positive colorectal carcinomas showed a better outcome than microsatellite stable cases (5-year survival rates, 0.81 and 0.60, respectively; logrank test, P = 0.006). CONCLUSIONS: Our results suggest that the prognosis of hereditary non-polyposis colorectal cancer-related colorectal carcinoma patients depends on the associated constitutional mismatch repair genotype.

Regulation of the human MSH6 gene by the Sp1 transcription factor and alteration of promoter activity and expression by polymorphisms.

Defects in human DNA mismatch repair have been reported to underlie a variety of hereditary and sporadic cancer cases. We characterized the structure of the MSH6 promoter region to examine the mechanisms of transcriptional regulation of the MSH6 gene. The 5'-flanking region of the MSH6 gene was found to contain seven functional Sp1 transcription factor binding sites that each bind Sp1 and Sp3 and contribute to promoter activity. Transcription did not appear to require a TATA box and resulted in multiple start sites, including two major start sites and at least nine minor start sites. Three common polymorphisms were identified in the promoter region (-557 T-->G, -448 G-->A, and -159 C-->T): the latter two were always associated, and each of these functionally inactivated a different Sp1 site. The polymorphic allele -448 A -159 T was demonstrated to be a common Caucasian polymorphism found in 16% of Caucasians and resulted in a five-Sp1-site promoter that had 50% less promoter activity and was more sensitive to inactivation by DNA methylation than the more common seven Sp1 site promoter allele, which was only partially inactivated by DNA methylation. In cell lines, this five-Sp1-site polymorphism resulted in reduced MSH6 expression at both the mRNA and protein level. An additional 2% of Caucasians contained another polymorphism, -210 C-->T, which inactivated a single Sp1 site that also contributes to promoter activity.

A hereditary nonpolyposis colorectal carcinoma case associated with hypermethylation of the MLH1 gene in normal tissue and loss of heterozygosity of the unmethylated allele in the resulting microsatellite instability-high tumor.

Fourteen suspected hereditary nonpolyposis colorectal carcinoma cases with microsatellite unstable(microsatellite instability-high; MSI-H) tumors but no germ-line MSH2, MSH6, or MLH1 mutations were examined for hypermethylation of CpG sites in the critical promoter region of MLH1. The methylation patterns were determined using methylation-specific PCR and by sequence analysis of sodium bisulfite-treated genomic DNA. In one case, DNA hypermethylation of one allele was detected in DNA isolated from blood. In the MSI-H tumor from this case, the unmethylated MLH1 allele was eliminated by loss of heterozygosity, and the methylated allele was retained. This biallelic inactivation resulted in loss of expression of MLH1 in the tumor as confirmed by immunohistochemistry. These results suggest a novel mode of germ-line inactivation of a cancer susceptibility gene.

Absence of α- and ß-dystroglycan is associated with Walker-Warburg syndrome.

OBJECTIVE: To identify the underlying genetic defect in 5 patients from a consanguineous family with a Walker-Warburg phenotype, together with intracranial calcifications. METHODS: Homozygosity mapping and exome sequencing, followed by Sanger sequencing of the obtained candidate gene, was performed. Expression of the candidate gene was tested by reverse transcription PCR. Patient fibroblasts were converted to myotubes, and the expression and function of dystroglycan was tested by Western blotting. RESULTS: We detected a homozygous loss-of-function frameshift mutation in the DAG1 gene and showed that this mutation results in a complete absence of both α- and ß-dystroglycan. CONCLUSIONS: A loss-of-function mutation in DAG1 can result in Walker-Warburg syndrome and is not embryonic lethal.

Comprehensive molecular analysis of mismatch repair gene defects in suspected Lynch syndrome (hereditary nonpolyposis colorectal cancer) cases.

An accurate algorithm is essential for effective molecular diagnosis of hereditary colorectal cancer (CRC). Here, we have extended the analysis of 71 CRC cases suspected to be Lynch syndrome cases for MSH2, MLH1, MSH6, and PMS2 gene defects. All cases were screened for mutations in MSH2, MLH1, and MSH6, and all cases where tumors were available were screened for microsatellite instability (MSI) and expression of MSH2 and MLH1. Subsequently, mutation-negative cases were screened for MLH1 methylation and mutations in PMS2. Of the MSI-high (MSI-H) cases, 96% had a mismatch repair (MMR) gene defect, mostly involving MSH2 or MLH1; one PMS2 mutation, one MLH1 epimutation, and no MSH6 mutations were found. Four of the 28 MSI-H cases, including one Amsterdam criteria case, had biallelic tumor MLH1 methylation, indicating that sporadic cases can be admixed in with Lynch syndrome cases, even those meeting the strongest criteria for Lynch syndrome. MMR gene defects were found in similar frequency in cases where tumors were and were not available. One MLH1 and one MSH2 deletion mutation were found in MSI-stable/low cases, indicating that MSI testing can exclude cases with pathogenic mutations. Our analysis supports a diagnostic algorithm where cases are selected for analysis based on clinical criteria or prediction models; isolated sporadic young-onset cases can be prescreened by tumor testing, whereas familial cases may be directly subjected to molecular analysis for mutations in MMR genes followed by MSI, protein expression, and DNA methylation analysis to aid in the resolution of mutation-negative cases.