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

Unraveling the Polysaccharide Biosynthesis Potential of Ganoderma lucidum: A Chromosome-Level Assembly Using Hi-C Sequencing

Ganoderma lucidum exhibits the ability to synthesize a diverse range of biologically active molecules with significant pharmaceutical potential, including xylomannan and fucogalactan, which have demonstrated antitumor activity. However, there exists considerable intra-species variability in the capacity to produce these metabolites at high concentrations, likely reflecting the high genomic diversity observed from a limited number of strains sequenced to date. We employed high-throughput shotgun sequencing to obtain the complete genome sequence of G. lucidum strain 5.1, which is distinguished by its remarkable xylomannan synthesis capabilities. Through the utilization of semi-automatic reordering based on conformation capture (Hi-C) data, we substantially enhanced the assembly process, resulting in the generation of 12 chromosome-level scaffolds with a cumulative length of 39 Mbp. By employing both de novo and homology-based approaches, we performed comprehensive annotation of the genome, thereby identifying a diverse repertoire of genes likely involved in polysaccharide biosynthesis. The genome sequence generated in this study serves as a valuable resource for elucidating the molecular mechanisms underlying the medicinal potential of Ganoderma species, discovering novel pharmaceutically valuable compounds, and elucidating the ecological mechanisms of the species. Furthermore, the chromosome contact map obtained for the first time for this species extends our understanding of 3D fungal genomics and provides insights into the functional and structural organization within the fungal kingdom.

De Novo Mining and Validating Novel Microsatellite Markers to Assess Genetic Diversity in Maruca vitrata (F.), a Legume Pod Borer

Maruca vitrata (Fabricius) is an invasive insect pest capable of causing enormous economic losses to a broad spectrum of leguminous crops. Microsatellites are valuable molecular markers for population genetic studies; however, an inadequate number of M. vitrata microsatellite loci are available to carry out population association studies. Thus, we utilized this insect's public domain databases for mining expressed sequence tags (EST)-derived microsatellite markers. In total, 234 microsatellite markers were identified from 10053 unigenes. We discovered that trinucleotide repeats were the most predominant microsatellite motifs (61.53%), followed by dinucleotide repeats (23.50%) and tetranucleotide repeats (14.95%). Based on the analysis, twenty-five markers were selected for validation in M. vitrata populations collected from various regions of India. The number of alleles (Na), observed heterozygosity (Ho), and expected heterozygosity (He) ranged from 2 to 5; 0.00 to 0.80; and 0.10 to 0.69, respectively. The polymorphic loci showed polymorphism information content (PIC), ranging from 0.09 to 0.72. Based on the genetic distance matrix, the unrooted neighbor-joining dendrogram differentiated the selected populations into two discrete groups. The SSR markers developed and validated in this study will be helpful in population-level investigations of M. vitrata to understand the gene flow, demography, dispersal patterns, biotype differentiation, and host dynamics.

Development of EST-SSR primers and genetic diversity analysis of the southern blight pathogen Sclerotium rolfsii using transcriptome data

Introduction

Sclerotium rolfsii Sacc. is a globally dispersed pathogenic fungus that causes southern blight disease in many crops and Chinese herbal medicine. The high degree of variation and diversity in the fungi altered population genetic structure. Therefore, the important factors of variation within the pathogen population should be considered during the development of management strategies for the disease.

Methods

In this study, S. rolfsii isolates from 13 hosts in 7 provinces of China were collected and analyzed to identify their morphological features and perform molecular characterization. To develop EST-SSR primers, transcriptome sequencing was performed on isolated CB1, and its SSR loci were comprehensively analyzed. In addition, we analyzed the polymorphisms among different populations based on screened EST-SSR primers.

Results

The results showed that all of these clean reads with total 36,165,475 assembled bases were clustered into 28,158 unigenes, ranged from 201 bp to 16,402 bp on the length, of which the average length was 1,284 bp. Of these, the SSR sequence appeared at an average interval of 15.43 kB, and the frequency of SSR was 0.0648 SSR/kB. Polymorphism of 9 primers was observed among 22 populations, and was verified by the Shannon's index (average = 1.414) and polymorphic information index (> 0.50). The genetic diversity analysis revealed diversity in all host populations and geographical populations. Further, molecular variance analysis (AMOVA) showed that the differences between groups were mainly related to geographical location. Based on cluster analysis, the 7 populations were roughly divided into 3 groups, and the results were highly consistent with those based on the geographical location, ultimately aligning with the results of STRUCTURE analysis.

