Statistical analysis of the exon-intron structure of higher and lower eukaryote genes

Selection for reduced translation costs at the intronic 5' end in fungi

It is generally believed that introns are not translated; therefore, the potential intronic features that may be related to the translation step (occurring after splicing) have yet to be thoroughly studied. Here, focusing on four fungi, we performed for the first time a comprehensive study aimed at characterizing how translation efficiency is encoded in introns and affects their evolution. By analysing their intronome we provide evidence of selection for STOP codons close to the intronic 5′ end, and show that the beginning of introns are selected for significantly high translation, presumably to reduce translation and metabolic costs in cases of non-spliced introns. Ribosomal profiling data analysis in Saccharomyces cerevisiae supports the conjecture that in this organism intron retention frequently occurs, introns are partially translated, and their translation efficiency affects organismal fitness. We show that the reported results are more significant in highly translated and highly spliced genes, but are not associated only with genes with a specific function. We also discuss the potential relation of the reported signals to efficient nonsense-mediated decay due to splicing errors. These new discoveries are supported by population-genetics considerations. In addition, they are contributory steps towards a broader understanding of intron evolution and the effect of silent mutations on gene expression and organismal fitness.

Nucleotide sequence composition adjacent to intronic splice sites improves splicing efficiency via its effect on pre-mRNA local folding in fungi

RNA splicing is the central process of intron removal in eukaryotes known to regulate various cellular functions such as growth, development, and response to external signals. The canonical sequences indicating the splicing sites needed for intronic boundary recognition are well known. However, the roles and evolution of the local folding of intronic and exonic sequence features adjacent to splice sites has yet to be thoroughly studied. Here, focusing on four fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, and Candida albicans), we performed for the first time a comprehensive high-resolution study aimed at characterizing the encoding of intronic splicing efficiency in pre-mRNA transcripts and its effect on intron evolution. Our analysis supports the conjecture that pre-mRNA local folding strength at intronic boundaries is under selective pressure, as it significantly affects splicing efficiency. Specifically, we show that in the immediate region of 12-30 nucleotides (nt) surrounding the intronic donor site there is a preference for weak pre-mRNA folding; similarly, in the region of 15-33 nt surrounding the acceptor and branch sites there is a preference for weak pre-mRNA folding. We also show that in most cases there is a preference for strong pre-mRNA folding further away from intronic splice sites. In addition, we demonstrate that these signals are not associated with gene-specific functions, and they correlate with splicing efficiency measurements (r = 0.77, P = 2.98 × 10(-21)) and with expression levels of the corresponding genes (P = 1.24 × 10(-19)). We suggest that pre-mRNA folding strength in the above-mentioned regions has a direct effect on splicing efficiency by improving the recognition of intronic boundaries. These new discoveries are contributory steps toward a broader understanding of splicing regulation and intronic/transcript evolution.

Origins and impacts of new mammalian exons

Mammalian genes are composed of exons, but the evolutionary origins and functions of new internal exons are poorly understood. Here, we analyzed patterns of exon gain using deep cDNA sequencing data from five mammals and one bird, identifying thousands of species- and lineage-specific exons. Most new exons derived from unique rather than repetitive intronic sequence. Unlike exons conserved across mammals, species-specific internal exons were mostly located in 5' UTRs and alternatively spliced. They were associated with upstream intronic deletions, increased nucleosome occupancy, and RNA polymerase II pausing. Genes containing new internal exons had increased gene expression, but only in tissues in which the exon was included. Increased expression correlated with the level of exon inclusion, promoter proximity, and signatures of cotranscriptional splicing. Altogether, these findings suggest that increased splicing at the 5' ends of genes enhances expression and that changes in 5' end splicing alter gene expression between tissues and between species.

Conservation/Mutation in the splice sites of cytokine receptor genes of mouse and human

