A Group I Intron in the 18S Ribosomal DNA from the Parasitic Fungus Isaria japonica

The nucleotide sequence of the 18S rDNA coding gene in the ascomycetes parasitic fungus Isaria japonica contains a group I intron with a length of 379 nucleotides. The identification of the DNA sequence as a group I intron is based on its sequence homology to other fungal group I introns. Its group I intron contained the highly conserved sequence elements P, Q, R, and S found in other group I introns. Surprisingly, the intron sequence of I. japonica is more similar to that of Ustilago maydis than to the one found in Sclerotinia sclerotiorum. This is in contrast to the sequence identity found on the neighboring rDNA. This is an interesting finding and suggests a horizontal transfer of group I intron sequences. Received: 19 September 1997 / Accepted: 10 September 1998     

Origin and Inheritance of Group I Introns in 26S rRNA Genes of Gaeumannomyces graminis

Studies of the distribution of the three group I introns (intron A, intron T, and intron AT) in the 26S rDNA of Gaeumannomyces graminis had suggested that they were transferred to a common ancestor of G. graminis var. avenae and var. tritici after it had branched off from var. graminis. Intron AT and intron A exhibited vertical inheritance and coevolved in concert with their hosts. Intron loss could occur after its acquisition. Loss of any one of the three introns could occur in var. tritici whereas only loss of intron T had been found in the majority of var. avenae isolates. The existence of isolates of var. tritici and var. avenae with three introns suggested that intron loss could be reversed by intron acquisition and that the whole process is a dynamic one. This process of intron acquisition and intron loss reached different equilibrium points for different varieties and subgroups, which explained the irregular distribution of these introns in G. graminis. Each of the three group I introns was more closely related to other intron sequences that share the same insertion point in the 26S rDNA than to each other. These introns in distantly related organisms appeared to have a common ancestry. This system had provided a good model for studies on both the lateral transfer and common ancestry of group I introns in the 26S rRNA genes. Received: 17 May 1996 / Accepted: 14 January 1997     

Group I Introns from Zygomycota: Evolutionary Implications for the Fungal IC1 Intron Subgroup

The origins of fungal group I introns within nuclear small-subunit (nSSU) rDNA are enigmatic. This is partly because they have never been reported in basal fungal phyla (Zygomycota and Chytridiomycota), which are hypothesized to be ancestral to derived phyla (Ascomycota and Basidiomycota). Here we report group I introns from the nSSU rDNA of two zygomycete fungi, Zoophagus insidians (Zoopagales) and Coemansia mojavensis (Kickxellales). Secondary structure analyses predicted that both introns belong to the IC1 subgroup and that they are distantly related to each other, which is also suggested by different insertion sites. Molecular phylogenetic analyses indicated that the IC1 intron of Z. insidians is closely related to the IC1 intron inserted in the LSU rDNA of the basidiomycete fungus Clavicorona taxophila, which strongly suggests interphylum horizontal transfer. The IC1 intron of C. mojavensis has a low phylogenetic affinity to other fungal IC1 introns inserted into site 943 of nSSU rDNA (relative to E. coli 16S rDNA). It is noteworthy that this intron contains a putative ORF containing a His–Cys box motif in the antisense strand, a hallmark for nuclear-encoded homing endonucleases. Overall, molecular phylogenetic analyses do not support the placement of these two introns in basal fungal IC1 intron lineages. This result leads to the suggestion that fungal IC1 introns might have invaded or been transferred laterally after the divergence of the four major fungal phyla. Received: 8 February 2001 / Accepted: 1 November 2001     

Multiple origins of fungal group I introns located in the same position of nuclear SSU rRNA gene

The archiascomycetous fungus Protomyces pachydermus has two group I introns within the nuclear small subunit (nSSU) rRNA gene. One of these introns has an internal open reading frame (ORF) that encodes a predicted protein of 228 amino acid residues. On the other hand, Protomyces macrosporus has two group I introns that insert at the same positions as P. pachydermus, which have no ORF. Each alignment was constructed with Protomyces group I introns located in the same position and other introns retrieved by the BLAST Search. Each phylogenetic tree based on the alignment shows that Protomyces introns are monophyletic but the relationships among fungal introns do not reflect on the fungal phylogeny. Therefore, it is suggested that two different horizontal transfers of group I introns occurred at the early stage of Protomyces species diversification. Received: 11 June 1997 / Accepted: 2 September 1997     

