Splicing Signals Research Articles

A novel mutation in EYA1 in a Chinese family with Branchio-oto-renal syndrome

BackgroundBranchio-oto-renal (BOR) syndrome is a dominant autosomal disorder characterized by phenotypes such as hearing loss, branchial fistulae, preauricular pits, and renal abnormalities. EYA1, the human homolog of the Drosophila “eye absent” gene on chromosome 8q13.3, is recognized as one of the most important genes associated with BOR syndrome.MethodsThe proposita in this study was a 5-year-old Chinese girl with hearing loss, bilateral otitis media with effusion, microtia, facial hypoplasia, palatoschisis, and bilateral branchial cleft fistulae. The girl’s family members, except two who were deceased, agreed to undergo clinical examination. We collected blood samples from 10 family members, including six who were affected by the syndrome. Genomic DNA was extracted and subjected to Sanger sequencing. A minigene assay was performed to confirm whether splicing signals were altered. In addition, we performed western blotting to determine alterations in protein levels of the wild-type and mutant gene.ResultsClinical tests showed that some of the family members met the criteria for BOR syndrome. The affected members harbored a novel heterozygous nonsense variation in exon 11 of EYA1, whereas no unaffected member carried the mutation at this position. Functional experiments did not detect abnormal splicing at the RNA level; however, western blotting showed that the mutated protein was truncated.ConclusionsThis study reports a novel mutation associated with BOR syndrome in a Chinese family. We highlight the usefulness of genetic testing in the diagnosis of BOR syndrome. Thus, we believe that this report would benefit clinicians in this field.

The genetic and molecular pathogenesis of myelodysplastic syndromes.

Myelodysplastic syndromes (MDS) comprise a diverse group of clonal and malignant myeloid disorders characterized by ineffective hematopoiesis, resultant peripheral cytopenias, and a meaningful increased risk of progression to acute myeloid leukemia. A wide array of recurring genetic mutations involved in RNA splicing, histone manipulation, DNA methylation, transcription factors, kinase signaling, DNA repair, cohesin proteins, and other signal transduction elements has been identified as important substrates for the development of MDS. Cytogenetic abnormalities, namely those characterized by loss of genetic material (including 5q- and 7q-), have also been strongly implicated and may influence the clonal architecture which predicts such mutations and may provoke an inflammatory bone marrow microenvironment as the substrate for clonal expansion. Other aspects of the molecular pathogenesis of MDS continue to be further elucidated, predicated upon advances in gene expression profiling and the development of new, and improved high-throughput techniques. More accurate understanding of the genetic and molecular basis for the development of MDS directly provides additional opportunity for treatment, which to date remains limited. In this comprehensive review, we examine the current understanding of the molecular pathogenesis and pathophysiology of MDS, as well as review future prospects which may enhance this understanding, treatment strategies, and hopefully outcomes.

Designing an effective drug combination for ER stress loading in cancer therapy using a real-time monitoring system

Excess stress caused by accumulation of misfolded proteins inside the endoplasmic reticulum (ER) lumen can cause cells to undergo apoptosis. Misfolded proteins exported from ER to cytoplasm are ubiquitinated and mostly degraded by the proteasome, but can also be destroyed by autophagy mediated by the docking proteins p62 and NBR1. When misfolded proteins accumulate beyond the capacity of these clearance systems, they are transported to the microtubule organization center to form aggresomes, which are also degraded by autophagy. Together, these phenomena suggest the existence of a coordinated intracellular network for coping with the accumulation of misfolded proteins. Thus, rational inhibition of this network system might enhance killing of cancer cells subjected to pronounced ER stress loading. Based on this rationale, we sought to establish a quantitative assay for monitoring ER stress loading. MDA-MB231 cells stably transfected with the ERAI-Venus vector exhibited a strong XBP1 splicing signal in response to ER stress. Using the IncuCyte cell imaging system, we monitored the fluorescence intensity of XBP1-Venus, normalized against cell density, as an ER stress indicator. This parameter correlated closely with other reporters of unfolded protein responses. Assessment of the XBP1-Venus signal during exposure to various drug combinations revealed that simultaneous inhibition of the proteasome, autophagy, and aggresome formation led to more effective ER stress loading and higher cytotoxicity than inhibition of only two components. Our data suggest that this monitoring system is a useful tool for designing effective drug combinations for ER stress loading in cancer therapy.

In vivo and In vitro methods to identify DNA sequence variants that alter RNA Splicing.

