Crossing Plasmopara viticola strains in controlled conditions to uncover the genomic bases of downy mildew resistance breakdown in grapevine

The biotrophic oomycete Plasmopara viticola is the causal agent of downy mildew, one of the major grapevine diseases. We report here a successful strategy to cross compatible strains of this pathogen and obtain a large viable progeny. We used this method to study the offspring between two P. viticola strains able to overcome two major grapevine resistance factors (Rpv3, Rpv10). Thanks to the genomic resources now available, we will genotype this progeny to build an unprecedented linkage map and uncover the genomic bases of grapevine resistance breakdown displayed by P. viticola virulent strains.

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

QTL mapping in Fusarium graminearum identified an allele of FgVe1 involved in reduced aggressiveness

Fusarium graminearum is one of the most frequent causal agents of the Fusarium Head Blight, a cereal disease spread throughout the world, reducing grain production and quality. F. graminearum isolates are genetically and phenotypically highly diverse. Notably, remarkable variations of aggressiveness between isolates have been observed, which could reflect an adaptive potential of this pathogen. In this study, we aimed to characterize the genetic basis of aggressiveness variation observed in an F1 population (n = 94), for which genome sequences of both parental strains are available. Aggressiveness was assessed by a panel of in planta and in vitro proxies during two phenotyping trials including, among others, disease severity and mycotoxin accumulation in wheat spike. One major and single QTL was mapped for all the traits measured, on chromosome I, that explained up to 90% of the variance for disease severity. The confidence interval at the QTL spanned 1.2 Mb and contained 428 genes on the reference genome. Of these, four candidates were selected based on the postulate that a non-synonymous mutation affecting protein function may be responsible for phenotypic differences. Finally, a new mutation was identified and functionally validated in the gene FgVe1, coding for a velvet protein known to be involved in pathogenicity and secondary metabolism production in several fungi.

Evolution of Fusarium tricinctum and Fusarium avenaceum mitochondrial genomes is driven by mobility of introns and of a new type of palindromic microsatellite repeats

Background

Increased contamination of European and Asian wheat and barley crops with “emerging” mycotoxins such as enniatins or beauvericin, produced by Fusarium avenaceum and Fusarium tricinctum, suggest that these phylogenetically close species could be involved in future food-safety crises.

Results

The mitochondrial genomes of F. tricinctum strain INRA104 and F. avenaceum strain FaLH27 have been annotated. A comparative analysis was carried out then extended to a set of 25 wild strains. Results show that they constitute two distinct species, easily distinguished by their mitochondrial sequences. The mitochondrial genetic variability is mainly located within the intergenic regions. Marks of variations show they have evolved (i) by Single Nucleotide Polymorphisms (SNPs), (ii) by length variations mediated by insertion/deletion sequences (Indels), and (iii) by length mutations generated by DNA sliding events occurring in mononucleotide (A)_n or (T)_n microsatellite type sequences arranged in a peculiar palindromic organization. The optionality of these palindromes between both species argues for their mobility. The presence of Indels and SNPs in palindrome neighbouring regions suggests their involvement in these observed variations. Moreover, the intraspecific and interspecific variations in the presence/absence of group I introns suggest a high mobility, resulting from several events of gain and loss during short evolution periods. Phylogenetic analyses of intron orthologous sequences suggest that most introns could have originated from lateral transfers from phylogenetically close or distant species belonging to various Ascomycota genera and even to the Basidiomycota fungal division.

Conclusions

Mitochondrial genome evolution between F. tricinctum and F. avenaceum is mostly driven by two types of mobile genetic elements, implicated in genome polymorphism. The first one is represented by group I introns. Indeed, both genomes harbour optional (inter- or intra-specifically) group I introns, all carrying putatively functional hegs, arguing for a high mobility of these introns during short evolution periods. The gain events were shown to involve, for most of them, lateral transfers between phylogenetically distant species. This study has also revealed a new type of mobile genetic element constituted by a palindromic arrangement of (A) n and (T) n microsatellite sequences whose presence was related to occurrence of SNPs and Indels in the neighbouring regions.

Landscape of genomic diversity and host adaptation in Fusarium graminearum

BACKGROUND