Conservation/mutation in the intronic initial and terminal hexanucleotides was studied in 26 orthologous cytokine receptor genes of Mouse and Human. Introns began and ended with the canonical dinucleotides GT and AG, respectively. Identical configurations were found in 57% of the 5′ hexanucleotides and 28% of the 3′ hexanucleotides. The actual conservation percentages of the individual variable nucleotides at each position in the hexanucleotides were determined, and the theoretical rates of conservation of groups of three nucleotides were calculated under the hypothesis of a mutual evolutionary independence of the neighboring nucleotides (random association). Analysis of the actual conservation of groups of variable nucleotides showed that, at 5′, GTGAGx was significantly more expressed and GTAAGx was significantly less expressed, as compared to the random association. At 3′, TTTxAG and xTGCAG were overexpressed as compared to a random association. Study of Mouse and Human transcript variants involving the splice sites showed that most variants were not inherited from the common ancestor but emerged during the process of speciation. In some variants the silencing of a terminal hexanucleotide determined skipping of the downstream exon; in other variants the constitutive splicing hexanucleotide was replaced by another potential, in-frame, splicing hexanucleotide, leading to alterations of exon lengths.

Aberrantly spliced HTT, a new player in Huntington's disease pathogenesis

Huntington's disease (HD) is an adult-onset neurodegenerative disorder caused by a mutated CAG repeat in the huntingtin gene that is translated into an expanded polyglutamine tract. The clinical manifestation of HD is a progressive physical, cognitive, and psychiatric deterioration that is eventually fatal. The mutant huntingtin protein is processed into several smaller fragments, which have been implicated as critical factors in HD pathogenesis. The search for proteases responsible for their production has led to the identification of several cleavage sites on the huntingtin protein. However, the origin of the small N-terminal fragments that are found in HD postmortem brains has remained elusive. Recent mapping of huntingtin fragments in a mouse model demonstrated that the smallest N-terminal fragment is an exon 1 protein. This discovery spurred our hypothesis that mis-splicing as opposed to proteolysis could be generating the smallest huntingtin fragment. We demonstrated that mis-splicing of mutant huntingtin intron 1 does indeed occur and results in a short polyadenylated mRNA, which is translated into an exon 1 protein. The exon 1 protein fragment is highly pathogenic. Transgenic mouse models containing just human huntingtin exon 1 develop a rapid onset of HD-like symptoms. Our finding that a small, mis-spliced HTT transcript and corresponding exon 1 protein are produced in the context of an expanded CAG repeat has unraveled a new molecular mechanism in HD pathogenesis. Here we present detailed models of how mis-splicing could be facilitated, what challenges remain in this model, and implications for therapeutic studies.

Genomic organization and promoter cloning of the human X11α gene APBA1

X11α is a brain specific multi-modular protein that interacts with the Alzheimer's disease amyloid precursor protein (APP). Aggregation of amyloid-β peptide (Aβ), an APP cleavage product, is believed to be central to the pathogenesis of Alzheimer's disease. Recently, overexpression of X11α has been shown to reduce Aβ generation and to ameliorate memory deficit in a transgenic mouse model of Alzheimer's disease. Therefore, manipulating the expression level of X11α may provide a novel route for the treatment of Alzheimer's disease. Human X11α is encoded by the gene APBA1. As evidence suggests that X11α expression can be regulated at transcription level, we have determined the gene structure and cloned the promoter of APBA1. APBA1 spans over 244 kb on chromosome 9 and is composed of 13 exons and has multiple transcription start sites. A putative APBA1 promoter has been identified upstream of exon 1 and functional analysis revealed that this is highly active in neurons. By deletion analysis, the minimal promoter was found to be located between -224 and +14, a GC-rich region that contains a functional Sp3 binding site. In neurons, overexpression of Sp3 stimulates the APBA1 promoter while an Sp3 inhibitor suppresses the promoter activity. Moreover, inhibition of Sp3 reduces endogenous X11α expression and promotes the generation of Aβ. Our findings reveal that Sp3 play an essential role in APBA1 transcription.

Longer first introns are a general property of eukaryotic gene structure

While many properties of eukaryotic gene structure are well characterized, differences in the form and function of introns that occur at different positions within a transcript are less well understood. In particular, the dynamics of intron length variation with respect to intron position has received relatively little attention. This study analyzes all available data on intron lengths in GenBank and finds a significant trend of increased length in first introns throughout a wide range of species. This trend was found to be even stronger when using high-confidence gene annotation data for three model organisms (Arabidopsis thaliana, Caenorhabditis elegans, and Drosophila melanogaster) which show that the first intron in the 5' UTR is--on average--significantly longer than all downstream introns within a gene. A partial explanation for increased first intron length in A. thaliana is suggested by the increased frequency of certain motifs that are present in first introns. The phenomenon of longer first introns can potentially be used to improve gene prediction software and also to detect errors in existing gene annotations.

