On the ancient nature of introns

Exon Elongation Added Intrinsically Disordered Regions to the Encoded Proteins and Facilitated the Emergence of the Last Eukaryotic Common Ancestor

Most prokaryotic proteins consist of a single structural domain (SD) with little intrinsically disordered regions (IDRs) that by themselves do not adopt stable structures, whereas the typical eukaryotic protein comprises multiple SDs and IDRs. How eukaryotic proteins evolved to differ from prokaryotic proteins has not been fully elucidated. Here, we found that the longer the internal exons are, the more frequently they encode IDRs in eight eukaryotes including vertebrates, invertebrates, a fungus, and plants. Based on this observation, we propose the "small bang" model from the proteomic viewpoint: the protoeukaryotic genes had no introns and mostly encoded one SD each, but a majority of them were subsequently divided into multiple exons (step 1). Many exons unconstrained by SDs elongated to encode IDRs (step 2). The elongated exons encoding IDRs frequently facilitated the acquisition of multiple SDs to make the last common ancestor of eukaryotes (step 3). One prediction of the model is that long internal exons are mostly unconstrained exons. Analytical results of the eight eukaryotes are consistent with this prediction. In support of the model, we identified cases of internal exons that elongated after the rat-mouse divergence and discovered that the expanded sections are mostly in unconstrained exons and preferentially encode IDRs. The model also predicts that SDs followed by long internal exons tend to have other SDs downstream. This prediction was also verified in all the eukaryotic species analyzed. Our model accounts for the dichotomy between prokaryotic and eukaryotic proteins and proposes a selective advantage conferred by IDRs.

Evidence for exon shuffling is sensitive to model choice

The exon shuffling theory posits that intronic recombination creates new domain combinations, facilitating the evolution of novel protein function. This theory predicts that introns will be preferentially situated near domain boundaries. Many studies have sought evidence for exon shuffling by testing the correspondence between introns and domain boundaries against chance intron positioning. Here, we present an empirical investigation of how the choice of null model influences significance. Although genome-wide studies have used a uniform null model, exclusively, more realistic null models have been proposed for single gene studies. We extended these models for genome-wide analyses and applied them to 21 metazoan and fungal genomes. Our results show that compared with the other two models, the uniform model does not recapitulate genuine exon lengths, dramatically underestimates the probability of chance agreement, and overestimates the significance of intron-domain correspondence by as much as 100 orders of magnitude. Model choice had much greater impact on the assessment of exon shuffling in fungal genomes than in metazoa, leading to different evolutionary conclusions in seven of the 16 fungal genomes tested. Genome-wide studies that use this overly permissive null model may exaggerate the importance of exon shuffling as a general mechanism of multidomain evolution.

The TIM Barrel Architecture Facilitated the Early Evolution of Protein-Mediated Metabolism

The triosephosphate isomerase (TIM) barrel protein fold is a structurally repetitive architecture that is present in approximately 10 % of all enzymes. It is generally assumed that this ubiquity in modern proteomes reflects an essential historical role in early protein-mediated metabolism. Here, we provide quantitative and comparative analyses to support several hypotheses about the early importance of the TIM barrel architecture. An information theoretical analysis of protein structures supports the hypothesis that the TIM barrel architecture could arise more easily by duplication and recombination compared to other mixed α/β structures. We show that TIM barrel enzymes corresponding to the most taxonomically broad superfamilies also have the broadest range of functions, often aided by metal and nucleotide-derived cofactors that are thought to reflect an earlier stage of metabolic evolution. By comparison to other putatively ancient protein architectures, we find that the functional diversity of TIM barrel proteins cannot be explained simply by their antiquity. Instead, the breadth of TIM barrel functions can be explained, in part, by the incorporation of a broad range of cofactors, a trend that does not appear to be shared by proteins in general. These results support the hypothesis that the simple and functionally general TIM barrel architecture may have arisen early in the evolution of protein biosynthesis and provided an ideal scaffold to facilitate the metabolic transition from ribozymes, peptides, and geochemical catalysts to modern protein enzymes.

Modeling the evolution dynamics of exon-intron structure with a general random fragmentation process

Background

Most eukaryotic genes are interrupted by spliceosomal introns. The evolution of exon-intron structure remains mysterious despite rapid advance in genome sequencing technique. In this work, a novel approach is taken based on the assumptions that the evolution of exon-intron structure is a stochastic process, and that the characteristics of this process can be understood by examining its historical outcome, the present-day size distribution of internal translated exons (exon). Through the combination of simulation and modeling the size distribution of exons in different species, we propose a general random fragmentation process (GRFP) to characterize the evolution dynamics of exon-intron structure. This model accurately predicts the probability that an exon will be split by a new intron and the distribution of novel insertions along the length of the exon.

Results

As the first observation from this model, we show that the chance for an exon to obtain an intron is proportional to its size to the 3rd power. We also show that such size dependence is nearly constant across gene, with the exception of the exons adjacent to the 5^′ UTR. As the second conclusion from the model, we show that intron insertion loci follow a normal distribution with a mean of 0.5 (center of the exon) and a standard deviation of 0.11. Finally, we show that intron insertions within a gene are independent of each other for vertebrates, but are more negatively correlated for non-vertebrate. We use simulation to demonstrate that the negative correlation might result from significant intron loss during evolution, which could be explained by selection against multi-intron genes in these organisms.

Conclusions

The GRFP model suggests that intron gain is dynamic with a higher chance for longer exons; introns are inserted into exons randomly with the highest probability at the center of the exon. GRFP estimates that there are 78 introns in every 10 kb coding sequences for vertebrate genomes, agreeing with empirical observations. GRFP also estimates that there are significant intron losses in the evolution of non-vertebrate genomes, with extreme cases of around 57% intron loss in Drosophila melanogaster, 28% in Caenorhabditis elegans, and 24% in Oryza sativa.

Whence genes in pieces: reconstruction of the exon-intron gene structures of the last eukaryotic common ancestor and other ancestral eukaryotes

