Activation of cryptic splice sites is a frequent splicing defect mechanism caused by mutations in exon and intron sequences of the OPA1 gene

Mutations in OPA1 are the most frequent cause underlying autosomal dominant optic atrophy (adOA). Until now only few putative splicing mutations in the OPA1 gene have been investigated at the mRNA level and all these result in exon skipping. Here, we report the identification and cDNA analysis of four intronic and three exonic OPA1 gene mutations that cause a variety of splicing defects including activation of cryptic splice sites in either flanking exon or intron sequences, and a leaky splicing mutation. Our results show that cDNA analysis is of prime importance for the full evaluation of the effect of putative splicing mutations in the OPA1 gene.     

Both U2 snRNA and U12 snRNA are required for accurate splicing of exon 5 of the rat calcitonin/CGRP gene

Two classes of spliceosome are present in eukaryotic cells. Most introns in nuclear pre-mRNAs are removed by a spliceosome that requires U1, U2, U4, U5, and U6 small nuclear ribonucleoprotein particles (snRNPs). A minor class of introns are removed by a spliceosome containing U11, U12, U5, U4atac, and U6 atac snRNPs. We describe experiments that demonstrate that splicing of exon 5 of the rat calcitonin/CGRP gene requires both U2 snRNA and U12 snRNA. In vitro, splicing to calcitonin/ CGRP exon 5 RNA was dependent on U2 snRNA, as preincubation of nuclear extract with an oligonucleotide complementary to U2 snRNA abolished exon 5 splicing. Addition of an oligonucleotide complementary to U12 snRNA increased splicing at a cryptic splice site in exon 5 from <5% to 50% of total spliced RNA. Point mutations in a candidate U12 branch sequence in calcitonin/CGRP intron 4, predicted to decrease U12-pre-mRNA base-pairing, also significantly increased cryptic splicing in vitro. Calcitonin/CGRP genes containing base changes disrupting the U12 branch sequence expressed significantly decreased CGRP mRNA levels when expressed in cultured cells. Coexpression of U12 snRNAs containing base changes predicted to restore U12-pre-mRNA base pairing increased CGRP mRNA synthesis to the level of the wild-type gene. These observations indicate that accurate, efficient splicing of calcitonin/CGRP exon 5 is dependent upon both U2 and U12 snRNAs.     

An Alu-derived intronic splicing enhancer facilitates intronic processing and modulates aberrant splicing in ATM

We have previously reported a natural GTAA deletion within an intronic splicing processing element (ISPE) of the ataxia telangiectasia mutated (ATM) gene that disrupts a non-canonical U1 snRNP interaction and activates the excision of the upstream portion of the intron. The resulting pre-mRNA splicing intermediate is then processed to a cryptic exon, whose aberrant inclusion in the final mRNA is responsible for ataxia telangiectasia. We show here that the last 40 bases of a downstream intronic antisense Alu repeat are required for the activation of the cryptic exon by the ISPE deletion. Evaluation of the pre-mRNA splicing intermediate by a hybrid minigene assay indicates that the identified intronic splicing enhancer represents a novel class of enhancers that facilitates processing of splicing intermediates possibly by recruiting U1 snRNP to defective donor sites. In the absence of this element, the splicing intermediate accumulates and is not further processed to generate the cryptic exon. Our results indicate that Alu-derived sequences can provide intronic splicing regulatory elements that facilitate pre-mRNA processing and potentially affect the severity of disease-causing splicing mutations.     

Optimizing high-resolution melting analysis for the detection of mutations of GPR30/GPER-1 in breast cancer

G protein-coupled receptor 30/G protein estrogen receptor-1 (GPR30/GPER-1) is a novel membrane receptor for estrogen whose mRNA is expressed at high levels in estrogen-dependent cells such as breast cancer cell lines. However, mutations in GRP30 related to diseases remain unreported. To detect unknown mutations in the GPR30 open reading frame (ORF) quickly, the experimental conditions for high-resolution melting (HRM) analysis were examined for PCR primers, Taq polymerases, saturation DNA binding dyes, Mg(2+) concentration, and normalized temperatures. Nine known SNPs and 13 artificial point mutations within the GPR30 ORF, as well as single nucleotide variants in DNA extracted from subjects with breast cancers were tested under the optimal experimental conditions. The combination of Expand High Fidelity(PLUS) and SYTO9 in the presence of 2.0 mM MgCl(2) produced the best separation in melting curves of mutations in all regions of the GPR30 ORF. Under these experimental conditions, the mutations were clearly detected in both heterozygotes and homozygotes. HRM analysis of GPR30 using genomic DNA from subjects with breast cancers showed a novel single nucleotide variant, 111C>T in GPR30 and 4 known SNPs. The experimental conditions determined in this study for HRM analysis are useful for high throughput assays to detect unknown mutations within the GPR30 ORF.     