Discussion

The findings build on current knowledge of the distribution of S. rolfsii in the southwest area of China, adding value to current knowledge base on the population structure and genetic diversity of S. rolfsii, specifically in the context of Chinese herbal medicine cultivation in China. Overall, our findings may provide valuable information for breeding of crops with enhanced resistance toward S. rolfsii.

Transcriptome analysis of Ganoderma lingzhi (Agaricomycetes) response to Trichoderma hengshanicum infection

Green mold caused by Trichoderma spp. has become one of the most serious diseases which threatening the production of Ganoderma lingzhi. To understand the possible resistance mechanism of the G. lingzhi response to T. hengshanicum infection, we examined the G. lingzhi transcript accumulation at 0, 12, and 24 h after T. hengshanicum inoculation. The gene expression analysis was conducted on the interaction between G. lingzhi and T. hengshanicum using RNA-seq and digital gene expression (DGE) profiling methods. Transcriptome sequencing indicated that there were 162 differentially expressed genes (DEGs) at three infection time points, containing 15 up-regulated DEGs and 147 down-regulated DEGs. Resistance-related genes thaumatin-like proteins (TLPs) (PR-5s), phenylalanine ammonia-lyase, and Beta-1,3-glucan binding protein were significantly up-regulated. At the three time points of infection, the heat shock proteins (HSPs) genes of G. lingzhi were down-regulated. The down-regulation of HSPs genes led to the inhibition of HSP function, which may compromise the HSP-mediated defense signaling transduction pathway, leading to G. lingzhi susceptibility. Pathway enrichment analyses showed that the main enriched pathways by G. lingzhi after infection were sphingolipid metabolism, ether lipid metabolism, and valine, leucine and isoleucine degradation pathway. Overall, the results described here improve fundamental knowledge of molecular responses to G. lingzhi defense and contribute to the design of strategies against Trichoderma spp.

Comparison of microsatellite distribution in the genomes of Pteropus vampyrus and Miniopterus natalensis (Chiroptera)

BACKGROUND

Microsatellites are a ubiquitous occurrence in prokaryotic and eukaryotic genomes. Microsatellites have become one of the most popular classes of genetic markers due to their high reproducibility, multi-allelic nature, co-dominant mode of inheritance, abundance and wide genome coverage. We characterised microsatellites in the genomes and genes of two bat species, Pteropus vampyrus and Miniopterus natalensis. This characterisation was used for gene ontology analysis and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment of coding sequences (CDS).

RESULTS

Compared to M. natalensis, the genome size of P. vampyrus is larger and contains more microsatellites, but the total diversity of both species is similar. Mononucleotide and dinucleotide repeats were the most diverse in the genome of the two species. In each bat species, the microsatellite bias was obvious. The microsatellites with the largest number of repeat motifs in P. vampyrus from mononucleotide to hexanucleotide were (A)_n, (AC)_n, (CAA)_n, (AAAC)_n, (AACAA)_n and (AAACAA)_n, with frequencies of 97.94%, 58.75%, 30.53%, 22.82%, 54.68% and 22.87%, respectively, while in M. natalensis were (A)_n, (AC)_n, (TAT)_n, (TTTA)_n, (AACAA)_n and (GAGAGG)_n, with of 92.00%, 34.08%, 40.36%, 21.83%, 25.42% and 12.79%, respectively. In both species, the diversity of microsatellites was highest in intergenic regions, followed by intronic, untranslated and exonic regions and lowest in coding regions. Location analysis indicated that microsatellites were mainly concentrated at both ends of the genes. Microsatellites in the CDS are thus subject to higher selective pressure. In the GO analysis, two unique GO terms were found only in P. vampyrus and M. natalensis, respectively. In KEGG enriched pathway, the biosynthesis of other secondary metabolites and metabolism of other amino acids in metabolism pathways were present only in M. natalensis. The combined biological process, cellular components and molecular function ontology are reflected in the GO analysis and six functional enrichments in KEGG annotation, suggesting advantageous mutations during species evolution.

CONCLUSIONS

Our study gives a comparative characterization of the genomes of microsatellites composition in the two bat species. And also allow further study on the effect of microsatellites on gene function as well as provide an insight into the molecular basis for species adaptation to new and changing environments.

Identification of novel SSR markers for predicting the geographic origin of fungus Schizophyllum commune Fr

