Original Research ARTICLE
Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity
- 1Department of Plant Sciences, University of California, Davis, Davis, CA, USA
- 2School of Life Sciences and Warwick Systems Biology Centre, University of Warwick, Coventry, UK
How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.
Assessing how genomic variation may play a role in assisting pathogens to adapt to the difficulties of their corresponding hosts has largely been driven by cross-species comparisons as illustrated by the efforts on generating 1000 or more genomes for entire Kingdoms (for example 1000 Fungal Genomes Project, 10,000 Microbial Genomes Project, 100 K Pathogen Genome Project). In fungi, published pathogen genomes include numerous potential plant and animal pathogen species such as, Aspergillus fumigatus (Nierman et al., 2005), Cryptococcus neoformans (Loftus et al., 2005), Pythium ultimum (Lévesque et al., 2010), Sclerotinia sclerotiorum (Amselem et al., 2011), Taphrina deformans (Cissé et al., 2013), and Botrytis cinerea (Amselem et al., 2011; Staats and van Kan, 2012; Blanco-Ulate et al., 2013). This cross-species comparative genomics has identified some key components of host/pathogen interactions such as plant pathogens typically having more carbohydrate metabolic genes than animal pathogens (Porcel et al., 2014).
Cross-species comparisons are largely constrained to look at genome changes that are old and likely fixed within the species being analyzed (Hudson et al., 1987). In contrast, host/pathogen interactions are continuously changing and evolving. To identify these younger changes requires whole-genome re-sequencing of numerous individuals within a host or pathogen species to find the standing genetic variation that could be contributing to variation within current host-pathogen interactions. While many re-sequencing projects are under way such as the 20,000 Global pneumococcal project and several Saccharomyces cerevisiae re-sequencing projects (Wilkening et al., 2013), currently published sequence is limited to the host of any host/pathogen interaction, i.e., from Humans (Altshuler et al., 2010) and Arabidopsis (Long et al., 2013). Available datasets investigating genetic diversity in pathogens have tended to follow gene specific approaches (Baltrus et al., 2011; Guo et al., 2014) rather than unbiased whole genome investigations. Similarly, these approaches have focused solely on nuclear genes with minimal assessment of mitochondrial genome variation. A number of studies are beginning to show that genetic variation within the organellar genomes of diverse species can epistatically interact with the variation in the nuclear genome to control adaptive traits suggesting that it is necessary to consider that organelle variation may alter host-pathogen interactions (Etterson et al., 2007; Tang et al., 2007, 2013; Wolf, 2009; Dowling et al., 2010; Tan et al., 2012; Joseph et al., 2013a,b, 2015). Thus, there is restricted genomic information available to aid in the identification of causal polymorphisms in the pathogens controlling the host-pathogen interaction. This indicates a stark need for studies providing detailed whole genome measurements of genetic diversity in pathogens for both the nuclear and mitochondrial genomes.
B. cinerea is a necrotrophic fungal plant pathogen that has been a focus of gene variation studies over the past several decades. It has a genome of approximately 41–42 Mbp spread across 16 chromosomes (Shirane et al., 1989; Amselem et al., 2011; Staats and van Kan, 2012; Blanco-Ulate et al., 2013). This fungus has the capacity to infect and cause disease in living tissue of diverse plants ranging from numerous dicots to gymnosperms and bryophytes (Darvas et al., 1978; Coley-Smith et al., 1980; Lorbeer, 1980; Ponce de León et al., 2007; Williamson et al., 2007; Ponce de León et al., 2012). This host range is similar to other fungi such as Sclerotinia scleortiorum, Verticillium dahlia, and Rhizoctonia solani but there is not a common understanding of the molecular basis of this host range or how this may impact the genome.
B. cinerea exists in diverse environmental conditions in an array of developmental forms such as mycelia, micro- and macro-conidia, chlamydospores, sclerotia, apothecia, and ascospores (Coley-Smith et al., 1980; Lorbeer, 1980). Sclerotia provide B. cinerea with the ability to survive within a soil reservoir for many years. The diverse array of sporulation forms enables numerous dispersal avenues. In addition to a wide range of life-styles, this pathogen is also striking in the lack of large-effect resistance loci found within tested plant hosts. This suite of characters generates a pathogen which causes endemic crop losses and has been highly recalcitrant to genetic methods of control and chemical control requires complex interchanging of fungicides to prevent buildup of resistant genotypes. While previous work has shown extensive genomic variation between two isolates, there is little understanding of the genomic frequency of variation within a collection of isolates across the species (Staats and van Kan, 2012). Thus, a genome wide survey of genetic diversity in this pathogen may help to better devise appropriate control methods by understanding its life style and host range.
Recent genetic studies have shown that B. cinerea contains considerable genetic and phenotypic diversity and populations rapidly reshuffle genetic material (Buttner et al., 1994; Giraud et al., 1997, 1999; Fournier et al., 2002, 2013; Munoz et al., 2002; Kliebenstein et al., 2005; Rowe and Kliebenstein, 2007, 2008; Fournier and Giraud, 2008; Rowe et al., 2010; Amselem et al., 2011; Estavillo et al., 2011; Aguileta et al., 2012; Schumacher et al., 2012). Some individual gene studies and marker based studies using microsatellites have identified variation that may contribute to the extensive host range (Giraud et al., 1997, 1999; Fournier et al., 2002; Munoz et al., 2002; Rowe and Kliebenstein, 2007; Fournier and Giraud, 2008). Other studies have however found no significant signature for host specificity using a genome wide survey with AFLPs (Ma and Michailides, 2005). Similar contradictions were present across the same studies about the role of biogeography in structuring the measured diversity. This analysis has been complicated by the use of dominant markers and microsatellites that have uncertainty in ascertaining identity by descent and suggests that genomic resequencing may be a better approach to survey the genetic diversity in this species (Simonsen et al., 1995; Kruglyak et al., 1998). Similarly both approaches only sample a small fraction of the genome and are unable to give a full genomic view of recombination or diversity. As such, by conducting a full genomic sequencing of a collection of isolates, we are working to provide a full genomic view of diversity and recombination within the species.
The presence of recombination in combination with a large pool of standing genetic variation may extend the genetic potential available to B. cinerea and contribute to the broad host range of B. cinerea the species. This recombination may arise from two potential sources within B. cinerea. First, shuffling may be via meiotic recombination which is supported by the identification that B. cinerea has a single mating type gene that has at least two alleles, and a sexual cycle requiring a partner of opposite mating type (Faretra et al., 1988; Faretra and Pollastro, 1996; Amselem et al., 2011). While the sexual cycle has been observed in the laboratory, no reported sexual morphs have been found within the field, potentially due to its small size leading to an ascertainment bias against its identification (Alfonso et al., 2000). Additionally, wild-collected isolates of B. cinerea have shown no heterozygosity at any measured SNPs suggesting that the field isolates are all haploid and that any diploid form must be rapidly resolved to the haploid if there is field meiosis (Rowe and Kliebenstein, 2007). Another possibility is that vegetative fusion may play a role in allowing genetic shuffling in the field via mitotic recombination (Glass et al., 2004; Beever and Weeds, 2007). Thus, while existing data shows that recombination contributes to the standing genetic diversity of B. cinerea it remains to be clarified if this recombination arises from sexual or parasexual recombination or to what extent recombination occurs. A genomic survey of breakpoint identification will provide insights into the level of recombination within B. cinerea. Further, an assessment of diversity at vegetative incompatibility loci may allow an insight into the potential role of parasexuality in this recombination.
To measure the genomic diversity found in the nuclear and mitochondrial genomes of B. cinerea and to begin mapping variation in recombination across the genome, we conducted a pilot study where we resequenced the genome of 13 isolates (Table S1). These 13 isolates were chosen to provide a representative sample of the breadth of diversity in a larger collection of isolates as assessed by the genetic variation at two cell-wall-degrading polygalacturonase loci, BcPG1 and BcPG2 (Rowe and Kliebenstein, 2007). This set of isolates provides us the ability to conduct a pilot study to estimate the level of recombination across the genome within the species and potentially to assess what the source of this recombination may be. In addition, we re-sequenced the genome of one of these isolates using four replicated cultures that had been separately maintained for a decade in two labs, to assess the potential rate of spontaneous mutation accumulation under standard lab practices. Our analysis of major effect polymorphisms within these isolates found 13 major effect polymorphism clusters with an enrichment in loci potentially controlling vegetative incompatibility. Sequencing these genomes focuses on the goal of determining the level of standing genetic variation in B. cinerea and how frequently it is shuffled by recombination to initially assess if this genomic architecture could in part contribute to the species broad host range.
Materials and Methods
Isolates were grown on standard potato dextrose agar (PDA; Difco, Becton, Dickson and Company, East Rutherford, NY) with cellophane on the surface to prevent hyphal growth into the media. Information on isolate collection and sources is presented in Table S1. Young hyphae were collected from the leading edge of the mycelial mat using a sterilized spatula, and ground using stainless steel ball-bearings in P1 buffer from a Qiagen Plant DNeasy kit (Venlo, Netherlands). High quality DNA was obtained by following the manufacturer's recommendations for fungal DNA extraction. Sequencing libraries were constructed using a unique in-line 5′ bp barcode and pooled in groups of eight isolates. Pooled and multiplexed DNA libraries were sequenced on a total of four lanes using either an Illumina GAIIx or HiSeq machine. Post-sequencing, DNA libraries were de-multiplexed by their unique barcode for downstream analysis.
Sequence Clean-up, Alignment, and Variant Calling
Raw sequence data was cleaned-up (quality threshold of 30) using Fastq-mcf from the ea-utils package (http://code.google.com/p/ea-utils), and 9 bp were trimmed from the start of each read. Sequences were aligned using STAMPY (Lunter and Goodson, 2011), with default parameters, to the T4 reference sequence (Staats and van Kan, 2012). Supercontigs were placed on 16 chromosomes for plotting purposes only by alignment to Sclerotinia sclerotiorum using PAGIT (Swain et al., 2012). Alignment and SNP calling of mitochondrial DNA was implemented as described for genomic DNA, except the reference sequence from BO5.10 was used (Staats and van Kan, 2012).
Samtools (Li et al., 2009) and Picard (http://sourceforge.net/projects/picard/) were used to prepare the alignments for determining polymorphisms using GATK's UnifiedGenotyper (Depristo et al., 2011). GATK and Picard software was used to remove duplicate sequence reads, recalibrate base quality scores, and locally realign regions around indels. All polymorphisms in both BO5.10 (82,212 bp) and the polygalacturonase genes (5492 bp × 13) were visually checked using SAVANT (Fiume et al., 2010).
Initially, the first round vcf file produced by GATK UnifiedGenotyper was filtered by removing anything with < 6X coverage and quality threshold of 30, which gave zero false positive SNPs when compared to known sequence of the BcPG genes (Rowe and Kliebenstein, 2007). However, there were ~25% false negative SNP calls for the highly diverse genes BcPG1 and BcPG3. After additional filtering by mapping quality (MQ) of < 85, the majority of SNPs were called but three false positives were found. These were removed by a second Indel realignment. However, three false negatives were still found in BcPG1, in a difficult to call region of the Grape isolate, caused by clustering of four polymorphisms (three SNPs and an Indel) within a 16 base region. This was out of a total of 62 polymorphisms (GenBank: EF195807.1) and only when Indels were included. Similar polymorphisms were also not called in the isolates GLO1, Acacia and A517. Calling all polymorphisms as haploid (ploidy 1 setting in UnifiedGenotyper) removed these false negatives and all polymorphisms in the BcPG genes were called without any false positives.
The majority of initial heterozygous calls were removed during filtering. Those that remained, upon visual inspection, were determined as ‘hard to validate’ for example due to location of more than three SNPs and an indel within 10 bp region. As such, and due to the improved call rate, polymorphism calling by UnifiedGenotyper was carried out with the option for haploid data.
Variant annotation and effect prediction was performed using SnpEff version 3.4 (Cingolani et al., 2012) based on the T4 gene models for genomic and BO5.10 for mitochondria DNA (http://www.broadinstitute.org, B. cinerea BO5.10 mito; Staats and van Kan, 2012).
Consensus sequences were created and converted to fasta format from called SNPs using GATK before BIONeighbor-joining radial phylogenetic trees were determined based on pairwise SNP differences using SEAVIEW (Gouy et al., 2010).
The four independently propagated Grape isolates were aligned to a consensus extracted via GATK from the Grape genomic and mitochondrial sequence. Although ~5000 variants were originally called, upon closer inspection these were revealed to be due to low coverage regions. When the coverage cut-off was reduced to three reads the remaining variants could be visually checked using SAVANT and were found to be caused by regions of no coverage where the program reverted to the assumed T4 sequence, as such no genuine variants were detected.
Diversity and Clustering Analysis
Non-synonymous and synonymous polymorphisms were extracted from the diversity data generated by GATK UnifiedGenotyper as vcfs using SnpEff, followed by conversion to fasta consensus sequence for each contig using GATK. Non-synonymous and synonymous diversity was calculated as S (segregating sites) per contig using 10 kb sliding windows and 200 bp steps in DnaSP v5 (Librado and Rozas, 2009).
Major effect polymorphisms were extracted from the diversity data as vcfs using SnpEff and summed per gene. Sliding window plots were made of major effect polymorphisms per 10 gene window. A binary version of this data, where all polymorphisms per gene were altered to one polymorphism per gene was also created to visualize areas of interest more easily and to minimize the impact of individual loci with multiple major effect polymorphisms.
Proteins with non-synonymous SNPs and the translated T4 reference were blasted using Blast2GOPro from the Blast2go suite (Conesa et al., 2005) to obtain GO annotations. The extracted GO terms from the two datasets were analyzed using WEGO (Ye et al., 2006), which uses the Pearson Chi-sqaure test to show significant relations between two input datasets.
Vegetative Incompatibility Assays
All isolates were grown from spores stored in 10% glycerol solution at –80°C. Stock solution was diluted in filter-sterilized ¼strength grape juice and plated on 39 g/L potato dextrose agar (PDA), and grown at 25°C, 12 h light/12 h dark. Plugs (~4 mm2) cut from the static region of the colonized PDA were transferred to new 15 cm PDA petri dishes. One plug was transferred per quadrant, for four paired isolate interactions per plate.
After 3–5 days of growth at 25°C, 12 h light/12 h dark, heterokaryon incompatibility was assessed as the formation of a melanic barrage (demarcation zone) at the mycelial interface of two genetically distinct isolates, due to accumulated dying cells. Interactions visible as growth arrest or faint, partial pigmentation were scored as a weak incompatibility phenotype. We examined 3–6 interactions per genotype pair, and reported the average observed phenotype.
Recombination Breakpoint Analysis
RDP4 (Martin et al., 2010) was used to detect recombination in the 13 isolates plus the T4 reference. Default parameters were used for the seven programs RDP (Martin and Rybicki, 2000), GENECONV (Padidam et al., 1999), Bootscan (Martin et al., 2005), Chimera (Posada and Crandall, 2001), SiScan (Gibbs et al., 2000), MaxChi (Smith, 1992), and 3Seq (Gibbs et al., 2000). Recombination was considered to be real if detected by four or more of these programs using a Bonferroni corrected P-value cut off of 0.05.
Genome Sequencing and Polymorphism Control
To survey genomic diversity in B. cinerea, 11 isolates were sequenced with GAIIx 100-bp Illumina paired-end reads to obtain a mean 45-fold coverage (31–56X) and five isolates were sequenced by Illumina HiSeq 2000 100-bp paired-end reads to obtain a mean 90-fold coverage (75–133X). Information on isolate collection and sources is presented in Table S1. Reference based assemblies to call polymorphisms were by alignment to the previously sequenced B. cinerea isolate T4 for nuclear genomic and BO5.10 for mitochondrial genomic analyses (http://www.broadinstitute.org; B. cinerea T4; http://www.broadinstitute.org; B. cinerea B05.10; Staats and van Kan, 2012). This reference was selected, over the alternative previously sequenced reference BO5.10, as it had fewer and longer super-scaffolds (56 and N50 of 1.71 Mbp vs. 82 with N50 of 970 kbp; Staats and van Kan, 2012). Additional variant calling filters were applied to maximize the quality of the genome sequences generated. We optimized the filters to maximize polymorphism accuracy by empirical comparison to the previously sequenced cell-wall-degrading polygalacturonase genes, BcPG1-3, that had been sequenced in the same isolates (Rowe and Kliebenstein, 2007). We also utilized the empirical information available from alignment and polymorphism calling when comparing the reference sample BO5.10 to our BO5.10 sample. Additionally, all sequences were visually checked using Savant (Fiume et al., 2010). Using this approach our final settings were set to produce no false positives or false negatives amongst the BcPG1-3 sequences. Analysis of the mitochondrial genome of our BO5.10 sequence against the published BO5.10 mitochondrial sequence found one false positive SNP and three false negatives using these settings across 153,000 bp (which were manually confirmed). All four of these polymorphisms were unique to the published BO5.10 genome sequence and not present in any other isolates suggesting that they are likely low coverage errors. Thus, when combining all of the empirically tested genes, our false positive rate has an upper bound of one per 153 kb while our false negative rates have a maximal bound of three per 153 kb. These error rates are acceptably low. All sequence polymorphism information is available in a vcf format at doi: 10.5061/dryad.4vd70.
Clonal Propagation and Spontaneous Mutation Accumulation
Previous studies have suggested that the study of natural variation in virulence amongst pathogens may be influenced by the potential accumulation of spontaneous mutations during lab cultivation as clonal propagates have been noted to vary for virulence (Jeon et al., 2013). However, recent work in several species has shown that clonal propagates are not showing elevated mutagenesis (Kohn et al., 2008; Park et al., 2010). To test if there was any evidence for clonal propagation in B. cinerea that could be linked to mutation accumulation, we resequenced the genomes of four clones that had been experimentally separated from one isolate (Grape) and independently propagated for a decade, in two separate labs, on different continents (US and UK). The genome resequenced samples were then realigned both to the T4 reference and also to the sequence of the progenitor Grape isolate as described for the other isolates. All potential variation between the four independent Grape genomes was manually investigated and found to be solely present in areas lacking sufficient coverage. Thus, we could not identify any de novo polymorphisms between the genomes of these four independently propagated Grape isolate clones after controlling for lower coverage thresholds. As such, the Grape strain of B. cinerea shows no evidence for the accumulation of spontaneous mutations over at least a decade of separate propagation under lab conditions. This suggests that B. cinerea has a similar mutation rate as Sclerotinia sclerotiorum and Magnaporthe oryzae which indicates that mutation rate is not the main driver in the level of genomic diversity we identify in our isolates as compared to these other fungi (Kohn et al., 2008; Park et al., 2010).
Elevated Nuclear Genomic Diversity
To identify genome-wide polymorphism levels in the nuclear genome, we aligned the genomes of the above isolates to the previously sequenced reference T4. This analysis identified 980,329 unique SNPs and 157,936 unique Indels across all isolates with between 144,109 and 403,608 polymorphisms per isolate (Table 1). There was no significant correlation between the level of sequence variation in an isolate and the depth of sequencing (Table 1, P = 0.09, R = 0.52, N = 12). Any appearance of a correlation was solely driven by the Grape isolate being both more sequenced and generally more diverged. Excluding this single outlier abolished any relationship between sequence variation and sequencing depth (Table 1, P = 0.99, R = 0.002, N =11). Aligning to the T4 reference also found extensive insertion and deletion polymorphisms. On average there were 12,560 deletions (7838–21,166) and 13,966 insertions (8492–25,137) per isolate in comparison to T4 (Table 1). Combining the polymorphisms (both SNP and indel) from all pairwise combinations of isolates yielded an average of at least one polymorphism between the isolates every 36 bp, giving an average density of 28 polymorphisms per kb within the population. Plotting the number of isolates that share a specific polymorphism in comparison to T4 showed that most alleles were of minor frequency with the major peak being a single isolate having each specific difference in comparison to T4 (Figure 1). These estimates for genomic diversity within B. cinerea are higher than that found so far in species with detailed whole genome resequencing projects such as Arabidopsis (Cao et al., 2011; Long et al., 2013), Humans (Altshuler et al., 2010), and Yeast (Schacherer et al., 2009).
Figure 1. Allele diversity in comparison to reference T4. The number of isolates that share a polymorphic allele in comparison to the reference T4 isolate was plotted against the total number of polymorphisms with that level of diversity.
To investigate the potential functional consequence of these polymorphisms, we focused on polymorphisms within coding regions of predicted genes. This found 98,400 missense (42.8%), 1255 nonsense (0.5%), and 130,295 (56.7%) silent polymorphisms. The missense/silent ratio of 0.76 is in line with that found using the three published B. cinerea genomes (Amselem et al., 2011; Staats and van Kan, 2012; Blanco-Ulate et al., 2013). Using SnpEff (Cingolani et al., 2012) we extracted polymorphisms predicted to have major effects (protein truncation, loss of function or triggering nonsense mediated decay) on gene function using the annotation from T4 (http://www.broadinstitute.org; B. cinerea T4). This found 5963 polymorphisms likely to have a high impact on the function of 1801 genes (out of a possible 10,401 genes annotated in the T4 reference). Thus, 17% of the annotated genes in B. cinerea have at least one major effect polymorphism (polymorphisms that lead to a predicted functional knockout) in this limited sample of isolates and the genetic diversity is likely to have significant phenotypic impact on the processes mediated by these genes.
Genome Wide Evidence of Recombinational Shuffling
To assess the genomic relationships amongst this collection of isolates, we used the genome wide diversity to plot the relationships of the isolates via a Neighbor-joining phylogenetic tree (Figure 2A). This was done using pairwise SNP differences to T4 to generate the alignments. The un-rooted radial phylogeny of the 17 isolates (including the reference T4) showed no statistically supported evidence of clustering indicating that there is no genetic isolation amongst the isolates. Of the isolates sampled, Katie Tomato showed the fewest and A517 the most SNP differences to the reference T4. The only genomes that were significantly associated were those from the repeated resequencing of independently propagated Grape isolate clones (Figure 2A).
Figure 2. Radial neighbor-joining phylogenetic trees of nuclear and mitochondrial diversity within B. cinerea. (A) Unrooted genomic phylogeny determined based on pairwise SNP differences in the alignments of 17 B. cinerea strains to the reference T4 (including T4 reference and three replicate lineages independently propagated from the Grape isolate). Branch lengths are proportional to the number of segregating sites that differentiate each pair of isolates. (B) Unrooted mitochondrial phylogeny determined based on pairwise SNP differences to B05.10 in the alignments of 17 B. cinerea strains (including BO5.10 reference and three replicate lineages independently propagated from the Grape isolate). Branch lengths are based on very few polymorphisms and are proportional to the number of segregating sites that differentiate each pair of isolates.
We also used the genomic variation to directly search for recombination breakpoints using the package RDP4 (Martin et al., 2010). Analyzing the data on a per contig basis allowed us to identify 7735 potential recombination breakpoints using a Bonferroni corrected P-value cut off of 0.05 (Table S2). This shows that within this collection of isolates, there is evidence for a recombination breakpoint every 5.4 kb. There was significant variation in local recombination rates with a low of one breakpoint per 9.7 kb on contig 11.0 and a high of a breakpoint every 0.8 kb on contig 15.3 (Table S2). This shows that there is significant evidence for recombination amongst isolates within this species. This analysis suggests extensive recombination amongst these B. cinerea isolates on a genome wide level and is in direct contrast to other published plant pathogens that show extensive clonality in their genome or phenotypic variation (Blandon-Diaz et al., 2012; Wicker et al., 2013).
Mitochondrial Genome Polymorphisms
Using the sequencing data, we next reconstructed the mitochondrial genomes in each of the isolates to compare the mitochondrial variation to the nuclear genome variation and obtained a mean 205-fold coverage (163–244X). Aligning the mitochondrial genome of the 13 isolates to the 82 kb of mitochondrial sequence from BO5.10 identified 53 SNPs, 52 insertions, and 32 deletions (Table 2). The BO5.10 mitochondrial genome was used because the T4 mitochondrial sequence was not published. This showed that mitochondrial genomic variation was approximately 20-fold lower than for the nuclear genome with an average polymorphism every 600 bp from the combined isolates. These SNPs included six missense and five silent polymorphisms within ORFs. Using SnpEff to identify major effect polymorphisms, nine mitochondrial polymorphisms were predicted to have large functional consequences in six genes (Table 3). These six genes include three genes within the NADH dehydrogenase complex, an ATP synthase subunit, LAGLIDADG endonuclease and NUMOD endonuclease. The LAGLIDADG and NUMOD endonuclease genes are self-mobile inteins/homing endonucleases that act like transposons suggesting that their major effect polymorphisms may affect the movement of these genes within the species (Stoddard, 2005). A larger collection of isolates will be required to test if there are different copy numbers and/or mobility of these mitochondrial transposons within the species. Thus, while the mitochondrial genome has less genetic diversity than the nuclear genome, there is still the potential for this genetic diversity to influence phenotypic diversity in the pathogen.
Within Fungi, mitochondria can either show uniparental inheritance or biparental inheritance with fusion and recombination (Basse, 2010). Phylogenetic analysis of the mitochondrial diversity suggested four clusters based on a total of 90 SNPs found amongst the isolates (Figure 2B). Although there was more apparent clustering in the mitochondrial phylogeny than found using the nuclear genome, it was not sufficient to support either of the potential patterns of mitochondrial inheritance. Interestingly, and without correlation to sequencing depth (Table 2), isolate A517 shows elevated mitochondrial variation in comparison to the other 12 isolates and accounts for 27% of the total mitochondrial variation, which is reflected in the phylogeny. This elevated mitochondrial diversity rate in A517 is significantly greater than the nuclear genome diversification of this isolate. It remains to be seen if there are any consequences of this elevated mitochondrial diversity. All of the independent propagates of Grape clustered given their lack of diversity. Thus, it is possible to identify significant levels of mitochondrial genomic diversity in B. cinerea and this may not reflect the diversity in the genome. Future studies are required to test if this mitochondrial diversity impacts the host/pathogen interaction.
Genomic Distribution of Polymorphisms
Given the high level of genomic diversity in B. cinerea, we sought to determine if this diversity was randomly distributed across the genome or showed evidence of clustering. The chromosomes of B. cinerea have not been delineated, as such the program PAGIT (Swain et al., 2012) was used to superficially overlay all contigs on chromosomes by comparison to Sclerotinia sclerotiorum (http://www.broadinstitute.org/sclerotinia_sclerotiorum/maps). This fungus is highly co-linear to B. cinerea with an identical chromosome number (Shirane et al., 1989; Amselem et al., 2011), and has an almost complete chromosome map. This created a visualization framework by placing the contigs on 15 chromosomes, with all unplaceable contigs being assigned into artificial chromosome 16 (as they were in the S. sclerotiorum used). The distribution of polymorphisms amongst contigs, placed on these superficial chromosomes, displayed regions or peaks of high density polymorphism, across all chromosomes (Figure 3A and Table S2). Change rates varied greatly per isolate and per contig, with the highest being one polymorphism every 40 bp in contig 15.3 of Grape and the lowest one change every 27,839 bp in contig 4.1 of Katie Tomato (Table S2).
Figure 3. Genomic distribution of polymorphism amongst the 14 isolates. (A) A sliding window analysis of genomic polymorphisms per 10 kb across the genome is shown. (B) Shown is a 10 gene sliding window analysis of the number of genes with major effect polymorphisms based on comparison to the isolate T4. (C) Again a 10 gene sliding window analysis of the number of genes with major effect polymorphisms based on comparison to the isolate T4, but with polymorphisms per gene given a binary value. This enabled clusters of seven or more genes in 10 gene windows affected by major effect polymorphisms to be visually determined more easily. Thousands of permutations of the major effect polymorphisms established a significance threshold of six genes per window, which is drawn for reference. No one isolate is key cause of these clusters.
To further dissect differences in the distribution of polymorphisms likely to affect function, the locations of major effect (frame shift, start lost, stop gained polymorphisms) were plotted within a 10-gene sliding window across the genome (Figure 3B). Plotting the major effect diversity as binary values rather than total number of major effect polymorphisms highlighted clearer clusters for gene dissection. This identified 13 windows where there were more major effect polymorphisms clustered than would be expected by random chance given this collection of isolates (Figure 3C, significance tested by a permutation test). Each of these clusters contains seven or more genes with high impact polymorphisms within a 10-gene window (Figure 3C).
An analysis of these polymorphisms showed that they were mostly moderate frequency polymorphisms (i.e., each allele occurring in more than one isolate) and not singleton polymorphisms where the alternate allele exists in only one particular isolate (Table S3). Gene families (determined from the T4 annotation) over-represented in these 13 major effect polymorphism clusters, showed statistically significant enrichments for heterokaryon incompatibility genes (5 of 67; P < 0.001 via χ2) and Nacht domain protein encoding genes (4 of 32; P < 0.001 via χ2) which have previously been associated with vegetative incompatibility in fungi (Koonin and Aravind, 2000; Paoletti et al., 2007; Zhao et al., 2015). Together, these five heterokaryon incompatibility genes and four Nacht domain protein encoding genes were present in six of the 13 major effect clusters. To test if the occurrence of these genes in the clusters was due to the particular gene families having higher diversity rates than the average gene, the rate of non-synonymous polymorphisms in all the vegetative incompatibility and Nacht domain genes was compared against the genome average. Neither gene family showed an elevated level of non-synonymous polymorphisms in comparison to the genomic average (8 of 32 NACHT genes and 11 of 67 HET proteins contained non-synonymous SNPs). Thus, the enriched presence of these genes within these major effect clusters is not being driven by elevated diversity of the gene family in relation to the rest of the genome (Table S4).
Five of the remaining seven clusters (that did not include heterokaryon incompatibility genes and or Nacht domain protein encoding genes) all contained one or more transporters of the Major Facilitator Superfamily (8 out of 287). Although this does not reflect enrichment of this family in the clusters, it suggests that transporters may be another target of genetic diversity. These transporters may play a role in controlling the pathogen's resistance to plant defense compounds (Stefanato et al., 2009). Alternatively, they may be required to allow the export of fungal phytotoxins that target the plant to enhance virulence (Stergiopoulos et al., 2002; Choquer et al., 2007). Future work with larger populations will be required to test the potential function of these transporters in any plant/pathogen interaction.
Biological Conformation of Heterokaryon Incompatible Loci Phenotype for Loci within Major Effect Polymorphism Clusters
To test if the major effect allelic diversity of the heterokaryon incompatible loci located within the major effect clusters could be controlling variation in vegetative incompatibility; the isolates were grouped based on their haplotypes at these loci (Figure 4 and replication information is provided in Table S6). We then conducted vegetative incompatibility assays between all pairs of isolates to test if the heterokaryon haplotype modulated vegetative fusion. In these tests, the incompatibility was significantly stronger when vegetative fusion was tested between isolates with different haplotypes than when they shared a haplotype (χ2 = 4.9, P = 0.02, N = 66, df = 2). All strong incompatibilities were found when vegetative fusion was tested between isolates that represent different haplotypes at the major effect clusters (Figure 4). In contrast, when vegetative fusion was tested within a haplotype group, there were no strong incompatibilities (Figure 4). This suggests that while we have likely identified some components controlling diversity in the vegetative fusion process, there are other vegetative incompatibility loci remaining to be identified potentially via sequencing a larger collection of isolates.
Figure 4. Pairwise Vegetative Incompatibility. All pairwise crosses were tested in at least triplicate and the incompatibility was scored as compatible (C), weakly incompatible (W), and incompatible (I) and results are shown above. χ2 was used to test if there was a difference in incompatibility when crossing within or between the Major Effect Haplotypes (χ2 = 4.9, P = 0.02, N = 66, df = 2).
Genomic Distribution of Non-synonymous and Synonymous Polymorphisms
To compare the distribution of potentially neutral and non-neutral polymorphisms the genomic diversity in non-synonymous and synonymous polymorphisms was plotted as the number of segregating sites (S) within 10 kb sliding windows (Figure 5). This identified a number of regions with an elevated rate of clustering of non-synonymous SNPs. Notably, non-synonymous SNPs are abundant in the regions of major effect polymorphism (Figures 3, 5). Highly polymorphic regions of non-synonymous and synonymous SNPs appeared to occur in the same 10 kb bins as each other when visualized in a mirrored plot (Figure 5). However, in the majority of these regions of clustered polymorphisms, the rate of non-synonymous SNPs is noticeably higher (Figures 5, 6A). In contrast, there were only three regions of clustered polymorphisms where the synonymous polymorphism level was higher than the non-synonymous rate (one on chromosome 1 and two on chromosome 13). In contrast to the regions of elevated diversity, synonymous SNPs were more abundant overall across the genome in the regions of moderate or lower polymorphism rates (Figure 6B). This suggests that these regions of high non-synonymous polymorphism and major effect polymorphisms are non-random. Fully testing this hypothesis and assessing the potential for diversifying selection will require the sequencing of more isolates to generate a larger sample size.
Figure 5. Genomic distribution of non-synonymous and synonymous SNPs. Polymorphisms were assigned as non-synonymous or synonymous amongst the 14 isolates (including T4) and then measured in a 10 gene sliding window as the number of segregating sites (S). The top mirrored plot shows the sliding window of non-synonymous polymorphisms while the bottom shows the corresponding synonymous analysis.
Figure 6. Relationship of Synonymous vs. Non-synonymous SNP frequencies. (A) An X/Y scatter plot showing the frequency of Synonymous vs. Non-synonymous SNPs per 10 kb bins. The x and y axis show the fraction of SNPs per kb for each window. No significant correlation was found using either pearson or spearman rank. (B) A zoomed in view of (A), focused on the regions up to four polymorphisms per kb per 10 kb bin.
To investigate non-synonymous SNPs further, as they are also likely to cause mutations which affect phenotype, we used Blast2GOPro from the Blast2go suite to test if there was enrichment for a particular function in genes with non-synonymous polymorphisms (Conesa et al., 2005). GO terms were extracted for both the translated T4 reference sequence and loci classified by presence or absence of non-synonymous SNPs. These two sets of GO terms were analyzed using WEGO (Ye et al., 2006), which uses the Pearson Chi-square test to show significant relations between two input datasets. GO enrichment analysis found that the proteins containing non-synonymous SNPs were involved in 46 of the ontologies that could be analyzed (Table S5 and Figure S1). Of interest amongst these (Table S5, highlighted in blue) were numerous non-synonymous SNP enrichments that were found in proteins associated with botcinin production. Botcinins are chemical toxins produced by B. cinerea that are virulence determinants (Dalmais et al., 2011) and as such are of interest in potentially determining isolate and crop specific interactions.
While previous work had shown significant genomic variation between two isolates, this had not assessed the potential frequency or level of diversity in the species (Staats and van Kan, 2012). Our genomic resequencing of 13 B. cinerea isolates revealed that the previous two isolates had only sampled a small range of the diversity found in this collection. Thus, we are still sampling new diversity at a high level (Amselem et al., 2011; Staats and van Kan, 2012; Blanco-Ulate et al., 2013). Importantly, there is a significant fraction of the polymorphisms that have moderate allele frequencies indicating that there is a large pool of standing variation. With the presence of extensive recombination in this species, this allows for this variation to be readily shuffled. Thus, we have only begun to sample the potential genetic diversity in this pathogen.
Even with this under-sampling, the level of diversity is higher than in any previously studied plant pathogen. For example, the sequencing of two Blumeria graminis ff. ssp. hordei isolates and comparison to the previously sequenced reference DH14 gave an approximate SNP change rate of one SNP per kb vs. 3–9 polymorphic SNPs per kb between individual isolates found within our study (Hacquard et al., 2013). Resequencing of Blumeria graminis ff. ssp tritici provided numbers of between 1.4 and 2.0 SNP per kb between genotypes (Wicker et al., 2013). This puts the estimates of genomic diversity in both plant pathogens lower than in B. cinerea. One potential reason for this level of genomic diversity within B. cinerea is that it may facilitate the broad host-range of the pathogen, including nearly every plant tested to date. Broader sampling of both B. cinerea as well as other broad and narrow host range plant pathogens is required to investigate if the level of genomic diversity may contribute to host-range diversity.
Mitochondrial Genomic Variation
In animals, mitochondrial DNA has been shown to accumulate polymorphisms more rapidly than genomic DNA (Birky and Walsh, 1992; Birky, 1995). In our B. cinerea analysis, this relationship was reversed, with the genomic DNA having a nearly 20-fold higher diversity rate than the mitochondria (Figure 2). A low level of mitochondrial diversity has also been found in other fungi (Zhan et al., 2003). For example, a study of four strains of Rhizophagus irregularis found zero SNPs (Formey et al., 2012) and a study of M. graminicola found a similar diversity rate as we detected in B. cinerea mitochondria (Torriani et al., 2008). Unfortunately, measures of nuclear genomic diversity were not available in these studies so it remains to be seen if most fungi have lower mitochondrial diversity than nuclear genome diversity. Even with this low level of diversity, there were major effect polymorphisms in specific genes within the mitochondria suggesting diversity could cause phenotypic consequences.
Observed Field Adaptation of B. cinerea is Likely via Selection of Standing Population Diversity Rather than Elevated Mutagenesis
Frequent studies focused on the agronomic control of B. cinerea infections in the field or within lab studies have suggested that the ability of B. cinerea to show resistance to fungicides or isolates to display variation might be caused by the accumulation of spontaneous mutations (Grindle, 1979; Kerssies et al., 1997; Yourman and Jeffers, 1999; Yourman et al., 2000, 2001). This mutation accumulation hypothesis has often been postulated to occur in fungi, viruses and prokaryotes (Drake et al., 1998; Lamb et al., 2008). However, recent work in in Sclerotinia sclerotiorum and Magnoparthe oryzae has shown that there does not appear to be an effect of lab culturing on mutation rate (Kohn et al., 2008; Park et al., 2010). Our genomic resequencing of four individual propagates detected no spontaneous polymorphisms in either the genomic or mitochondrial alignments. Each line was propagated on average at least four times per year but given the continuous formation of conidia in B. cinerea it is not possible to generate an accurate estimate of generation number. As such, it would appear that Sclerotinia sclerotiorum, Magnoparthe oryzae and B. cinerea do not have an elevated spontaneous mutation rate (Kohn et al., 2008; Park et al., 2010). This suggests that there is a large level of standing diversity within B. cinerea populations and that the observed field “evolution” against a fungicide or other agronomic control is likely a result of selection on massive standing diversity rather than de novo evolution.
Vegetative Incompatibility Genes were found in the Major Effect Polymorphism Clusters
Our analysis of major effect polymorphisms within these isolates found 13 major effect polymorphism clusters with an enrichment in loci potentially controlling vegetative incompatibility. Filamentous ascomycete fungi can fuse their hyphae, fuse their nuclei, recombine and resolve back into haploid isolates in a process called parasexuality. This process is prevented by incompatibility loci called heterokaryon incompatibility loci which stimulate programmed cell death when incompatible hyphae fuse (Glass et al., 2004). These het (heterokaryon) loci have previously been shown to be under selection for diversity in other fungal species (Hall et al., 2010; Zhao et al., 2015). Further, they can prevent transmission of mycoviruses by hyphal fusion. Loss-of-function mutations at het genes in N. crassa do not affect the vegetative growth phenotype of mutants but these mutants fail to distinguish self from non-self and will form vigorous heterokaryons with strains with which they were formerly incompatible (Shiu and Glass, 1999). Noticeably, all genes containing HET-E-1 domains that are present within major effect polymorphism clusters contained major effect polymorphisms, i.e., variants that should abolish their function. However, these polymorphisms at independent loci were shuffled amongst the isolates with different genes being polymorphic in different isolates.
Biological confirmation of the major effect polymorphism within these heterokaryon incompatible loci showed that we could predict the vegetative fusion interaction using the haplotype at these het loci. This observation is in agreement with previous reports that suggested a potential role for parasexuality in B. cinerea diversity (Beever and Weeds, 2007). What this means in regards to reproductive barriers between isolates and how this may affect the adaptation to new plant hosts should become clearer as more isolates are resequenced.
Our genomic resequencing has shown that the broad-host range plant pathogen B. cinerea contains a high level of genomic diversity. Further, there appears to be shuffling of this polymorphism within the nuclear genome as evidenced by the lack of clonality. That B. cinerea has a large pool of shuffling standing genetic diversity may enhance its ability to adapt to new hosts. This high level of standing variation suggests that there is significant utility in the creation of a large B. cinerea mapping population. Future efforts will focus on conducting genome-wide association mapping within the pathogen to determine how this genomic diversity leads to phenotypic variation in the host-pathogen interaction. As such it will be interesting to assess how larger diversity studies in B. cinerea identify the link between genotype and phenotype in altering the plant/pathogen interaction.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the NSF IOS grant 1021861 to DK, the NSF IOS grant 1339125 to DK, the USDA National Institute of Food and Agriculture, Hatch project number CA-D-PLS-7033-H to DK and by the Danish National Research Foundation (DNRF99) grant to DK.
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2015.00996
Aguileta, G., Lengelle, J., Chiapello, H., Giraud, T., Viaud, M., Fournier, E., et al. (2012). Genes under positive selection in a model plant pathogenic fungus, Botrytis. Infect. Genet. Evol. 12, 987–996. doi: 10.1016/j.meegid.2012.02.012
Alfonso, C., Raposo, R., and Melgarejo, P. (2000). Genetic diversity in Botrytis cinerea populations on vegetable crops in greenhouses in south-eastern Spain. Plant Pathol. 49, 243–251. doi: 10.1046/j.1365-3059.2000.00452.x
Altshuler, D., Durbin, R. M., Abecasis, G. R., Bentley, D. R., Chakravarti, A., Clark, A. G., et al. (2010). A map of human genome variation from population-scale sequencing. Nature 467, 1061–1073. doi: 10.1038/nature09534
Amselem, J., Cuomo, C. A., van Kan, J. A., Viaud, M., Benito, E. P., Couloux, A., et al. (2011). Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. Plos Genet. 7:e1002230. doi: 10.1371/journal.pgen.1002230
Baltrus, D. A., Nishimura, M. T., Romanchuk, A., Chang, J. H., Mukhtar, M. S., Cherkis, K., et al. (2011). Dynamic evolution of pathogenicity revealed by sequencing and comparative genomics of 19 Pseudomonas syringae isolates. Plos Pathog. 7:e1002132. doi: 10.1371/journal.ppat.1002132
Beever, R. E., and Weeds, P. L. (2007). “Taxonomy and genetic variation of Botrytis and Botryotinia,” in Botrytis: Biology, Pathology and Control, eds Y. Elad, B. Williamson, P. Tudzynski, and N. Delen (Heidelberg: Springer), 29–52.
Blanco-Ulate, B., Allen, G., Powell, A. L., and Cantu, D. (2013). Draft genome sequence of Botrytis cinerea BcDW1, inoculum for noble rot of grape berries. Genome Announc. 1:e00252-13. doi: 10.1128/genomeA.00252-13
Blandon-Diaz, J. U., Widmark, A. K., Hannukkala, A., Andersson, B., Hogberg, N., and Yuen, J. E. (2012). phenotypic variation within a clonal lineage of Phytophthora infestans infecting both tomato and potato in Nicaragua. Phytopathology 102, 323–330. doi: 10.1094/PHYTO-02-11-0033
Buttner, P., Koch, F., Voigt, K., Quidde, T., Risch, S., Blaich, R., et al. (1994). Variations in ploidy among isolates of Botrytis cinerea - implications for genetic and molecular analyses. Curr. Genet. 25, 445–450. doi: 10.1007/BF00351784
Cao, J., Schneeberger, K., Ossowski, S., Guenther, T., Bender, S., Fitz, J., et al. (2011). Whole-genome sequencing of multiple Arabidopsis thaliana populations. Nat. Genet. 43, 956–963. doi: 10.1038/ng.911
Choquer, M., Fournier, E., Kunz, C., Levis, C., Pradier, J. M., Simon, A., et al. (2007). Botrytis cinerea virulence factors: new insights into a necrotrophic and polyphageous pathogen. FEMS Microbiol. Lett. 277, 1–10. doi: 10.1111/j.1574-6968.2007.00930.x
Cingolani, P., Platts, A., Wang, L. L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w(1118); iso-2; iso-3. Fly (Austin) 6, 80–92. doi: 10.4161/fly.19695
Cissé, O. H., Almeida, J. M., Fonseca, A., Kumar, A. A., Salojarvi, J., Overmyer, K., et al. (2013). Genome sequencing of the plant pathogen Taphrina deformans, the causal agent of peach leaf curl. Mbio 4:e00055-13. doi: 10.1128/mBio.00055-13
Conesa, A., Gotz, S., Garcia-Gómez, J. M., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Dalmais, B., Schumacher, J., Moraga, J., le Pecheur, P., Tudzynski, B., Collado, I. G., et al. (2011). The Botrytis cinerea phytotoxin botcinic acid requires two polyketide synthases for production and has a redundant role in virulence with botrydial. Mol. Plant Pathol. 12, 564–579. doi: 10.1111/j.1364-3703.2010.00692.x
Darvas, J. M., Scott, D. B., and Kotzé, J. M. (1978). Fungi associated with Damping-off in coniferous seedlings in South African nurseries. South Afr. For. J. 104, 15–19. doi: 10.1080/00382167.1978.9629481
Depristo, M. A., Banks, E., Poplin, R., Garimella, K. V., Maguire, J. R., Hartl, C., et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498. doi: 10.1038/ng.806
Estavillo, G. M., Crisp, P. A., Pornsiriwong, W., Wirtz, M., Collinge, D., Carrie, C., et al. (2011). Evidence for a SAL1-PAP chloroplast retrograde pathway that functions in drought and high light signaling in Arabidopsis. Plant Cell 23, 3992–4012. doi: 10.1105/tpc.111.091033
Etterson, J. R., Keller, S. R., and Galloway, L. F. (2007). Epistatic and cytonuclear interactions govern outbreeding depression in the autotetraploid Campanulastrum americanum. Evolution 61, 2671–2683. doi: 10.1111/j.1558-5646.2007.00234.x
Faretra, F., and Pollastro, S. (1996). Genetic studies of the phytopathogenic fungus Botryotinia fuckeliana (Botrytis cinerea) by analysis of ordered tetrads. Mycol. Res. 100, 620–624. doi: 10.1016/S0953-7562(96)80018-9
Formey, D., Molés, M., Haouy, A., Savelli, B., Bouchez, O., Bécard, G., et al. (2012). Comparative analysis of mitochondrial genomes of Rhizophagus irregularis - syn. Glomus irregulare - reveals a polymorphism induced by variability generating elements. New Phytol. 196, 1217–1227. doi: 10.1111/j.1469-8137.2012.04283.x
Fournier, E., and Giraud, T. (2008). Sympatric genetic differentiation of a generalist pathogenic fungus, Botrytis cinerea, on two different host plants, grapevine and bramble. J. Evol. Biol. 21, 122–132. doi: 10.1111/j.1420-9101.2007.01462.x
Fournier, E., Giraud, T., Loiseau, A., Vautrin, D., Estoup, A., Solignac, M., et al. (2002). Characterization of nine polymorphic microsatellite loci in the fungus Botrytis cinerea (Ascomycota). Mol. Ecol. Notes 2, 253–255. doi: 10.1046/j.1471-8286.2002.00207.x
Fournier, E., Gladieux, P., and Giraud, T. (2013). The ‘Dr Jekyll and Mr Hyde fungus’: noble rot versus gray mold symptoms of Botrytis cinerea on grapes. Evol. Appl. 6, 960–969. doi: 10.1111/eva.12079
Gibbs, M. J., Armstrong, J. S., and Gibbs, A. J. (2000). Sister-Scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics 16, 573–582. doi: 10.1093/bioinformatics/16.7.573
Giraud, T., Fortini, D., Levis, C., Lamarque, C., Leroux, P., Lobuglio, K., et al. (1999). Two sibling species of the Botrytis cinerea complex, transposa and vacuma, are found in sympatry on numerous host plants. Phytopathology 89, 967–973. doi: 10.1094/PHYTO.19126.96.36.1997
Giraud, T., Fortini, D., Levis, C., Leroux, P., and Brygoo, Y. (1997). RFLP markers show genetic recombination in Botryotinia fuckeliana (Botrytis cinerea) and transposable elements reveal two sympatric species. Mol. Biol. Evol. 14, 1177–1185. doi: 10.1093/oxfordjournals.molbev.a025727
Gouy, M., Guindon, S., and Gascuel, O. (2010). SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol. 27, 221–224. doi: 10.1093/molbev/msp259
Guo, L, Han, L., Yang, L., Zeng, H., Fan, D., Zhu, Y., et al. (2014). Genome and transcriptome analysis of the fungal pathogen Fusarium oxysporum f. sp cubense causing banana vascular wilt disease. PLoS ONE 9:e95543. doi: 10.1371/journal.pone.0095543
Hacquard, S., Kracher, B., Maekawa, T., Vernaldi, S., Schulze-Lefert, P., and Van Themaat, E. (2013). Mosaic genome structure of the barley powdery mildew pathogen and conservation of transcriptional programs in divergent hosts. Proc. Natl. Acad. Sci. U.S.A. 110, E2219–E2228. doi: 10.1073/pnas.1306807110
Jeon, J., Choi, J., Lee, G. W., Dean, R. A., and Lee, Y. H. (2013). Experimental evolution reveals genome-wide spectrum and dynamics of mutations in the rice blast fungus, Magnaporthe oryzae. PLoS ONE 8:e65416. doi: 10.1371/journal.pone.0065416
Joseph, B., Corwin, J. A., and Kliebenstein, D. J. (2015). Genetic variation in the nuclear and organellar genomes modulates stochastic variation in the metabolome, growth, and defense. PLoS Genet. 11:e1004779. doi: 10.1371/journal.pgen.1004779
Joseph, B., Corwin, J. A., Li, B., Atwell, S., and Kliebenstein, D. J. (2013a). Cytoplasmic genetic variation and extensive cytonuclear interactions influence natural variation in the metabolome. Elife 2:e00776. doi: 10.7554/elife.00776
Joseph, B., Corwin, J. A., Züest, T., Li, B., Iravani, M., Schaepman-Strub, G., et al. (2013b). Hierarchical nuclear and cytoplasmic genetic architectures for plant growth and defense within Arabidopsis. Plant Cell 25, 1929–1945. doi: 10.1105/tpc.113.112615
Kerssies, A., Boskervanzessen, A. I., Wagemakers, C. A. M., and Vankan, J. A. L. (1997). Variation in pathogenicity and DNA polymorphism among Botrytis cinerea isolates sampled inside and outside a glasshouse. Plant Dis. 81, 781–786. doi: 10.1094/PDIS.19188.8.131.521
Kliebenstein, D. J., Rowe, H. C., and Denby, K. J. (2005). Secondary metabolites influence Arabidopsis/Botrytis interactions: variation in host production and pathogen sensitivity. Plant J. 44, 25–36. doi: 10.1111/j.1365-313X.2005.02508.x
Kohn, L. M., Schaffer, M. R., Anderson, J. B., and Grunwald, N. J. (2008). Marker stability throughout 400 days of in vitro hyphal growth in the filamentous ascomycete, Sclerotinia sclerotiorum. Fungal Genet. Biol. 45, 613–617. doi: 10.1016/j.fgb.2007.09.013
Koonin, E. V., and Aravind, L. (2000). The NACHT family - a new group of predicted NTPases implicated in apoptosis and MHC transcription activation. Trends Biochem. Sci. 25, 223–224. doi: 10.1016/S0968-0004(00)01577-2
Kruglyak, S., Durrett, R. T., Schug, M. D., and Aquadro, C. F. (1998). Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations. Proc. Natl. Acad. Sci. U.S.A. 95, 10774–10778. doi: 10.1073/pnas.95.18.10774
Lamb, B. C., Mandaokar, S., Bahsoun, B., Grishkan, I., and Nevo, E. (2008). Differences in spontaneous mutation frequencies as a function of environmental stress in soil fungi at “Evolution Canyon,” Israel. Proc. Natl. Acad. Sci. U.S.A. 105, 5792–5796. doi: 10.1073/pnas.0801995105
Lévesque, C. A., Brouwer, H., Cano, L., Hamilton, J. P., Holt, C., Huitema, E., et al. (2010). Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire. Genome Biol. 11:R73. doi: 10.1186/gb-2010-11-7-r73. Available online at: http://www.genomebiology.com/2010/11/7/R73/about
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The Sequence alignment/map (SAM) format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352
Loftus, B. J., Fung, E., Roncaglia, P., Rowley, D., Amedeo, P., Bruno, D., et al. (2005). The genome of the basidiomycetous yeast and human pathogen Cryptococcus neoformans. Science 307, 1321–1324. doi: 10.1126/science.1103773
Long, Q., Rabanal, F. A., Meng, D., Huber, C. D., Farlow, A., Platzer, A., et al. (2013). Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden. Nat. Genet. 45, 884–890. doi: 10.1038/ng.2678
Martin, D. P., Lemey, P., Lott, M., Moulton, V., Posada, D., and Lefeuvre, P. (2010). RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics 26, 2462–2463. doi: 10.1093/bioinformatics/btq467
Martin, D. P., Posada, D., Crandall, K. A., and Williamson, C. (2005). A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res. Hum. Retroviruses 21, 98–102. doi: 10.1089/aid.2005.21.98
Nierman, W. C., Pain, A., Anderson, M. J., Wortman, J. R., Kim, H. S., Arroyo, J., et al. (2005). Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus. Nature 438, 1151–1156. doi: 10.1038/nature04332
Park, S. Y., Chi, M. H., Milgroom, M. G., Kim, H., Han, S. S., Kang, S., et al. (2010). Genetic stability of Magnaporthe oryzae during successive passages through rice plants and on artificial medium. Plant Pathol. J. 26, 313–320. doi: 10.5423/PPJ.2010.26.4.313
Ponce de León, I., Oliver, J. P., Castro, A., Gaggero, C., Bentancor, M., and Vidal, S. (2007). Erwinia carotovora elicitors and Botrytis cinerea activate defense responses in Physcomitrella patens. BMC Plant Biol. 7:52. doi: 10.1186/1471-2229-7-52
Ponce de León, I., Schmelz, E. A., Gaggero, C., Castro, A., Álvarez, A., and Montesano, M. (2012). Physcomitrella patens activates reinforcement of the cell wall, programmed cell death and accumulation of evolutionary conserved defence signals, such as salicylic acid and 12-oxo-phytodienoic acid, but not jasmonic acid, upon Botrytis cinerea infection. Mol. Plant Pathol. 13, 960–974. doi: 10.1111/j.1364-3703.2012.00806.x
Porcel, B. M., Denoeud, F., Opperdoes, F., Noel, B., Madoui, M. A., Hammarton, T. C., et al. (2014). The streamlined genome of Phytomonas spp. relative to human pathogenic kinetoplastids reveals a parasite tailored for plants. Plos Genet. 10:e1004007. doi: 10.1371/journal.pgen.1004007
Posada, D., and Crandall, K. A. (2001). Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc. Natl. Acad. Sci. U.S.A. 98, 13757–13762. doi: 10.1073/pnas.241370698
Rowe, H. C., and Kliebenstein, D. J. (2007). Elevated genetic variation within virulence-associated Botrytis cinerea polygalacturonase loci. Mol. Plant Microbe Interact. 20, 1126–1137. doi: 10.1094/MPMI-20-9-1126
Rowe, H. C., and Kliebenstein, D. J. (2008). Complex genetics control natural variation in Arabidopsis thaliana resistance to Botrytis cinerea. Genetics 180, 2237–2250. doi: 10.1534/genetics.108.091439
Rowe, H. C., Walley, J. W., Corwin, J., Chan, E. K., Dehesh, K., and Kliebenstein, D. J. (2010). Deficiencies in jasmonate-mediated plant defense reveal quantitative variation in Botrytis cinerea pathogenesis. PLoS Pathog. 6:e1000861. doi: 10.1371/journal.ppat.1000861
Schacherer, J., Shapiro, J. A., Ruderfer, D. M., and Kruglyak, L. (2009). Comprehensive polymorphism survey elucidates population structure of Saccharomyces cerevisiae. Nature 458, 342–345. doi: 10.1038/nature07670
Schumacher, J., Pradier, J. M., Simon, A., Traeger, S., Moraga, J., Collado, I. G., et al. (2012). Natural variation in the VELVET gene bcvel1 affects virulence and light-dependent differentiation in Botrytis cinerea. PLoS ONE 7:e47840. doi: 10.1371/journal.pone.0047840
Stefanato, F. L., Abou-Mansour, E., Buchala, A., Kretschmer, M., Mosbach, A., Hahn, M., et al. (2009). The ABC transporter BcatrB from Botrytis cinerea exports camalexin and is a virulence factor on Arabidopsis thaliana. Plant J. 58, 499–510. doi: 10.1111/j.1365-313X.2009.03794.x
Stergiopoulos, I., Zwiers, L. H., and De Waard, M. A. (2002). Secretion of natural and synthetic toxic compounds from filamentous fungi by membrane transporters of the ATP-binding cassette and major facilitator superfamily. European Journal of Plant Pathology 108, 719–734. doi: 10.1023/A:1020604716500
Swain, M. T., Tsai, I. J., Assefa, S. A., Newbold, C., Berriman, M., and Otto, T. D. (2012). A post-assembly genome-improvement toolkit (PAGIT) to obtain annotated genomes from contigs. Nat. Protoc. 7, 1260–1284. doi: 10.1038/nprot.2012.068
Tan, Z., Hu, Z., Yang, Z., Yu, B., and Xu, C. (2012). Framework for dissection of complex cytonuclear epistasis by a two-dimensional genome scan. Chin. Sci. Bull. 57, 2675–2680. doi: 10.1007/s11434-012-5116-0
Tang, Z., Yang, Z., Hu, Z., Zhang, D., Lu, X., Deng, D., et al. (2013). Cytonuclear epistatic quantitative trait locus mapping for plant height and ear height in maize. Mol. Breed. 31, 1–14. doi: 10.1007/s11032-012-9762-3
Torriani, S. F., Goodwin, S. B., Kema, G. H., Pangilinan, J. L., and McDonald, B. A. (2008). Intraspecific comparison and annotation of two complete mitochondrial genome sequences from the plant pathogenic fungus Mycosphaerella graminicola. Fungal Genet. Biol. 45, 628–637. doi: 10.1016/j.fgb.2007.12.005
Wicker, T., Oberhaensli, S., Parlange, F., Buchmann, J. P., Shatalina, M., Roffler, S., et al. (2013). The wheat powdery mildew genome shows the unique evolution of an obligate biotroph. Nat. Genet. 45, 1092–1096. doi: 10.1038/ng.2704
Wilkening, S., Tekkedil, M. M., Lin, G., Fritsch, E. S., Wei, W., Gagneur, J., et al. (2013). Genotyping 1000 yeast strains by next-generation sequencing. BMC Genomics 14:90. doi: 10.1186/1471-2164-14-90
Yourman, L. F., Jeffers, S. N., and Dean, R. A. (2000). Genetic analysis of isolates of Botrytis cinerea sensitive and resistant to benzimidazole and dicarboximide fungicides. Phytopathology 90, 851–859. doi: 10.1094/PHYTO.2000.90.8.851
Yourman, L. F., Jeffers, S. N., and Dean, R. A. (2001). Phenotype instability in Botrytis cinerea in the absence of benzimidazole and dicarboximide fungicides. Phytopathology 91, 307–315. doi: 10.1094/PHYTO.2001.91.3.307
Zhan, J., Pettway, R. E., and McDonald, B. A. (2003). The global genetic structure of the wheat pathogen Mycosphaerella graminicola is characterized by high nuclear diversity, low mitochondrial diversity, regular recombination, and gene flow. Fungal Genet. Biol. 38, 286–297. doi: 10.1016/S1087-1845(02)00538-8
Zhao, J., Gladieux, P., Hutchison, E., Bueche, J., Hall, C., Perraudeau, F., et al. (2015). Identification of allorecognition loci in Neurospora crassa by genomics and evolutionary approaches. Mol. Biol. Evol. 2015:msv125. doi: 10.1093/molbev/msv125
Keywords: Botrytis cinerea, genome resequencing, standing genetic variation, diversity, vegetative incompatibility loci
Citation: Atwell S, Corwin JA, Soltis NE, Subedy A, Denby KJ and Kliebenstein DJ (2015) Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity. Front. Microbiol. 6:996. doi: 10.3389/fmicb.2015.00996
Received: 20 May 2015; Accepted: 07 September 2015;
Published: 24 September 2015.
Edited by:Michael H. Perlin, University of Louisville, USA
Reviewed by:Ernesto P. Benito, University of Salamanca, Spain
Dominik Begerow, Ruhr-Universität Bochum, Germany
Copyright © 2015 Atwell, Corwin, Soltis, Subedy, Denby and Kliebenstein. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel J. Kliebenstein, Department of Plant Sciences, University of California, Davis, One Shields Ave, Davis, CA 95616, USA, email@example.com