Editor's choice: Functional classification of 15 million SNPs detected from diverse chicken populations

Next-generation sequencing has prompted a surge of discovery of millions of genetic variants from vertebrate genomes. Besides applications in genetic association and linkage studies, a fraction of these variants will have functional consequences. This study describes detection and characterization of 15 million SNPs from chicken genome with the goal to predict variants with potential functional implications (pfVars) from both coding and non-coding regions. The study reports: 183K amino acid-altering SNPs of which 48% predicted as evolutionary intolerant, 13K splicing variants, 51K likely to alter RNA secondary structures, 500K within most conserved elements and 3K from non-coding RNAs. Regions of local fixation within commercial broiler and layer lines were investigated as potential selective sweeps using genome-wide SNP data. Relationships with phenotypes, if any, of the pfVars were explored by overlaying the sweep regions with known QTLs. Based on this, the candidate genes and/or causal mutations for a number of important traits are discussed. Although the fixed variants within sweep regions were enriched with non-coding SNPs, some non-synonymous-intolerant mutations reached fixation, suggesting their possible adaptive advantage. The results presented in this study are expected to have important implications for future genomic research to identify candidate causal mutations and in poultry breeding.     

Exon skipping in IVD RNA processing in isovaleric acidemia caused by point mutations in the coding region of the IVD gene

Isovaleric acidemia (IVA) is a recessive disorder caused by a deficiency of isovaleryl-CoA dehydrogenase (IVD). We have reported elsewhere nine point mutations in the IVD gene in fibroblasts of patients with IVA, which lead to abnormalities in IVD protein processing and activity. In this report, we describe eight IVD gene mutations identified in seven IVA patients that result in abnormal splicing of IVD RNA. Four mutations in the coding region lead to aberrantly spliced mRNA species in patient fibroblasts. Three of these are amino acid altering point mutations, whereas one is a single-base insertion that leads to a shift in the reading frame of the mRNA. Two of the coding mutations strengthen pre-existing cryptic splice acceptors adjacent to the natural splice junctions and apparently interfere with exon recognition, resulting in exon skipping. This mechanism for missplicing has not been reported elsewhere. Four other mutations alter either the conserved gt or ag dinucleotide splice sites in the IVD gene. Exon skipping and cryptic splicing were confirmed by transfection of these mutations into a Cos-7 cell line model splicing system. Several of the mutations were predicted by individual information analysis to inactivate or significantly weaken adjacent donor or acceptor sites. The high frequency of splicing mutations identified in these patients is unusual, as is the finding of missplicing associated with missense mutations in exons. These results may lead to a better understanding of the phenotypic complexity of IVA, as well as provide insight into those factors important in defining intron/exon boundaries in vivo.     

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5′ splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5′ and 3′ splice sites, and the activation of cryptic U12-dependent branch/3′ splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3′ splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.     

Identification of Novel GRM1 Mutations and Single Nucleotide Polymorphisms in Prostate Cancer Cell Lines and Tissues

Metabotropic glutamate receptor 1 (GRM1) signaling has been implicated in benign and malignant disorders including prostate cancer (PCa). To further explore the role of genetic alterations of GRM1 in PCa, we screened the entire human GRM1 gene including coding sequence, exon-intron junctions, and flanking untranslated regions (UTRs) for the presence of mutations and single nucleotide polymorphisms (SNPs) in several PCa cell lines and matched tumor-normal tissues from Caucasian Americans (CAs) and African Americans (AAs). We used bidirectional sequencing, allele-specific PCR, and bioinformatics to identify the genetic changes in GRM1 and to predict their functional role. A novel missense mutation identified at C1744T (582 Pro > Ser) position of GRM1 gene in a primary AA-PCa cell line (E006AA) was predicted to affect the protein stability and functions. Another novel mutation identified at exon-intron junction of exon-8 in C4-2B cell line resulted in alteration of the GRM1 splicing donor site. In addition, we found missense SNP at T2977C (993 Ser > Pro) position and multiple non-coding mutations and SNPs in 3′-UTR of GRM1 gene in PCa cell lines and tissues. These novel mutations may contribute to the disease by alterations in GRM1 gene splicing, receptor activation, and post-receptor downstream signaling.     

A Bayesian ensemble approach with a disease gene network predicts damaging effects of missense variants of human cancers

