Deep Sequencing, Nested PCR, and Denaturing Gradient Gel Electrophoresis Reveal a Wider Distribution of Chinese Caterpillar Mushroom, Ophiocordyceps sinensis (Ascomycetes), in Native Soil Types.

Ophiocordyceps sinensis appears as stroma emerging from underground sclerotium enclosed by the skeleton of Thitarodes moth larvae. However, the actual distribution of the fungus in soil still remains unclarified. In this study, 40 soil samples were used for detection of O. sinensis to confirm its distribution in native habitats using denaturing gradient gel electrophoresis, nested internal transcribed spacer (ITS) PCR, and 454 pyrosequencing methods. The soil samples included six types: Os, where both stromata and host moth larvae were found; NL, representing no signs of stromata, but where moth larvae were found; NOs, where neither stroma nor moth larvae were found; BS, with bare soil without the presence of stroma of O. sinensis or moth larvae; AF, from soil surrounding the stroma; and MP, soil particles firmly wrapping the sclerotium of O. sinensis. Of 40 samples tested, 36 showed positive detection of O. sinensis by at least one of the three detection methods, with positive detection in all six sample types at all five sites. The results showed that traces of O. sinensis can be detected in locations with no macroscopically visible evidence of the fungus or its host and at least 100 m away from such locations.

The Reliability of DNA Sequences in Public Databases Belonging to the Most Economically Important Shiitake Culinary-Medicinal Mushroom Lentinus edodes (Agaricomycetes) in Asia.

Large numbers of DNA sequences deposited in the International Nucleotide Sequence Databases (INSD) are erroneously annotated. The erroneous information may lead to misleading conclusions or cause great economic losses to farmers. Lentinus edodes (= Lentinula edodes (Berk.) Pegler) is one of the most important and popular culinary-medicinal mushrooms with a high nutritional value. In this study, experimental and in silico methods were used to correct the sequences annotated as L. edodes in the INSD. A total of 3,426 nucleotide entries were retrieved from public databases, including 140 different types of genetic sequences. Excluding 1,893 genome sequences, the most abundant signatures represented by ITS (258) and IGS1 (259) sequences accounted for 33.23% of the total entries. A total of 3,058 sequences were annotated correctly, 350 were indeterminate, and 18 were annotated erroneously based on the two methods. Correction of sequences will be beneficial for species identification and annotation. Phylogenic analysis based on ITS sequences suggested that L. edodes segregate in four clades in the tree based on ITS sequences. The isolates from China were distributed into two clades. In L. edodes, the intraspecific variation of the ITS2 sequences was much higher than that of the ITS1 sequences. In addition, the genetic diversity of the L. edodes sequences from China was much higher than that of any other regions included in this study. The northwest and southwest regions of China were L. edodes diversity centers.

Comparison of different sequencing and assembly strategies for a repeat-rich fungal genome, Ophiocordyceps sinensis.

Ophiocordyceps sinensis is one of the most expensive medicinal fungi world-wide, and has been used as a traditional Chinese medicine for centuries. In a recent report, the genome of this fungus was found to be expanded by extensive repetitive elements after assembly of Roche 454 (223Mb) and Illumina HiSeq (10.6Gb) sequencing data, producing a genome of 87.7Mb with an N50 scaffold length of 12kb and 6972 predicted genes. To test whether the assembly could be improved by deeper sequencing and to assess the amount of data needed for optimal assembly, genomic sequencing was run several times on genomic DNA extractions of a single ascospore isolate (strain 1229) on an Illumina HiSeq platform (25Gb total data). Assemblies were produced using different data types (raw vs. trimmed) and data amounts, and using three freely available assembly programs (ABySS, SOAP and Velvet). In nearly all cases, trimming the data for low quality base calls did not provide assemblies with higher N50 values compared to the non-trimmed data, and increasing the amount of input data (i.e. sequence reads) did not always lead to higher N50 values. Depending on the assembly program and data type, the maximal N50 was reached with between 50% to 90% of the total read data, equivalent to 100× to 200× coverage. The draft genome assembly was improved over the previously published version resulting in a 114Mb assembly, scaffold N50 of 70kb and 9610 predicted genes. Among the predicted genes, 9213 were validated by RNA-Seq analysis in this study, of which 8896 were found to be singletons. Evidence from genome and transcriptome analyses indicated that species assemblies could be improved with defined input material (e.g. haploid mono-ascospore isolate) without the requirement of multiple sequencing technologies, multiple library sizes or data trimming for low quality base calls, and with genome coverages between 100× and 200×.

Bacterial diversity in native habitats of the medicinal fungus Ophiocordyceps sinensis on Tibetan Plateau as determined using Illumina sequencing data.

Ophiocordyceps sinensis is one of the most well-known traditional Chinese medicinal fungi. In this study, bacterial diversity in the soils of native habitats of O. sinensis was investigated using Illumina sequencing data. A total of 525,000 sequences of V6-16S rRNA were analyzed. The number of OTUs from each sample ranged from 13,858 to 15,978 at 97% sequence similarity cut-off. The results demonstrated that the deep sequencing approach provides improved access to rare genotypes. Richness indices and Shannon's diversity index did not differ significantly between samples collected from locations where O. sinensis was present (Os1-3) and not present (NOs1-3). Classified bacterial sequences were grouped into 23 phyla including Proteobacteria, Actinobacteria, Acidobacteria, Verrucomicrobia, etc. The Venn diagram revealed that 7183 OTUs belonging to 14 phyla were shared by Os, NOs and MP (mycelial pellicle wrapping the sclerotium of O. sinensis) samples, possibly representing a core microbiome existing in native habitats of O. sinensis, and that 863 belonging to 12 phyla were shared by Os and MP samples, possibly related to the occurrence of O. sinensis. Overall, the results revealed a high bacterial diversity in the soil samples and the relationships between the bacterial diversity and O. sinensis merit further investigation.