In eukaryotes, protein-coding sequences are interrupted by non-coding sequences known as introns. During mRNA maturation, introns are excised by the spliceosome and the coding regions, exons, are spliced to form the mature coding region. The intron densities widely differ between eukaryotic lineages, from 6 to 7 introns per kb of coding sequence in vertebrates, some invertebrates and green plants, to only a few introns across the entire genome in many unicellular eukaryotes. Evolutionary reconstructions using maximum likelihood methods suggest intron-rich ancestors for each major group of eukaryotes. For the last common ancestor of animals, the highest intron density of all extant and extinct eukaryotes was inferred, at 120-130% of the human intron density. Furthermore, an intron density within 53-74% of the human values was inferred for the last eukaryotic common ancestor. Accordingly, evolution of eukaryotic genes in all lines of descent involved primarily intron loss, with substantial gain only at the bases of several branches including plants and animals. These conclusions have substantial biological implications indicating that the common ancestor of all modern eukaryotes was a complex organism with a gene architecture resembling those in multicellular organisms. Alternative splicing most likely initially appeared as an inevitable result of splicing errors and only later was employed to generate structural and functional diversification of proteins.

DsHsp90 is involved in the early response of Dunaliella salina to environmental stress

Heat shock protein 90 (Hsp90) is a molecular chaperone highly conserved across the species from prokaryotes to eukaryotes. Hsp90 is essential for cell viability under all growth conditions and is proposed to act as a hub of the signaling network and protein homeostasis of the eukaryotic cells. By interacting with various client proteins, Hsp90 is involved in diverse physiological processes such as signal transduction, cell mobility, heat shock response and osmotic stress response. In this research, we cloned the dshsp90 gene encoding a polypeptide composed of 696 amino acids from the halotolerant unicellular green algae Dunaliella salina. Sequence alignment indicated that DsHsp90 belonged to the cytosolic Hsp90A family. Further biophysical and biochemical studies of the recombinant protein revealed that DsHsp90 possessed ATPase activity and existed as a dimer with similar percentages of secondary structures to those well-studied Hsp90As. Analysis of the nucleotide sequence of the cloned genomic DNA fragment indicated that dshsp90 contained 21 exons interrupted by 20 introns, which is much more complicated than the other plant hsp90 genes. The promoter region of dshsp90 contained putative cis-acting stress responsive elements and binding sites of transcriptional factors that respond to heat shock and salt stress. Further experimental research confirmed that dshsp90 was upregulated quickly by heat and salt shock in the D. salina cells. These findings suggested that dshsp90 might serve as a component of the early response system of the D. salina cells against environmental stresses.

Origin and evolution of spliceosomal introns

Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded ‘introns first’ held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers’ Reports section.

Split genes: another surprise from Giardia

Some genes in the candidate early-branching eukaryote Giardia lamblia occur in separate pieces, transcribed from non-contiguous chromosomal locations. The pre-mRNAs from the separate pieces apparently find each other by regions of complementarity and are subsequently spliced together by the spliceosome. Could genes in pieces, transcribed into separate pre-mRNAs, have been an early feature of spliceosomal evolution?

Evolution of spliceosomal introns following endosymbiotic gene transfer

Background

Spliceosomal introns are an ancient, widespread hallmark of eukaryotic genomes. Despite much research, many questions regarding the origin and evolution of spliceosomal introns remain unsolved, partly due to the difficulty of inferring ancestral gene structures. We circumvent this problem by using genes originated by endosymbiotic gene transfer, in which an intron-less structure at the time of the transfer can be assumed.

Results

By comparing the exon-intron structures of 64 mitochondrial-derived genes that were transferred to the nucleus at different evolutionary periods, we can trace the history of intron gains in different eukaryotic lineages. Our results show that the intron density of genes transferred relatively recently to the nuclear genome is similar to that of genes originated by more ancient transfers, indicating that gene structure can be rapidly shaped by intron gain after the integration of the gene into the genome and that this process is mainly determined by forces acting specifically on each lineage. We analyze 12 cases of mitochondrial-derived genes that have been transferred to the nucleus independently in more than one lineage.

Conclusions

Remarkably, the proportion of shared intron positions that were gained independently in homologous genes is similar to that proportion observed in genes that were transferred prior to the speciation event and whose shared intron positions might be due to vertical inheritance. A particular case of parallel intron gain in the nad7 gene is discussed in more detail.

Intron-dominated genomes of early ancestors of eukaryotes

Evolutionary reconstructions using maximum likelihood methods point to unexpectedly high densities of introns in protein-coding genes of ancestral eukaryotic forms including the last common ancestor of all extant eukaryotes. Combined with the evidence of the origin of spliceosomal introns from invading Group II self-splicing introns, these results suggest that early ancestral eukaryotic genomes consisted of up to 80% sequences derived from Group II introns, a much greater contribution of introns than that seen in any extant genome. An organism with such an unusual genome architecture could survive only under conditions of a severe population bottleneck.

Domain mobility in proteins: functional and evolutionary implications

A substantial fraction of eukaryotic proteins contains multiple domains, some of which show a tendency to occur in diverse domain architectures and can be considered mobile (or 'promiscuous'). These promiscuous domains are typically involved in protein-protein interactions and play crucial roles in interaction networks, particularly those contributing to signal transduction. They also play a major role in creating diversity of protein domain architecture in the proteome. It is now apparent that promiscuity is a volatile and relatively fast-changing feature in evolution, and that only a few domains retain their promiscuity status throughout evolution. Many such domains attained their promiscuity status independently in different lineages. Only recently, we have begun to understand the diversity of protein domain architectures and the role the promiscuous domains play in evolution of this diversity. However, many of the biological mechanisms of protein domain mobility remain shrouded in mystery. In this review, we discuss our present understanding of protein domain promiscuity, its evolution and its role in cellular function.

Molecular cloning, partial genomic structure and functional characterization of succinic semialdehyde dehydrogenase genes from the parasitic insects Lucilia cuprina and Ctenocephalides felis