The fungus Schizophyllum commune is a cosmopolitan basidiomycete, which is popular as an edible, medical mushroom. It causes wood rot and often used as a model object in research. In this study, we analyzed thirty-two genomes of S. commune strains from the NCBI database and designed forty-seven unique SSR DNA markers. The detailed analysis revealed the enrichment of the S. commune genome for CG, GC, CTC, GAG, and TCG motifs. Principal components analysis confirmed the effectiveness of novel SSR DNA markers that preserve the initial heterogeneity of populations. The construction of a network between strains showed single one at a maximum similarity of 38%, and increasing the similarity to 55% breaks the linkage between large groups while separating two new groups containing strains of the population Ru and test cultures S. commune. The amplicons' presence was identified as a sufficient sign of relation of the culture to a specific population. Testing the novel SSR markers allowed to establish a clear delimitation of all groups by geographic location and to differentiate the H4-8 (GCA_000143185.1) strain from the USA population. This research is the basis for the further analysis of S. commune populations and improvement of approaches to determine its genetic diversity.

Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops

Weevils, classified in the family Curculionidae (true weevils), constitute a group of phytophagous insects of which many species are considered significant pests of crops. Within this family, the red palm weevil (RPW), Rhynchophorus ferrugineus, has an integral role in destroying crops and has invaded all countries of the Middle East and many in North Africa, Southern Europe, Southeast Asia, Oceania, and the Caribbean Islands. Simple sequence repeats (SSRs), also termed microsatellites, have become the DNA marker technology most applied to study population structure, evolution, and genetic diversity. Although these markers have been widely examined in many mammalian and plant species, and draft genome assemblies are available for many species of true weevils, very little is yet known about SSRs in weevil genomes. Here we carried out a comparative analysis examining and comparing the relative abundance, relative density, and GC content of SSRs in previously sequenced draft genomes of nine true weevils, with an emphasis on R. ferrugineus. We also used Illumina paired-end sequencing to generate draft sequence for adult female RPW and characterized it in terms of perfect SSRs with 1–6 bp nucleotide motifs. Among weevil genomes, mono- to trinucleotide SSRs were the most frequent, and mono-, di-, and hexanucleotide SSRs exhibited the highest GC content. In these draft genomes, SSR number and genome size were significantly correlated. This work will aid our understanding of the genome architecture and evolution of Curculionidae weevils and facilitate exploring SSR molecular marker development in these species.

Identification of a Novel Metabolic Target for Bioactive Triterpenoids Biosynthesis in Ganoderma lucidum

Triterpenoids are crucial active ingredients of Ganoderma lucidum (G. lucidum) with various health benefits. Yet the low biosynthesis greatly restricts their industrial applications, novel metabolic engineering strategies are needed to further enhance Ganoderma triterpenoids production. Transcription factors play vital roles in the metabolic regulation of terpenoids, which are still scarce to study in G. lucidum. Herein, a transcription factor GlbHLH5 (GenBank No. MZ436906.1) potential for metabolic regulation of Ganoderma triterpenes was identified for the first time. MeJA could increase Ganoderma triterpenoids biosynthesis, and GlbHLH5 significantly responded to MeJA induction, suggesting GlbHLH5 is a new target for Ganoderma triterpenoids overproduction. The regulatory effect of the newly identified target was further validated by homologous gene overexpression and silence in G. lucidum. It’s demonstrated that overexpression of GlbHLH5 significantly increased triterpenoids accumulation and the key enzyme genes transcription in the biosynthetic pathway, while silencing it displayed the opposite effect, indicating GlbHLH5 could positively regulate the triterpenoids biosynthesis by activating the synergistic expression of key enzyme genes in the biosynthetic pathway. Consequently, GlbHLH5 was identified as a positive regulator and novel metabolic target for Ganoderma triterpenoids biosynthesis, it sheds new lights on the regulatory effect regulation and synthetic biology of Ganoderma triterpenoids.

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.

Distinguishing Homokaryons and Heterokaryons in Medicinal Polypore Mushroom Wolfiporia cocos (Agaricomycetes) Based on Cultural and Genetic Characteristics

