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

A SLAF-based high-density genetic map construction and genetic architecture of thermotolerant traits in maize (Zea mays L.)

The leaf scorching trait at flowering is a crucial thermosensitive phenotype in maize under high temperature stress (HS), yet the genetic basis of this trait remains poorly understood. In this study, we genotyped a 254 RIL-F2:8 population, derived from the leaf scorch-free parental inbred line Abe2 and the leaf scorching maternal inbred line B73, using the specific-locus amplified fragment sequencing (SLAF-seq) method. A total of 10,112 polymorphic SLAF markers were developed, and a high-density genetic map with a total length of 1,475.88 cM was constructed. The average sequencing depth of the parents was 55.23X, and that of the progeny was 12.53X. Then, we identified a total of 16 QTLs associated with thermotolerant traits at flowering, of which four QTLs of leaf scorching damage (LS) were distributed on chromosomes 1 (qLS1), 2 (qLS2.1, qLS2.2) and 3 (qLS3), which could explain 19.73% of phenotypic variation. Combining one qLS1 locus with QTL-seq results led to the identification of 6 candidate genes. Expression experiments and sequence variation indicated that Zm00001d033328, encoding N-acetyl-gamma-glutamyl-phosphate reductase, was the most likely candidate gene controlling thermotolerant traits at flowering. In summary, the high-density genetic map and genetic basis of thermotolerant traits lay a critical foundation for mapping other complex traits and identifying the genes associated with thermotolerant traits in maize.

Genetic Linkage and Physical Mapping for an Oyster Mushroom (Pleurotus cornucopiae) and Quantitative Trait Locus Analysis for Cap Color

Oyster mushrooms are grown commercially worldwide, especially in many developing countries, for their easy cultivation and high biological efficiency. Pleurotus cornucopiae is one of the main oyster mushroom species because of its gastronomic value and nutraceutical properties. Cap color is an important trait, since consumers prefer dark mushrooms, which are now represented by only a small portion of the commercial varieties. Breeding efforts are required to improve quality-related traits to satisfy various demands of consumers. Here, we present a saturated genetic linkage map of P. cornucopiae constructed by using a segregating population of 122 monokaryons and 3,449 single nucleotide polymorphism (SNP) markers generated by the 2b-RAD approach. The map contains 11 linkage groups covering 961.6 centimorgans (cM), with an average marker spacing of 0.27 cM. The genome of P. cornucopiae was de novo sequenced, resulting in 425 scaffolds (>1,000 bp) with a total genome size of 35.1 Mb. The scaffolds were assembled to the pseudochromosome level with the assistance of the genetic linkage map. A total of 97% SNP markers (3,357) were physically localized on 140 scaffolds that were assigned to 11 pseudochromosomes, with a total of 32.5 Mb, representing 92.5% of the whole genome. Six quantitative trait loci (QTL) controlling cap color of P. cornucopiae were detected, accounting for a total phenotypic variation of 65.6%, with the highest value for the QTL on pseudochromosome 5 (18%). The results of our study provide a solid base for marker-assisted breeding for agronomic traits and especially for studies on biological mechanisms controlling cap color in oyster mushrooms. IMPORTANCE Oyster mushrooms are produced and consumed all over the world. Pleurotus cornucopiae is one of the main oyster mushroom species. Dark-cap oyster mushrooms are becoming more and more popular with consumers, but dark varieties are rare on the market. Prerequisites for efficient breeding programs are the availability of high-quality whole genomes and genetic linkage maps. Genetic studies to fulfill some of these prerequisites have hardly been done for P. cornucopiae. In this study, we de novo sequenced the genome and constructed a saturated genetic linkage map for P. cornucopiae. The genetic linkage map was effectively used to assist the genome assembly and identify QTL that genetically control the trait cap color. As well, the genome characteristics of P. cornucopiae were compared to the closely related species Pleurotus ostreatus. The results provided a basis for understanding the genetic background and marker-assisted breeding of this economically important mushroom species.

Tolerance to oxidative stress is associated with both oxidative stress response and inherent growth in a fungal wheat pathogen