The enzyme succinic semialdehyde dehydrogenase (SSADH; EC1.2.1.24) is a component of the gamma-aminobutyric acid degradation pathway in mammals and is essential for development and function of the nervous system. Here we report the identification, cDNA cloning and functional expression of SSADH from the parasitic insects Lucilia cuprina and Ctenocephalides felis. The recombinant proteins possess potent NAD+-dependent SSADH activity, while their catalytic efficiency for other aldehyde substrates is lower. A genomic copy of the L. cuprina SSADH gene contains two introns, while a genomic gene version of C. felis is devoid of introns. In contrast to the single copy SSADH genes in Drosophila melanogaster and mammals, in L. cuprina and C. felis, multiple SSADH gene copies are present in the genome.

Mechanisms of intron gain and loss in Cryptococcus

BACKGROUND

Genome comparisons across deep phylogenetic divergences have revealed that spliceosomal intron gain and loss are common evolutionary events. However, because of the deep divergences involved in these comparisons, little is understood about how these changes occur, particularly in the case of intron gain. To ascertain mechanisms of intron gain and loss, we compared five relatively closely related genomes from the yeast Cryptococcus.

RESULTS

We observe a predominance of intron loss over gain and identify a relatively slow intron loss rate in Cryptococcus. Some genes preferentially lose introns and a large proportion of intron losses occur in the middle of genes (so called internal intron loss). Finally, we identify a gene that displays a differential number of introns in a repetitive DNA region.

CONCLUSION

Based the observed patterns of intron loss and gain, population resequencing and population genetic analysis, it appears that recombination causes the widely observed but poorly understood phenomenon of internal intron loss and that DNA repeat expansion can create new introns in a population.

Conserved intron positions in ancient protein modules

Background

The timing of the origin of introns is of crucial importance for an understanding of early genome architecture. The Exon theory of genes proposed a role for introns in the formation of multi-exon proteins by exon shuffling and predicts the presence of conserved splice sites in ancient genes. In this study, large-scale analysis of potential conserved splice sites was performed using an intron-exon database (ExInt) derived from GenBank.

Results

A set of conserved intron positions was found by matching identical splice sites sequences from distantly-related eukaryotic kingdoms. Most amino acid sequences with conserved introns were homologous to consensus sequences of functional domains from conserved proteins including kinases, phosphatases, small GTPases, transporters and matrix proteins. These included ancient proteins that originated before the eukaryote-prokaryote split, for instance the catalytic domain of protein phosphatase 2A where a total of eleven conserved introns were found. Using an experimental setup in which the relation between a splice site and the ancientness of its surrounding sequence could be studied, it was found that the presence of an intron was positively correlated to the ancientness of its surrounding sequence. Intron phase conservation was linked to the conservation of the gene sequence and not to the splice site sequence itself. However, no apparent differences in phase distribution were found between introns in conserved versus non-conserved sequences.

Conclusion

The data confirm an origin of introns deep in the eukaryotic branch and is in concordance with the presence of introns in the first functional protein modules in an 'Exon theory of genes' scenario. A model is proposed in which shuffling of primordial short exonic sequences led to the formation of the first functional protein modules, in line with hypotheses that see the formation of introns integral to the origins of genome evolution.

Reviewers

This article was reviewed by Scott Roy (nominated by Anthony Poole), Sandro de Souza (nominated by Manyuan Long), and Gáspár Jékely.

Exon size distribution and the origin of introns

Since it was first recognised that eukaryotic genes are fragmented into coding segments (exons) separated by non-coding segments (introns), the reason for this phenomenon has been debated. There are two dominant theories: that the piecewise arrangement of genes allows functional protein domains, represented by exons, to recombine by shuffling to form novel proteins with combinations of functions; or that introns represent parasitic DNA that can infest the eukaryotic genome because it does not interfere grossly with the fitness of its host. Differing distributions of exon lengths are predicted by these two theories. In this paper we examine distributions of exon lengths for six different organisms and find that they offer empirical evidence that both theories may in part be correct.

A glimpse of a putative pre-intron phase of eukaryotic evolution

Comparison of the exon-intron structures of ancient eukaryotic paralogs reveals the absence of conserved intron positions in these genes. This is in contrast to the conservation of intron positions in orthologous genes from even the most evolutionarily distant eukaryotes and in more recent paralogs. The lack of conserved intron positions in ancient paralogs probably reflects the origination of these genes during the earliest phase of eukaryotic evolution, which was characterized by concomitant invasion of genes by group II self-splicing elements (which were to become introns in the future) and extensive duplication of genes.

Evolution of secretin family GPCR members in the metazoa

Background

Comparative approaches using protostome and deuterostome data have greatly contributed to understanding gene function and organismal complexity. The family 2 G-protein coupled receptors (GPCRs) are one of the largest and best studied hormone and neuropeptide receptor families. They are suggested to have arisen from a single ancestral gene via duplication events. Despite the recent identification of receptor members in protostome and early deuterostome genomes, relatively little is known about their function or origin during metazoan divergence. In this study a comprehensive description of family 2 GPCR evolution is given based on in silico and expression analyses of the invertebrate receptor genes.

Results

Family 2 GPCR members were identified in the invertebrate genomes of the nematodes C. elegans and C. briggsae, the arthropods D. melanogaster and A. gambiae (mosquito) and in the tunicate C. intestinalis. This suggests that they are of ancient origin and have evolved through gene/genome duplication events. Sequence comparisons and phylogenetic analyses have demonstrated that the immediate gene environment, with regard to gene content, is conserved between the protostome and deuterostome receptor genomic regions. Also that the protostome genes are more like the deuterostome Corticotrophin Releasing Factor (CRF) and Calcitonin/Calcitonin Gene-Related Peptide (CAL/CGRP) receptors members than the other family 2 GPCR members. The evolution of family 2 GPCRs in deuterostomes is characterised by acquisition of new family members, with SCT (Secretin) receptors only present in tetrapods. Gene structure is characterised by an increase in intron number with organismal complexity with the exception of the vertebrate CAL/CGRP receptors.

Conclusion

The family 2 GPCR members provide a good example of gene duplication events occurring in tandem with increasing organismal complexity during metazoan evolution. The putative ancestral receptors are proposed to be more like the deuterostome CAL/CGRP and CRF receptors and this may be associated with their fundamental role in calcium regulation and the stress response, both of which are essential for survival.

Origins and Evolution of Spliceosomal Introns