Fusarium graminearum is one of the main causal agents of the Fusarium Head Blight, a worldwide disease affecting cereal cultures, whose presence can lead to contaminated grains with chemically stable and harmful mycotoxins. Resistant cultivars and fungicides are frequently used to control this pathogen, and several observations suggest an adaptation of F. graminearum that raises concerns regarding the future of current plant disease management strategies. To understand the genetic basis as well as the extent of its adaptive potential, we investigated the landscape of genomic diversity among six French isolates of F. graminearum, at single-nucleotide resolution using whole-genome re-sequencing.

RESULTS

A total of 242,756 high-confidence genetic variants were detected when compared to the reference genome, among which 96% are single nucleotides polymorphisms. One third of these variants were observed in all isolates. Seventy-seven percent of the total polymorphism is located in 32% of the total length of the genome, comprising telomeric/subtelomeric regions as well as discrete interstitial sections, delineating clear variant enriched genomic regions- 7.5 times in average. About 80% of all the F. graminearum protein-coding genes were found polymorphic. Biological functions are not equally affected: genes potentially involved in host adaptation are preferentially located within polymorphic islands and show greater diversification rate than genes fulfilling basal functions. We further identified 29 putative effector genes enriched with non-synonymous effect mutation.

CONCLUSIONS

Our results highlight a remarkable level of polymorphism in the genome of F. graminearum distributed in a specific pattern. Indeed, the landscape of genomic diversity follows a bi-partite organization of the genome according to polymorphism and biological functions. We measured, for the first time, the level of sequence diversity for the entire gene repertoire of F. graminearum and revealed that the majority are polymorphic. Those assumed to play a role in host-pathogen interaction are discussed, in the light of the subsequent consequences for host adaptation. The annotated genetic variants discovered for this major pathogen are valuable resources for further genetic and genomic studies.

A detailed analysis of the recombination landscape of the button mushroom Agaricus bisporus var. bisporus

The button mushroom (Agaricus bisporus) is one of the world's most cultivated mushroom species, but in spite of its economic importance generation of new cultivars by outbreeding is exceptional. Previous genetic analyses of the white bisporus variety, including all cultivars and most wild isolates revealed that crossing over frequencies are low, which might explain the lack of introducing novel traits into existing cultivars. By generating two high quality whole genome sequence assemblies (one de novo and the other by improving the existing reference genome) of the first commercial white hybrid Horst U1, a detailed study of the crossover (CO) landscape was initiated. Using a set of 626 SNPs in a haploid offspring of 139 single spore isolates and whole genome sequencing on a limited number of homo- and heterokaryotic single spore isolates, we precisely mapped all COs showing that they are almost exclusively restricted to regions of about 100kb at the chromosome ends. Most basidia of A. bisporus var. bisporus produce two spores and pair preferentially via non-sister nuclei. Combined with the COs restricted to the chromosome ends, these spores retain most of the heterozygosity of the parent thus explaining how present-day white cultivars are genetically so close to the first hybrid marketed in 1980. To our knowledge this is the first example of an organism which displays such specific CO landscape.

Genetic Analyses of the Internal Transcribed Spacer Sequences Suggest Introgression and Duplication in the Medicinal Mushroom Agaricus subrufescens

The internal transcribed spacer (ITS) region of the nuclear ribosomal RNA gene cluster is widely used in fungal taxonomy and phylogeographic studies. The medicinal and edible mushroom Agaricus subrufescens has a worldwide distribution with a high level of polymorphism in the ITS region. A previous analysis suggested notable ITS sequence heterogeneity within the wild French isolate CA487. The objective of this study was to investigate the pattern and potential mechanism of ITS sequence heterogeneity within this strain. Using PCR, cloning, and sequencing, we identified three types of ITS sequences, A, B, and C with a balanced distribution, which differed from each other at 13 polymorphic positions. The phylogenetic comparisons with samples from different continents revealed that the type C sequence was similar to those found in Oceanian and Asian specimens of A. subrufescens while types A and B sequences were close to those found in the Americas or in Europe. We further investigated the inheritance of these three ITS sequence types by analyzing their distribution among single-spore isolates from CA487. In this analysis, three co-dominant markers were used firstly to distinguish the homokaryotic offspring from the heterokaryotic offspring. The homokaryotic offspring were then analyzed for their ITS types. Our genetic analyses revealed that types A and B were two alleles segregating at one locus ITSI, while type C was not allelic with types A and B but was located at another unlinked locus ITSII. Furthermore, type C was present in only one of the two constitutive haploid nuclei (n) of the heterokaryotic (n+n) parent CA487. These data suggest that there was a relatively recent introduction of the type C sequence and a duplication of the ITS locus in this strain. Whether other genes were also transferred and duplicated and their impacts on genome structure and stability remain to be investigated.

