Complete mitochondrial genome of the medicinal fungus Ophiocordyceps sinensis

Characterization and Phylogenetic Analysis of the Complete Mitogenomes of Valsa mali and Valsa pyri

Apple Valsa canker, caused by Valsa mali and Valsa pyri, is a devastating disease of apple trees and poses a severe threat to the sustainable development of apple production. Although the two species' whole genomes have been sequenced, their mitochondrial genomes are still uncharacterized. In this study, the complete mitochondrial genomes of V. mali and V. pyri were assembled, annotated, and compared by bioinformatic methods. The results indicate that the mitogenomes are both circular DNA molecules with sizes of 213,406 bp and 128,022 bp, respectively. The AT skew values of the two Valsa species' mitogenomes were positive, while the GC skew values were negative. Comparative mitogenome analysis revealed that the length and base composition of protein-coding genes (PCGs), rRNA genes, and tRNA genes differed between the two Valsa species. It was found that the expansion of V. mali was primarily attributable to the intronic regions. There are large numbers of interspersed repetitive sequences (IRS) in both Valsa mitogenomes; however, the proportion of IRS in V. mali (43.56%) was much higher than that in V. pyri (2.41%). The alignment of large fragments between the mitochondrial and nuclear genomes of both V. mali (1.73 kb) and V. pyri (5.17 kb) indicates that gene transfer between mitochondrial and nuclear genomes occurred during evolution. The ka/ks ratios for 15 core PCGs were below one, suggesting that these genes were subjected to purifying selection pressure. Comparative mitogenomics revealed that the two fungi had significant mitogenomic collinearity and large-scale gene rearrangements. The results of phylogenetic analysis based on Bayesian inference (BI) and maximum likelihood (ML) using a combined mitochondrial gene set confirmed that V. mali and V. pyri were fully independent taxa with a high bootstrap value of 100 (ML) and a high posterior probability of 1.0 (BI). This is the first report on the mitogenomes within the genus Valsa. These results will pave the way to understanding the evolution and differentiation of mitogenomes in the genus Valsa.

Complete mitochondrial genome of Paraisaria gracilioides (Hypocreales: Ophiocordycipitaceae) and phylogenetic analysis

Paraisaria gracilioides (Kobayasi) Luangsa-ard, Mongkolsamrit & Samson, originally described as Ophiocordyceps gracilioides, is a member of the Ophiocordycipitaceae, within the order Hypocreales. This study, presented the mitochondrial genome of P. gracilioides structured as a circular molecule spanning 181,140 bp, was larger than those of most Ophiocordycipitaceae species. Despite this size variation, encoding 49 genes, including 15 PCGs, 26 tRNAs, two rRNAs, and six ORFs. Phylogenetic analysis based on nucleotide sequences of 14 PCGs revealed that P. gracilioides clustered with Paraisaria gracilis, forming a subgroup within Ophiocordyceps species. The results contribute to refining the taxonomic framework of Ophiocordycipitaceae.

Mitogenomics, phylogeny and morphology reveal two new entomopathogenic species of Ophiocordyceps (Ophiocordycipitaceae, Hypocreales) from south-western China

Ophiocordyceps encompasses over 300 species, demonstrating a wide range of morphological features, hosts and habitats within its species diversity. In this study, two novel species in Ophiocordyceps were revealed parasitising Hepialidae larva buried in soil. Ophiocordycepsalbastroma was morphologically characterised by white stromata, solitary and cylindrical conidiogenous cells and smooth ovoid or ellipsoidal conidia. Ophiocordycepsnigristroma was characterised by woody and dark brown stromata, monophialidic, swollen base and lageniform conidiogenous cells and smooth fusiform or oval conidia. The two new species formed a separate clade, respectively, based on the phylogenetic analyses of a combined dataset including nrSSU, nrLSU, rpb1, rpb2, and tef-1α, as well as a dataset of mitochondrial 14 protein coding genes (PCGs). They were all closely grouped with O.sinensis. The mitochondrial genomes of them were first reported. Their mitogenomes were all typical of circular molecules, with positive AT and GC skew, similar GC content, similar genetic composition, similar codon usage and conservative gene positions. However, the length of the mitogenomes varied. Changes in the length of the genes were the leading cause of changes in the length of mitochondrial genome of Ophiocordyceps. The discovery and identification of new Ophiocordyceps species and analysis their mitochondrial genomes may serve as foundations for phylogeny and diversity research within the genus Ophiocordyceps.

Comparative Analysis of the Mitochondrial Genome Sequences of Diaporthe longicolla (syn. Phomopsis longicolla) Isolates Causing Phomopsis Seed Decay in Soybean

Diaporthe longicolla (syn. Phomopsis longicolla) is an important seed-borne fungal pathogen and the primary cause of Phomopsis seed decay (PSD) in soybean. PSD is one of the most devastating seed diseases, reducing soybean seed quality and yield worldwide. As part of a genome sequencing project on the fungal Diaporthe-Phomopsis complex, draft genomes of eight D. longicolla isolates were sequenced and assembled. Sequences of mitochondrial genomes were extracted and analyzed. The circular mitochondrial genomes ranged from 52,534 bp to 58,280 bp long, with a mean GC content of 34%. A total of 14 core protein-coding genes, 23 tRNA, and 2 rRNA genes were identified. Introns were detected in the genes of atp6, cob, cox1, cox2, cox3, nad1, nad2, nad5, and rnl. Three isolates (PL7, PL10, and PL185E) had more introns than other isolates. Approximately 6.4% of the mitochondrial genomes consist of repetitive elements. Moreover, 48 single-nucleotide polymorphisms (SNPs) and were identified. The mitochondrial genome sequences of D. longicolla will be useful to further study the molecular basis of seed-borne pathogens causing seed diseases, investigate genetic variation among isolates, and develop improved control strategies for Phomopsis seed decay of soybean.

Tracing the geographical origin of endangered fungus Ophiocordyceps sinensis, especially from Nagqu, using UPLC-Q-TOF-MS

Ophiocordyceps sinensis (OS), known as "soft gold", played an important role in local economic development. OS from different producing areas was difficult to be discriminated by the appearance. Nagqu OS, a distinguished and safeguarded geographical indication product, commands a premium price in market. The real claim of OS geographical origins is urgently required. Here, 81 OS samples were collected from Tibetan Plateau in China to explore markers for tracing origins. OS from Xigazê can be distinguished by dark color of head of caterpillar. Then 57 samples, a fully representative training-sample set, were used to set up OPLS-DA models by nontargeted metabolomics from UPLC-QTOF-MS. Certain markers were successfully identified and validation using 21 blind test samples confirmed that the markers can trace the geographical origin of OS, especially Nagqu samples. It was affirmed that UPLC-QTOF-MS-based untargeted metabolomics coupled with OPLS-DA was a reliable strategy to trace the geographical origins of OS.

Comparative Mitogenomics Analysis Revealed Evolutionary Divergence among Neopestalotiopsis Species Complex (Fungi: Xylariales)

