Distribution and localization of microsatellites in the Perigord black truffle genome and identification of new molecular markers

Comparative plastome genomics and phylogenetic relationships of the genus Trollius

Trollius, a genus in the Ranunculaceae family, has significant medicinal and ornamental value. It is widely distributed in China with 16 different species accepted. However, due to the lack of enough samples and information sites, the molecular phylogenetic relationships of Trollius have been unresolved till now. Here we sequenced, assembled and annotated the plastomes of 16 Trollius species to investigate their genomic characteristics, inverted repeat (IR) boundaries, sequence repeats, and hypervariable loci. In addition, the phylogenetic relationships of this genus was reconstructed based on the whole plastomes and the protein-coding sequences data-sets. The plastomes of Trollius ranged between 159,597 bp and 160,202 bp in length, and contained 113 unique genes, including 79 protein coding, 30 tRNA, and 4 rRNA. The IR boundaries were relatively conserved within the genus Trollius. 959 simple sequence repeats and 657 long sequence repeats were detected in the Trollius plastomes. We identified 12 highly polymorphic loci (Pi > 0.0115) that can be used as plastid markers in molecular identification and phylogenetic investigation of the genus. Besides, Trollius was a monophyletic group with the earliest divergence clade being Trollius lilacinus Bunge, and the remaining species were divided into two strongly-supported clades. The phylogeny in our study supported the traditional classification systems based on the color of sepal, but not the previous classification system based on the types and relative lengths of the nectaries, and distribution. The genomic resources provided in our study can be used in the taxonomy of the genus Trollius, promoting the development and utilization of this genus.

Plastid phylogenomics and species discrimination in the "Chinese" clade of Roscoea (Zingiberaceae)

Roscoea is an alpine or subalpine genus from the pan-tropical family Zingiberaceae, which consists of two disjunct groups in geography, namely the "Chinese" clade and the "Himalayan" clade. Despite extensive research on the genus, Roscoea species remain poorly defined and relationships between these species are not well resolved. In this study, we used plastid genomes of nine species and one variety to resolve phylogenetic relationships within the "Chinese" clade of Roscoea and as DNA super barcodes for species discrimination. We found that Roscoea plastid genomes ranged in length from 163,063 to 163,796 bp, and encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, four rRNA genes. In addition, expansion and contraction of the IR regions showed obvious infraspecific conservatism and interspecific differentiation. Plastid phylogenomics revealed that species belonging to the "Chinese" clade of Roscoea can be divided into four distinct subclades. Furthermore, our analysis supported the independence of R. cautleoides var. pubescens, the recovery of Roscoea pubescens Z.Y. Zhu, and a close relationship between R. humeana and R. cautloides. When we used the plastid genome as a super barcode, we found that it possessed strong discriminatory power (90%) with high support values. Intergenic regions provided similar resolution, which was much better than that of protein-coding regions, hypervariable regions, and DNA universal barcodes. However, plastid genomes could not completely resolve Roscoea phylogeny or definitively discriminate species. These limitations are likely related to the complex history of Roscoea speciation, poorly defined species within the genus, and the maternal inheritance of plastid genomes.

Analysis of the Genome Sequence of Strain GiC-126 of Gloeostereum incarnatum with Genetic Linkage Map

Gloeostereum incarnatum has edible and medicinal value and was first cultivated and domesticated in China. We sequenced the G. incarnatum monokaryotic strain GiC-126 on an Illumina HiSeq X Ten system and obtained a 34.52-Mb genome assembly sequence that encoded 16,895 predicted genes. We combined the GiC-126 genome with the published genome of G. incarnatum strain CCMJ2665 to construct a genetic linkage map (GiC-126 genome) that had 10 linkage groups (LGs), and the 15 assembly sequences of CCMJ2665 were integrated into 8 LGs. We identified 1912 simple sequence repeat (SSR) loci and detected 700 genes containing 768 SSRs in the genome; 65 and 100 of them were annotated with gene ontology (GO) terms and KEGG pathways, respectively. Carbohydrate-active enzymes (CAZymes) were identified in 20 fungal genomes and annotated; among them, 144 CAZymes were annotated in the GiC-126 genome. The A mating-type locus (MAT-A) of G. incarnatum was located on scaffold885 at 38.9 cM of LG1 and was flanked by two homeodomain (HD1) genes, mip and beta-fg. Fourteen segregation distortion markers were detected in the genetic linkage map, all of which were skewed toward the parent GiC-126. They formed three segregation distortion regions (SDR1-SDR3), and 22 predictive genes were found in scaffold1920 where three segregation distortion markers were located in SDR1. In this study, we corrected and updated the genomic information of G. incarnatum. Our results will provide a theoretical basis for fine gene mapping, functional gene cloning, and genetic breeding the follow-up of G. incarnatum.

Genome-wide characterization of microsatellite DNA in fishes: survey and analysis of their abundance and frequency in genome-specific regions

BACKGROUND

Microsatellite repeats are ubiquitous in organism genomes and play an important role in the chromatin organization, regulation of gene activity, recombination and DNA replication. Although microsatellite distribution patterns have been studied in most phylogenetic lineages, they are unclear in fish species.

RESULTS

Here, we present the first systematic examination of microsatellite distribution in coding and non-coding regions of 14 fish genomes. Our study showed that the number and type of microsatellites displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation and DNA replication slippage theories alone were insufficient to explain the distribution patterns. Our results showed that microsatellites are dominant in non-coding regions. The total number of microsatellites ranged from 78,378 to 1,012,084, and the relative density varied from 4925.76 bp/Mb to 25,401.97 bp/Mb. Overall, (A + T)-rich repeats were dominant. The dependence of repeat abundance on the length of the repeated unit (1-6 nt) showed a great similarity decrease, whereas more tri-nucleotide repeats were found in exonic regions than tetra-nucleotide repeats of most species. Moreover, the incidence of different repeated types appeared species- and genomic-specific. These results highlight potential mechanisms for maintaining microsatellite distribution, such as selective forces and mismatch repair systems.