Reactive oxygen species are toxic byproducts of aerobic respiration that are also important in mediating a diversity of cellular functions. Reactive oxygen species form an important component of plant defenses to inhibit microbial pathogens during pathogen-plant interactions. Tolerance to oxidative stress is likely to make a significant contribution to the viability and pathogenicity of plant pathogens, but the complex network of oxidative stress responses hinders identification of the genes contributing to this trait. Here, we employed a forward genetic approach to investigate the genetic architecture of oxidative stress tolerance in the fungal wheat pathogen Zymoseptoria tritici. We used quantitative trait locus (QTL) mapping of growth and melanization under axenic conditions in two cross-populations to identify genomic regions associated with tolerance to oxidative stress. We found that QTLs associated with growth under oxidative stress as well as inherent growth can affect oxidative stress tolerance, and we identified two uncharacterized genes in a major QTL associated with this trait. Our data suggest that melanization does not affect tolerance to oxidative stress, which differs from what was found for animal pathogens. This study provides a whole-genome perspective on the genetic basis of oxidative stress tolerance in a plant pathogen.

The Assembled and Annotated Genome of the Fairy-Ring Fungus Marasmius oreades

Marasmius oreades is a basidiomycete fungus that grows in so called "fairy rings," which are circular, underground mycelia common in lawns across temperate areas of the world. Fairy rings can be thought of as natural, long-term evolutionary experiments. As each ring has a common origin and expands radially outwards over many years, different sectors will independently accumulate mutations during growth. The genotype can be followed to the next generation, as mushrooms producing the sexual spores are formed seasonally at the edge of the ring. Here, we present new genomic data from 95 single-spore isolates of the species, which we used to construct a genetic linkage map and an updated version of the genome assembly. The 44-Mb assembly was anchored to 11 linkage groups, producing chromosome-length scaffolds. Gene annotation revealed 13,891 genes, 55% of which contained a pfam domain. The repetitive fraction of the genome was 22%, and dominated by retrotransposons and DNA elements of the KDZ and Plavaka groups. The level of assembly contiguity we present is so far rare in mushroom-forming fungi, and we expect studies of genomics, transposons, phylogenetics, and evolution to be facilitated by the data we present here of the iconic fairy-ring mushroom.

Genetic linkage map construction and quantitative trait loci mapping of agronomic traits in Gloeostereum incarnatum

Gloeostereum incarnatum is an edible medicinal mushroom widely grown in China. Using the whole genome of G. incarnatum, simple sequence repeat (SSR) markers were developed and synthetic primers were designed to construct its first genetic linkage map. The 1,048.6 cm map is composed of 10 linkage groups and contains 183 SSR markers. In total, 112 genome assembly sequences were anchored, representing 16.43 Mb and covering 46.41% of the genome. Selfing populations were used for quantitative trait loci (QTL) targeting, and the composite interval mapping method was used to co-localize the mycelium growth rate (potato dextrose agar and sawdust), growth period, yield and fruiting body length, and width and thickness. The 14 QTLs of agronomic traits had LOD values of 3.20-6.51 and contribution rates of 2.22-13.18%. No linkage relationship was found between the mycelium growth rate and the growth period, but a linkage relationship was observed among the length, width and thickness of the fruiting bodies. Using NCBI's BLAST alignment, the genomic sequences corresponding to the QTL regions were compared, and a TPR-like protein candidate gene was selected. Using whole-genome data, 138 candidate genes were found in four sequence fragments of two SSR markers located in the same scaffold. The genetic map and QTLs established in this study will aid in developing selective markers for agronomic traits and identifying corresponding genes, thereby providing a scientific basis for the further gene mapping of quantitative traits and the marker-assisted selection of functional genes in G. incarnatum breeding programs.

Chromosomal assembly and analyses of genome-wide recombination rates in the forest pathogenic fungus Armillaria ostoyae

