Adaptation for Horizontal Transfer in a Homing Endonuclease

Vassiliki Koufopanou, Matthew R. Goddard and Austin Burt

*Department of Biology
{dagger}NERC Centre for Population Biology, Imperial College, Silwood Park, Ascot, Berkshire, U.K


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 
Selfish genes of no function other than self-propagation are susceptible to degeneration if they become fixed in a population, and regular transfer to new species may be the only means for their long-term persistence. To test this idea we surveyed 24 species of yeast for VDE, a nuclear, intein-associated homing endonuclease gene (HEG) originally discovered in Saccharomyces cerevisiae. Phylogenetic analyses show that horizontal transmission has been a regular occurrence in its evolutionary history. Moreover, VDE appears to be specifically adapted for horizontal transmission. Its 31-bp recognition sequence is an unusually well-conserved region in an unusually well-conserved gene. In addition, the nine nucleotide sites most critical for homing are also unusually well conserved. Such adaptation for horizontal transmission presumably arose as a consequence of selection, both among HEGs at different locations in the genome and among variants at the same location. The frequency of horizontal transmission must therefore be a key feature constraining the distribution and abundance of these genes.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 
Homing endonuclease genes (HEG) are nonessential genes with no known function, at least in eukaryotes (Mueller et al. 1993Citation ; Belfort and Roberts 1997Citation ; Jurica and Stoddard 1999Citation ; Gimble 2000Citation ). Rather, they are thought to be selfish or parasitic genes that spread and persist in populations because their catalytic activity results in super-Mendelian inheritance. Frequent horizontal transmission has previously been demonstrated for mitochondrial HEGs associated with self-splicing introns in both plants and fungi (Cho et al. 1998Citation ; Goddard and Burt 1999Citation ). HEGs also exist in the nucleus (Dalgaard et al. 1997aCitation ; Perler, Olsen, and Adam 1997Citation ), where horizontal transmission might be expected to be more difficult. VDE (also known as PI-SceI) is an apparently selfish element inside the VMA1 gene of Saccharomyces cerevisiae (Gimble and Thorner 1992Citation ). It encodes a compound protein consisting of a self-splicing intein domain (analogous to a self-splicing intron but operating at the protein level) and an endonuclease domain with two LAGLIDADG motifs. These two domains are functionally independent and have had separate evolutionary origins (Dalgaard et al. 1997aCitation ; Derbyshire et al. 1997Citation ). VMA1 genes that contain VDE are translated into a polypeptide, and then VDE is spliced out in an autocatalytic reaction which ligates the two flanking regions, leaving a functional VMA1 protein. In individuals heterozygous for the presence of VDE, the excised VDE protein recognizes a 31-bp sequence of the VDE- allele and cuts it; the VDE+ allele is immune to cutting as the VDE gene interrupts the 31-bp sequence. The presence of the cut allele induces the cell's recombinational repair system, which uses the VDE-containing homologue as a template for repair, after which both alleles are VDE+. This "homing" mechanism results in VDE being inherited in 90% of meiotic products, rather than the Mendelian 50% (Gimble and Thorner 1992Citation ), which should allow it to increase in frequency and go to fixation within a population. However, once fixed, there would be little or no selection maintaining endonuclease function, as there would be no chromosomes left to cut (except for rare deletion mutations), and one might expect the endonuclease domain to degenerate over evolutionary time because of mutation pressure, aided perhaps by natural selection if there is any cost to the host cell of producing a functional endonuclease (Goddard and Burt 1999Citation ; Gimble 2000Citation ). The splicing domain, however, should remain functional, regardless of saturation, as there will always be strong selection for a functional VMA1 gene.

If fixation is followed by degeneration, then the only way for VDE to persist as a HEG over long evolutionary time spans would be to occasionally move to a new location, on average at least once before degeneration. In principle, movement might be to a new place in the same genome (transposition) or to a different species (horizontal transmission). For mitochondrial HEGs associated with self-splicing introns, horizontal transmission is much more common than transposition (Cho et al. 1998Citation ; Goddard and Burt 1999Citation ). In no case is the mechanism of horizontal transmission known, and unlike transposable elements, there is no extrachromosomal phase in the normal propagation of a HEG which might facilitate horizontal transmission. For yeasts, perhaps the simplest scenario is that there is occasionally some leakage of mtDNA during the preliminary phase of interspecific matings that are subsequently aborted. Such direct transfer might be more difficult for a nuclear HEG, as it would seem to require full nuclear fusion, and hybrids of even close relatives are usually sterile (Naumov 1987Citation ; Marinoni et al. 1999Citation ). It is for these reasons that the study of a nuclear HEG is of interest. Furthermore, note that, whatever the mechanism of crossing the species barrier, a key prerequisite for the successful transfer and spread of a HEG is that the recognition sequence exists in the recipient species. If the ability to undergo horizontal transmission is an important factor in the persistence of a HEG, then one expects there to have been selection against HEGs with poorly conserved recognition sequences. VDE, which has persisted, should therefore have an unusually well-conserved recognition sequence.