CONCLUSIONS

Our data could be beneficial for the studies of genome evolution and microsatellite DNA evolutionary dynamics, and facilitate the exploration of microsatellites structural, function, composition mode and molecular markers development in these species.

Genome-Wide Characterization and Comparative Analyses of Simple Sequence Repeats among Four Miniature Pig Breeds

Simple Summary

Simple sequence repeats (SSRs) are present at high densities in regulatory elements, suggesting that they may affect gene function and phenotypic traits. Therefore, SSRs can be exploited in marker-assisted selection. In addition, they can be widely used as molecular markers to study genetic diversity, population structure, and evolution. While SSRs have been widely studied in many mammalian species, very little research has focused on genome-wide SSRs of miniature pigs, a small but special group of pigs that express the dwarf phenotype. Based on the SSR-enriched library building and sequencing, about 30,000 novel polymorphic SSRs for four miniature pig breeds were mapped to the Duroc pig reference genome. The four miniature pig breeds had different numbers and types of SSRs and distributions of repeat units. There were 2518 polymorphic SSRs in the intron or exon regions that were common to all four breeds and functional analyses revealed 17 genes that were associated with body size and other genes that were associated with growth and development. In conclusion, the SSRs detected in the miniature pigs in this study may provide useful genetic markers for the selection of farm animals and the polymorphic SSRs provide valuable insights into the determination of mature body size, as well as the immunity, growth and development of animals.

Abstract

Simple sequence repeats (SSRs) are commonly used as molecular markers in research on genetic diversity and discrimination among taxa or breeds because polymorphisms in these regions contribute to gene function and phenotypically important traits. In this study, we investigated genome-wide characteristics, repeat units, and polymorphisms of SSRs using sequencing data from SSR-enriched libraries created from Wuzhishan (WZS), Bama (BM), inbred Luchuan (LC) and Zangxiang (ZX) miniature pig breeds. The numbers and types of SSRs, distributions of repeat units and polymorphic SSRs varied among the four breeds. Compared to the Duroc pig reference genome, 2518 polymorphic SSRs were unique and common to all four breeds and functional annotation revealed that they may affect the coding and regulatory regions of genes. Several examples, such as FGF23, MYF6, IGF1R, and LEPROT, are associated with growth and development in pigs. Three of the polymorphic SSRs were selected to confirm the polymorphism and the corresponding alleles through fluorescence polymerase chain reaction (PCR) and capillary electrophoresis. Together, this study provides useful insights into the discovery, characteristics and distribution of SSRs in four pig breeds. The polymorphic SSRs, especially those common and unique to all four pig breeds, might affect associated genes and play important roles in growth and development.

Two ectomycorrhizal truffles, Tuber melanosporum and T. aestivum, endophytically colonise roots of non-ectomycorrhizal plants in natural environments

Serendipitous findings and studies on Tuber species suggest that some ectomycorrhizal fungi, beyond their complex interaction with ectomycorrhizal hosts, also colonise roots of nonectomycorrhizal plants in a loose way called endophytism. Here, we investigate endophytism of T. melanosporum and T. aestivum. We visualised endophytic T. melanosporum hyphae by fluorescent in situ hybridisation on nonectomycorrhizal plants. For the two Tuber species, microsatellite genotyping investigated the endophytic presence of the individuals whose mating produced nearby ascocarps. We quantified the expression of four T. aestivum genes in roots of endophyted, non-ectomycorrhizal plants. Tuber melanosporum hyphae colonised the apoplast of healthy roots, confirming endophytism. Endophytic Tuber melanosporum and T. aestivum contributed to nearby ascocarps, but only as maternal parents (forming the flesh). Paternal individuals (giving only genes found in meiotic spores of ascocarps) were not detected. Gene expression of T. aestivum in non-ectomycorrhizal plants confirmed a living status. Tuber species, and likely other ectomycorrhizal fungi found in nonectomycorrhizal plant roots in this study, can be root endophytes. This is relevant for the ecology (brûlé formation) and commercial production of truffles. Evolutionarily speaking, endophytism may be an ancestral trait in some ectomycorrhizal fungi that evolved from root endophytes.

Tuberomics: a molecular profiling for the adaption of edible fungi (Tuber magnatum Pico) to different natural environments

BACKGROUND

Truffles are symbiotic fungi that develop underground in association with plant roots, forming ectomycorrhizae. They are primarily known for the organoleptic qualities of their hypogeous fruiting bodies. Primarily, Tuber magnatum Pico is a greatly appreciated truffle species mainly distributed in Italy and Balkans. Its price and features are mostly depending on its geographical origin. However, the genetic variation within T. magnatum has been only partially investigated as well as its adaptation to several environments.

RESULTS

Here, we applied an integrated omic strategy to T. magnatum fruiting bodies collected during several seasons from three different areas located in the North, Center and South of Italy, with the aim to distinguish them according to molecular and biochemical traits and to verify the impact of several environments on these properties. With the proteomic approach based on two-dimensional electrophoresis (2-DE) followed by mass spectrometry, we were able to identify proteins specifically linked to the sample origin. We further associated the proteomic results to an RNA-seq profiling, which confirmed the possibility to differentiate samples according to their source and provided a basis for the detailed analysis of genes involved in sulfur metabolism. Finally, geographical specificities were associated with the set of volatile compounds produced by the fruiting bodies, as quantitatively and qualitatively determined through proton transfer reaction-mass spectrometry (PTR-MS) and gas-chromatography-mass spectrometry (GC-MS). In particular, a partial least squares-discriminant analysis (PLS-DA) model built from the latter data was able to return high confidence predictions of sample source.

CONCLUSIONS

Results provide a characterization of white fruiting bodies by a wide range of different molecules, suggesting the role for specific compounds in the responses and adaptation to distinct environments.

Genetic Structure and Phylogeography of Tuber magnatum Populations