Evolution of Mhc–DRB Introns: Implications for the Origin of Primates

Introns are generally believed to evolve too rapidly and too erratically to be of much use in phylogenetic reconstructions. Few phylogenetically informative intron sequences are available, however, to ascertain the validity of this supposition. In the present study the supposition was tested on the example of the mammalian class II major histocompatibility complex (Mhc) genes of the DRB family. Since the Mhc genes evolve under balancing selection and are believed to recombine or rearrange frequently, the evolution of their introns could be expected to be particularly rapid and subject to scrambling. Sequences of intron 4 and 5 DRB genes were obtained from polymerase chain reaction-amplified fragments of genomic DNA from representatives of six eutherian orders—Primates, Scandentia, Chiroptera, Dermoptera, Lagomorpha, and Insectivora. Although short stretches of the introns have indeed proved to be unalignable, the bulk of the intron sequences from all six orders, spanning >85 million years (my) of evolution, could be aligned and used in a study of the tempo and mode of intron evolution. The analysis has revealed the Mhc introns to evolve at a rate similar to that of other genes and of synonymous sites of non-Mhc genes. No evidence of homogenization or large-scale scrambling of the intron sequences could be found. The Mhc introns apparently evolve largely by point mutations and insertions/deletions. The phylogenetic signals contained in the intron sequences could be used to identify Scandentia as the sister group of Primates, to support the existence of the Archonta superorder, and to confirm the monophyly of the Chiroptera. Received: 26 October 1998 / Accepted: 21 December 1998     

Evolutionary History of MHC Class I Genes in the Mammalian Order Perissodactyla

We carried out an analysis of partial sequences from expressed major histocompatibility complex (MHC) class I genes isolated from a range of equid species and more distantly related members of the mammalian order Perissodactyla. Phylogenetic analysis revealed a minimum of six groups, five of which contained genes and alleles that are found in equid species and one group specific to the rhinoceros. Four of the groups contained only one, or very few sequences, indicating the presence of relatively nonpolymorphic loci, while another group contained the majority of the equid sequences identified. These data suggest that a diversification of MHC genes took place after the split between the Equidae and the Rhinocerotidae yet before the speciation events within the genus Equus. Received: 17 November 1998 / Accepted: 7 April 1999     

The Peperomia Mitochondrial coxI Group I Intron: Timing of Horizontal Transfer and Subsequent Evolution of the Intron

The Peperomia polybotrya coxI gene intron is the only currently reported group I intron in a vascular plant mitochondrial genome and it likely originated by horizontal transfer from a fungal donor. We provide a clearer picture of the horizontal transfer and a portrayal of the evolution of the group I intron since it was gained by the Peperomia mitochondrial genome. The intron was transferred recently in terms of plant evolution, being restricted to the single genus Peperomia among the order Piperales. Additional support is presented for the suggestion that a recombination/repair mechanism was used by the intron for integration into the Peperomia mitochondrial genome, as a perfect 1:1 correspondence exists between the intron's presence in a species and the presence of divergent nucleotide markers flanking the intron insertion site. Sequencing of coxI introns from additional Peperomia species revealed that several mutations have occurred in the intron since the horizontal transfer, but sequence alterations have not caused frameshifts or created stop codons in the intronic open reading frame. In addition, two coxI pseudogenes in Peperomia cubensis were discovered that lack a large region of coxI exon 2 and contain a truncated version of the group I intron that likely cannot be spliced out. Received: 29 May 1997 / Accepted: 1 November 1997     

Regular Spliceosomal Introns Are Invasive in Chlamydomonas reinhardtii: 15 Introns in the Recently Relocated Mitochondrial cox2 and cox3 Genes