The sclerotia of Wolfiporia cocos are a kind of traditional medicine and food with excellent benefits and are widely used in China, Japan, and other Asian countries. The mating system of fungi is not only of practical importance for breeding but also has profound effects on genetic variability and molecular evolution. However, the lack of clamp connections in W. cocos increases the difficulty of research on mating systems. In this study, homokaryons and heterokaryons were distinguished by comparing the characteristics of culture, fruiting tests, and molecular markers, which was further demonstrated by k-mer analysis based on Illumina sequencing. Uninucleate, binucleate, and nuclei-free condition basidiospores of W. cocos were observed, and binucleate basidiospores were the most predominant. Brown-type colonies, slow growth rates in both PDA medium and sawdust substrate, and neutral pH after the growth of mycelia and unfruiting were found to be the morphological and growth characteristics of homokaryotic strains. Primers SSR37 and 38 were screened to identify homokaryons. K-mer analysis based on Illumina sequencing exhibited different heterozygous ratios for homokaryons and heterokaryons. The results revealed that pseudo-homothallism was the predominant mode of reproduction in the Chinese population of W. cocos, and heterothallism also existed in all probability. This study will be helpful for the cross-breeding of this precious medicinal mushroom and for understanding its evolution and population structure.

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.

Polysaccharide biosynthetic pathway profiling and putative gene mining of Dendrobium moniliforme using RNA-Seq in different tissues

BACKGROUND

Dendrobium moniliforme (Linnaeus) Swartz is a well-known plant used in traditional Chinese medicine due to bioactive constituents. Polysaccharides are the main medicinal ingredients, yet no studies have been published on polysaccharide biosynthesis in D. moniliforme. To comprehensively investigate the polysaccharide at the transcription level, we performed de novo transcriptome sequencing for the first time to produce a comprehensive transcriptome of D. moniliforme.

RESULTS

In our study, a database of 562,580 unigenes (average length = 1115.67 bases) was generated by performing transcriptome sequencing. Based on the gene annotation of the transcriptome, we identified 1204 carbohydrate-active related unigenes against CAZy database, including 417 glycosyltransferase genes (GTs), 780 glycoside hydrolases (GHs), 19 carbohydrate esterases (CEs), 75 carbohydrate-binding modules (CBMs), and 44 polysaccharide lyases (PLs). In the cellulose synthase family, 21 differential expression genes (DEGs) related to polysaccharide were identified. Subsequently, the tissue-specific expression patterns of the genes involved in polysaccharide pathway were investigated, which provide understanding of the biosynthesis and regulation of DMP at the molecular level. The two key enzyme genes (Susy and SPS) involved in the polysaccharide pathway were identified, and their expression patterns in different tissues were further analyzed using quantitative real-time PCR.

CONCLUSIONS

We determined the content of polysaccharides from Dendrobium moniliforme under different tissues, and we obtained a large number of differential genes by transcriptome sequencing. This database provides a pool of candidate genes involved in biosynthesis of polysaccharides in D. moniliforme. Furthermore, the comprehensive analysis and characterization of the significant pathways are expected to give a better insight regarding the diversity of chemical composition, synthetic characteristics, and the regulatory mechanism which operate in this medical herb.

Genome-wide mining of microsatellites in king cobra (Ophiophagus hannah) and cross-species development of tetranucleotide SSR markers in Chinese cobra (Naja atra)

The complete genome sequence provides the opportunity for genome-wide and coding region analysis of SSRs in the king cobra and for cross-species identification of microsatellite markers in the Chinese cobra. In the Ophiophagus hannah genome, tetranucleotide repeats (38.03%) were the most abundant category, followed by dinucleotides (23.03%), pentanucleotides (13.07%), mononucleotides (11.78%), trinucleotides (11.49%) and hexanucleotides (2.6%). Twenty predominant motifs in the O. hannah genome were (A)n (C)n, (AC)n, (AG)n, (AT)n, (AGG)n, (AAT)n, (AAG)n, (AAC)n, (ATG)n, (ATAG)n, (AAGG)n, (ATCT)n, (CCTT)n, (ATTT)n, (AAAT)n, (AATAG)n, (ATTCT)n, (ATATGT)n, (AGATAT)n. In total, 4344 SSRs were found in coding sequences (CDSs). Tetranucleotides (52.79%) were the most abundant microsatellite type in CDS, followed by trinucleotides (28.50%), dinucleotides (11.02%), pentanucleotides (4.42%), mononucleotides (1.77%), and hexanucleotides (1.50%). A total of 984 CDSs containing microsatellites were assigned 11152 Gene Ontology (GO) functional terms. Gene Ontology (GO) analysis demonstrated that cellular process, cell and binding were the most frequent GO terms in biological process, cellular component and molecular function, respectively. Thirty-two novel highly polymorphic (PIC > 0.5) SSR markers for Naja atra were developed from cross-species amplification based on the tetranucleotide microsatellite sequences in the king cobra genome. The number of alleles (N_A) per locus had between 3 and 11 alleles with an average of 6.5, the polymorphism information content (PIC) value ranged from 0.521 to 0.858 (average = 0.707), the observed heterozygosity (Ho) of 32 microsatellite loci ranged from 0.292 to 0.875 (mean = 0.678), the expected heterozygosity (H_E) ranged from 0.561 to 0.889 (average = 0.761), and 3 microsatellite loci exhibited statistically significant departure from Hardy-Weinberg equilibrium (HWE) after Bonferroni correction (p < 0.003).

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.