The ectomycorrhizal fungus Tuber magnatum produces the white truffle appreciated worldwide for its unique aroma. With respect to other Tuber spp. of economic interest, T. magnatum presents a narrower geographical range. This species has, in fact, long been considered endemic to Italy. However, over the last few decades several reports have documented the presence of white truffles in different Mediterranean countries and in particular in various areas of south-east Europe. In this study, samples from several Pannonian and Balkan countries such as Hungary, Serbia, Romania, Bulgaria and Greece have been collected and genotyped with microsatellite markers and the data merged with those available for Italian populations. Our objectives were to test whether Italian and south-east European populations are differentiated and to evaluate the genetic diversity of T. magnatum all over its distributional range. We show the genetic structure of T. magnatum populations with the differentiation of four main groups: northern Italy, central-northern Italy, southern Italy and the Balkan/Pannonian region. The present study allowed us to refine the evolutionary history of T. magnatum and track the possible post-glacial expansion route of this species. The assessment of T. magnatum’s genetic structure is not only of scientific relevance, but it is also important for the conservation and market traceability of this prestigious fungus.

Phylogenetic analysis of Fritillaria cirrhosa D. Don and its closely related species based on complete chloroplast genomes

Fritillaria cirrhosa D. Don, whose bulb is used in a well-known traditional Chinese medicine to relieve cough and eliminate phlegm, is one of the most important medicinal plants of Fritillaria L. The species is widely distributed among the alpine regions in southwestern China and possesses complex morphological variations in different distributions. A series of newly related species were reported, based on obscure morphological differences. As a result, F. cirrhosa and its closely related species constitute a taxonomically complex group. However, it is difficult to accurately identify these species and reveal their phylogenetic relationships using traditional taxonomy. Molecular markers and gene fragments have been adopted but they are not able to afford sufficient phylogenetic resolution in the genus. Here, we report the complete chloroplast genome sequences of F. cirrhosa and its closely related species using next generation sequencing (NGS) technology. Eight plastid genomes ranged from 151,058 bp to 152,064 bp in length and consisted of 115 genes. Gene content, gene order, GC content, and IR/SC boundary structures were highly similar among these genomes. SSRs and five large repeat sequences were identified and the total number of them ranged from 73 to 79 and 63 to 75, respectively. Six highly divergent regions were successfully identified that could be used as potential genetic markers of Fritillaria. Phylogenetic analyses revealed that eight Fritillaria species were clustered into three clades with strong supports and F. cirrhosa was closely related to F. przewalskii and F. sinica. Overall, this study indicated that the complete chloroplast genome sequence was an efficient tool for identifying species in taxonomically complex groups and exploring their phylogenetic relationships.

Ascoma genotyping and mating type analyses of mycorrhizas and soil mycelia of Tuber borchii in a truffle orchard established by mycelial inoculated plants

Tuber borchii (the Bianchetto truffle) is a heterothallic Ascomycete living in symbiotic association with trees and shrubs. Maternal and paternal genotype dynamics have already been studied for the black truffles Tuber melanosporum and Tuber aestivum but not yet for T. borchii. In this study, we analysed maternal and paternal genotypes in the first truffle orchard realized with plants inoculated with five different T. borchii mycelia. Our aims were to test the persistence of the inoculated mycelia, if maternal and/or paternal genotypes correspond to inoculated mycelia and to assess the hermaphroditism of T. borchii. The mating type of each isolate as well as those of mycorrhizas, ascomata and extraradical soil mycelia was determined. Moreover, simple sequence repeat (SSR) profiles of maternal and paternal genotypes were assessed in 18 fruiting bodies to investigate the sexual behaviour of this truffle. The maternal genotypes of the fruiting bodies corresponded to those of the inoculated mycelia with only two exceptions. This confirmed that the inoculated mycelia persisted 9 years after plantation. As regards paternal partner, only two had the same genotype as those of the inoculated mycelia, suggesting hermaphroditism. Most of the new paternal genotypes originated from a recombination of those of inoculated mycelia.

Hybridization, characterization and transferability of SSRs in the genus Morchella

Recently, Morchella importuna, M. sextelata, M. eximia, M. exuberans, Mel-13, and Mel-21 have been successfully cultivated in China and planting areas rapidly expanded because of their economic importance. Effective molecular markers are urgently needed for accurately identifying morel cultivars. Microsatellites are widely utilized for strain authentication in many fungal groups. To our knowledge, for the first time we characterized the distribution of microsatellites (simple sequence repeats, SSRs) in the M. importuna genome with 12902 SSRs and reported the first set of SSRs developed for Morchella species. Mono-nucleotides (66.2 %) were the most frequent motifs, followed by tri- (15.4 %), di- (12.1 %), tetra- (3.7 %), penta- (1.3 %) and hexa-nucleotides (1.3 %). We tested the cross-species amplification of 180 SSRs on 24 samples from the six species and high cross-species transferability of SSRs (87.7 %) was found. Among twenty-two microsatellites selected for genetic diversity analysis on 127 samples from the six species, fifteen to twenty polymorphic loci were identified in M. importuna, M. sextelata, M. eximia, M. exuberans, Mel-13 and Mel-21. Interspecific hybridization events were detected among morel species, indicating the potential application of morel crossbreeding. Ninety-one cultivated samples were characterized as new cultivars with different genotypes, but cultivar names used for these by farmers was confusing, with misnaming, synonyms and homonyms. Our results are not only helpful for cultivar identification and morel breeding programs in China, but also provide molecular tools for genetic studies in morels.

Ribosomal DNA polymorphisms reveal genetic structure and a phylogeographic pattern in the Burgundy truffle Tuber aestivum Vittad