Selection for the miniaturization of highly expressed genes

Most widely expressed genes are also highly expressed. Based on high or wide expression, different models were proposed to explain the small sizes of highly/widely expressed genes. We found that housekeeping genes are not more compact than narrowly expressed genes with similar expression levels, but compactness and expression level are correlated in housekeeping genes (except that highly expressed Arabidopsis HK genes have longer intron length). Meanwhile, we found evidence that genes with high functional/regulatory complexity do not have longer introns and longer proteins. The genome design hypothesis is thus not supported. Furthermore, we found that housekeeping genes are not more compact than the narrowly expressed somatic genes with similar average expression levels. Because housekeeping genes are expected to have much higher germline expression levels than narrowly expressed somatic genes, transcription-associated deletion bias is not supported. Selection of the compactness of highly expressed genes for economy is supported.

HybGFS: a hybrid method for genome-fingerprint scanning

Background

Protein identification based on mass spectrometry (MS) has previously been performed using peptide mass fingerprinting (PMF) or tandem MS (MS/MS) database searching. However, these methods cannot identify proteins that are not already listed in existing databases. Moreover, the alternative approach of de novo sequencing requires costly equipment and the interpretation of complex MS/MS spectra. Thus, there is a need for novel high-throughput protein-identification methods that are independent of existing predefined protein databases.

Results

Here, we present a hybrid method for genome-fingerprint scanning, known as HybGFS. This technique combines genome sequence-based peptide MS/MS ion searching with liquid-chromatography elution-time (LC-ET) prediction, to improve the reliability of identification. The hybrid method allows the simultaneous identification and mapping of proteins without a priori information about their coding sequences. The current study used standard LC-MS/MS data to query an in silico-generated six-reading-frame translation and the enzymatic digest of an entire genome. Used in conjunction with precursor/product ion-mass searching, the LC-ETs increased confidence in the peptide-identification process and reduced the number of false-positive matches. The power of this method was demonstrated using recombinant proteins from the Escherichia coli K12 strain.

Conclusion

The novel hybrid method described in this study will be useful for the large-scale experimental confirmation of genome coding sequences, without the need for transcriptome-level expression analysis or costly MS database searching.

Effect of 5'UTR introns on gene expression in Arabidopsis thaliana

BACKGROUND

The majority of introns in gene transcripts are found within the coding sequences (CDSs). A small but significant fraction of introns are also found to reside within the untranslated regions (5'UTRs and 3'UTRs) of expressed sequences. Alignment of the whole genome and expressed sequence tags (ESTs) of the model plant Arabidopsis thaliana has identified introns residing in both coding and non-coding regions of the genome.

RESULTS

A bioinformatic analysis revealed some interesting observations: (1) the density of introns in 5'UTRs is similar to that in CDSs but much higher than that in 3'UTRs; (2) the 5'UTR introns are preferentially located close to the initiating ATG codon; (3) introns in the 5'UTRs are, on average, longer than introns in the CDSs and 3'UTRs; and (4) 5'UTR introns have a different nucleotide composition to that of CDS and 3'UTR introns. Furthermore, we show that the 5'UTR intron of the A. thaliana EF1alpha-A3 gene affects the gene expression and the size of the 5'UTR intron influences the level of gene expression.

CONCLUSION

Introns within the 5'UTR show specific features that distinguish them from introns that reside within the coding sequence and the 3'UTR. In the EF1alpha-A3 gene, the presence of a long intron in the 5'UTR is sufficient to enhance gene expression in plants in a size dependent manner.

The presence of GC-AG introns in Neurospora crassa and other euascomycetes determined from analyses of complete genomes: implications for automated gene prediction