Recombination shapes the evolutionary trajectory of populations and plays an important role in the faithful transmission of chromosomes during meiosis. Levels of sexual reproduction and recombination are important properties of host-pathogen interactions because the speed of antagonistic co-evolution depends on the ability of hosts and pathogens to generate genetic variation. However, our understanding of the importance of recombination is limited because large taxonomic groups remain poorly investigated. Here, we analyze recombination rate variation in the basidiomycete fungus Armillaria ostoyae, which is an aggressive pathogen on a broad range of conifers and other trees. We analyzed a previously constructed, dense genetic map based on 198 single basidiospore progeny from a cross. Progeny were genotyped at a genome-wide set of single-nucleotide polymorphism (SNP) markers using double digest restriction site associated DNA sequencing. Based on a linkage map of on 11,700 SNPs spanning 1007.5 cM, we assembled genomic scaffolds into 11 putative chromosomes of a total genome size of 56.6 Mb. We identified 1984 crossover events among all progeny and found that recombination rates were highly variable along chromosomes. Recombination hotspots tended to be in regions close to the telomeres and were more gene-poor than the genomic background. Genes in proximity to recombination hotspots were encoding on average shorter proteins and were enriched for pectin degrading enzymes. Our analyses enable more powerful population and genome-scale studies of a major tree pathogen.

A Resequencing-Based Ultradense Genetic Map of Hericium erinaceus for Anchoring Genome Sequences and Identifying Genetic Loci Associated With Monokaryon Growth

Hericium erinaceus has attracted tremendous interest owing to its compelling health-promoting properties. However, breeding of elite cultivars of H. erinaceus is hindered by the lack of a genetic and molecular toolbox. Here, we performed resequencing analysis of 127 F1 single-spore isolates and constructed the first high-resolution genetic map of H. erinaceus. With the use of recombination bins as markers, an ultradense genetic map consisting of 1,174 bins (including 37,082 single-nucleotide polymorphisms) was generated. This newly developed genetic map covered 1,096.5 cM, with an average bin spacing of 0.95 cM. High collinearity between genetic map and H. erinaceus genome assembly was revealed by aligning scaffolds to this genetic map using bin markers as anchors. The application of this newly developed genetic map in quantitative trait locus (QTL) mapping was also elucidated, and four QTLs for monokaryon growth were recovered. One QTL, mgr1, which contributes 12.1% of growth variations, was located near the mating type A (MAT-A) loci. Overall, this newly constructed high-resolution genetic map (or bin map) could be used as reference in future genetic, genomic, and breeding studies on H. erinaceus.

Construction of a genetic linkage map of Lentinula edodes based on SSR, SRAP and TRAP markers

Genetic mapping is a basic tool for eukaryotic genomic research. It allows the localization of genes or quantitative trait loci (QTLs) and map-based cloning. In this study, we constructed a linkage map based on DNA samples from a commercial strain L808, including two parental monokaryons and 93 single spore isolates considered with segregating to 1:1:1:1 at four mating types (A₁B₁, A₁B₂, A₂B₁ and A₂B₂). Using Simple Sequence Repeats (SSR), Sequence Related Amplified Polymorphism (SRAP), Target Region Amplified Polymorphism (TRAP) molecular markers, 182 molecular markers and two mating factors were located on 11 linkage groups (LGs). The total length of the map was 948.083 centimorgan (cM), with an average marker interval distance of 4.817 cM. Only two gaps spanning more than 20 cM was observed. The probability of 20 cM, 10 cM, 5 cM genetic distance cover one marker was 99.68%, 94.36%, 76.43% in our genetic linkage map, respectively. This is the first linkage map of Lentinula edodes using SSR markers, which provides essential information for quantitative trait analyses and improvement of genome assembly.

Comparative Methods for Molecular Determination of Host-Specificity Factors in Plant-Pathogenic Fungi

Many plant-pathogenic fungi are highly host-specific. In most cases, host-specific interactions evolved at the time of speciation of the respective host plants. However, host jumps have occurred quite frequently, and still today the greatest threat for the emergence of new fungal diseases is the acquisition of infection capability of a new host by an existing plant pathogen. Understanding the mechanisms underlying host-switching events requires knowledge of the factors determining host-specificity. In this review, we highlight molecular methods that use a comparative approach for the identification of host-specificity factors. These cover a wide range of experimental set-ups, such as characterization of the pathosystem, genotyping of host-specific strains, comparative genomics, transcriptomics and proteomics, as well as gene prediction and functional gene validation. The methods are described and evaluated in view of their success in the identification of host-specificity factors and the understanding of their functional mechanisms. In addition, potential methods for the future identification of host-specificity factors are discussed.