Identification of sequence variants that create or eliminate splice sites has proven to be a significant challenge and represents one of many roadblocks in the clinical interpretation of rare genetic variation. Current methods of identifying splice altering sequence variants exist, however, these are limited by an imperfect understanding of splice signals and cumbersome functional assays. We have recently developed a computational tool that prioritizes putative splice-altering sequence variants, and a moderate-throughput minigene assay that confirms the variants which alter splicing. This bioinformatic strategy represents a substantial increase in accuracy and efficiency of historical in vitro splicing assays. In this unit we give detailed instructions on how to organize, run, and interpret various features of this protocol. We expect that splice-altering variants revealed through this protocol can be reliably carried forward for further clinical and biological analyses.

SpliceVec: Distributed feature representations for splice junction prediction

Identification of intron boundaries, called splice junctions, is an important part of delineating gene structure and functions. This also provides valuable insights into the role of alternative splicing in increasing functional diversity of genes. Identification of splice junctions through RNA-seq is by mapping short reads to the reference genome which is prone to errors due to random sequence matches. This encourages identification of splicing junctions through computational methods based on machine learning. Existing models are dependent on feature extraction and selection for capturing splicing signals lying in the vicinity of splice junctions. But such manually extracted features are not exhaustive. We introduce distributed feature representation, SpliceVec, to avoid explicit and biased feature extraction generally adopted for such tasks. SpliceVec is based on two widely used distributed representation models in natural language processing. Learned feature representation in form of SpliceVec is fed to multilayer perceptron for splice junction classification task. An intrinsic evaluation of SpliceVec indicates that it is able to group true and false sites distinctly. Our study on optimal context to be considered for feature extraction indicates inclusion of entire intronic sequence to be better than flanking upstream and downstream region around splice junctions. Further, SpliceVec is invariant to canonical and non-canonical splice junction detection. The proposed model is consistent in its performance even with reduced dataset and class-imbalanced dataset. SpliceVec is computationally efficient and can be trained with user-defined data as well.

Requirement for intron structures in activating the Cd8a locus

During differentiation of CD4+CD8+ double-positive (DP) thymocytes into the CD4-CD8+ single-positive (CD8SP) thymocytes committed to the cytotoxic T cell lineage, Cd8a transcription is temporally terminated after positive selection and is subsequently reinitiated, a process known as coreceptor reversal. Despite the identification of a transcriptional enhancer in the Cd8a gene that directs reporter transgene expression specifically in CD8SP thymocytes, the molecular mechanisms controlling reactivation of the Cd8a gene are not fully understood. Here, we show that, after positive selection, hCD2 reporter expression from the Cd8a locus, which was generated by insertion of hCD2 cDNA into the first exon of the Cd8a gene, requires the incorporation of intron sequences into the hCD2 transcript. The presence of polyadenylation signals after hCD2 cDNA inhibited hCD2 expression in mature CD8+ T cells, whereas hCD2 expression in DP thymocytes recapitulated the Cd8a expression. Incorporation of the endogenous short intron structure and heterologous intron structure of the Cd4 locus restored hCD2 expression in mature CD8+ T cells in a variegated manner. Interestingly, stage-specific DNA demethylation was impaired in Cd8a reporter alleles that failed to express hCD2 in CD8+ T cells, and intron sequences lacking RNA splicing signals still restored hCD2 expression. These observations indicate that "intron-mediated enhancement" is involved in a stage-specific reactivation of the Cd8a locus harboring hCD2 cDNA. However, the Cd8a gene was transcribed in mature CD8+ T cells, albeit at a lower level, from a mutant Cd8a locus lacking intron structures, suggesting that protein-coding sequences in transcripts affect sensitivity to intron-mediated enhancement.

Enhancing heterologous expression in Chlamydomonas reinhardtii by transcript sequence optimization.