Research into the origins of introns is at a critical juncture in the resolution of theories on the evolution of early life (which came first, RNA or DNA?), the identity of LUCA (the last universal common ancestor, was it prokaryotic- or eukaryotic-like?), and the significance of noncoding nucleotide variation. One early notion was that introns would have evolved as a component of an efficient mechanism for the origin of genes. But alternative theories emerged as well. From the debate between the "introns-early" and "introns-late" theories came the proposal that introns arose before the origin of genetically encoded proteins and DNA, and the more recent "introns-first" theory, which postulates the presence of introns at that early evolutionary stage from a reconstruction of the "RNA world." Here we review seminal and recent ideas about intron origins. Recent discoveries about the patterns and causes of intron evolution make this one of the most hotly debated and exciting topics in molecular evolutionary biology today.

The origin of introns and their role in eukaryogenesis: a compromise solution to the introns-early versus introns-late debate?

BACKGROUND

Ever since the discovery of 'genes in pieces' and mRNA splicing in eukaryotes, origin and evolution of spliceosomal introns have been considered within the conceptual framework of the 'introns early' versus 'introns late' debate. The 'introns early' hypothesis, which is closely linked to the so-called exon theory of gene evolution, posits that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. Under this scenario, the absence of spliceosomal introns in prokaryotes is considered to be a result of "genome streamlining". The 'introns late' hypothesis counters that spliceosomal introns emerged only in eukaryotes, and moreover, have been inserted into protein-coding genes continuously throughout the evolution of eukaryotes. Beyond the formal dilemma, the more substantial side of this debate has to do with possible roles of introns in the evolution of eukaryotes.

RESULTS

I argue that several lines of evidence now suggest a coherent solution to the introns-early versus introns-late debate, and the emerging picture of intron evolution integrates aspects of both views although, formally, there seems to be no support for the original version of introns-early. Firstly, there is growing evidence that spliceosomal introns evolved from group II self-splicing introns which are present, usually, in small numbers, in many bacteria, and probably, moved into the evolving eukaryotic genome from the alpha-proteobacterial progenitor of the mitochondria. Secondly, the concept of a primordial pool of 'virus-like' genetic elements implies that self-splicing introns are among the most ancient genetic entities. Thirdly, reconstructions of the ancestral state of eukaryotic genes suggest that the last common ancestor of extant eukaryotes had an intron-rich genome. Thus, it appears that ancestors of spliceosomal introns, indeed, have existed since the earliest stages of life's evolution, in a formal agreement with the introns-early scenario. However, there is no evidence that these ancient introns ever became widespread before the emergence of eukaryotes, hence, the central tenet of introns-early, the role of introns in early evolution of proteins, has no support. However, the demonstration that numerous introns invaded eukaryotic genes at the outset of eukaryotic evolution and that subsequent intron gain has been limited in many eukaryotic lineages implicates introns as an ancestral feature of eukaryotic genomes and refutes radical versions of introns-late. Perhaps, most importantly, I argue that the intron invasion triggered other pivotal events of eukaryogenesis, including the emergence of the spliceosome, the nucleus, the linear chromosomes, the telomerase, and the ubiquitin signaling system. This concept of eukaryogenesis, in a sense, revives some tenets of the exon hypothesis, by assigning to introns crucial roles in eukaryotic evolutionary innovation.

CONCLUSION

The scenario of the origin and evolution of introns that is best compatible with the results of comparative genomics and theoretical considerations goes as follows: self-splicing introns since the earliest stages of life's evolution--numerous spliceosomal introns invading genes of the emerging eukaryote during eukaryogenesis--subsequent lineage-specific loss and gain of introns. The intron invasion, probably, spawned by the mitochondrial endosymbiont, might have critically contributed to the emergence of the principal features of the eukaryotic cell. This scenario combines aspects of the introns-early and introns-late views.

REVIEWERS

this article was reviewed by W. Ford Doolittle, James Darnell (nominated by W. Ford Doolittle), William Martin, and Anthony Poole.

Intron gain and loss in segmentally duplicated genes in rice

BACKGROUND

Introns are under less selection pressure than exons, and consequently, intronic sequences have a higher rate of gain and loss than exons. In a number of plant species, a large portion of the genome has been segmentally duplicated, giving rise to a large set of duplicated genes. The recent completion of the rice genome in which segmental duplication has been documented has allowed us to investigate intron evolution within rice, a diploid monocotyledonous species.

RESULTS

Analysis of segmental duplication in rice revealed that 159 Mb of the 371 Mb genome and 21,570 of the 43,719 non-transposable element-related genes were contained within a duplicated region. In these duplicated regions, 3,101 collinear paired genes were present. Using this set of segmentally duplicated genes, we investigated intron evolution from full-length cDNA-supported non-transposable element-related gene models of rice. Using gene pairs that have an ortholog in the dicotyledonous model species Arabidopsis thaliana, we identified more intron loss (49 introns within 35 gene pairs) than intron gain (5 introns within 5 gene pairs) following segmental duplication. We were unable to demonstrate preferential intron loss at the 3' end of genes as previously reported in mammalian genomes. However, we did find that the four nucleotides of exons that flank lost introns had less frequently used 4-mers.

CONCLUSION

We observed that intron evolution within rice following segmental duplication is largely dominated by intron loss. In two of the five cases of intron gain within segmentally duplicated genes, the gained sequences were similar to transposable elements.

An analysis of intron positions in relation to nucleotides, amino acids, and protein secondary structure

We present an analysis of intron positions in relation to nucleotides, amino acid residues, and protein secondary structure. Previous work has shown that intron sites in proteins are not randomly distributed with respect to secondary structures. Here we show that this preference can be almost totally explained by the nucleotide bias of splice site machinery, and may well not relate to protein stability or conformation at all. Each intron phase is preferentially associated with its own set of residues: phase 0 introns with lysine, glutamine, and glutamic acid before the intron, and valine after; phase 1 introns with glycine, alanine, valine, aspartic acid, and glutamic acid; and phase 2 introns with arginine, serine, lysine, and tryptophan. These preferences can be explained principally on the basis of nucleotide bias at intron locations, which is in accordance with previous literature. Although this work does not prove that introns are inserted into genomes at specific proto-splice sites, it shows that the nucleotide bias surrounding introns, however it originally occurred, explains the observed correlations between introns and protein secondary structure.

