Phenotypic consequences of branch point substitutions

Exploiting public databases of genomic variation to quantify evolutionary constraint on the branch point sequence in 30 plant and animal species

The branch point sequence is a degenerate intronic heptamer required for the assembly of the spliceosome during pre-mRNA splicing. Disruption of this motif may promote alternative splicing and eventually cause phenotype variation. Despite its functional relevance, the branch point sequence is not included in most genome annotations. Here, we predict branch point sequences in 30 plant and animal species and attempt to quantify their evolutionary constraints using public variant databases. We find an implausible variant distribution in the databases from 16 of 30 examined species. Comparative analysis of variants from whole-genome sequencing shows that variants submitted from exome sequencing or false positive variants are widespread in public databases and cause these irregularities. We then investigate evolutionary constraint with largely unbiased public variant databases in 14 species and find that the fourth and sixth position of the branch point sequence are more constrained than coding nucleotides. Our findings show that public variant databases should be scrutinized for possible biases before they qualify to analyze evolutionary constraint.

Cytokine gene polymorphism and parasite susceptibility in free-living rodents: Importance of non-coding variants

Associations between genetic variants and susceptibility to infections have long been studied in free-living hosts so as to infer the contemporary evolutionary forces that shape the genetic polymorphisms of immunity genes. Despite extensive studies of proteins interacting with pathogen-derived ligands, such as MHC (major histocompatilbility complex) or TLR (Toll-like receptors), little is known about the efferent arm of the immune system. Cytokines are signalling molecules that trigger and modulate the immune response, acting as a crucial link between innate and adaptive immunity. In the present study we investigated how genetic variation in cytokines in bank voles Myodes glareolus affects their susceptibility to infection by parasites (nematodes: Aspiculuris tianjensis, Heligmosomum mixtum, Heligmosomoides glareoli) and microparasites (Cryptosporidium sp, Babesia microti, Bartonella sp.). We focused on three cytokines: tumour necrosis factor (TNF), lymphotoxin alpha (LTα), and interferon beta (IFNβ1). Overall, we identified four single nucleotide polymorphisms (SNPs) associated with susceptibility to nematodes: two located in LTα and two in IFNβ1. One of those variants was synonymous, another located in an intron. Each SNP associated with parasite load was located in or next to a codon under selection, three codons displayed signatures of positive selection, and one of purifying selection. Our results indicate that cytokines are prone to parasite-driven selection and that non-coding variants, although commonly disregarded in studies of the genetic background of host-parasite co-evolution, may play a role in susceptibility to infections in wild systems.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Alternatively Spliced Isoforms of the P2X7 Receptor: Structure, Function and Disease Associations

The P2X7 receptor (P2X7R) is an ATP-gated membrane ion channel that is expressed by multiple cell types. Following activation by extracellular ATP, the P2X7R mediates a broad range of cellular responses including cytokine and chemokine release, cell survival and differentiation, the activation of transcription factors, and apoptosis. The P2X7R is made up of three P2X7 subunits that contain specific domains essential for the receptor’s varied functions. Alternative splicing produces P2X7 isoforms that exclude one or more of these domains and assemble in combinations that alter P2X7R function. The modification of the structure and function of the P2X7R may adversely affect cellular responses to carcinogens and pathogens, and alternatively spliced (AS) P2X7 isoforms have been associated with several cancers. This review summarizes recent advances in understanding the structure and function of AS P2X7 isoforms and their associations with cancer and potential role in modulating the inflammatory response.

The intronic branch point sequence is under strong evolutionary constraint in the bovine and human genome

The branch point sequence is a cis-acting intronic motif required for mRNA splicing. Despite their functional importance, branch point sequences are not routinely annotated. Here we predict branch point sequences in 179,476 bovine introns and investigate their variability using a catalogue of 29.4 million variants detected in 266 cattle genomes. We localize the bovine branch point within a degenerate heptamer "nnyTrAy". An adenine residue at position 6, that acts as branch point, and a thymine residue at position 4 of the heptamer are more strongly depleted for mutations than coding sequences suggesting extreme purifying selection. We provide evidence that mutations affecting these evolutionarily constrained residues lead to alternative splicing. We confirm evolutionary constraints on branch point sequences using a catalogue of 115 million SNPs established from 3,942 human genomes of the gnomAD database.

