A Meiotic Drive Element in the Maize Pathogen Fusarium verticillioides Is Located Within a 102 kb Region of Chromosome V

Fusarium verticillioides is an agriculturally important fungus because of its association with maize and its propensity to contaminate grain with toxic compounds. Some isolates of the fungus harbor a meiotic drive element known as Spore killer (SkK) that causes nearly all surviving meiotic progeny from an SkK × Spore killer-susceptible (SkS) cross to inherit the SkK allele. SkK has been mapped to chromosome V but the genetic element responsible for meiotic drive has yet to be identified. In this study, we used cleaved amplified polymorphic sequence markers to genotype individual progeny from an SkK × SkS mapping population. We also sequenced the genomes of three progeny from the mapping population to determine their single nucleotide polymorphisms. These techniques allowed us to refine the location of SkK to a contiguous 102 kb interval of chromosome V, herein referred to as the Sk region. Relative to SkS genotypes, SkK genotypes have one extra gene within this region for a total of 42 genes. The additional gene in SkK genotypes, herein named SKC1 for Spore Killer Candidate 1, is the most highly expressed gene from the Sk region during early stages of sexual development. The Sk region also has three hyper-variable regions, the longest of which includes SKC1. The possibility that SKC1, or another gene from the Sk region, is an essential component of meiotic drive and spore killing is discussed.

