A new pseudoexon activation due to ultrarare branch point formation in Duchenne muscular dystrophy

Deep-intronic variants that create or enhance a splice site are increasingly reported as a significant cause of monogenic diseases. However, deep-intronic variants that activate pseudoexons by affecting a branch point are extremely rare in monogenic diseases. Here, we describe a novel deep-intronic DMD variant that created a branch point in a Duchenne muscular dystrophy (DMD) patient. A 7.0-year-old boy was enrolled because he was suspected of DMD based on his clinical, muscle imaging, and pathological features. Routine genetic testing did not discover a pathogenic DMD variant. We then performed muscle-derived dystrophin mRNA analysis and detected an aberrant pseudoexon-containing transcript. Further genomic Sanger sequencing and bioinformatic analyses revealed a novel deep-intronic splicing variant in DMD (NM_004006.2:c.5325+1759G>T), which created a new branch point sequence and thus activated a new dystrophin pseudoexon (NM_004006.2:r.5325_5326ins5325+1779_5325+1855). Our study highlights the significant role of branch point alterations in the pathogenesis of monogenic diseases.

From the comparative study of a circRNA originating from an mammalian ATXN2L intron to understanding the genesis of intron lariat-derived circRNAs

Circular intronic RNAs (ciRNAs) are still unexplored regarding mechanisms for their emergence. We considered the ATXN2L intron lariat-derived circular RNA (ciRNA-ATXN2L) as an opportunity to conduct a cross-species examination of ciRNA genesis. To this end, we investigated 207 datasets from 4 tissues and from 13 mammalian species. While in eight species, ciRNA-ATXN2L was never detected, in pigs and rabbits, ciRNA-ATXN2L was expressed in all tissues and sometimes at very high levels. Bovine tissues were an intermediate case and in macaques and cats, only ciRNA-ATXN2L traces were detected. The pattern of ciRNA-ATXN2L restricted to only five species is not related to a particular evolution of intronic sequences. To empower our analysis, we considered 221 additional introns including 80 introns where a lariat-derived ciRNA was previously described. The primary driver of micro-ciRNA genesis (< 155 nt as ciRNA-ATXN2L) appears to be the absence of a canonical "A" (i.e. a "tnA" located in the usual branching region) to build the lariat around this adenosine. The balance between available "non canonical-A" (no ciRNA genesis) and "non-A" (ciRNA genesis) for use as a branch point to build the lariat could modify the expression level of ciRNA-ATXN2L. In addition, the rare localization of the 2'-5' bond in an open RNA secondary structure could also negatively affect the lifetime of ciRNAs (macaque ciRNA-ATXN2L). Our analyses suggest that ciRNA-ATXN2L is likely a functionless splice remnant. This study provides a better understanding of the ciRNAs origin, especially drivers for micro ciRNA genesis.

Identification of essential intron sequences that enhance gene expression independently of splicing in the yeast Saccharomyces cerevisiae

Gene expression in eukaryotes is enhanced by the presence of introns in a process known as intron-mediated enhancement (IME), but its mechanism remains unclear. In Saccharomyces cerevisiae, sequences at the 5'-splice sites (SS) and branch point sites (BPS) are highly conserved compared with other higher eukaryotes. Here, the minimum intron sequence essential for IME was investigated using various short introns and a yeast codon-optimized luciferase gene as an IME model. Mutations at the 5'-SS conserved sequence and branch point in the QCR10 intron caused splicing deficiency with either a complete loss or a marked decrease in IME. By contrast, however, the 3'-AG to tG mutant was spliced and retained IME function. Moreover, heterologous introns, which did not show IME in S. cerevisiae, gained splicing competency and IME ability by substitutions to the S. cerevisiae-type 5'-SS and BPS sequences. Intriguingly, several deletion mutants between the 5'-SS and BPS in introns exhibited high levels of IME despite a loss in splicing competency. In most cases, further deletions or substitutions did not recover splicing competency and were found to decrease IME. However, a 16-nt variant consisting of the conserved 5'-SS and BPS sequences and 3'-CAG showed an IME level comparable with that of the wild-type intron. These results indicate that IME can be independent of splicing in S. cerevisiae while intron sequences at the 5'-SS and BPS play an essential role in IME.