In the unicellular green alga, Chlamydomonas reinhardtii, cytochrome oxidase subunit 2 (cox2) and 3 (cox3) genes are missing from the mitochondrial genome. We isolated and sequenced a BAC clone that carries the whole cox3 gene and its corresponding cDNA. Almost the entire cox2 gene and its cDNA were also determined. Comparison of the genomic and the corresponding cDNA sequences revealed that the cox3 gene contains as many as nine spliceosomal introns and that cox2 bears six introns. Putative mitochondria targeting signals were predicted at each N terminal of the cox genes. These spliceosomal introns were typical GT–AG-type introns, which are very common not only in Chlamydomonas nuclear genes but also in diverse eukaryotic taxa. We found no particular distinguishing features in the cox introns. Comparative analysis of these genes with the various mitochondrial genes showed that 8 of the 15 introns were interrupting the conserved mature protein coding segments, while the other 7 introns were located in the N-terminal target peptide regions. Phylogenetic analysis of the evolutionary position of C. reinhardtii in Chlorophyta was carried out and the existence of the cox2 and cox3 genes in the mitochondrial genome was superimposed in the tree. This analysis clearly shows that these cox genes were relocated during the evolution of Chlorophyceae. It is apparent that long before the estimated period of relocation of these mitochondrial genes, the cytosol had lost the splicing ability for group II introns. Therefore, at least eight introns located in the mature protein coding region cannot be the direct descendant of group II introns. Here, we conclude that the presence of these introns is due to the invasion of spliceosomal introns, which occurred during the evolution of Chlorophyceae. This finding provides concrete evidence supporting the ``intron-late' model, which rests largely on the mobility of spliceosomal introns. Received: 22 August 2000 / Accepted: 28 February 2001     

Identification of Sequences Encoding the Detoxification Metalloisomerase Glyoxalase I in Microbial Genomes from Several Pathogenic Organisms

The ubiquitous glyoxalase system, which is composed of two enzymes, removes cellular cytotoxic methylglyoxal (MG). In an effort to identify critical residues conserved in the evolution of the first enzyme in this system, glyoxalase I (GlxI), as well as the structural implications of sequence alterations in this enzyme, a search of the National Center for Biotechnology Information (NCBI) database of unfinished genomes was undertaken. Eleven putative GlxI sequences from pathogenic organisms were identified and analyses of these sequences in relation to the known and previously identified GlxI enzymes were performed. Several of these sequences show a very high similarity to the Escherichia coli GlxI sequence, most notably the 79% identity of the sequence identified from Yersinia pestis, the causative agent of bubonic plague. In addition to the conservation of residues critical to binding the catalytic metal in all of the proposed GlxI enzymes, four regions in the Homo sapiens GlxI enzyme are absent in all of the bacterial GlxI sequences, with the exception of Pseudomonas putida. Removal of these regions may alter the active-site conformation of the bacterial enzymes in relation to that of the H. sapiens. These differences may be targeted for the development of inhibitors selective to the bacterial enzymes. Received: 13 October 1999 / Accepted: 17 January 2000     

Identification of Protein-Tyrosine Phosphatases in Archaea

Protein-tyrosine dephosphorylation is a major mechanism in cellular regulation. A large number of protein-tyrosine phosphatases is known from Eukarya, and more recently bacterial homologues have also been identified. By employing conserved sequence patterns from both eukaryotic and bacterial protein-tyrosine phosphatases, we have identified three homologous sequences in two of the four complete archaeal genomes. Two hypothetical open reading frames in the genome of Methanococcus jannaschii (MJ0215 and MJECL20) and one in the genome of Pyrococcus horikoshii (PH1732) clearly bear all the conserved residues of this family. No homologues were found in the genomes of Archaeoglobus fulgidus and Methanobacterium thermoautotrophicum. This is the first report of protein-tyrosine phosphatase sequences in Archaea. Received: 29 April 1998 / Accepted: 27 November 1998     

Locus-specific conservation of the HLA class I introns by intra-locus homogenization

 The protein-coding sequences of the major histocompatibility complex (MHC) genes are characterized by extraordinarily high polymorphism, apparently maintained by balancing selection, which favors diversity in the peptide-binding domains of the MHC glycoproteins. Here we report that the introns flanking the polymorphic exons of the human MHC class I loci HLA-A, -B, and -C genes have been relatively conserved and have become locus-specific apparently as a result of recombination and subsequent genetic drift, leading to homogenization within loci over evolutionary time. Thus, HLA class I genes have been shaped by contrasting evolutionary forces maintaining polymorphism in the exons and leading to conservation in the introns. This study provides the first extensive analysis of the introns of a highly polymorphic gene family. Received: 10 April 1997 / Revised: 15 July 1997     