Ectomycorrhizal ascomycetes belonging to the genus Tuber produce edible fruiting bodies known as truffles. Tuber aestivum, in particular, is a fungus appreciated worldwide and has a natural distribution throughout Europe. Most of the molecular studies conducted on this species have been focused on the question as to whether or not T. aestivum and the morphologically similar T. uncinatum are conspecific. Conversely, only a handful of studies have assessed the level and distribution of genetic diversity and occurrence of phylogeographic patterns in this species. Here, we analyzed the genetic diversity of T. aestivum over a wide geographic range, performing an extensive sampling of specimens from Turkey, which is novel, to the best of our knowledge. We compared the internal transcribed spacer (ITS) profiles of 45 samples collected in different Turkish areas with those of 144 samples from all over Europe. We identified 63 haplotypes, 32 of which were exclusively present in Turkey. The majority of these haplotyes were also population specific. Haplotype network analysis and statistical tests highlighted the presence of a genetic structure and phylogeographic pattern, with three spatially distinct genetic clusters (northeastern Europe, southern Europe, and Turkey), with Turkey representing a diversity hotspot. Based on these results, we hypothesize the presence of glacial refugia for T. aestivum in Turkey, whereas European populations likely experienced a population bottleneck. The possible occurrence of cryptic species among Turkish T. aestivum samples also emerged. Our results are of practical relevance for the marketing of T. aestivum truffles and mycorrhizal seedlings and the preservation of the biodiversity of this species.

De novo assembly and comparative transcriptome analysis: novel insights into terpenoid biosynthesis in Chamaemelum nobile L

Analysis of terpenoids content, transcriptome from Chamaemelum nobile showed that the content of the terpenoids in the roots was the highest and key genes involved in the terpenoids synthesis pathway were identified. Chamaemelum nobile is a widely used herbaceous medicinal plant rich in volatile oils, mainly composed of terpenoids. It is widely used in food, cosmetics, medicine, and other fields. In this study, we analyzed the transcriptome and the content and chemical composition of the terpenoids in different organs of C. nobile. Gas chromatography-mass spectrometry analysis showed that the total content of the terpenoids among C. nobile organs was highest in the roots, followed by the flowers. Illumina HiSeq 2500 high-throughput sequencing technology was used to sequence the transcripts of roots, stems, leaves, and flowers of C. nobile. We obtained 139,757 unigenes using the Trinity software assembly. A total of 887 unigenes were annotated to secondary metabolism. In total, 55,711 differentially expressed genes were screened among different organs of C. nobile. We identified 16 candidate genes that may be involved in the terpenoid biosynthesis from C. nobile and analyzed their expression patterns using real-time PCR. Results showed that the expression pattern of these genes was tissue-specific and had significant differential expression levels in different organs of C. nobile. Among these genes, 13 were expressed in roots with the highest levels. Furthermore, the transcript levels of these 13 genes were positively correlated with the content of α-pinene, β-phellandrene, 1,8-cineole, camphor, α-terpineol, carvacrol, (E,E)-farnesol and chamazulene, suggesting that these 13 genes may be involved in the regulation of the synthesis of the volatile terpenoids. These results laid the foundation for the subsequent improvement of C. nobile quality through genetic engineering.

The Uneven Distribution of Mating Type Genes in Natural and Cultivated Truffle Orchards Contributes to the Fructification of Tuber indicum

The aim of this study was to investigate the pattern of distribution of mating type (MAT) genes of Tuber indicum in ectomycorhizosphere soils from natural T. indicum-producing areas and cultivated truffle orchards and ascocarp samples from different regions. Quantitative real-time PCR and multiplex PCR were used to weight the copy numbers of MAT1-1-1 and MAT1-2-1 in natural truffle soils and cultivated orchard soils. The effect of limestone on the pattern of truffle MAT genes and the correlation between soil properties and the proportion of MAT genes were also assessed. These results indicated that an uneven and nonrandom distribution of MAT genes was common in truffle-producing areas, cultivated truffle orchards, and ascocarps gleba. The competition between the two mating type genes and the expansion of unbalanced distribution was found to be closely related to truffle fructification. Limestone treatments failed to alter the proportion of the two mating type genes in the soil. The content of available phosphorus in soil was significantly correlated with the value of MAT1-1-1/MAT1-2-1 in cultivated and natural ectomycorhizosphere soils. The application of real-time quantitative PCR can provide reference for monitoring the dynamic changes of mating type genes in soil. This study investigates the distributional pattern of T. indicum MAT genes in the ectomycorhizosphere soil and ascocarp gleba from different regions, which may provide a foundation for the cultivation of T. indicum.

Distinct patterns of simple sequence repeats and GC distribution in intragenic and intergenic regions of primate genomes

As the first systematic examination of simple sequence repeats (SSRs) and guanine-cytosine (GC) distribution in intragenic and intergenic regions of ten primates, our study showed that SSRs and GC displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation. Our results suggest that the majority of SSRs are distributed in non-coding regions, such as the introns, TEs, and intergenic regions. In these primates, trinucleotide perfect (P) SSRs were the most abundant repeats type in the 5'UTRs and CDSs, whereas, mononucleotide P-SSRs were the most in the intron, 3'UTRs, TEs, and intergenic regions. The GC-contents varied greatly among different intragenic and intergenic regions: 5'UTRs > CDSs > 3'UTRs > TEs > introns > intergenic regions, and high GC-content was frequently distributed in exon-rich regions. Our results also showed that in the same intragenic and intergenic regions, the distribution of GC-contents were great similarity in the different primates. Tri- and hexanucleotide P-SSRs had the most GC-contents in the 5'UTRs and CDSs, whereas mononucleotide P-SSRs had the least GC-contents in the six genomic regions of these primates. The most frequent motifs for different length varied obviously with the different genomic regions.

Survey and analysis of simple sequence repeats in the Ustilaginoidea virens genome and the development of microsatellite markers

Ustilaginoidea virens is the causal agent of rice false smut, causing quantitative and qualitative losses in rice industry. However, the development and application of simple sequence repeat (SSR) markers for genetic diversity studies in U. virens were limited. This study is the first to perform large-scale development of SSR markers of this pathogen at the genome level, to (1) compare these SSR markers with those of other fungi, (2) analyze the pattern of the SSRs, and (3) obtain more informative genetic markers. U. virens is rich in SSRs, and 13,778 SSRs were identified with a relative abundance of 349.7SSRs/Mb. The most common motifs in the genome or in noncoding regions were mononucleotides, whereas trinucleotides in coding sequences. A total of 6 out of 127 primers were randomly selected to be used to analyze 115 isolates, and these 6 primers showed high polymorphism in U. virens. This study may serve as an important resource for molecular genetic studies in U. virens.

