The mitochondrial genome of the fission yeast Schizosaccharomyces pombe. The cytochrome b gene has an intron closely related to the first two introns in the Saccharomyces cerevisiae cox1 gene

The DNA sequence of the cob region of the Schizosaccharomyces pombe mitochondrial DNA has been determined. The cytochrome b structural gene is interrupted by an intron of 2526 base-pairs, which has an open reading frame of 2421 base-pairs in phase with the upstream exon. The position of the intron differs from those found in the cob genes of Saccharomyces cerevisiae, Aspergillus nidulans or Neurospora crassa. The Sch. pombe cob intron has the potential of assuming an RNA secondary structure almost identical to that proposed for the first two cox1 introns (group II) in S. cerevisiae and the p1-cox1 intron in Podospora anserina. It has most of the consensus nucleotides in the central core structure described for this group of introns and its comparison with other group II introns allows the identification of an additional conserved nucleotide stretch. A comparison of the predicted protein sequences of group II intronic coding regions reveals three highly conserved blocks showing pairwise amino acid identities of 34 to 53%. These regions comprise over 50% of the coding length of the intron but do not include the 5' region, which has strong secondary structural features. In addition to the potential intron folding, long helical structures involving repetitive sequences can be formed in the flanking cob exon regions. A comparison of the Sch. pombe cytochrome b sequence with those available from other organisms indicates that Sch. pombe is evolutionarily distant from both budding yeasts and filamentous fungi. As was seen for the Sch. pombe cox1 gene (Lang, 1984), the cob exons are translated using the universal genetic code and this distinguishes Sch. pombe mitochondria from all other fungal and animal mitochondrial systems.     

Gene expression profile of human bone marrow stromal cells: high-throughput expressed sequence tag sequencing analysis

Human bone marrow stromal cells (HBMSC) are pluripotent cells with the potential to differentiate into osteoblasts, chondrocytes, myelosupportive stroma, and marrow adipocytes. We used high-throughput DNA sequencing analysis to generate 4258 single-pass sequencing reactions (known as expressed sequence tags, or ESTs) obtained from the 5' (97) and 3' (4161) ends of human cDNA clones from a HBMSC cDNA library. Our goal was to obtain tag sequences from the maximum number of possible genes and to deposit them in the publicly accessible database for ESTs (dbEST of the National Center for Biotechnology Information). Comparisons of our EST sequencing data with nonredundant human mRNA and protein databases showed that the ESTs represent 1860 gene clusters. The EST sequencing data analysis showed 60 novel genes found only in this cDNA library after BLAST analysis against 3.0 million ESTs in NCBI's dbEST database. The BLAST search also showed the identified ESTs that have close homology to known genes, which suggests that these may be newly recognized members of known gene families. The gene expression profile of this cell type is revealed by analyzing both the frequency with which a message is encountered and the functional categorization of expressed sequences. Comparing an EST sequence with the human genomic sequence database enables assignment of an EST to a specific chromosomal region (a process called digital gene localization) and often enables immediate partial determination of intron/exon boundaries within the genomic structure. It is expected that high-throughput EST sequencing and data mining analysis will greatly promote our understanding of gene expression in these cells and of growth and development of the skeleton.     

Evidence for splice site pairing via intron definition in Schizosaccharomyces pombe

Schizosaccharomyces pombe pre-mRNAs are generally multi-intronic and share certain features with pre-mRNAs from Drosophila melanogaster, in which initial splice site pairing can occur via either exon or intron definition. Here, we present three lines of evidence suggesting that, despite these similarities, fission yeast splicing is most likely restricted to intron definition. First, mutating either or both splice sites flanking an internal exon in the S. pombe cdc2 gene produced almost exclusively intron retention, in contrast to the exon skipping observed in vertebrates. Second, we were unable to induce skipping of the internal microexon in fission yeast cgs2, whereas the default splicing pathway excludes extremely small exons in mammals. Because nearly quantitative removal of the downstream intron in cgs2 could be achieved by expanding the microexon, we propose that its retention is due to steric occlusion. Third, several cryptic 5' junctions in the second intron of fission yeast cdc2 are located within the intron, in contrast to their generally exonic locations in metazoa. The effects of expanding and contracting this intron are as predicted by intron definition; in fact, even highly deviant 5' junctions can compete effectively with the standard 5' splice site if they are closer to the 3' splicing signals. Taken together, our data suggest that pairing of splice sites in S. pombe most likely occurs exclusively across introns in a manner that favors excision of the smallest segment possible.     