A Pharmacogenetic Study of CYP2C19 in Acute Coronary Syndrome Patients of Colombian Origin Reveals New Polymorphisms Potentially Related to Clopidogrel Therapy

Clopidogrel, an oral platelet P2Y₁₂ receptor blocker, is used in the treatment of acute coronary syndrome. Interindividual variability in treatment response and the occurrence of adverse effects has been attributed to genetic variants in CYP2C19. The analysis of relevant pharmacogenes in ethnically heterogeneous and poorly studied populations contributes to the implementation of personalized medicine. We analyzed the coding and regulatory regions of CYP2C19 in 166 patients with acute coronary syndrome (ACS) treated with clopidogrel. The allele frequencies of CYP2C19 alleles *1, *2, *4, *17, *27 and *33 alleles were 86.1%, 7.2%, 0.3%, 10.2%, 0.3% and 0.3%, respectively. A new potentially pathogenic mutation (p.L15H) and five intronic variants with potential splicing effects were detected. In 14.4% of the patients, a new haplotype in strong linkage disequilibrium was identified. The clinical outcome indicated that 13.5% of the patients presented adverse drugs reactions with a predominance of bleeding while 25% of these patients were carriers of at least one polymorphic allele. We propose that new regulatory single-nucleotide variants (SNVs) might potentially influence the response to clopidogrel in Colombian individuals.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

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.

Impact of matrix metalloproteinase-11 gene polymorphisms on development and clinicopathologcial variables of uterine cervical cancer in Taiwanese women