Metabolism of ganoderic acids by a Ganoderma lucidum cytochrome P450 and the 3-keto sterol reductase ERG27 from yeast

Ganoderic acids, a group of oxygenated lanostane-type triterpenoids, are the major bioactive compounds produced by the well-known medicinal macro fungus Ganoderma lucidum. More than 150 ganoderic acids have been identified, and the genome of G. lucidum has been sequenced recently. However, the biosynthetic pathways of ganoderic acids have not yet been elucidated. Here, we report the functional characterization of a cytochrome P450 gene CYP512U6 from G. lucidum, which is involved in the ganoderic acid biosynthesis. CYP512U6 hydroxylates the ganoderic acids DM and TR at the C-23 position to produce hainanic acid A and ganoderic acid Jc, respectively. In addition, CYP512U6 can also hydroxylate a modified ganoderic acid DM in which the C-3 ketone has been reduced to hydroxyl by the sterol reductase ERG27 from Saccharomyces cerevisiae. An NADPH-dependent cytochrome P450 reductase from G. lucidum was also isolated and characterized. These results will help elucidate the biosynthetic pathways of ganoderic acids.

Genome microsatellite diversity within the Apicomplexa phylum

The Apicomplexa phylum groups include unicellular and obligate intracellular protozoan parasites with an apical complex used for attachment and invasion to host cells. In this study, we analyze single sequence repeats (SSRs) in the whole genome of 20 apicomplexan organisms that represent four different lineages within the phylum. Only perfect SSRs with at least 12 nucleotides and composed of 2-6 mers were included. To better understand the association of SSR types with the genomic regions, the SSRs were classified accordingly with the genomic location into exon, intron and intergenic categories. Our results showed heterogeneous SSRs density within the studied genomes. However, the most frequent SSRs types were di- and tri-nucleotide repeats. The former was associated with intergenic regions, while the latter was associated with exon 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.

A comprehensive characterization of simple sequence repeats in pepper genomes provides valuable resources for marker development in Capsicum

The sequences of the full set of pepper genomes including nuclear, mitochondrial and chloroplast are now available for use. However, the overall of simple sequence repeats (SSR) distribution in these genomes and their practical implications for molecular marker development in Capsicum have not yet been described. Here, an average of 868,047.50, 45.50 and 30.00 SSR loci were identified in the nuclear, mitochondrial and chloroplast genomes of pepper, respectively. Subsequently, systematic comparisons of various species, genome types, motif lengths, repeat numbers and classified types were executed and discussed. In addition, a local database composed of 113,500 in silico unique SSR primer pairs was built using a homemade bioinformatics workflow. As a pilot study, 65 polymorphic markers were validated among a wide collection of 21 Capsicum genotypes with allele number and polymorphic information content value per marker raging from 2 to 6 and 0.05 to 0.64, respectively. Finally, a comparison of the clustering results with those of a previous study indicated the usability of the newly developed SSR markers. In summary, this first report on the comprehensive characterization of SSR motifs in pepper genomes and the very large set of SSR primer pairs will benefit various genetic studies in Capsicum.

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.

Genome-Wide Survey and Analysis of Microsatellite Sequences in Bovid Species