A high-resolution genetic map of the cereal crown rot pathogen Fusarium pseudograminearum provides a near-complete genome assembly

Fusarium pseudograminearum is an important pathogen of wheat and barley, particularly in semi-arid environments. Previous genome assemblies for this organism were based entirely on short read data and are highly fragmented. In this work, a genetic map of F. pseudograminearum has been constructed for the first time based on a mapping population of 178 individuals. The genetic map, together with long read scaffolding of a short read-based genome assembly, was used to give a near-complete assembly of the four F. pseudograminearum chromosomes. Large regions of synteny between F. pseudograminearum and F. graminearum, the related pathogen that is the primary causal agent of cereal head blight disease, were previously proposed in the core conserved genome, but the construction of a genetic map to order and orient contigs is critical to the validation of synteny and the placing of species-specific regions. Indeed, our comparative analyses of the genomes of these two related pathogens suggest that rearrangements in the F. pseudograminearum genome have occurred in the chromosome ends. One of these rearrangements includes the transposition of an entire gene cluster involved in the detoxification of the benzoxazolinone (BOA) class of plant phytoalexins. This work provides an important genomic and genetic resource for F. pseudograminearum, which is less well characterized than F. graminearum. In addition, this study provides new insights into a better understanding of the sexual reproduction process in F. pseudograminearum, which informs us of the potential of this pathogen to evolve.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

Construction of a genetic linkage map and QTL mapping of agronomic traits in Auricularia auricula-judae

Auricularia auricula-judae is a traditional edible fungus that is cultivated widely in China. In this study, a genetic linkage map for A. auricula-judae was constructed using a mapping population consisting of 138 monokaryons derived from a hybrid strain (A119-5). The monokaryotic parent strains A14-5 and A18-119 were derived from two cultivated varieties, A14 (Qihei No. 1) and A18 (Qihei No. 2), respectively. In total, 130 simple sequence repeat markers were mapped. These markers were developed using the whole genome sequence of A. auricula-judae and amplified in A14-5, A18- 119, and the mapping population. The map consisted of 11 linkage groups (LGs) spanning 854 cM, with an average interval length of 6.57 cM. A testcross population was derived from crossing between the monokaryon A184-57 (from the wild strain A184 as a tester strain) and the mapping population. Important agronomic trait-related QTLs, including mycelium growth rate on potato dextrose agar for the mapping population, mycelium growth rate on potato dextrose agar and sawdust for the testcross population, growth period (days from inoculation to fruiting body harvesting), and yield for the testcross population, were identified using the composite interval mapping method. Six mycelium growth raterelated QTLs were identified on LG1 and LG4, two growth period-related QTLs were identified on LG2, and three yieldrelated QTLs were identified on LG2 and LG6. The results showed no linkage relationship between mycelium growth rate and growth period. The present study provides a foundation for locating genes for important agronomic characteristics in A. auricula-judae in the future.

High-density genetic mapping identifies the genetic basis of a natural colony morphology mutant in the root rot pathogen Armillaria ostoyae

Filamentous fungi exhibit a broad spectrum of heritable growth patterns and morphological variations reflecting the adaptation of the different species to distinct ecological niches. But also within species, isolates show considerable variation in growth rates and other morphological characteristics. The genetic basis of this intraspecific variation in mycelial growth and morphology is currently poorly understood. By chance, a growth mutant in the root rot pathogen Armillaria ostoyae was discovered. The mutant phenotype was characterized by extremely compact and slow growth, as well as shorter aerial hyphae and hyphal compartments in comparison to the wildtype phenotype. Genetic analysis revealed that the abnormal phenotype is caused by a recessive mutation, which segregates asa single locus in sexual crosses. In order to identify the genetic basis of the mutant phenotype, we performed a quantitative trait locus (QTL) analysis. A mapping population of 198 haploid progeny was genotyped at 11,700 genome-wide single nucleotide polymorphisms (SNPs) making use of double digest restriction site associated DNA sequencing (ddRADseq). In accordance with the genetic analysis, a single significant QTL was identified for the abnormal growth phenotype. The QTL confidence interval spans a narrow, gene dense region of 87kb in the A. ostoyae genome which contains 37 genes. Overall, our study reports the first high-density genetic map for an Armillaria species and shows its successful application in forward genetics by resolving the genetic basis of a mutant phenotype with a severe defect in hyphal growth.