The purposes of this study were to examine whether there were associations among matrix metalloproteinase-11 (MMP-11) gene polymorphisms, development and clinicopathological characteristics of uterine cervical cancer as well as patient survival or not. Five single-nucleotide polymorphisms (SNPs) of the MMP-11 gene rs738791, rs738792, rs2267029, rs28382575, and rs131451 from one hundred and thirty patients with invasive cancer, 99 patients with high-grade cervical intraepithelial neoplasia (CIN) of uterine and 335 normal controls were analyzed using real-time polymerase chain reaction. Our results revealed that genotypic frequencies of CT/TT in MMP-11 SNP rs738791, with CC as a reference, tended to exhibit significantly different distributions (p=0.044, AOR: 0.63, 95% CI: 0.41-0.99) between patients with cervical invasive cancer and normal control women when controlling age. After multiple significance adjustment, the tendency becomes insignificant (Holm's adjusted p 0.176). Although CT/TT genotype of MMP-11 gene rs738791 tended to increase the risk of developing stage II disease at least (p=0.035; OR: 2.16, 95% CI: 1.05-4.44) and deep stromal invasion more than 10 mm (p=0.043; OR: 2.08, 95% CI: 1.02-4.26) with CC as a reference in patients with uterine cervical cancer. They became insignificant after multiple significance adjustment and the Holm's adjusted p values would become as 0.245 and 0.258, respectively. However, lymph node metastasis exhibited significant worse recurrence-free survival (p=0.033; HR: 2.83, 95% CI: 1.09-7.35), and overall survival (p=0.001; HR: 4.80, 95% CI: 1.82-12.62) compared to those without pelvic lymph node metastasis. In conclusion, it indicates no impact of the MMP-11 SNPs on uterine cervical cancer in Taiwanese women.

Novel insights into FAS defects underlying autoimmune lymphoproliferative syndrome revealed by studies in consanguineous patients

Autoimmune lymphoproliferative syndrome (ALPS) is a primary immunodeficiency disease due to impaired Fas-Fas ligand apoptotic pathway. It is characterized by chronic nonmalignant, noninfectious lymphadenopathy and/or splenomegaly associated with autoimmune manifestations primarily directed against blood cells. Herein, we review the heterogeneous ALPS molecular bases and discuss recent findings revealed by the study of consanguineous patients. Indeed, this peculiar genetic background favored the identification of a novel form of AR ALPS-FAS associated with normal or residual protein expression, expanding the spectrum of ALPS types. In addition, rare mutational mechanisms underlying the splicing defects of FAS exon 6 have been identified in AR ALPS-FAS with lack of protein expression. These findings will help decipher critical regions required for the tight regulation of FAS exon 6 splicing. We also discuss the genotype-phenotype correlation and disease severity in AR ALPS-FAS. Altogether, the study of ALPS molecular bases in endogamous populations helps to better classify the disease subgroups and to unravel the Fas pathway functioning.

Identification of functional single nucleotide polymorphisms in the branchpoint site

Background

The human genome contains millions of single nucleotide polymorphisms (SNPs); many of these SNPs are intronic and have unknown functional significance. SNPs occurring within intron branchpoint sites, especially at the adenine (A), would presumably affect splicing; however, this has not been systematically studied. We employed a splicing prediction tool to identify human intron branchpoint sites and screened dbSNP for identifying SNPs located in the predicted sites to generate a genome-wide branchpoint site SNP database.

Results

We identified 600 SNPs located within branchpoint sites; among which, 216 showed a change in A. After scoring the SNPs by counting the As in the ± 10 nucleotide region, only four SNPs were identified without additional As (rs13296170, rs12769205, rs75434223, and rs67785924). Using minigene constructs, we examined the effects of these SNPs on splicing. The three SNPs (rs13296170, rs12769205, and rs75434223) with nucleotide substitution at the A position resulted in abnormal splicing (exon skipping and/or intron inclusion). However, rs67785924, a 5-bp deletion that abolished the branchpoint A nucleotide, exhibited normal RNA splicing pattern, presumably using two of the downstream As as alternative branchpoints. The influence of additional As on splicing was further confirmed by studying rs2733532, which contains three additional As in the ± 10 nucleotide region.

Conclusions

We generated a high-confidence genome-wide branchpoint site SNP database, experimentally verified the importance of A in the branchpoint, and suggested that other nearby As can protect branchpoint A substitution from abnormal splicing.

Electronic supplementary material

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

Rare splicing defects of FAS underly severe recessive autoimmune lymphoproliferative syndrome

Autoimmune lymphoproliferative syndrome (ALPS) is a prototypic disorder of impaired apoptosis characterized by autoimmune features and lymphoproliferation. Heterozygous germline or somatic FAS mutations associated with preserved protein expression have been described. Very rare cases of homozygous germline FAS mutations causing severe autosomal recessive form of ALPS with a complete defect of Fas expression have been reported. We report two unrelated patients from highly inbred North African population showing a severe ALPS phenotype and an undetectable Fas surface expression. Two novel homozygous mutations have been identified underlying rare splicing defects mechanisms. The first mutation breaks a branch point sequence and the second alters a regulatory exonic splicing site. These splicing defects induce the skipping of exon 6 encoding the transmembrane domain of CD95. Our findings highlight the requirement of tight regulation of FAS exon 6 splicing for balanced alternative splicing and illustrate the importance of such studies in highly consanguineous populations.

Alport syndrome cold cases: Missing mutations identified by exome sequencing and functional analysis

Alport syndrome (AS) is an inherited progressive renal disease caused by mutations in COL4A3, COL4A4, and COL4A5 genes. Despite simultaneous screening of these genes being widely available, mutation detection still remains incomplete in a non-marginal portion of patients. Here, we applied whole-exome sequencing (WES) in 3 Italian families negative after candidate-gene analyses. In Family 1, we identified a novel heterozygous intronic variant (c.2245-40A>G) -outside the conventionally screened candidate region for diagnosis- potentially disrupting COL4A5 exon29 splicing. Using a minigene-based approach in HEK293 cells we demonstrated that this variant abolishes exon29 branch site, causing exon skipping. Moreover, skewed X-inactivation of the c.2245-40A>G allele correlated with disease severity in heterozygous females. In Family 2, WES highlighted a novel COL4A5 hemizygous missense mutation (p.Gly491Asp), which segregates with the phenotype and impacts on a highly-conserved residue. Finally, in Family 3, we detected a homozygous 24-bp in-frame deletion in COL4A3 exon1 (NM_000091.4:c.30_53del:p.Val11_Leu18del or c.40_63del24:p.Leu14_Leu21del), which is ambiguously annotated in databases, although it corresponds to a recurrent AS mutation. Functional analyses showed that this deletion disrupts COL4A3 signal peptide, possibly altering protein secretion. In conclusion, WES -together with functional studies- was fundamental for molecular diagnosis in 3 AS families, highlighting pathogenic variants that escaped previous screenings.

Limb girdle myasthenia with digenic RAPSN and a novel disease gene AK9 mutations

Though dysfunction of neuromuscular junction (NMJ) is associated with congenital myasthenic syndrome (CMS), the proteins involved in neuromuscular transmission have not been completely identified. In this study, we aimed to identify a novel CMS gene in a consanguineous family with limb-girdle type CMS. Homozygosity mapping of the novel CMS gene was performed using high-density single-nucleotide polymorphism microarrays. The variants in CMS gene were identified by whole-exome sequencing (WES) and Sanger sequencing. A 20 MB-region of homozygosity (ROH) was mapped on chromosome 6q15-21. This was the only ROH that present in all clinically affected siblings and absent in all clinically unaffected siblings. WES showed a novel variant of AK9 gene located in this ROH. This variant was a start-gain mutation and introduced a cryptic 5'-UTR signal in intron 5 of the AK9 gene. The normal splicing signal would be interfered by the cryptic translation signal leading to defective splicing. Another 25 MB-ROH was found on chromosome 11p13-q12 in all siblings. WES showed a homozygous RAPSN pathogenic variant in this ROH. Since RAPSN-associated limb-girdle type CMS was only manifested in AK9 homozygous variant carriers, the disease phenotype was of digenic inheritance, and was determined by the novel disease modifier AK9 which provides NTPs for N-glycosylation. This is the first time that this specific genotype-phenotype correlation is reported. Importantly, the AK9-associated nucleotide deficiency may replete by dietary supplements. Since AK9 is a disease modifier, enhancing N-glycosylation by increasing dietary nucleotides may be a new therapeutic option for CMS patients.

Identification of new branch points and unconventional introns in Saccharomyces cerevisiae

Spliced messages constitute one-fourth of expressed mRNAs in the yeast Saccharomyces cerevisiae, and most mRNAs in metazoans. Splicing requires 5' splice site (5'SS), branch point (BP), and 3' splice site (3'SS) elements, but the role of the BP in splicing control is poorly understood because BP identification remains difficult. We developed a high-throughput method, Branch-seq, to map BPs and 5'SSs of isolated RNA lariats. Applied to S. cerevisiae, Branch-seq detected 76% of expressed, annotated BPs and identified a comparable number of novel BPs. We performed RNA-seq to confirm associated 3'SS locations, identifying some 200 novel splice junctions, including an AT-AC intron. We show that several yeast introns use two or even three different BPs, with effects on 3'SS choice, protein coding potential, or RNA stability, and identify novel introns whose splicing changes during meiosis or in response to stress. Together, these findings show unanticipated complexity of splicing in yeast.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