The genus Neopestalotiopsis consists of obligate parasites that cause ring spot, scab, and leaf blight diseases in higher plant species. We assembled the three complete mitogenomes for the guava fruit ring spot pathogen, Neopestalotiopsis cubana. The mitogenomes are circular, with sizes of 38,666 bp, 33,846 bp, and 32,593 bp. The comparative analyses with Pestalotiopsis fici showed that N. cubana differs greatly from it in the length of the mitogenomes and the number of introns. Moreover, they showed significant differences in the gene content and tRNAs. The two genera showed little difference in gene skewness and codon preference for core protein-coding genes (PCGs). We compared gene sequencing in the mitogenomes of the order Xylariales and found large-scale gene rearrangement events, such as gene translocations and the duplication of tRNAs. N. cubana shows a unique evolutionary position in the phylum Ascomycota constructed in phylogenetic analyses. We also found a more concentrated distribution of evolutionary pressures on the PCGs of Neopestalotiopsis in the phylum Ascomycota and that they are under little selective pressure compared to other species and are subjected to purifying selection. This study explores the evolutionary dynamics of the mitogenomes of Neopestalotiopsis and provides important support for genetic and taxonomic studies.

Systematic analyses with genomic and metabolomic insights reveal a new species, Ophiocordyceps indica sp. nov. from treeline area of Indian Western Himalayan region

Ophiocordyceps is a species-rich genus in the order Hypocreales (Sordariomycetes, Ascomycota) depicting a fascinating relationship between microbes and insects. In the present study, a new species, Ophiocordyceps indica sp. nov., is discovered infecting lepidopteran larvae from tree line locations (2,202-2,653 m AMSL) of the Kullu District, Himachal Pradesh, Indian Western Himalayan region, using combinations of morphological and molecular phylogenetic analyses. A phylogeny for Ophiocordyceps based on a combined multigene (nrSSU, nrLSU, tef-1α, and RPB1) dataset is provided, and its taxonomic status within Ophiocordycipitaceae is briefly discussed. Its genome size (~59 Mb) revealed 94% genetic similarity with O. sinensis; however, it differs from other extant Ophiocordyceps species based on morphological characteristics, molecular phylogenetic relationships, and genetic distance. O. indica is identified as the second homothallic species in the family Ophiocordycipitaceae, after O. sinensis. The presence of targeted marker components, viz. nucleosides (2,303.25 μg/g), amino acids (6.15%), mannitol (10.13%), and biological activity data, suggests it to be a new potential source of nutraceutical importance. Data generated around this economically important species will expand our understanding regarding the diversity of Ophiocordyceps-like taxa from new locations, thus providing new research avenues.

Mitochondrial genome of Cordyceps blackwelliae: organization, transcription, and evolutionary insights into Cordyceps

Cordyceps is a diverse genus of insect pathogenic fungi, with about 180 accepted species, including some well-known ones used as ethnic medicine and/or functional food. Nevertheless, mitogenomes are only available for four members of the genus. The current study reports the mitogenome of Cordyceps blackwelliae, a newly described entomopathogenic fungus. The 42,257-bp mitogenome of the fungus encoded genes typically found in fungal mitogenomes, and a total of 14 introns inserted into seven genes, including cob (1 intron), cox1 (4), cox3 (3), nad1 (1), nad4 (1), nad5 (1), and rnl (3). RNA-Seq analysis revealed differential expression of mitochondrial genes and supported annotations resulting from in silico analysis. There was clear evidence for polycistronic transcription and alternative splicing of mitochondrial genes. Comparison among mitogenomes of five different Cordyceps species (i.e., C. blackwelliae, C. chanhua, C. militaris, C. pruinosa, and C. tenuipes) revealed a high synteny, with mitogenome size expansion correlating with intron insertions. Different mitochondrial protein-coding genes showed variable degrees of genetic differentiation among these species, but they were all under purifying selection. Mitochondrial phylogeny based on either nucleotide or amino acid sequences confirmed the taxonomic position of C. blackwelliae in Cordycipitaceae, clustering together with C. chanhua. This study promotes our understanding of fungal evolution in Cordyceps.

Comparative analysis of the complete mitochondrial genomes of four cordyceps fungi

Cordyceps is a large group of entomogenous, medicinally important fungi. In this study, we sequenced, assembled, and annotated the entire mitochondrial genome of Ophiocordyceps xuefengensis, in addition to comparing it against other three complete cordyceps mitogenomes that were previously published. Comparative analysis indicated that the four complete mitogenomes are all composed of circular DNA molecules, although their sizes significantly differ due to high variability in intron and intergenic region sizes in the Ophiocordyceps sinensis and O. xuefengensis mitogenomes. All mitogenomes contain 14 conserved genes and two ribosomal RNA genes, but varying numbers of tRNA introns. The Ka/Ks ratios for all 14 PCGs and rps3 were all less than 1, indicating that these genes have been subject to purifying selection. Phylogenetic analysis was conducted using concatenated amino acid and nucleotide sequences of the 14 PCGs and rps3 using two different methods (Maximum Likelihood and Bayesian analysis), revealing highly supported relationships between O. xuefengensis and other Ophiocordyceps species, in addition to a close relationship with O. sinensis. Further, the analyses indicated that cox1 and rps3 play important roles in population differentiation. These mitogenomes will allow further study of the population genetics, taxonomy, and evolutionary biology of medicinally important cordyceps species.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the pathogenic fungus Ilyonectria destructans

Ilyonectria destructans is a pathogenic fungus causing root rot and other symptoms on trees and many crops. This paper analyses the mitochondrial genome of I. destructans and compares it with other published Nectriaceae mitogenomes. The I. destructans mitogenome appears as a circular DNA molecule of 42,895 bp and an overall GC content of 28.23%. It contains 28 protein-coding genes (15 core protein genes and 13 free-standing ORFs), two rRNAs and 27 tRNAs. The gene content and order were found to be conserved in the mitogenome of I. destructans and other Nectriaceae, although the genome size varies because of the variation in the number and length of intergenic regions and introns. For most core protein-coding genes in Nectriaceae species, Ka/Ks < 1 indicates purifying selection. Among some Nectriaceae representatives, only the rps3 gene was found under positive selection. Phylogenetic analyses based on nucleotide sequences of 15 protein-coding genes divided 45 Hypocreales species into six major clades matching the families Bionectriaceae, Cordycipitaceae, Clavicipitaceae, Ophiocordycipitaceae, Hypocreaceae and Nectriaceae. I. destructans appeared as a sister species to unidentified Ilyonectia sp., closely related to C. ilicicola, N. cinnabarina and a clad of ten Fusarium species and G. moniliformis. The complete mitogenome of I. destructans reported in the current paper will facilitate the study of epidemiology, biology, genetic diversity of the species and the evolution of family Nectriace and the Hypocreales order.

The complete mitochondrial genome of Ophiocordyceps gracilis and its comparison with related species