Developments in breeding of Agaricus bisporus var. bisporus: progress made and technical and legal hurdles to take

True breeding of button mushrooms has hardly been done in the last decades, despite this species being one of the most cultivated mushrooms worldwide. Research done in the last 20 years has identified and characterised new germplasm and improved our understanding of the genetic base for some traits. A substantial collection of wild-collected strains is now available and partly characterised for a number of important traits such as disease resistance and yield. Most of the variations found in a number of important agronomic traits have a considerable heritability and are thus useful for breeding. Genetic marker technology has also developed considerably for this mushrooms in the last decade and used to identify quantitative trait loci (QTL) for important agronomic traits. This progress has, except for one example, not resulted so far into new commercially varieties. One of the reasons lies in the typical life cycle of the button mushroom Agaricus bisporus var. bisporus which hampers breeding. Joint investment is needed to solve technical problems in breeding. Special attention is needed for the protection of new varieties. Due to its typical life cycle, it is very easy to generate so called "look-a-likes" from protected cultivars by screening fertile single spore cultures. A consensus has been reached within the mushroom (breeding) industry to consider this method as the generation of essentially derived varieties as defined in plant breeding.

Association Mapping Reveals Genetic Loci Associated with Important Agronomic Traits in Lentinula edodes, Shiitake Mushroom

Association mapping is a robust approach for the detection of quantitative trait loci (QTLs). Here, by genotyping 297 genome-wide molecular markers of 89 Lentinula edodes cultivars in China, the genetic diversity, population structure and genetic loci associated with 11 agronomic traits were examined. A total of 873 alleles were detected in the tested strains with a mean of 2.939 alleles per locus, and the Shannon's information index was 0.734. Population structure analysis revealed two robustly differentiated groups among the Chinese L. edodes cultivars (F_ST = 0.247). Using the mixed linear model, a total of 43 markers were detected to be significantly associated with four traits. The number of markers associated with traits ranged from 9 to 26, and the phenotypic variations explained by each marker varied from 12.07% to 31.32%. Apart from five previously reported markers, the remaining 38 markers were newly reported here. Twenty-one markers were identified as simultaneously linked to two to four traits, and five markers were associated with the same traits in cultivation tests performed in two consecutive years. The 43 traits-associated markers were related to 97 genes, and 24 of them were related to 10 traits-associated markers detected in both years or identified previously, 13 of which had a >2-fold expression change between the mycelium and primordium stages. Our study has provided candidate markers for marker-assisted selection (MAS) and useful clues for understanding the genetic architecture of agronomic traits in the shiitake mushroom.

Multi-trait QTL analysis for agronomic and quality characters of Agaricus bisporus (button mushrooms)

The demand for button mushrooms of high quality is increasing. Superior button mushroom varieties require the combination of multiple traits to maximize productivity and quality. Very often these traits are correlated and should, therefore, be evaluated together rather than as single traits. In order to unravel the genetic architecture of multiple traits of Agaricus bisporus and the genetic correlations among traits, we have investigated a total of six agronomic and quality traits through multi-trait QTL analyses in a mixed-model. Traits were evaluated in three heterokaryon sets. Significant phenotypic correlations were observed among traits. For instance, earliness (ER) correlated to firmness (FM), cap color, and compost colonization, and FM correlated to scales (SC). QTLs of different traits located on the same chromosomes genetically explains the phenotypic correlations. QTL detected on chromosome 10 mainly affects three traits, i.e., ER, FM and SC. It explained 31.4 % phenotypic variation of SC on mushroom cap (heterokaryon Set 1), 14.9 % that of the FM (heterokaryon Set 3), and 14.2 % that of ER (heterokaryon Set 3). High value alleles from the wild parental line showed beneficial effects for several traits, suggesting that the wild germplasm is a valuable donor in terms of those traits. Due to the limitations of recombination pattern, we only made a start at understanding the genetic base for several agronomic and quality traits in button mushrooms.

Experimental evolution in fungi: An untapped resource

