Analysis of the exon-intron structures of fish, amphibian, bird and mammalian hatching enzyme genes, with special reference to the intron loss evolution of hatching enzyme genes in Teleostei

Using gene cloning and in silico cloning, we analyzed the structures of hatching enzyme gene orthologs of vertebrates. Comparison led to a hypothesis that hatching enzyme genes of Japanese eel conserve an ancestral structure of the genes of fishes, amphibians, birds and mammals. However, the exon-intron structure of the genes was different from species to species in Teleostei: Japanese eel hatching enzyme genes were 9-exon-8-intron genes, and zebrafish genes were 5-exon-4-intron genes. In the present study, we further analyzed the gene structures of fishes belonging to Acanthopterygii. In the species of Teleostei we examined, diversification of hatching enzyme gene into two paralogous genes for HCE (high choriolytic enzyme) and LCE (low choriolytic enzyme) was found only in the acanthopterygian fishes such as medaka Oryzias latipes, Fundulus heteroclitus, Takifugu rubripes and Tetraodon nigroviridis. In addition, the HCE gene had no intron, while the LCE gene consisted of 8 exons and 7 introns. Phylogenetic analysis revealed that HCE and LCE genes were paralogous to each other, and diverged during the evolutionary lineage to Acanthopterygii. Analysis of gene synteny and cluster structure showed that the syntenic genes around the HCE and LCE genes were highly conserved between medaka and Teraodon, but such synteny was not found around the zebrafish hatching enzyme genes. We hypothesize that the zebrafish hatching enzyme genes were translocated from chromosome to chromosome, and lost some of their introns during evolution.     

Evolution of teleostean hatching enzyme genes and their paralogous genes

We isolated genes for hatching enzymes and their paralogs having two cysteine residues at their N-terminal regions in addition to four cysteines conserved in all the astacin family proteases. Genes for such six-cysteine-containing astacin proteases (C6AST) were searched out in the medaka genome database. Five genes for MC6AST1 to 5 were found in addition to embryo-specific hatching enzyme genes. RT-PCR and whole-mount in situ hybridization evidenced that MC6AST1 was expressed in embryos and epidermis of almost all adult tissues examined, while MC6AST2 and 3 were in mesenterium, intestine, and testis. MC6AST4 and 5 were specifically expressed in jaw. In addition, we cloned C6AST cDNA homologs from zebrafish, ayu, and fugu. The MC6AST1 to 5 genes were classified into three groups in the phylogenetic positions, and the expression patterns and hatching enzymes were clearly discriminated from other C6ASTs. Analysis of the exon–intron structures clarified that genes for hatching enzymes MHCE and MAHCE were intron-less, while other MC6AST genes were basically the same as the gene for another hatching enzyme MLCE. In the basal Teleost, the C6AST genes having the ancestral exon–intron structure (nine exon/eight intron structure) first appeared by duplication and chromosomal translocation. Thereafter, maintaining such ancestral exon–intron structure, the LCE gene was newly diversified in Euteleostei, and the MC6AST1 to 5 gene orthologs were duplicated and diversified independently in respective fish lineages. The HCE gene lost all introns in Euteleostei, whereas in the lineage to zebrafish, it was translocated from chromosome to chromosome and lost some of its introns.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.The nucleotide sequence data reported in the present paper will appear in the DDBJ/EMBL/GenBank nucleotide sequence databases with accession numbers from AB256940 to AB256952.     

Spliceosomal intron insertions in genome compacted ray-finned fishes as evident from phylogeny of MC receptors, also supported by a few other GPCRs

Background

Insertions of spliceosomal introns are very rare events during evolution of vertebrates and the mechanisms governing creation of novel intron(s) remain obscure. Largely, gene structures of melanocortin (MC) receptors are characterized by intron-less architecture. However, recently a few exceptions have been reported in some fishes. This warrants a systematic survey of MC receptors for understanding intron insertion events during vertebrate evolution.

Methodology/Principal Findings

We have compiled an extended list of MC receptors from different vertebrate genomes with variations in fishes. Notably, the closely linked MC2Rs and MC5Rs from a group of ray-finned fishes have three and one intron insertion(s), respectively, with conserved positions and intron phase. In both genes, one novel insertion was in the highly conserved DRY motif at the end of helix TM3. Further, the proto-splice site MAG↑R is maintained at intron insertion sites in these two genes. However, the orthologs of these receptors from zebrafish and tetrapods are intron-less, suggesting these introns are simultaneously created in selected fishes. Surprisingly, these novel introns are traceable only in four fish genomes. We found that these fish genomes are severely compacted after the separation from zebrafish. Furthermore, we also report novel intron insertions in P2Y receptors and in CHRM3. Finally, we report ultrasmall introns in MC2R genes from selected fishes.