In this study, the complete mitochondrial genome of O. gracilis was sequenced and assembled before being compared with related species. As the second largest mitogenome reported in the family Ophiocordycipitaceae, the mitogenome of O. gracilis (voucher OG201301) is a circular DNA molecule of 134,288 bp that contains numerous introns and longer intergenomic regions. UCA was detected as anticodon in tRNA-Sec of O. gracilis, while comparative mitogenome analysis of nine Ophiocordycipitaceae fungi indicated that the order and contents of PCGs and rRNA genes were considerably conserved and could descend from a common ancestor in Ophiocordycipitaceae. In addition, the expansion of mitochondrial organization, introns, gene length, and order of O. gracilis were determined to be similar to those of O. sinensis, which indicated common mechanisms underlying adaptive evolution in O. gracilis and O. sinensis. Based on the mitochondrial gene dataset (15 PCGs and 2 RNA genes), a close genetic relationship between O. gracilis and O. sinensis was revealed through phylogenetic analysis. This study is the first to investigate the molecular evolution, phylogenetic pattern, and genetic structure characteristics of mitogenome in O. gracilis. Based on the obtained results, the mitogenome of O. gracilis can increase understanding of the genetic diversity and evolution of cordycipitoid fungi.

Comparative Mitogenomic Analysis and the Evolution of Rhizoctonia solani Anastomosis Groups

Mitochondria are the major energy source for cell functions. However, for the plant fungal pathogens, mitogenome variations and their roles during the host infection processes remain largely unknown. Rhizoctonia solani, an important soil-borne pathogen, forms different anastomosis groups (AGs) and adapts to a broad range of hosts in nature. Here, we reported three complete mitogenomes of AG1-IA RSIA1, AG1-IB RSIB1, and AG1-IC, and performed a comparative analysis with nine published Rhizoctonia mitogenomes (AG1-IA XN, AG1-IB 7/3/14, AG3, AG4, and five Rhizoctonia sp. mitogenomes). These mitogenomes encoded 15 typical proteins (cox1-3, cob, atp6, atp8-9, nad1-6, nad4L, and rps3) and several LAGLIDADG/GIY-YIG endonucleases with sizes ranging from 109,017 bp (Rhizoctonia sp. SM) to 235,849 bp (AG3). We found that their large sizes were mainly contributed by repeat sequences and genes encoding endonucleases. We identified the complete sequence of the rps3 gene in 10 Rhizoctonia mitogenomes, which contained 14 positively selected sites. Moreover, we inferred a robust maximum-likelihood phylogeny of 32 Basidiomycota mitogenomes, representing that seven R. solani and other five Rhizoctonia sp. lineages formed two parallel branches in Agaricomycotina. The comparative analysis showed that mitogenomes of Basidiomycota pathogens had high GC content and mitogenomes of R. solani had high repeat content. Compared to other strains, the AG1-IC strain had low substitution rates, which may affect its mitochondrial phylogenetic placement in the R. solani clade. Additionally, with the published RNA-seq data, we investigated gene expression patterns from different AGs during host infection stages. The expressed genes from AG1-IA (host: rice) and AG3 (host: potato) mainly formed four groups by k-mean partitioning analysis. However, conserved genes represented varied expression patterns, and only the patterns of rps3-nad2 and nad1-m3g18/mag28 (an LAGLIDADG endonuclease) were conserved in AG1-IA and AG3 as shown by the correlation coefficient analysis, suggesting regulation of gene repertoires adapting to infect varied hosts. The results of variations in mitogenome characteristics and the gene substitution rates and expression patterns may provide insights into the evolution of R. solani mitogenomes.

A Comparative Analyses of the Complete Mitochondrial Genomes of Fungal Endosymbionts in Sogatella furcifera, White-Backed Planthoppers

Sogatella furcifera Horvath, commonly known as the white-backed planthoppers (WBPH), is an important pest in East Asian rice fields. Fungal endosymbiosis is widespread among planthoppers in the infraorder Fulgoromorpha and suborder Auchenorrhyncha. We successfully obtained complete mitogenome of five WBPH fungal endosymbionts, belonging to the Ophiocordycipitaceae family, from next-generation sequencing (NGS) reads obtained from S. furcifera samples. These five mitogenomes range in length from 55,390 bp to 55,406 bp, which is shorter than the mitogenome of the fungal endosymbiont found in Ricania speculum, black planthoppers. Twenty-eight protein-coding genes (PCGs), 12 tRNAs, and 2 rRNAs were found in the mitogenomes. Two single-nucleotide polymorphisms, two insertions, and three deletions were identified among the five mitogenomes, which were fewer in number than those of four species of Ophiocordycipitaceae, Ophiocordyceps sinensis, Hirsutella thompsonii, Hirsutella rhossiliensis, and Tolypocladium inflatum. Noticeably short lengths (up to 18 bp) of simple sequence repeats were identified in the five WBPH fungal endosymbiont mitogenomes. Phylogenetic analysis based on conserved PCGs across 25 Ophiocordycipitaceae mitogenomes revealed that the five mitogenomes were clustered with that of R. speculum, forming an independent clade. In addition to providing the full mitogenome sequences, obtaining complete mitogenomes of WBPH endosymbionts can provide insights into their phylogenetic positions without needing to isolate the mtDNA from the host. This advantage is of value to future studies involving fungal endosymbiont mitogenomes.

Revealing the high variability on nonconserved core and mobile elements of Austropuccinia psidii and other rust mitochondrial genomes

Mitochondrial genomes are highly conserved in many fungal groups, and they can help characterize the phylogenetic relationships and evolutionary biology of plant pathogenic fungi. Rust fungi are among the most devastating diseases for economically important crops around the world. Here, we report the complete sequence and annotation of the mitochondrial genome of Austropuccinia psidii (syn. Puccinia psidii), the causal agent of myrtle rust. We performed a phylogenomic analysis including the complete mitochondrial sequences from other rust fungi. The genome composed of 93.299 bp has 73 predicted genes, 33 of which encoded nonconserved proteins (ncORFs), representing almost 45% of all predicted genes. A. psidii mtDNA is one of the largest rust mtDNA sequenced to date, most likely due to the abundance of ncORFs. Among them, 33% were within intronic regions of diverse intron groups. Mobile genetic elements invading intron sequences may have played significant roles in size but not shaping of the rust mitochondrial genome structure. The mtDNAs from rust fungi are highly syntenic. Phylogenetic inferences with 14 concatenated mitochondrial proteins encoded by the core genes placed A. psidii according to phylogenetic analysis based on 18S rDNA. Interestingly, cox1, the gene with the greatest number of introns, provided phylogenies not congruent with the core set. For the first time, we identified the proteins encoded by three A. psidii ncORFs using proteomics analyses. Also, the orf208 encoded a transmembrane protein repressed during in vitro morphogenesis. To the best of our knowledge, we presented the first report of a complete mtDNA sequence of a member of the family Sphaerophragmiacea.

The Ubiquity and Development-Related Abundance Dynamics of Ophiocordyceps Fungi in Soft Scale Insects