Optimal antisense target reducing INS intron 1 retention is adjacent to a parallel G quadruplex

Splice-switching oligonucleotides (SSOs) have been widely used to inhibit exon usage but antisense strategies that promote removal of entire introns to increase splicing-mediated gene expression have not been developed. Here we show reduction of INS intron 1 retention by SSOs that bind transcripts derived from a human haplotype expressing low levels of proinsulin. This haplotype is tagged by a polypyrimidine tract variant rs689 that decreases the efficiency of intron 1 splicing and increases the relative abundance of mRNAs with extended 5' untranslated region (5' UTR), which curtails translation. Co-expression of haplotype-specific reporter constructs with SSOs bound to splicing regulatory motifs and decoy splice sites in primary transcripts revealed a motif that significantly reduced intron 1-containing mRNAs. Using an antisense microwalk at a single nucleotide resolution, the optimal target was mapped to a splicing silencer containing two pseudoacceptor sites sandwiched between predicted RNA guanine (G) quadruplex structures. Circular dichroism spectroscopy and nuclear magnetic resonance of synthetic G-rich oligoribonucleotide tracts derived from this region showed formation of a stable parallel 2-quartet G-quadruplex on the 3' side of the antisense retention target and an equilibrium between quadruplexes and stable hairpin-loop structures bound by optimal SSOs. This region interacts with heterogeneous nuclear ribonucleoproteins F and H that may interfere with conformational transitions involving the antisense target. The SSO-assisted promotion of weak intron removal from the 5' UTR through competing noncanonical and canonical RNA structures may facilitate development of novel strategies to enhance gene expression.

