Genomic Comparative Analysis of Cordyceps pseudotenuipes with Other Species from Cordyceps

Abstract

The whole genome of Cordyceps pseudotenuipes was sequenced, annotated, and compared with three related species to characterize the genome. The antibiotics and Secondary Metabolites Analysis Shell (antiSMASH) and local BLAST analysis were used to explore the secondary metabolites (SMs) and biosynthesis gene clusters (BGCs) of the genus Cordyceps. The genome-wide basic characteristics of C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris revealed unequal genome size, with C. cicadae as the largest (34.11 Mb), followed by C. militaris (32.27 Mb). However, the total gene lengths of C. pseudotenuipes and C. tenuipes were similar (30.1 Mb and 30.06 Mb). The GC contents of C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris genomes differed slightly (51.40% to 54.11%). AntiSMASH and local BLAST analysis showed that C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris had 31, 28, 31, and 29 putative SM BGCs, respectively. The SM BGCs contained different quantities of polyketide synthetase (PKS), nonribosomal peptide synthetase (NRPS), terpene, hybrid PKS + NRPS, and hybrid NRPS + Other. Moreover, C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris had BGCs for the synthesis of dimethylcoprogen. C. pseudotenuipes, C. tenuipes, and C. cicadae had BGCs for the synthesis of leucinostatin A/B, neosartorin, dimethylcoprogen, wortmanamide A/B, and beauvericin. In addition, the SM BGCs unique to C. pseudotenuipes were clavaric acid, communesin, and deoxynivalenol. Synteny analysis indicated that the scaffolds where the SM BGC was located were divided into more than 70 collinear blocks, and there might be rearrangements. Altogether, these findings improved our understanding of the molecular biology of the genus Cordyceps and will facilitate the discovery of new biologically active SMs from the genus Cordyceps using heterologous expression and gene knockdown methods.