To test whether VDE shows horizontal transmission, we surveyed 22 species of saccharomycete yeasts from four closely related genera (Saccharomyces, Torulaspora, Zygosaccharomyces, and Kluyveromyces) for the presence or absence of VDE in the VMA1 gene or elsewhere in the genome and combined our results with published data from two more distantly related yeasts, Candida tropicalis (Gu et al. 1993Citation ) and C. albicans (Stanford DNA Sequencing and Technology Centre; http://www-sequence.Stanford.edu/group/candida). We then analyzed the distribution of VDE in the 24 species and compared the phylogenetic relationships among all VDEs recovered with those of the host species. To examine whether VDE might be specifically adapted for horizontal transmission, we tested whether the VDE recognition site is an unusually well-conserved region of the genome.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 
Yeast Strains and Molecular Methods
All strains were obtained from the Centraalbureau Voor Schimmelcultures (CBS; Delft, The Netherlands), except S. paradoxus, S. cariocanus, and S. kudriavzevii, which were obtained from Dr. Edward Louis, Leicester. Strains were grown overnight and DNA extracted (Strathern and Higgins 1991Citation ). Presence or absence of VDE in the VMA1 gene was assayed by PCR using degenerate primers flanking the VDE insertion site, 83-bp upstream and 291-bp downstream, based on VMA1 sequences from S. cerevisiae, C. tropicalis, Schizosaccharomyces pombe, and Neurospora crassa (primer sequences, 5' to 3': VR1: GGI GCI TTY GGI TGY GGI AA; VR2: GAA ICC TTG RTC IGC IAG CAT, where I is inosine, Y is C + T and R is A + G). All PCR products, both with and without VDE, were sequenced directly using an ABI 373 automated sequencer (EMBL accession numbers AJ307891–902 and AJ312425–33). To screen for VDE elsewhere in the genome, total genomic DNA from the 22 species was digested with HaeIII, electrophoresed, and transferred to a nylon membrane. It was then probed with PCR products from VDE sequences obtained using primers internal and near the two ends of VDE. For control, the membrane was also probed separately with PCR products from empty (i.e., VDE-) VMA1 alleles, amplified with primers VR1 and VR2.

Phylogenetic Analyses
DNA sequences were translated using the ABI Sequence Navigator program and amino acid sequences aligned with Clustal W (Thompson, Higgins, and Gibson 1994Citation ). To reconstruct the host phylogeny, nuclear 18S, 26S, and ITS1-5.8S-ITS2 data were obtained from GenBank (Cai, Roberts, and Collins 1996Citation ; James et al. 1997Citation ; Kurtzman and Robnet 1998Citation ; Goddard and Burt 1999Citation ). Sequence alignments can be obtained from the authors. Phylogenetic analyses were performed with PAUP* (Version 4.0d64) (Swofford 1999Citation ), MacClade (Version 3.07) (Maddison and Maddison 1992Citation ), MultiLocus (Version 1.0b; http://www.bio.ic.ac.uk/evolve/software/multilocus), and Mathematica (Version 3.0.1) (Wolfram 1996Citation ).


    Results and Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 
Detection and Isolation
All species gave a single, unambiguous PCR product, the size of which indicated whether it was VDE+ or VDE-, except Z. rouxii, which gave both VDE+ and VDE- amplicons. These latter bands were reconfirmed by multiple isolations and PCR-amplifications of VMA1 from genomic DNA initially derived from single cell cultures. Both amplicons were sequenced, and the VMA1 fragments differed at 12% of nucleotide sites (all changes but one were silent). This is substantially more than would be expected for alleles at a single locus and is greater than the divergence between some pairs of species (e.g., 3.6% nucleotide divergence for S. cerevisiae-S. cariocanus VMA1 fragments), suggesting that there has been a gene duplication event in the lineage leading to Z. rouxii. Note that the high divergence of the empty VMA1 allele might prevent it from being repaired by the other copy of the gene were it cut by the endonuclease. Only the VDE+ gene is included in our analyses. Of the 24 species, 14 had a VDE insertion in the VMA1 gene (fig. 1 ). VDE genes varied in length from 1,230 to 1,551 bp, and for the alignable regions, average amino acid divergence was 40%, ranging 7%–56%. As expected, and unlike the {omega} endonuclease (Goddard and Burt 1999Citation ), none of the 14 sequences contained mutations disrupting the reading frame, as these would interrupt translation of the VMA1 gene. Also, there is no correlation between the presence of VDE and {omega} (G = 0.01, n = 20 species, P > 0.5) and thus no evidence that some species are particularly susceptible or resistant to HEGs in general.



View larger version (32K):
[in this window]
[in a new window]
 
Fig. 1.—Distribution of VDE in the Saccharomycetaceae: shaded boxes indicate taxa with a VDE insertion in the VMA1 gene. The host phylogeny is a maximum parsimony tree (heuristic search, with gaps scored as 5th bases) of 18S, 26S, and ITS1-5.8S-ITS2 rDNA data (3,392 sites in total, of which 2,731 were aligned unambiguously; 313 informative sites). Values along branches are bootstrap values from 100 replicates, lack of number indicating <50% bootstrap support. Relationships within the Saccharomyces sensu stricto group are based on a separate alignment of the ITS region to improve resolution (marked by a break in the tree; note that only a very small part of the ITS region could be aligned unambiguously and thus be included in the global analysis of all 24 species). Bootstrap support for the branch leading to the S. sensu stricto clade was derived from the global analysis; values for branches within that clade are from Goddard and Burt (1999)Citation . Numbers besides species names give the CBS strain identifications, "T" indicating a "Type" strain for that species. Below is a cyclical model for the evolution of VDE within host lineages, with maximum average waiting times

 
To test whether VDE occurs at other locations in the genome, DNA from the 22 species was probed with 7 of the 14 VDE sequences chosen to span the full range of VDE diversity, in two batches. The first batch included VDE probes from S. castellii, S. unisporus, K. lactis, and K. polysporus and gave hybridization signals from these four species plus S. exiguus and S. dairenensis. The second batch included probes from S. cerevisiae, Z. bailii, and T. pretoriensis and gave signals only from these three species. As control, the membrane was also probed separately with the VDE- sequences from S. kudriavzevii and K. thermotolerans, and all species gave a strong hybridization signal. The lack of crosshybridization for VDE is consistent with its high divergence among the species assayed. With the exception of S. exiguus, the position, number, and size of all VDE hybridization bands were as expected from the VDE- hybridizations and the presence of HaeIII restriction sites, i.e., either they were identical to the VDE- hybridization bands or their size was as expected by the presence of additional HaeIII sites inside the VDE sequences. These results indicate that VDE only exists inside the VMA1 gene, if at all. For S. exiguus, one extra band was obtained for both VMA1 and VDE hybridizations, though neither fragment included a HaeIII restriction site, perhaps indicating that the strain is heterozygous for a site beyond the fragments considered here or there has been a gene duplication not detected by PCR. Thus, VDE does not appear to exist at any location in the genome of these yeast species other than the VMA1 gene.

Horizontal Transmission
The distribution of VDE in the Saccharomycetaceae appears to be random, without significant clustering along specific clades, indicating multiple independent gains or losses (or both) (fig. 1 ). Significant clustering of states on the tree would have suggested that gains and losses are relatively infrequent events. The most parsimonious character reconstruction for the observed distribution of VDE on the host tree infers eight steps and is no shorter than the length obtained when presence/absence is randomized on the tree (P = 0.66, from 100 randomizations). We also failed to find significant clustering if gains or losses are modeled as irreversible (P = 0.94 and 0.59, from 100 randomizations, respectively). To estimate the rates of VDE gain and loss, we used a maximum likelihood approach, combining information from the observed distribution of VDE and the maximum parsimony host tree (fig. 1 ; Pagel 1994Citation ). Branch lengths were estimated by maximum likelihood while enforcing a molecular clock (Swofford 1999Citation ) and were calibrated following Berbee and Taylor (1993)Citation . The fit of a two-parameter model, in which rates of gain and loss are allowed to differ, is not significantly better than the one-parameter model with equal rates ({Delta}lnL = 0.33, df = 1, NS [not significant]) but is significantly better than unidirectional models in which only gains or only losses are allowed ({Delta}lnL = 19.46 and 16.62, respectively, df = 1, P < 0.001). Therefore, to estimate the rates of gain and loss we fitted a one-parameter model with equal rates. The estimate maximizing the likelihood of the data is {infty} (infinity), with a lower bound of 0.17 events per million years (from the 2-unit support limit, i.e., the value at which the lnL is 2 units less than the maximum, which is equivalent to the 95% confidence limit; Edwards 1992Citation ; note that this does not take into account error in the phylogeny; fig. 1 ). This is equivalent to a maximum average waiting time of 6 Myr for each transition. As the total time represented in the entire phylogeny is estimated to be 766 Myr, this implies at least 128 gains and losses in the ancestry of our 24 species. Note this last estimate is independent of the temporal calibration of the host tree. It indicates that VDE has been an extremely dynamic element, continuously moving in and out of species throughout its evolutionary history.

A history of horizontal transmission should mean that the phylogenies of host and VDE are significantly different. Before testing this prediction, we first determined whether the self-splicing and endonuclease domains of VDE themselves have different phylogenetic histories (Dalgaard et al. 1997aCitation ). Separate phylogenies for the two domains have no well-supported branches (i.e., with bootstrap values >70%) that are incompatible between trees (not shown), and the data sets do not differ by the partition homogeneity test (P = 0.11; Farris et al. 1995Citation ). However, the host and VDE phylogenies are strikingly different (fig. 2 ). First, each tree has well-supported branches that are incompatible with the other tree. Second, each of the two data sets fits its own tree significantly better than the alternative tree (winning sites test; 83 of 95 sites, P < 0.0001, and 63 of 71 sites, P < 0.0001, for host and VDE data sets, respectively; Swofford 1999Citation ). Finally, the two data sets are significantly different by the partition-homogeneity test (P << 0.001; fig. 2 ; Farris et al. 1995Citation ; Swofford 1999Citation ), despite having similar amounts of homoplasy (consistency index [CI] = Lmin/Lshortest: 0.73 = 475/649 and 0.78 = 472/605, for the host and VDE trees, respectively, with character changes distributed over 153 and 126 informative sites). These incompatibilities provide compelling evidence for horizontal transmission of VDE across yeast species. An absolute lower bound on the number of horizontal transmission events can be obtained by calculating how many branches of one tree have to be cut off and reconnected in a new position to transform it into the other tree (Hein 1993Citation ). Starting with the better resolved VDE tree, it takes four such rearrangements to get the host tree or three rearrangements to get a tree that is not significantly worse than the best host tree (table 1 ). However, such an approach is likely to substantially underestimate the actual number of rearrangements (Goddard and Burt 1999Citation ).



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 2.—Comparison of host and VDE gene genealogies for the 14 species containing VDE. The host genealogy is a maximum parsimony tree, as in figure 1 (153 informative sites). The VDE genealogy is based on an amino acid alignment (maximum parsimony tree, heuristic search; 126 informative sites). Resolution within groups, indicated by breaks on the tree, is based on separate alignments, to maximize the number of informative sites, with branches to the broken clades derived from the global alignment and analysis and those within the clades from local analyses. Numbers along branches are bootstrap support values from 100 replicates, lack of number indicating <50% bootstrap support. The VDE tree was rooted with the Chlamydomonas eugametos clpP gene intein (Huang et al. 1994Citation ; Dalgaard et al. 199Citation 7a). The histogram below shows the distribution of summed tree lengths for 1,000 random partitions of the data, with the arrow indicating the observed partition (partition homogeneity test, P << 0.001; Swofford 1999Citation )

 

View this table:
[in this window]
[in a new window]
 
Table 1 Heuristic Search for the Minimum Number of Horizontal Transfers (clade relocations) Required to Derive the Host Tree from the VDE Genealogy

 
To test whether there is still an association remaining between the host and VDE phylogenies, despite the horizontal transmission events, we compared the length of the shortest tree fitted to the combined host-VDE data set to the length of trees fitted to data sets in which the host-VDE association was randomized (Koufopanou, Burt, and Taylor 1997Citation ). The observed association gave significantly shorter tree lengths than the randomized associations (1,208 vs. 1,220–1,290 for the observed and 100 random associations, respectively; P < 0.01). One possible explanation for such an association is that transfers tend to be between close relatives. Consistent with this possibility, deep nodes in the host tree (leading to C. tropicalis and K. lactis) correspond to deep nodes in the VDE tree (fig. 2 ). Moreover, TBLASTN searches of GenBank with our VDE sequences failed to retrieve any close homologue except the yeast mating-type switching (HO) gene, implying that transfers among yeasts are more common than transfers from the more distantly related taxa known to have intein-associated HEGs (bacteriophages, eubacteria, archaea, and algal chloroplasts; Dalgaard et al. 1997aCitation ; Perler, Olsen, and Adam 1997Citation ).

Adaptation for Horizontal Transfer
Three lines of evidence indicate that the recognition sequence of VDE is particularly well conserved. First, the host gene, VMA1, is unusually well conserved among yeasts. Complete or virtually complete genomes are now available for two yeasts, S. cerevisiae and C. albicans, and we used these to compute an average amino acid divergence of genes: 20 genes of known function were chosen randomly from the S. cerevisiae genome directory (Goffeau et al. 1997Citation ) and used as queries to search the C. albicans genome data base (Stanford DNA sequencing and technology centre; http://www-sequence.Stanford.edu/group/candida). VMA1 is 12% divergent, significantly better conserved than random genes in the two species (one-tailed t-test on arcsine-transformed values, t = 2.76, P < 0.01; fig. 3a ). Second, an alignment of VMA1 genes from six ascomycete fungi shows that the 31-bp recognition sequence of VDE is significantly better conserved than random 31-bp regions within the gene (P = 0.03; fig. 3b ). This last result could be due in part to the homogenization of the recognition sequence, through coconversion of flanking regions during homing (Dujon 1989Citation ; Cho et al. 1998Citation ), though the high divergence of VDE among species suggests that this cannot be a strong effect. Finally, within VDEs 31-bp recognition sequence, there are nine sites which have been identified by site-directed mutagenesis as being critical for effective homing (Gimble and Wang 1996Citation ). Our VMA1 sequences include the recognition site for all 24 species, and these nine sites are perfectly conserved in all 24 species and are significantly better conserved than the other 22 sites (G = 9.83, P < 0.005, fig. 3c ). These results indicate that VDE has evolved to use as its homing site a region that is very highly conserved and thus likely to be present in a large variety of species.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 3.—Conservation of the VDE homing site. a, Frequency distribution of S. cerevisiae-C. albicans amino acid divergence of 20 random S. cerevisiae genes of known function. Arrow indicates divergence of VMA1 (12%; t-test on arcsine-transformed values, t = 2.76, P < 0.01 [one-tailed]). b, Sliding window analysis of nucleotide diversity (probability that two sequences chosen at random are different at a particular site) in VMA1. Full length VMA1 sequences were available for six ascomycete fungi (S. cerevisiae, C. tropicalis, C. albicans, Ashbya gossypii, Neurospora crassa, and Schizosaccharomyces pombe); window size was 31 bp. The position of the VDE recognition sequence is also shown. c, The VDE recognition sequence. Stars mark nucleotide sites critical for effective homing (Gimble and Wang 1996Citation ) and those variable among the 24 species of yeast. Arrow indicates VDE insertion site (Gimble 2001)

 

    Conclusions
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 
Simple evolutionary reasoning suggests that HEGs may be unable to persist over long time spans without regular movement, either by transposition to a new location in the genome or by horizontal transmission to a new species. Similar logic has also been applied to other selfish genetic elements (Hurst, Hurst, and Majerus 1992Citation ; Lohe et al. 1995Citation ; Hurst and McVean 1996Citation ). Motivated by this idea, we surveyed 24 species of yeasts for VDE and found it to be widespread, present in 14 species. Various phylogenetic analyses give strong evidence for regular horizontal transmission, with some suggestion that it is more common between more closely related species; transposition to new sites was not detected. A specific reliance on horizontal transmission might also plausibly explain the absence of HEGs in animals (except in the mitochondrion of the sea anemone Metridium senile; Beagley, Okada, and Wolstenholme 1996Citation ), where access to the germ line is more difficult.

If horizontal transmission is critical in determining the distribution of HEGs, then within taxa there ought to be selection among HEGs for those that are better able to transfer successfully, which in turn would produce adaptation for horizontal transmission. For transfer to be successful (i.e., the transferred gene spreads through the recipient population), the recipient population must contain the HEG recognition site. Therefore, HEGs with highly conserved recognition sequences should transfer successfully more frequently and become more widespread. This line of reasoning suggests that extant HEGs should target unusually well-conserved regions of the genome. The fact that VDE is in a gene instead of free standing in noncoding DNA supports this idea, but our analyses go further, indicating that VDE is in an unusually well-conserved region of an unusually well-conserved gene. This positioning presumably reflects selection among HEGs targeting different regions of the genome. Such selection probably also explains why HEGs are typically found in functionally important parts of genes (Garret et al. 1991Citation ; Dalgaard et al. 1997bCitation ; Edgell, Belfort, and Shub 2000Citation ). Moreover, within the recognition site, VDE targets particularly well-conserved base pairs, and presumably this adaptation arose and is maintained by selection among variant VDEs at the same site, favoring those with increased binding affinity for highly-conserved nonsynonymous sites. These features of VDE biology we take to be adaptations for successful horizontal transmission and evidence that horizontal transmission is a key parameter in the distribution of these selfish genes. There are no adaptations for horizontal transmission, of which we are aware, other than the operon structure of prokaryotic genomes (Lawrence and Roth 1996Citation ).

The mechanism of horizontal transfer is unclear as infectious viruses and plasmids appear to be absent from yeasts. Possibilities include interspecific hybridization (Marinoni et al. 1999Citation ), vectoring by predacious yeasts (Lachance and Pang 1997Citation ), and uptake of naked DNA from the environment (Nevoigt, Fassbender, and Stahl 2000)Citation . Once inside the cell, fragments of DNA are readily integrated into the yeast genome by homologous recombination. Note that this last step could also select for HEGs in well-conserved regions of the genome. Also unclear is the mechanism of intein loss, as imprecise deletions would disrupt the host gene and be highly deleterious.

More generally, the results also show that horizontal transmission in yeasts is not restricted to mitochondrial DNA but that nuclear DNA is also transferred intact from one species to another. HEGs have no extrachromosomal phase in their life cycle (unlike, say, transposable elements), and therefore normal host-benefiting genes probably also get transferred occasionally. The extent to which these show evidence of past horizontal transmission should then be determined by selection, rather than opportunity. Some degree of horizontal transmission of nuclear genes has been shown for filamentous fungi (Rosewich and Kistler 2000)Citation . For yeasts we are not yet aware of any compelling examples, but interestingly some ecologically contingent genes occur in clusters, for example, those involved in biotin synthesis (Phalip et al. 1999Citation ), galactose metabolism (Greger and Proudfoot 1998Citation ), and arsenic tolerance (Bobrowicz et al. 1997Citation ). By analogy with prokaryotes, these would be excellent candidates for horizontal transmission.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 
We thank C. Godfray for useful comments on the manuscript. This work was supported by NERC grant GR3/10626.


    Footnotes
 
Richard Thomas, Reviewing Editor

Abbreviations: VDE, VMA1-derived endonuclease; HEG, homing endonuclease gene; VMA1, vacuolar membrane H+-ATPase. Back

Keywords: horizontal (lateral) transfer (transmission) homing endonuclease adaptation yeast selfish gene VDE distribution Back

Address for correspondence and reprints: Vassiliki Koufopanou, Department of Biology, Imperial College, Silwood Park, Ascot, Berkshire SL5 7PY, UK. v.koufopanou{at}ic.ac.uk . Back


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 References
 

    Beagley C. T., N. A. Okada, D. R. Wolstenholme, 1996 Two mitochondrial group I introns in a metazoan, the sea anemone Metridium senile: one intron contains genes for subunits 1 and 3 of NADH dehydrogenase Proc. Natl. Acad. Sci. USA 93:5619-5623[Abstract/Free Full Text]

    Belfort M., R. J. Roberts, 1997 Homing endonucleases: keeping the house in order Nucleic Acids Res 25:3379-3388[Abstract/Free Full Text]

    Berbee M. L., J. W. Taylor, 1993 Dating the evolutionary radiations of the true fungi Can. J. Bot 71:1114-1127[ISI]

    Bobrowicz P., R. Wysocki, G. Owsianik, A. Goffeau, S. Ulaszewski, 1997 Isolation of three contiguous genes, ACR1,ACR2,ACR3, involved in resistance to arsenic compounds in the yeast Saccharomyces cerevisiae Yeast 13:819-828[ISI][Medline]

    Cai J., I. N. Roberts, M. D. Collins, 1996 Phylogenetic relationships among members of the Ascomycetous yeast genera Brettanomyces, Debaromyces, Dekkera and Kluyveromyces deduced by small-subunit rRNA gene sequences Int. J. Syst. Bacteriol 46:542-549[Abstract]

    Cho Y., Y.-L. Qiu, P. Kuhlman, J. D. Palmer, 1998 Explosive invasion of plant mitochondria by a group I intron Proc. Natl. Acad. Sci. USA 95:14244-14249[Abstract/Free Full Text]

    Dalgaard J. Z., A. J. Klar, M. J. Moser, W. R. Holley, A. Chatterjee, I. S. Mian, 1997a. Statistical modeling and analysis of the LAGLIDADG family of site-specific endonucleases and identification of an intein that encodes a site-specific endonuclease of the HNH family Nucleic Acids Res 25:4626-4638[Abstract/Free Full Text]

    Dalgaard J. Z., M. J. Moser, R. Hughey, I. S. Mian, 1997b. Statistical modeling, phylogenetic analysis and structure prediction of a protein splicing domain common to inteins and hedgehog proteins J. Comput. Biol 4:193-214[ISI][Medline]

    Derbyshire V., D. W. Wood, W. Wu, J. T. Dansereau, J. Z. Dalgaard, M. Belfort, 1997 Genetic definition of a protein-splicing domain: functional mini-inteins support structure predictions and a model for intein evolution Proc. Natl. Acad. Sci. USA 94:11466-11471[Abstract/Free Full Text]

    Dujon B., 1989 Group I introns as mobile genetic elements: facts and mechanistic speculations—a review Gene 82:91-114[ISI][Medline]

    Edgell D. R., M. Belfort, D. A. Shub, 2000 Barriers to intron promiscuity in bacteria J. Bacteriol 182:5281-5289[Free Full Text]

    Edwards A. W. F., 1992 Likelihood Expanded edition. The Johns Hopkins University Press, Baltimore

    Farris J. S., M. Kallersjo, A. G. Kluge, C. Bult, 1995 Testing significance of incongruence Cladistics 10:315-319[ISI]

    Garret R. A., J. Dalgaard, N. Larsen, J. Kjems, A. S. Mankins, 1991 Archaeal rRNA operons Trends Biochem. Sci 16:22-26[ISI][Medline]

    Gimble F. S., 2000 Invasion of a multitude of genetic niches by mobile endonuclease genes FEMS Microbiol. Lett 185:99-107[ISI][Medline]

    Gimble F. S., 2001 Degeneration of a homing endonuclease and its target sequence in a wild yeast strain Nucleic Acids Res 29:4215-4223[Abstract/Free Full Text]

    Gimble F. S., J. Thorner, 1992 Homing of a DNA endonuclease gene by meiotic gene conversion in Saccharomyces cerevisiae Nature 357:301-305[ISI][Medline]

    Gimble F. S., J. Wang, 1996 Substrate recognition and induced DNA distortion by the PI-SceI endonuclease, an enzyme generated by protein splicing J. Mol. Biol 263:163-180[ISI][Medline]

    Goddard M. R., A. Burt, 1999 Recurrent invasion and extinction of a selfish gene Proc. Natl. Acad. Sci. USA 96:1380-1385

    Goffeau A., A. Jones, B. Jones, C. Jones, D. Jones, E. Jones, F. Jones, H. Jones, 1997 The Yeast genome directory Nature 387: (Suppl.) 1-105

    Greger I. H., N. J. Proudfoot, 1998 Poly(A) signals control both transcriptional termination and initiation between the tandem GAL10 and GAL7 genes of Saccharomyces cerevisiae EMBO J 17:4771-4779[Abstract/Free Full Text]

    Gu H. H., J. Xu, M. Gallagher, G. E. Dean, 1993 Peptide splicing in the vacuolar ATPase subunit a from Candida tropicalis J. Biol. Chem 268:7372-7381[Abstract/Free Full Text]

    Hein J., 1993 A heuristic method to reconstruct the history of sequences subject to recombination J. Mol. Evol 36:396-405[ISI]

    Huang C., S. Wang, L. Chen, C. Lemieux, C. Otis, M. Turmel, X. Q. Liu, 1994 The Chlamydomonas chloroplast clpP gene contains translated large insertion sequences and is essential for cell growth Mol. Gen. Genet 244:151-159[ISI][Medline]

    Hurst G. D. D., L. D. Hurst, M. E. N. Majerus, 1992 Selfish genes move sideways Nature 356:659-660[ISI][Medline]

    Hurst L. D., G. T. McVean, 1996 Clade selection, reversible evolution and the persistence of selfish elements: the evolutionary dynamics of cytoplasmic incompatibility Proc. R. Soc. Lond. B 263:97-104[ISI]

    James S. A., J. Cai, I. N. Roberts, M. D. Collins, 1997 A phylogenetic analysis of the genus Saccharomyces based on 18S rRNA gene sequences: description of Saccharomyces kunashirensis sp. nov. and Saccharomyces martiniae sp. nov. Int J. Syst. Bacteriol 47:453-460[Abstract/Free Full Text]

    Jurica M. S., B. L. Stoddard, 1999 Homing endonucleases: structure, function and evolution Cell. Mol. Life Sci 55:1304-1326[ISI][Medline]

    Kishino H., M. Hasegawa, 1989 Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA-sequence data and the branching order in Hominoidea J. Mol. Evol 29:170-179[ISI][Medline]

    Koufopanou V., A. Burt, J. W. Taylor, 1997 Concordance of gene genealogies reveals reproductive isolation in the pathogenic fungus Coccidioides immitis Proc. Natl. Acad. Sci. USA 94:5478-5482 (see also 95:8414) [Abstract/Free Full Text]

    Kurtzman C. P., C. J. Robnet, 1998 Identification and phylogeny of ascomycetous yeasts and analysis of nuclear large subunit (26S) ribosomal DNA partial sequences Antonie Leeuwenhoek 73:331-371[ISI]

    Lachance M. A., W. M. Pang, 1997 Predacious yeasts Yeast 13:225-232[ISI][Medline]

    Lawrence J. G., J. R. Roth, 1996 Selfish operons: horizontal transfer may drive the evolution of gene clusters Genetics 143:1843-1860[Abstract/Free Full Text]

    Lohe A. R., E. N. Moriyama, D. A. Lidholm, D. L. Hartl, 1995 Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements Mol. Biol. Evol 12:62-72[Abstract]

    Maddison W. P., D. R. Maddison, 1992 MacClade 3.07 ed. Sinauer, Sunderland, Mass

    Marinoni G., M. Manuel, R. Fons Petersen, J. Hvidtfeldt, P. Sulo, J. Piskur, 1999 Horizontal transfer of genetic material among Saccharomyces yeasts J. Bacteriol 181:6488-6496[Abstract/Free Full Text]

    Mueller J. E., M. Bryk, N. Loizos, M. Belfort, 1993 Homing endonucleases Pp. 111–143 in S. M. Linn, R. S. Lloyd, and R. J. Roberts, eds. Nucleases, Vol. 2. Cold Spring Harbor Laboratory Press, Plainview, New York

    Naumov G. I., 1987 Genetic basis for classification and identification of the ascomycetous yeasts Stud. Mycol 30:469-475

    Nevoigt E., A. Fassbender, U. Stahl, 2000 Cells of the yeast Saccharomyces cerevisiae are transformable by DNA under non-artificial conditions Yeast 16:1107-1110[ISI][Medline]

    Pagel M., 1994 Detecting correlated evolution on phylogenies: a general method for the comparative analysis of discrete characters Proc. R. Soc. Lond. B 255:37-45[ISI]

    Perler P. B., G. J. Olsen, E. Adam, 1997 Compilation and analysis of intein sequences Nucleic Acids Res 25:1087-1093[Abstract/Free Full Text]

    Phalip V., I. Kuhn, Y. Lemoine, J.-M. Jeltsch, 1999 Characterization of the biotin biosynthesis pathway in Saccharomyces cerevisiae and evidence for a gene cluster containing BIO5, a novel gene involved in vitamer uptake Gene 232:43-51[ISI][Medline]

    Rosewich U. L., H. C. Kistler, 2000 Role of horizontal gene transfer in the evolution of fungi Annu. Rev. Phytopathol 38:325-363[ISI][Medline]

    Strathern J., D. Higgins, 1991 Recovery of plasmids from yeast into E. coli: shuttle vectors Pp. 322–323 in C. Guthrie and G. Fink, eds. Guide to yeast genetics and molecular biology. Academic Press, San Diego, Calif

    Swofford D. L., 1999 PAUP* phylogenetic analysis using parsimony 4.0d 64 ed. Sinauer, Sunderland, Mass

    Thompson J. D., D. G. Higgins, T. J. Gibson, 1994 Clustal W—improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice Nucleic Acids Res 22:4673-4680[Abstract]

    Wolfram S., 1996 The mathematica book Wolfram Media, Champaign

Accepted for publication October 9, 2001.