First deep intronic mutation in the NOTCH3 gene in a family with late-onset CADASIL

CADASIL is the most prominent inherited form of vascular dementia. The main clinical features include migraine with aura, stroke, mood disturbances, and cognitive decline, with a mid-life (30s-60s) adult onset. Genetic testing is the gold standard for the diagnosis. CADASIL is caused mostly by missense mutations in the NOTCH3 gene, invariably involving a cysteine residue. Only a couple of splice site mutations have been reported. In a few pathologically defined patients, genetic mutations remain unidentified. We report a family with late-onset CADASIL phenotype carrying a novel intronic deletion in the NOTCH3 gene (c.341-26_24delAAC). Transcript analysis revealed a splicing alteration, with the complete intron 3 retention. The insertion was in-frame and encoded an extra 25 amino acids, including 1 cysteine. This is the first report of an aberrant splicing event of the NOTCH3 gene associated with a mutation far away from the canonical splice site. Our finding suggests that the assays used to evaluate splicing should be mandatory in the diagnostic setting of genetically undefined CADASIL cases.

Exome sequencing identifies a branch point variant in Aarskog-Scott syndrome

Aarskog-Scott syndrome (ASS) is a rare disorder with characteristic facial, skeletal, and genital abnormalities. Mutations in the FGD1 gene (Xp11.21) are responsible for ASS. However, mutation detection rates are low. Here, we report a family with ASS where conventional Sanger sequencing failed to detect a pathogenic change in FGD1. To identify the causative gene, we performed whole-exome sequencing in two patients. An initial analysis did not reveal a likely candidate gene. After relaxing our filtering criteria, accepting larger intronic segments, we unexpectedly identified a branch point (BP) variant in FGD1. Analysis of patient-derived RNA showed complete skipping of exon 13, leading to premature translation termination. The BP variant detected is one of very few reported so far proven to affect splicing. Our results show that besides digging deeper to reveal nonobvious variants, isolation and analysis of RNA provides a valuable but under-appreciated tool to resolve cases with unknown genetic defects.

Allele-specific recognition of the 3' splice site of INS intron 1

Genetic predisposition to type 1 diabetes (T1D) has been associated with a chromosome 11 locus centered on the proinsulin gene (INS) and with differential steady-state levels of INS RNA from T1D-predisposing and -protective haplotypes. Here, we show that the haplotype-specific expression is determined by INS variants that control the splicing efficiency of intron 1. The adenine allele at IVS1-6 (rs689), which rapidly expanded in modern humans, renders the 3' splice site of this intron more dependent on the auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF). This interaction required both zinc fingers of the 35-kD U2AF subunit (U2AF35) and was associated with repression of a competing 3' splice site in INS exon 2. Systematic mutagenesis of reporter constructs showed that intron 1 removal was facilitated by conserved guanosine-rich enhancers and identified additional splicing regulatory motifs in exon 2. Sequencing of intron 1 in primates revealed that relaxation of its 3' splice site in Hominidae coevolved with the introduction of a short upstream open reading frame, providing a more efficient coupled splicing and translation control. Depletion of SR proteins 9G8 and transformer-2 by RNA interference was associated with exon 2 skipping whereas depletion of SRp20 with increased representation of transcripts containing a cryptic 3' splice site in the last exon. Together, these findings reveal critical interactions underlying the allele-dependent INS expression and INS-mediated risk of T1D and suggest that the increased requirement for U2AF35 in higher primates may hinder thymic presentation of autoantigens encoded by transcripts with weak 3' splice sites.

Four exons of the serotonin receptor 4 gene are associated with multiple distant branch points

