Role of the 3' splice site consensus sequence in mammalian pre-mRNA splicing

Purifying selection against spurious splicing signals contributes to the base composition evolution of the polypyrimidine tract

Among eukaryotes, the major spliceosomal pathway is highly conserved. While long introns may contain additional regulatory sequences, the ones in short introns seem to be nearly exclusively related to splicing. Although these regulatory sequences involved in splicing are well-characterized, little is known about their evolution. At the 3' end of introns, the splice signal nearly universally contains the dimer AG, which consists of purines, and the polypyrimidine tract upstream of this 3' splice signal is characterized by over-representation of pyrimidines. If the over-representation of pyrimidines in the polypyrimidine tract is also due to avoidance of a premature splicing signal, we hypothesize that AG should be the most under-represented dimer. Through the use of DNA-strand asymmetry patterns, we confirm this prediction in fruit flies of the genus Drosophila and by comparing the asymmetry patterns to a presumably neutrally evolving region, we quantify the selection strength acting on each motif. Moreover, our inference and simulation method revealed that the best explanation for the base composition evolution of the polypyrimidine tract is the joint action of purifying selection against a spurious 3' splice signal and the selection for pyrimidines. Patterns of asymmetry in other eukaryotes indicate that avoidance of premature splicing similarly affects the nucleotide composition in their polypyrimidine tracts.

Regulation of Arp5 expression by alternative splicing coupled to nonsense-mediated RNA decay

Actin-related protein 5 (ARP5) inhibits the differentiation of skeletal, smooth, and cardiac muscle tissues, and ARP5 expression increases or decreases according to physiological and pathological changes in the muscle differentiation status. However, the regulatory mechanisms of ARP5 expression are largely unknown. Here, we identified a novel Arp5 mRNA isoform that contains premature termination codons in alternative exon 7b and is thus targeted by nonsense-mediated mRNA decay (NMD). In mouse skeletal muscle cells, switching from the canonical Arp5 isoform, i.e., Arp5(7a), to the NMD-targeted isoform Arp5(7b) occurred during differentiation, suggesting that Arp5 expression is regulated by alternative splicing coupled to NMD (AS-NMD). We developed an original method to accurately quantify the proportion of both Arp5 isoforms and measured higher levels of Arp5(7b) in muscle and brain tissues, where ARP5 is less expressed. The 3' splice site in Arp5 exon 7 has an unusual acceptor sequence that often leads to the skip of the authentic splice site and the use of the cryptic splice site localized 16 bases downstream. When the unusual acceptor sequence was mutated to the usual one, the Arp5(7b) isoform was barely detectable. The expression of several splicing factors involved in 3' splice site recognition was reduced after muscle differentiation. Additionally, knockdown of splicing factors increased the levels of Arp5(7b) and decreased the expression of Arp5(7a). Furthermore, strong positive correlations were found between Arp5 expression and the levels of these splicing factors in human skeletal and cardiac muscle tissues. Thus, Arp5 expression in muscle tissues is most likely regulated by the AS-NMD pathway.

An LMNA synonymous variant associated with severe dilated cardiomyopathy: Case report

Dilated cardiomyopathy (DCM) is one of the most common cardiac phenotypes caused by mutations of lamin A/C (LMNA) gene in humans. In our study, a cohort of 57 patients who underwent heart transplant for dilated cardiomyopathy was screened for variants in LMNA. We identified a synonymous variant c.936G>A in the last nucleotide of exon 5 of LMNA in a DCM family. Clinically, the LMNA variant carriers presented with severe familial DCM, conduction disease, and high creatine-kinase level. The LMNA c.936G>A variant is novel and has not been reported in current genetic variant databases. Sanger sequencing results showed the presence of LMNA c.936G>A variant in the genomic DNA but not in the cDNA derived from one family member's heart tissue. Real-time quantitative polymerase chain reaction showed significantly lower LMNA mRNA levels in the patient's heart compared to the controls, suggesting that the c.936G>A LMNA variant resulted in reduced mRNA and possibly lower protein expression of LMNA. These findings expand the understanding on the association between synonymous variant of LMNA and the molecular pathogenesis in DCM patients.

Dynamic stacking of an expected branch point adenosine in duplexes containing pseudouridine-modified or unmodified U2 snRNA sites

The pre-mRNA branch point sequence (BPS) anneals with a pseudouridine-modified region of the U2 small nuclear (sn)RNA, and offers a 2' hydroxyl group of a bulged adenosine as the nucleophile for the first catalytic step of pre-mRNA splicing. To increase our structural understanding of branch site selection, we characterized a duplex containing a BPS sequence that is common among multicellular eukaryotes (5'-UACUGAC-3') and the complementary U2 snRNA site using NMR. A major conformation of the expected branch site adenosine stacked within the duplex and paired with the conserved pseudouridine of the U2 snRNA strand. In contrast, the guanosine preceding the branch site appeared flexible and had weak contacts with the surrounding nucleotides. Pseudouridine-modified and unmodified U2 snRNA-BPS-containing duplexes remained structurally similar. These results highlight the importance of auxiliary factors to achieve the active bulged conformation of the branch site nucleophile for the first step of pre-mRNA splicing.

BPP: a sequence-based algorithm for branch point prediction

Motivation

Although high-throughput sequencing methods have been proposed to identify splicing branch points in the human genome, these methods can only detect a small fraction of the branch points subject to the sequencing depth, experimental cost and the expression level of the mRNA. An accurate computational model for branch point prediction is therefore an ongoing objective in human genome research.

Results

We here propose a novel branch point prediction algorithm that utilizes information on the branch point sequence and the polypyrimidine tract. Using experimentally validated data, we demonstrate that our proposed method outperforms existing methods. Availability and implementation: https://github.com/zhqingit/BPP.

Contact

djguo@cuhk.edu.hk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Splicing and transcription touch base: co-transcriptional spliceosome assembly and function

Several macromolecular machines collaborate to produce eukaryotic messenger RNA. RNA polymerase II (Pol II) translocates along genes that are up to millions of base pairs in length and generates a flexible RNA copy of the DNA template. This nascent RNA harbours introns that are removed by the spliceosome, which is a megadalton ribonucleoprotein complex that positions the distant ends of the intron into its catalytic centre. Emerging evidence that the catalytic spliceosome is physically close to Pol II in vivo implies that transcription and splicing occur on similar timescales and that the transcription and splicing machineries may be spatially constrained. In this Review, we discuss aspects of spliceosome assembly, transcription elongation and other co-transcriptional events that allow the temporal coordination of co-transcriptional splicing.

Identification of human short introns

Canonical pre-mRNA splicing requires snRNPs and associated splicing factors to excise conserved intronic sequences, with a minimum intron length required for efficient splicing. Non-canonical splicing-intron excision without the spliceosome-has been documented; most notably, some tRNAs and the XBP1 mRNA contain short introns that are not removed by the spliceosome. There have been some efforts to identify additional short introns, but little is known about how many short introns are processed from mRNAs. Here, we report an approach to identify RNA short introns from RNA-Seq data, discriminating against small genomic deletions. We identify hundreds of short introns conserved among multiple human cell lines. These short introns are often alternatively spliced and are found in a variety of RNAs-both mRNAs and lncRNAs. Short intron splicing efficiency is increased by secondary structure, and we detect both canonical and non-canonical short introns. In many cases, splicing of these short introns from mRNAs is predicted to alter the reading frame and change protein output. Our findings imply that standard gene prediction models which often assume a lower limit for intron size fail to predict short introns effectively. We conclude that short introns are abundant in the human transcriptome, and short intron splicing represents an added layer to mRNA regulation.