Conclusions/Significance

The current repository of MC receptors illustrates that fishes have no MC3R ortholog. MC2R, MC5R, P2Y receptors and CHRM3 have novel intron insertions only in ray-finned fishes that underwent genome compaction. These receptors share one intron at an identical position suggestive of being inserted contemporaneously. In addition to repetitive elements, genome compaction is now believed to be a new hallmark that promotes intron insertions, as it requires rapid DNA breakage and subsequent repair processes to gain back normal functionality.     

Sense-antisense gene overlap is probably a cause for retaining the few introns in <Emphasis Type="Italic">Giardia</Emphasis> genome and the implications

Background

It is widely accepted that the last eukaryotic common ancestor and early eukaryotes were intron-rich and intron loss dominated subsequent evolution, thus the presence of only very few introns in some modern eukaryotes must be the consequence of massive loss. But it is striking that few eukaryotes were found to have completely lost introns. Despite extensive research, the causes of massive intron losses remain elusive. Actually the reverse question -- how the few introns can be retained under the evolutionary selection pressure of intron loss -- is equally significant but was rarely studied, except that it was conjectured that the essential functions of some introns prevent their loss. The situation that extremely few (eight) spliceosome-mediated cis-spliced introns present in the relatively simple genome of Giardia lamblia provides an excellent opportunity to explore this question.

Results

Our investigation found three types of distribution patterns of the few introns in the intron-containing genes: ancient intron in ancient gene, later-evolved intron in ancient gene, and later-evolved intron in later-evolved gene, which can reflect to some extent the dynamic evolution of introns in Giardia. Without finding any special features or functional importance of these introns responsible for their retention, we noticed and experimentally verified that some intron-containing genes form sense-antisense gene pairs with transcribable genes on their complementary strands, and that the introns just reside in the overlapping regions.

Conclusions

In Giardia’s evolution, despite constant evolutionary selection pressure of intron loss, intron gain can still occur in both ancient and later-evolved genes, but only a few introns are retained; at least the evolutionary retention of some of the introns might not be due to the functional constraint of the introns themselves but the causes outside of introns, such as the constraints imposed by other genomic functional elements overlapping with the introns. These findings can not only provide some clues to find new genomic functional elements -- in the areas overlapping with introns, but suggest that “functional constraint” of introns may not be necessarily directly associated with intron loss and gain, and that the real functions are probably still outside of our current knowledge.

Reviewers

This article was reviewed by Mikhail Gelfand, Michael Gray, and Igor Rogozin.

     

The exon/intron organization of the globin gene of Scapharca inaequivalvis homodimeric hemoglobin: unusual intron homology with other bivalve mollusc globin genes