Conserved modularity and potential for alternate splicing in mouse and human Slit genes

The vertebrate Slit gene family currently consists of three members; Slit1, Slit2 and Slit3. Each gene encodes a protein containing multiple epidermal growth factor and leucine rich repeat motifs, which are likely to have importance in cell-cell interactions. In this study, we sought to fully define and characterise the vertebrate Slit gene family. Using long distance PCR coupled with in silico mapping, we determined the genomic structure of all three Slit genes in mouse and man. Analysis of EST and genomic databases revealed no evidence of further Slit family members in either organism. All three Slit genes were encoded by 36 (Slit3) or 37 (Slit1 and Slit2) exons covering at least 143 kb or 183 kb of mouse or human genomic DNA respectively. Two additional potential leucine-rich repeat encoding exons were identified within intron 12 of Slit2. These could be inserted in frame, suggesting that alternate splicing may occur in Slit2. A search for STS sequences within human Slit3 anchored this gene to D5S2075 at the 5' end (exon 4) and SGC32449 within the 3' UTR, suggesting that Slit3 may cover greater than 693 kb. The genomic structure of all Slit genes demonstrated considerable modularity in the placement of exon-intron boundaries such that individual leucine-rich repeat motifs were encoded by individual 72 bp exons. This further implies the potential generation of multiple Slit protein isoforms varying in their number of repeat units. cDNA library screening and EST database searching verified that such alternate splicing does occur.     

Introns and splicing elements of five diverse fungi

Genomic sequences and expressed sequence tag data for a diverse group of fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, Neurospora crassa, and Cryptococcus neoformans) provided the opportunity to accurately characterize conserved intronic elements. An examination of large intron data sets revealed that fungal introns in general are short, that 98% or more of them belong to the canonical splice site (ss) class (5'GU...AG3'), and that they have polypyrimidine tracts predominantly in the region between the 5' ss and the branch point. Information content is high in the 5' ss, branch site, and 3' ss regions of the introns but low in the exon regions adjacent to the introns in the fungi examined. The two yeasts have broader intron length ranges and correspondingly higher intron information content than the other fungi. Generally, as intron length increases in the fungi, so does intron information content. Homologs of U2AF spliceosomal proteins were found in all species except for S. cerevisiae, suggesting a nonconventional role for U2AF in the absence of canonical polypyrimidine tracts in the majority of introns. Our observations imply that splicing in fungi may be different from that in vertebrates and may require additional proteins that interact with polypyrimidine tracts upstream of the branch point. Theoretical protein homologs for Nam8p and TIA-1, two proteins that require U-rich regions upstream of the branch point to function, were found. There appear to be sufficient differences between S. cerevisiae and S. pombe introns and the introns of two filamentous members of the Ascomycota and one member of the Basidiomycota to warrant the development of new model organisms for studying the splicing mechanisms of fungi.     

The mitochondrial genome of the fission yeast Schizosaccharomyces pombe: highly homologous introns are inserted at the same position of the otherwise less conserved cox1 genes in Schizosaccharomyces pombe and Aspergillus nidulans.