A combination of experimental and computational approaches was employed to identify introns with noncanonical GC-AG splice sites (GC-AG introns) within euascomycete genomes. Evaluation of 2335 cDNA-confirmed introns from Neurospora crassa revealed 27 such introns (1.2%). A similar frequency (1.0%) of GC-AG introns was identified in Fusarium graminearum, in which 3 of 292 cDNA-confirmed introns contained GC-AG splice sites. Computational analyses of the N. crassa genome using a GC-AG intron consensus sequence identified an additional 20 probable GC-AG introns in this fungus. For 8 of the 47 GC-AG introns identified in N. crassa a GC donor site is also present in a homolog from Magnaporthe grisea, F. graminearum, or Aspergillus nidulans. In most cases, however, homologs in these fungi contain a GT-AG intron or no intron at the corresponding position. These findings have important implications for fungal genome annotation, as the automated annotations of euascomycete genomes incorrectly identified intron boundaries for all of the confirmed and probable GC-AG introns reported here.

The Drosophila homologue of the dystrophin gene - introns containing promoters are the major contributors to the large size of the gene

We show that the drosophila gene encoding the dystrophin-like protein (DLP) is as complex as the mammalian dystrophin gene. Three 5' promoters and three internal promoters regulate the expression of three full-length and three truncated products, respectively. The existence of this complex gene structure in such evolutionary remote organisms suggests that both types of products have diverse important functions. The promoters of both the DLP gene and the mammalian dystrophin gene are located in very large introns. These introns contribute significantly to the large size of the genes. The possible relevance of the conservation of the large size of introns containing promoters to the regulation of promoter activity is discussed.

Gene organization features in A/T-rich organisms

Several species have genomes in which the four nucleotides are not equally represented (Glöckner 2000). Interestingly, shifts to very high A/T or G/C levels can occur in several distinct branches of the tree of life. The underlying reasons for these shifts therefore may be of different origin. Now entire chromosome sequences from two different A/T-rich genomes, Dictyostelium discoideum and Plasmodium falciparum, are available (Bowman et al. 1999; Gardner et al. 2002; Glöckner et al. 2002). This gives us the opportunity to investigate how a high A/T content may influence the signals that are the landmarks for gene specification. We found that, in contrast with most known metazoan and plant genomes, splice signals contain, little information other than the canonical GT-AG dinucleotides. Intron lengths in A/T rich organisms, on the other hand, are comparable to those of other lower eukaryotes. Intergenic regions show, dependent on the orientation of adjacent genes, a size pattern with a ratio of 1 (3'-3') to 2 (3'-5') to 3 (5'-5'). Overall, gene organization patterns seem not to be influenced by the A/T bias. Surprisingly, the slightly higher A/T content of the P. falciparum genome compared to that of D. discoideum (80.1 versus 77.4%) is not achieved by increased A/T richness in intergenic regions. Instead both the shift of the nucleotide usage in coding regions to A/T-rich codons and the longer intergenic regions make an equal contribution to the higher A/T content in this organism.

Genomic organization of the cadmium-inducible tandem repeat 25-kDa metallothionein of the oligochaete worm Enchytraeus buchholzi

The terrestric oligochaete worm Enchytraeus buchholzi survives in cadmium (Cd)-polluted environments by aid of its Cd-inducible 25 kDa cysteine-rich protein (CRP). Here, we analyze promoter and structure of the crp gene and compare its relationship to MT genes. The crp gene, approximately 12 kbp long, consists of 10 exons with exons 2 to 9 encoding eight almost identical repeats of predominantly 31 amino acids of the CRP. The introns of the crp gene contain various repetitive elements including retrotransposon-like sequences. The 683-bp promoter of the non-constitutive crp gene exhibits a much higher basal activity than the mouse MT-II promoter in HepG2 cells. Essential for crp promoter activity is the distal region (-683/-521) with a GC box and the proximal region (-308/-8) with the four MREa, b, c, d and AP-1, -2, -3 elements, whereas the central portion (-521/-309) with CAAT box, CRE and a XRE causes promoter repression. The TATA box-, MREc- and the AP-2, -3-containing region are required for high crp promoter activity. Our data support the view that the crp gene is a unique MT-gene and has evolved by exon duplications from a MT-like ancestral gene.