Mutual relationships with symbionts play a crucial role in the evolution and ecology of plant-feeding hemipteran insects. However, there was no specific dominant bacterium observed in soft scales (Coccidae) in the previous studies, it is still unclear whether soft scales have specific primary symbionts. In this study, a nuclear ribosomal internal transcribed spacer (ITS)gene fragment was used to analyze the diversity of fungal communities in 28 Coccidae species based on next-generation sequencing (NGS). Furthermore, samples from different developmental stages of Ceroplastes japonicus were sequenced to illustrate the dynamics of fungal community. Our results showed that Coccidae-associated Ophiocordyceps fungi (COF) were prevalent in all 28 tested species with high relative abundance. Meanwhile, the first and second instars of C. japonicus, two important stages for growth and development, had high relative abundance of COF, while the relative abundances in other stages were low, ranging from 0.68% to 2.07%. The result of fluorescent in situ hybridization showed that the COF were widely present in hemolymph and vertically transmitted from mother to offspring. Our study confirms that the COF have intimate associations with the growth and development of soft scales, and provides new evidence to support that COF are primary fungal symbionts for Coccidae.

Chinese caterpillar fungus (Ophiocordyceps sinensis) in China: Current distribution, trading, and futures under climate change and overexploitation

Chinese caterpillar fungus (Ophiocordyceps sinensis) is a precious traditional medicine which is mostly distributed on the Qinghai-Tibetan Plateau (QTP). Due to its medicinal values, it has become one of the most valuable biological commodities and widely traded in recent years worldwide. However, its habitat has changed profoundly in recent years under global warming as well as anthropogenic pressures, resulting in a sharp decline in its wild population in recent years. Based on the occurrence samples, this paper estimates the potential distribution of caterpillar fungus using MaxEnt model. The model simulates potential geographical distribution of the species under current climate conditions, and examine future distributions under different climatic change scenarios (i.e., RCP 2.6, RCP 4.5, RCP 6.0 and RCP 8.5 have been modelled in 2050s and 2070s, respectively). For examining the impacts of climate change in future, the integrated effects of climatic impact, trading, and overexploitation had been analyzed in detailed routes. The results show that: 1) The distribution patterns of caterpillar fungus under scenario RCP 2.6 have been predicted without obvious changes. However, range shift has been observed with significant shrinks across all classes of suitable areas in Tianshan, Kunlun Mountains, and the southwestern QTP in 2050s and 2070s under RCP 4.5, RCP 6.0 and RCP 8.5 scenarios, respectively. 2) The exports were decreasing drastically in recent years. Guangzhou and Hongkong are two international super import and consumption centres of caterpillar fungus in the world. 3) Both ecological and economic sustainable utilization of the caterpillar fungus has been threatened by the combined pressures of climate change and overexploitation. A strict but effective regulation and protection system, even a systematic management plan not just on the collectors but the whole explore process are urgently needed and has to be issued in the QTP.

Genomic analyses reveal evolutionary and geologic context for the plateau fungus Ophiocordyceps sinensis

Background

Ophiocordyceps sinensis, which is only naturally found in the high-elevation extreme environment of the Tibetan Plateau, has been used in traditional Chinese medicine. Information concerning the evolutionary and geologic context of O. sinensis remains limited, however.

Methods

We constructed the high-quality genome of O. sinensis and provided insight into the evolution and ecology of O. sinensis using comparative genomics.

Results

We mapped the whole genome of the anamorph/asexual form Hirsutella of O. sinensis using Illumina and PacBio sequencing technologies and obtained a well assembled genome of 119.2 Mbp size. Long-read Single Molecule Real Time (SMRT) sequencing technology generated an assembly with more accurate representation of repeat sequence abundances and placement. Evolutionary analyses indicated that O. sinensis diverged from other fungi 65.9 Mya in the Upper Cretaceous, during the uplift of the Tibetan Plateau. Gene family expansions and contractions in addition to genome inflation via long terminal repeat (LTR) retrotransposon insertions were implicated as an important driver of O. sinensis divergence. The insertion rate of LTR sequences into the O. sinensis genome peaked ~ 30-40 Mya, when the Tibetan Plateau rose rapidly. Gene Ontology (GO) enrichment analysis suggested that O. sinensis contained more genes related to ice binding compared to other closely related fungi, which may aid in their adaptability to the cold Tibetan Plateau. Further, heavy metal resistance genes were in low abundance in the O. sinensis genome, which may help to explain previous observations that O. sinensis tissues contain high levels of heavy metals.

Conclusions

Our results reveal the evolutionary, geological, and ecological context for the evolution of the O. sinensis genome and the factors that have contributed to the environmental adaptability of this valuable fungus. These findings suggest that genome inflation via LTR retrotransposon insertions in O. sinensis coincided with the uplift of the Tibetan Plateau. LTRs and the specific genetic mechanisms of O. sinensis contributed to its adaptation to the environment on the plateau.

Characterization of Two Mitochondrial Genomes and Gene Expression Analysis Reveal Clues for Variations, Evolution, and Large-Sclerotium Formation in Medical Fungus Wolfiporia cocos

Wolfiporia cocos, a precious mushroom with a long history as an edible food and Asian traditional medicine, remains unclear in the genetic mechanism underlying the formation of large sclerotia. Here, two complete circular mitogenomes (BL16, 135,686 bp and MD-104 SS10, 124,842 bp, respectively) were presented in detail first. The salient features in the mitogenomes of W. cocos include an intron in the tRNA (trnQ-UUG²), and an obvious gene rearrangement identified between the two mitogenomes from the widely geographically separated W. cocos strains. Genome comparison and phylogenetic analyses reveal some variations and evolutional characteristics in W. cocos. Whether the mitochondrion is functional in W. cocos sclerotium development was investigated by analyzing the mitogenome synteny of 10 sclerotium-forming fungi and mitochondrial gene expression patterns in different W. cocos sclerotium-developmental stages. Three common homologous genes identified across ten sclerotium-forming fungi were also found to exhibit significant differential expression levels during W. cocos sclerotium development. Most of the mitogenomic genes are not expressed in the mycelial stage but highly expressed in the sclerotium initial or developmental stage. These results indicate that some of mitochondrial genes may play a role in the development of sclerotium in W. cocos, which needs to be further elucidated in future studies. This study will stimulate new ideas on cytoplasmic inheritance of W. cocos and facilitate the research on the role of mitochondria in large sclerotium formation.

Intra-specific comparison of mitochondrial genomes reveals host gene fragment exchange via intron mobility in Tremella fuciformis

Background

Mitochondrial genomic sequences are known to be variable. Comparative analyses of mitochondrial genomes can reveal the nature and extent of their variation.

Results