Splicing of vertebrate introns involves recognition of three consensus elements at the 3' end. The branch point (BP) and polypyrimidine tract (PPT) are usually located within 40 nucleotides (nt) of the 3' splice site (3' ss), AG, but can be much more distant. A characteristic of the region between distant BPs (dBPs) and the 3' ss is the absence of intervening AG dinucleotides, leading to its designation as the "AG exclusion zone" (AGEZ). The human HTR4 gene, which encodes serotonin receptor 4 and has been associated with schizophrenia, bipolar disease, and gastrointestinal disorders, has four exons with extensive AGEZs. We have mapped the BPs for HTR4 exons 3, 4, 5, and g generated by in vitro splicing, and validated them by mutagenesis in exon-trapping vectors. All exons used dBPs up to 273 nt upstream of the exon. Strikingly, exons 4 and 5 used combinations of both distant and conventionally located BPs, suggesting that successful splicing of these exons can occur by distinct pathways. Our results emphasize the importance for single nucleotide polymorphism resequencing projects to take account of potential dBPs, as the extended AGEZs are vulnerable to mutations that could affect splicing itself or regulation of alternative splicing.

Mutation mechanisms

A mutation is an event that produces heritable changes in the DNA. There are many different types of mutations, including point mutations (changes that imply loss, duplication, or alterations of small DNA segments, often involving a single or a few nucleotides) and major DNA changes (loss, duplication, or rearrangements of entire genes or of gene segments). This article reviews how different types of mutation may result in defective gene expression.

A KCNH2 branch point mutation causing aberrant splicing contributes to an explanation of genotype-negative long QT syndrome

BACKGROUND

Genetic screening of long QT syndrome (LQTS) fails to identify disease-causing mutations in about 30% of patients. So far, molecular screening has focused mainly on coding sequence mutations or on substitutions at canonical splice sites.

OBJECTIVE

The purpose of this study was to explore the possibility that intronic variants not at canonical splice sites might affect splicing regulatory elements, lead to aberrant transcripts, and cause LQTS.

METHOD

Molecular screening was performed through DHPLC and sequence analysis. The role of the intronic mutation identified was assessed with a hybrid minigene splicing assay.

RESULTS

A three-generation LQTS family was investigated. Molecular screening failed to identify an obvious disease-causing mutation in the coding sequences of the major LQTS genes but revealed an intronic A-to-G substitution in KCNH2 (IVS9-28A/G) cosegregating with the clinical phenotype in family members. In vitro analysis proved that the mutation disrupts the acceptor splice site definition by affecting the branch point (BP) sequence and promoting intron retention. We further demonstrated a tight functional relationship between the BP and the polypyrimidine tract, whose weakness is responsible for the pathological effect of the IVS9-28A/G mutation.

CONCLUSIONS

We identified a novel BP mutation in KCNH2 that disrupts the intron 9 acceptor splice site definition and causes LQT2. The present finding demonstrates that intronic mutations affecting pre-mRNA processing may contribute to the failure of traditional molecular screening in identifying disease-causing mutations in LQTS subjects and offers a rationale strategy for the reduction of genotype-negative cases.

Genome-wide survey of allele-specific splicing in humans

Background

Accurate mRNA splicing depends on multiple regulatory signals encoded in the transcribed RNA sequence. Many examples of mutations within human splice regulatory regions that alter splicing qualitatively or quantitatively have been reported and allelic differences in mRNA splicing are likely to be a common and important source of phenotypic diversity at the molecular level, in addition to their contribution to genetic disease susceptibility. However, because the effect of a mutation on the efficiency of mRNA splicing is often difficult to predict, many mutations that cause disease through an effect on splicing are likely to remain undiscovered.

Results

We have combined a genome-wide scan for sequence polymorphisms likely to affect mRNA splicing with analysis of publicly available Expressed Sequence Tag (EST) and exon array data. The genome-wide scan uses published tools and identified 30,977 SNPs located within donor and acceptor splice sites, branch points and exonic splicing enhancer elements. For 1,185 candidate splicing polymorphisms the difference in splicing between alternative alleles was corroborated by publicly available exon array data from 166 lymphoblastoid cell lines. We developed a novel probabilistic method to infer allele-specific splicing from EST data. The method uses SNPs and alternative mRNA isoforms mapped to EST sequences and models both regulated alternative splicing as well as allele-specific splicing. We have also estimated heritability of splicing and report that a greater proportion of genes show evidence of splicing heritability than show heritability of overall gene expression level. Our results provide an extensive resource that can be used to assess the possible effect on splicing of human polymorphisms in putative splice-regulatory sites.

Conclusion

We report a set of genes showing evidence of allele-specific splicing from an integrated analysis of genomic polymorphisms, EST data and exon array data, including several examples for which there is experimental evidence of polymorphisms affecting splicing in the literature. We also present a set of novel allele-specific splicing candidates and discuss the strengths and weaknesses of alternative technologies for inferring the effect of sequence variants on mRNA splicing.

Patterns of missplicing caused by RB1 gene mutations in patients with retinoblastoma and association with phenotypic expression

We have analyzed RNA from retinoblastoma patients and unaffected carriers with various RB1 gene mutations to determine the patterns of missplicing and associations with phenotypic expression. Most sequence alterations in or in the neighborhood of conserved splice signals that we tested resulted in simple exon skipping (15 mutations) or intron inclusion (new acceptor AG-sites, four mutations) as expected. Two mutations resulted in skipping of a neighboring exon (exon 11), a complex pattern indicating competition for correct lariat formation. We observed no activation of a cryptic splice site but found that a recurrent missense mutation in exon 7 creates a new splice site (two families). RT-PCR analysis enabled us to confirm the presence and to characterize the transcriptional consequences of gross insertions and deletions in the RB1 gene in six patients, including two patients with mutational mosaicism. We also used RT-PCR analysis to search for unknown mutations in 15 patients and identified three oncogenic point mutations deep in introns. Two of these mutations are recurrent thus indicating that, despite the vast extent of the introns of the RB1 gene, few bases are effective targets for oncogenic mutations. When analyzing associations between phenotypic expression (16 families) and mutational consequences we observed no link to the presence or absence of a premature termination codon in the mutant transcript. However, the location of a mutation relative to the splice sequence has a strong and consistent influence on phenotypic expression.

Human branch point consensus sequence is yUnAy

Yeast carries a strictly conserved branch point sequence (BPS) of UACUAAC, whereas the human BPS is degenerative and is less well characterized. The human consensus BPS has never been extensively explored in vitro to date. Here, we sequenced 367 clones of lariat RT-PCR products arising from 52 introns of 20 human housekeeping genes. Among the 367 clones, a misincorporated nucleotide at the branch point was observed in 181 clones, for which we can precisely pinpoint the branch point. The branch points were comprised of 92.3% A, 3.3% C, 1.7% G and 2.8% U. Our analysis revealed that the human consensus BPS is simply yUnAy, where the underlined is the branch point at position zero and the lowercase pyrimidines ('y') are not as well conserved as the uppercase U and A. We found that the branch points are located 21-34 nucleotides upstream of the 3' end of an intron in 83% clones. We also found that the polypyrimidine tract spans 4-24 nucleotides downstream of the branch point. Our analysis demonstrates that the human BPSs are more degenerative than we have expected and that the human BPSs are likely to be recognized in combination with the polypyrimidine tract and/or the other splicing cis-elements.

Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Alternative splicing involving the 3' tandem splice site NAGNAG sequence may play a role in the structure-function diversity of proteins. However, how 3' tandem splice site utilization is determined is not well understood. We previously demonstrated that 3' NAGNAG-based wobble splicing occurs mostly in a tissue- and developmental stage-independent manner. Bioinformatic analysis reveals that the nucleotide preceding the AG dinucleotide may influence 3' splice site utilization; this is also supported by an in vivo splicing assay. Moreover, we found that the intron sequence plays an important role in 3' splice site selection for NAGNAG wobble splicing. Mutations of the region between the branch site and the NAGNAG 3' splice site, indeed, affected the ratio of the distal/proximal AG selection. Finally, we found that single nucleotide polymorphisms around the NAGNAG motif could affect the splice site choice, which may lead to a change in mRNA patterns and influence protein function. We conclude that the NAGNAG motif and its upstream region to the branch point sequence are required for 3' tandem splice site selection.