Gopinath G, Graudens E, Hahn Y, Han M, Han ZG, Hanada K, Hanaoka H, Harada E, Hashimoto K, Hinz U, Hirai M, Hishiki T, Hopkinson I, Imbeaud S, Inoko H, Kanapin A, Kaneko Y, Kasukawa T, Kelso J, Kersey P, Kikuno R, Kimura K, Korn B, Kuryshev V, Makalowska I, Makino T, Mano S, Mariage-Samson R, Mashima J, Matsuda H, Mewes HW, Minoshima S, Nagai K, Nagasaki H, Nagata N, Nigam R, Ogasawara O, Ohara O, Ohtsubo M, Okada N, Okido T, Oota S, Ota M, Ota T, Otsuki T, Piatier-Tonneau D, Poustka A, Ren SX, Saitou N, Sakai K, Sakamoto S, Sakate R, Schupp I, Servant F, Sherry S, Shiba R, Shimizu N, Shimoyama M, Simpson AJ, Soares B, Steward C, Suwa M, Suzuki M, Takahashi A, Tamiya G, Tanaka H, Taylor T, Terwilliger JD, Unneberg P, Veeramachaneni V, Watanabe S, Wilming L, Yasuda N, Yoo HS, Stodolsky M, Makalowski W, Go M, Nakai K, Takagi T, Kanehisa M, Sakaki Y, Quackenbush J, Okazaki Y, Hayashizaki Y, Hide W, Chakraborty R, Nishikawa K, Sugawara H, Tateno Y, Chen Z, Oishi M, Tonellato P, Apweiler R, Okubo K, Wagner L, Wiemann S, Strausberg RL, Isogai T, Auffray C, Nomura N, Gojobori T, Sugano S. Integrative annotation of 21,037 human genes validated by full-length cDNA clones

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.

Xpro: database of eukaryotic protein-encoding genes

Xpro is a relational database that contains all the eukaryotic protein-encoding DNA sequences contained in GenBank with associated data required for the analysis of eukaryotic gene architecture. In addition to the information found in the GenBank records, which includes properties such as sequence, position, length and description about introns, exons and protein-coding regions, Xpro provides annotations on the splice sites and intron phases. Furthermore, Xpro validates intron positions using alignment information between the record's sequence and EST sequences found in dbEST. In the process of validation, alternative splicing information is also obtained and can be found in the database. The intron-containing genes in the Xpro are also classified as experimental or predicted based on the intron position validation and specific keywords in the GenBank records that are present in predicted genes. An Entrez-like query system, which is familiar to most biologists, is provided for accessing the information present in the database system. A non-redundant set of Xpro database contents is also obtained by cross-referencing to the Swiss-Prot/TrEMBL and Pfam databases. The database currently contains information for 493,983 genes--351,918 intron- containing genes and 142,065 intron-less genes. Xpro is updated for each new GenBank release and is freely available via the internet at http://origin.bic. nus.edu.sg/xpro.

DNA analysis servers: plot.it, bend.it, model.it and IS

The WWW servers at http://www.icgeb.trieste.it/dna/ are dedicated to the analysis of user-submitted DNA sequences; plot.it creates parametric plots of 45 physicochemical, as well as statistical, parameters; bend.it calculates DNA curvature according to various methods. Both programs provide 1D as well as 2D plots that allow localisation of peculiar segments within the query. The server model.it creates 3D models of canonical or bent DNA starting from sequence data and presents the results in the form of a standard PDB file, directly viewable on the user's PC using any molecule manipulation program. The recently established introns server allows statistical evaluation of introns in various taxonomic groups and the comparison of taxonomic groups in terms of length, base composition, intron type etc. The options include the analysis of splice sites and a probability test for exon-shuffling.

On biased distribution of introns in various eukaryotes