Draft mitochondrial genomes of 16 Tremella fuciformis isolates (TF01-TF16) were assembled from Illumina and PacBio sequencing data. Mitochondrial DNA contigs were extracted and assembled into complete circular molecules, ranging from 35,104 bp to 49,044 bp in size. All mtDNAs contained the same set of 41 conserved genes with identical gene order. Comparative analyses revealed that introns and intergenic regions were variable, whereas genic regions (including coding sequences, tRNA, and rRNA genes) were conserved. Among 24 introns detected, 11 were in protein-coding genes, 3 in tRNA genes, and the other 10 in rRNA genes. In addition, two mobile fragments were found in intergenic regions. Interestingly, six introns containing N-terminal duplication of the host genes were found in five conserved protein-coding gene sequences. Comparison of genes with and without these introns gave rise to the following proposed model: gene fragment exchange with other species can occur via gain or loss of introns with N-terminal duplication of the host genes.

Conclusions

Our findings suggest a novel mechanism of fungal mitochondrial gene evolution: partial foreign gene replacement though intron mobility.

Fungal Mitogenomes: Relevant Features to Planning Plant Disease Management

Mitochondrial genomes (mt-genomes) are characterized by a distinct codon usage and their autonomous replication. Mt-genomes encode highly conserved genes (mt-genes), like proteins involved in electron transport and oxidative phosphorylation but they also carry highly variable regions that are in part responsible for their high plasticity. The degree of conservation of their genes is such that they allow the establishment of phylogenetic relationships even across distantly related species. Here, we describe the mechanisms that generate changes along mt-genomes, which play key roles at enlarging the ability of fungi to adapt to changing environments. Within mt-genomes of fungal pathogens, there are dispensable as well as indispensable genes for survival, virulence and/or pathogenicity. We also describe the different complexes or mechanisms targeted by fungicides, thus addressing a relevant issue regarding disease management. Despite the controversial origin and evolution of fungal mt-genomes, the intrinsic mechanisms and molecular biology involved in their evolution will help to understand, at the molecular level, the strategies for fungal disease management.

Exploring the Relationship Among Divergence Time and Coding and Non-coding Elements in the Shaping of Fungal Mitochondrial Genomes

The order Hypocreales (Ascomycota) is composed of ubiquitous and ecologically diverse fungi such as saprobes, biotrophs, and pathogens. Despite their phylogenetic relationship, these species exhibit high variability in biomolecules production, lifestyle, and fitness. The mitochondria play an important role in the fungal biology, providing energy to the cells and regulating diverse processes, such as immune response. In spite of its importance, the mechanisms that shape fungal mitogenomes are still poorly understood. Herein, we investigated the variability and evolution of mitogenomes and its relationship with the divergence time using the order Hypocreales as a study model. We sequenced and annotated for the first time Trichoderma harzianum mitochondrial genome (mtDNA), which was compared to other 34 mtDNAs species that were publicly available. Comparative analysis revealed a substantial structural and size variation on non-coding mtDNA regions, despite the conservation of copy number, length, and structure of protein-coding elements. Interestingly, we observed a highly significant correlation between mitogenome length, and the number and size of non-coding sequences in mitochondrial genome. Among the non-coding elements, group I and II introns and homing endonucleases genes (HEGs) were the main contributors to discrepancies in mitogenomes structure and length. Several intronic sequences displayed sequence similarity among species, and some of them are conserved even at gene position, and were present in the majority of mitogenomes, indicating its origin in a common ancestor. On the other hand, we also identified species-specific introns that advocate for the origin by different mechanisms. Investigation of mitochondrial gene transfer to the nuclear genome revealed that nuclear copies of the nad5 are the most frequent while atp8, atp9, and cox3 could not be identified in any of the nuclear genomes analyzed. Moreover, we also estimated the divergence time of each species and investigated its relationship with coding and non-coding elements as well as with the length of mitogenomes. Altogether, our results demonstrated that introns and HEGs are key elements on mitogenome shaping and its presence on fast-evolving mtDNAs could be mostly explained by its divergence time, although the intron sharing profile suggests the involvement of other mechanisms on the mitochondrial genome evolution, such as horizontal transference.

Complete Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma atroviride: Genomic Features, Comparative Analysis and Insight Into the Mitochondrial Evolution in Trichoderma

The improvement of biopesticides for use in the agriculture industry requires an understanding of the biological- and ecological principles underlying their behavior in natural environments. The nuclear genomes of members of the genus Trichoderma, which are representative fungal biocontrol agents, have been actively studied in relation to the unique characteristics of these species as effective producers of CAZymes/secondary metabolites and biopesticides, but their mitochondrial genomes have received much less attention. In this study, the mitochondrial genome of Trichoderma atroviride (Hypocreales, Sordariomycetes), which targets wood-decaying fungal pathogens and has the ability to degrade chemical fungicides, was assembled de novo. A 32,758 bp circular DNA molecule was revealed with specific features, such as a few more protein CDS and trn genes, two homing endonucleases (LAGLIDADG-/GIY-YIG-type), and even a putative overlapping tRNA gene, on a closer phylogenetic relationship with T. gamsii among hypocrealean fungi. Particularly, introns were observed with several footprints likely to be evolutionarily associated with the intron dynamics of the Trichoderma mitochondrial genomes. This study is the first to report the complete de novo mitochondrial genome of T. atroviride, while comparative analyses of Trichoderma mitochondrial genomes were also conducted from the perspective of mitochondrial evolution for the first time.

Mitochondrial genome and diverse inheritance patterns in Pleurotus pulmonarius

Pleurotus pulmonarius, a member of the Pleurotaceae family in Basidiomycota, is an edible, economically important mushroom in most Asian countries. In this study, the complete mitochondrial genomes (mtDNA) of three P. pulmonarius strains - two monokaryotic commercial (J1-13 and ZA3) and one wild (X1-15) - were sequenced and analyzed. In ZA3 and X1-15, the mtDNA molecule was found to be a single circle of 68,305 bp and 73,435 bp, respectively. Both strains contain 14 core protein-coding genes and two ribosomal RNA (rRNA) subunit genes. The ZA3 strain has 22 transfer RNA (tRNA) genes and nine introns: eight in cytochrome c oxidase subunit 1 (coxl), and one in the rRNA large subunit (rnl). Monokaryotic J1-13 and ZA3 mtDNAs were found to be similar in their structure. However, the wild strain X1-15 contains 25 tRNA genes and only seven introns in coxl. Open reading frames (ORFs) of ZA3/J1-13 and X1-15 encode LAGLIDADG, ribosomal protein S3, and DNA polymerase II. In addition, mtDNA inheritance in J1-13, ZA3, and X1-15 was also studied. Results showed that the mtDNA inheritance pattern was uniparental and closely related to dikaryotic hyphal location with respect to the parent. Results also show that mtDNA inheritance is influenced by both the parental nuclear genome and mitogenome in the zone of contact between two compatible parents. In summary, this analysis provides valuable information and a basis for further studies to improve our understanding of the inheritance of fungal mtDNA.

Subchromosome-Scale Nuclear and Complete Mitochondrial Genome Characteristics of Morchella crassipes