Assessment of branch point prediction tools to predict physiological branch points and their alteration by variants

BACKGROUND

Branch points (BPs) map within short motifs upstream of acceptor splice sites (3'ss) and are essential for splicing of pre-mature mRNA. Several BP-dedicated bioinformatics tools, including HSF, SVM-BPfinder, BPP, Branchpointer, LaBranchoR and RNABPS were developed during the last decade. Here, we evaluated their capability to detect the position of BPs, and also to predict the impact on splicing of variants occurring upstream of 3'ss.

RESULTS

We used a large set of constitutive and alternative human 3'ss collected from Ensembl (n = 264,787 3'ss) and from in-house RNAseq experiments (n = 51,986 3'ss). We also gathered an unprecedented collection of functional splicing data for 120 variants (62 unpublished) occurring in BP areas of disease-causing genes. Branchpointer showed the best performance to detect the relevant BPs upstream of constitutive and alternative 3'ss (99.48 and 65.84% accuracies, respectively). For variants occurring in a BP area, BPP emerged as having the best performance to predict effects on mRNA splicing, with an accuracy of 89.17%.

CONCLUSIONS

Our investigations revealed that Branchpointer was optimal to detect BPs upstream of 3'ss, and that BPP was most relevant to predict splicing alteration due to variants in the BP area.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene

BACKGROUND

The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates.

RESULTS

The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination.

CONCLUSIONS

We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

HPV18 Utilizes Two Alternative Branch Sites for E6*I Splicing to Produce E7 Protein

Human papillomavirus 18 (HPV18) E6 and E7 oncogenes are transcribed as a single bicistronic E6E7 pre-mRNA. The E6 ORF region in the bicistronic E6E7 pre-mRNA contains an intron. Splicing of this intron disrupts the E6 ORF integrity and produces a spliced E6*I RNA for efficient E7 translation. Here we report that the E6 intron has two overlapped branch point sequences (BPS) upstream of its 3′ splice site, with an identical heptamer AACUAAC, for E6*I splicing. One heptamer has a branch site adenosine (underlined) at nt 384 and the other at nt 388. E6*I splicing efficiency correlates to the expression level of E6 and E7 proteins and depends on the selection of which branch site. In general, E6*I splicing prefers the 3′ss-proximal branch site at nt 388 over the distal branch site at nt 384. Inactivation of the nt 388 branch site was found to activate a cryptic acceptor site at nt 636 for aberrant RNA splicing. Together, these data suggest that HPV18 modulates its production ratio of E6 and E7 proteins by alternative selection of the two mapped branch sites for the E6*I splicing, which could be beneficial in its productive or oncogenic infection according to the host cell environment.

A purine-rich element in foamy virus pol regulates env splicing and gag/pol expression

Background

The foamy viral genome encodes four central purine-rich elements localized in the integrase-coding region of pol. Previously, we have shown that the first two of these RNA elements (A and B) are required for protease dimerization and activation. The D element functions as internal polypurine tract during reverse transcription. Peters et al., described the third element (C) as essential for gag expression suggesting that it might serve as an RNA export element for the unspliced genomic transcript.

Results

Here, we analysed env splicing and demonstrate that the described C element composed of three GAA repeats known to bind SR proteins regulates env splicing, thus balancing the amount of gag/pol mRNAs. Deletion of the C element effectively promotes a splice site switch from a newly identified env splice acceptor to the intrinsically strong downstream localised env 3′ splice acceptor permitting complete splicing of almost all LTR derived transcripts. We provide evidence that repression of this env splice acceptor is a prerequisite for gag expression. This repression is achieved by the C element, resulting in impaired branch point recognition and SF1/mBBP binding. Separating the branch point from the overlapping purine-rich C element, by insertion of only 20 nucleotides, liberated repression and fully restored splicing to the intrinsically strong env 3′ splice site. This indicated that the cis-acting element might repress splicing by blocking the recognition of essential splice site signals.

Conclusions

The foamy viral purine-rich C element regulates splicing by suppressing the branch point recognition of the strongest env splice acceptor. It is essential for the formation of unspliced gag and singly spliced pol transcripts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-017-0337-6) contains supplementary material, which is available to authorized users.