was correlated with an outbreak of equine leukoencephalomalacia (ELEM) in South Africa (Kellerman et al. 1972). Since then, fumonisins have been confirmed to cause ELEM in horses Kellerman et al. 1990), pulmonary edema in pigs (Harrison et al. 1990), and liver cancer in laboratory rats (Gelderblom et al. , 1991. Fumonisin-contaminated grain has also been linked to neural tube defects in humans (Missmer et al. 2006). Thus, the potential presence of fumonisins, even in apparently healthy grain, makes it necessary to aggressively screen grain for these mycotoxins. Contaminated grain must be destroyed, resulting in crop losses worth millions of dollars per yr (Wu 2007).
The primary sources of inoculum for endophytic and pathogenic growth of F. verticillioides on maize are asexual spores known as macroconidia and microconidia (Munkvold 2003). Sexual spores, called ascospores, may also be important to the F. verticillioides life cycle (Reynoso et al. 2009). Sexual reproduction in this heterothallic fungus begins when the immature fruiting bodies of one strain are fertilized by a strain of opposite mating type (Martin et al. 2011). Fruiting bodies are called perithecia and meiosis occurs in spore sacs, called asci, within the perithecia. At the start of meiosis, a single diploid nucleus is formed from the haploid genomes of both parents. This diploid nucleus separates into four haploid nuclei through the standard reductional and equational division stages of meiosis, and each of the four haploid nuclei undergoes a postmeiotic mitosis to form eight nuclei. Each nucleus is then incorporated into an ascospore through a process known as ascospore delimitation. A phenotypically normal ascus should thus contain eight ascospores at maturity in F. verticillioides.
During normal development, asci break down and ascospores are exuded from perithecia in gelatinous masses. It is possible, however, to examine ascospores within intact asci by perithecial dissection. Kathariou and Spieth (1982) used this technique during their discovery of an allele named Spore killer (Sk K ). In Sk K · Sk S crosses, where the latter stands for Spore killer-susceptible, asci contain four instead of eight ascospores (Kathariou and Spieth 1982). While early stages of ascus development are cytologically normal in these crosses, four ascospores degenerate shortly after ascospore delimitation (Raju 1994). The four surviving ascospores are almost always of the Sk K genotype (Kathariou and Spieth 1982;Xu and Leslie 1996). Unlike Sk K · Sk S crosses, Sk K · Sk K crosses produce asci with eight ascospores (Kathariou and Spieth 1982). The Sk K allele is thus sufficient for both spore killing and resistance to killing.
The ability of Sk K -ascospores to kill Sk S -ascospores suggests that Sk K could be the predominant allele in some F. verticillioides populations. Although the literature is limited on this subject, a screen of 225 F. verticillioides isolates from 24 fields across Europe and North America found the Sk K allele to be present in 81% of isolates (Kathariou and Spieth 1982). Shortly after the discovery of Sk K in F. verticillioides, a spore killing allele was also discovered in F. subglutinans (Sidhu 1984). While this allele was also named Spore killer, it was given the slightly different notation of SK k . Sidhu (1984) presented evidence that SK k was present in 10 of 15 F. subglutinans isolates obtained from a single maize field. Therefore, Sk K and SK k appear to be the predominant alleles in at least some populations of F. verticillioides and F. subglutinans.
The relationship between Sk K and SK k is still unclear. To our knowledge, the only other significant research involving either allele was performed by Xu and Leslie (1996) during construction of an F. verticillioides genetic map. During this study, Sk K was mapped between two restriction fragment length polymorphism (RFLP) markers on chromosome V. These markers, named RFLP1 and 11p18, are located 2.5 cM and 8.6 cM from Sk K , respectively (Xu and Leslie 1996). Xu and Leslie's mapping population was derived from an Sk K · Sk S cross, and, interestingly, only one of more than 100 progeny did not inherit the Sk K allele (Xu and Leslie 1996), demonstrating that Sk K can achieve transmission rates of over 99% in laboratory crosses.
An F. verticillioides reference genome, derived from an Sk S strain known as Fv149-Sk S , was published 14 yr after the initial mapping of Sk K (Ma et al. 2010). Here, we advance Fusarium Spore killer research by delineating physical borders for Sk K with respect to this reference genome. Our data place Sk K within a 102 kb contiguous sequence of DNA, which we refer to as the Sk region. Notable differences exist between this region in the Sk K and Sk S strains examined in this study. These are described and discussed below.

MATERIALS AND METHODS
Strains, media, and culture conditions Key strains are listed in Table 1. Vegetative propagation was performed on V8 juice agar (VJA) (Tuite 1969) in 16 mm test tubes at room temperature on a laboratory benchtop. Carrot agar (CA), which was originally described by Klittich and Leslie (1988), was prepared as follows: 200 g of organic peeled baby cut carrots were autoclaved in 200 ml of water, pureed with a blender, then adjusted to 500 ml with sterile water to create a 1· stock. Diluted CA (e.g., 0.1· and 0.25·) was prepared by mixing the appropriate volumes of 1· CA stock and sterile water before adding agar to a final concentration of 2% and autoclaving. Liquid Vogel's Minimal medium (VMM, Vogel 1956) or GYP medium (2% glucose, 1% peptone, and 0.3% yeast extract) were used to produce mycelia for genomic DNA isolation. Liquid cultures for genomic DNA isolation were incubated at 28°in the dark without agitation.

Sexual crosses
Crosses were performed on CA in a manner similar to previously described methods (Klittich and Leslie 1988). Asexual spores (conidia) from the female parent were transferred to the center of a 60 mm petri dish containing 20 ml of CA or diluted CA. The female parent was then cultured for 10214 d before fertilization with a suspension of conidia from the male parent. Conidial suspensions were prepared by adding 2.0 ml of 0.001% Tween-20 to a 10-to 14-d-old test-tube culture of the male parent and dislodging the conidia with a pipette tip. Fertilization was performed by transferring 1.0 ml of this conidial suspension to the surface of a culture of the female parent. The conidial suspension was spread over the surface of the female culture with a glass rod. Crosses were incubated in a culture chamber that alternated between 23.0°n a Total, the total length of each genome assembly. b Avg. length, the average length of each read in the assembly. c Coverage, calculated by multiplying the number of reads by the average read length and dividing the product by 41.9 Mb, the genome size of Fv149-Sk S (NCBI ASM14955v1).
The Fv999-Sk K 3 Fv149-Sk S mapping population In F. verticillioides, ascospores are exuded from mature perithecia in a hair-like structure called a cirrus. To obtain an Sk K mapping population, Fv999-Sk K was crossed with Fv149-Sk S and cirri were isolated from the tops of a few perithecia with a sterile needle, dispersed in sterile water, and spread onto a plate of 4% water agar. Germinating ascospores were transferred to VJA in 16 mm test tubes.

Microscopy
Asci were dissected from perithecia in 25% glycerol under magnification. A Vanguard 1433Phi light microscope and attached digital camera (Amscope MU1000) were used for imaging. The condenser and aperture diaphragm were set for high contrast, which allowed for the number of ascospores in mature asci to be determined without tissue staining.

DNA methods
Genomic DNA was prepared using one of three methods. In method one, strains were cultured in 25 ml liquid VMM at 28°in the dark for 3 d. Mycelia were washed with 0.9% NaCl and dried by lyophilization before extraction with IBI Scientific's Genomic DNA Mini Kit for Plants. Method two, based on Henderson et al. (2005), was used as an inexpensive and time-efficient alternative to method one. A 6-inch plain-tipped wood applicator was used to transfer #10 mg of conidia to 200 ml of TE buffer (10 mM Tris, 1 mM EDTA, pH 8.0). The suspension was boiled at 105°in a heat block for 12 min, incubated on ice for 2 min, then vortexed for 5 sec. Insoluble material was pelleted at 15,000 · g for 10 min at room temperature, after which 25 ml of supernatant was transferred to a new vial and frozen at 220°for storage. In our hands, this method works well for the amplification of polymerase chain reaction (PCR) products shorter than 1 kb. Method three was used for the preparation of genomic DNA for high-throughput sequencing. Strains were cultured in 25 ml of liquid GYP at 28°for 48 hr. Mycelia were washed with water and DNA was extracted with the Zymo Research Fungal DNA Miniprep Kit.
PCR assays were performed with MidSci Bullseye Taq DNA Polymerase or New England Biolabs Phusion High-Fidelity DNA Polymerase.
Genome sequencing DNA libraries for MiSeq sequencing (Illumina) were constructed from 1 ng of genomic DNA using the Illumina Nextera XT DNA Library Preparation Kit. Sequencing was performed with MiSeq Reagent Kit Version 3. Adapters were removed and low-quality reads were trimmed with CLC Genomics Workbench (Version 8.0). The datasets were deposited in the National Center for Biotechnology Information's (NCBI) Sequence Read Archive (Leinonen et al. 2011). They can be obtained with the following accession numbers: SRR3271586 (Fv999-Sk K ), SRR3273544 (JP98.75-Sk K ), SRR3273545 (JP98.111-Sk K ), and SRR3273546 (JP98.118-Sk K ). The datasets are of high quality. For example, draft genomes were assembled with CLC Genomics Workbench and all assemblies had N50 values over 90 kb with coverage levels between 49-and 73-fold (Table 2).

CAPS markers
Cleaved amplified polymorphic sequence (CAPS) markers (Konieczny and Ausubel 1993) were used to help refine the location of Sk K (Table 3). CAPS markers were identified by first aligning genome sequencing reads from strain Fv999-Sk K to the Fv149-Sk S reference genome (NCBI, ASM14955v1) with Bowtie 2 (Langmead and Salzberg 2012) and then visually scanning aligned reads with Tablet (Milne et al. 2010) for polymorphisms in GGCC sites. This four-base sequence is cleaved by the restriction endonuclease HaeIII. Six CAPS markers on chromosome V were chosen for this study, along with one CAPS marker on each of chromosomes I, VII, and XI (Table 3). PCR primers for each CAPS marker were designed to amplify a short product (,500 bp) from both Fv999-Sk K and Fv149-Sk S sequences. The PCR primers used for amplification of each CAPS marker are described in Supplemental Material, Table S1.
The following protocol was used to analyze the segregation patterns of CAPS markers in each individual of the Fv999-Sk K · Fv149-Sk S mapping population. Genomic DNA was isolated from each progeny, CAPS markers were amplified by PCR, and PCR products were digested with HaeIII. The digested products were then examined for Fv999-Sk K or Fv149-Sk S cleavage patterns by gel electrophoresis on 2% agarose-TAE (40 mM Tris, 20 mM acetic acid, and 1 mM EDTA) gels. The HaeIII-based DNA digest patterns of each CAPS marker for Fv999-Sk K and Fv149-Sk S alleles are listed in Table S2.

Single nucleotide polymorphisms (SNPs)
Reads from each MiSeq dataset were aligned to the reference genome of strain Fv149-Sk S with Bowtie 2 (Langmead and Salzberg 2012). SAMtools 1.3 and BCFtools 1.3 were then used to report SNPs in variant call format (VCF) (Li et al. 2009;Danecek et al. 2011). Only reads that aligned to a single region of Fv149-Sk S with ,10 n  mismatches were used to produce the VCF files. Custom Perl scripts were then used to extract significant SNPs from the VCF files. Significant SNPs were defined as those which were supported by at least 90% of reads at each position. Only positions covered by more than nine reads were considered. Polymorphisms caused by insertions or deletions were ignored.

Gene predictions
Protein-coding genes within the Sk region of Fv999-Sk K were predicted as follows: first, the sequence of the Sk region in Fv999-Sk K was obtained by de novo assembly of MiSeq reads (Table 2); second, the coding sequences of annotated genes from the Sk region of Fv149-Sk S were obtained from GenBank (CM000582.1); third, the Fv149-Sk S coding sequences were aligned to the Sk region of Fv999-Sk K with Clustal Omega (Sievers et al. 2011;Li et al. 2015); and fourth, the alignments were used to manually annotate proteincoding sequences within the Fv999-Sk K Sk region. These steps identified 41 putative protein-coding genes within the Sk region of Fv999-Sk K . Augustus 3.2.1 (Stanke and Waack 2003) was then used for de novo prediction of protein-coding genes. Augustus predictions matched our manual annotation with the notable exception of an additional gene, SKC1 (described below), within the Sk region of Fv999-Sk K . The complete sequence and annotation of the 102 kb Sk K region from Fv999-Sk K can be downloaded from GenBank with accession number KU963213. Sikhakolli et al. (2012) analyzed transcriptional changes in F. verticillioides during fruiting body development in Fv999-Sk K · Fv149-Sk S crosses by RNA sequencing (RNAseq) and deposited the datasets in NCBI's SRA database (Leinonen et al. 2011). The following datasets were downloaded from SRA: 2 hr post fertilization (hpf) (SRR1592416), 24 hpf (SRR1592417), 48 hpf (SRR1592418), 72 hpf (SRR1592419), 96 hpf (SRR1592420), and 144 hpf (SRR1592421). Reads were aligned to coding sequences of each predicted gene from the Sk region, plus flanking genes FVEG_03199 and FVEG_03163 , with Bowtie 2 (Langmead and Salzberg 2012). Reads per kilobase exon model per thousand mapped reads (RPKK), a variation upon RPKM as described by Mortazavi et al. (2008), were calculated for each coding sequence. Because of sequence differences between Sk K -linked and Sk S -linked alleles of genes in the Sk region, alignments and RPKK calculations were performed separately for Sk K -linked and Sk S -linked alleles. Thus, reported RPKK values are averages except for SKC1, which is only found in Sk K genotypes.

Data availability
Data deposition: genome sequencing data are available through NCBI's Sequence Read Archive under accession numbers SRR3271586, SRR3273544, SRR3273545, and SRR3273546. Other sequence data are available from GenBank (KU963213). The annotated chromosome V sequence was downloaded from NCBI (CM000582.1). The positions of the first (FVEG_03430 ) and last (FVEG_16592 ) genes on the chromosome are indicated in the diagram. The approximate location of RFLP marker 11p18 is also shown. Marker 11p18 overlaps most of gene FVEG_02851. (B) A total of 60 progeny were isolated from a cross between Fv999-Sk K and Fv149-Sk S . These progeny were genotyped with nine different CAPS markers. Markers CAPS-1 through CAPS-6 are located on chromosome V. Markers CAPS-9, CAPS-10, and CAPS-11 are found on chromosomes XI, VII, and I, respectively. Each bar represents the percentage of Fv999-Sk K genotypes recovered for each marker. An asterisk is placed above all markers whose recovery deviated significantly from a Mendelian ratio of 1:1 according to a x 2 test with P , 0.01. The biased transmission of Fv999-Sk K sequences was detected for all markers on chromosome V, except for CAPS-5. (C) The marker patterns for progeny JP98.75-Sk K , JP98.111-Sk K , and JP98.118-Sk K suggest that cross-over events occurred between CAPS-3 and CAPS-4 (for JP98.75-Sk K ), CAPS-4 and CAPS-5 (for JP98.111-Sk K ), and CAPS-2 and CAPS-3 (for JP98.118-Sk K ) (Table S3). These cross-overs did not identify an Fv999-Sk K -derived region of chromosome V common to all three progeny. The predicted Fv999-Sk K -and Fv149-Sk S -inherited regions are shown in black and white, respectively. Irregular borders are used to indicate that exact cross-over positions could not be determined from the CAPS marker data. (D) SNP analysis was used to more accurately define cross-over positions for progeny JP98.75-Sk K , JP98.111-Sk K , and JP98.118-Sk K . In addition to refining the location of the cross-overs identified by CAPS analysis, two additional cross-overs were identified between CAPS-3 and CAPS-4 on chromosome V of JP98.111-Sk K . These findings helped delineate a single contiguous interval between position 665.7 kb and 767.4 kb on chromosome V as the only Fv999-Sk K -derived chromosome V sequence common to all three progeny (red dotted lines). Therefore, this interval is referred to as the Sk region.

RESULTS
Increased production of fruiting bodies by strain Fv999-Sk K on diluted CA Directional crosses are often used to investigate sexual reproduction in heterothallic ascomycete fungi. In this type of cross, one parent is designated as the female and the other is designated as the male. The female provides protoperithecia, i.e., immature fruiting bodies, which are fertilized by conidia from the male. After fertilization, protoperithecia develop into perithecia.
With F. verticillioides, CA is commonly used as a growth medium to study sexual processes. Its preparation essentially involves purchasing carrots from a local market, followed by peeling, blending, and autoclaving a predetermined weight of carrot within a specified volume of water (Leslie and Summerell 2006). Because we were curious about how the amount of carrot in CA influences the productivity of Fv999-Sk K · Fv149-Sk S crosses, Fv999-Sk K was cultured on 0.001·, 0.01·, 0.1·, 0.25·, 0.5·, and 1.0· CA, then fertilized with Fv149-Sk S conidia. Crosses performed on 0.1· and 0.25· CA resulted in three-to six-fold more perithecia than crosses performed on 0.5· and 1.0· CA (Figure 1). No perithecia were produced when crosses were performed on 0.001· and 0.01· CA (Figure 1). Additionally, less conidia were produced on 0.1· CA than on 0.25· CA (data not shown). Since conidia can be a source of cross-contamination, 0.1· CA was used as the crossing medium for the remainder of this study.
RFLP 11p18 overlaps with FVEG_02851 on chromosome V Although the sequence of RFLP1 was not available, we were able to obtain the sequence of RFLP 11p18 from the Fusarium Comparative Database. We used this sequence as the query in a Basic Local Alignment Search Tool (BLAST) (Altschul et al. 1990) search of the F. verticillioides reference genome. This search identified positions 1,672,753 to 1,673,602 on chromosome V (CM000582.1) as a match to 11p18. This region overlaps with most of the predicted coding sequence of gene FVEG_02851 , which spans positions 1,672,886 to 1,673,739 and encodes a protein of unknown function. This position on chromosome V was used as a reference point to design CAPS markers near the Sk locus ( Figure 2A).

Loci linked to Sk K drive through meiosis
To begin refining the location of Sk K on chromosome V, Fv999-Sk K was crossed with Fv149-Sk S and the segregation patterns of nine CAPS markers were examined in 60 progeny. Both Fv999-Sk K and Fv149-Sk S patterns were inherited in a 1:1 ratio for chromosomes XI, VII, and I ( Figure 2B and Table S3). In contrast, Fv999-Sk K patterns of CAPS markers on chromosome V were inherited more frequently than expected ($79.6%; Figure 2B and Table S3). The only exception was CAPS-5, which is located 20 kb from a telomere (Figure 2A and Table  3). Overall, these results mark the first independent confirmation of Xu and Leslie's (1996) original findings on the biased recovery of molecular markers linked to Sk K in Sk K · Sk S crosses.
Genome sequencing refines Sk to a 102 kb region of chromosome V To delineate the position of Sk K on chromosome V, we narrowed our focus to progeny JP98.75-Sk K , JP98.111-Sk K , and JP98.118-Sk K . Even though CAPS marker analysis failed to identify an Fv999-Sk K region of chromosome V that was inherited by all three of these progeny ( Figure  2C, black shaded regions), all were determined to carry the Sk K allele by phenotypic analysis (Figure 3 and data not shown). Therefore, their genomes were sequenced (Table 2) and their SNPs along chromosome V were examined. This analysis allowed us to identify two additional recombination events. Both of these recombination events occurred between CAPS-3 and CAPS-4 in the ascus that produced progeny JP98.111-Sk K (Figure 2D), explaining why they were not identified by our CAPS marker analysis. More importantly, the SNP profiles of all three progeny require that Sk K be located between positions 665,669 and 767,411 ( Figure 2D, red lines) with respect to the reference sequence (GenBank, CM000582.1). We refer to this 102-kb contiguous sequence as the Sk region (GenBank, KU963213).

Sk K strains carry a unique gene in the Sk region
The Sk region spans 102,256 bases in Fv999-Sk K , but only 101,743 bases in Fv149-Sk S . A Clustal W-based alignment (Thompson et al. 1994) of the sequences from both strains is 102,557 bases long, with 301 gaps in the Fv999-Sk K sequence and 814 gaps in the Fv149-Sk S sequence (Figure 4). To identify specific regions in Sk K and Sk S strains with high levels of SNPs or gaps, the number of SNPs or gaps were calculated for each Figure 3 Representative images of spore killing and resistance to spore killing in F. verticillioides asci. (A) Asci from a cross of progeny JP98.111-Sk K with Fv149-Sk S . Only four ascospores are detected in most of the mature asci, confirming that JP98.111-Sk K carries the Sk K allele. The asterisk denotes an ascus where there appears to be five ascospores. Two possible explanations are that this 5th ascospore escaped killing, or it strayed from another ascus during the dissection and imaging process. (B) Asci from a cross of progeny JP98.111-Sk K with BM43.31-Sk K . Eight ascospores can be detected in most of the mature asci, again confirming that progeny JP98.111-Sk K carries the Sk K allele.
100-base window across the alignment. A qualitative analysis of these results finds at least three hyper-variable intervals (Figure 4). One of the short hyper-variable intervals spans genes FVEG_03180 to FVEG_03175 , while another spans genes FVEG_03174 and FVEG_03173 (Figure 4). The longest hyper-variable interval is 14 kb long, spans gene FVEG_03167 , and extends to the right border of the Sk region (Figure 4).
The Sk region in strain Fv149-Sk S includes 41 putative protein-coding genes ( Figure 4 and Table 4). The Sk region in Fv999-Sk K carries the same 41 genes as well as one additional protein-coding gene, SKC1 (for Spore Killer Candidate 1), which is located between FVEG_03165 and FVEG_03164 (Figure 4). A BLASTP search (Altschul et al. 1997) of the NCBI nonredundant protein database with the predicted 70 amino acid sequence of Skc1 identified hypothetical proteins with Expect values ranging from 2 · 10 249 to 8 · 10 204 from various formae speciales of the F. oxysporum species complex ( Figure S1). No significant hits were identified in other species. A search of the NCBI Conserved Domain Database (Marchler-Bauer et al. 2010) failed to identify a domain within Skc1.
To confirm that SKC1 is absent in Fv149-Sk S , the sequences corresponding to the FVEG_03165 -FVEG_03164 intergenic region from Fv999-Sk K and Fv149-Sk S were aligned and examined. This intergenic region consists of 1234 and 755 bases in Fv999-Sk K and Fv149-Sk S respectively ( Figure S2). The alignment revealed that this difference in length was due to the presence of SKC1 in Fv999-Sk K and its absence from Fv149-Sk S ( Figure S2). To confirm that the absence of SKC1 in Fv149-Sk S was not due to an error in the reference genome sequence, we examined the lengths of the FVEG_03165 -FVEG_03164 intergenic region in our laboratory stocks of Fv999-SkK and Fv149-SkS by PCR. The PCR product lengths were consistent with the presence of SKC1 in Fv999-Sk K and its absence from Fv149-Sk S (Figure 5).
Clustal W alignments were used to investigate the level of polymorphism between proteins encoded by the Sk region in Fv999-Sk K and Fv149-Sk S strains (Table 4). Analysis of these pairwise alignments revealed that Fveg_15199 is the most polymorphic protein of the region. Only 80% of the 133 amino acids of Fveg_15199 are identical between the two strains. This is a remarkably high level of polymorphism given that 34 of 41 proteins within the Sk region of Fv999-Sk K and Fv149-Sk S are .94% identical, and 40 of 41 are .88% identical (Table 4). A BLAST search of the NCBI nonredundant protein database identified homologs of Fveg_15199 in several Fusarium, Neonectria, Scedosporium, and Aspergillus species ( Figure S3). A search of the NCBI Conserved Domain Database failed to identify a domain within Fveg_15199.
To shed light upon the transcriptional profile of protein-coding genes within the Sk region, we analyzed five F. verticillioides RNAseq datasets from NCBI's Sequence Read Archive. These datasets were for five time points following the induction of sexual development in a cross of Fv999-Sk K · Fv149-Sk S (Sikhakolli et al. 2012). Of the genes within the Sk region, FVEG_03171 exhibited the greatest increase in expression from the first time point to the last time point (Table 5, 2.95 and 94.58 RPKK), suggesting that this gene may have an important role during later stages of sexual development. The corresponding protein, Fveg_03171 , includes a WSC-domain (pfam01822), which has been linked to carbohydrate-binding, cell wall integrity, and stress response. The gene FVEG_03194 exhibited the greatest fold change when comparing differences between all five time points (.1050·) ( Table 5). The Fveg_03194 protein does not contain recognized domains, although putative homologs can be found in Fusarium, Nectria, Acremonium, and Trichoderma fungi ( Figure S4). Interestingly, SKC1 reached the highest level of expression of all examined genes 2 hr after fertilization  (Hall 1999) and aligned with Clustal W. Custom Perl scripts were used to examine base mismatches and gap positions, and to generate the diagram. The total number of mismatches (SNPs) between the two sequences was calculated for each 100-base nonoverlapping window of the alignment. The total number of gaps was also calculated for each 100-base nonoverlapping window. A gap position was not considered an SNP. By this definition, a window with 100 gaps cannot have SNPs. Black rectangles represent the coding regions of predicted genes. Gene names have been abbreviated according to their identification tags in the Fv149-Sk S annotation. For example: FVEG_03199 was shortened to 3199. A striking feature of the alignment is the presence of three hyper-variable intervals; spanning FVEG_03180 to FVEG_03175 , FVEG_03174 to FVEG_03173 , and FVEG_03167 to the right border of the Sk region. A second striking feature is the presence of a unique gene, SKC1, in Fv999-Sk K only. (Table 5). SKC1 transcript levels were also second only to FVEG_03197 at the last time point (144 hr, Table 5). Fveg_03197 appears to be from a family of proteins widely conserved among bacteria and eukaryotes ( Figure S5), but a function for it or its homologs is unknown.

DISCUSSION
The segregation of alternate alleles into separate gametes during meiosis encourages genetic conflict (Burt and Trivers 2008). Evidence for this is found in meiotic drive, which occurs when an allele is transmitted through meiosis in a biased manner. Meiotic drive elements are found in a diverse range of fungi, where they achieve biased transmission through sexual reproduction by killing meiospores carrying an alter-nate allele (Turner and Perkins 1979;Kathariou and Spieth 1982;Raju 1994;Dalstra et al. 2003;Grognet et al. 2014). Therefore, fungal meiotic drive elements are often referred to as Spore killers. The molecular basis of meiotic drive by most Spore killers is unknown.
The existence of meiotic drive by spore killing in Fusarium was first recognized in 1982. More than three decades later, we still do not understand the molecular mechanism that mediates this process. This is similar to the situation in Neurospora, where three distinct Spore killers, namely Sk-1, Sk-2, and Sk-3, were identified in 1979 (Turner and Perkins 1979). A breakthrough in Sk-2 research was recently made by refining the location of an Sk-2 resistance gene to a 52 kb sequence of DNA (Hammond et al. 2012a). This provided the necessary foundation to clone and characterize this resistance gene (Hammond et al. 2012b), which in turn allowed for the isolation of killer-less Sk-2 mutants (Harvey et al. 2014). Here, we have taken a similar approach toward identifying the genetic basis of Sk K in F. verticillioides by refining the position of Sk K to a 102 kb region of chromosome V.
As with previous work on Sk K (Xu and Leslie 1996), we observed a biased transmission of Sk K -linked molecular markers ( Figure 2B). This result was expected since the driving ability of Sk K should also affect the transmission of alleles linked to Sk K . This phenomenon is referred to as genetic hitchhiking (Lyttle 1991). The closer a hitchhiker is to a meiotic driver, the more likely it is to be transmitted to the next generation through the sexual cycle. In our study, patterns of CAPS markers from the parent Fv999-Sk K were inherited at a higher frequency than could be attributed to chance alone for five of six markers on chromosome V, and the transmission bias generally decreased with increasing distance from the Sk region ( Figure 2B). Surprisingly, CAPS-1 and CAPS-2 demonstrated higher levels of hitchhiking than CAPS-4 ( Figure 2B), despite the latter being closer to the Sk region. One explanation for this could be the existence of a recombination hotspot between CAPS-4 and the Sk region; however, we do not have additional data to support this hypothesis.
While a recombination hotspot could explain the relatively low level of hitchhiking by CAPS-4, the opposite phenomenon of recombinationsuppression is often associated with meiotic drive elements (reviewed by Lyttle 1991). For example, Neurospora Sk-2 requires specific alleles of at least two genes to mediate drive, a resistance gene called rsk and a killer gene called rfk (Hammond et al. 2012b;Harvey et al. 2014). Each gene is located on a different arm of chromosome III (Harvey et al. 2014). For Sk-2 to succeed as a meiotic driver, it is imperative that an Sk-2 ascospore inherit both rsk and rfk alleles because separation can lead to n Similarity (%) Figure 5 The Sk region in Fv149-Sk S is missing SKC1. To confirm that SKC1 is missing from the Sk region in Fv149-Sk S , genomic DNA was isolated from two liquid cultures of Fv149-Sk S and two liquid cultures of Fv999-Sk K . All four genomic DNA samples were used as templates in standard PCR reactions with the following oligonucleotide primers: 59 CGAATGACCTGGGGAGCCATAA 39 and 59 TCTCTCCAC CACCTCCATCAGC 39, which amplify the FVEG_03165 -FVEG_03164 intergenic regions from both Fv999-Sk K and Fv149-Sk S . PCR products were visualized by ethidium bromide staining after electrophoresis through a 1% agarose-TAE gel. We observed PCR products with lengths that were consistent with absence of SKC1 from Fv149-Sk S (287 bp) and presence of SKC1 in Fv999-Sk K (767 bp).
a self-killing genotype (Hammond et al. 2012b). This helps explain why Sk-2 is associated with a 30 cM "recombination-blocked" interval of chromosome III (Campbell and Turner 1987;Harvey et al. 2014). By suppressing recombination between rsk and rfk, Sk-2 can prevent these critical components of drive from separating during meiosis. With respect to F. verticillioides Sk K , recombination-suppression, if it exists at all, does not appear to be a significant phenomenon. For example, our CAPS marker analysis identified 12 recombination events between CAPS-3 and CAPS-4 (n = 59, Table S3), and two recombination events between CAPS-2 and CAPS-3 (n = 59, Table S3). This relatively high number of recombination events near the Sk region argues against an Sk-2-like recombination block for Sk K . If F. verticillioides Sk K requires multiple genes to function as a meiotic drive element, all should be found within the 102 kb Sk region defined by this study. For example, SNP-profiling of progeny JP98.75-Sk K , JP98.111-Sk K , and JP98.118-Sk K indicates that only this region of chromosome V is common between these three progeny and their Fv999-Sk K parent ( Figure 2D and Figure 4). Assuming killing is mediated by two or more distinct genes within the Sk region, the close proximity of these genes may negate the requirement for a Neurospora Sk-2-like recombination block. Alternatively, meiotic drive and spore killing may be mediated by a single gene. A precedent for this is seen in the het-s spore killing mechanism of Podospora anserina (Coustou et al. 1997;Saupe 2011). This system is controlled by two alternate alleles of a single gene, named het-s and het-S. The former allele, het-s, encodes the HET-s prion and is the meiotic driver, while the latter allele, het-S, encodes the HET-S protein. Interaction of the HET-s prion with the HET-S protein causes HET-S to relocate to cell membranes, resulting in cell death, presumably through loss of membrane integrity (Seuring et al. 2012). Ascospores with het-s genotypes escape cell death because n b The maximum RPKK in the table. c Fold change was calculated by dividing the maximum RPKK by the minimum RPKK. If the minimum RPKK for a gene was 0 it was arbitrarily assigned a minimum RPKK of 0.1 to approximate a fold change value. Unrounded maximum and minimum RPKK values (not shown) were used to calculate fold change values. they do not produce the HET-S protein, thus the het-s prion is not toxic to them. This is just one example of how meiotic drive and spore killing can be mediated by alternate alleles of a single gene in fungi. Grognet et al. (2014) have recently identified spok1 and spok2, two additional spore killing genes in P. anserina, both of which also appear to be single-gene-based meiotic drive systems.
In the current study, comparative sequence analysis revealed multiple differences between the 102 kb Sk region in Sk K and Sk S strains of F. verticillioides. Presumably, one or more of these differences is the genetic basis for the different Sk phenotypes exhibited by the strains. An analysis of SNPs and gaps identified at least three hyper-variable intervals within the Sk region (Figure 4). The interspersed nature of these intervals is somewhat surprising. The most striking difference, however, is the presence of the SKC1 gene in Sk K strains. Differences between the hypervariable intervals and SKC1 in Sk K vs. Sk S strains raise at least two key questions: (1) are the hyper-variable intervals and/or SKC1 responsible for killing, and (2) do they correspond to the Sk K allele that was previously defined by phenotypic analysis (Kathariou and Spieth 1982). The high level of expression of SKC1 throughout sexual development is consistent with its involvement in sexual reproduction. However, in addition to SKC1 and the hyper-variable intervals, there are less dramatic sequence differences between genes and intergenic sequences of the Sk region in Sk K and Sk S strains. The current data do not rule out the possibility that one or more of these differences is responsible for the Sk K and Sk S phenotypes. As a result, our efforts are now focused on targeted deletions of SKC1 and other sequences within the Sk region to determine the genetic basis of the spore killing phenotype in F. verticillioides.