Morchella crassipes (Vent.) Pers., a typical yellow morel species with high economic value, is mainly distributed in the low altitude plains of Eurasia. However, rare research has been performed on its genomics and polarity, thus limiting its research and development. Here, we reported a fine physical map of the nuclear genome at the subchromosomal-scale and the complete mitochondrial genome of M. crassipes. The complete size of the nuclear genome was 56.7 Mb, and 23 scaffolds were assembled, with eight of them being complete chromosomes. A total of 11,565 encoding proteins were predicted. The divergence time analysis showed that M. crassipes representing yellow morels differentiated with black morels at ~33.98 Mya (million years), with 150 gene families contracted and expanded in M. crassipes versus the two black morels (M. snyderi and M. importuna). Furthermore, 409 CAZYme genes were annotated in M. crassipes, containing almost all plant cell wall degrading enzymes compared with the mycorrhizal fungi (truffles). Genomic annotation of mating type loci and amplification of the mating genes in the monospore population was conducted, the results indicated that M. crassipes is a heterothallic fungus. Additionally, a complete circular mitochondrial genome of M. crassipes was assembled, the size reached as large as 531,195 bp. It can be observed that the strikingly large size was the biggest up till now, coupled with 14 core conserved mitochondrial protein-coding genes, two rRNAs, 31 tRNAs, 51 introns, and 412 ncORFs. The total length of intron sequences accounted for 53.67% of the mitochondrial genome, with 19 introns having a length over 5 kb. Particularly, 221 of 412 ncORFs were distributed within 51 introns, and the total length of the ncORFs sequence accounted for 40.83% of the mitochondrial genome, and 297 ncORFs had expression activity in the mycelium stage, suggesting their potential functions in M. crassipes. Meanwhile, there was a high degree of repetition (51.31%) in the mitochondria of M. crassipes. Thus, the large number of introns, ncORFs and internal repeat sequences may contribute jointly to the largest fungal mitochondrial genome to date. The fine physical maps of nuclear genome and mitochondrial genome obtained in this study will open a new door for better understanding of the mysterious species of M. crassipes.

Revealing mitogenome-wide DNA methylation and RNA editing of three Ascomycotina fungi using SMRT sequencing

Beauveria bassiana, Cordyceps militaris and Ophiocordyceps sinensis (Ascomycotina) are traditional Chinese medicines. Here, mitogenomes of these three Ascomycotina fungi were sequenced and de-novo assembled using single-molecule real-time sequencing. The results showed that their complete mitogenomes were 31,258, 31,854 and 157,584 bp, respectively, with sequencing depth approximately 278,760×, 326,283× and 69,385×. Types of repeat sequences were mainly (AA)n, (AAT)n, (TA)n and (TATT)n. DNA methylation motifs were revealed in DNA modifications of these three fungi. We discovered new models of RNA editing through analysis of transcriptomes from B. bassiana and C. militaris. These data lay a solid foundation for further genetic and biological studies about these three fungi, especially for elucidating the mitogenome evolution and exploring the regulatory mechanism of adapting environment.

Artificial Cultivation of the Chinese Cordyceps From Injected Ghost Moth Larvae

The Chinese cordyceps, regarded as the 'Himalayan Viagra', is highly valued for its medicinal benefits. The decline of its yield due to over-exploitation and increased market demand have stimulated efforts to artificially cultivate Chinese cordyceps for over half a century. However, successful cultivation of Chinese cordyceps through caterpillar infection by the fungus Ophiocordyceps sinensis (Berk.) and the induction of the fruiting body from each mummified cadaver remains difficult for its complex life cycle. Herein, we report the developmental dynamics of hyphal bodies in hemolymph of injected Thitarodes xiaojinensis (Tu, Ma & Zhang) larvae and the success in artificial cultivation of sexual fruiting bodies from the mummified cadavers in the low-altitude area. We find that not only the numbers of hyphal bodies but also the conversion of hyphal bodies into hyphae played important roles in the mummification of the injected larvae. This cultivation will be beneficial for sustainable utilization of natural resources and provides the possibility for further research on the mechanism of the interaction between pathogenic fungus and host insect.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

Background

Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.

Results

Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.

Conclusions

Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Electronic supplementary material

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

Comparative Characterization of the Complete Mitochondrial Genomes of the Three Apple Snails (Gastropoda: Ampullariidae) and the Phylogenetic Analyses

The apple snails Pomacea canaliculata, Pomacea diffusa and Pomacea maculate (Gastropoda: Caenogastropoda: Ampullariidae) are invasive pests causing massive economic losses and ecological damage. We sequenced and characterized the complete mitochondrial genomes of these snails to conduct phylogenetic analyses based on comparisons with the mitochondrial protein coding sequences of 47 Caenogastropoda species. The gene arrangements, distribution and content were canonically identical and consistent with typical Mollusca except for the tRNA-Gln absent in P. diffusa. An identifiable control region (d-loop) was absent. Bayesian phylogenetic analysis indicated that all the Ampullariidae species clustered on the same branch. The genus Pomacea clustered together and then with the genus Marisa. The orders Architaenioglossa and Sorbeoconcha clustered together and then with the order Hypsogastropoda. Furthermore, the intergenic and interspecific taxonomic positions were defined. Unexpectedly, Ceraesignum maximum, Dendropoma gregarium, Eualetes tulipa and Thylacodes squamigerus, traditionally classified in order Hypsogastropoda, were isolated from the order Hypsogastropoda in the most external branch of the Bayesian inference tree. The divergence times of the Caenogastropoda indicated that their evolutionary process covered four geological epochs that included the Quaternary, Neogene, Paleogene and Cretaceous periods. This study will facilitate further investigation of species identification to aid in the implementation of effective management and control strategies of these invasive species.

A breakthrough in the artificial cultivation of Chinese cordyceps on a large-scale and its impact on science, the economy, and industry

Chinese cordyceps, an entity of the Chinese caterpillar fungus (Ophiocordyceps sinensis, syn. Cordyceps sinensis) that parasitizes ghost moth larvae, is one of the best known traditional Chinese medicines and is found exclusively on the Tibetan Plateau with limited natural resources. Although the fungus O. sinensis can grow on artificial substrates and the ghost moth has been successfully reared, the large-scale artificial cultivation of Chinese cordyceps has only recently been accomplished after several decades of efforts and attempts. In this article, research progress related to this breakthrough from living habitats, the life history of the fungus, its host insect, fungal isolation and culture, host larvae rearing, infection cycle of the fungus to the host, primordium induction, and fruiting body development have been reviewed. An understanding of the basic biology of O. sinensis, its host insect and the simulation of the Tibetan alpine environment resulted in the success of artificial cultivation on a large scale. Practical workshop production has reached annual yields of 2.5, 5, and 10 tons in 2014, 2015, and 2016, respectively. There was no difference in the chemical components detected between the cultivated and natural Chinese cordyceps. However, the artificial cultivation system can be controlled to avoid heavy metal contamination and results in high-quality products. Although omics studies, including genomic, transcriptomic, proteomic, and metabolomic studies, have helped to understand the biology of the fungus, the success of the artificial cultivation of the Chinese cordyceps is clearly a milestone and provides the possibility for research on the in-depth mechanisms of the interaction between the fungus and host insects and their adaptation to the harsh habitats. This cultivation will not only result in a large industry to alleviate the pressure of human demand but also protect the limited natural resources for sustainable utilization.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