SSR-based identification of genetic groups within European populations of Tuber aestivum Vittad

Tuber species are ectomycorrhizal ascomycetes establishing relationships with different host trees and forming hypogeous fruiting bodies known as truffles. Among Tuber species, Tuber aestivum Vittad. has a wide distributional range being found naturally all over Europe. Here, we performed large-scale population genetic analyses in T. aestivum to (i) investigate its genetic diversity at the European scale, (ii) characterize its genetic structure and test for the presence of ecotypes and (iii) shed light into its demographic history. To reach these goals, 230 ascocarps from different populations were genotyped using 15 polymorphic simple sequence repeat markers. We identified 181 multilocus genotypes and four genetic groups which did not show a clear geographical separation; although, one of them was present exclusively in Southeast France, Italy and Spain. Fixation index values between pairs of genetic groups were generally high and ranged from 0.29 to 0.45. A significant deficit of heterozygosity indicated a population expansion instead of a recent population bottleneck, suggesting that T. aestivum is not endangered in Europe, not even in Mediterranean regions. Our study based on a large-scale population genetic analysis suggests that genetically distinct populations and likely ecotypes within T. aestivum are present. In turn, this study paves the way to future investigations aimed at addressing the biological and/or ecological factors that have concurred in shaping the population genetic structure of this species. Present results should also have implications for the truffle market since defining genetic markers are now possible at least for some specific T. aestivum genetic groups.

Selection processes in simple sequence repeats suggest a correlation with their genomic location: insights from a fungal model system

BACKGROUND

Adaptive processes shape the evolution of genomes and the diverse functions of different genomic regions are likely to have an impact on the trajectory and outcome of this evolution. The main underlying hypothesis of this study is that the evolution of Simple Sequence Repeats (SSRs) is correlated with the evolution of the genomic region in which they are located, resulting in differences of motif size, number of repeats, and levels of polymorphisms. These differences should be clearly detectable when analyzing the frequency and type of SSRs within the genome of a species, when studying populations within a species, and when comparing closely related sister taxa. By coupling a genome-wide SSR survey in the genome of the plant pathogenic fungus Heterobasidion irregulare with an analysis of intra- and interspecific variability of 39 SSR markers in five populations of the two sibling species H. irregulare and H. annosum, we investigated mechanisms of evolution of SSRs.

RESULTS

Results showed a clear dominance of trirepeats and a selection against other repeat number, i.e. di- and tetranucleotides, both in regions inside Open Reading Frames (ORFs) and upstream 5' untranslated region (5'UTR). Locus per locus AMOVA showed SSRs both inside ORFs and upstream 5'UTR were more conserved within species compared to SSRs in other genomic regions, suggesting their evolution is constrained by the functions of the regions they are in. Principal coordinates analysis (PCoA) indicated that even if SSRs inside ORFs were less polymorphic than those in intergenic regions, they were more powerful in differentiating species. These findings indicate SSRs evolution undergoes a directional selection pressure comparable to that of the ORFs they interrupt and to that of regions involved in regulatory functions.

CONCLUSIONS

Our work linked the variation and the type of SSRs with regions upstream 5'UTR, putatively harbouring regulatory elements, and shows that the evolution of SSRs might be affected by their location in the genome. Additionally, this study provides a first glimpse on a possible molecular basis for fast adaptation to the environment mediated by SSRs.

Expanding genomics of mycorrhizal symbiosis

The mycorrhizal symbiosis between soil fungi and plant roots is a ubiquitous mutualism that plays key roles in plant nutrition, soil health, and carbon cycling. The symbiosis evolved repeatedly and independently as multiple morphotypes [e.g., arbuscular mycorrhizae (AM), ectomycorrhizal (ECM)] in multiple fungal clades (e.g., phyla Glomeromycota, Ascomycota, Basidiomycota). The accessibility and cultivability of many mycorrhizal partners make them ideal models for symbiosis studies. Alongside molecular, physiological, and ecological investigations, sequencing led to the first three mycorrhizal fungal genomes, representing two morphotypes and three phyla. The genome of the ECM basidiomycete Laccaria bicolor showed that the mycorrhizal lifestyle can evolve through loss of plant cell wall-degrading enzymes (PCWDEs) and expansion of lineage-specific gene families such as short secreted protein (SSP) effectors. The genome of the ECM ascomycete Tuber melanosporum showed that the ECM type can evolve without expansion of families as in Laccaria, and thus a different set of symbiosis genes. The genome of the AM glomeromycete Rhizophagus irregularis showed that despite enormous phylogenetic distance and morphological difference from the other two fungi, symbiosis can involve similar solutions as symbiosis-induced SSPs and loss of PCWDEs. The three genomes provide a solid base for addressing fundamental questions about the nature and role of a vital mutualism.

Genome wide characterization of short tandem repeat markers in sweet orange (Citrus sinensis)

Sweet orange (Citrus sinensis) is one of the major cultivated and most-consumed citrus species. With the goal of enhancing the genomic resources in citrus, we surveyed, developed and characterized microsatellite markers in the ≈347 Mb sequence assembly of the sweet orange genome. A total of 50,846 SSRs were identified with a frequency of 146.4 SSRs/Mbp. Dinucleotide repeats are the most frequent repeat class and the highest density of SSRs was found in chromosome 4. SSRs are non-randomly distributed in the genome and most of the SSRs (62.02%) are located in the intergenic regions. We found that AT-rich SSRs are more frequent than GC-rich SSRs. A total number of 21,248 SSR primers were successfully developed, which represents 89 SSR markers per Mb of the genome. A subset of 950 developed SSR primer pairs were synthesized and tested by wet lab experiments on a set of 16 citrus accessions. In total we identified 534 (56.21%) polymorphic SSR markers that will be useful in citrus improvement. The number of amplified alleles ranges from 2 to 12 with an average of 4 alleles per marker and an average PIC value of 0.75. The newly developed sweet orange primer sequences, their in silico PCR products, exact position in the genome assembly and putative function are made publicly available. We present the largest number of SSR markers ever developed for a citrus species. Almost two thirds of the markers are transferable to 16 citrus relatives and may be used for constructing a high density linkage map. In addition, they are valuable for marker-assisted selection studies, population structure analyses and comparative genomic studies of C. sinensis with other citrus related species. Altogether, these markers provide a significant contribution to the citrus research community.