Organisation of trypsin genes in the salmon louse (Lepeophtheirus salmonis, Crustacea, copepoda) genome

Trypsins constitute a subclass of the S1A family of serine peptidases found in all groups of animal and some bacteria. At present, no information about the genomic organisation of trypsins is available for copepods. The only data of copepod trypsins indicate several different trypsins in the marine parasitic copepod Lepeophtheirus salmonis. In the present study, 31.7 kbp of genomic DNA surrounding the previously described LsTryp1-5 sequences was sequenced. The sequenced regions contain nine full-length and three partial trypsin genes. A conservative estimate based on PCR analysis and genomic sequence indicated at least 22 different trypsin genes in L. salmonis, of which 18 are most similar to the previously described LsTryp1 and -2 cDNA sequences. Four of these genes are putative pseudogenes. In addition, a putative mariner like transposase gene was identified. The genomic sequences suggest that the L. salmonis trypsin genes reside within one or more gene clusters. Three different LsTryp intron exon structures were identified, and all three are different from the intron exon organisation previously reported for other S1A peptidases. This implies several intron loss and gain events in the evolution of the L. salmonis trypsin genes.

Heterogeneity of intron presence or absence in rDNA genes of the lichen species Physcia aipolia and P. stellaris

Intron origin and evolution are of high interest, yet the rates of insertion and loss are unclear. To investigate their spread, we studied ribosomal (r)DNA introns from the closely related lichens Physcia aipolia and P. stellaris. Both taxa are replete with rDNA spliceosomal introns and autocatalytic group I introns, many of which show presence/absence polymorphism when screened with the PCR approach. This initially suggested that Physcia could be a model for studying intron retention and loss. However, during the course of a population-level analysis, we discovered widespread intron presence/absence heterogeneity within lichen thalli. To address this result, we sequenced multiple clones encoding nuclear rDNA and the single-copy elongation factor-1alpha (EF-1alpha) from individual thalli. These data showed extensive rDNA heterogeneity within individuals, rather than the presence of multiple fungi within a thallus. Our results suggest that considerable care must be taken when interpreting intron presence/absence in lichen rDNA, an observation that has general implications for the study of rDNA intron evolution.

Preferential loss and gain of introns in 3' portions of genes suggests a reverse-transcription mechanism of intron insertion

In an attempt to gain insight into the dynamics of intron evolution in eukaryotic protein-coding genes, the distributions of old introns, that are conserved between distant phylogenetic lineages, and new, lineage-specific introns along the gene length, were examined. A significant excess of old introns in 5'-regions of genes was detected. New introns, when analyzed in bulk, showed a nearly flat distribution from the 5'- to the 3'-end. However, analysis of new intron distributions in individual genomes revealed notable lineage-specific features. While in intron-poor genomes, particularly yeast Schizosaccharomyces pombe (Sp), the 5'-portions of genes contain a significantly greater number of new introns than the 3'-portions, the intron-rich genomes of humans and Arabidopsis show the opposite trend. These observations seem to be compatible with the view that introns are both lost and inserted in 3'-terminal portions of genes more often than in 5'-portions. Overrepresentation of 3'-terminal sequences among cDNAs that mediate intron loss appears to be the most likely explanation for the apparent preferential loss of introns in the distal parts of genes. Preferential insertion of introns in the 3'-portions suggests that introns might be inserted via a reverse-transcription-mediated pathway similar to that implicated in intron loss. This mechanism could involve duplication of a portion of the coding region during reverse transcription followed by homologous recombination and subsequent rapid sequence divergence in the copy that becomes a new intron.

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.

Remarkable interkingdom conservation of intron positions and massive, lineage-specific intron loss and gain in eukaryotic evolution

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.

The unusual phylogenetic distribution of retrotransposons: a hypothesis

Retrotransposons have proliferated extensively in eukaryotic lineages; the genomes of many animals and plants comprise 50% or more retrotransposon sequences by weight. There are several persuasive arguments that the enzymatic lynchpin of retrotransposon replication, reverse transcriptase (RT), is an ancient enzyme. Moreover, the direct progenitors of retrotransposons are thought to be mobile self-splicing introns that actively propagate themselves via reverse transcription, the group II introns, also known as retrointrons. Retrointrons are represented in modern genomes in very modest numbers, and thus far, only in certain eubacterial and organellar genomes. Archaeal genomes are nearly devoid of RT in any form. In this study, I propose a model to explain this unusual distribution, and rationalize it with the proposed ancient origin of the RT gene. A cap and tail hypothesis is proposed. By this hypothesis, the specialized terminal structures of eukaryotic mRNA provide the ideal molecular environment for the lengthening, evolution, and subsequent massive expansion of highly mobile retrotransposons, leading directly to the retrotransposon-cluttered structure that typifies modern metazoan genomes and the eventual emergence of retroviruses.

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.

An intron-less betagamma-crystallin-type gene from the sponge Geodia cydonium

We report the cloning of a gene encoding a betagamma-crystallin-type protein from a porifera, the Geodia cydonium sponge. The data provide direct, conclusive evidence of the existence of such a gene in the genome of an early diverged metazoan. The cloned gene is found to contain no introns, while proto-splice sites are identified in the nucleotide sequence at positions where introns are located in homologous, very recently diverged vertebrate genes. These findings are discussed in the light of the debate between the introns-late and introns-early theories.

Conservation of the positions of metazoan introns from sponges to humans