Validation of three splice donor and three splice acceptor sites for regulating four novel low-abundance spliced transcripts of human cytomegalovirus UL21.5 gene locus

In a previous study, one spliced transcript of human cytomegalovirus (HCMV), named UL21.5 was identified. UL21.5 has been found to be one of the viral transcripts packaged within HCMV particles. The UL21.5 mRNA is translated into a secreted glycoprotein, which is a viral chemokine decoy receptor specifically interacting with regulated upon activation normal T cell expressed and secreted (RANTES). In the present study, four novel low-abundance 3'-coterminal spliced transcripts were identified to be transcribed from the UL21.5 gene region of a low-passage HCMV strain during the late infection phase by cDNA library screening, northern blot hybridization, reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE)-PCR. Three splicing donor and three splicing acceptor sites found in the UL21.5 gene region were validated to be functional in an in vitro expression system. In addition, the determinant regulatory region that is necessary for the splice donor site at nucleotide (nt) 25533 was located in a 9-bp sequence around the site; the regulatory regions for the splice acceptor sites at nt 26597 and nt 26633 were located in a 20-bp sequence upstream of the site at nt 26597 and in a 10-bp sequence from nt 26641 to nt 26650 downstream of the site at nt 26633, respectively.

Origin of spliceosomal introns and alternative splicing

In this work we review the current knowledge on the prehistory, origins, and evolution of spliceosomal introns. First, we briefly outline the major features of the different types of introns, with particular emphasis on the nonspliceosomal self-splicing group II introns, which are widely thought to be the ancestors of spliceosomal introns. Next, we discuss the main scenarios proposed for the origin and proliferation of spliceosomal introns, an event intimately linked to eukaryogenesis. We then summarize the evidence that suggests that the last eukaryotic common ancestor (LECA) had remarkably high intron densities and many associated characteristics resembling modern intron-rich genomes. From this intron-rich LECA, the different eukaryotic lineages have taken very distinct evolutionary paths leading to profoundly diverged modern genome structures. Finally, we discuss the origins of alternative splicing and the qualitative differences in alternative splicing forms and functions across lineages.

U2AF65 adapts to diverse pre-mRNA splice sites through conformational selection of specific and promiscuous RNA recognition motifs

Degenerate splice site sequences mark the intron boundaries of pre-mRNA transcripts in multicellular eukaryotes. The essential pre-mRNA splicing factor U2AF⁶⁵ is faced with the paradoxical tasks of accurately targeting polypyrimidine (Py) tracts preceding 3′ splice sites while adapting to both cytidine and uridine nucleotides with nearly equivalent frequencies. To understand how U2AF⁶⁵ recognizes degenerate Py tracts, we determined six crystal structures of human U2AF⁶⁵ bound to cytidine-containing Py tracts. As deoxy-ribose backbones were required for co-crystallization with these Py tracts, we also determined two baseline structures of U2AF⁶⁵ bound to the deoxy-uridine counterparts and compared the original, RNA-bound structure. Local structural changes suggest that the N-terminal RNA recognition motif 1 (RRM1) is more promiscuous for cytosine-containing Py tracts than the C-terminal RRM2. These structural differences between the RRMs were reinforced by the specificities of wild-type and site-directed mutant U2AF⁶⁵ for region-dependent cytosine- and uracil-containing RNA sites. Small-angle X-ray scattering analyses further demonstrated that Py tract variations select distinct inter-RRM spacings from a pre-existing ensemble of U2AF⁶⁵ conformations. Our results highlight both local and global conformational selection as a means for universal 3′ splice site recognition by U2AF⁶⁵.

Establishment of an in vitro trans-splicing system in Trypanosoma brucei that requires endogenous spliced leader RNA

In trypanosomes a 39 nucleotide exon, the spliced leader (SL) is donated to all mRNAs from a small RNA, the SL RNA, by trans-splicing. Since the discovery of trans-splicing in trypanosomes two decades ago, numerous attempts failed to reconstitute the reaction in vitro. In this study, a crude whole-cell extract utilizing the endogenous SL RNA and synthetic tubulin pre-mRNA were used to reconstitute the trans-splicing reaction. An RNase protection assay was used to detect the trans-spliced product. The reaction was optimized and shown to depend on ATP and intact U2 and U6 snRNPs. Mutations introduced at the polypyrimidine tract and the AG splice site reduced the reaction efficiency. To simplify the assay, RT–PCR and quantitative real-time PCR assays were established. The system was used to examine the structural requirements for SL RNA as a substrate in the reaction. Interestingly, synthetic SL RNA assembled poorly to its cognate particle and was not utilized in the reaction. However, SL RNA synthesized in cells lacking Sm proteins, which is defective in cap-4 modification, was active in the reaction. This study is the first step towards further elucidating the mechanism of trans-splicing, an essential reaction which determines the trypanosome transcriptome.

A 5'UTR-spliced mRNA isoform is specialized for enhanced HIV-2 gag translation

Full-length unspliced genomic RNA plays critical roles in HIV replication, serving both as mRNA for the synthesis of the key viral polyproteins Gag and Gag-Pol and as genomic RNA for encapsidation into assembling viral particles. We show that a second gag mRNA species that differs from the genomic RNA molecule by the absence of an intron in the 5' untranslated region (5'UTR) is produced during HIV-2 replication in cell culture and in infected patients. We developed a cotransfection system in which epitopically tagged Gag proteins can be traced back to their mRNA origins in the translation pool. We show that a disproportionate amount of Gag is translated from 5'UTR intron-spliced mRNAs, demonstrating a role for the 5'UTR intron in the regulation of gag translation. To further characterize the effects of the HIV-2 5'UTR on translation, we fused wild-type, spliced, or mutant leader RNA constructs to a luciferase reporter gene and assayed their translation in reticulocyte lysates. These assays confirmed that leaders lacking the 5'UTR intron increased translational efficiency compared to that of the unspliced leader. In addition, we found that removal or mutagenesis of the C-box, a pyrimidine-rich sequence located in the 5'UTR intron and previously shown to affect RNA dimerization, also strongly influenced translational efficiency. These results suggest that the splicing of both the 5'UTR intron and the C-box element have key roles in regulation of HIV-2 gag translation in vitro and in vivo.

The cathepsin L first intron stimulates gene expression in rat sertoli cells

Large amounts of cathepsin L (CTSL), a cysteine protease required for quantitatively normal spermatogenesis, are synthesized by mouse and rat Sertoli cells during stages VI to VII of the cycle of the seminiferous epithelium. We previously demonstrated that all of the regulatory elements required in vivo for both Sertoli cell- and stage-specific expression of the Ctsl gene are present within a ~3-kb genomic fragment that contains 2065 nucleotides upstream of the transcription start site and 977 nucleotides of downstream sequence. Most of the downstream region encodes the first intron. In this study, transient transfection assays using primary Sertoli cell cultures and the TM4 Sertoli cell line established that the Ctsl first intron increased reporter gene activity by ~5-fold. While the intron-mediated enhancement in reporter gene activity was not restricted to the Ctsl promoter, positioning the first intron upstream of the Ctsl promoter in either orientation abolished its stimulatory activity, suggesting that it does not contain a typical enhancer. Mutating the 5'-splice site of the Ctsl first intron or replacing the first intron by the Ctsl fourth intron abolished the stimulatory effect. Finally, the intron-dependent increase in reporter gene activity could be explained in part by an increase in the amounts of total RNA and transcript polyadenylation. Results from this study suggest that the stimulatory effect mediated by the Ctsl first intron may explain in part why Sertoli cells in seminiferous tubules at stages VI to VII produce high levels of CTSL.