Historically, evolutionary biology has been considered an observational science. Examining populations and inferring evolutionary histories mold evolutionary theories. In contrast, laboratory evolution experiments make use of the amenability of traditional model organisms to study fundamental processes underlying evolution in real time in simple, but well-controlled, environments. With advances in high-throughput biology and next generation sequencing, it is now possible to propagate hundreds of parallel populations over thousands of generations and to quantify precisely the frequencies of various mutations over time. Experimental evolution combines the ability to simultaneously monitor replicate populations with the power to vary individual parameters to test specific evolutionary hypotheses, something that is impractical or infeasible in natural populations. Many labs are now conducting laboratory evolution experiments in nearly all model systems including viruses, bacteria, yeast, nematodes, and fruit flies. Among these systems, fungi occupy a unique niche: with a short generation time, small compact genomes, and sexual cycles, fungi are a particularly valuable and largely untapped resource for propelling future growth in the field of experimental evolution. Here, we describe the current state of fungal experimental evolution and why fungi are uniquely positioned to answer many of the outstanding questions in the field. We also review which fungal species are most well suited for experimental evolution.

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.

Genetic dissection of fruiting body-related traits using quantitative trait loci mapping in Lentinula edodes

To provide a better understanding of the genetic architecture of fruiting body formation of Lentinula edodes, quantitative trait loci (QTLs) mapping was employed to uncover the loci underlying seven fruiting body-related traits (FBRTs). An improved L. edodes genetic linkage map, comprising 572 markers on 12 linkage groups with a total map length of 983.7 cM, was constructed by integrating 82 genomic sequence-based insertion-deletion (InDel) markers into a previously published map. We then detected a total of 62 QTLs for seven target traits across two segregating testcross populations, with individual QTLs contributing 5.5 %-30.2 % of the phenotypic variation. Fifty-three out of the 62 QTLs were clustered in six QTL hotspots, suggesting the existence of main genomic regions regulating the morphological characteristics of fruiting bodies in L. edodes. A stable QTL hotspot on MLG2, containing QTLs for all investigated traits, was identified in both testcross populations. QTLs for related traits were frequently co-located on the linkage groups, demonstrating the genetic basis for phenotypic correlation of traits. Meta-QTL (mQTL) analysis was performed and identified 16 mQTLs with refined positions and narrow confidence intervals (CIs). Nine genes, including those encoding MAP kinase, blue-light photoreceptor, riboflavin-aldehyde-forming enzyme and cyclopropane-fatty-acyl-phospholipid synthase, and cytochrome P450s, were likely to be candidate genes controlling the shape of fruiting bodies. The study has improved our understanding of the genetic architecture of fruiting body formation in L. edodes. To our knowledge, this is the first genome-wide QTL detection of FBRTs in L. edodes. The improved genetic map, InDel markers and QTL hotspot regions revealed here will assist considerably in the conduct of future genetic and breeding studies of L. edodes.

Quantitative trait locus mapping of melanization in the plant pathogenic fungus Zymoseptoria tritici

Melanin plays an important role in virulence and antimicrobial resistance in several fungal pathogens. The wheat pathogen Zymoseptoria tritici is important worldwide, but little is known about the genetic architecture of pathogenicity, including the production of melanin. Because melanin production can exhibit complex inheritance, we used quantitative trait locus (QTL) mapping in two crosses to identify the underlying genes. Restriction site−associated DNA sequencing was used to genotype 263 (cross 1) and 261 (cross 2) progeny at ~8500 single-nucleotide polymorphisms and construct two dense linkage maps. We measured gray values, representing degrees of melanization, for single-spore colonies growing on Petri dishes by using a novel image-processing approach that enabled high-throughput phenotyping. Because melanin production can be affected by stress, each offspring was grown in two stressful environments and one control environment. We detected six significant QTL in cross 1 and nine in cross 2, with three QTL shared between the crosses. Different QTL were identified in different environments and at different colony ages. By obtaining complete genome sequences for the four parents and analyzing sequence variation in the QTL confidence intervals, we identified 16 candidate genes likely to affect melanization. One of these candidates was PKS1, a polyketide synthase gene known to play a role in the synthesis of dihydroxynaphthalene melanin. Three candidate quantitative trait nucleotides were identified in PKS1. Many of the other candidate genes were not previously associated with melanization.