The DNA sequence of the second intron in the mitochondrial gene for subunit 1 of cytochrome oxidase (cox1), and the 3'' part of the structural gene have been determined in Schizosaccharomyces pombe. Comparing the presumptive amino acid sequence of the 3'' regions of the cox1 genes in fungi reveals similarly large evolutionary distances between Aspergillus nidulans, Saccharomyces cerevisiae and S. pombe. The comparison of exon sequences also reveals a stretch of only low homology and of general size variation among the fungal and mammalian genes, close to the 3'' ends of the cox1 genes. The second intron in the cox1 gene of S. pombe contains an open reading frame, which is contiguous with the upstream exon and displays all characteristics common to class I introns. Three findings suggest a recent horizontal gene transfer of this intron from an Aspergillus type fungus to S. pombe. (i) The intron is inserted at exactly the same position of the cox1 gene, where an intron is also found in A. nidulans. (ii) Both introns contain the highest amino acid homology between the intronic unassigned reading frames of all fungi identified so far (70% identity over a stretch of 253 amino acids). However, in the most homologous region, a GC-rich sequence is inserted in the A. nidulans intron, flanked by two direct repeats of 5 bp. The 37-bp insert plus 5 bp of direct repeat amounts to an extra 42 bp in the A. nidulans intron. (iii) TGA codons are the preferred tryptophan codons compared with TGG in all mitochondrial protein coding sequences of fungi and mammalia.(ABSTRACT TRUNCATED AT 250 WORDS)     

IsoSVM – Distinguishing isoforms and paralogs on the protein level

Background  

Recent progress in cDNA and EST sequencing is yielding a deluge of sequence data. Like database search results and proteome databases, this data gives rise to inferred protein sequences without ready access to the underlying genomic data. Analysis of this information (e.g. for EST clustering or phylogenetic reconstruction from proteome data) is hampered because it is not known if two protein sequences are isoforms (splice variants) or not (i.e. paralogs/orthologs). However, even without knowing the intron/exon structure, visual analysis of the pattern of similarity across the alignment of the two protein sequences is usually helpful since paralogs and orthologs feature substitutions with respect to each other, as opposed to isoforms, which do not.     

Structure of the murine complement factor H gene

Factor H is a regulatory protein of the alternative pathway of complement activation comprised of 20 tandem repeating units of 60 amino acids each. A factor H cDNA clone was used to identify 17 genomic clones from a cosmid library. Four clones were selected for analysis of intron/exon junctions and 5' and 3' regions of the gene and for mapping of the exons. The factor H gene was found to be comprised of 22 exons. Each repeating unit is encoded by one exon, except the second repeat, which is coded by two exons; the leader sequence is encoded by a separate exon. The exons range in size from 77 to 210 base pairs (bp) and average 178 bp. They span a region of approximately 100 kilobases (kb) on chromosome 1. The leader sequence exon is 26 kb upstream of the first repeat exon, representing the largest intron. The other introns range in size from 86 bp to 12.9 kb, and the average intron size is 4.7 kb. Analysis of the genomic organization of the factor H gene has provided insight into the protein structure and will enable the construction of deletion mutants for functional studies.     

The length of the downstream exon and the substitution of specific sequences affect pre-mRNA splicing in vitro.

We have shown previously that truncation of the human beta-globin pre-mRNA in the second exon, 14 nucleotides downstream from the 3' splice site, leads to inhibition of splicing but not cleavage at the 5' splice site. We now show that several nonglobin sequences substituted at this site can restore splicing and that the efficiency of splicing depends on the length of the second (downstream) exon and not a specific sequence. Deletions in the first exon have no effect on the efficiency of in vitro splicing. Surprisingly, an intron fragment from the 5' region of the human or rabbit beta-globin intron 2, when placed 14 nucleotides downstream from the 3' splice site, inhibited all the steps in splicing beginning with cleavage at the 5' splice site. This result suggests that the intron 2 fragment carries a "poison" sequence that can inhibit the splicing of an upstream intron.     

Partial characterization of bovine elastin gene; comparison with the gene for human elastin

A 14 kilobase (kb) genomic clone of the gene for bovine elastin, containing exons 1 and 2, has been characterized. This clone extends about 6.5 kb in the 5' direction from the initiation codon and 978 nucleotides in the 3' direction from exon 2. The size of the first intron is about 6.4 kb. The sequence immediately 5' to the initiation codon is highly conserved between the genes for bovine and human elastins and contains a TATA box consensus sequence (ATAAA), CAAT, and Sp1 binding sites. Several putative AP-2 binding sites are also present. Comparative analysis of the sequences flanking the first exon in the genes for bovine and human elastins identified conserved sequences that may be regulatory control elements. A putative enhancer core sequence is present in the first intron of the genes for bovine and human elastins.