Looking for organization patterns of highly expressed genes: purine-pyrimidine composition of precursor mRNAs

We analyzed precursor messenger RNAs (pre-mRNAs) of 12 eukaryotic species. In each species, three groups of highly expressed genes, ribosomal proteins, heat shock proteins, and amino-acyl tRNA synthetases, were compared with a control group (randomly selected genes). The purine-pyrimidine (R-Y) composition of pre-mRNAs of the three targeted gene groups proved to differ significantly from the control. The exons of the three groups tested have higher purine contents and R-tract abundance and lower abundance of Y-tracts compared to the control (R-tract-tract of sequential purines with Rn>or=5; Y-tract-tract of sequential pyrimidines with Yn>or=5). In species widely employing "intron definition" in the splicing process, the Y content of introns of the three targeted groups appeared to be higher compared to the control group. Furthermore, in all examined species, the introns of the targeted genes have a lower abundance of R-tracts compared to the control. We hypothesized that the R-Y composition of the targeted gene groups contributes to high rate and efficiency of both splicing and translation, in addition to the mRNA coding role. This is presumably achieved by (1) reducing the possibility of the formation of secondary structures in the mRNA, (2) using the R-tracts and R-biased sequences as exonic splicing enhancers, (3) lowering the amount of targets for pyrimidine tract binding protein in the exons, and (4) reducing the amount of target sequences for binding of serine/arginine-rich (SR) proteins in the introns, thereby allowing SR proteins to bind to proper (exonic) targets.

Polypyrimidine tract binding protein modulates efficiency of polyadenylation

Polypyrimidine tract binding protein (PTB) is a major hnRNP protein with multiple roles in mRNA metabolism, including regulation of alternative splicing and internal ribosome entry site-driven translation. We show here that a fourfold overexpression of PTB results in a 75% reduction of mRNA levels produced from transfected gene constructs with different polyadenylation signals (pA signals). This effect is due to the reduced efficiency of mRNA 3' end cleavage, and in vitro analysis reveals that PTB competes with CstF for recognition of the pA signal's pyrimidine-rich downstream sequence element. This may be analogous to its role in alternative splicing, where PTB competes with U2AF for binding to pyrimidine-rich intronic sequences. The pA signal of the C2 complement gene unusually possesses a PTB-dependent upstream sequence, so that knockdown of PTB expression by RNA interference reduces C2 mRNA expression even though PTB overexpression still inhibits polyadenylation. Consequently, we show that PTB can act as a regulator of mRNA expression through both its negative and positive effects on mRNA 3' end processing.

SNTG1, the gene encoding gamma1-syntrophin: a candidate gene for idiopathic scoliosis

Idiopathic scoliosis (IS) affects approximately 2%-3% of the population and has a heritable component. The genetics of this disorder are complex. Here, we describe a family in which a pericentric inversion of chromosome 8 co-segregates with IS. We have used fluorescence in situ hybridization to identify cloned DNAs that span the breakpoints on the two arms of the chromosome. We have identified a bacterial artificial chromosome (BAC) of 150 kb that crosses the q-arm breakpoint and a BAC of 120 kb that crosses the p-arm breakpoint. The complete genomic DNA sequence of these BACs has been analyzed to identify candidate genes and to localize further the precise breakpoints. This has revealed that the p-arm break does not interrupt any known gene and occurs in a region of highly repetitive sequence elements. On the q-arm, the break occurs between exons 10 and 11 of the gamma-1 syntrophin (SNTG1) gene. Syntrophins are a group of cytoplasmic peripheral membrane proteins that associate directly with dystrophin, the Duchenne muscular dystrophy gene; gamma1-syntrophin has been shown to be a neuronal cell-specific protein. Mutational analysis of SNTG1 exons in 152 sporadic IS patients has revealed a 6-bp deletion in exon 10 of SNTG1 in one patient and a 2-bp insertion/deletion mutation occurring in a polypyrimidine tract of intronic sequence 20 bases upstream of the SNTG1 exon 5 splice site in two patients. These changes were not seen in a screen of 480 control chromosomes. Genomic DNAs from seven affected individuals within the family of a patient carrying the 6-bp deletion were typed to determine whether the alteration co-segregated with IS. The deletion was only observed in five out of these seven individuals. Thus, although genetic heterogeneity or multiple alleles cannot be ruled out, the 6-bp deletion does not consistently co-segregate with the disease in this family.

Human cytomegalovirus UL37 immediate early target minigene RNAs are accurately spliced and polyadenylated

The human cytomegalovirus (HCMV) UL36-38 immediate early (IE) locus encodes proteins required for virus growth. The UL37 IE promoter drives production of differentially spliced and unspliced RNAs. To study their post-transcriptional processing, we generated target minigenes encoding each UL37 RNA splicing substrate. Target 1 RNA, spanning UL37 exon 1 (x1) donor and 2 (x2) acceptor as well as adjacent intronic sequences, but not the UL38 gene, accurately reproduced UL37 x1/x2 RNA splicing in transfected permissive cells. Surprisingly, deletion of distal intronic sequences nt -82 to -143 from the UL37x2 acceptor resulted in aberrant splicing to an upstream non-consensus exonic donor. Target 1 RNAs carry the UL37x1 polyadenylation (PA) signal and site as well as a downstream SV40 early PA signal. Both the UL37x1 and SV40 PA signals are used in wild-type target 1 RNAs but inhibited in UL37x1 PA signal mutants. Alternative RNA splicing of UL37 exons 2 to 3 or 3A as well as exons 3 to 4, observed in HCMV mature UL37 and UL36 spliced RNAs, is accurately reproduced with target minigene RNAs carrying the corresponding UL37 exonic and intronic sequences. Moreover, alternative splicing using two novel UL37 exon 3 consensus splice donors (di and dii) was found in target and in HCMV-infected cell RNA. These results demonstrate that: (i) target minigene RNAs accurately recapitulate the processing of UL37 IE RNAs in the HCMV-infected cell; (ii) precise UL37x1 donor selection is modulated by 3'-distal UL37 intronic sequences; and (iii) UL37 exon 3 contains multiple alternative consensus splice donors.

Mutation analysis of SYNJ1: a possible candidate gene for chromosome 21q22-linked bipolar disorder

Genes involved in the regulation of synaptic vesicle function are potential candidates for the development of psychiatric disorders. In addition to experimental and theoretical considerations, a number of genes involved in synaptic vesicle function map to regions of the genome that have been linked to bipolar disorder (BPD) and schizophrenia (SZ). One is synaptojanin 1 (SYNJ1) which maps to 21q22.2, a chromosomal region that has been linked to BPD in a subset of families in several studies. Synaptojanin 1 is an inositol 5-phosphatase that has an important role in synaptic vesicle endocytosis. Mutation screening of 32 exons, intron--exon junctions, and 839 bases of 5'-flanking DNA resulted in the identification of 11 mutations of which four were very common and seven were very rare. Of the 11 mutations identified, several may have functional significance including two coding variants, two that may affect the binding of a transcription factor, and two that involve known splicing regulatory domains. Five bipolar patients out of 149 analyzed were found who have one of the four rare variants that were most likely to have functional significance compared with 0/148 controls. The allele frequencies for three of the four common variants were very similar in bipolar patients and controls. A slight difference in allele frequency was found for an interesting mutation we detected in intron 12 in which two non-adjacent thymidine residues are deleted in a poly-AT tract located near the exon 12 splice donor site (chi(2) = 2.45, P = 0.12, 2-tailed). Although we failed to unequivocally identify a specific SYNJ1 allele that could be responsible for putative chromosome 21q22-linked BPD, several interesting variants were found to be increased in bipolar subjects and should be further investigated.