Analysis of simple sequence repeats in the Gaeumannomyces graminis var. tritici genome and the development of microsatellite markers

Understanding the genetic structure of Gaeumannomyces graminis var. tritici is essential for the establishment of efficient disease control strategies. It is becoming clear that microsatellites, or simple sequence repeats (SSRs), play an important role in genome organization and phenotypic diversity, and are a large source of genetic markers for population genetics and meiotic maps. In this study, we examined the G. graminis var. tritici genome (1) to analyze its pattern of SSRs, (2) to compare it with other plant pathogenic filamentous fungi, such as Magnaporthe oryzae and M. poae, and (3) to identify new polymorphic SSR markers for genetic diversity. The G. graminis var. tritici genome was rich in SSRs; a total 13,650 SSRs have been identified with mononucleotides being the most common motifs. In coding regions, the densities of tri- and hexanucleotides were significantly higher than in noncoding regions. The di-, tri-, tetra, penta, and hexanucleotide repeats in the G. graminis var. tritici genome were more abundant than the same repeats in M. oryzae and M. poae. From 115 devised primers, 39 SSRs are polymorphic with G. graminis var. tritici isolates, and 8 primers were randomly selected to analyze 116 isolates from China. The number of alleles varied from 2 to 7 and the expected heterozygosity (He) from 0.499 to 0.837. In conclusion, SSRs developed in this study were highly polymorphic, and our analysis indicated that G. graminis var. tritici is a species with high genetic diversity. The results provide a pioneering report for several applications, such as the assessment of population structure and genetic diversity of G. graminis var. tritici.

Impact of the competition between mating types on the cultivation of Tuber melanosporum: Romeo and Juliet and the matter of space and time

Major breakthroughs in our understanding of the life cycles of the symbiotic ascomycetes belonging to the genus Tuber have occurred over the last several years. A number of Tuber species produce edible fruiting bodies, known as truffles, that are marketed worldwide. A better understanding of the basic biological characteristics of Tuber spp. is likely to have tremendous practical relevance for their cultivation. Tuber melanosporum produces the most valuable black truffles and its genome has been recently sequenced. This species is now serving as a model for studying the biology of truffles. Here, we review recent progress in the understanding of sexual reproduction modalities in T. melanosporum. The practical relevance of these findings is outlined. In particular, the discoveries that T. melanosporum is heterothallic and that strains of different mating types compete to persist on the roots of host plants suggest that the spatial and temporal distributional patterns of strains of different mating types are key determinants of truffle fructification. The spatial segregation of the two mating types in areas where T. melanosporum occurs likely limits truffle production. Thus, host plant inoculation techniques and agronomic practices that might be pursued to manage T. melanosporum orchards with a balanced presence of the two mating partners are described.

Population genetics of the westernmost distribution of the glaciations-surviving black truffle Tuber melanosporum

The black truffle (Tuber melanosporum Vittad.) is an important natural resource due to its relevance as a delicacy in gastronomy. Different aspects of this hypogeous fungus species have been studied, including population genetics of French and Italian distribution ranges. Although those studies include some Spanish populations, this is the first time that the genetic diversity and genetic structure of the wide geographical range of the natural Spanish populations have been analysed. To achieve this goal, 23 natural populations were sampled across the Spanish geographical distribution. ISSR technique demonstrated its reliability and capability to detect high levels of polymorphism in the species. Studied populations showed high levels of genetic diversity (h N  = 0.393, h S  = 0.678, Hs = 0.418), indicating a non threatened genetic conservation status. These high levels may be a consequence of the wide distribution range of the species, of its spore dispersion by animals, and by its evolutionary history. AMOVA analysis showed a high degree of genetic structure among populations (47.89%) and other partitions as geographical ranges. Bayesian genetic structure analyses differentiated two main Spanish groups separated by the Iberian Mountain System, and showed the genetic uniqueness of some populations. Our results suggest the survival of some of these populations during the last glaciation, the Spanish southern distribution range perhaps surviving as had occurred in France and Italy, but it is also likely that specific northern areas may have acted as a refugia for the later dispersion to other calcareous areas in the Iberian Peninsula and probably France.

Microsatellites in the genome of the edible mushroom, Volvariella volvacea

Using bioinformatics software and database, we have characterized the microsatellite pattern in the V. volvacea genome and compared it with microsatellite patterns found in the genomes of four other edible fungi: Coprinopsis cinerea, Schizophyllum commune, Agaricus bisporus, and Pleurotus ostreatus. A total of 1346 microsatellites have been identified, with mono-nucleotides being the most frequent motif. The relative abundance of microsatellites was lower in coding regions with 21 No./Mb. However, the microsatellites in the V. volvacea gene models showed a greater tendency to be located in the CDS regions. There was also a higher preponderance of trinucleotide repeats, especially in the kinase genes, which implied a possible role in phenotypic variation. Among the five fungal genomes, microsatellite abundance appeared to be unrelated to genome size. Furthermore, the short motifs (mono- to tri-nucleotides) outnumbered other categories although these differed in proportion. Data analysis indicated a possible relationship between the most frequent microsatellite types and the genetic distance between the five fungal genomes.

Fine-scale spatial genetic structure of the black truffle (Tuber melanosporum) investigated with neutral microsatellites and functional mating type genes

