Distribution of introns in fungal histone genes

Saccharomycotina and Taphrinomycotina lack intron in their histone genes, except for an intron in one of histone H4 genes of Yarrowia lipolytica. On the other hand, Basidiomycota and Perizomycotina have introns in their histone genes. We compared the distributions of 81, 47, 79, and 98 introns in the fungal histone H2A, H2B, H3, and H4 genes, respectively. Based on the multiple alignments of the amino acid sequences of histones, we identified 19, 13, 31, and 22 intron insertion sites in the histone H2A, H2B, H3, and H4 genes, respectively. Surprisingly only one hot spot of introns in the histone H2A gene is shared between Basidiomycota and Perizomycotina, suggesting that most of introns of Basidiomycota and Perizomycotina were acquired independently. Our findings suggest that the common ancestor of Ascomycota and Basidiomycota maybe had a few introns in the histone genes. In the course of fungal evolution, Saccharomycotina and Taphrinomycotina lost the histone introns; Basidiomycota and Perizomycotina acquired other introns independently. In addition, most of the introns have sequence similarity among introns of phylogenetically close species, strongly suggesting that horizontal intron transfer events between phylogenetically distant species have not occurred recently in the fungal histone genes.     

Molecular characterization of Yarrowia lipolytica and Candida zeylanoides isolated from poultry

Yeast isolates from raw and processed poultry products were characterized using PCR amplification of the internally transcribed spacer (ITS) 5.8S ribosomal DNA region (ITS-PCR), restriction analysis of amplified products, randomly amplified polymorphic DNA (RAPD) analysis, and pulsed-field gel electrophoresis (PFGE). ITS-PCR resulted in single fragments of 350 and 650 bp, respectively, from eight strains of Yarrowia lipolytica and seven strains of Candida zeylanoides. Digestion of amplicons with HinfI and HaeIII produced two fragments of 200 and 150 bp from Y. lipolytica and three fragments of 350, 150, and 100 bp from C. zeylanoides, respectively. Although these fragments showed species-specific patterns and confirmed species identification, characterization did not enable intraspecies typing. Contour-clamped heterogeneous electric field PFGE separated chromosomal DNA of Y. lipolytica into three to five bands, most larger than 2 Mbp, whereas six to eight bands in the range of 750 to 2,200 bp were obtained from C. zeylanoides. Karyotypes of both yeasts showed different polymorphic patterns among strains. RAPD analysis, using enterobacterial repetitive intergenic sequences as primers, discriminated between strains within the same species. Cluster analysis of patterns formed groups that correlated with the source of isolation. For ITS-PCR, extraction of DNA by boiling yeast cells was successfully used.     

Promiscuous DNA in the nuclear genomes of hemiascomycetous yeasts

Transfer of fragments of mtDNA to the nuclear genome is a general phenomenon that gives rise to NUMTs (NUclear sequences of MiTochondrial origin). We present here the first comparative analysis of the NUMT content of entirely sequenced species belonging to a monophyletic group, the hemiascomycetous yeasts ( Candida glabrata, Kluyveromyces lactis, Kluyveromyces thermotolerans, Debaryomyces hansenii and Yarrowia lipolytica , along with the updated NUMT content of Saccharomyces cerevisiae ). This study revealed a huge diversity in NUMT number and organization across the six species. Debaryomyces hansenii harbors the highest number of NUMTs (145), half of which are distributed in numerous large mosaics of up to eight NUMTs arising from multiple noncontiguous mtDNA fragments inserted at the same chromosomal locus. Most NUMTs, in all species, are found within intergenic regions including seven NUMTs in pseudogenes. However, five NUMTs overlap a gene, suggesting a positive impact of NUMTs on protein evolution. Contrary to the other species, K. lactis and K. thermotolerans harbor only a few diverged NUMTs, suggesting that mitochondrial transfer to the nuclear genome has decreased or ceased in these phylogenetic branches. The dynamics of NUMT acquisition and loss are illustrated here by their species-specific distribution.     

Palindromic repetitive elements in the mitochondrial genome of Volvox