Large-scale sequencing of cancer genomes has revealed many novel mutations and inter-tumoral heterogeneity. Therefore, prioritizing variants according to their potential deleterious effects has become essential. We constructed a disease gene network and proposed a Bayesian ensemble approach that integrates diverse sources to predict the functional effects of missense variants. We analyzed 23,336 missense disease mutations and 36,232 neutral polymorphisms of 12,039 human proteins. The results showed successful improvement of prediction accuracy in both sensitivity and specificity, and we demonstrated the utility of the method by applying it to somatic mutations obtained from colorectal and breast cancer cell lines. The candidate genes with predicted deleterious mutations as well as known cancer genes were significantly enriched in many KEGG pathways related to carcinogenesis, supporting genetic homogeneity of cancer at the pathway level. The breast cancer-specific network increased the prediction accuracy for breast cancer mutations. This study provides a ranked list of deleterious mutations and candidate cancer genes and suggests that mutations affecting cancer may occur in important pathways and should be interpreted on the phenotype-related network or pathway. A disease gene network may be of value in predicting functional effects of novel disease-specific mutations.     

Identification of a novel splice variant of X-linked inhibitor of apoptosis-associated factor 1

XAF1 (XIAP-associated factor 1) binds to XIAP and blocks its anti-apoptotic activity. It has been reported that XAF1 is mainly expressed in normal tissues but is missing or present at low levels in most cancer cell lines, which implies a tumor-suppressing function. In the present study we describe the identification of a novel splice variant of human XAF1, designated XAF1C, which contains a cryptic exon. Incorporation of this exon (exon 4b) into the mRNA introduces an in-frame stop codon, resulting in a shortened open-reading frame (ORF) of 495 nucleotides. This ORF is predicted to encode a 164 amino acid (AA) protein lacking the C-terminal domain of the previously described XAF1(A), but containing a unique 24 AA carboxy terminus. Like XAF1(A), XAF1C mRNA expression was detected in a variety of human cancer cell lines and also in normal human tissues. The ratio of XAF1(A) and XAF1C mRNA expression differs amongst the cell lines tested, suggesting differential mRNA stabilities and/or the existence of tissue- or cell type-specific splicing regulation. In transfected cells, xaf1c encodes a truncated protein of 18kDa, which is distributed primarily in the nucleus.     

The Alternative Splicing of Cytoplasmic Polyadenylation Element Binding Protein 2 Drives Anoikis Resistance and the Metastasis of Triple Negative Breast Cancer

Triple negative breast cancer (TNBC) represents an anomalous subset of breast cancer with a greatly reduced (30%) 5-year survival rate. The enhanced mortality and morbidity of TNBC arises from the high metastatic rate, which requires the acquisition of AnR, a process whereby anchorage-dependent cells become resistant to cell death induced by detachment. In this study TNBC cell lines were selected for AnR, and these cell lines demonstrated dramatic enhancement in the formation of lung metastases as compared with parental cells. Genetic analysis of the AnR subclones versus parental cells via next generation sequencing and analysis of global alternative RNA splicing identified that the mRNA splicing of cytoplasmic polyadenylation element binding 2 (CPEB2), a translational regulator, was altered in AnR TNBC cells. Specifically, increased inclusion of exon 4 into the mature mRNA to produce the CPEB2B isoform was observed in AnR cell lines. Molecular manipulations of CPEB2 splice variants demonstrated a key role for this RNA splicing event in the resistance of cells to anoikis. Specifically, down-regulation of the CPEB2B isoform using siRNA re-sensitized the AnR cell lines to detachment-induced cell death. The ectopic expression of CPEB2B in parental TNBC cell lines induced AnR and dramatically increased metastatic potential. Importantly, alterations in the alternative splicing of CPEB2 were also observed in human TNBC and additional subtypes of human breast cancer tumors linked to a high metastatic rate. Our findings demonstrate that the regulation of CPEB2 mRNA splicing is a key mechanism in AnR and a driving force in TNBC metastasis.     

Single nucleotide polymorphism-based validation of exonic splicing enhancers

Because deleterious alleles arising from mutation are filtered by natural selection, mutations that create such alleles will be underrepresented in the set of common genetic variation existing in a population at any given time. Here, we describe an approach based on this idea called VERIFY (variant elimination reinforces functionality), which can be used to assess the extent of natural selection acting on an oligonucleotide motif or set of motifs predicted to have biological activity. As an application of this approach, we analyzed a set of 238 hexanucleotides previously predicted to have exonic splicing enhancer (ESE) activity in human exons using the relative enhancer and silencer classification by unanimous enrichment (RESCUE)-ESE method. Aligning the single nucleotide polymorphisms (SNPs) from the public human SNP database to the chimpanzee genome allowed inference of the direction of the mutations that created present-day SNPs. Analyzing the set of SNPs that overlap RESCUE-ESE hexamers, we conclude that nearly one-fifth of the mutations that disrupt predicted ESEs have been eliminated by natural selection (odds ratio = 0.82 +/- 0.05). This selection is strongest for the predicted ESEs that are located near splice sites. Our results demonstrate a novel approach for quantifying the extent of natural selection acting on candidate functional motifs and also suggest certain features of mutations/SNPs, such as proximity to the splice site and disruption or alteration of predicted ESEs, that should be useful in identifying variants that might cause a biological phenotype.