Analysis of splicing patterns by pyrosequencing

Capitalizing on the heterogeneous effects of CFTR nonsense and frameshift variants to inform therapeutic strategy for cystic fibrosis

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic 'nulls' as some may allow generation of protein that can be targeted to achieve clinical benefit.

Experimental assessment of splicing variants using expression minigenes and comparison with in silico predictions

Assessment of the functional consequences of variants near splice sites is a major challenge in the diagnostic laboratory. To address this issue, we created expression minigenes (EMGs) to determine the RNA and protein products generated by splice site variants (n = 10) implicated in cystic fibrosis (CF). Experimental results were compared with the splicing predictions of eight in silico tools. EMGs containing the full-length Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) coding sequence and flanking intron sequences generated wild-type transcript and fully processed protein in Human Embryonic Kidney (HEK293) and CF bronchial epithelial (CFBE41o-) cells. Quantification of variant induced aberrant mRNA isoforms was concordant using fragment analysis and pyrosequencing. The splicing patterns of c.1585-1G>A and c.2657+5G>A were comparable to those reported in primary cells from individuals bearing these variants. Bioinformatics predictions were consistent with experimental results for 9/10 variants (MES), 8/10 variants (NNSplice), and 7/10 variants (SSAT and Sroogle). Programs that estimate the consequences of mis-splicing predicted 11/16 (HSF and ASSEDA) and 10/16 (Fsplice and SplicePort) experimentally observed mRNA isoforms. EMGs provide a robust experimental approach for clinical interpretation of splice site variants and refinement of in silico tools.

Pyrosequencing: powerful and quantitative sequencing technology

Pyrosequencing is a sequencing-by-synthesis method for DNA analysis that has emerged as a platform not only for de novo sequencing applications, but also for quantitative analysis of genomic methylation, single-nucleotide polymorphisms, and allele quantification. In this unit, we describe a complete workflow from sample to result that is suitable for each of these applications. As cytosine conversion is a key element of successful methylation analysis using pyrosequencing, a support protocol for bisulfite treatment is also included.

Quantitative evaluation of all hexamers as exonic splicing elements

We describe a comprehensive quantitative measure of the splicing impact of a complete set of RNA 6-mer sequences by deep sequencing successfully spliced transcripts. All 4096 6-mers were substituted at five positions within two different internal exons in a 3-exon minigene, and millions of successfully spliced transcripts were sequenced after transfection of human cells. The results allowed the assignment of a relative splicing strength score to each mutant molecule. The effect of 6-mers on splicing often depended on their location; much of this context effect could be ascribed to the creation of different overlapping sequences at each site. Taking these overlaps into account, the splicing effect of each 6-mer could be quantified, and 6-mers could be designated as enhancers (ESEseqs) and silencers (ESSseqs), with an ESRseq score indicating their strength. Some 6-mers exhibited positional bias relative to the two splice sites. The distribution and conservation of these ESRseqs in and around human exons supported their classification. Predicted RNA secondary structure effects were also seen: Effective enhancers, silencers and 3' splice sites tend to be single stranded, and effective 5' splice sites tend to be double stranded. 6-mers that may form positive or negative synergy with another were also identified. Chromatin structure may also influence the splicing enhancement observed, as a good correspondence was found between splicing performance and the predicted nucleosome occupancy scores of 6-mers. This approach may prove of general use in defining nucleic acid regulatory motifs, substitute for functional SELEX in most cases, and provide insights about splicing mechanisms.

Reverse genetics in eukaryotes