Various species of microalgae have recently emerged as promising host-organisms for use in biotechnology industries due to their unique properties. These include efficient conversion of sunlight into organic compounds, the ability to grow in extreme conditions and the occurrence of numerous post-translational modification pathways. However, the inability to obtain high levels of nuclear heterologous gene expression in microalgae hinders the development of the entire field. To overcome this limitation, we analyzed different sequence optimization algorithms while studying the effect of transcript sequence features on heterologous expression in the model microalga Chlamydomonas reinhardtii, whose genome consists of rare features such as a high GC content. Based on the analysis of genomic data, we created eight unique sequences coding for a synthetic ferredoxin-hydrogenase enzyme, used here as a reporter gene. Following in silico design, these synthetic genes were transformed into the C. reinhardtii nucleus, after which gene expression levels were measured. The empirical data, measured invivo show a discrepancy of up to 65-fold between the different constructs. In this work we demonstrate how the combination of computational methods and our empirical results enable us to learn about the way gene expression is encoded in the C. reinhardtii transcripts. We describe the deleterious effect on overall expression of codons encoding for splicing signals. Subsequently, our analysis shows that utilization of a frequent subset of preferred codons results in elevated transcript levels, and that mRNA folding energy in the vicinity of translation initiation significantly affects gene expression.

Intron specificity in pre-mRNA splicing.

The occurrence of spliceosomal introns in eukaryotic genomes is highly diverse and ranges from few introns in an organism to multiple introns per gene. Introns vary with respect to their lengths, strengths of splicing signals, and position in resident genes. Higher intronic density and diversity in genetically complex organisms relies on increased efficiency and accuracy of spliceosomes for pre-mRNA splicing. Since intron diversity is critical for functions in RNA stability, regulation of gene expression and alternative splicing, RNA-binding proteins, spliceosomal regulatory factors and post-translational modifications of splicing factors ought to make the splicing process intron-specific. We recently reported function and regulation of a ubiquitin fold harboring splicing regulator, Sde2, which following activation by ubiquitin-specific proteases facilitates excision of selected introns from a subset of multi-intronic genes in Schizosaccharomyces pombe (Thakran et al. EMBO J, https://doi.org/10.15252/embj.201796751 , 2017). By reviewing our findings with understandings of intron functions and regulated splicing processes, we propose possible functions and mechanism of intron-specific pre-mRNA splicing and suggest that this process is crucial to highlight importance of introns in eukaryotic genomes.

FUNCTIONAL ANALYSIS OF INTRONS FOR THEIR EFFECT ON TRANSGENE EXPRESSION IN PLANT

Most plant genes contain intervening sequences known as introns that are found in nearly all plant genes and they are transcribed into pre-mRNA, the later will be removed by splicing machinery. Proteins are encoded by exons which separated by intron segments while the non-coding exons are conserved in the mature transcript. In this study, level of spliced and non-spliced introns was measured in transgenic Nicotiana tabacum lines expressing gusint-gfpint fusion gene controlled by the CaMV 35S promoter. Transcription analysis of the transgenes was performed by (RT)-PCR. All transcripts of transgenic lines generated two different length of fragments for gfp and gus genes. The results indicated that the gene expression of transcripts (without intron) in which their intron was spliced out, was higher than those of non-spliced transcript (with intron) in which the intron was maintained i.e. the both introns inside gus and gfp were partially and not completely spliced out of the transcripts. Furthermore, the relative quantifications of band intensity of gfp transcripts revealed that the percentage of non-spliced transcripts compared to the housekeeping gene, elongation factor 1-alpha (ef1α) using different amplification cycles (25, 30, 35 and 40) was 23.3 %, 28.3 %, 21.7 % and 20 % respectively and the percentage of spliced transcripts was 76.7 %, 71.7%, 78.3 % and 80 % in that order. While, the proportion of gus transcripts revealed that the non-spliced transcripts were 34.2 %, 15.7%, 12.8 % and 16.3 % and the spliced was 65.8 %, 84.3 %, 87.2 %, and 80 % at 25, 30, 35 and 40 cycles respectively. Conclusively, in the present study, the expression of gus and gfp genes might linked tothe process of splicing and not to a particular intron sequence or the partial splicing may due to splicing signals within the sequences of the both introns. In conclusion, in spite of the incomplete splicing of the IV2 second intron of the potato ST-LS-1 gene and fristintron of a castor bean gene for catalase, the expression of gus and gfp genes was not affected. Where, it was ranged from 71.7 % to 80 % and from 65.8 to 87.2 for the gus and gfp, respectively in the spliced introns. Therefore, the expression of the non-spliced gus and gfp introns was much lower than those of splicedones.

Identification and characterization of circular RNAs in the silkworm midgut following Bombyx mori cytoplasmic polyhedrosis virus infection

ABSTRACTThe pathogenesis of Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) infection is unclear, although accumulating evidence indicates that circular RNAs (circRNAs), which act as competing endogenous RNAs or positive regulators, play important roles in regulating gene expression in eukaryotes and, thus, may play a role in BmCPV infections. To explore the expression and biological functions of circRNAs in the silkworm midgut following BmCPV infection, silkworm circRNA expression profiles of normal midgut tissue (control) and BmCPV-infected midgut tissue (test) were determined using high-through sequencing. A total of 9,753 and 7,475circRNAs were detected from the control and test samples, respectively. The two samples shared 6,085 circRNAs, while 646 and 737 circRNAs were expressed specifically in the control and test samples, respectively. A total of 3,638 circRNAs were shown to be differentially expressed, and 400 circRNAs were substantially differentially expressed with a fold-change ≥ 2.0 (p< 0.05 and a false discover rate < 0.05), of which 294 were up-regulated and 106 were down-regulated following infection. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted to determine the principal functions of the substantially differentially regulated genes. circRNA-miRNA interaction networks were constructed based on a correlation analysis between the differentially expressed circRNAs and the nature of their microRNA (miRNA) binding sites. The network inferred that 13 miRNAs interacting with 193 circRNAs were among the 300 most abundant relationships. bmo-miR-3389-5p, bmo-miR-745-3p, and bmo-miR-3262 were related to 30, 34, and 34 circRNAs, respectively. circRNA_8115, circRNA_9444, circRNA_4553, circRNA_0827, and circRNA_6649 contained six, five, four, four, and four miRNA binding sites, respectively. We further found that alternative circularization of circRNAs is a common feature in silkworms and that the junction sites of many silkworm circRNAs are flanked by canonical GT/AG splicing signals. Our study is the first to show the circRNA response to virus infection. Thus, it provides a novel perspective on circRNA-miRNA interactions during BmCPV pathogenesis, and it lays the foundation for future research of the potential roles of circRNAs in BmCPV pathogenesis.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing.

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3’ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5’ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5’ss in dtd1+ intron 1 and of an upstream alternative 3’ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5’ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5’ ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3’ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