Identification of an alternative transcript of ABCA1 gene in different human cell types

We have observed two ABCA1 gene transcripts in human skin fibroblasts. The RT-PCR amplification of the exon 3-exon 8 region generated a normal fragment (740 bp) and an abnormal fragment (600 bp) in a ratio ranging from 3:1 to 8/9:1. These two transcripts were present in other cells such as leukemia T-cells, endothelial and smooth muscle cells as well human hepatoma cells (HepG2). Restriction enzyme analysis and sequencing indicated that in the abnormal fragment exon 3 was followed by exon 5. The complete skipping of exon 4 leads to a premature stop and a predicted translation product of 74 amino acids. The ratio between the normal and alternative transcript is not affected by variation in ABCA1 gene expression induced by incubating cells in serum-free medium and in the presence of cholesterol. It is possible that this alternative splicing represents as mechanism that regulates the ABCA1 content in tissues.

Role of the 3' splice site in U12-dependent intron splicing

U12-dependent introns containing alterations of the 3' splice site AC dinucleotide or alterations in the spacing between the branch site and the 3' splice site were examined for their effects on splice site selection in vivo and in vitro. Using an intron with a 5' splice site AU dinucleotide, any nucleotide could serve as the 3'-terminal nucleotide, although a C residue was most active, while a U residue was least active. The penultimate A residue, by contrast, was essential for 3' splice site function. A branch site-to-3' splice site spacing of less than 10 or more than 20 nucleotides strongly activated alternative 3' splice sites. A strong preference for a spacing of about 12 nucleotides was observed. The combined in vivo and in vitro results suggest that the branch site is recognized in the absence of an active 3' splice site but that formation of the prespliceosomal complex A requires an active 3' splice site. Furthermore, the U12-type spliceosome appears to be unable to scan for a distal 3' splice site.

Identification of a novel PTEN intronic deletion in Li-Fraumeni syndrome and its effect on RNA processing

Germline mutations of the TP53 tumor suppressor gene account for the predisposition to cancer observed in many Li-Fraumeni syndrome (LFS) families. A causative genetic factor in the remaining families that harbor no TP53 mutations remains to be elucidated. The PTEN phosphatase tumor suppressor gene is mutated in human cancers observed in LFS. There also exists some phenotypic overlap in the occurrence of cancers in LFS and Cowden's disease (CD), for which germline PTEN mutations are believed to be responsible. We hypothesized that PTEN may be altered in some TP53-wild-type LFS families. We examined LFS primary patient lymphocytes for PTEN alterations using SSCP and sequence analysis. A novel intronic deletion was found in two unrelated individuals, adjacent to the splice acceptor site of PTEN exon 4. Based on an in vitro mRNA processing assay this alteration is predicted to be a polymorphism. The in vivo effects of this proximal splice site deletion are unknown and a genetic cause for the cancers in these families remains to be elucidated. Germline mutations of PTEN were not detected in other families, suggesting that alterations of this tumor suppressor gene do not account for the cancers observed in the subset of LFS individuals with wild-type germline TP53.

Analysis of the pseudoautosomal X-linked gene SYBL1in bipolar affective disorder: description of a new candidate allele for psychiatric disorders

The absence of father-to-son transmission has been observed in a subset of families with bipolar disorder (BPD), suggestive of a susceptibility gene on the sex-linked portion of the X chromosome. This is supported by some genetic linkage studies that have provided evidence for a susceptibility locus near Xq28. We have analyzed one candidate gene on Xq28, SYBL1, which maps to the Xq pseudoautosomal region (PAR). SYBL1 encodes a member of the synaptobrevin family of proteins that is involved in synaptic vesicle docking and membrane transport. Genes in the PAR generally escape X-chromosome inactivation and have an active homolog on the Y chromosome, which would result in an increase in same-sex concordance in paternal transmitted traits. However, SYBL1 is neither expressed on the Y chromosome nor the inactive X chromosome and would therefore be expected to show typical sex-linked transmission. We have screened SYBL1 for mutations that could be tested as candidate alleles in the development of BPD. Following single-strand conformation polymorphism (SSCP) analysis and DNA sequencing, four single nucleotide polymorphisms were detected: a silent mutation at codon 108, two intron mutations without any obvious biological significance, and a G-->C transversion in the polypyrimidine tract at the 3' splice acceptor site preceding exon 8. This polymorphism, which creates a perfect 16/16 stretch of pyrimidines, was analyzed in 110 patients with BPD not selected for sex-linked transmission and 119 control subjects. The results show a statistical trend toward an increase in the frequency of the C allele in males with BPD but not females. Males: chi(2) = 3.46, 1 df, p =.06; Females: chi(2) =.20, 1 df, p =.66.

Beta-thalassaemia in Cubans: novel allele increases the genetic diversity at the HBB locus in the Caribbean

In order to establish the molecular basis of beta-thalassaemia in Cubans, a total of 75 unrelated individuals, with beta-thalassaemia major (7), Hb S-beta-thalassaemia (28), Hb C-beta-thalassaemia (1), and beta-thalassaemia trait (39) yielding 82 beta-thalassaemia alleles, were analyzed. Seventeen different point mutations were identified accounting for 93% of the beta-thalassaemia alleles studied, revealing a high genetic heterogeneity at the HBB locus in this population. The more prevalent mutations, namely, CD 39 (C --> T) (30.5%), -29 (A --> G) (13.4%), IVS-I-110 (G --> A) (8.5%), and IVS-II-1 (G --> A) (8.5%), reflect the Mediterranean and African predominant ancestry of the extant Cuban population. We also report the identification of a novel allele, IVS-I-108 (T --> C), that possibly activates a cryptic branch site, in a beta-thalassaemia carrier with no other molecular defect within the beta-globin gene and its proximal promoter. This study shows that prenatal diagnosis of beta-thalassaemia should be feasible in about 60% of at-risk pregnancies by direct detection of selected point mutations. However, due to the wide spectrum of mutations, and in order to offer fully informative prenatal diagnosis to more than 87% of at-risk couples, the screening for beta-thalassaemia mutations in Cubans should be performed by using a general point mutation detection method, such as DGGE (denaturing gradient gel electrophoresis).

Role of adenine functional groups in the recognition of the 3'-splice-site AG during the second step of pre-mRNA splicing

The AG dinucleotide at the 3' splice sites of metazoan nuclear pre-mRNAs plays a critical role in catalytic step II of the splicing reaction. Previous studies have shown that replacement of the guanine by adenine in the AG (AG --> GG) inhibits this step. We find that the second step was even more severely inhibited by cytosine (AG --> CG) or uracil (AG --> UG) substitutions at this position. By contrast, a relatively moderate inhibition was observed with a hypoxanthine substitution (AG --> HG). When adenine was replaced by a purine base (AG --> PG) or by 7-deazaadenine (AG --> c(7)AG), little effect on the second step was observed, suggesting that the 6-NH(2) and N(7) groups do not play a critical role in adenine recognition. Finally, replacement of adenine by 2-aminopurine (AG --> 2-APG) had no effect on the second step. Taken together, our results suggest that the N(1) group of adenine functions as an essential determinant in adenine recognition during the second step of pre-mRNA splicing.

Allele-specific genetic interactions between Prp8 and RNA active site residues suggest a function for Prp8 at the catalytic core of the spliceosome

