The prediction of exons through an analysis of spliceable open reading frames.

We have developed a computer program which predicts internal exons from naive genomic sequence data and which will run on any IBM-compatible 80286 (or higher) computer. The algorithm searches a sequence for 'spliceable open reading frames' (SORFs), which are open reading frames bracketed by suitable splice-recognition sequences, and then analyzes the region for codon usage. Potential exons are stratified according to the reliability of their prediction, from confidence levels 1 to 5. The program is designed to predict internal exons of length greater than 60 nucleotides. In an analysis of 116 genes of a training set, 384 out of 441 such exons (87.1%) are identified, with 280 (63.5%) of predictions matching the true exon exactly (at both 5' and 3' splice junctions and in the correct reading frame), and with 104 (23.6%) exons matching partially. In a similar analysis of 14 genes in a test set unrelated to the genes used to generate the parameters of the program, 70 out of 80 internal exons greater than 60 bp in length are identified (87.5%), with 47 completely and 23 partially matched. SORFs that partially match true internal exons share at least one splice junction with the exon, or share both splice junctions but are interpreted in an incorrect reading frame. Specificity (the percentage of SORFs that correspond to true exons) varies from 91% at confidence level 1 to 16% at confidence level 5, with an overall specificity of 35-40%. The output displays nucleotide position, confidence level, reading frame phase at the 5' and 3' ends, acceptor and donor sequences and scoring statistics and also gives an amino acid translation of the potential exon. SORFIND compares favourably with other programs currently used to predict protein-coding regions.     

Epstein-Barr virus mRNAs produced by alternative splicing.

The structure of Epstein-Barr virus mRNAs transcribed in B95-8 cells has been studied by cDNA cloning and sequencing. We present here the analysis of four cDNAs. The corresponding mRNAs are probably transcribed from a single promoter located in the US region. They are produced by alternative splicing of exons transcribed from the US, IR and UL regions. The exons are spread over 100 kbp. The exons from the IR region constitute a unit which is repeated several times. The cDNAs share the exons from the US and IR regions. Some of the cDNAs also share some of the exons from the UL region. Each cDNA contains a long open reading frame or the 5' end of a long open reading frame which ends several hundred nucleotides downstream on the viral genome. The 5' untranslated regions are unusually long. Three mRNA species differing in their 5' untranslated regions may encode for the nuclear antigen EBNA-1. The other mRNAs encode for polypeptides which may not have any common region.     

Analysis and recognition of 5' UTR intron splice sites in human pre-mRNA

Prediction of splice sites in non-coding regions of genes is one of the most challenging aspects of gene structure recognition. We perform a rigorous analysis of such splice sites embedded in human 5' untranslated regions (UTRs), and investigate correlations between this class of splice sites and other features found in the adjacent exons and introns. By restricting the training of neural network algorithms to 'pure' UTRs (not extending partially into protein coding regions), we for the first time investigate the predictive power of the splicing signal proper, in contrast to conventional splice site prediction, which typically relies on the change in sequence at the transition from protein coding to non-coding. By doing so, the algorithms were able to pick up subtler splicing signals that were otherwise masked by 'coding' noise, thus enhancing significantly the prediction of 5' UTR splice sites. For example, the non-coding splice site predicting networks pick up compositional and positional bias in the 3' ends of non-coding exons and 5' non-coding intron ends, where cytosine and guanine are over-represented. This compositional bias at the true UTR donor sites is also visible in the synaptic weights of the neural networks trained to identify UTR donor sites. Conventional splice site prediction methods perform poorly in UTRs because the reading frame pattern is absent. The NetUTR method presented here performs 2-3-fold better compared with NetGene2 and GenScan in 5' UTRs. We also tested the 5' UTR trained method on protein coding regions, and discovered, surprisingly, that it works quite well (although it cannot compete with NetGene2). This indicates that the local splicing pattern in UTRs and coding regions is largely the same. The NetUTR method is made publicly available at www.cbs.dtu.dk/services/NetUTR.     

Organization and chromosomal localization of the human interleukin 5 receptor alpha-chain gene.

The gene for the hIL-5R alpha subunit, which is present in a single copy in the human genome, has been analysed in detail. It is located on chromosome 3 in the region 3p26. The gene organization reflects its relationship to the cytokine/haematopoietin receptor superfamily. Three introns are located in the 5' untranslated region. The subsequent exons determine the functional domains of the hIL-5R alpha protein: the signal peptide, three fibronectin type III-like (FN-like) modules, each built up by two exons, the membrane anchor and two exons forming the cytoplasmic tail, the first of which contains the proline cluster region. In addition, a specific exon generating a soluble isoform is located before the membrane anchor exon. This specific exon contains an in frame TAA stop codon, followed by a polyadenylation signal. Hence, a normal splicing event leads to a soluble IL-5R alpha variant, whereas alternative splicing is required for cell membrane anchoring. A second area of alternative splicing is found in the 5' leader sequence, and possibly relates to the presence of short open reading frames preceding the main ATG. All intron-exon junctions meet the GT-AG rule. The gene structures of all cytokine/haematopoietin receptors documented so far have also been compared with respect to intron phasing. This shows that all introns between the FN-III-like modules are of the +1 type, but in addition, splice sites within the Cys-module and WS-WS-module are invariably of the +2 and 0 type, respectively.     

Characterization of multiple mRNA species of simian immunodeficiency virus from macaques in a CD4+ lymphoid cell line.

Cytoplasmic poly(A)+ RNA was isolated from CEMX721.174 cells 5 to 10 days after infection with molecularly cloned simian immunodeficiency virus SIVmac239. Expression of SIV RNA was analyzed by Northern (RNA) blot hybridization and by sequencing of cDNA clones. As expected, a splice donor site was demonstrated in the untranslated leader sequence outside the left long terminal repeat. The region between pol and env was found to contain at least two splice donor and six splice acceptor sites. Splice acceptor and donor sites in the intergenic region were suitably positioned for expression of vpx, vpr, tat, and rev. Splice acceptor sites at nucleotides 8802 and 8805 were demonstrated in singly and doubly spliced RNAs with the potential of expressing nef and the second exons of tat and rev. Our results demonstrate a complex pattern of alternative splicing of SIV mRNAs. The results are very similar to what has been observed in human immunodeficiency virus type 1-infected cells, suggesting that both human and simian immunodeficiency viruses are subject to multiple levels of regulation.