In this study, we have investigated the positions of introns in the globin gene of Scapharca inaequivalvis homodimeric hemoglobin. We found the three exon/two intron organization typical of vertebrate globin genes, with the two introns in highly conserved positions, as it occurs in the A and B globin genes of the tetrameric hemoglobin from the same organism, confirming the absence of the so-called `central intron' found in the globin genes of plants and of some invertebrates. We identified two homodimeric globin genes (3207 and 2723 bp) that differ only with respect to the size of the first intron. Sequence analysis of the two first introns (1668 and 1364 bp) has revealed that they are highly homologous, except for a 569- and 296-bp insertion in each intron I. Interestingly, the two first introns contain regions with an unusually high identity (∼80%) with regions of the first intron of the congeneric clam Anadara trapezia and the related clam Barbatia reveana globin genes, suggesting that these uncoding regions may have played a regulatory role that has subsequently been lost during the course of the evolution.     

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.     

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.     

Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns

As part of the exploratory sequencing program Génolevures, visual scrutinisation and bioinformatic tools were used to detect spliceosomal introns in seven hemiascomycetous yeast species. A total of 153 putative novel introns were identified. Introns are rare in yeast nuclear genes (<5% have an intron), mainly located at the 5′ end of ORFs, and not highly conserved in sequence. They all share a clear non-random vocabulary: conserved splice sites and conserved nucleotide contexts around splice sites. Homologues of metazoan snRNAs and putative homologues of SR splicing factors were identified, confirming that the spliceosomal machinery is highly conserved in eukaryotes. Several introns’ features were tested as possible markers for phylogenetic analysis. We found that intron sizes vary widely within each genome, and according to the phylogenetic position of the yeast species. The evolutionary origin of spliceosomal introns was examined by analysing the degree of conservation of intron positions in homologous yeast genes. Most introns appeared to exist in the last common ancestor of present day yeast species, and then to have been differentially lost during speciation. However, in some cases, it is difficult to exclude a possible sliding event affecting a pre-existing intron or a gain of a novel intron. Taken together, our results indicate that the origin of spliceosomal introns is complex within a given genome, and that present day introns may have resulted from a dynamic flux between intron conservation, intron loss and intron gain during the evolution of hemiascomycetous yeasts.     

Analysis of Fungal Genomes Reveals Commonalities of Intron Gain or Loss and Functions in Intron-Poor Species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.     

Diversity of sponge mitochondrial introns revealed by <Emphasis Type="Italic">cox 1</Emphasis>sequences of Tetillidae

Background  

Animal mitochondrial introns are rare. In sponges and cnidarians they have been found in the cox 1 gene of some spirophorid and homosclerophorid sponges, as well as in the cox 1 and nad 5 genes of some Hexacorallia. Their sporadic distribution has raised a debate as to whether these mobile elements have been vertically or horizontally transmitted among their hosts. The first sponge found to possess a mitochondrial intron was a spirophorid sponge from the Tetillidae family. To better understand the mode of transmission of mitochondrial introns in sponges, we studied cox 1 intron distribution among representatives of this family.     

Endogenous retroviruses of the chicken genome

 

We report that the positions of minor, U12 introns are conserved in orthologous genes from human and Arabidopsis to an even greater extent than the positions of the major, U2 introns. The U12 introns, especially, conserved ones are concentrated in 5'-portions of plant and animal genes, where the U12 to U2 conversions occurs preferentially in the 3'-portions of genes. These results are compatible with the hypothesis that the high level of conservation of U12 intron positions and their persistence in genomes despite the unidirectional U12 to U2 conversion are explained by the role of the slowly excised U12 introns in down-regulation of gene expression.

Reviewers

This article was reviewed by John Logsdon and Manyuan Long. For the full reviews, please go to the Reviewers' Reports section.     

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

Statistics of the exon-intron structure and splicing sites of several diverse eukaryotes was studied. The yeast exon-intron structures have a number of unique features. A yeast gene usually have at most one intron. The branch site is strongly conserved, whereas the polypirimidine tract is short. Long yeast introns tend to have stronger acceptor sites. In other species the branch site is less conserved and often cannot be determined. In non-yeast samples there is an almost universal correlation between lengths of neighboring exons (all samples excluding protists) and correlation between lengths of neighboring introns (human, drosophila, protists). On the average first introns are longer, and anomalously long introns are usually first introns in a gene. There is a universal preference for exons and exon pairs with the (total) length divisible by 3. Introns positioned between codons are preferred, whereas those positioned between the first and second positions in codon are avoided. The choice of A or G at the third position of intron (the donor splice sites generally prefer purines at this position) is correlated with the overall GC-composition of the gene. In all samples dinucleotide AG is avoided in the region preceding the acceptor site.     

ConservedPrimers 2.0: A high-throughput pipeline for comparative genome referenced intron-flanking PCR primer design and its application in wheat SNP discovery

Background  

In some genomic applications it is necessary to design large numbers of PCR primers in exons flanking one or several introns on the basis of orthologous gene sequences in related species. The primer pairs designed by this target gene approach are called "intron-flanking primers" or because they are located in exonic sequences which are usually conserved between related species, "conserved primers". They are useful for large-scale single nucleotide polymorphism (SNP) discovery and marker development, especially in species, such as wheat, for which a large number of ESTs are available but for which genome sequences and intron/exon boundaries are not available. To date, no suitable high-throughput tool is available for this purpose.     

Hard tick calreticulin (CRT) gene coding regions have only one intron with conserved positions and variable sizes

Calreticulin (CRT) is a unique eukaryotic gene. The CRT gene product, calreticulin, was first identified as a calcium binding protein in 1974, but further investigations have indicated that CRT protein performs many functions in cells, including involvement in evading the host's immune system by parasites. Many studies of CRT have been published since the molecule was first discovered; however, the CRT gene exon-intron structure is only known for a limited number of ectoparasite species. In this study, we compared tick CRT genomic sequences to the corresponding cDNA from 28 species and found that 2 exons and 1 intron are present in the tick CRT gene. The intron position is conserved in 28 hard ticks, but intron size and nucleotide sequences vary. Three tick introns possess duplicated fragments and are twice as long as other introns. All tick CRT introns obey the GT-AG rule in the splice-site junctions and are phase 1 introns. By comparing tick CRT introns to those of fruit fly, mouse, and human, we conclude that tick CRT introns belong to the intron-late type. The number and size of CRT introns have increased through the evolution of eukaryotes.