Group I introns were found in the cob and cox I genes of Volvox carteri. These introns contain tandem arrays of short palindromic sequences that are related to each other. Inspection of other regions in the mtDNA revealed that similar palindromic repetitive sequences are dispersed in the non-protein coding regions of the mitochondrial genome. Analysis of the group I intron in the cob gene of another member of Volvocaceae, Volvox aureus, has shown that its sequence is highly homologous to its counterpart in V. carteri with the exception of a cluster of palindromic sequences not found in V. carteri. This indicates that the palindromic clusters were inserted into the introns after divergence of the two species, presumably due to frequent insertions of the palindromic elements during evolution of the Volvocaceae. Possible involvement of the palindromic repetitive elements in the molecular evolution of functional RNAs is discussed.     

Evolutionary relationships among pathogenic Candida species and relatives.

Small subunit rRNA sequences have been determined for 10 of the most clinically important pathogenic species of the yeast genus Candida (including Torulopsis [Candida] glabrata and Yarrowia [Candida] lipolytica) and for Hansenula polymorpha. Phylogenetic analyses of these sequences and those of Saccharomyces cerevisiae, Kluyveromyces marxianus var. lactis, and Aspergillus fumigatus indicate that Candida albicans, C. tropicalis, C. parapsilosis, and C. viswanathii form a subgroup within the genus. The remaining significant pathogen, T. glabrata, falls into a second, distinct subgroup and is specifically related to S. cerevisiae and more distantly related to C. kefyr (psuedotropicalis) and K. marxianus var. lactis. The 18S rRNA sequence of Y. lipolytica has evolved rapidly in relation to the other Candida sequences examined and appears to be only distantly related to them. As anticipated, species of several other genera appear to bear specific relationships to members of the genus Candida.     

Yarrowia lipolytica mutants devoid of pyruvate carboxylase activity show an unusual growth phenotype

We have cloned and characterized the gene PYC1, encoding the unique pyruvate carboxylase in the dimorphic yeast Yarrowia lipolytica. The protein putatively encoded by the cDNA has a length of 1,192 amino acids and shows around 70% identity with pyruvate carboxylases from other organisms. The corresponding genomic DNA possesses an intron of 269 bp located 133 bp downstream of the starting ATG. In the branch motif of the intron, the sequence CCCTAAC, not previously found at this place in spliceosomal introns of Y. lipolytica, was uncovered. Disruption of the PYC1 gene from Y. lipolytica did not abolish growth in glucose-ammonium medium, as is the case in other eukaryotic microorganisms. This unusual growth phenotype was due to an incomplete glucose repression of the function of the glyoxylate cycle, as shown by the lack of growth in that medium of double pyc1 icl1 mutants lacking both pyruvate carboxylase and isocitrate lyase activity. These mutants grew when glutamate, aspartate, or Casamino Acids were added to the glucose-ammonium medium. The cDNA from the Y. lipolytica PYC1 gene complemented the growth defect of a Saccharomyces cerevisiae pyc1 pyc2 mutant, but introduction of either the S. cerevisiae PYC1 or PYC2 gene into Y. lipolytica did not result in detectable pyruvate carboxylase activity or in growth on glucose-ammonium of a Y. lipolytica pyc1 icl1 double mutant.     

Complete sequence of a bovine type I cytokeratin gene: conserved and variable intron positions in genes of polypeptides of the same cytokeratin subfamily

The complete sequence of a bovine gene encoding an epidermal cytokeratin of mol. wt. 54 500 (No VIb) of the acidic (type I) subfamily is presented, including an extended 5' upstream region. The gene (4377 bp, seven introns) which codes for a representative of the glycine-rich subtype of cytokeratins of this subfamily, is compared with genes coding for: another subtype of type I cytokeratin; a basic (type II) cytokeratin gene; and vimentin, a representative of another intermediate filament (IF) protein class. The positions of the five introns located within the highly homologous alpha-helix-rich rod domain are identical or equivalent, i.e., within the same triplet, in the two cytokeratin I genes. Four of these intron positions are also identical with intron sites in the vimentin gene, and three of these intron positions are identical or similar in the type I and type II cytokeratin subfamilies. On the other hand, the gene organization of both type I cytokeratins differs from that of the type II cytokeratin in the rod region in five intron positions and in the introns located in the carboxy-terminal tail region, with the exception of one position at the rod-tail junction. Remarkably, the two type I cytokeratins also differ from each other in the positions of two introns located at and in the region coding for the hypervariable, carboxy-terminal portion. The introns and the 5' upstream regions of the cytokeratin VIb gene do not display notable sequence homologies with the other IF protein genes, but sequences identical with--or very similar to--certain viral and immunoglobulin enhancers have been identified.(ABSTRACT TRUNCATED AT 250 WORDS)     

