During evolution, new species emerge when individuals from different populations of similar organisms no longer breed with each other, or when the offspring produced if they do breed are sterile. This process is known as “reproductive isolation” and, for over 100 years, evolutionary biologists have tried to better understand how this process happens.

Animals, plants and fungi produce sex cells – known as gametes – when they are preparing to reproduce. These cells are made when cells containing two copies of every gene in the organism divide to produce new cells that each only have one copy of each gene. Therefore, a particular gene copy usually has a 50% chance of being carried by an individual gamete. There are genes that selfishly increase their chances of being transmitted to the next generation by destroying the gametes that do not carry them. These “gamete killer” genes reduce the fertility of the organism and lead to reproductive isolation.

Fission yeast is a fungus that is widely used in research. There are different strains of fission yeast that are reproductively isolated from each other, but it is not known whether gamete killers are responsible for this isolation. To address this question, Hu et al. investigated the causes of reproductive isolation in fission yeast.

The experiments identified two gamete killers, referred to as cw9 and cw27. Both genes belong to the wtf gene family. Each gene is believed to encode two different proteins, one that acts as a poison and one that acts as an antidote. The poison is capable of killing all gametes, but the antidote protects the cells that contain the gamete killer gene. Further experiments show that the antidote produced by one of the gamete killer genes cannot protect cells against the poison produced by the other gene.

A separate study by Nuckolls et al. found that another member of the wtf family also acts as a gamete killer in fission yeast. Together, these findings shed new light on the causes of reproductive isolation, and will contribute to deeper understanding of speciation and evolution in general.

Mendel’s law of equal segregation stipulates that paternal and maternal alleles of a gene should have an equal chance of being transmitted to progenies. This law guarantees a fair competition between different alleles and enables beneficial ones to prevail during natural selection. Meiotic drivers, a type of selfish genetic element, break Mendel’s law by skewing transmission ratios to their advantage, and thus can spread in a population even when having a deleterious effect on organismal fitness (Lindholm et al., 2016; Werren, 2011). The term ‘meiotic drive’ was initially coined to describe segregation distortion resulting from preferential inclusion in the gamate during asymmetric female meiosis but has now been used more broadly to include biased transmission caused by postmeiotic mechanisms. In fact, some of the best-known meiotic drivers, such as Segregation Distorter in Drosophila and the t haplotype in mouse (Lyon, 2003; Larracuente and Presgraves, 2012), act postmeiotically to disable male gametes (sperms) that do not inherit them. This type of meiotic driver, called gamete killer, exists in animals, plants, and fungi (Burt and Trivers, 2006). Fungal gamete killers, or spore killers, have been found in several filamentous ascomycetes, most notably Neurospora and Podospora anserina (Dalstra et al., 2003; Grognet et al., 2014; Hammond et al., 2012; Turner and Perkins, 1979). It is unclear how widespread spore killers are among fungal species.