Characterization of the Complete Mitochondrial Genome Sequences of Three Croakers (Perciformes, Sciaenidae) and Novel Insights into the Phylogenetics

The three croakers (Nibea coibor, Protonibea diacanthus and Argyrosomus amoyensis, Perciformes, Sciaenidae) are important commercial species inhabiting the Eastern Indian Ocean and Western Pacific. Molecular data employed in previous research on phylogenetic reconstruction have not been adequate and complete, and systematic and comprehensive phylogenetic relationships for these fish are unresolved. We sequenced the complete mitochondrial genomes of the three croakers using next-generation sequencing for the first time. We analyzed the composition and phylogenies between 19 species in the family Sciaenidae using the mitochondrial protein coding sequences of 204 species in the Series Eupercaria. We present the characterization of the complete mitochondrial genome sequences of the three croakers. Gene arrangement and distribution of the three croakers are canonically identical and consistent with other vertebrates. We found that the family Sciaenidae is an independent branch that is isolated from the order Perciformes and does not belong to any extant classification. Therefore, this family is expected to belong to a new classification at the order level and needs further analysis. The evolution of Sciaenidae has lagged far behind the Perciformes differentiation. This study presents a novel insight into the phylogenetics of the family Sciaenidae from the order Perciformes and facilitates additional studies on the evolution and phylogeny of Series Eupercaria.

Comparative analysis of the mitochondrial genome of the fungus Colletotrichum lindemuthianum, the causal agent of anthracnose in common beans

Fungi of the genus Colletotrichum are economically important and are used as models in plant-pathogen interaction studies. In this study, the complete mitochondrial genomes of two Colletotrichum lindemuthianum isolates were sequenced and compared with the mitochondrial genomes of seven species of Colletotrichum. The mitochondrial genome of C. lindemuthianum is a typical circular molecule 37,446 bp (isolate 89 A₂ 2-3) and 37,440 bp (isolate 83.501) in length. The difference of six nucleotides between the two genomes is the result of a deletion in the ribosomal protein S3 (rps3) gene in the 83.501 isolate. In addition, substitution of adenine for guanine within the rps3 gene in the mitochondrial genome of the 83.501 isolate was observed. Compared to the previously sequenced C. lindemuthianum mitochondrial genome, an exon no annotated in the cytochrome c oxidase I (cox1) gene and a non-conserved open reading frame (ncORF) were observed. The size of the mitochondrial genomes of the seven species of Colletotrichum was highly variable, being attributed mainly to the ncORF, ranging from one to 10 and also from introns ranging from one to 11 and which encode a total of up to nine homing endonucleases. This paper reports for the first time by means of transcriptome that then ncORFs are transcribed in Colletotrichum spp. Phylogeny data revealed that core mitochondrial genes could be used as an alternative in phylogenetic relationship studies in Colletotrichum spp. This work contributes to the genetic and biological knowledge of Colletotrichum spp., which is of great economic and scientific importance.

The complete mitochondrial genome of Engyodontium album and comparative analyses with Ascomycota mitogenomes

Engyodontium album is a widespread pathogen that causes different kinds of dermatoses and respiratory tract diseases in humans and animals. In spite of its perniciousness, the basic genetic and molecular background of this species remains poorly understood. In this study, the mitochondrial genome sequence of E. album was determined using a high-throughput sequencing platform. The circular mitogenome was found to be 28,081 nucleotides in length and comprised of 17 protein-coding genes, 24 tRNA genes, and 2 rRNA genes. The nucleotide composition of the genome was A+T-biased (74.13%). Group-II introns were found in the nad1, nad5, and cob genes. The most frequently used codon of protein-coding genes was UAU. Isoleucine was identified as the most common amino acid, while proline was the least common amino acid in protein-coding genes. The gene-arrangement order is nearly the same when compared with other Ascomycota mitogenomes. Phylogenetic relationships based on the shared protein-coding genes revealed that E. album is closely related to the Cordycipitaceae family, with a high-confidence support value (100%). The availability of the mitogenome of E. album will shed light on the molecular systematic and genetic differentiation of this species.

Comparative genomics of maize ear rot pathogens reveals expansion of carbohydrate-active enzymes and secondary metabolism backbone genes in Stenocarpella maydis

Stenocarpella maydis is a plant pathogenic fungus that causes Diplodia ear rot, one of the most destructive diseases of maize. To date, little information is available regarding the molecular basis of pathogenesis in this organism, in part due to limited genomic resources. In this study, a 54.8 Mb draft genome assembly of S. maydis was obtained with Illumina and PacBio sequencing technologies, and analyzed. Comparative genomic analyses with the predominant maize ear rot pathogens Aspergillus flavus, Fusarium verticillioides, and Fusarium graminearum revealed an expanded set of carbohydrate-active enzymes for cellulose and hemicellulose degradation in S. maydis. Analyses of predicted genes involved in starch degradation revealed six putative α-amylases, four extracellular and two intracellular, and two putative γ-amylases, one of which appears to have been acquired from bacteria via horizontal transfer. Additionally, 87 backbone genes involved in secondary metabolism were identified, which represents one of the largest known assemblages among Pezizomycotina species. Numerous secondary metabolite gene clusters were identified, including two clusters likely involved in the biosynthesis of diplodiatoxin and chaetoglobosins. The draft genome of S. maydis presented here will serve as a useful resource for molecular genetics, functional genomics, and analyses of population diversity in this organism.

SMRT Sequencing Revealed Mitogenome Characteristics and Mitogenome-Wide DNA Modification Pattern in Ophiocordyceps sinensis

Single molecule, real-time (SMRT) sequencing was used to characterize mitochondrial (mt) genome of Ophiocordyceps sinensis and to analyze the mt genome-wide pattern of epigenetic DNA modification. The complete mt genome of O. sinensis, with a size of 157,539 bp, is the fourth largest Ascomycota mt genome sequenced to date. It contained 14 conserved protein-coding genes (PCGs), 1 intronic protein rps3, 27 tRNAs and 2 rRNA subunits, which are common characteristics of the known mt genomes in Hypocreales. A phylogenetic tree inferred from 14 PCGs in Pezizomycotina fungi supports O. sinensis as most closely related to Hirsutella rhossiliensis in Ophiocordycipitaceae. A total of 36 sequence sites in rps3 were under positive selection, with dN/dS >1 in the 20 compared fungi. Among them, 16 sites were statistically significant. In addition, the mt genome-wide base modification pattern of O. sinensis was determined in this study, especially DNA methylation. The methylations were located in coding and uncoding regions of mt PCGs in O. sinensis, and might be closely related to the expression of PCGs or the binding affinity of transcription factor A to mtDNA. Consequently, these methylations may affect the enzymatic activity of oxidative phosphorylation and then the mt respiratory rate; or they may influence mt biogenesis. Therefore, methylations in the mitogenome of O. sinensis might be a genetic feature to adapt to the cold and low PO2 environment at high altitude, where O. sinensis is endemic. This is the first report on epigenetic modifications in a fungal mt genome.