The highly conserved spliceosomal protein Prp8 is known to cross-link the critical sequences at both the 5' (GU) and 3' (YAG) ends of the intron. We have identified prp8 mutants with the remarkable property of suppressing exon ligation defects due to mutations in position 2 of the 5' GU, and all positions of the 3' YAG. The prp8 mutants also suppress mutations in position A51 of the critical ACAGAG motif in U6 snRNA, which has been observed previously to cross-link position 2 of the 5' GU. Other mutations in the 5' splice site, branchpoint, and neighboring residues of the U6 ACAGAG motif are not suppressed. Notably, the suppressed residues are specifically conserved from yeast to man, and from U2- to U12-dependent spliceosomes. We propose that Prp8 participates in a previously unrecognized tertiary interaction between U6 snRNA and both the 5' and 3' ends of the intron. This model suggests a mechanism for positioning the 3' splice site for catalysis, and assigns a fundamental role for Prp8 in pre-mRNA splicing.

Functional interactions of Prp8 with both splice sites at the spliceosomal catalytic center

A U5 snRNP protein, hPrp8, interacts closely with the GU dinucleotide at the 5' splice site (5'SS), forming a specific UV-inducible cross-link. To test if this physical contact between the 5'SS and the carboxy-terminal region of Prp8 reflects a functional recognition of the 5'SS during spliceosome assembly, we mutagenized the corresponding region of yeast Prp8 and screened the resulting mutants for suppression of 5'SS mutations in vivo. All of the isolated prp8 alleles not only suppress 5'SS but also 3'SS mutations, affecting the second catalytic step. Suppression of the 5'SS mutations by prp8 alleles was also tested in the presence of U1-7U snRNA, a predicted suppressor of the U+2A mutation. As expected, U1-7U efficiently suppresses prespliceosome formation, and the first, but not the second, step of U+2A pre-mRNA splicing. Independently, Prp8 functionally interacts with both splice sites at the later stage of splicing, affecting the efficiency of the second catalytic step. The striking proximity of two of the prp8 suppressor mutations to the site of the 5'SS:hPrp8 cross-link suggests that some protein:5'SS contacts made before the first step may be subsequently extended to accommodate the 3'SS for the second catalytic step. Together, these results strongly implicate Prp8 in specific interactions at the catalytic center of the spliceosome.

A cooperative interaction between U2AF65 and mBBP/SF1 facilitates branchpoint region recognition

During the early events of pre-mRNA splicing, intronic cis-acting sequences are recognized and interact through a network of RNA-RNA, RNA-protein, and protein-protein contacts. Recently, we identified a branchpoint sequence binding protein in yeast (BBP). The mammalian ortholog (mBBP/SF1) also binds specifically to branchpoint sequences and interacts with the well studied mammalian splicing factor U2AF65, which binds to the adjacent polypyrimidine (PY) tract. In this paper we demonstrate that the mBBP/SF1-U2AF65 interaction promotes cooperative binding to a branchpoint sequence-polypyrimidine tract-containing RNA, and we suggest that this cooperative RNA binding contributes to initial recognition of the branchpoint sequence (BPS) during pre-mRNA splicing. We also demonstrate the essential nature of the third RBD of U2AF65 for the interaction between the two proteins, both in the presence and absence of RNA.

Both the polypyrimidine tract and the 3' splice site function prior to the first step of splicing in fission yeast

While it is known that several trans -acting splicing factors are highly conserved between Schizosaccharomyces pombe and mammals, the roles of cis -acting signals have received comparatively little attention. In Saccharomyces cerevisiae, sequences downstream from the branch point are not required prior to the first transesterification reaction, whereas in mammals the polypyrimidine tract and, in some introns, the 3' AG dinucleotide are critical for initial recognition of an intron. We have investigated the contribution of these two sequence elements to splicing in S.pombe. To determine the stage at which the polypyrimidine tract functions, we analyzed the second intron of the cdc2 gene (cdc 2-Int2), in which pyrimidines span the entire interval between the branch point and 3' splice site. Our data indicate that substitution of a polypurine tract results in accumulation of linear pre-mRNA, while expanding the polypyrimidine tract enhances splicing efficiency, as in mammals. To examine the role of the AG dinucleotide in cdc 2-Int2 splicing, we mutated the 3' splice junction in both the wild-type and pyrimidine tract variant RNAs. These changes block the first transesterification reaction, as in a subset of mammalian introns. However, in contrast to the situation in mammals, we were unable to rescue the first step of splicing in a 3' splice site mutant by expanding the polypyrimidine tract. Mutating the terminal G in the third intron of the nda 3 gene (nda 3-Int3) also blocks the first transesterification reaction, suggesting that early recognition of the 3' splice site is a general property of fission yeast introns. Counter to earlier work with an artificial intron, it is not possible to restore the first step of splicing in cdc 2-Int2 and nda 3-Int3 3' splice site mutants by introducing compensatory changes in U1 snRNA. These results highlight the diversity and probable redundancy of mechanisms for identifying the 3' ends of introns.

Higher proportion of intact exon 9 CFTR mRNA in nasal epithelium compared with vas deferens

The 5-thymidine (5T) variant of the cystic fibrosis transmembrane conductance regulator (CFTR) intron 8 polypyrimidine tract (IVS8-T tract) is the most frequent CFTR gene alteration identified in men with congenital bilateral absence of vas deferens (CBAVD). This alternative splicing variant gives rise to two transcripts, one normal with exon 9 intact and the other with in-frame deletion of exon 9. That CBAVD men usually have none of the other clinical signs of classical cystic fibrosis (CF) suggests less functional CFTR is produced in the reproductive tract than in other CF-associated organs. Nasal epithelia and segments of vas deferens were obtained from healthy, previously vasectomized men who presented for vasectomy reversal. Quantitative RT-PCR was performed on these specimens, with the region of CFTR cDNA spanning exon 9 amplified. For both nasal and vasal tissues, a strong positive correlation was found between the length of the IVS8-T tract and the proportion of mRNA with exon 9 intact. In addition, within the same subject, a significantly higher level of transcripts lacking exon 9 was found in vas deferens than nasal epithelia, regardless of the IVS8-T genotype. These findings suggest that the splicing of CFTR precursor mRNA is less efficient in vasal epithelia compared with respiratory epithelia. Thus, differential splicing efficiency between the various tissues which express CFTR provides one possible explanation for the reproductive tract abnormalities observed in infertile men with CFTR gene alterations but without other clinical manifestations of CF.

Gene structure and functional analysis of the human Na+/phosphate co-transporter

Three lambda phage clones encompassing the Na+/phosphate co-transporter (NaPi-3) gene and its 5' flanking region were isolated from a human genomic DNA library. The gene comprises 13 exons and 12 introns and spans approx. 14 kb. All exon-intron junctions conform to the GT/AG rule. The major transcription-initiation site was determined by primer-extension analysis and is an adenosine residue 57 bp upstream of the 3' end of the first exon. There is a typical TATA box 28 bp upstream of the major transcription-initiation site and various cis-acting elements, including a cAMP-responsive element, AP-1, AP-2 and SP-1 sites in the 5' flanking region. This region also contains three direct-repeat-like sequences that resemble the consensus binding sequence for members of the steroid-thyroid hormone receptor superfamily, including vitamin D. Deletion analysis suggests that the region from nt-2409 to nt-1259 in the 5' flanking region may be involved in kidney-specific gene expression. Vitamin D responsiveness of the NaPi-3 promoter was also detected in COS-7 cells co-transfected with a human vitamin D receptor expression vector. The presence of the three vitamin D receptor- responsive elements in the NaPi-3 promoter may be important in mediating the enhanced expression of the gene by 1,25-dihydroxyvitamin D3.

Bimolecular exon ligation by the human spliceosome