Sponges (phylum Porifera) are the phylogenetic oldest Metazoa still extant. They can be considered as reference animals (Urmetazoa) for the understanding of the evolutionary processes resulting in the creation of Metazoa in general and also for the metazoan gene organization in particular. In the marine sponge Suberites domuncula, genes encoding p38 and JNK kinases contain nine and twelve introns, respectively. Eight introns in both genes share the same positions and the identical phases. One p38 intron slipped for six bases and the JNK gene has three more introns. However, the sequences of the introns are not conserved and the introns in JNK gene are generally much longer. Introns interrupt most of the conserved kinase subdomains I-XI and are found in all three phases (0, 1 and 2). We analyzed in details p38 and JNK genes from human, Caenorhabditis elegans and Drosophila melanogaster and found in most genes introns at the positions identical to those in sponge genes. The exceptions are two p38 genes from D. melanogaster that have lost all introns in the coding sequence. The positions of 11 introns in each of four human p38 genes are fully conserved and ten introns occupy identical positions as the introns in sponge p38 or JNK genes. The same is true for nine, out of ten introns in the human JNK-1 gene. The introns in human p38 and JNK genes are on average more than ten times longer than corresponding introns in sponges. It was proposed that yeast HOG1-like kinases (from i.e. Saccharomyces cerevisiae and Emericella nidulans) and metazoan p38 and JNK kinases are orthologues. p38 and JNK genes were created after the split from fungi by the duplication and diversification of the HOG1-like progenitor gene. Our results further support the common origin of p38 and JNK genes and speak in favor of a very early time of duplication. The ancestral gene contained at least ten introns, which are still present at the very conserved positions in p38 and JNK genes of extant animals. Four of these introns are present at the same positions in the HOG-like gene in the fungus E. nidulans. The others probably entered the ancestral gene after the split of fungi, but before the duplication of the gene and before the creation of the common, urmetazoan progenitor of all multicellular animals. A second gene coding for an immune molecule is described, the allograft inflammatory factor, which likewise showed a highly conserved exon/intron structure in S. domuncula and in human. These data show that the intron/exon borders are highly conserved in genes from sponges to human.

Anticipatory evolution and DNA shuffling

DNA shuffling has proven to be a powerful technique for the directed evolution of proteins. A mix of theoretical and applied research has now provided insights into how recombination can be guided to more efficiently generate proteins and even organisms with altered functions.

Comparative genomics and evolution of proteins involved in RNA metabolism

RNA metabolism, broadly defined as the compendium of all processes that involve RNA, including transcription, processing and modification of transcripts, translation, RNA degradation and its regulation, is the central and most evolutionarily conserved part of cell physiology. A comprehensive, genome-wide census of all enzymatic and non-enzymatic protein domains involved in RNA metabolism was conducted by using sequence profile analysis and structural comparisons. Proteins related to RNA metabolism comprise from 3 to 11% of the complete protein repertoire in bacteria, archaea and eukaryotes, with the greatest fraction seen in parasitic bacteria with small genomes. Approximately one-half of protein domains involved in RNA metabolism are present in most, if not all, species from all three primary kingdoms and are traceable to the last universal common ancestor (LUCA). The principal features of LUCA's RNA metabolism system were reconstructed by parsimony-based evolutionary analysis of all relevant groups of orthologous proteins. This reconstruction shows that LUCA possessed not only the basal translation system, but also the principal forms of RNA modification, such as methylation, pseudouridylation and thiouridylation, as well as simple mechanisms for polyadenylation and RNA degradation. Some of these ancient domains form paralogous groups whose evolution can be traced back in time beyond LUCA, towards low-specificity proteins, which probably functioned as cofactors for ribozymes within the RNA world framework. The main lineage-specific innovations of RNA metabolism systems were identified. The most notable phase of innovation in RNA metabolism coincides with the advent of eukaryotes and was brought about by the merge of the archaeal and bacterial systems via mitochondrial endosymbiosis, but also involved emergence of several new, eukaryote-specific RNA-binding domains. Subsequent, vast expansions of these domains mark the origin of alternative splicing in animals and probably in plants. In addition to the reconstruction of the evolutionary history of RNA metabolism, this analysis produced numerous functional predictions, e.g. of previously undetected enzymes of RNA modification.

Intron distribution difference for 276 ancient and 131 modern genes suggests the existence of ancient introns

o introns delineate elements of protein tertiary structure? This issue is crucial to the debate about the role and origin of introns. We present an analysis of the full set of proteins with known three-dimensional structures that have homologs with intron positions recorded in GenBank. A computer program was generated that maps on a reference sequence the positions of all introns in homologous genes. We have applied this program to a set of 665 nonredundant protein sequences with defined three-dimensional structures in the Protein Data Bank (PDB), which yielded 8,217 introns in 407 proteins. For the subset of proteins corresponding to ancient conserved regions (ACR), we find that there is a correlation of phase-zero introns with the boundary regions of modules and no correlation for the phase-one and phase-two positions. However, for a subset of proteins without prokaryotic counterparts (131 non-ACR proteins), a set of presumably modern proteins (or proteins that have diverged extremely far from any ancestral form), we do not find any correlation of phase-zero intron positions with three-dimensional structure. Furthermore, we find an anticorrelation of phase-one intron positions with module boundaries: they actually have a preference for the interior of modules. This finding is explicable as a preference for phase-one introns to lie in glycines, between G/G sequences, the preference for glycines being anticorrelated with the three-dimensional modules. We interpret this anticorrelation as a sign that a number of phase-one introns, and hence many modern introns, have been inserted into G/G "protosplice" sequences.

Scale-free behavior in protein domain networks

Several technical, social, and biological networks were recently found to demonstrate scale-free and small-world behavior instead of random graph characteristics. In this work, the topology of protein domain networks generated with data from the ProDom, Pfam, and Prosite domain databases was studied. It was found that these networks exhibited small-world and scale-free topologies with a high degree of local clustering accompanied by a few long-distance connections. Moreover, these observations apply not only to the complete databases, but also to the domain distributions in proteomes of different organisms. The extent of connectivity among domains reflects the evolutionary complexity of the organisms considered.

Class I release factors in ciliates with variant genetic codes