Microsatellites or simple sequence repeats (SSRs) have become the most popular source of genetic markers, which are ubiquitously distributed in many eukaryotic and prokaryotic genomes. This is the first study examining and comparing SSRs in completely sequenced genomes of the Bovidae. We analyzed and compared the number of SSRs, relative abundance, relative density, guanine-cytosine (GC) content and proportion of SSRs in six taxonomically different bovid species: Bos taurus, Bubalus bubalis, Bos mutus, Ovis aries, Capra hircus, and Pantholops hodgsonii. Our analysis revealed that, based on our search criteria, the total number of perfect SSRs found ranged from 663,079 to 806,907 and covered from 0.44% to 0.48% of the bovid genomes. Relative abundance and density of SSRs in these Bovinae genomes were non-significantly correlated with genome size (Pearson, r < 0.420, p > 0.05). Perfect mononucleotide SSRs were the most abundant, followed by the pattern: perfect di- > tri- > penta- > tetra- > hexanucleotide SSRs. Generally, the number of SSRs, relative abundance, and relative density of SSRs decreased as the motif repeat length increased in each species of Bovidae. The most GC-content was in trinucleotide SSRs and the least was in the mononucleotide SSRs in the six bovid genomes. The GC-contents of tri- and pentanucleotide SSRs showed a great deal of similarity among different chromosomes of B. taurus, O. aries, and C. hircus. SSR number of all chromosomes in the B. taurus, O.aries, and C. hircus is closely positively correlated with chromosome sequence size (Pearson, r > 0.980, p < 0.01) and significantly negatively correlated with GC-content (Pearson, r < -0.638, p < 0.01). Relative abundance and density of SSRs in all chromosomes of the three species were significantly negatively correlated with GC-content (Pearson, r < -0.333, P < 0.05) but not significantly correlated with chromosome sequence size (Pearson, r < -0.185, P > 0.05). Relative abundances of the same nucleotide SSR type showed great similarity among different chromosomes of B. taurus, O. aries, and C. hircus.

Genome-wide survey and analysis of microsatellites in giant panda (Ailuropoda melanoleuca), with a focus on the applications of a novel microsatellite marker system

Background

The giant panda (Ailuropoda melanoleuca) is a critically endangered species endemic to China. Microsatellites have been preferred as the most popular molecular markers and proven effective in estimating population size, paternity test, genetic diversity for the critically endangered species. The availability of the giant panda complete genome sequences provided the opportunity to carry out genome-wide scans for all types of microsatellites markers, which now opens the way for the analysis and development of microsatellites in giant panda.

Results

By screening the whole genome sequence of giant panda in silico mining, we identified microsatellites in the genome of giant panda and analyzed their frequency and distribution in different genomic regions. Based on our search criteria, a repertoire of 855,058 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. A total of 160 primer pairs were designed to screen for polymorphic microsatellites using the selected tetranucleotide microsatellite sequences. The 51 novel polymorphic tetranucleotide microsatellite loci were discovered based on genotyping blood DNA from 22 captive giant pandas in this study. Finally, a total of 15 markers, which showed good polymorphism, stability, and repetition in faecal samples, were used to establish the novel microsatellite marker system for giant panda. Meanwhile, a genotyping database for Chengdu captive giant pandas (n = 57) were set up using this standardized system. What’s more, a universal individual identification method was established and the genetic diversity were analysed in this study as the applications of this marker system.

Conclusion

The microsatellite abundance and diversity were characterized in giant panda genomes. A total of 154,677 tetranucleotide microsatellites were identified and 15 of them were discovered as the polymorphic and stable loci. The individual identification method and the genetic diversity analysis method in this study provided adequate material for the future study of giant panda.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1268-z) contains supplementary material, which is available to authorized users.

De novo assembly of Auricularia polytricha transcriptome using Illumina sequencing for gene discovery and SSR marker identification

Auricularia polytricha (Mont.) Sacc., a type of edible black-brown mushroom with a gelatinous and modality-specific fruiting body, is in high demand in Asia due to its nutritional and medicinal properties. Illumina Solexa sequenceing technology was used to generate very large transcript sequences from the mycelium and the mature fruiting body of A. polytricha for gene discovery and molecular marker development. De novo assembly generated 36,483 ESTs with an N50 length of 636 bp. A total of 28,108 ESTs demonstrated significant hits with known proteins in the nr database, and 94.03% of the annotated ESTs showed the greatest similarity to A. delicata, a related species of A. polytricha. Functional categorization of the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways revealed the conservation of genes involved in various biological processes in A. polytricha. Gene expression profile analysis indicated that a total of 2,057 ESTs were differentially expressed, including 1,020 ESTs that were up-regulated in the mycelium and 1,037 up-regulated in the fruiting body. Functional enrichment showed that the ESTs associated with biosynthesis, metabolism and assembly of proteins were more active in fruiting body development. The expression patterns of homologous transcription factors indicated that the molecular mechanisms of fruiting body formation and development were not exactly the same as for other agarics. Interestingly, an EST encoding tyrosinase was significantly up-regulated in the fruiting body, indicating that melanins accumulated during the processes of the formation of the black-brown color of the fruiting body in A. polytricha development. In addition, a total of 1,715 potential SSRs were detected in this transcriptome. The transcriptome analysis of A. polytricha provides valuable sequence resources and numerous molecular markers to facilitate further functional genomics studies and genetic researches on this fungus.

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.