Cas9/gRNA targeted excision of cystic fibrosis-causing deep-intronic splicing mutations restores normal splicing of CFTR mRNA.

Cystic Fibrosis is an autosomal recessive disorder caused by mutations in the CFTR gene. CRISPR mediated, template-dependent homology-directed gene editing has been used to correct the most common mutation, c.1521_1523delCTT / p.Phe508del (F508del) which affects ~70% of individuals, but the efficiency was relatively low. Here, we describe a high efficiency strategy for editing of three different rare CFTR mutations which together account for about 3% of individuals with Cystic Fibrosis. The mutations cause aberrant splicing of CFTR mRNA due to the creation of cryptic splice signals that result in the formation of pseudoexons containing premature stop codons c.1679+1634A>G (1811+1.6kbA>G) and c.3718-2477C>T (3849+10kbC>T), or an out-of-frame 5’ extension to an existing exon c.3140-26A>G (3272-26A>G). We designed pairs of Cas9 guide RNAs to create targeted double-stranded breaks in CFTR either side of each mutation which resulted in high efficiency excision of the target genomic regions via non-homologous end-joining repair. When evaluated in a mini-gene splicing assay, we showed that targeted excision restored normal splicing for all three mutations. This approach could be used to correct aberrant splicing signals or remove disruptive transcription regulatory motifs caused by deep-intronic mutations in a range of other genetic disorders.

Large-scale prediction of ADAR-mediated effective human A-to-I RNA editing.