Genome-wide patterns of segregation and linkage disequilibrium: the construction of a linkage genetic map of the poplar rust fungus Melampsora larici-populina

The poplar rust fungus Melampsora larici-populina causes significant yield reduction and severe economic losses in commercial poplar plantations. After several decades of breeding for qualitative resistance and subsequent breakdown of the released resistance genes, breeders now focus on quantitative resistance, perceived to be more durable. But quantitative resistance also can be challenged by an increase of aggressiveness in the pathogen. Thus, it is of primary importance to better understand the genetic architecture of aggressiveness traits. To this aim, our goal is to build a genetic linkage map for M. larici-populina in order to map quantitative trait loci related to aggressiveness. First, a large progeny of M. larici-populina was generated through selfing of the reference strain 98AG31 (which genome sequence is available) on larch plants, the alternate host of the poplar rust fungus. The progeny's meiotic origin was validated through a segregation analysis of 115 offspring with 14 polymorphic microsatellite markers, of which 12 segregated in the expected 1:2:1 Mendelian ratio. A microsatellite-based linkage disequilibrium analysis allowed us to identify one potential linkage group comprising two scaffolds. The whole genome of a subset of 47 offspring was resequenced using the Illumina HiSeq 2000 technology at a mean sequencing depth of 6X. The reads were mapped onto the reference genome of the parental strain and 144,566 SNPs were identified across the genome. Analysis of distribution and polymorphism of the SNPs along the genome led to the identification of 2580 recombination blocks. A second linkage disequilibrium analysis, using the recombination blocks as markers, allowed us to group 81 scaffolds into 23 potential linkage groups. These preliminary results showed that a high-density linkage map could be constructed by using high-quality SNPs based on low-coverage resequencing of a larger number of M. larici-populina offspring.

Constructing a new integrated genetic linkage map and mapping quantitative trait loci for vegetative mycelium growth rate in Lentinula edodes

The most saturated linkage map for Lentinula edodes to date was constructed based on a monokaryotic population of 146 single spore isolates (SSIs) using sequence-related amplified polymorphism (SRAP), target region amplification polymorphism (TRAP), insertion-deletion (InDel) markers, and the mating-type loci. Five hundred and twenty-four markers were located on 13 linkage groups (LGs). The map spanned a total length of 1006.1 cM, with an average marker spacing of 2.0 cM. Quantitative trait loci (QTLs) mapping was utilized to uncover the loci regulating and controlling the vegetative mycelium growth rate on various synthetic media, and complex medium for commercial cultivation of L. edodes. Two and 13 putative QTLs, identified respectively in the monokaryotic population and two testcross dikaryotic populations, were mapped on seven different LGs. Several vegetative mycelium growth rate-related QTLs uncovered here were clustered on LG4 (Qmgr1, Qdgr1, Qdgr2 and Qdgr9) and LG6 (Qdgr3, Qdgr4 and Qdgr5), implying the presence of main genomic areas responsible for growth rate regulation and control. The QTL hotspot region on LG4 was found to be in close proximity to the region containing the mating-type A (MAT-A) locus. Moreover, Qdgr2 on LG4 was detected on different media, contributing 8.07 %-23.71 % of the phenotypic variation. The present study provides essential information for QTL mapping and marker-assisted selection (MAS) in L. edodes.

Genetic variation and combining ability analysis of bruising sensitivity in Agaricus bisporus

Advanced button mushroom cultivars that are less sensitive to mechanical bruising are required by the mushroom industry, where automated harvesting still cannot be used for the fresh mushroom market. The genetic variation in bruising sensitivity (BS) of Agaricus bisporus was studied through an incomplete set of diallel crosses to get insight in the heritability of BS and the combining ability of the parental lines used and, in this way, to estimate their breeding value. To this end nineteen homokaryotic lines recovered from wild strains and cultivars were inter-crossed in a diallel scheme. Fifty-one successful hybrids were grown under controlled conditions, and the BS of these hybrids was assessed. BS was shown to be a trait with a very high heritability. The results also showed that brown hybrids were generally less sensitive to bruising than white hybrids. The diallel scheme allowed to estimate the general combining ability (GCA) for each homokaryotic parental line and to estimate the specific combining ability (SCA) of each hybrid. The line with the lowest GCA is seen as the most attractive donor for improving resistance to bruising. The line gave rise to hybrids sensitive to bruising having the highest GCA value. The highest negative SCA possibly indicates heterosis effects for resistance to bruising. This study provides a foundation for estimating breeding value of parental lines to further study the genetic factors underlying bruising sensitivity and other quality-related traits, and to select potential parental lines for further heterosis breeding. The approach of studying combining ability in a diallel scheme was used for the first time in button mushroom breeding.