We conducted comprehensive analyses on intron positions in the Mus musculus genome by comparing genomic sequences in the GenBank database and cDNA sequences in the mouse cDNA library recently developed by Riken Genomic Sciences Center. Our results confirm that introns have a tendency to be located toward the 5' end of the gene. The same type of analysis was conducted in the coding region of seven eukaryotes (Saccharomyces cerevisiae, Plasmodium falciparum, Caenorhabditis elegans, Drosophila melanogaster, M. musculus, Homo sapiens, Arabidopsis thaliana). Introns in genes with a single intron have a locational bias toward the 5' end in all species except A. thaliana. We also measured the distance from the start codon to the position of the intron, and found that single introns prefer the location immediately after the start codon in S. cerevisiae and P. falciparum. We discuss three possible explanations for these findings: (1) they are the consequence of intron loss by reverse-transcriptase; (2) they are necessary to accommodate the function; and (3) they are concerned with the mechanism of pre-mRNA splicing.

Polymorphism Patterns in Two Tightly Linked Developmental Genes, Idgf1 and Idgf3, of Drosophila melanogaster

A new developmental gene family, recently identified in D. melanogaster, has been called imaginal disc growth factors (IDGF) because the proteins promote growth of cell lineages derived from imaginal discs. These are the first genes reported that encode polypeptide factors with mitotic activity in invertebrates. Characteristics such as similar arrangement of introns and exons, small size, and different cytological localization make this family an excellent candidate for evolutionary studies. We focus on the loci Idgf1 and Idgf3, two genes that possess the most distinctive features. We examine the pattern of intra- and interspecific nucleotide variation in the sequences from 20 isogenic lines of D. melanogaster and sequences from D. simulans and D. yakuba. While MK, HKA, and Tajima’s tests of neutrality fail to reject a neutral model of molecular evolution, Fu and Li’s test with outgroup and McDonald’s test suggest that balancing selection is modulating the evolution of the Idgf1 locus. The rate of recombination between the two loci is high enough to uncouple any linkage disequilibrium arising between Idgf1 and Idgf3, despite their close physical proximity.

An intron-containing glycoside hydrolase family 9 cellulase gene encodes the dominant 90 kDa component of the cellulosome of the anaerobic fungus Piromyces sp. strain E2

The cellulosome produced by Piromyces sp. strain E2 during growth on filter paper was purified by using an optimized cellulose-affinity method consisting of steps of EDTA washing of the cellulose-bound protein followed by elution with water. Three dominant proteins were identified in the cellulosome preparation, with molecular masses of 55, 80 and 90 kDa. Treatment of cellulose-bound cellulosome with a number of denaturing agents was also tested. Incubation with 0.5% (w/v) SDS or 8 M urea released most cellulosomal proteins, while leaving the greater fraction of the 80, 90 and 170 kDa components. To investigate the major 90 kDa cellulosome protein further, the corresponding gene, cel9A, was isolated, using immunoscreening and N-terminal sequencing. Inspection of the cel9A genomic organization revealed the presence of four introns, allowing the construction of a consensus for introns in anaerobic fungi. The 2800 bp cDNA clone contained an open reading frame of 2334 bp encoding a 757-residue extracellular protein. Cel9A includes a 445-residue glycoside hydrolase family 9 catalytic domain, and so is the first fungal representative of this large family. Both modelling of the catalytic domain as well as the activity measured with low level expression in Escherichia coli indicated that Cel9A is an endoglucanase. The catalytic domain is succeeded by a putative beta-sheet module of 160 amino acids with unknown function, followed by a threonine-rich linker and three fungal docking domains. Homology modelling of the Cel9A dockerins suggested that the cysteine residues present are all involved in disulphide bridges. The results presented here are used to discuss evolution of glycoside hydrolase family 9 enzymes.

SURVEY AND SUMMARY: exon-intron organization of genes in the slime mold Physarum polycephalum

The slime mold Physarum polycephalum is a morphologically simple organism with a large and complex genome. The exon-intron organization of its genes exhibits features typical for protists and fungi as well as those characteristic for the evolutionarily more advanced species. This indicates that both the taxonomic position as well as the size of the genome shape the exon-intron organization of an organism. The average gene has 3.7 introns which are on average 138 bp, with a rather narrow size distribution. Introns are enriched in AT base pairs by 13% relative to exons. The consensus sequences at exon-intron boundaries resemble those found for other species, with minor differences between short and long introns. A unique feature of P.polycephalum introns is the strong preference for pyrimidines in the coding strand throughout their length, without a particular enrichment at the 3'-ends.