The fission yeast Schizosaccharomyces pombe is a prominent model organism for molecular and cell biology and has been increasingly used to study natural variation and genome evolution (Brown et al., 2011; Hu et al., 2015; Jeffares et al., 2015; Rhind et al., 2011). S. pombe natural isolates, which are nearly all haploids and have pair-wise nucleotide differences of less than 1% (Jeffares et al., 2015; Rhind et al., 2011), can readily mate with each other to form hybrid diploids, but the viability of spores derived from inter-isolate crosses is often below 5% and in many instances under 1% (Gutz and Doe, 1975; Kondrat’'eva and Naumov, 2001; Naumov et al., 2015). This is in stark contrast to Saccharomyces cerevisiae, whose natural isolates have a similar level of nucleotide diversity but much better spore viability when inter-crossed (Hou et al., 2014). One explanation for the within-species hybrid sterility in S. pombe is chromosomal rearrangement (Avelar et al., 2013; Brown et al., 2011; Zanders et al., 2014). However, because one rearrangement reduces spore viability at most by half (Avelar et al., 2013; Hou et al., 2014), other factors are likely in play. It was shown recently that when the S. pombe laboratory strain, which was isolated from French grape juice (Hu et al., 2015; Osterwalder, 1924), was crossed to a strain isolated from fermented tea (initially called Schizosaccharomyces kambucha [Singh and Klar, 2003], later renamed S. pombe var. kambucha [Rhind et al., 2011]), at least three spore killers contributed to hybrid sterility (Bomblies, 2014; Zanders et al., 2014), but the identities of the killer genes were unknown.

In this study, through investigating the causes of intraspecific hybrid sterility of fission yeast, we uncovered the molecular identities of two active spore killers, which both belong to the wtf gene family. We analyzed their killing behaviors in both native and non-native genomic contexts. Furthermore, we performed comparative genomic analysis of wtf genes and revealed interesting patterns of divergence among them.

We have previously studied CBS5557, an S. pombe strain isolated from grapes in Spain (Hu et al., 2015; Rankine and Fornachon, 1964). Like most S. pombe isolates, CBS5557 differs from the laboratory strain, whose genome is the S. pombe reference genome, in a 2.23 Mb pericentric inversion on chromosome I (Brown et al., 2011). After introducing the same inversion into the laboratory strain (Hu et al., 2015), the viability of spores derived from the cross between the laboratory strain and CBS5557 increased from 14% to 22% (Figure 1—figure supplement 1 and Supplementary file 1). To uncover the causes of the remaining sterility, we performed next-generation sequencing (NGS)-assisted bulk segregant analysis (Hu et al., 2015) and found that among the viable progenies from the cross, SNP alleles of the laboratory strain (reference alleles) on chromosome III were markedly under-represented (Figure 1A). Two broad transmission distortion peaks, one on the left arm (hereafter referred to as the left peak) and one on the right arm (hereafter referred to as the right peak), were observed, suggesting that there are at least two spore killer genes on chromosome III of CBS5557.

Figure 1 with 3 supplements see all

Download asset Open asset

To narrow down the genomic regions containing the spore killer genes, we performed multiple rounds of backcross and obtained strains with only a portion of chromosome III originating from CBS5557 and the rest of the genome coming from the laboratory strain (see 'Materials and methods'). We selected two backcrossed strains that retained the ability to skew inheritance of markers inserted in the left and right peak regions, respectively. These two strains, referred to as backcrossed-1 and backcrossed-2, respectively, were crossed to the laboratory strain and viable progenies were subjected to NGS-assisted bulk segregant analysis, which revealed the locations of the CBS5557 genomic segments present in the two backcrossed strains and confirmed the allele transmission bias (left-side plots in Figure 1B and C).

Through inspecting the Illumina and PacBio sequencing data of the CBS5557 genome, we discovered that, within the CBS5557 genomic segments in the two backcrossed strains, sequences most divergent from the reference genome are those containing wtf family genes (Figure 1D and E). wtf (for with Tf LTRs) genes were so named because they are often in close proximity to solo long terminal repeats (LTRs) derived from Tf retrotransposons (Bowen et al., 2003; Wood et al., 2002). wtf genes are predicted to contain multiple introns and encode multi-transmembrane proteins and have been shown to be highly up-regulated transcriptionally during meiosis and sporulation by large-scale analyses (Bowen et al., 2003). No homologs of wtf genes have been found in any other species, including other Schizosaccharomyces species (Rhind et al., 2011). We deleted from the backcrossed-1 strain a wtf gene locating at the position of wtf9 in the reference genome but exhibiting a low sequence identity (42%) to wtf9 in the 5’ portion (Figure 1D and Figure 1—figure supplement 2) and found that the deletion notably reduced the segregation distortion (right-side plot in Figure 1B), suggesting that this wtf gene, which we named cw9 (the nomenclature of wtf genes in the CBS5557 genome will be explained later in this paper), is a spore killer. The residual allele frequency bias is probably due to additional spore killer gene(s) in the CBS5557 genomic segments present in the backcrossed-1 strain. Within the CBS5557 genomic segment in the backcrossed-2 strain, we found a wtf gene absent in the reference genome, possibly owing to a deletion event mediated by homologous recombination between LTRs (Figure 1E and Figure 1—figure supplement 3). Removing this wtf gene, which we named cw27, from the backcrossed-2 strain completely abolished segregation distortion (right-side plot in Figure 1C).

To determine how cw9 and cw27 behave in their native strain background, we generated h⁺/h^- CBS5557 diploid strains containing homozygous or heterozygous deletion of cw9 or cw27, and performed tetrad analysis to evaluate spore viability (Figure 2A and Supplementary file 1). Homozygous deletion of cw9 or cw27 had no effect on spore viability (Figure 2A and Figure 2—figure supplement 1). In contrast, heterozygous deletion of either cw9 or cw27 caused a significant reduction of spore viability, which resulted mainly from the death of spores carrying the deletion (Figure 2B). When cw9 and cw27 were both in the heterozygous state, a more severe loss of spore viability was observed, with the double deletion mutant spores suffering a particularly high level of death (Figure 2A and B, and Supplementary file 1). These data indicate that cw9 and cw27 can act as spore killers in the self-cross of CBS5557 when either of them is in the heterozygous state, and inheriting one killer can prevent its killing effect but not that of the other killer.

Figure 2 with 1 supplement see all

Download asset Open asset

To determine whether cw9 and cw27 can cause spore killing in a non-native genomic context, we constructed integrating plasmids carrying cw9 or cw27, integrated the plasmids at various genomic locations in the laboratory strain, and subjected the resulting strains to crosses. In crosses where only one of the parental strains harbored cw9 or cw27, spores inheriting the killer had normal viability, whereas spores not inheriting the killer suffered significant loss of viability, regardless of the genomic location of the killer (Figure 3A and Figure 3—figure supplement 1). We then performed crosses in which both parental strains harbored a killer at the leu1 locus (Figure 3B and Figure 3—figure supplement 2). When both parental strains contained the same killer, spore viability was normal. In contrast, when one parental strain had cw9 and the other parental strain had cw27, spores inheriting either cw9 or cw27 experienced a dramatic loss of viability. These results demonstrate that cw9 and cw27 can act independently of genomic context and confirm that they do not exhibit mutual resistance.

Figure 3 with 2 supplements see all

Download asset Open asset

To observe the ultrastructural details of the spore killing process, we performed electron microscopy analysis on laboratory-background h⁺/h^- diploid cells undergoing synchronous meiosis and sporulation induced by a shift to a nitrogen-free synthetic sporulation medium (SSL-N medium) (Figure 3C). Two diploid strains were used: a strain heterozygous for a cw27-containing plasmid integrated at the his3 locus and a control strain heterozygous for an empty vector integrated at the same locus. 12 hr after the shift, cw27-containing cells formed four-spored asci indistinguishable from those formed by the control cells. At this time point, organelles inside spores were clearly visible. 18 hr after the shift, the cytoplasm of spores derived from the control cells became electron-dense and organelles were no longer discernable, owing to a little-understood spore maturation process (Yoo et al., 1973). Interestingly, at the same time point, asci derived from the cw27-containing cells often exhibited a 2:2 pattern of spore morphology: two spores in an ascus looked like the matured spores in the control sample, whereas the other two spores exhibited cytoplasm with a lower electron density and visible organelles, an electron-lucent spore wall, and frequently a non-spherical shape. These data suggest that killer-affected spores fail to undergo a proper maturation process.

cw9 and cw27 are both flanked by solo LTRs, and upstream of their predicted start codons they share a nearly identical 288-bp-long sequence (hereafter referred to as the conserved_up sequence) (Figure 4—figure supplement 1), which is also conserved upstream of the coding sequences of many other wtf genes (Bowen et al., 2003). To determine whether these and other sequence features are important for spore killing, we performed truncation analysis (Figure 4A and Figure 4—figure supplement 2). The truncated versions were named Ta, Tb, Tc, Td, and Te. Their effects on spore viability in a heterozygous cross of laboratory strains were examined. The Ta versions, which lack the LTRs, caused spore killing phenotypes similar to those caused by longer constructs containing the flanking LTRs, indicating that the LTRs are not important. Interestingly, the Tb versions missing the sequence between the upstream LTR and the conserved_up sequence caused spores containing the killer to suffer a moderate viability loss, suggesting that these truncated killer genes partially lost the ability to protect killer-containing spores from being harmed. Removing from the Ta versions, the sequence downstream of the predicted stop codon (Tc versions) did not have a notable effect. Removing from the Tc versions, the last 10 amino-acid-coding codons (Td versions) greatly diminished spore killing, whereas removing from the Tc versions the sequences upstream of the start codons (Te versions) resulted in more indiscriminate killing than that caused by the Tb versions. These results suggest that intact C termini are important for the killing activity, and the upstream sequences are important for the protecting activity.

Figure 4 with 3 supplements see all

Download asset Open asset

To determine whether the Td versions, which are largely devoid of the killing activity, can still protect, we performed crosses in which a strain with a Td-version-containing plasmid integrated at the ars1 locus was crossed to a strain with a Tc-version-containing plasmid integrated at the leu1 locus, and found that the Td versions of both cw9 and cw27 can effectively protect spores harboring them from being killed by the respective Tc-version killers (Figure 4B and Figure 4—figure supplement 3). Thus, perturbing the C termini resulted in separation-of-function mutants that retain the protecting activity but can no longer kill. Such kind of mutants, if present in natural populations, would behave similarly to the resistant alleles of other previously studied gamete killers, and may also represent an intermediate state toward the extinction of a killer (Burt and Trivers, 2006).

In the reference genome, there are 25 wtf genes at 20 locations (15 singletons and five tandem pairs) (Bowen et al., 2003) (Figure 5A and B). Using PacBio sequencing, we found that in the CBS5557 genome, these 20 locations are occupied by 29 wtf genes (12 singletons, seven pairs, and one triplet), with four locations containing an extra wtf gene compared to the reference genome (Figure 5A and B). We systematically named these 29 genes according to their syntenic relationship with the 25 wtf genes in the reference genome, with each name consisting of the prefix cw, a number (from 1 to 25), and a suffix for genes at locations where extra genes exist in the CBS5557 genome (Figure 5A and B). In addition, there are three CBS5557 singleton wtf genes without syntenic counterparts in the reference genome and we named them cw26, cw27, and cw28 according to their order in the genome (Figure 5A). Like cw27, cw26 and cw28 also appeared to have been lost in the laboratory strain through LTR-mediated recombination (Figure 5—figure supplement 1), suggesting that this is a common mode of wtf gene turnover.

Figure 5 with 5 supplements see all

Download asset Open asset

Figure 5—source data 1

MAFFT-aligned DNA sequences of 57 wtf genes of the reference genome and the CBS5557 genome.

The sequences include the conserved_up regions, predicted coding sequences, and associated introns.

https://doi.org/10.7554/eLife.26057.017

Download elife-26057-fig5-data1-v1.fasta

Figure 5—source data 2

Gene structure predictions for 57 wtf genes of the reference genome and the CBS5557 genome.

https://doi.org/10.7554/eLife.26057.018

Download elife-26057-fig5-data2-v1.txt

Figure 5—source data 3

MAFFT-aligned amino acid sequences of the predicted protein products of 57 wtf genes of the reference genome and the CBS5557 genome.

https://doi.org/10.7554/eLife.26057.019

Download elife-26057-fig5-data3-v1.fasta

Among the 12 CBS5557 singleton wtf genes with syntenic counterparts in the reference genome, cw9 and five others exhibit exceptionally high DNA sequence diversity (<85% identity) from their counterparts in the reference genome (Figure 1D, Figure 1—figure supplement 2, and Figure 5—figure supplement 2). Given that the genome-wide nucleotide difference between CBS5557 and the laboratory strain is only 3.1 SNPs/kb (Hu et al., 2015), wtf genes must have diverged at a much faster pace than average genes.

To examine the phylogenetic relationship among all 57 wtf genes of the two genomes, we performed a maximum likelihood analysis (Figure 5C). In the resulting phylogenetic tree, we used colored rectangles to highlight 16 pairs of genes, with each pair consisting of phylogenetic neighbors, defined as two genes separated by a single internal node that has an ultrafast bootstrap support value >= 95% (Figure 5C). Interestingly, only 11 of these 16 pairs conform to typical orthologous relationship where two members of a pair are syntenic (highlighted by green rectangles in Figure 5C), whereas the five other pairs each consist of two genes that are non-syntenic (highlighted by yellow rectangles in Figure 5C). This pattern suggests that ectopic gene conversion (also called non-allelic or interlocus gene conversion) may contribute to the inter-isolate difference of wtf genes (Chen et al., 2007; Petes and Hill, 1988). Interestingly, genes belonging to the same gene cluster (pair or triplet), despite being subject more strongly to the homogenizing effect of gene conversion due to their physical proximity, are without exception very distinct from each other. One possible explanation is that diversifying selection may have prevented sequence homogenization within gene clusters.

The substantial difference in the numbers and sequences of wtf gene between the laboratory strain and CBS5557 suggests that rapid sequence divergence and independent loss/gain of wtf genes have resulted in different S. pombe natural isolates harboring distinct sets of active spore killers, which become reproductive barriers during inter-isolate crosses and contribute to the extensive and severe hybrid sterility within this species.

Only a small number of gamete killers have been characterized at the molecular level. Two well-studied gamete killers in animals, Segregation Distorter in Drosophila and the t haplotype in mouse, both contain at least two key loci, a killer locus and a target locus (Lyon, 2003; Burt and Trivers, 2006; Larracuente and Presgraves, 2012; Lindholm et al., 2016; Bauer et al., 2012). In fungi, spore killers of Neurospora also appear to comprise of two genes, rsk (resistant to Spore killer) and a yet-to-be-identified killer gene, which are proposed to act together in a fashion similar to the toxin-antitoxin (TA) systems of bacterial plasmids, in which a stable toxin and a labile antitoxin together ensure that bacterial cells cannot lose the plasmid expressing them (Hammond et al., 2012; McLaughlin and Malik, 2017; Saupe, 2012). The involvement of more than one gene renders these gamete killer systems vulnerable to recombination that may separate the component genes, and as a consequence, they have all evolved recombination-suppressing features such as chromosomal inversions (Burt and Trivers, 2006; Harvey et al., 2014). On the other hand, two types of single-gene spore killers, [Het-s] prion and the Spok genes, have been found in the filamentous fungus Podospora anserina (Dalstra et al., 2003; Grognet et al., 2014). The killing activity of the [Het-s] prion requires the presence of the HET-S protein, whereas the Spok genes can act autonomously. Fission yeast wtf killers share two prominent features with the Spok killers: being composed of a single predicted coding sequence and acting in a self-sufficient manner. These features arguably make wtf and spok killers the most compact and efficient forms of gamete killers.

Despite the similarities they share, wtf and spok killers are very different entities. Firstly, wtf genes only exist in a single fission yeast species S. pombe, indicating that they probably have arisen after the divergence of S. pombe from its closest relatives S. octosporus and S. cryophilus around 120 million years ago (Rhind et al., 2011), and have not jumped across the species barriers. In contrast, Spok family genes have a patchy but widespread distribution among Pezizomycotina species (Grognet et al., 2014), whose last common ancestor dates to 400 million years ago (Prieto and Wedin, 2013), suggesting that Spok killers have originated much earlier than wtf killers or have been especially prone to transfer horizontally across species. Secondly, spok family genes encode soluble proteins with a C-terminal kinase-like domain (InterPro:IPR011009), whereas wtf genes are predicted to encode multi-transmembrane proteins, suggesting that the detailed molecular mechanisms of these two types of killers are likely to be quite distinct.

The single-gene nature of wtf killers begs the question: how does a single gene encode both the killing and the protecting activities? Our findings on the mutations that can disrupt one but not the other activity suggest the possibility that more than one type of expression products can be generated from a wtf killer gene. Consistent with this idea, Nuckolls et al. showed that the wtf4 gene from the S. pombe var. kambucha strain (Sk wtf4) acts as a spore killer by expressing two overlapping transcripts, with the shorter transcript encoding a spore-killing poison and the longer transcript encoding an antidote that protects the spores from the killing effect of the poison (Nuckolls et al., 2017). This dual-transcript model can satisfactorily explain why the Tb and Te forms of cw9 and cw27 lost protecting activities but not killing activities (Figure 4A). Presumably, the transcription of the longer isoforms, which encode antidote proteins initiated from a start codon in the predicted exon 1, partially relies on the sequence between the upstream LTR and the conserved_up sequence and absolutely requires the conserved_up sequence; in contrast, the transcription of the shorter isoforms, which encode poison proteins initiated from a downstream alternative start codon, is not affected by removing the upstream sequences.

We hypothesized that the transcription of the poison isoforms may depend on a special promoter residing within the predicted intron 1 of cw9 and cw27. Indeed, through inspecting the DNA sequence alignment (Figure 5—source data 1), we found that cw9, cw27, and 10 other wtf genes in the reference genome and the CBS5557 genome share a 150 bp conserved sequence in the predicted intron 1, and this sequence is closely upstream of the transcription start site of a shorter-than-predicted wtf23 transcript isoform uncovered by the Iso-Seq long-read RNA sequencing method (Figure 5—figure supplement 3) (Kuang et al., 2017). These 12 genes can be classified into two subtypes, which we call Intron-1-ATG genes (wtf4, wtf13, cw4, cw18, and cw23) and Exon-2-ATG genes (wtf19, wtf23, cw8a, cw9, cw21, cw22, and cw27), based on the locations of the first ATG codons downstream of the 150 bp sequence (Figure 5C and Figure 5—figure supplement 3). These ATG codons are all in frame with the predicted coding sequences. The within-intron-1 ATG codons in the 5 Intron-1-ATG genes are only 7 bp away from the predicted exon 2, and if used as start codons, would result in proteins beginning with three intron-1-coded amino acids (Met, Leu, and Ser). The two active spore killer genes found in the S. pombe var. kambucha strain by Nuckolls et al., Sk wtf4 and Sk wtf28, also possess this 150 bp sequence and can be classified as an Intron-1-ATG gene and an Exon-2-ATG gene, respectively. We propose that this 150 bp sequence is a hallmark of the ‘poison-and-antidote’ genes that can, at least potentially, generate both poison and antidote isoforms according to the dual-transcript model proposed by Nuckolls et al.

The other 45 wtf genes in the reference genome and the CBS5557 genome lack this 150 bp sequence. In addition, except for the eight most divergent genes (wtf7, cw7, wtf11, cw11, wtf14, cw14, wtf15, and cw15), they do not have an in-frame ATG codon in intron one or in the 5’-section of exon 2, suggesting that most of these genes are not able to generate a shorter isoform according to the dual-transcript model, and are thus ‘antidote-only’ genes. Alternatively, it is possible that they may generate two different protein products using other mechanisms, such as alternative splicing (Kuang et al., 2017), overlapping proteins translated in different frames (Duncan and Mata, 2014), and post-translational processing and modification.

Five wtf genes in the reference genome, wtf1, wtf2, wtf3, wtf22, and wtf24, are annotated by PomBase as pseudogenes (McDowall et al., 2015). Our analysis of the gene structures of the 57 wtf genes in the reference genome and the CBS5557 genome showed that in addition to these five genes, four other wtf genes in the reference genome (wtf6, wtf8, wtf12, and wtf17), and three wtf genes in the CBS5557 genome (cw3, cw12, and cw16), appear to have suffered pseudogenizing mutations (see Materials and methods and Figure 5—source data 2), which include loss of start codons (wtf3 and cw3), premature stop codons (wtf1, wtf12, cw12, and cw16), frameshifting indels (wtf6, wtf8, wtf22, and wtf24), and segmental deletions (wtf2 and wtf17). These ‘dead’ wtf genes presumably are not able to generate either antidote or poison proteins.

Most of the wtf genes in the reference genome and the CBS5557 genome are predicted to contain either 5 or 6 exons. Interestingly, the 6-exon genes, including cw9, cw27, and all but one of the other ten 150-bp-sequence-containing genes, clustered together in the DNA sequence-based phylogenetic tree (Figure 5C), and their protein products are predicted to be longer and contain two or three more transmembrane helices than the protein products of 5-exon genes (Figure 5—source data 3, Figure 5—figure supplement 4, and Figure 5—figure supplement 5). For each of the 12 150-bp-sequence-containing genes, the antidote protein product is approximately 50 amino acids longer than the poison protein product at the N-terminal soluble region. The 10 C-terminal amino acids missing from the protein products of the Td forms of cw9 and cw27 are downstream of the last transmembrane helix, and are strongly conserved among the predicted protein products of 6-exon genes (Figure 5—figure supplement 4), indicating that Td-like mutations do not occur naturally to a significant extent.

How do transmenbrane proteins act as killers? In the prokaryotic toxin-antitoxin (TA) systems, many toxins are transmembrane proteins and kill cells by damaging cellular membranes (Schuster and Bertram, 2013). The best-known examples include the type I TA toxins Hok and TisB, and the type V TA toxin GhoT (Cheng et al., 2014; Gerdes et al., 1986; Unoson and Wagner, 2008). We speculate that Wtf proteins may also possess abilities to disrupt the integrity of certain cellular membranes.

wtf genes were named after their association with the solo LTRs of the Tf retrotransposons (Wood et al., 2002). It was initially hypothesized that retrotransposition may aid the expansion of the wtf gene family (Wood et al., 2002). However, the presence of introns in wtf genes and the lack of sequence similarity between wtf-flanking LTRs indicate that wtf genes have not been retrotransposed (Bowen et al., 2003). Bowen et al. have put forward an explanation for the enrichment of LTRs near wtf genes: the regions flanking wtf genes may have been favored insertion targets of Tf retrotransposons, perhaps owing to the high transcription levels of wtf genes during meiosis and sporulation (Bowen et al., 2003). Our findings that LTRs flanking cw9 and cw27 are dispensable for spore killing support the idea that the wtf-LTR association does not benefit wtf genes but may rather reflect an exploitation of wtf genes by Tf retrotransposons for their own expression and/or inheritance benefit. Furthermore, we show that three singleton wtf genes (cw26, cw27, and cw28) appear to have been lost in the laboratory strain due to LTR-mediated recombination, suggesting that the wtf-LTR association promotes wtf gene loss and is detrimental to the long-term persistence of wtf genes. In addition, it may not be coincidental that both cw9 and cw27 require an approximately 90-bp-long sequence between the upstream LTR and the conserved_up sequence for the full protecting activity, raising the possibility that some wtf genes that lack such an intervening sequence (e.g. wtf4 and cw4) may have suffered functional loss due to a Tf insertion too close to the conserved_up sequence or may have become functionally dependent on the upstream LTR.

In terms of the number of genes per genome, wtf genes are probably the most evolutionarily successful gamete killers known to date. During a relatively short evolutionary time span, wtf gene family has emerged and become the largest species-specific gene family in S. pombe (Lespinet et al., 2002). The compactness (minimal sizes of fully functional cw9 and cw27 are both around 1.9 kb) and the self-sufficiency of the wtf killers likely have been advantageous to its family size expansion. Furthermore, given the dramatic divergence of wtf genes between the reference genome and the CBS5557 genome, and the mutual-killing behaviors of cw9 and cw27, it is possible that two identical copies of a wtf gene generated by a gene duplication event can undergo relatively rapid sequence diversification to become mutually non-resistant, so that their selective advantage is not diminished by resistance. Ectopic gene conversion is known to play a causal role in the sequence diversification of animal major histocompatibility complex (MHC) genes (Ohta, 1991; Takuno et al., 2008), plant disease-resistance genes (R genes) (Mondragon-Palomino and Gaut, 2005), and Plasmodium falciparum var genes (Claessens et al., 2014). Our phylogenetic analysis suggests that ectopic gene conversion may also contribute to the historical expansion and the ongoing dynamic change of the wtf gene family.

It is intriguing that most of the wtf genes locate on chromosome III, the smallest of the three fission yeast chromosomes. To explain this phenomenon, Bowen et al. proposed that the laboratory strain may have inherited only chromosome III but not the other two chromosomes from an ancestor strain in which wtf genes were abundant on all three chromosomes (Bowen et al., 2003). Our finding here of the same kind of chromosomal distribution bias in CBS5557, whose nucleotide diversity from the laboratory strain is higher on chromosome III than on the other two chromosomes (Hu et al., 2015), makes this explanation less likely. The only type of viable aneuploidy in fission yeast is chromosome III disomy (Niwa et al., 2006; Niwa and Yanagida, 1985), which occurs at a frequency of >1% during normal meiosis (Molnar et al., 1995). We speculate that the tolerance of chromosome III disomy may in some way be advantageous to the expansion of the wtf gene family on this chromosome.

Spore killers and other meiotic drivers, despite having only been found in a small number of taxa, are believed to exist in a wide range of species, where they are an important force influencing genome organization, gametogenesis, and speciation (Burt and Trivers, 2006; Lindholm et al., 2016; Werren, 2011). The molecular basis of their action and the evolutionary mechanisms behind their birth, spread, and extinction are scarcely known. The unveiling of the molecular identities of spore killers in S. pombe by this study and that of Nuckolls et al. (Nuckolls et al., 2017), promises to bring the power of a highly tractable model organism to the pursuit of the many mysteries of these evolutionary wonders.

1. 1. Hu W
   2. Suo F
   3. Du LL

   (2017) Schizosaccharomyces pombe Genome sequencing

   Publicly available at the NCBI BioProject (accession no: PRJNA358837).

   https://www.ncbi.nlm.nih.gov/bioproject/PRJNA358837

1. 1. Arcangioli B
   2. Klar AJ

   (1991)

   A novel switch-activating site (SAS1) and its cognate binding factor (SAP1) required for efficient mat1 switching in Schizosaccharomyces pombe

   The EMBO Journal 10:3025–3032.

   • PubMed
   • Google Scholar
2. 1. Avelar AT
   2. Perfeito L
   3. Gordo I
   4. Ferreira MG

   (2013) Genome architecture is a selectable trait that can be maintained by antagonistic pleiotropy

   Nature Communications 4:2235.

   https://doi.org/10.1038/ncomms3235

   • PubMed
   • Google Scholar
3. 1. Bauer H
   2. Schindler S
   3. Charron Y
   4. Willert J
   5. Kusecek B
   6. Herrmann BG

   (2012) The nucleoside diphosphate kinase gene Nme3 acts as quantitative trait locus promoting non-mendelian inheritance

   PLoS Genetics 8:e1002567.

   https://doi.org/10.1371/journal.pgen.1002567

   • PubMed
   • Google Scholar
4. 1. Bomblies K

   (2014) Cheaters divide and conquer

   eLife 3:e03371.

   https://doi.org/10.7554/eLife.03371

   • PubMed
   • Google Scholar
5. 1. Bowen NJ
   2. Jordan IK
   3. Epstein JA
   4. Wood V
   5. Levin HL

   (2003) Retrotransposons and their recognition of pol II promoters: a comprehensive survey of the transposable elements from the complete genome sequence of Schizosaccharomyces pombe

   Genome Research 13:1984–1997.

   https://doi.org/10.1101/gr.1191603

   • PubMed
   • Google Scholar
6. 1. Brown WR
   2. Liti G
   3. Rosa C
   4. James S
   5. Roberts I
   6. Robert V
   7. Jolly N
   8. Tang W
   9. Baumann P
   10. Green C
   11. Schlegel K
   12. Young J
   13. Hirchaud F
   14. Leek S
   15. Thomas G
   16. Blomberg A
   17. Warringer J

   (2011) A geographically diverse collection of schizosaccharomyces pombe isolates shows limited phenotypic variation but extensive karyotypic diversity

   G3 (Bethesda) 1:615–626.

   https://doi.org/10.1534/g3.111.001123

   • PubMed
   • Google Scholar
7. Book

   1. Burt A
   2. Trivers R

   (2006)

   Genes in Conflict: The Biology of Selfish Genetic Elements

   Cambridge, Mass: Belknap Press of Harvard University Press.

   • Google Scholar
8. 1. Chen JM
   2. Cooper DN
   3. Chuzhanova N
   4. Férec C
   5. Patrinos GP

   (2007) Gene conversion: mechanisms, evolution and human disease

   Nature Reviews Genetics 8:762–775.

   https://doi.org/10.1038/nrg2193

   • PubMed
   • Google Scholar
9. 1. Cheng HY
   2. Soo VW
   3. Islam S
   4. McAnulty MJ
   5. Benedik MJ
   6. Wood TK

   (2014) Toxin GhoT of the GhoT/GhoS toxin/antitoxin system damages the cell membrane to reduce adenosine triphosphate and to reduce growth under stress

   Environmental Microbiology 16:1741–1754.

   https://doi.org/10.1111/1462-2920.12373

   • PubMed
   • Google Scholar
10. 1. Claessens A
    2. Hamilton WL
    3. Kekre M
    4. Otto TD
    5. Faizullabhoy A
    6. Rayner JC
    7. Kwiatkowski D

    (2014) Generation of antigenic diversity in plasmodium falciparum by structured rearrangement of Var genes during mitosis

    PLoS Genetics 10:e1004812.

    https://doi.org/10.1371/journal.pgen.1004812

    • PubMed
    • Google Scholar
11. 1. Dalstra HJ
    2. Swart K
    3. Debets AJ
    4. Saupe SJ
    5. Hoekstra RF

    (2003) Sexual transmission of the [Het-S] prion leads to meiotic drive in Podospora anserina

    PNAS 100:6616–6621.

    https://doi.org/10.1073/pnas.1030058100

    • PubMed
    • Google Scholar
12. 1. Duncan CD
    2. Mata J

    (2014) The translational landscape of fission-yeast meiosis and sporulation

    Nature Structural and Molecular Biology 21:641–647.

    https://doi.org/10.1038/nsmb.2843

    • PubMed
    • Google Scholar
13. 1. Egel R
    2. Egel-Mitani M

    (1974) Premeiotic DNA synthesis in fission yeast

    Experimental Cell Research 88:127–134.

    https://doi.org/10.1016/0014-4827(74)90626-0

    • PubMed
    • Google Scholar
14. 1. Forsburg SL
    2. Rhind N

    (2006) Basic methods for fission yeast

    Yeast 23:173–183.

    https://doi.org/10.1002/yea.1347

    • PubMed
    • Google Scholar
15. 1. Gerdes K
    2. Bech FW
    3. Jørgensen ST
    4. Løbner-Olesen A
    5. Rasmussen PB
    6. Atlung T
    7. Boe L
    8. Karlstrom O
    9. Molin S
    10. von Meyenburg K

    (1986)

    Mechanism of postsegregational killing by the hok gene product of the parB system of plasmid R1 and its homology with the relF gene product of the E. coli relB operon

    The EMBO Journal 5:2023–2029.

    • PubMed
    • Google Scholar
16. 1. Grognet P
    2. Lalucque H
    3. Malagnac F
    4. Silar P

    (2014) Genes that bias mendelian segregation

    PLoS Genetics 10:e1004387.

    https://doi.org/10.1371/journal.pgen.1004387

    • PubMed
    • Google Scholar
17. 1. Gutz H
    2. Doe FJ

    (1975) On Homo- and heterothallism in Schizosaccharomyces pombe

    Mycologia 67:748–759.

    https://doi.org/10.2307/3758336

    • PubMed
    • Google Scholar
18. 1. Hammond TM
    2. Rehard DG
    3. Xiao H
    4. Shiu PK

    (2012) Molecular dissection of Neurospora spore killer meiotic drive elements

    PNAS 109:12093–12098.

    https://doi.org/10.1073/pnas.1203267109

    • PubMed
    • Google Scholar
19. 1. Harvey AM
    2. Rehard DG
    3. Groskreutz KM
    4. Kuntz DR
    5. Sharp KJ
    6. Shiu PK
    7. Hammond TM

    (2014) A critical component of meiotic drive in Neurospora is located near a chromosome rearrangement

    Genetics 197:1165–1174.

    https://doi.org/10.1534/genetics.114.167007

    • PubMed
    • Google Scholar
20. 1. Hess PN
    2. De Moraes Russo CA

    (2007) An empirical test of the midpoint rooting method

    Biological Journal of the Linnean Society 92:669–674.

    https://doi.org/10.1111/j.1095-8312.2007.00864.x

    • Google Scholar
21. 1. Hou J
    2. Friedrich A
    3. de Montigny J
    4. Schacherer J

    (2014) Chromosomal rearrangements as a major mechanism in the onset of reproductive isolation in Saccharomyces cerevisiae

    Current Biology 24:1153–1159.

    https://doi.org/10.1016/j.cub.2014.03.063

    • PubMed
    • Google Scholar
22. 1. Hu W
    2. Suo F
    3. Du LL

    (2015) Bulk Segregant Analysis reveals the genetic basis of a natural trait variation in Fission yeast

    Genome Biology and Evolution 7:3496–3510.

    https://doi.org/10.1093/gbe/evv238

    • PubMed
    • Google Scholar
23. Data

    1. Hu W
    2. Suo F
    3. Du L-L

    (authors) (2017) Schizosaccharomyces pombeGenome sequencing

    NCBI - BioProject. ID PRJNA358837.

    https://www.ncbi.nlm.nih.gov/bioproject/PRJNA358837
24. 1. Jeffares DC
    2. Rallis C
    3. Rieux A
    4. Speed D
    5. Převorovský M
    6. Mourier T
    7. Marsellach FX
    8. Iqbal Z
    9. Lau W
    10. Cheng TM
    11. Pracana R
    12. Mülleder M
    13. Lawson JL
    14. Chessel A
    15. Bala S
    16. Hellenthal G
    17. O'Fallon B
    18. Keane T
    19. Simpson JT
    20. Bischof L
    21. Tomiczek B
    22. Bitton DA
    23. Sideri T
    24. Codlin S
    25. Hellberg JE
    26. van Trigt L
    27. Jeffery L
    28. Li JJ
    29. Atkinson S
    30. Thodberg M
    31. Febrer M
    32. McLay K
    33. Drou N
    34. Brown W
    35. Hayles J
    36. Carazo Salas RE
    37. Ralser M
    38. Maniatis N
    39. Balding DJ
    40. Balloux F
    41. Durbin R
    42. Bähler J

    (2015) The genomic and phenotypic diversity of Schizosaccharomyces pombe

    Nature Genetics 47:235–241.

    https://doi.org/10.1038/ng.3215

    • PubMed
    • Google Scholar
25. 1. Katoh K
    2. Misawa K
    3. Kuma K
    4. Miyata T

    (2002) MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transform

    Nucleic Acids Research 30:3059–3066.

    https://doi.org/10.1093/nar/gkf436

    • PubMed
    • Google Scholar
26. 1. Katoh K
    2. Standley DM

    (2014) MAFFT: iterative refinement and additional methods

    Methods in Molecular Biology 1079:131–146.

    https://doi.org/10.1007/978-1-62703-646-7_8

    • PubMed
    • Google Scholar
27. 1. Kondrat'eva VI
    2. Naumov GI

    (2001) The phenomenon of spore killing in Schizosaccharomyces pombe hybrids

    Doklady Biological Sciences 379:385–388.

    https://doi.org/10.1023/A:1011624918673

    • PubMed
    • Google Scholar
28. 1. Kuang Z
    2. Boeke JD
    3. Canzar S

    (2017) The dynamic landscape of fission yeast meiosis alternative-splice isoforms

    Genome Research 27:145–156.

    https://doi.org/10.1101/gr.208041.116

    • PubMed
    • Google Scholar
29. 1. Käll L
    2. Krogh A
    3. Sonnhammer EL

    (2005) An HMM posterior decoder for sequence feature prediction that includes homology information

    Bioinformatics 21 Suppl 1:i251–i257.

    https://doi.org/10.1093/bioinformatics/bti1014

    • PubMed
    • Google Scholar
30. 1. Larracuente AM
    2. Presgraves DC

    (2012) The selfish segregation distorter gene complex of Drosophila melanogaster

    Genetics 192:33–53.

    https://doi.org/10.1534/genetics.112.141390

    • PubMed
    • Google Scholar
31. 1. Lespinet O
    2. Wolf YI
    3. Koonin EV
    4. Aravind L

    (2002) The role of lineage-specific gene family expansion in the evolution of eukaryotes

    Genome Research 12:1048–1059.

    https://doi.org/10.1101/gr.174302

    • PubMed
    • Google Scholar
32. 1. Li H
    2. Hou J
    3. Bai L
    4. Hu C
    5. Tong P
    6. Kang Y
    7. Zhao X
    8. Shao Z

    (2015) Genome-wide analysis of core promoter structures in Schizosaccharomyces pombe with DeepCAGE

    RNA Biology 12:525–537.

    https://doi.org/10.1080/15476286.2015.1022704

    • PubMed
    • Google Scholar
33. 1. Lindholm AK
    2. Dyer KA
    3. Firman RC
    4. Fishman L
    5. Forstmeier W
    6. Holman L
    7. Johannesson H
    8. Knief U
    9. Kokko H
    10. Larracuente AM
    11. Manser A
    12. Montchamp-Moreau C
    13. Petrosyan VG
    14. Pomiankowski A
    15. Presgraves DC
    16. Safronova LD
    17. Sutter A
    18. Unckless RL
    19. Verspoor RL
    20. Wedell N
    21. Wilkinson GS
    22. Price TA

    (2016) The ecology and evolutionary dynamics of meiotic drive

    Trends in Ecology and Evolution 31:315–326.

    https://doi.org/10.1016/j.tree.2016.02.001

    • PubMed
    • Google Scholar
34. 1. Lyon MF

    (2003) Transmission ratio distortion in mice

    Annual Review of Genetics 37:393–408.

    https://doi.org/10.1146/annurev.genet.37.110801.143030

    • PubMed
    • Google Scholar
35. 1. Matsuyama A
    2. Shirai A
    3. Yashiroda Y
    4. Kamata A
    5. Horinouchi S
    6. Yoshida M

    (2004) pDUAL, a multipurpose, multicopy vector capable of chromosomal integration in fission yeast

    Yeast 21:1289–1305.

    https://doi.org/10.1002/yea.1181

    • PubMed
    • Google Scholar
36. 1. Matsuyama A
    2. Shirai A
    3. Yoshida M

    (2008) A novel series of vectors for chromosomal integration in fission yeast

    Biochemical and Biophysical Research Communications 374:315–319.

    https://doi.org/10.1016/j.bbrc.2008.07.015

    • PubMed
    • Google Scholar
37. Book

    1. McDonald JH

    (2014)

    Handbook of Biological Statistics (3rd edn)

    Sparky House Publishing.

    • Google Scholar
38. 1. McDowall MD
    2. Harris MA
    3. Lock A
    4. Rutherford K
    5. Staines DM
    6. Bähler J
    7. Kersey PJ
    8. Oliver SG
    9. Wood V

    (2015) PomBase 2015: updates to the fission yeast database

    Nucleic Acids Research 43:D656–D661.

    https://doi.org/10.1093/nar/gku1040

    • PubMed
    • Google Scholar
39. 1. McLaughlin RN
    2. Malik HS

    (2017) Genetic conflicts: the usual suspects and beyond

    The Journal of Experimental Biology 220:6–17.

    https://doi.org/10.1242/jeb.148148

    • PubMed
    • Google Scholar
40. 1. Molnar M
    2. Bähler J
    3. Sipiczki M
    4. Kohli J

    (1995)

    The rec8 gene of Schizosaccharomyces pombe is involved in linear element formation, chromosome pairing and sister-chromatid cohesion during meiosis

    Genetics 141:61–73.

    • PubMed
    • Google Scholar
41. 1. Mondragon-Palomino M
    2. Gaut BS

    (2005) Gene conversion and the evolution of three leucine-rich repeat gene families in Arabidopsis thaliana

    Molecular Biology and Evolution 22:2444–2456.

    https://doi.org/10.1093/molbev/msi241

    • PubMed
    • Google Scholar
42. 1. Naumov GI
    2. Kondratieva VI
    3. Naumova ES

    (2015) [Hybrid sterility of the yeast Schizosaccharomyces pombe: genetic genus and many species in statu nascendi?]

    Microbiology 84:159–169.

    https://doi.org/10.1134/S0026261715010099

    • PubMed
    • Google Scholar
43. 1. Nguyen LT
    2. Schmidt HA
    3. von Haeseler A
    4. Minh BQ

    (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies

    Molecular Biology and Evolution 32:268–274.

    https://doi.org/10.1093/molbev/msu300

    • PubMed
    • Google Scholar
44. 1. Niwa O
    2. Yanagida M

    (1985) Triploid meiosis and aneuploidy in Schizosaccharomyces pombe: an unstable aneuploid disomic for chromosome III

    Current Genetics 9:463–470.

    https://doi.org/10.1007/BF00434051

    • Google Scholar
45. 1. Niwa O
    2. Tange Y
    3. Kurabayashi A

    (2006) Growth arrest and chromosome instability in aneuploid yeast

    Yeast 23:937–950.

    https://doi.org/10.1002/yea.1411

    • PubMed
    • Google Scholar
46. 1. Nuckolls NL
    2. Núñez MAB
    3. Eickbush MT
    4. Young JM
    5. Lange JJ
    6. Yu JS
    7. Zanders SE

    (2017) Wtf genes are prolific dual poison-antidote meiotic drivers

    eLife 6:e26033.

    https://doi.org/10.7554/eLife.26033

    • Google Scholar
47. 1. Ohta T

    (1991) Role of diversifying selection and gene conversion in evolution of Major histocompatibility complex loci

    PNAS 88:6716–6720.

    https://doi.org/10.1073/pnas.88.15.6716

    • PubMed
    • Google Scholar
48. 1. Olson LW
    2. Eden U
    3. Egel-Mitani M
    4. Egel R

    (1978) Asynaptic meiosis in fission yeast?

    Hereditas 89:189–199.

    https://doi.org/10.1111/j.1601-5223.1978.tb01275.x

    • Google Scholar
49. 1. Osterwalder A

    (1924)

    Schizosaccharomyces liquefaciens n. sp., Eine Gegen freie schweflige säure widerstandsfähige gärhefe

    Mitt Gebiete Lebensmittelunters Hyg 15:5–28.

    • Google Scholar
50. 1. Petes TD
    2. Hill CW

    (1988) Recombination between repeated genes in microorganisms

    Annual Review of Genetics 22:147–168.

    https://doi.org/10.1146/annurev.ge.22.120188.001051

    • PubMed
    • Google Scholar
51. 1. Prieto M
    2. Wedin M

    (2013) Dating the diversification of the major lineages of ascomycota (Fungi)

    PLoS One 8:e65576.

    https://doi.org/10.1371/journal.pone.0065576

    • PubMed
    • Google Scholar
52. 1. Rankine BC
    2. Fornachon JC

    (1964) Schizosaccharomyces malidevorans SP.N., a yeast decomposing l-malic acid

    Antonie Van Leeuwenhoek 30:73–75.

    https://doi.org/10.1007/BF02046704

    • PubMed
    • Google Scholar
53. 1. Rhind N
    2. Chen Z
    3. Yassour M
    4. Thompson DA
    5. Haas BJ
    6. Habib N
    7. Wapinski I
    8. Roy S
    9. Lin MF
    10. Heiman DI
    11. Young SK
    12. Furuya K
    13. Guo Y
    14. Pidoux A
    15. Chen HM
    16. Robbertse B
    17. Goldberg JM
    18. Aoki K
    19. Bayne EH
    20. Berlin AM
    21. Desjardins CA
    22. Dobbs E
    23. Dukaj L
    24. Fan L
    25. FitzGerald MG
    26. French C
    27. Gujja S
    28. Hansen K
    29. Keifenheim D
    30. Levin JZ
    31. Mosher RA
    32. Müller CA
    33. Pfiffner J
    34. Priest M
    35. Russ C
    36. Smialowska A
    37. Swoboda P
    38. Sykes SM
    39. Vaughn M
    40. Vengrova S
    41. Yoder R
    42. Zeng Q
    43. Allshire R
    44. Baulcombe D
    45. Birren BW
    46. Brown W
    47. Ekwall K
    48. Kellis M
    49. Leatherwood J
    50. Levin H
    51. Margalit H
    52. Martienssen R
    53. Nieduszynski CA
    54. Spatafora JW
    55. Friedman N
    56. Dalgaard JZ
    57. Baumann P
    58. Niki H
    59. Regev A
    60. Nusbaum C

    (2011) Comparative functional genomics of the fission yeasts

    Science 332:930–936.

    https://doi.org/10.1126/science.1203357

    • PubMed
    • Google Scholar
54. 1. Robinson O
    2. Dylus D
    3. Dessimoz C

    (2016) Phylo.io: interactive viewing and comparison of large phylogenetic Trees on the web

    Molecular Biology and Evolution 33:2163–2166.

    https://doi.org/10.1093/molbev/msw080

    • PubMed
    • Google Scholar
55. 1. Saupe SJ

    (2012) A fungal gene reinforces Mendel's laws by counteracting genetic cheating

    PNAS 109:11900–11901.

    https://doi.org/10.1073/pnas.1209748109

    • PubMed
    • Google Scholar
56. 1. Schuster CF
    2. Bertram R

    (2013) Toxin-antitoxin systems are ubiquitous and versatile modulators of prokaryotic cell fate

    FEMS Microbiology Letters 340:73–85.

    https://doi.org/10.1111/1574-6968.12074

    • PubMed
    • Google Scholar
57. 1. Singh G
    2. Klar AJ

    (2003) DNA sequence of the mat2,3 region of Schizosaccharomyces kambucha shares high homology with the corresponding sequence from Sz. pombe

    Yeast 20:1273–1278.

    https://doi.org/10.1002/yea.1042

    • PubMed
    • Google Scholar
58. 1. Stanke M
    2. Diekhans M
    3. Baertsch R
    4. Haussler D

    (2008) Using native and syntenically mapped cDNA alignments to improve de novo gene finding

    Bioinformatics 24:637–644.

    https://doi.org/10.1093/bioinformatics/btn013

    • PubMed
    • Google Scholar
59. 1. Stothard P

    (2000)

    The sequence manipulation suite: javascript programs for analyzing and formatting protein and DNA sequences

    BioTechniques 28:1102–1104.

    • PubMed
    • Google Scholar
60. 1. Sun LL
    2. Li M
    3. Suo F
    4. Liu XM
    5. Shen EZ
    6. Yang B
    7. Dong MQ
    8. He WZ
    9. Du LL

    (2013) Global analysis of fission yeast mating genes reveals new autophagy factors

    PLoS Genetics 9:e1003715.

    https://doi.org/10.1371/journal.pgen.1003715

    • PubMed
    • Google Scholar
61. 1. Takuno S
    2. Nishio T
    3. Satta Y
    4. Innan H

    (2008) Preservation of a pseudogene by gene conversion and diversifying selection

    Genetics 180:517–531.

    https://doi.org/10.1534/genetics.108.091918

    • PubMed
    • Google Scholar
62. 1. Turner BC
    2. Perkins DD

    (1979)

    Spore killer, a chromosomal factor in neurospora that kills meiotic products not containing it

    Genetics 93:587–606.

    • PubMed
    • Google Scholar
63. 1. Unoson C
    2. Wagner EG

    (2008) A small SOS-induced toxin is targeted against the inner membrane in Escherichia coli

    Molecular Microbiology 70:258–270.

    https://doi.org/10.1111/j.1365-2958.2008.06416.x

    • PubMed
    • Google Scholar
64. 1. Walker BJ
    2. Abeel T
    3. Shea T
    4. Priest M
    5. Abouelliel A
    6. Sakthikumar S
    7. Cuomo CA
    8. Zeng Q
    9. Wortman J
    10. Young SK
    11. Earl AM

    (2014) Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement

    PLoS One 9:e112963.

    https://doi.org/10.1371/journal.pone.0112963

    • PubMed
    • Google Scholar
65. 1. Waterhouse AM
    2. Procter JB
    3. Martin DM
    4. Clamp M
    5. Barton GJ

    (2009) Jalview Version 2--a multiple sequence alignment editor and analysis workbench

    Bioinformatics 25:1189–1191.

    https://doi.org/10.1093/bioinformatics/btp033

    • PubMed
    • Google Scholar
66. 1. Werren JH

    (2011) Selfish genetic elements, genetic conflict, and evolutionary innovation

    PNAS 108 Suppl 2:10863–10870.

    https://doi.org/10.1073/pnas.1102343108

    • PubMed
    • Google Scholar
67. 1. Wood V
    2. Gwilliam R
    3. Rajandream MA
    4. Lyne M
    5. Lyne R
    6. Stewart A
    7. Sgouros J
    8. Peat N
    9. Hayles J
    10. Baker S
    11. Basham D
    12. Bowman S
    13. Brooks K
    14. Brown D
    15. Brown S
    16. Chillingworth T
    17. Churcher C
    18. Collins M
    19. Connor R
    20. Cronin A
    21. Davis P
    22. Feltwell T
    23. Fraser A
    24. Gentles S
    25. Goble A
    26. Hamlin N
    27. Harris D
    28. Hidalgo J
    29. Hodgson G
    30. Holroyd S
    31. Hornsby T
    32. Howarth S
    33. Huckle EJ
    34. Hunt S
    35. Jagels K
    36. James K
    37. Jones L
    38. Jones M
    39. Leather S
    40. McDonald S
    41. McLean J
    42. Mooney P
    43. Moule S
    44. Mungall K
    45. Murphy L
    46. Niblett D
    47. Odell C
    48. Oliver K
    49. O'Neil S
    50. Pearson D
    51. Quail MA
    52. Rabbinowitsch E
    53. Rutherford K
    54. Rutter S
    55. Saunders D
    56. Seeger K
    57. Sharp S
    58. Skelton J
    59. Simmonds M
    60. Squares R
    61. Squares S
    62. Stevens K
    63. Taylor K
    64. Taylor RG
    65. Tivey A
    66. Walsh S
    67. Warren T
    68. Whitehead S
    69. Woodward J
    70. Volckaert G
    71. Aert R
    72. Robben J
    73. Grymonprez B
    74. Weltjens I
    75. Vanstreels E
    76. Rieger M
    77. Schäfer M
    78. Müller-Auer S
    79. Gabel C
    80. Fuchs M
    81. Düsterhöft A
    82. Fritzc C
    83. Holzer E
    84. Moestl D
    85. Hilbert H
    86. Borzym K
    87. Langer I
    88. Beck A
    89. Lehrach H
    90. Reinhardt R
    91. Pohl TM
    92. Eger P
    93. Zimmermann W
    94. Wedler H
    95. Wambutt R
    96. Purnelle B
    97. Goffeau A
    98. Cadieu E
    99. Dréano S
    100. Gloux S
    101. Lelaure V
    102. Mottier S
    103. Galibert F
    104. Aves SJ
    105. Xiang Z
    106. Hunt C
    107. Moore K
    108. Hurst SM
    109. Lucas M
    110. Rochet M
    111. Gaillardin C
    112. Tallada VA
    113. Garzon A
    114. Thode G
    115. Daga RR
    116. Cruzado L
    117. Jimenez J
    118. Sánchez M
    119. del Rey F
    120. Benito J
    121. Domínguez A
    122. Revuelta JL
    123. Moreno S
    124. Armstrong J
    125. Forsburg SL
    126. Cerutti L
    127. Lowe T
    128. McCombie WR
    129. Paulsen I
    130. Potashkin J
    131. Shpakovski GV
    132. Ussery D
    133. Barrell BG
    134. Nurse P
    135. Cerrutti L

    (2002) The genome sequence of Schizosaccharomyces pombe

    Nature 415:871–880.

    https://doi.org/10.1038/nature724

    • PubMed
    • Google Scholar
68. 1. Yoo BY
    2. Calleja GB
    3. Johnson BF

    (1973) Ultrastructural changes of the fission yeast (Schizosaccharomyces pombe) during ascospore formation

    Archiv FuR Mikrobiologie 91:1–10.

    https://doi.org/10.1007/BF00409533

    • PubMed
    • Google Scholar
69. 1. Zanders SE
    2. Eickbush MT
    3. Yu JS
    4. Kang JW
    5. Fowler KR
    6. Smith GR
    7. Malik HS

    (2014) Genome rearrangements and pervasive meiotic drive cause hybrid infertility in fission yeast

    eLife 3:e02630.

    https://doi.org/10.7554/eLife.02630

    • PubMed
    • Google Scholar

Author details

1. Wen Hu

   National Institute of Biological Sciences, Beijing, China

   Contribution

   WH, Conceptualization, Investigation, Writing—original draft

   Competing interests

   WH: A patent application related to this work has been filed

   "This ORCID iD identifies the author of this article:" 0000-0002-0604-2119

2. Zhao-Di Jiang

   1. National Institute of Biological Sciences, Beijing, China
   2. PTN Graduate Program, School of Life Sciences, Tsinghua University, Beijing, China

   Contribution

   Z-DJ, Investigation

   Competing interests

   No competing interests declared.

3. Fang Suo

   National Institute of Biological Sciences, Beijing, China

   Contribution

   FS, Formal analysis, Writing—original draft

   Competing interests

   FS: A patent application related to this work has been filed

4. Jin-Xin Zheng

   National Institute of Biological Sciences, Beijing, China

   Contribution

   J-XZ, Investigation

   Competing interests

   No competing interests declared.

5. Wan-Zhong He

   National Institute of Biological Sciences, Beijing, China

   Contribution

   W-ZH, Formal analysis

   Competing interests

   No competing interests declared.

6. Li-Lin Du

   National Institute of Biological Sciences, Beijing, China

   Contribution

   L-LD, Conceptualization, Funding acquisition, Writing—original draft, Writing—review and editing

   For correspondence

   dulilin@nibs.ac.cn

   Competing interests

   L-LD: A patent application related to this work has been filed

   "This ORCID iD identifies the author of this article:" 0000-0002-1028-7397

• 4,096

  views

• 645

  downloads

• 82

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.