The genetic structure of ectomycorrhizal (ECM) fungal populations results from both vegetative and sexual propagation. In this study, we have analysed the spatial genetic structure of Tuber melanosporum populations, a heterothallic ascomycete that produces edible fruit bodies. Ectomycorrhizas from oaks and hazels from two orchards were mapped and genotyped using simple sequence repeat markers and the mating type locus. The distribution of the two T. melanosporum mating types was also monitored in the soil. In one orchard, the genetic profiles of the ascocarps were compared with those of the underlying mycorrhizas. A pronounced spatial genetic structure was found. The maximum genet sizes were 2.35 and 4.70 m in the two orchards, with most manifesting a size < 1 m. Few genets persisted throughout two seasons. A nonrandom distribution pattern of the T. melanosporum was observed, resulting in field patches colonized by genets that shared the same mating types. Our findings suggest that competition occurs between genets and provide basic information on T. melanosporum propagation patterns that are relevant for the management of productive truffle orchards.

Analysis of the Dendrobium officinale transcriptome reveals putative alkaloid biosynthetic genes and genetic markers

Dendrobium officinale Kimura et Migo (Orchidaceae) is a traditional Chinese medicinal plant. The stem contains an alkaloid that is the primary bioactive component. However, the details of alkaloid biosynthesis have not been effectively explored because of the limited number of expressed sequence tags (ESTs) available in GenBank. In this study, we analyzed RNA isolated from the stem of D. officinale using a single half-run on the Roche 454 GS FLX Titanium platform to generate 553,084 ESTs with an average length of 417 bases. The ESTs were assembled into 36,407 unique putative transcripts. A total of 69.97% of the unique sequences were annotated, and a detailed view of alkaloid biosynthesis was obtained. Functional assignment based on Kyoto Encyclopedia of Genes and Genomes (KEGG) terms revealed 69 unique sequences representing 25 genes involved in alkaloid backbone biosynthesis. A series of qRT-PCR experiments confirmed that the expression levels of 5 key enzyme-encoding genes involved in alkaloid biosynthesis are greater in the leaves of D. officinale than in the stems. Cytochrome P450s, aminotransferases, methyltransferases, multidrug resistance protein (MDR) transporters and transcription factors were screened for possible involvement in alkaloid biosynthesis. Furthermore, a total of 1061 simple sequence repeat motifs (SSR) were detected from 36,407 unigenes. Dinucleotide repeats were the most abundant repeat type. Of these, 179 genes were associated with a metabolic pathway in KEGG. This study is the first to produce a large volume of transcriptome data from D. officinale. It extends the foundation to facilitate gene discovery in D. officinale and provides an important resource for the molecular genetic and functional genomic studies in this species.

First identification of polymorphic microsatellite markers in the Burgundy truffle, Tuber aestivum (Tuberaceae)

PREMISE OF THE STUDY

Tuber aestivum, the most common truffle in Europe, plays an important role in the commercial truffle market. For the first time, microsatellite primers were developed to investigate polymorphism within this species. •

METHODS AND RESULTS

Using direct shotgun pyrosequencing, 15 polymorphic microsatellites were identified out of the 7784 perfect microsatellites present in the 534620 reads obtained. Tested on 75 samples, these microsatellites were highly polymorphic. The number of alleles varied from four to 15, and the expected heterozygosity ranged from 0.266 to 0.620. A multilocus analysis allowed the identification of 63 genotypes over the 75 samples analyzed. •

CONCLUSIONS

Direct shotgun pyrosequencing is a fast and relatively low-cost technique allowing identification of microsatellites in nonmodel species. The microsatellites developed in this study will be useful in population genetic studies to infer the evolutionary history of this species.

Genome-wide analysis of simple sequence repeats in the model medicinal mushroom Ganoderma lucidum

Simple sequence repeats (SSRs) or microsatellites are one of the most popular sources of genetic markers and play a significant role in gene function and genome organization. We identified SSRs in the genome of Ganoderma lucidum and analyzed their frequency and distribution in different genomic regions. We also compared the SSRs in G. lucidum with six other Agaricomycetes genomes: Coprinopsis cinerea, Laccaria bicolor, Phanerochaete chrysosporium, Postia placenta, Schizophyllum commune and Serpula lacrymans. Based on our search criteria, the total number of SSRs found ranged from 1206 to 6104 and covered from 0.04% to 0.15% of the fungal genomes. The SSR abundance was not correlated with the genome size, and mono- to tri-nucleotide repeats outnumbered other SSR categories in all of the species examined. In G. lucidum, a repertoire of 2674 SSRs was detected, with mono-nucleotides being the most abundant. SSRs were found in all genomic regions and were more abundant in non-coding regions than coding regions. The highest SSR relative abundance was found in introns (108 SSRs/Mb), followed by intergenic regions (84 SSRs/Mb). A total of 684 SSRs were found in the protein-coding sequences (CDSs) of 588 gene models, with 81.4% of them being tri- or hexa-nucleotides. After scanning for InterPro domains, 280 of these genes were successfully annotated, and 215 of them could be assigned to Gene Ontology (GO) terms. SSRs were also identified in 28 bioactive compound synthesis-related gene models, including one 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), three polysaccharide biosynthesis genes and 24 cytochrome P450 monooxygenases (CYPs). Primers were designed for the identified SSR loci, providing the basis for the future development of SSR markers of this medicinal fungus.

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.

The importance of individuals: intraspecific diversity of mycorrhizal plants and fungi in ecosystems

A key component of biodiversity is the number and abundance of individuals (i.e. genotypes), and yet such intraspecific diversity is rarely considered when investigating the effects of biodiversity of mycorrhizal plants and fungi on ecosystem processes. Within a species, individuals vary considerably in important reproductive and functional attributes, including carbon fixation, mycelial growth and nutrient utilization, but this is driven by both genetic and environmental (including climatic) factors. The interactions between individual plants and mycorrhizal fungi can have important consequences for the maintenance of biodiversity and regulation of resource transfers in ecosystems. There is also emerging evidence that assemblages of genotypes may affect ecosystem processes to a similar extent as assemblages of species. The application of whole-genome sequencing and population genomics to mycorrhizal plants and fungi will be crucial to determine the extent to which individual variation in key functional attributes is genetically based. We argue the need to unravel the importance of the diversity (especially assemblages of different evenness and richness) of individuals of both mycorrhizal plants and fungi, and the need to take a 'community genetics' approach to better understand the functional significance of the biodiversity of mycorrhizal symbioses.