Intron excision is an essential step in eukaryotic gene expression, but the molecular mechanisms by which the spliceosome accurately identifies splice sites in nuclear precursors to messenger RNAs (pre-mRNAs) are not well understood. A bimolecular assay for the second step of splicing has now revealed that exon ligation by the human spliceosome does not require covalent attachment of a 3' splice site to the branch site. Furthermore, accurate definition of the 3' splice site in this system is independent of either a covalently attached polypyrimidine tract or specific 3' exon sequences. Rather, in this system 3' splice site selection apparently occurs with a 5' --> 3' directionality.

Functional analysis of the polypyrimidine tract in pre-mRNA splicing

The polypyrimidine tract is one of the important cis-acting sequence elements directing intron removal in pre-mRNA splicing. Progressive deletions of the polypyrimidine tract have been found to abolish correct lariat formation, spliceosome assembly and splicing. In addition, the polypyrimidine tract can alter 3'-splice site selection by promoting alternative branch site selection. However, there appears to be great flexibility in the specific sequence of a given tract. Not only the optimal composition of the polypyrimidine tract, but also the role of the tract in introns with no apparent polypyrimidine tracts or where changes in the tract are apparently harmless are uncertain. Accordingly, we have designed a series of cis-competition splicing constructs to test the functional competitive efficiency of a variety of systematically mutated polypyrimidine tracts. An RT/PCR assay was used to detect spliced product formation as a result of differential branch point selection dependent on direct competition between two opposing polypyrimidine tracts. We found that pyrimidine tracts containing 11 continuous uridines are the strongest pyrimidine tracts. In such cases, the position of the uridine stretch between the branch point and 3'-splice site AG is unimportant. In contrast, decreasing the continuous uridine stretch to five or six residues requires that the tract be located immediately adjacent to the AG for optimal competitive efficiency. The block to splicing with decreasing polypyrimidine tract strength is primarily prior to the first step of splicing. While lengthy continuous uridine tracts are the most competitive, tracts with decreased numbers of consecutive uridines and even tracts with alternating purine/pyrimidine residues can still function to promote branch point selection, but are far less effective competitors in 3'-splice site selection assays.

Species-specific alternative splicing of the epidermal growth factor-like domain 1 of cartilage aggrecan

Aggrecan transcripts contain two alternatively spliced exons that code for two epidermal growth factor-like domains (EGF1 and EGF2). Whereas the EGF2 sequence is expressed at a uniform level among different species, the EGF1 sequence has been detected only in human aggrecan transcripts. In this study we have used the nested primer reverse transcription-PCR (RT-PCR) method to compare the expression of the EGF1 exon in human, bovine and dog aggrecan transcripts. Our results indicate that this exon is expressed in a species-specific manner. In addition to its significant expression level in human transcripts, the EGF1 sequence can be detected in a small portion of bovine aggrecan transcripts as shown with nested primer RT-PCR. In contrast, the same module is not detectable in dog aggrecan transcripts, although an EGF1 sequence is present in the dog aggrecan gene. The expression level of the EGF1 exon in the aggrecan transcripts correlates with the strength of the polypyrimidine tract upstream of the exon. The EGF1 sequence also shows much less conservation between the species than the EGF2 sequence. The species-specific expression and high sequence variation of the EGF1 exon imply that this sequence is likely to code for an aggrecan domain having no cartilage-specific function.

Sequence of the polypyrimidine tract of the 3'-terminal 3' splicing signal can affect intron-dependent pre-mRNA processing in vivo

Most pre-mRNAs require an intron for efficient processing in higher eukaryotes. However, not all introns can provide this function. For example, transcripts synthesized from a variant of the human beta-globin gene lacking its second intervening sequence (IVS2), yet retaining its first intervening sequence (IVS1), exhibit multiple defects in mRNA biogenesis. To investigate why, we transfected into monkey cells plasmids containing the human beta-globin gene and variants of it altered in (i) IVS1, (ii) the 3'-terminal exon, and (iii) the polyadenylation signal. The beta-globin RNAs accumulated in these cells were analyzed by quantitative S1 nuclease mapping for nuclear accumulation, intron excision, polyadenylation and cytoplasmic accumulation. We found that the 3' splicing signal of IVS1, with multiple purines interrupting its polypyrimidine tract, could efficiently function as an internal 3' splicing signal; however, it could not efficiently function as the 3'-terminal 3' splicing signal for any of these steps in intron-dependent mRNA biogenesis unless (i) its polypyrimidine tract was made uninterrupted in pyrimidines, or (ii) specific sequences were deleted from the 3'-terminal exon. We conclude that whether an intron can provide the function necessary for efficient processing of intron-dependent pre-mRNA is dependent upon the ability of its 3' splicing signal to define the 3'-terminal exon. On the practical side, this finding means one needs to consider both the sequence and location of the intron to be included in an intron-dependent gene to obtain efficient expression in vivo.

Polymorphism of the mouse E2A gene due to an intronic deletion of 536 bp

Examination of genetic polymorphism of the transcription factor-encoding gene E2A in laboratory and wild mice by Southern blotting has revealed the presence of two alleles. The most frequent allele is found in Mus musculus (Mm) musculus, as well as Mm domesticus. The less common allele is restricted to the Mm domesticus subspecies. Characterisation of these alleles has shown that the less common allele contains a deletion of approx. 500 bp located within a 1.8-kb intron immediately upstream from the E12 basic helix-loop-helix exon. DNA sequencing determined the deletion to span 536 bp including nucleotides 1045-1580 of the intron within the common allele. The deleted region includes several sequences with similarity to gene regulatory motifs; however, expression of E12 and intron splicing appeat unaltered. The occurrence of an identical deletion in mice from different geographical regions suggests that the uncommon allele may have a long evolutionary history.

Human adenovirus type 5 variants with sequence alterations flanking the E2A gene: effects on E2 expression and DNA replication

The human adenovirus type 5 (Ad5) E2 transcription unit is divided into a promoter-proximal region, E2A, and a distal region, E2B, each with its own polyadenylation site. Together these regions encode the three virus-derived proteins necessary for genome replication. Ad5 variants were produced that carried linker insertion mutations immediately 5' and/or 3' to the coding sequence for the E2A gene DNA binding protein (DBP). Two variants carrying solely a 5' lesion showed decreased usage of the adjacent 3' splice site, via which the DBP mRNA is produced, and an increased usage of the alternative downstream splice sites in the E2B region, wherein viral DNA polymerase and terminal protein precursor are encoded; these viruses showed somewhat reduced growth. A variant carrying a 3' lesion showed a marginal increase in DBP expression and slightly accelerated growth. When lesions 5' and 3' to the DBP coding sequence were combined in cis, the resulting virus was severely defective for growth and expressed E2B products to the virtual exclusion of E2A DBP. These data indicate that interactions must occur between the E2A 3' splice site and polyadenylation site before this region can be treated as an exon by the RNA processing machinery, and that a sequence alteration at the polyadenylation site that alone has only minor effects on the pattern of RNA processing can drastically affect terminal exon usage when placed in cis with a mutation that reduces splicing efficiency at the upstream 3' splice site. The data further indicate that, in vivo, Ad5 DNA replication is limited by prevailing DBP levels rather than by levels of polymerase or terminal protein precursor.

A T to G mutation in the polypyrimidine tract of the second intron of the human beta-globin gene reduces in vitro splicing efficiency: evidence for an increased hnRNP C interaction