Mitochondrial genome of the nematode endoparasitic fungus Hirsutella vermicola reveals a high level of synteny in the family Ophiocordycipitaceae

Ophiocordycipitaceae is a diverse fungal family comprising multiple ecologically, economically, medicinally, and culturally important fungal species; however, only four species of the family have available mitochondrial genomes (mitogenomes). In this study, the complete mitogenome of the nematode endoparasitic fungus Hirsutella vermicola in Ophiocordycipitaceae was sequenced, and a comparative mitogenomic analysis of Ophiocordycipitaceae was performed. We found that the 53,793-bp circular mitogenome of H. vermicola, except for standard fungal mitochondrial genes, harbors seven introns acquired possibly through lateral transfer from other fungi and three free-standing open reading frames (ORFs) coding for hypothetical proteins. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed its placement in Ophiocordycipitaceae. Comparison on five mitogenomes of Ophiocordycipitaceae revealed great variation on their sizes, from 35.2 kb in Tolypocladium ophioglossoides to 157.5 kb in Ophiocordyceps sinensis, mainly due to variable numbers of introns (from 7 to 54) as well as variable lengths of intergenic regions. The five mitogenomes, however, are highly syntenic to each other in terms of gene order, the presence of an intronic ORF encoding ribosomal protein S3 within rnl, and the nad2/nad3 joining pattern. Our study is the first report of the mitogenome of H. vermicola and has facilitated the understanding of mitogenome evolution of Ophiocordycipitaceae.

Systems Pharmacology-based strategy to screen new adjuvant for hepatitis B vaccine from Traditional Chinese Medicine Ophiocordyceps sinensis

Adjuvants are common component for many vaccines but there are still few licensed for human use due to low efficiency or side effects. The present work adopted Systems Pharmacology analysis as a new strategy to screen adjuvants from traditional Chinese medicine. Ophiocordyceps sinensis has been used for many years in China and other Asian countries with many biological properties, but the pharmacological mechanism has not been fully elucidated. First in this study, 190 putative targets for 17 active compounds in Ophiocordyceps sinensis were retrieved and a systems pharmacology-based approach was applied to provide new insights into the pharmacological actions of the drug. Pathway enrichment analysis found that the targets participated in several immunological processes. Based on this, we selected cordycepin as a target compound to serve as an adjuvant of the hepatitis B vaccine because the existing vaccine often fails to induce an effective immune response in many subjects. Animal and cellular experiments finally validated that the new vaccine simultaneously improves the humoral and cellular immunity of BALB/c mice without side effects. All this results demonstrate that cordycepin could work as adjuvant to hepatitis b vaccine and systems-pharmacology analysis could be used as a new method to select adjuvants.

Evolution of group I introns in Porifera: new evidence for intron mobility and implications for DNA barcoding

BACKGROUND

Mitochondrial introns intermit coding regions of genes and feature characteristic secondary structures and splicing mechanisms. In metazoans, mitochondrial introns have only been detected in sponges, cnidarians, placozoans and one annelid species. Within demosponges, group I and group II introns are present in six families. Based on different insertion sites within the cox1 gene and secondary structures, four types of group I and two types of group II introns are known, which can harbor up to three encoding homing endonuclease genes (HEG) of the LAGLIDADG family (group I) and/or reverse transcriptase (group II). However, only little is known about sponge intron mobility, transmission, and origin due to the lack of a comprehensive dataset. We analyzed the largest dataset on sponge mitochondrial group I introns to date: 95 specimens, from 11 different sponge genera which provided novel insights into the evolution of group I introns.

RESULTS

For the first time group I introns were detected in four genera of the sponge family Scleritodermidae (Scleritoderma, Microscleroderma, Aciculites, Setidium). We demonstrated that group I introns in sponges aggregate in the most conserved regions of cox1. We showed that co-occurrence of two introns in cox1 is unique among metazoans, but not uncommon in sponges. However, this combination always associates an active intron with a degenerating one. Earlier hypotheses of HGT were confirmed and for the first time VGT and secondary losses of introns conclusively demonstrated.

CONCLUSION

This study validates the subclass Spirophorina (Tetractinellida) as an intron hotspot in sponges. Our analyses confirm that most sponge group I introns probably originated from fungi. DNA barcoding is discussed and the application of alternative primers suggested.

Cordyceps collected from Bhutan, an appropriate alternative of Cordyceps sinensis

Natural Cordyceps collected in Bhutan has been widely used as natural Cordyceps sinensis, an official species of Cordyceps used as Chinese medicines, around the world in recent years. However, whether Cordyceps from Bhutan could be really used as natural C. sinensis remains unknown. Therefore, DNA sequence, bioactive components including nucleosides and polysaccharides in twelve batches of Cordyceps from Bhutan were firstly investigated, and compared with natural C. sinensis. Results showed that the fungus of Cordyceps from Bhutan was C. sinensis and the host insect belonged to Hepialidae sp. In addition, nucleosides and their bases such as guanine, guanosine, hypoxanthine, uridine, inosine, thymidine, adenine, and adenosine, as well as compositional monosaccharides, partial acid or enzymatic hydrolysates, molecular weights and contents of polysaccharides in Cordyceps from Bhutan were all similar to those of natural C. sinensis. All data suggest that Cordyceps from Bhutan is a rational alternative of natural C. sinensis, which is beneficial for the improvement of their performance in health and medicinal food areas.

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×.

Complete Mitochondrial Genome Sequence of the Pezizomycete Pyronema confluens

The complete mitochondrial genome of the ascomycete Pyronema confluens has been sequenced. The circular genome has a size of 191 kb and contains 48 protein-coding genes, 26 tRNA genes, and two rRNA genes. Of the protein-coding genes, 14 encode conserved mitochondrial proteins, and 31 encode predicted homing endonuclease genes.

The mitochondrial genome of the nematode endoparasitic fungus Hirsutella rhossiliensis

In this study, we report the complete mitochondrial genome of Hirsutella rhossiliensis (Ophiocordycipitaceae, Hypocreales, Ascomycota). We construct the mitochondrial DNA genome organization of 62 483 bp in length of H. rhossiliensis by using the whole-genome resequencing method. Conserved genes including the large and small rRNA subunits, 26 tRNA and 14 protein-coding genes are identified. These protein-coding genes utilize ATG, GTG or TTG as initiation codons and TAA or TAG as termination codons. Moreover, we detect 10 group I introns and one unclassified intron in six genes (rnl, cob, cox1, cox3, nad1 and nad5) encoding ORFs of ribosomal protein S3 and GIY-YIG/LAGLIDADG endonucleases or hypothetical proteins. This mitochondrial genome will be useful in understanding the distribution and genetic diversity of this species.