Phylogenetics of wigeons and allies (Anatidae: Anas): the importance of sampling multiple loci and multiple individuals

Species-level DNA phylogenies frequently suffer from two shortcomings--gene trees usually are constructed from a single locus, and often species are represented by only one individual. To evaluate the effect of these two shortcomings, we tested phylogenetic hypotheses within the wigeons and allies, a clade of Anas ducks (Anatidae) composed of five species. We sequenced two nuclear introns from the Z-chromosome-linked chromo-helicase binding protein gene (CHD1Zb and CHD1Za) and the mitochondrial DNA (mtDNA) control region for multiple individuals sampled from widespread geographic locations. We compared these phylogenies to previously published phylogenies constructed from morphology and protein coding regions of mtDNA. Relative to other nuclear introns, CHD showed remarkable phylogenetic utility. Of the 26 CHD1Zb alleles identified, only one was shared between two species, and the combined CHD datasets revealed that four of the five species were consistent with monophyly. Several species shared mtDNA haplotypes, which probably was a result of interspecific hybridization. Overall, the nuclear CHD tree and the mtDNA tree were more congruent with coding regions of mtDNA than they were with morphology.     

Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms

Sequences of spacers and group I introns in plant chloroplast genomes have recently been shown to be very effective in phylogenetic reconstruction at higher taxonomic levels and not only for inferring relationships among species. Group II introns, being more frequent in those genomes than group I introns, may be further promising markers. Because group II introns are structurally constrained, we assumed that sequences of a group II intron should be alignable across seed plants. We designed universal amplification primers for the petD intron and sequenced this intron in a representative selection of 47 angiosperms and three gymnosperms. Our sampling of taxa is the most representative of major seed plant lineages to date for group II introns. Through differential analysis of structural partitions, we studied patterns of molecular evolution and their contribution to phylogenetic signal. Nonpairing stretches (loops, bulges, and interhelical nucleotides) were considerably more variable in both substitutions and indels than in helical elements. Differences among the domains are basically a function of their structural composition. After the exclusion of four mutational hotspots accounting for less than 18% of sequence length, which are located in loops of domains I and IV, all sequences could be aligned unambiguously across seed plants. Microstructural changes predominantly occurred in loop regions and are mostly simple sequence repeats. An indel matrix comprising 241 characters revealed microstructural changes to be of lower homoplasy than are substitutions. In showing Amborella first branching and providing support for a magnoliid clade through a synapomorphic indel, the petD data set proved effective in testing between alternative hypotheses on the basal nodes of the angiosperm tree. Within angiosperms, group II introns offer phylogenetic signal that is intermediate in information content between that of spacers and group I introns on the one hand and coding sequences on the other.     

Evolution of rbcL group IA introns and intron open reading frames within the colonial Volvocales (Chlorophyceae)

Mobile group I introns sometimes contain an open reading frame (ORF) possibly encoding a site-specific DNA endonuclease. However, previous phylogenetic studies have not clearly deduced the evolutionary roles of the group I intron ORFs. In this paper, we examined the phylogeny of group IA2 introns inserted in the position identical to that of the chloroplast-encoded rbcL coding region (rbcL-462 introns) and their ORFs from 13 strains of five genera (Volvox, Pleodorina, Volvulina, Astrephomene, and Gonium) of the colonial Volvocales (Chlorophyceae) and a related unicellular green alga, Vitreochlamys. The rbcL-462 introns contained an intact or degenerate ORF of various sizes except for the Gonium multicoccum rbcL-462 intron. Partial amino acid sequences of some rbcL-462 intron ORFs exhibited possible homology to the endo/excinuclease amino acid terminal domain. The distribution of the rbcL-462 introns is sporadic in the phylogenetic trees of the colonial Volvocales based on the five chloroplast exon sequences (6021 bp). Phylogenetic analyses of the conserved intron sequences resolved that the G. multicoccum rbcL-462 intron had a phylogenetic position separate from those of other colonial volvocalean rbcL-462 introns, indicating the recent horizontal transmission of the intron in the G. multicoccum lineage. However, the combined data set from conserved intron sequences and ORFs from most of the rbcL-462 introns resolved robust phylogenetic relationships of the introns that were consistent with those of the host organisms. Therefore, most of the extant rbcL-462 introns may have been vertically inherited from the common ancestor of their host organisms, whereas such introns may have been lost in other lineages during evolution of the colonial Volvocales. In addition, apparently higher synonymous substitutions than nonsynonymous substitutions in the rbcL-462 intron ORFs indicated that the ORFs might evolve under functional constraint, which could result in homing of the rbcL-462 intron in cases of spontaneous intron loss. On the other hand, the presence of intact to largely degenerate ORFs of the rbcL-462 introns within the three isolates of Gonium viridistellatum and the rare occurrence of the ORF-lacking rbcL-462 intron suggested that the ORFs might degenerate to result in the spontaneous intron loss during a very short evolutionary time following the loss of the ORF function. Thus, the sporadic distribution of the rbcL-462 introns within the colonial Volvocales can be largely explained by an equilibrium between maintenance of the introns by the intron ORF and spontaneous loss of introns when the introns do not have a functional ORF.     

