Evolutionary Maintenance of Selfish Homing Endonuclease Genes in the Absence of Horizontal Transfer
Total Page:16
File Type:pdf, Size:1020Kb
Evolutionary maintenance of selfish homing endonuclease genes in the absence of horizontal transfer Koji Yaharaa,b,c, Masaki Fukuyod, Akira Sasakie,f,g, and Ichizo Kobayashid,1 Divisions of aBiostatistics and bInfectious Diseases, Kurume University School of Medicine, Fukuoka 830-0011, Japan; cLife Science Systems Department, Fujitsu Kyushu Systems, Fukuoka 814-8589, Japan; dLaboratory of Social Genome Sciences, Department of Medical Genome Sciences, Graduate School of Frontier Science and Institute of Medical Science, University of Tokyo, Tokyo 108-8639, Japan; eDepartment of Evolutionary Studies of Biosystems (Sokendai-Hayama), The Graduate University for Advanced Studies (Sokendai), Hayama, Kanagawa 240-0193, Japan; fPRESTO, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan; and gEvolution and Ecology Program, International Institute for Applied Systems Analysis, A-2361 Laxenburg, Austria Edited by Simon A. Levin, Princeton University, Princeton, NJ, and approved August 25, 2009 (received for review July 26, 2009) Homing endonuclease genes are ‘‘selfish’’ mobile genetic elements particular site (3, 4). Specifically, they have been hypothesized to whose endonuclease promotes the spread of its own gene by persist over a long evolutionary period only if they can transmit creating a break at a specific target site and using the host to other species before extinction. For this reason, active homing machinery to repair the break by copying and inserting the gene at endonuclease genes are often used as an indicator of DNA this site. Horizontal transfer across the boundary of a species or transfer events. This interpretation is consistent with phyloge- population within which mating takes place has been thought to netic analysis of homing endonuclease genes such as and be necessary for their evolutionary persistence. This is based on the PI-SceI (VDE) of Saccharomyces cerevisiae (5), the group I assumption that they will become fixed in a host population, intron of the mitochondrial cox1 gene (6), and the group I intron where opportunities of homing will disappear, and become sus- of the 23S rRNA genes of hyper-thermophilic bacteria. ceptible to degeneration. To test this hypothesis, we modeled Although these interpretations are widely accepted, a func- behavior of a homing endonuclease gene that moves during tional homing intein endonuclease in the PRP8 gene of euas- meiosis through double-strand break repair. We mathematically comycetes (Pezizomycota) has been found that can survive within explored conditions for persistence of the homing endonuclease a species for several hundred million years, in the absence of gene and elucidated their parameter dependence as phase dia- interspecies horizontal transfer (4). Because the conditions that grams. We found that, if the cost of the pseudogene is lower than allow the evolutionary persistence of homing endonuclease that of the homing endonuclease gene, the 2 forms can persist in genes without horizontal transfer are unknown, this process must a population through autonomous periodic oscillation. If the cost be explored. In this work, we investigate if horizontal gene of the pseudogene is higher, 2 types of dynamics appear that transfer across the boundary of a species or population within BIOLOGY enable evolutionary persistence: bistability dependent on initial which mating takes place is necessary for the evolutionary POPULATION frequency or fixation irrespective of initial frequency. The predic- maintenance of homing endonuclease genes. We built a model tion of long persistence in the absence of horizontal transfer was of a hypothetical homing endonuclease gene that spreads confirmed by stochastic simulations in finite populations. The through meiotic double-strand break repair, similar to VDE (7). average time to extinction of the endonuclease gene was found to We analyzed its allele frequency over evolutionary periods in be thousands of meiotic generations or more based on realistic which an alternation of meiotic generations was interspersed parameter values. These results provide a solid theoretical basis for with stages of haploid and diploid mitosis, followed by degen- an understanding of these and other extremely selfish elements. erate mutation, resulting in a pseudogene, which would be lost. We mathematically explored conditions that allowed persistence ͉ ͉ ͉ ͉ horizontal gene transfer selfish gene VDE meiosis of the homing endonuclease gene, and numerically confirmed haploid-diploid cycle the evolutionary dynamics. We discovered conditions for persistence of the homing oming endonucleases are a large class of proteins that are endonuclease gene, elucidated their dependence on parameters Hcommon in simple eukaryotes and in prokaryotes (1). Their such as its cost, and generated phase diagrams. Persistence was coding regions (called genes here for convenience), are inserted possible through autonomous periodic oscillation, or through 2 into a specific DNA sequence and spliced out of either the other types of evolutionary dynamics, even in the absence of transcribed mRNA as an intron, or out of the translated protein horizontal transfer. These results were confirmed by stochastic as an intein. Homing endonuclease genes act as ‘‘selfish’’ mobile simulations in finite populations. genetic elements, because their gene product introduces a break at an ‘‘empty’’ target site and forces the host machinery to repair Model. We constructed a model of a haploid-diploid cycle in a the break by copying the endonuclease gene. These genes are not Mendelian population of unicellular eukaryotes, in which the essential to their hosts. Currently, much is known about the haploid cells correspond to gametes in higher eukaryotes. We molecular structure and action of homing endonucleases and focused on a homing endonuclease locus that moves by homing their genes, but little information is available on their population in the early phase of meiosis. At each round of meiosis, which was biology (2). From a population genetics standpoint, a homing endonucle- ase allele could easily become fixed in a host population and Author contributions: I.K. designed research; K.Y., M.F., and A.S. performed research; and susceptible to degeneration because homing opportunities dis- K.Y., A.S., and I.K. wrote the paper. appear. Therefore, the long-term survival of a homing endonu- The authors declare no conflict of interest. clease allele at a specific site is proposed to depend on horizontal This article is a PNAS Direct Submission. transfer across the boundary of a species or population, within 1To whom correspondence should be addressed. E-mail: [email protected]. which mating takes place, to another species or population that This article contains supporting information online at www.pnas.org/cgi/content/full/ do not yet contain the homing endonuclease allele at that 0908404106/DCSupplemental. www.pnas.org͞cgi͞doi͞10.1073͞pnas.0908404106 PNAS ͉ November 3, 2009 ͉ vol. 106 ͉ no. 44 ͉ 18861–18866 Downloaded by guest on September 25, 2021 x y z probability of the formation of diploid Hϩ/HϪ was 2xy, and that - + A H H Hd of Hϩ/Hd was 2yz. We assumed that the haploid cells could freely switch their mating type (8) or that the homing endonuclease syngamy alleles were randomly associated with the mating types. ϩ Ϫ 2xy 2yz …other 4 patterns During meiosis of an H /H diploid cell, the homing endo- nuclease produced from the Hϩ allele cleaves the unique empty target site (HϪ) of the homologous chromosome. Homing H+ H+ endonucleases are known to have a single target site in a haploid genome (9). Although they can evolve new target specificities H- Hd (10, 11), they also tolerate individual bp variation at their mutation recognition sequences (12). We therefore ignored the possibility r endonuclease that mutations at their target prevented recognition. The DNA double-strand break generated by endonuclease cleavage is efficiently repaired by the recombination repair system of the H+ host cell that is induced at the meiosis stage. Repair occurs with H+ no deleterious effects (13) through the copying of the homolo- meiosis gous Hϩ allele, which results in gene conversion from HϪ to Hϩ. This homing process takes place at the rate r. The homing endonuclease will not cleave the Hϩ or Hd allele because of insertion at the recognition site. The pseudogene Hd and the B empty target site HϪ do not produce the functional homing + - H u Hd v H endonuclease, so the homing process and the accompanying ϩ Ϫ Degeneration Loss of pseudogene change in allele frequency can occur only in H /H diploids. (pseudogene formation) After the homing step, each diploid cell that results from Fig. 1. Model of homing endonucleases. (A) Life cycle of homing endonu- production of a zygote was assumed to undergo meiosis to clease alleles: allele Hϩ produces a homing endonuclease that cleaves the regenerate haploid cells. The frequency of each allele was empty allele, HϪ, on its homologous chromosome, in a diploid cell during calculated at this stage. We assumed a cost to the host cell of ϩ meiosis. The DNA double-strand break is repaired by the host with the H carrying a homing endonuclease gene, c1, manifested as a allele as the template. This process of homing results in gene conversion from reduced growth rate. The relative fecundity of a cell carrying an HϪ to Hϩ. The homing endonuclease does not cleave Hd, which represents a ϩ Ϫ ϩ H allele to one carrying only an H empty target site depended pseudogene degenerated from H . Diploid cells are formed by syngamy from on the cost expressed as ␣ ϭ eϪc1. Also assumed was the cost of all pair-wise combinations of haploid cells, whose probabilities depend on the frequency of each haploid cell type. (B) Mutation: the Hϩ allele becomes carrying a pseudogene, c2, which was also represented by a d d d reduced growth rate. The relative fecundity of a cell carrying H pseudogene H by degenerate mutation at the rate u. Pseudogene H is Ϫ occasionally lost from its target site with the rate v, which is smaller than u.