The Genetic Linkage Map of the Medicinal Mushroom Agaricus subrufescens Reveals Highly Conserved Macrosynteny with the Congeneric Species Agaricus bisporus

Comparative linkage mapping can rapidly facilitate the transfer of genetic information from model species to orphan species. This macrosynteny analysis approach has been extensively used in plant species, but few example are available in fungi, and even fewer in mushroom crop species. Among the latter, the Agaricus genus comprises the most cultivable or potentially cultivable species. Agaricus bisporus, the button mushroom, is the model for edible and cultivable mushrooms. We have developed the first genetic linkage map for the basidiomycete A. subrufescens, an emerging mushroom crop known for its therapeutic properties and potential medicinal applications. The map includes 202 markers distributed over 16 linkage groups (LG), and covers a total length of 1701 cM, with an average marker spacing of 8.2 cM. Using 96 homologous loci, we also demonstrated the high level of macrosynteny with the genome of A. bisporus The 13 main LG of A. subrufescens were syntenic to the 13 A. bisporus chromosomes. A disrupted synteny was observed for the three remaining A. subrufescens LG. Electronic mapping of a collection of A. subrufescens expressed sequence tags on A. bisporus genome showed that the homologous loci were evenly spread, with the exception of a few local hot or cold spots of homology. Our results were discussed in the light of Agaricus species evolution process. The map provides a framework for future genetic or genomic studies of the medicinal mushroom A. subrufescens.

Evidence for amphithallism and broad geographical hybridization potential among Agaricus subrufescens isolates from Brazil, France, and Thailand

Agaricus subrufescens is a cultivated edible and medicinal mushroom. Its known geographical distribution encompasses the Americas, Europe, Oceania, and Asia. The objective of this study was to assess mating compatibility and interfertility of strains originating from Brazil, France, and Thailand. Progeny of each strain were analyed with codominant molecular markers. Multilocus genotype tests revealed that the three strains were amphithallic with percentages of heterokaryotic single spore progenies of 75% for the Thai strain and around 40% for the Brazilian and French strains. In mating tests A. subrufescens had a multiallelic unifactorial system of sexual incompatibility. The three parent strains were interfertile based on experimental pairings of single-spore isolates, the recovery of hybrid heterokaryons from compatible matings, and the ability of hybrids to produce mushrooms and fertile spores. This biological approach supports the inclusion of the European strains within the species and the extension of the geographical distribution range to Asia. Our data should help to develop breeding strategies and to better manage and exploit the diversity existing in A. subrufescens.

Development of polymorphic microsatellite markers issued from pyrosequencing technology for the medicinal mushroom Agaricus subrufescens

The recently described procedure of microsatellite-enriched library pyrosequencing was used to isolate microsatellite loci in the gourmet and medicinal mushroom Agaricus subrufescens. Three hundred and five candidate loci containing at least one simple sequence repeats (SSR) locus and for which primers design was successful, were obtained. From a subset of 95 loci, 35 operational and polymorphic SSR markers were developed and characterized on a sample of 14 A. subrufescens genotypes from diverse origins. These SubSSR markers each displayed from two to 10 alleles with an average of 4.66 alleles per locus. The observed heterozygosity ranged from 0 to 0.71. Several multiplex combinations can be set up, making it possible to genotype up to six markers easily and simultaneously. Cross-amplification in some closely congeneric species was successful for a subset of loci. The 35 microsatellite markers developed here provide a highly valuable molecular tool to study genetic diversity and reproductive biology of A. subrufescens.

Genetic linkage mapping in fungi: current state, applications, and future trends

Genetic mapping is a basic tool for eukaryotic genomic research. Linkage maps provide insights into genome organization and can be used for genetic studies of traits of interest. A genetic linkage map is a suitable support for the anchoring of whole genome sequences. It allows the localization of genes of interest or quantitative trait loci (QTL) and map-based cloning. While genetic mapping has been extensively used in plant or animal models, this discipline is more recent in fungi. The present article reviews the current status of genetic linkage map research in fungal species. The process of linkage mapping is detailed, from the development of mapping populations to the construction of the final linkage map, and illustrated based on practical examples. The range of specific applications in fungi is browsed, such as the mapping of virulence genes in pathogenic species or the mapping of agronomically relevant QTL in cultivated edible mushrooms. Future prospects are finally discussed in the context of the most recent advances in molecular techniques and the release of numerous fungal genome sequences.