Intron Conservation in a UV-Specific DNA Repair Gene Encoded by Chlorella Viruses

Large dsDNA-containing chlorella viruses encode a pyrimidine dimer-specific glycosylase (PDG) that initiates repair of UV-induced pyrimidine dimers. The PDG enzyme is a homologue of the bacteriophage T4-encoded endonuclease V. The pdg gene was cloned and sequenced from 42 chlorella viruses isolated over a 12-year period from diverse geographic regions. Surprisingly, the pdg gene from 15 of these 42 viruses contain a 98-nucleotide intron that is 100% conserved among the viruses and another 4 viruses contain an 81-nucleotide intron, in the same position, that is nearly 100% identical (one virus differed by one base). In contrast, the nucleotides in the pdg coding regions (exons) from the intron-containing viruses are 84 to 100% identical. The introns in the pdg gene have 5′-AG/GTATGT and 3′-TTGCAG/AA splice site sequences which are characteristic of nuclear-located, spliceosomal processed pre-mRNA introns. The 100% identity of the 98-nucleotide intron sequence in the 15 viruses and the near-perfect identity of an 81-nucleotide intron sequence in another 4 viruses imply strong selective pressure to maintain the DNA sequence of the intron when it is in the pdg gene. However, the ability of intron-plus and intron-minus viruses to repair UV-damaged DNA in the dark was nearly identical. These findings contradict the widely accepted dogma that intron sequences are more variable than exon sequences. Received: 13 May 1999 / Accepted: 20 August 1999     

Predicting the secondary structures and tertiary interactions of 211 group I introns in IE subgroup

The large number of currently available group I intron sequences in the public databases provides opportunity for studying this large family of structurally complex catalytic RNA by large-scale comparative sequence analysis. In this study, the detailed secondary structures of 211 group I introns in the IE subgroup were manually predicted. The secondary structure-favored alignments showed that IE introns contain 14 conserved stems. The P13 stem formed by long-range base-pairing between P2.1 and P9.1 is conserved among IE introns. Sequence variations in the conserved core divide IE introns into three distinct minor subgroups, namely IE1, IE2 and IE3. Co-variation of the peripheral structural motifs with core sequences supports that the peripheral elements function in assisting the core structure folding. Interestingly, host-specific structural motifs were found in IE2 introns inserted at S516 position. Competitive base-pairing is found to be conserved at the junctions of all long-range paired regions, suggesting a possible mechanism of establishing long-range base-pairing during large RNA folding. These findings extend our knowledge of IE introns, indicating that comparative analysis can be a very good complement for deepening our understanding of RNA structure and function in the genomic era.     

Comparison and Evolutionary Analysis of the Glycosomal Glyceraldehyde-3-Phosphate Dehydrogenase from Different Kinetoplastida

In this work, we present the sequences and a comparison of the glycosomal GAPDHs from a number of Kinetoplastida. The complete gene sequences have been determined for some species (Crithidia fasciculata, Herpetomonas samuelpessoai, Leptomonas seymouri, and Phytomonas sp), whereas for other species (Trypanosoma brucei gambiense, Trypanosoma congolense, Trypanosoma vivax, and Leishmania major), only partial sequences have been obtained by PCR amplification. The structure of all available glycosomal GAPDH genes was analyzed in detail. Considerable variations were observed in both their nucleotide composition and their codon usage. The GC content varies between 64.4% in L. seymouri and 49.5% in the previously sequenced GAPDH gene from Trypanoplasma borreli. A highly biased codon usage was found in C. fasciculata, with only 34 triplets used, whereas in T. borreli 57 codons were employed. No obvious correlation could be observed between the codon usage and either the nucleotide composition or the level of gene expression. The glycosomal GAPDH is a very well-conserved enzyme. The maximal overall difference observed in the amino acid sequences is only 25%. Specific insertions and extensions are retained in all sequences. The residues involved in catalysis, substrate, and inorganic phosphate binding are fully conserved, whereas some variability is observed in the cofactor-binding pocket. The implications of these data for the design of new trypanocidal drugs targeted against GAPDH are discussed. All available gene and amino acid sequences of glycosomal GAPDHs were used for a phylogenetic analysis. The division of the Kinetoplastida into two suborders, Bodonina and Trypanosomatina, was well supported. Within the letter group, the Trypanosoma species appeared to be monophyletic, whereas the other trypanosomatids form a second clade. Received: 23 February 1998/Accepted: 26 March 1998     