Discovery of a group I intron in the SSU rDNA of Ploeotia costata (Euglenozoa)

The gene coding for the small ribosomal subunit RNA of Ploeotia costata contains an actively splicing group I intron (Pco.S516) which is unique among euglenozoans. Secondary structure predictions indicate that paired segments P1-P10 as well as several conserved elements typical of group I introns and of subclass IC1 in particular are present. Phylogenetic analyses of SSU rDNA sequences demonstrate a well-supported placement of Ploeotia costata within the Euglenozoa; whereas, analyses of intron data sets uncover a close phylogenetic relation of Pco.S516 to S-516 introns from Acanthamoeba, Aureoumbra lagunensis (Stramenopila) and red algae of the order Bangiales. Discrepancies between SSU rDNA and intron phylogenies suggest horizontal spread of the group I intron. Monophyly of IC1 516 introns from Ploeotia costata, A. lagunensis and rhodophytes is supported by a unique secondary structure element: helix P5b possesses an insertion of 19 nt length with a highly conserved tetraloop which is supposed to take part in tertiary interactions. Neither functional nor degenerated ORFs coding for homing endonucleases can be identified in Pco.S516. Nevertheless, degenerated ORFs with His-Cys box motifs in closely related intron sequences indicate that homing may have occurred during evolution of the investigated intron group.     

The RNA polymerase III-dependent family of genes in hemiascomycetes: comparative RNomics, decoding strategies, transcription and evolutionary implications

We present the first comprehensive analysis of RNA polymerase III (Pol III) transcribed genes in ten yeast genomes. This set includes all tRNA genes (tDNA) and genes coding for SNR6 (U6), SNR52, SCR1 and RPR1 RNA in the nine hemiascomycetes Saccharomyces cerevisiae, Saccharomyces castellii, Candida glabrata, Kluyveromyces waltii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and the archiascomycete Schizosaccharomyces pombe. We systematically analysed sequence specificities of tRNA genes, polymorphism, variability of introns, gene redundancy and gene clustering. Analysis of decoding strategies showed that yeasts close to S.cerevisiae use bacterial decoding rules to read the Leu CUN and Arg CGN codons, in contrast to all other known Eukaryotes. In D.hansenii and C.albicans, we identified a novel tDNA-Leu (AAG), reading the Leu CUU/CUC/CUA codons with an unusual G at position 32. A systematic 'p-distance tree' using the 60 variable positions of the tRNA molecule revealed that most tDNAs cluster into amino acid-specific sub-trees, suggesting that, within hemiascomycetes, orthologous tDNAs are more closely related than paralogs. We finally determined the bipartite A- and B-box sequences recognized by TFIIIC. These minimal sequences are nearly conserved throughout hemiascomycetes and were satisfactorily retrieved at appropriate locations in other Pol III genes.     

A novel group I intron in Candida dubliniensis is homologous to a Candida albicans intron

In the present study, we determined the sequence of group I self-splicing introns found in the large ribosomal RNA subunit of Candida albicans, Candida stellatoidea and the recently-described species Candida dubliniensis. It was found that both the intron and ribosomal RNA nucleotide sequences are almost perfectly identical between different C. albicans strains as well as between C. albicans and C. stellatoidea strains. Comparisons of ribosomal RNA sequences suggest that local isolates of atypical C. albicans from individuals infected with human immunodeficiency virus can be assigned to the C. dubliniensis species. C. dubliniensis strains also harbor a group I intron in their ribosomal RNA, as observed in about 40% of C. albicans strains and all C. stellatoidea strains. This novel C. dubliniensis group I intron is identical to the C. albicans and C. stellatoidea intron, except for two widely divergent stem-loop regions. Despite these differences, the C. dubliniensis intron possesses self-splicing ability in an in vitro assay. Taken together, these data support the idea that C. albicans and C. stellatoidea should be joined together as variants of the same species while C. dubliniensis is a distinct but closely related microorganism. To our knowledge, the C. albicans and C. dubliniensis introns are the first example of a pair of homologous group I introns differing only by the presence of apparently facultative sequences in some stem-loops suspected to be involved in stabilization of tertiary structure.     