In eukaryotes with the universal genetic code a single class I release factor (eRF1) most probably recognizes all stop codons (UAA, UAG and UGA) and is essential for termination of nascent peptide synthesis. It is well established that stop codons have been reassigned to amino acid codons at least three times among ciliates. The codon specificities of ciliate eRF1s must have been modified to accommodate the variant codes. In this study we have amplified, cloned and sequenced eRF1 genes of two hypotrichous ciliates, Oxytricha trifallax (UAA and UAG for Gln) and Euplotes aediculatus (UGA for Cys). We also sequenced/identified three protist and two archaeal class I RF genes to enlarge the database of eRF1/aRF1s with the universal code. Extensive comparisons between universal code eRF1s and those of Oxytricha, Euplotes, and Tetrahymena which represent three lineages that acquired variant codes independently, provide important clues to identify stop codon-binding regions in eRF1. Domain 1 in the five ciliate eRF1s, particularly the TASNIKS heptapeptide and its adjacent region, differs significantly from domain 1 in universal code eRF1s. This observation suggests that domain 1 contains the codon recognition site, but that the mechanism of eRF1 codon recognition may be more complex than proposed by Nakamura et al. or Knight and Landweber.

About the last common ancestor, the universal life-tree and lateral gene transfer: a reappraisal

An organismal tree rooted in the bacterial branch and derived from a hyperthermophilic last common ancestor (LCA) is still widely assumed to represent the path followed by evolution from the most primeval cells to the three domains recognized among contemporary organisms: Bacteria, Archaea and Eucarya. In the past few years, however, more and more discrepancies between this pattern and individual protein trees have been brought to light. There has been an overall tendency to attribute these incongruities to widespread lateral gene transfer. However, recent developments, a reappraisal of earlier evidence and considerations of our own lead us to a quite different view. It would appear (i) that the role of lateral gene transfer was overemphasized in recent discussions of molecular phylogenies; (ii) that the LCA was probably a non-thermophilic protoeukaryote from which both Archaea and Bacteria emerged by reductive evolution but not as sister groups, in keeping with a current evolutionary scheme for the biosynthesis of membrane lipids; and (iii) that thermophilic Archaea may have been the first branch to diverge from the ancestral line.

Intron-exon structure and gene copy number of a gene encoding for a membrane-intrinsic light-harvesting polypeptide of the red alga Galdieria sulphuraria

Genes for light-harvesting proteins (lhc genes) of higher plants are well examined. However, little is known about the corresponding genes of algae, although this knowledge might give valuable information about the evolution of photosynthetic antennae. In the case of rhodophytes only two cDNA sequences from a single organism, Porphyridium cruentum, have been published. Here we describe an additional sequence from another species, the thermo-acidophilic red alga Galdieria sulphuraria. For the first time also a genomic sequence for a red algal lhc gene is presented. From a cDNA library of G. sulphuraria we isolated a clone containing an open reading frame for a protein of 302 amino acids with a deduced molecular mass of 33.86kDa. It shares major structural features with eukaryotic light-harvesting polypeptides. A proposed cleavage site between transit peptide and mature protein gives rise to a transit peptide of 119 amino acids and a mature protein of 183 residues. Hydropathy analysis suggests that the mature protein consists of three transmembrane helices. Several amino acid residues supposed to bind chlorophyll a and chlorophyll b in higher plants are conserved. The protein shows up to 69% identity and 81% similarity to the Porphyridium polypeptides in the transmembrane helices 1 and 3. Using oligonucleotides annealing in the regions of the start and stop codons of the gene as primers, a DNA sequence was amplified from nuclear G. sulphuraria DNA by PCR. Compared with the cDNA clone, this sequence contains five additional intervening DNA strings of 50-74bp length. Four of them show typical features of spliceosomal introns with GT-AG borders, and the fifth differs by starting with GC. Three of the supposed introns are located in similar positions as introns of higher plant light-harvesting proteins. Southern blotting and hybridization experiments indicate that G. sulphuraria contains at least three copies of this gene.

Early evolution: prokaryotes, the new kids on the block

Prokaryotes are generally assumed to be the oldest existing form of life on earth. This assumption, however, makes it difficult to understand certain aspects of the transition from earlier stages in the origin of life to more complex ones, and it does not account for many apparently ancient features in the eukaryotes. From a model of the RNA world, based on relic RNA species in modern organisms, one can infer that there was an absolute requirement for a high-accuracy RNA replicase even before proteins evolved. In addition, we argue here that the ribosome (together with the RNAs involved in its assembly) is so large that it must have had a prior function before protein synthesis. A model that connects and equates these two requirements (high-accuracy RNA replicase and prior function of the ribosome) can explain many steps in the origin of life while accounting for the observation that eukaryotes have retained more vestiges of the RNA world. The later derivation of prokaryote RNA metabolism and genome structure can be accounted for by the two complementary mechanisms of r-selection and thermoreduction.

Genome evolution and the evolution of exon-shuffling--a review

Recent studies on the genomes of protists, plants, fungi and animals confirm that the increase in genome size and gene number in different eukaryotic lineages is paralleled by a general decrease in genome compactness and an increase in the number and size of introns. It may thus be predicted that exon-shuffling has become increasingly significant with the evolution of larger, less compact genomes. To test the validity of this prediction, we have analyzed the evolutionary distribution of modular proteins that have clearly evolved by intronic recombination. The results of this analysis indicate that modular multidomain proteins produced by exon-shuffling are restricted in their evolutionary distribution. Although such proteins are present in all major groups of metazoa from sponges to chordates, there is practically no evidence for the presence of related modular proteins in other groups of eukaryotes. The biological significance of this difference in the composition of the proteomes of animals, fungi, plants and protists is best appreciated when these modular proteins are classified with respect to their biological function. The majority of these proteins can be assigned to functional categories that are inextricably linked to multicellularity of animals, and are of absolute importance in permitting animals to function in an integrated fashion: constituents of the extracellular matrix, proteases involved in tissue remodelling processes, various proteins of body fluids, membrane-associated proteins mediating cell-cell and cell-matrix interactions, membrane associated receptor proteins regulating cell cell communications, etc. Although some basic types of modular proteins seem to be shared by all major groups of metazoa, there are also groups of modular proteins that appear to be restricted to certain evolutionary lineages. In summary, the results suggest that exon-shuffling acquired major significance at the time of metazoan radiation. It is interesting to note that the rise of exon-shuffling coincides with a spectacular burst of evolutionary creativity: the Big Bang of metazoan radiation. It seems probable that modular protein evolution by exon-shuffling has contributed significantly to this accelerated evolution of metazoa, since it facilitated the rapid construction of multidomain extracellular and cell surface proteins that are indispensable for multicellularity.