Relationship between yield components and partial resistance to Lecanicillium fungicola in the button mushroom, Agaricus bisporus, assessed by quantitative trait locus mapping

Dry bubble, caused by Lecanicillium fungicola, is one of the most detrimental diseases affecting button mushroom cultivation. In a previous study, we demonstrated that breeding for resistance to this pathogen is quite challenging due to its quantitative inheritance. A second-generation hybrid progeny derived from an intervarietal cross between a wild strain and a commercial cultivar was characterized for L. fungicola resistance under artificial inoculation in three independent experiments. Analysis of quantitative trait loci (QTL) was used to determine the locations, numbers, and effects of genomic regions associated with dry-bubble resistance. Four traits related to resistance were analyzed. Two to four QTL were detected per trait, depending on the experiment. Two genomic regions, on linkage group X (LGX) and LGVIII, were consistently detected in the three experiments. The genomic region on LGX was detected for three of the four variables studied. The total phenotypic variance accounted for by all QTL ranged from 19.3% to 42.1% over all traits in all experiments. For most of the QTL, the favorable allele for resistance came from the wild parent, but for some QTL, the allele that contributed to a higher level of resistance was carried by the cultivar. Comparative mapping with QTL for yield-related traits revealed five colocations between resistance and yield component loci, suggesting that the resistance results from both genetic factors and fitness expression. The consequences for mushroom breeding programs are discussed.

Quantitative genetics to dissect the fungal-fungal interaction between Lecanicillium verticillium and the white button mushroom Agaricus bisporus

Lecanicillium fungicola (formerly Verticillium fungicola) is responsible for dry bubble disease in the white button mushroom Agaricus bisporus. Selection for resistance to this pathogen raises an important challenge for mushroom breeders. We have investigated the inheritance of resistance to dry bubble under artificial inoculation in three independent experiments, using a progeny of 89 hybrids derived from an intervarietal A. bisporus var. bisporus×A. bisporus var. burnettii cross. Overall, phenotypic correlations were highly significant between the different experiments. Principal component analysis, together with analysis of variance results stated that the disease reactions were accurately assessed using the percentage of bubbles (PB) and the percentage of spotty cap mushrooms (PS) separately rather than with the combination of both. An original contribution of this study lies in the effective use of area under the disease-progress curve (AUDPC) to describe the dry bubble resistance. The continuous phenotypic distribution observed for the resistance traits suggested that tolerance to dry bubble was under polygenic control. Heritability estimates for either PB or AUDPC were high (0.67-0.86) while it was inconsistent for PS (0.33-0.68) suggesting a strong impact of the environment on this latter trait. Earliness and latent period were found highly correlated with disease incidence. The earliest strains appeared to be the most resistant ones. These results contribute to disentangle the complex fungal-fungal A. bisporus / L. fungicola interaction and to provide genetic basis as a prerequisite for mushroom breeding program.

An expanded genetic linkage map of an intervarietal Agaricus bisporus var. bisporusxA. bisporus var. burnettii hybrid based on AFLP, SSR and CAPS markers sheds light on the recombination behaviour of the species

A genetic linkage map for the edible basidiomycete Agaricus bisporus was constructed from 118 haploid homokaryons derived from an intervarietal A. bisporus var. bisporus x A. bisporus var. burnettii hybrid. Two hundred and thirty-one AFLP, 21 SSR, 68 CAPS markers together with the MAT, BSN, PPC1 loci and one allozyme locus (ADH) were evenly spread over 13 linkage groups corresponding to the chromosomes of A. bisporus. The map covers 1156cM, with an average marker spacing of 3.9cM and encompasses nearly the whole genome. The average number of crossovers per chromosome per individual is 0.86. Normal recombination over the entire genome occurs in the heterothallic variety, burnettii, contrary to the homothallic variety, bisporus, which showed adaptive genome-wide suppressed recombination. This first comprehensive genetic linkage map for A. bisporus provides foundations for quantitative trait analyses and breeding programme monitoring, as well as genome organisation studies.