In a patient with a beta-thalassemia intermedia, a mutation was identified in the second intron of the human beta-globin gene. The U-->G mutation is located within the polypyrimidine tract at position -8 upstream of the 3' splice site. In vivo, this mutation leads to decreased levels of the hemoglobin protein. Because of the location of the mutation and the role of the polypyrimidine tract in the splicing process, we performed in vitro splicing assays on the pre-messenger RNA (pre-mRNA). We found that the splicing efficiency of the mutant pre-mRNA is reduced compared to the wild type and that no cryptic splice sites are activated. Analysis of splicing complex formation shows that the U-->G mutation affects predominantly the progression of the H complex towards the pre-spliceosome complex. By cross-linking and immunoprecipitation assays, we show that the hnRNP C protein interacts more efficiently with the mutant precursor than with the wild-type. This stronger interaction could play a role, directly or indirectly, in the decreased splicing efficiency.

A novel role for a U5 snRNP protein in 3' splice site selection

The choice of a 3' splice site in Saccharomyces cerevisiae introns involves recognition of a uridine-rich tract upstream of the AG dinucleotide splice junction. By isolating mutants that eliminate the normal preference for uridine-containing 3' splice sites in a cis-competition, we identified a mutation that is an allele of PRP8, prp8-101. This was unexpected because previous analysis has demonstrated that the U5 snRNP protein encoded by PRP8 is required for spliceosome assembly prior to the first catalytic step of splicing. In contrast, the uridine recognition defect caused by the prp8-101 mutation selectively inhibits the second catalytic step of splicing. This defect is seen not only in 3' splice site cis-competitions but also in the splicing of an unusual intron in the TUB3 gene and in the ACT1 intron when utilization of its 3' splice site is rate limiting for splicing. Consistent with a direct role in 3' splice site selection, Prp8 can be cross-linked to the 3' splice site during the splicing reaction. These data demonstrate a novel function for Prp8 in 3' splice site recognition and utilization.

Molecular characterization of the mouse ribosomal protein S24 multigene family: a uniquely expressed intron-containing gene with cell-specific expression of three alternatively spliced mRNAs

A family of 16 genes encoding the mouse ribosomal protein S24 was identified, and four members from this family were cloned. A single expressed intron-containing S24 gene (termed mrpS24) and one pseudogene (mrpS24p) were completely sequenced and characterized. The mrpS24 gene has seven exons and six introns spanning over 5.1 x 10(3) nucleotides (nt). The cap site of S24 was mapped to a G residue four nt upstream of a polypyrimidine tract and 15 nt downstream of a TATA-like (TATGA) element. The 5' region (-325 to +33) of the mrpS24 gene has a functional promoter that was able to express the fused chloramphenicol acetyltransferase (CAT) reporter gene. Two different forms of mouse S24 cDNA clones were previously isolated. Sequence analysis showed that one of these cDNA clones (termed S24a) lacks the entire exon V sequence (18 nt), and the deduced amino acid sequence is missing a C-terminal lysine residue encoded by the other cDNA (S24b). The pseudogene mrpS24p is flanked by an 11-bp direct repeat, and its sequence is almost identical to the S24 cDNA sequence, but it lacks two mini-exons, V and VI (20 nt), as in the cases of the human and rat S24 cDNAs. RT-PCR experiments demonstrated the existence of a third form (S24c) that similarly lacks both of the mini-exons, and suggested that different species of S24 mRNA might arise from alternative splicing of the mini-exons V and VI. Northern blot analysis showed that S24 expression is down- and up-regulated during adipocyte differentiation and in cellular transformation, respectively. RNase protection assays and RT-PCR experiments suggested that these cell-specific changes of S24 mRNA levels are mainly due to fluctuations in S24c mRNA level. Our results provide the first indication that a ribosomal protein gene is regulated by alternative usage of two mini-exons in a cell-specific manner.

Acceptor splice site mutation in the invariant AG of intron 5 of the protein C gene, causing type I protein C deficiency

An acceptor splice-site mutation (3318, A-->G) in the invariant AG of intron 5 of the human protein C gene has been identified in a Spanish family with heterozygous type I protein C (PC) deficiency and thromboembolic disease. Family studies confirmed cosegregation of the mutation with type I PC deficiency. Computer analysis of the mutated sequence predicted the normal splicing site to be abolished by this mutation, whereas a cryptic splice site located two nucleotides downstream, in exon 6, is probably activated. According to this, 3318, A-->G should result in a frameshift with a stop at codon 119, in agreement with the presence of a type I or quantitative PC deficient phenotype in the affected members of the family.

Analytical expression of the purine/pyrimidine codon probability after and before random mutations

Recently, we proposed a new model of DNA sequence evolution (Arquès and Michel. 1990b. Bull. math. Biol. 52, 741-772) according to which actual genes on the purine/pyrimidine (R/Y) alphabet (R = purine = adenine or guanine, Y = pyrimidine = cytosine or thymine) are the result of two successive evolutionary genetic processes: (i) a mixing (independent) process of non-random oligonucleotides (words of base length less than 10: YRY(N)6, YRYRYR and YRYYRY are so far identified; N = R or Y) leading to primitive genes (words of several hundreds of base length) and followed by (ii) a random mutation process, i.e., transformations of a base R (respectively Y) into the base Y (respectively R) at random sites in these primitive genes. Following this model the problem investigated here is the study of the variation of the 8 R/Y codon probabilities RRR, ..., YYY under random mutations. Two analytical expressions solved here allow analysis of this variation in the classical evolutionary sense (from the past to the present, i.e., after random mutations), but also in the inverted evolutionary sense (from the present to the past, i.e., before random mutations). Different properties are also derived from these formulae. Finally, a few applications of these formulae are presented. They prove the proposition in Arquès and Michel (1990b. Bull. math. Biol. 52, 741-772), Section 3.3.2, with the existence of a maximal mean number of random mutations per base of the order 0.3 in the protein coding genes. They also confirm the mixing process of oligonucleotides by excluding the purine/pyrimidine contiguous and alternating tracts from the formation process of primitive genes.

Solution structure of lariat RNA by 500 MHz NMR spectroscopy and molecular dynamics studies in water

A 500 MHz NMR study of the lariat RNA tetramer 1 and pentamer 2 mimicking the naturally occurring lariat RNA is reported. The conformational properties of 1 and 2 were compared with those of a linear branched RNA tetramer 3, which show that the conformational features of the two lariat RNAs, 1 and 2, are quite constrained and significantly different from those observed for the linear branched RNA tetramer 3. The conformation of all sugar residues forming the lariat ring in 1 and 2 are locked in a rigid South-type conformation. All residues in both lariat RNAs have a high population of gamma+ (67-85%) and beta t (95-100%) rotamers except guanosine where the gamma+ population is low. The conformation around the glycosidic bond is anti for all residues except for guanosine where NOE data indicates an equilibrium of syn<-->anti. In both lariat RNAs, 1 and 2, the temperature dependent 1H and 31P chemical shifts as well as the oligomerization shifts, with respect to adenosine 2',3',5'-triethyl-phosphate (Sund et al., 1992, Tetrahedron 48, 695) suggests that the 3'-->5' linked U4 or C4 residue is stacked on guanosine. Subsequently, 1H-1H, 1H-31P and 13C-31P coupling constants derived torsional constraints were used for molecular dynamics study in water with counter sodium ions for a total of 226 ps. The MD simulations were first carried out with harmonic torsional constraints which were derived from J couplings (0-86 ps) and then completely without constraints (96-226 ps). The lack of any major changes in the conformation of the two lariat-RNA structures upon releasing the NMR constraints indicate that the conformers generated in the MD simulation in water agree well with the structural features suggested by experimental observables. This means that the ensemble of conformers generated during the MD trajectory of 226 ps are not artificially held in these conformations due to the NMR constraints, suggesting that these conformers can be considered to be good representatives of the actual NMR observed solution structures.

Scanning and competition between AGs are involved in 3' splice site selection in mammalian introns