A mini-protein designed by removing a module from barnase: molecular modeling and NMR measurements of the conformation

A globular domain can be decomposed into compact modules consisting of contiguous 10-30 amino acid residues. The correlation between modules and exons observed in different proteins suggests that each module was encoded by an ancestral exon and that modules were combined into globular domains by exon fusion. Barnase is a single domain RNase consisting of 110 amino acid residues and was decomposed into six modules. We designed a mini-protein by removing the second module, M2, from barnase in order to gain an insight into the structural and functional roles of the module. In the molecular modeling of the mini-protein, we evaluated thermodynamic stability and aqueous solubility together with mechanical stability of the model. We chemically synthesized a mini-barnase with (15)N-labeling at 10 residues, whose corresponding residues in barnase are all found in the region around the hydrophobic core. Circular dichroism and NMR measurements revealed that mini-barnase takes a non-random specific conformation that has a similar hydrophobic core structure to that of barnase. This result, that a module could be deleted without altering the structure of core region of barnase, supports the view that modules act as the building blocks of protein design.

Structural, functional, and genetic characterization of Gastrophilus hemoglobin

Hemoglobin of Gastrophilus intestinalis (Insecta, Diptera), was purified and characterized. At least two isoforms have been identified by isoelectrofocusing, mass spectrometry, and genomic Southern blotting. Functional studies show a high oxygen affinity due to a low ligand dissociation rate (koff = 2.4 s-1) and a relatively high autoxidation rate (t1/2 = 1.6/h). The globins were separated under denaturing conditions, and the sequence of Hb1 (Mr = 17,965 +/- 2) was determined at the protein and DNA level. The open reading frame codes for a polypeptide of 150 amino acids. Although the globin is distantly related to globins from other species, it has a low penalty score against globin templates. Freshly isolated hemoglobin was crystallized from polyethylene glycol. Crystals contain two hemoglobin molecules per asymmetric unit. Solution of the three-dimensional structure by molecular replacement could not be achieved, possibly due to the presence of three protein isoforms in the crystals. In order to determine its three-dimensional structure, G. intestinalis Hb1 was overexpressed in Escherichia coli, resulting in a fully functional molecule as confirmed by ligand binding affinity. The globin gene contains two introns at positions D7.0 and G7.0. The D7.0 intron is unprecedented, suggesting that globin gene evolution is much more complex than originally thought.

The recent origins of spliceosomal introns revisited

Does the intron/exon structure of eukaryotic genes belie their ancient assembly by exon-shuffling or have introns been inserted into preformed genes during eukaryotic evolution? These are the central questions in the ongoing 'introns-early' versus 'introns-late' controversy. The phylogenetic distribution of spliceosomal introns continues to strongly favor the intronslate theory. The introns-early theory, however, has claimed support from intron phase and protein structure correlations.

Characterization of the myoglobin and its coding gene of the mollusc Biomphalaria glabrata

A cDNA clone isolated from a Biomphalaria glabrata (Mollusca, Gastropoda) neural cDNA library was identified as encoding a myoglobin-like protein of 148 amino acids with a single domain and a calculated mass of 16,049.29. Alignment with globin sequences with known tertiary structure confirms its overall globin nature. The expressed myoglobin was identified in the radular muscle and isolated. Oxygen equilibrium measurements on the protein reveal a high oxygen affinity. Val-B10 and Gln-E7, important residues for the determination of the oxygen affinity, are strikingly different from the standard molluscan pattern (Conti, E., Moser, C., Rizzi, M., Mattevi, A., Lionetti, C., Coda, A., Ascenzi, P., Brunori, M., Bolognesi, M. (1993) J. Mol. Biol. 233, 498-508). The single gene encoding the globin chain is interrupted by three introns at positions A3.2, B12.2, and G7.0. Comparison with other nonvertebrate globin genes reveals on the one hand conservation (B12.2 and G7.0) and on the other hand variability of the insertion positions (A3.2). The Biomphalaria myoglobin sequence was used together with all other molluscan globin sequences available to assess the origin and phylogeny of the phylum. Our results confirm the doubts raised about monophyletic origin of the Mollusca, which was first observed using SSU rRNA as a molecular marker.

Structure and chromosomal localization of mouse G protein subunit gamma 4 gene

The G protein gamma subunits are members of a multigene family and are implicated in determining the specificity of receptor-G protein interaction. The gene structures for many of the gamma subunits remain to be determined. Here, we report the gene structure for the brain-specific gamma 4 subunit and its map position on a mouse chromosome. The gene (Gng4) comprises at least three exons spanning over 20 kb. The 225-bp coding region, which spans two exons, is interrupted by a large 18.2-kb intron whose position is conserved in other gamma subunit genes. There is a putative additional intron in the 5' untranslated region just upstream of the translation initiation codon. Introns are present in most of the other gamma subunits at this position. The mouse Gng4 gene is mapped to chromosome 13.

De novo insertion of an intron into the mammalian sex determining gene, SRY

Two theories have been proposed to explain the evolution of introns within eukaryotic genes. The introns early theory, or "exon theory of genes," proposes that introns are ancient and that recombination within introns provided new exon structure, and thus new genes. The introns late theory, or "insertional theory of introns," proposes that ancient genes existed as uninterrupted exons and that introns have been introduced during the course of evolution. There is still controversy as to how intron-exon structure evolved and whether the majority of introns are ancient or novel. Although there is extensive evidence in support of the introns early theory, phylogenetic comparisons of several genes indicate recent gain and loss of introns within these genes. However, no example has been shown of a protein coding gene, intronless in its ancestral form, which has acquired an intron in a derived form. The mammalian sex determining gene, SRY, is intronless in all mammals studied to date, as is the gene from which it recently evolved. However, we report here comparisons of genomic and cDNA sequences that now provide evidence of a de novo insertion of an intron into the SRY gene of dasyurid marsupials. This recently (approximately 45 million years ago) inserted sequence is not homologous with known transposable elements. Our data demonstrate that introns may be inserted as spliced units within a developmentally crucial gene without disrupting its function.