Reverse genetics consists in the modification of the activity of a target gene to analyse the phenotypic consequences. Four main approaches are used towards this goal and will be explained in this review. Two of them are centred on genome alterations. Mutations produced by random chemical or insertional mutagenesis can be screened to recover only mutants in a specific gene of interest. Alternatively, these alterations may be specifically targeted on a gene of interest by HR (homologous recombination). The other two approaches are centred on mRNA. RNA interference is a powerful method to reduce the level of gene products, while MO (morpholino) antisense oligonucleotides alter mRNA metabolism or translation. Some model species, such as Drosophila, are amenable to most of these approaches, whereas other model species are restricted to one of them. For example, in mice and yeasts, gene targeting by HR is prevalent, whereas in Xenopus and zebrafish MO oligonucleotides are mainly used. Genome-wide collections of mutants or inactivated models obtained in several species by these approaches have been made and will help decipher gene functions in the post-genomic era.

MapSplice: accurate mapping of RNA-seq reads for splice junction discovery

The accurate mapping of reads that span splice junctions is a critical component of all analytic techniques that work with RNA-seq data. We introduce a second generation splice detection algorithm, MapSplice, whose focus is high sensitivity and specificity in the detection of splices as well as CPU and memory efficiency. MapSplice can be applied to both short (<75 bp) and long reads (≥75 bp). MapSplice is not dependent on splice site features or intron length, consequently it can detect novel canonical as well as non-canonical splices. MapSplice leverages the quality and diversity of read alignments of a given splice to increase accuracy. We demonstrate that MapSplice achieves higher sensitivity and specificity than TopHat and SpliceMap on a set of simulated RNA-seq data. Experimental studies also support the accuracy of the algorithm. Splice junctions derived from eight breast cancer RNA-seq datasets recapitulated the extensiveness of alternative splicing on a global level as well as the differences between molecular subtypes of breast cancer. These combined results indicate that MapSplice is a highly accurate algorithm for the alignment of RNA-seq reads to splice junctions. Software download URL: http://www.netlab.uky.edu/p/bioinfo/MapSplice.

Integration of molecular functions at the ecosystemic level: breakthroughs and future goals of environmental genomics and post-genomics

Environmental genomics and genome-wide expression approaches deal with large-scale sequence-based information obtained from environmental samples, at organismal, population or community levels. To date, environmental genomics, transcriptomics and proteomics are arguably the most powerful approaches to discover completely novel ecological functions and to link organismal capabilities, organism-environment interactions, functional diversity, ecosystem processes, evolution and Earth history. Thus, environmental genomics is not merely a toolbox of new technologies but also a source of novel ecological concepts and hypotheses. By removing previous dichotomies between ecophysiology, population ecology, community ecology and ecosystem functioning, environmental genomics enables the integration of sequence-based information into higher ecological and evolutionary levels. However, environmental genomics, along with transcriptomics and proteomics, must involve pluridisciplinary research, such as new developments in bioinformatics, in order to integrate high-throughput molecular biology techniques into ecology. In this review, the validity of environmental genomics and post-genomics for studying ecosystem functioning is discussed in terms of major advances and expectations, as well as in terms of potential hurdles and limitations. Novel avenues for improving the use of these approaches to test theory-driven ecological hypotheses are also explored.

ERK-regulated differential expression of the Mitf 6a/b splicing isoforms in melanoma

The master regulator of the melanocyte lineage Mitf is intimately involved in development as well as melanoma, controlling cell survival, differentiation, proliferation and metastasis/migration. Consistent with its central role, Mitf expression and Mitf post-translational modifications are tightly regulated. An additional potential level of regulation is afforded by differential splicing of Mitf exon-6 leading to the generation of two isoforms that differ by the presence of six amino-acids in the Mitf (+) isoform and which have differential effects on cell cycle progression. However, whether the ratio of the two isoforms is regulated and whether their expression correlates with melanoma progression is not known. Here, we show that the differential expression of the Mitf 6a/b isoforms is dependent on the MAPKinase signalling, being linked to the activation of MEK1-ERK2, but not to N-RAS/B-RAF mutation status. In addition, quantification of Mitf 6a/b splicing forms in 86 melanoma samples revealed substantially increased levels of the Mitf (-) form in a subset of metastatic melanomas. The results suggest that differential expression of the Mitf 6a/b isoforms may represent an additional mechanism for regulating Mitf function and melanoma biology.