Using exon and intron sequences of the gene Mp20 to resolve basal relationships in Cicindela (Coleoptera:Cicindelidae)

The genus Cicindela (Coleoptera: Cicindelidae) is a species-rich cosmopolitan group of tiger beetles useful for comparing clade diversification worldwide. Knowledge about relationships of major groups is important for this analysis but basal nodes in Cicindela have been difficult to resolve with standard mtDNA markers. Here we developed the Mp20 gene, a single-copy nuclear marker coding for a muscle-associated protein in insects, for phylogenetic analysis of basal groups of Cicindela. Nearly full-length sequences were obtained for 51 cicindelids, including major taxonomic groups from all continents. Sequences of Mp20 were between 1.2 and 1.7 kb and spanning three introns. Phylogenetic signal of exon and intron sequences was compared with that from four gene regions of mtDNA (COI, COIII, Cytb, 16S rRNA; 2.4 kb total). Because introns differed in length, sequence alignment was conducted using various procedures of phenetic and parsimony-based character coding of indels to assess their phylogenetic information content, but major nodes were recovered consistently. Mp20 sequences contributed two thirds of the total support of the combined analysis, with most signal from the introns. We found major clades of Cicindela to be geographically largely coincident with continental regions, confined to Australasia, the Holarctic, the Indian subcontinent, Africa, and South and Central America. Clock estimates using various maximum-likelihood (ML) branch length calculations resulted in roughly similar divergence times whether Mp20 exon, introns, or mtDNA were used, and they were not greatly affected by different procedures for coding and optimizing indel characters. Based on existing clock calibrations in Cicindela, basal splits of continental lineages occurred in the mid-Miocene, placing the radiation of basal groups of Cicindela to a period when their open-vegetation habitats expanded globally.     

The chloroplast chlL gene of the green alga Chlorella vulgaris C-27 contains a self-splicing group I intron

The chlL gene product is involved in the light-independent synthesis of chlorophyll in photosynthetic bacteria, green algae and non-flowering plants. The chloroplast genome of Chlorella vulgaris strain C-27 contains the first example of a split chlL gene, which is interrupted by a 951?bp group I intron in the coding region. In vitro synthesized pre-mRNA containing the entire intron and parts of the flanking exon sequences is able to efficiently self-splice in vitro in the presence of a divalent and a monovalent cation and GTP, to yield the ligated exons and other splicing intermediates characteristic of self-splicing group I introns. The 5′ and 3′ splice sites were confirmed by cDNA sequencing and the products of the splicing reaction were characterized by primer extension analysis. The absence of a significant ORF in the long P9 region (522?nt), separating the catalytic core from the 3′ splice site, makes this intron different from the other known examples of group I introns. Guanosine-mediated attack at the 3′ splice site and the presence of G-exchange reaction sites internal to the intron are some other properties demonstrated for the first time by an intron of a protein-coding plastid gene.     

Nucleolar introns from Physarum flavicomum contain insertion elements that may explain how mobile group I introns gained their open reading frames.

Comparison of two group I intron sequences in the nucleolar genome of the myxomycete Physarum flavicomum to their homologs in the closely related Physarum polycephalum revealed insertion-like elements. One of the insertion-like elements consists of two repetitive sequence motifs of 11 and 101 bp in five and three copies, respectively. The smaller motif, which flanks the larger, resembles a target duplication and indicates a relationship to transposons or retroelements. The insertion-like elements are found in the peripheral loops of the RNA structure; the positions occupied by the ORFs of mobile nucleolar group I introns. The P. flavicomum introns are 1184 and 637 bp in size, located in the large subunit ribosomal RNA gene, and can be folded into group I intron structures at the RNA level. However, the intron 2s from both P. flavicomum and P. polycephalum contain an unusual core region that lacks the P8 segment. None of the introns are able to self-splice in vitro. Southern analysis of different isolates indicates that the introns are not optional in myxomycetes.