Molecular Evolution of P Transposable Elements in the Genus Drosophila. II. The obscura Species Group

A phylogenetic analysis of P transposable elements in the Drosophila obscura species group is described. Multiple P sequences from each of 10 species were obtained using PCR primers that flank a conserved region of exon 2 of the transposase gene. In general, the P element phylogeny is congruent with the species phylogeny, indicating that the dominant mode of transmission has been vertical, from generation to generation. One manifestation of this is the distinction of P elements from the Old World obscura and subobscura subgroups from those of the New World affinis subgroup. However, the overall distribution of elements within the obscura species group is not congruent with the phylogenetic relationships of the species themselves. There are at least four distinct subfamilies of P elements, which differ in sequence from each other by as much as 34%, and some individual species carry sequences belonging to different subfamilies. P sequences from D. bifasciata are particularly interesting. These sequences belong to two subfamilies and both are distinct from all other P elements identified in this survey. Several mechanisms are postulated to be involved in determining phylogenetic relationships among P elements in the obscura group. In addition to vertical transmission, these include retention of ancestral polymorphisms and horizontal transfer by an unknown mating-independent mechanism.     

Sequences and Evolution of Human and Squirrel Monkey Blue Opsin Genes

The sequences of the entire blue opsin gene in the squirrel monkey (Saimiri boliviensis) and the five introns of the human blue opsin gene were obtained. Intron 3 of these genes contains an Alu sequence and intron 4 contains a partial mer13 sequence. A comparison of the squirrel monkey opsin sequence with published mammalian opsin sequences shows that features believed to be functionally critical are all conserved. However, the blue opsin has evolved twice as fast as rhodopsin and is only as conservative as the β globin, which has evolved at the average rate of mammalian proteins. Interestingly, the interhelical loops are, on average, actually more conservative than the transmembrane α helical regions. The introns of the blue opsin gene have evolved at the average rate of introns in primate genes. Received: 5 August 1996 / Accepted: 2 October 1996     

Heat Shock Protein 70 Family: Multiple Sequence Comparisons, Function, and Evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus ``matches' 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins. Received: 29 January 1998 / Accepted: 14 May 1998     

Sequence Peculiarity of Gnetalean Legumin-Like Seed Storage Proteins

The development of seeds as a specialized organ for the nutrition, protection, and dispersal of the next generation was an important step in the evolution of land plants. Seed maturation is accompanied by massive synthesis of storage compounds such as proteins, starch, and lipids. To study the processes of seed storage protein evolution we have partially sequenced storage proteins from maturing seeds of representatives from the gymnosperm genera Gnetum, Ephedra, and Welwitschia—morphologically diverse and unusual taxa that are grouped in most formal systems into the common order Gnetales. Based on partial N-terminal amino acid sequences, oligonucleotide primers were derived and used for PCR amplification and cloning of the corresponding cDNAs. We also describe the structure of the nuclear gene for legumin of Welwitschia mirabilis. This first gnetalean nuclear gene structure contains introns in only two of the four conserved positions previously characterized in other spermatophyte legumin genes. The distinct phylogenetic status of the gnetalean taxa is also reflected in a sequence peculiarity of their legumin genes. A comparative analysis of exon/intron sequences leads to the hypothesis that legumin genes from Gnetales belong to a monophyletic evolutionary branch clearly distinct from that of legumin genes of extant Ginkgoales and Coniferales as well as from all angiosperms. Received: 5 June 1997 / Accepted: 31 March 1998