Intraspecific genotypic variability determines concentrations of key truffle volatiles

• Aroma variability in truffles has been attributed to maturation (Tuber borchii), linked to environmental factors (Tuber magnatum), but the involvement of genetic factors has been ignored. We investigated aroma variability in Tuber uncinatum, a species with wide distribution. Our aim was to assess aroma variability at different spatial scales (i.e. trees, countries) and to quantify how aroma was affected by genotype, fruiting body maturity, and geographical origin. • A volatile fingerprinting method was used to analyze the aroma of 223 T. uncinatum fruiting bodies from seven European countries. Maturity was estimated from spore melanization. Genotypic fingerprinting was performed by amplified fragment length polymorphism (AFLP). • Discriminant analysis revealed that, regardless of the geographical origin of the truffles, most of the aroma variability was caused by eight-carbon-containing volatiles (C8-VOCs). In an orchard of T. uncinatum, truffles producing different concentrations of C8-VOCs clustered around distinct host trees. This clustering was not associated with maturity, but was associated with fungal genotype. • These results indicate that the variation in C8-VOCs in truffles is most likely under genetic control. They exemplify that understanding the factors behind aroma variability requires a holistic approach. Furthermore, they also raise new questions regarding the ecological role of 1-octen-3-ol in truffles.

Mutualistic interactions on a knife-edge between saprotrophy and pathogenesis

Saprophytic, ectomycorrhizal (ECM) and pathogenic fungi play a key role in carbon and nutrient cycling in forest ecosystems. Whereas more than 50 genomes of saprotrophic and pathogenic fungi have been published, only two genomes of ECM fungi, Laccaria bicolor and Tuber melanosporum, have been released. Comparative analysis of the genomes of biotrophic species highlighted convergent evolution. Mutualistic and pathogenic biotrophic fungi share expansion of genome size through transposon proliferation and common strategies to avoid plant detection. Differences mainly rely on nutritional strategies. Such analyses also pinpointed how blurred the molecular boundaries are between saprotrophism, symbiosis and pathogenesis. Sequencing of additional ECM species, as well as soil saprotrophic fungi, will facilitate the identification of conserved traits for ECM symbiosis and those leading to the transition from white-rotting and brown-rotting to the ECM lifestyle.

Population genetics of ectomycorrhizal fungi: from current knowledge to emerging directions

Ectomycorrhizal (EM) fungi are major microbial components of boreal, temperate and Mediterranean forests, as well as some tropical forest ecosystems. Nearly two decades of studies have clarified many aspects of their population biology, based on several model species from diverse lineages of fungi where the EM symbiosis evolved, i.e. among Hymenomycetes and, to a lesser extent, among Ascomycetes. In this review, we show how tools for individual recognition have changed, shifting from the use of somatic incompatibility reactions to dominant and non-specific markers (such as random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP)) and, more recently, to co-dominant and specific markers (such as microsatellites and single nucleotide polymorphisms (SNPs)). At the same time, the theoretical focus has also changed. In earlier studies, a major aim was the description of genet size and popul/ation strategy. For example, we show how some studies supported or challenged the simple, classical model of colonization of new forest stands by ruderal (R) species, propagating by spores and forming small genets, progressively replaced in older forests by more competitive (C) species, propagating by mycelial growth and forming larger genets. By contrast, more recent studies give insights into some genetic traits, such as partners' assortment (allo- versus autogamy), genetic structure of populations and gene flow that turn out to depend both on distance and on whether spores are animal- or wind-dispersed. We discuss the rising awareness that (i) many morphospecies contain cryptic biological species (often sympatric) and (ii) trans- and inter-continental species may often contain several biological species isolated by distance. Finally, we show the emergence of biogeographic approaches and call for some aspects to be developed, such as fine-scale and long-term population monitoring, analyses of subterranean populations of extra-radical mycelia, or more model species from the tropics, as well as from the Ascomycetes (whose genetic idiosyncrasies are discussed). With the rise of the '-omics' sciences, analysis of population structure for non-neutral genes is expected to develop, and forest management and conservation biology will probably profit from published and expected work.

On the road to understanding truffles in the underground

The genome of the ectomycorrhizal ascomycete Tubermelanosporum has recently been published and this has given researchers unique opportunities to learn more about the biology of this precious edible fungus. The epigeous ascomycete lives in Mediterranean countries in symbiotic interaction with roots of broad-leaf trees such as oaks and hazel. A most important new finding was the single mating type locus in the genome that occurs with two alleles in natural populations. The life cycle is now confirmed to be heterothallic and the species is outcrossing. Unlike sexual development in the soil, mycorrhization of the roots by homokaryotic haploid mycelia is mating-type-independent. Gene regulation during mycorrhization and fruiting and environmental influences on it is now genome-wide addressed. Genome profiling for functions in specific metabolic pathways is undertaken. Insights in most enthralling features of tubers such as on odor formation are thus gained.

Deep RNA sequencing improved the structural annotation of the Tuber melanosporum transcriptome

• The functional complexity of the Tuber melanosporum transcriptome has not yet been fully elucidated. Here, we applied high-throughput Illumina RNA-sequencing (RNA-Seq) to the transcriptome of T. melanosporum at different major developmental stages, that is free-living mycelium, fruiting body and ectomycorrhiza. • Sequencing of cDNA libraries generated a total of c. 24 million sequence reads representing > 882 Mb of sequence data. To construct a coverage signal profile across the genome, all reads were then aligned to the reference genome assembly of T. melanosporum Mel28. • We were able to identify a substantial number of novel transcripts, antisense transcripts, new exons, untranslated regions (UTRs), alternative upstream initiation codons and upstream open reading frames. • This RNA-Seq analysis allowed us to improve the genome annotation. It also provided us with a genome-wide view of the transcriptional and post-transcriptional mechanisms generating an increased number of transcript isoforms during major developmental transitions in T. melanosporum.