Adenosine-to-inosine (A-to-I) editing by adenosine deaminase acting on the RNA (ADAR) proteins is one of the most frequent modifications during post- and co-transcription. To facilitate the assignment of biological functions to specific editing sites, we designed an automatic online platform to annotate A-to-I RNA editing sites in pre-mRNA splicing signals, microRNAs (miRNAs) and miRNA target untranslated regions (3' UTRs) from human (Homo sapiens) high-throughput sequencing data and predict their effects based on large-scale bioinformatic analysis. After analysing plenty of previously reported RNA editing events and human normal tissues RNA high-seq data, >60 000 potentially effective RNA editing events on functional genes were found. The RNA Editing Plus platform is available for free at https://www.rnaeditplus.org/, and we believe our platform governing multiple optimized methods will improve further studies of A-to-I-induced editing post-transcriptional regulation.

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.

Neurofibromatosis type 1 alternative splicing is a key regulator of Ras/ERK signaling and learning behaviors in mice.

Appropriate activation of the Ras/extracellular signal-regulated kinase (ERK) protein signaling cascade within the brain is crucial for optimal learning and memory. One key regulator of this cascade is the Nf1 Ras GTPase activating protein (RasGAP), which attenuates Ras/ERK signaling by converting active Ras is bound to guanosine triphosphate, activating Ras into inactive Ras is bound to guanosine diphosphate, inactivating Ras. A previous study using embryonic stem cells and embryonic stem cell-derived neurons indicated that Nf1 RasGAP activity is modulated by the highly regulated alternative splicing of Nf1 exon 23a. In this study, we generated Nf123aIN/23aIN mice, in which the splicing signals surrounding Nf1 exon 23a were manipulated to increase exon inclusion. Nf123aIN/23aIN mice are viable and exon 23a inclusion approaches 100% in all tissues, including the brain, where the exon is normally almost completely skipped. Ras activation and phosphorylation of ERK1/2 downstream of Ras are both greatly increased in Nf123aIN/23aIN mouse brain lysates, confirming that exon 23a inclusion inhibits Nf1 RasGAP activity in vivo as it does in cultured cells. Consistent with the finding of altered Ras/ERK signaling in the brain, Nf123aIN/23aIN mice showed specific deficits in learning and memory compared with Nf1+/+ mice. Nf123aIN/23aIN mice performed poorly on the T-maze and Morris water maze tests, which measure short- and long-term spatial memory, respectively. In addition, Nf123aIN/23aIN mice showed abnormally elevated context-dependent fear and a diminished ability to extinguish a cued fear response, indicating defective associative fear learning. Therefore, the regulated alternative splicing of Nf1 is an important mechanism for fine-tuning Ras/ERK signaling as well as learning and memory in mice.

Identification of pathogenic gene mutations in LMNA and MYBPC3 that alter RNA splicing

Genetic variants that cause haploinsufficiency account for many autosomal dominant (AD) disorders. Gene-based diagnosis classifies variants that alter canonical splice signals as pathogenic, but due to imperfect understanding of RNA splice signals other variants that may create or eliminate splice sites are often clinically classified as variants of unknown significance (VUS). To improve recognition of pathogenic splice-altering variants in AD disorders, we used computational tools to prioritize VUS and developed a cell-based minigene splicing assay to confirm aberrant splicing. Using this two-step procedure we evaluated all rare variants in two AD cardiomyopathy genes, lamin A/C (LMNA) and myosin binding protein C (MYBPC3). We demonstrate that 13 LMNA and 35 MYBPC3 variants identified in cardiomyopathy patients alter RNA splicing, representing a 50% increase in the numbers of established damaging splice variants in these genes. Over half of these variants are annotated as VUS by clinical diagnostic laboratories. Familial analyses of one variant, a synonymous LMNA VUS, demonstrated segregation with cardiomyopathy affection status and altered cardiac LMNA splicing. Application of this strategy should improve diagnostic accuracy and variant classification in other haploinsufficient AD disorders.

Antisense Oligonucleotide-Mediated Removal of the Polyglutamine Repeat in Spinocerebellar Ataxia Type 3 Mice.

Spinocerebellar ataxia type 3 (SCA3) is a currently incurable neurodegenerative disorder caused by a CAG triplet expansion in exon 10 of the ATXN3 gene. The resultant expanded polyglutamine stretch in the mutant ataxin-3 protein causes a gain of toxic function, which eventually leads to neurodegeneration. One important function of ataxin-3 is its involvement in the proteasomal protein degradation pathway, and long-term downregulation of the protein may therefore not be desirable. In the current study, we made use of antisense oligonucleotides to mask predicted exonic splicing signals, resulting in exon 10 skipping from ATXN3 pre-mRNA. This led to formation of a truncated ataxin-3 protein lacking the toxic polyglutamine expansion, but retaining its ubiquitin binding and cleavage function. Repeated intracerebroventricular injections of the antisense oligonucleotides in a SCA3 mouse model led to exon skipping and formation of the modified ataxin-3 protein throughout the mouse brain. Exon skipping was long lasting, with the modified protein being detectable for at least 2.5 months after antisense oligonucleotide injection. A reduction in insoluble ataxin-3 and nuclear accumulation was observed following antisense oligonucleotide treatment, indicating a beneficial effect on pathogenicity. Together, these data suggest that exon 10 skipping is a promising therapeutic approach for SCA3.

ISVASE: identification of sequence variant associated with splicing event using RNA-seq data

BackgroundExon recognition and splicing precisely and efficiently by spliceosome is the key to generate mature mRNAs. About one third or a half of disease-related mutations affect RNA splicing. Software PVAAS has been developed to identify variants associated with aberrant splicing by directly using RNA-seq data. However, it bases on the assumption that annotated splicing site is normal splicing, which is not true in fact.ResultsWe develop the ISVASE, a tool for specifically identifying sequence variants associated with splicing events (SVASE) by using RNA-seq data. Comparing with PVAAS, our tool has several advantages, such as multi-pass stringent rule-dependent filters and statistical filters, only using split-reads, independent sequence variant identification in each part of splicing (junction), sequence variant detection for both of known and novel splicing event, additional exon-exon junction shift event detection if known splicing events provided, splicing signal evaluation, known DNA mutation and/or RNA editing data supported, higher precision and consistency, and short running time. Using a realistic RNA-seq dataset, we performed a case study to illustrate the functionality and effectiveness of our method. Moreover, the output of SVASEs can be used for downstream analysis such as splicing regulatory element study and sequence variant functional analysis.ConclusionsISVASE is useful for researchers interested in sequence variants (DNA mutation and/or RNA editing) associated with splicing events. The package is freely available at https://sourceforge.net/projects/isvase/.

Potentially pathogenic germline CHEK2 c.319+2T>A among multiple early-onset cancer families.

To study the potential contribution of genes other than BRCA1/2, PTEN, and TP53 to the biological and clinical characteristics of multiple early-onset cancers in Norwegian families, including early-onset breast cancer, Cowden-like and Li-Fraumeni-like syndromes (BC, CSL and LFL, respectively). The Hereditary Cancer Biobank from the Norwegian Radium Hospital was used to identify early-onset BC, CSL or LFL for whom no pathogenic variants in BRCA1/2, PTEN, or TP53 had been found in routine diagnostic DNA sequencing. Forty-four cancer susceptibility genes were selected and analyzed by our in-house designed TruSeq amplicon-based assay for targeted sequencing. Protein- and RNA splicing-dedicated in silico analyses were performed for all variants of unknown significance (VUS). Variants predicted as the more likely to affect splicing were experimentally analyzed by minigene assay. We identified a CSL individual carrying a variant in CHEK2 (c.319+2T>A, IVS2), here considered as likely pathogenic. Out of the five VUS (BRCA2, CDH1, CHEK2, MAP3K1, NOTCH3) tested in the minigene splicing assay, only NOTCH3 c.14090C>T (p.Ser497Leu) showed a significant effect on RNA splicing, notably by inducing partial skipping of exon 9. Among 13 early-onset BC, CSL and LFL patients, gene panel sequencing identified a potentially pathogenic variant in CHEK2 that affects a canonical RNA splicing signal. Our study provides new information on genetic loci that may affect the risk of developing cancer in these patients and their families, demonstrating that genes presently not routinely tested in molecular diagnostic settings may be important for capturing cancer predisposition in these families.

Influence of low frequency PSEN1 variants on familial Alzheimer's disease risk in Brazil.

About 30-70% of familial Alzheimer's disease (AD) cases are related to mutations in presenilin-1 gene (PSEN1). Although the role of mutations and common variants in AD had been extensively investigated, the contribution of rare or low frequency PSEN1 variants on AD risk remains unclear. In the current study, we performed a mutational screening of PSEN1 coding exons and flanking intronic sequences among 53 index cases with familial history of AD from Rio de Janeiro (Brazil). Two missense variants (rs63750592; rs17125721), one rare and a low frequency variant, and two intronic variants (rs3025786; rs165932) were identified. In silico tools were used to predict the functional impact of the variants, revealing no changes in protein functionality by exonic variants. Otherwise, all variants were predicted to alter splicing signals. Prediction results, together with previous reports, suggest a correlation between rs17125721 and AD. So, a subsequent case-control study to evaluate the role of rs1712572 on AD risk was performed in an additional sample of 120 AD sporadic cases and in 149 elderly healthy controls by TaqMan Genotyping Assay. Our data indicates a risk association for rs17125721 in familial AD cases (OR=6.0; IC95%=1.06-33.79; p=0.042). In addition, we tested the multiplicative interaction between allele ε4 of the apolipoprotein E (APOE) and rs17125721 and no statistical association was found. Taken together, our findings provide new insight about the genetic relevance of low frequency PSEN1 variants for familial AD development.