In mammalian intron splicing, the mechanism by which the 3' splice site AG is accurately and efficiently identified has remained unresolved. We have previously proposed that the 3' splice site in mammalian introns is located by a scanning mechanism for the first AG downstream of the branch point-polypyrimidine tract. We now present experiments that lend further support to this model while identifying conditions under which competition can occur between adjacent AGs. The data show that the 3' splice site is identified as the first AG downstream from the branch point by a mechanism that has all the characteristics expected of a 5'-to-3' scanning process that starts from the branch point rather than the pyrimidine tract. Failure to recognize the proximal AG may arise, however, from extreme proximity to the branch point or sequestration within a hairpin. Once an AG has been encountered, the spliceosome can still see a limited stretch of downstream RNA within which an AG more competitive than the proximal one may be selected. Proximity to the branch point is a major determinant of competition, although steric effects render an AG less competitive in close proximity (approximately 12 nucleotides). In addition, the nucleotide preceding the AG has a striking influence upon competition between closely spaced AGs. The order of competitiveness, CAG congruent to UAG > AAG > GAG, is similar to the nucleotide preference at this position in wild-type 3' splice sites. Thus, 3' splice site selection displays properties of both a scanning process and competition between AGs based on immediate sequence context. This refined scanning model, incorporating elements of competition, is the simplest interpretation that is consistent with all of the available data.

3' splice site selection in dicot plant nuclei is position dependent

In contrast to mammalian and yeast systems, the mechanism for intron recognition and splice site selection in plant pre-mRNAs is poorly understood. Splice site sequences and putative branchpoint sequences are loosely conserved in plant introns compared with other eukaryotes. Perhaps to compensate for these variations, plant introns are significantly richer in adenosine and uridine residues than are their adjacent exons. To define elements critical for 3' splice site selection in dicotyledonous plant nuclei, pre-mRNA transcripts containing intron 3 of the maize Adh1 gene were expressed in Nicotiana benthamiana nuclei by using an autonomously replicating plant expression vector. Using a series of intron rearrangements which reposition the 3' intron-exon border, we demonstrate that the normal 3' splice site is defined in a position-dependent manner and that cryptic 3' splice sites within the intron are masked by the presence of a functional downstream 3' splice site. Disruption of the AU-rich elements upstream from the normal 3' splice site indicates that multiple AU elements between -66 and -6 cooperatively define the 3' boundary of the intron. These results are consistent with a model for plant intron recognition in which AU-rich elements spread throughout the length of the intron roughly define the intron boundaries by generating strong AU transition points. Functional 3' splice sites located downstream from these AU-rich sequences are preferentially selected over sites embedded within them.

Scanning and competition between AGs are involved in 3' splice site selection in mammalian introns

In mammalian intron splicing, the mechanism by which the 3' splice site AG is accurately and efficiently identified has remained unresolved. We have previously proposed that the 3' splice site in mammalian introns is located by a scanning mechanism for the first AG downstream of the branch point-polypyrimidine tract. We now present experiments that lend further support to this model while identifying conditions under which competition can occur between adjacent AGs. The data show that the 3' splice site is identified as the first AG downstream from the branch point by a mechanism that has all the characteristics expected of a 5'-to-3' scanning process that starts from the branch point rather than the pyrimidine tract. Failure to recognize the proximal AG may arise, however, from extreme proximity to the branch point or sequestration within a hairpin. Once an AG has been encountered, the spliceosome can still see a limited stretch of downstream RNA within which an AG more competitive than the proximal one may be selected. Proximity to the branch point is a major determinant of competition, although steric effects render an AG less competitive in close proximity (approximately 12 nucleotides). In addition, the nucleotide preceding the AG has a striking influence upon competition between closely spaced AGs. The order of competitiveness, CAG congruent to UAG > AAG > GAG, is similar to the nucleotide preference at this position in wild-type 3' splice sites. Thus, 3' splice site selection displays properties of both a scanning process and competition between AGs based on immediate sequence context. This refined scanning model, incorporating elements of competition, is the simplest interpretation that is consistent with all of the available data.

3' splice site selection in dicot plant nuclei is position dependent

In contrast to mammalian and yeast systems, the mechanism for intron recognition and splice site selection in plant pre-mRNAs is poorly understood. Splice site sequences and putative branchpoint sequences are loosely conserved in plant introns compared with other eukaryotes. Perhaps to compensate for these variations, plant introns are significantly richer in adenosine and uridine residues than are their adjacent exons. To define elements critical for 3' splice site selection in dicotyledonous plant nuclei, pre-mRNA transcripts containing intron 3 of the maize Adh1 gene were expressed in Nicotiana benthamiana nuclei by using an autonomously replicating plant expression vector. Using a series of intron rearrangements which reposition the 3' intron-exon border, we demonstrate that the normal 3' splice site is defined in a position-dependent manner and that cryptic 3' splice sites within the intron are masked by the presence of a functional downstream 3' splice site. Disruption of the AU-rich elements upstream from the normal 3' splice site indicates that multiple AU elements between -66 and -6 cooperatively define the 3' boundary of the intron. These results are consistent with a model for plant intron recognition in which AU-rich elements spread throughout the length of the intron roughly define the intron boundaries by generating strong AU transition points. Functional 3' splice sites located downstream from these AU-rich sequences are preferentially selected over sites embedded within them.

A functional association between the 5' and 3' splice site is established in the earliest prespliceosome complex (E) in mammals

The earliest detectable mammalian prespliceosome complex (E) contains the non-snRNP splicing factor U2AF, U1 snRNP, and several spliceosome-associated proteins (SAPs). We show that specific complexes, designated E3' and E5', assemble independently on RNAs containing only a 3' or 5' splice site, respectively. U2AF is enriched in E3', whereas U1 snRNP is enriched in E5'. Using a highly sensitive substrate-competition assay, we show that both the 5' splice site and the pyrimidine tract at the 3' splice site are required for efficient E complex assembly on intact pre-mRNA. We conclude that the 5' and 3' splice sites are associated functionally as early as E complex by either direct or indirect interactions between U1 snRNP and U2AF. Our observations predict that E complex assembly is a major control point for establishing splice site selection in both constitutively and alternatively spliced pre-mRNAs.

Expression of maize Adh1 intron mutants in tobacco nuclei

In vivo and in vitro gene transfer experiments have suggested that the elements mediating intron recognition differ in mammalian, yeast and plant nuclei. Differences in the sequence dependencies, which also exist between dicotyledonous and monocotyledonous nuclei, have prevented some monocot introns from being spliced in dicot nuclei. To locate elements which modulate efficient recognition of introns in dicot nuclei, the maize Adh1 gene has been expressed in full-length and single intron constructs in Nicotiana benthamiana nuclei using an autonomously replicating plant expression vector. Quantitative PCR-Southern analyses indicate that the inefficient splicing of the maize Adh1 intron 1 (57% AU) in these dicot nuclei can be dramatically enhanced by increasing the degree of U1 snRNA complementarity at the 5' splice site. This indicates that the 5' splice site plays a significant role in defining the splicing efficiency of an intron in dicot nuclei and that, most importantly, the remainder of this monocot intron contains no elements which inhibit its accurate recognition in dicot nuclei. Deletions in intron 3 (66% AU) which effectively move the 3' boundary between AU-rich intron and GC-rich exon sequences strongly activate a cryptic upstream splice site; those which do not reposition this boundary activate a downstream cryptic splice site. This suggests that 3' splice site selection in dicot nuclei is extremely flexible and not dependent on strict sequence requirements but rather on the transition points between introns and exons. Our results are consistent with a model in which potential splice sites are selected if they are located upstream (5' splice site) or downstream (3' splice site) of AU transition points and not if they are embedded within AU-rich sequences.