Genome wide association identifies novel loci involved in fungal communication

Understanding how genomes encode complex cellular and organismal behaviors has become the outstanding challenge of modern genetics. Unlike classical screening methods, analysis of genetic variation that occurs naturally in wild populations can enable rapid, genome-scale mapping of genotype to phenotype with a medium-throughput experimental design. Here we describe the results of the first genome-wide association study (GWAS) used to identify novel loci underlying trait variation in a microbial eukaryote, harnessing wild isolates of the filamentous fungus Neurospora crassa. We genotyped each of a population of wild Louisiana strains at 1 million genetic loci genome-wide, and we used these genotypes to map genetic determinants of microbial communication. In N. crassa, germinated asexual spores (germlings) sense the presence of other germlings, grow toward them in a coordinated fashion, and fuse. We evaluated germlings of each strain for their ability to chemically sense, chemotropically seek, and undergo cell fusion, and we subjected these trait measurements to GWAS. This analysis identified one gene, NCU04379 (cse-1, encoding a homolog of a neuronal calcium sensor), at which inheritance was strongly associated with the efficiency of germling communication. Deletion of cse-1 significantly impaired germling communication and fusion, and two genes encoding predicted interaction partners of CSE1 were also required for the communication trait. Additionally, mining our association results for signaling and secretion genes with a potential role in germling communication, we validated six more previously unknown molecular players, including a secreted protease and two other genes whose deletion conferred a novel phenotype of increased communication and multi-germling fusion. Our results establish protein secretion as a linchpin of germling communication in N. crassa and shed light on the regulation of communication molecules in this fungus. Our study demonstrates the power of population-genetic analyses for the rapid identification of genes contributing to complex traits in microbial species.

Many phenotypes of interest are controlled by multiple loci, and in biological systems identifying determinants of such complex traits is challenging. Here, we genotyped 112 wild isolates of Neurospora crassa and used this resource to identify genes that mediate a fundamental but poorly-understood attribute of this filamentous fungus: the ability of germinating spores to sense each other at a distance, extend projections toward one another, and fuse. Inheritance at a secretion gene, cse-1, was associated strongly with germling communication across wild strains; this association was validated in experiments showing reduced communication in a cse-1 deletion strain. By testing interacting partners of CSE1, and by assessing additional secretion and signaling factors whose inheritance associated more modestly with germling communication in wild strains, we identified eight other novel determinants of this phenotype. Our population of genotyped wild isolates provides a flexible and powerful community resource for the rapid identification of any varying, complex phenotype in N. crassa. The success of our approach, which used a phenotyping scheme far more tractable than would be required in a screen of the entire N. crassa gene deletion collection, serves as a proof of concept for association studies of wild populations for any organism.

Linkage mapping identifies the sex determining region as a single locus in the Pennate diatom Seminavis robusta

The pennate diatom Seminavis robusta, characterized by an archetypical diatom life cycle including a heterothallic mating system, is emerging as a model system for studying the molecular regulation of the diatom cell and life cycle. One of its main advantages compared with other diatom model systems is that sexual crosses can be made routinely, offering unprecedented possibilities for forward genetics. To date, nothing is known about the genetic basis of sex determination in diatoms. Here, we report on the construction of mating type-specific linkage maps for S. robusta, and use them to identify a single locus sex determination system in this diatom. We identified 13 mating type plus and 15 mating type minus linkage groups obtained from the analysis of 463 AFLP markers segregating in a full-sib family, covering 963.7 and 972.2 cM, respectively. Five linkage group pairs could be identified as putative homologues. The mating type phenotype mapped as a monogenic trait, disclosing the mating type plus as the heterogametic sex. This study provides the first evidence for a genetic sex determining mechanism in a diatom.