The complete mitochondrial genome of the nematode endoparasitic fungus Hirsutella minnesotensis

The complete mitochondrial genome of Linnemannia amoeboidea (W. Gams) Vandepol & Bonito (Mortierellales: Mortierellaceae)

The complete mitochondrial genome of Linnemannia amoeboidea (W. Gams) Vandepol & Bonito 2020 (Strain no.: CBS 889.72) was sequenced under the next-generation sequencing platform. It was the second one in the family Mortierellaceae Luerss. 1877. The circular genome was 49,702 bp in size, with a GC content of 20.86%. Gene prediction revealed 15 PCGs, two rRNA genes, 26 tRNA genes, one rnpB gene and seven ORFs. Phylogenetic analyses showed that L. amoeboidea was closely related to Podila verticillate (Linnem.) Vandepol & Bonito 2020.

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.

Mitogenomics, Phylogeny and Morphology Reveal Ophiocordyceps pingbianensis Sp. Nov., an Entomopathogenic Fungus from China

The new entomopathogenic fungus Ophiocordyceps pingbianensis, collected from Southeast China, was described by mitogenomic, morphological, and phylogenetic evidence. The systematic position of O. pingbianensis was determined by phylogenetic analyses based on six nuclear gene (ITS, tef1-α, nrSSU, nrLSU, rpb1 and rpb2) and 14 mitochondrial protein-coding gene (PCGs) (cox1, cox2, cox3, atp6, atp8, atp9, cob, nad1, nad2, nad3, nad4, nad5, nad6 and nad4L) data. Phylogenetic analyses reveal that O. pingbianensis was belonged to the Hirsutella nodulosa clade in the genus Ophiocordyceps of Ophiocordycipiaceae. This fungus exhibits distinctive characteristics which differed from other related Ophiocordyceps species with slender and geminate stromata, monophialidic conidiogenous cells with an inflated awl-shaped base, a twisty and warty phialide neck and a fusiform or oval conidia, as well as being found on a tiger beetle of Coleoptera buried in moss at the cave. The complete mitochondrial genome of O. pingbianensis was a circular DNA molecule 80,359 bp in length, containing 15 PCGs, 24 open reading frames genes (ORFs), 25 transfer RNA genes (tRNAs) and 27 introns. Ophiocordyceps pingbianensis, containing 27 introns, has the second largest mitogenome in Ophiocordycipiaceae and was next to O. sinensis. To our knowledge, this is the first report of the mitogenome from a new entomopathogenic fungus, and thus provides an important foundation for future studies on taxonomy, genetics and evolutionary biology of Ophiocordycipiaceae.

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.

Characterization of the Complete Mitochondrial Genome of Basidiomycete Yeast Hannaella oryzae: Intron Evolution, Gene Rearrangement, and Its Phylogeny

In this study, the mitogenome of Hannaella oryzae was sequenced by next-generation sequencing (NGS) and successfully assembled. The H. oryzae mitogenome comprised circular DNA molecules with a total size of 26,444 bp. We found that the mitogenome of H. oryzae partially deleted the tRNA gene transferring cysteine. Comparative mitogenomic analyses showed that intronic regions were the main factors contributing to the size variations of mitogenomes in Tremellales. Introns of the cox1 gene in Tremellales species were found to have undergone intron loss/gain events, and introns of the H. oryzae cox1 gene may have different origins. Gene arrangement analysis revealed that H. oryzae contained a unique gene order different from other Tremellales species. Phylogenetic analysis based on a combined mitochondrial gene set resulted in identical and well-supported topologies, wherein H. oryzae was closely related to Tremella fuciformis. This study represents the first report of mitogenome for the Hannaella genus, which will allow further study of the population genetics, taxonomy, and evolutionary biology of this important phylloplane yeast and other related species.

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.

Natural diversity in the predatory behavior facilitates the establishment of a robust model strain for nematode-trapping fungi

Nematode-trapping fungi (NTF) are a group of specialized microbial predators that consume nematodes when food sources are limited. Predation is initiated when conserved nematode ascaroside pheromones are sensed, followed by the development of complex trapping devices. To gain insights into the coevolution of this interkingdom predator-prey relationship, we investigated natural populations of nematodes and NTF that we found to be ubiquitous in soils. Arthrobotrys species were sympatric with various nematode species and behaved as generalist predators. The ability to sense prey among wild isolates of Arthrobotrys oligospora varied greatly, as determined by the number of traps after exposure to Caenorhabditis elegans While some strains were highly sensitive to C. elegans and the nematode pheromone ascarosides, others responded only weakly. Furthermore, strains that were highly sensitive to the nematode prey also developed traps faster. The polymorphic nature of trap formation correlated with competency in prey killing, as well as with the phylogeny of A. oligospora natural strains, calculated after assembly and annotation of the genomes of 20 isolates. A chromosome-level genome assembly and annotation were established for one of the most sensitive wild isolates, and deletion of the only G-protein β-subunit-encoding gene of A. oligospora nearly abolished trap formation. In summary, our study establishes a highly responsive A. oligospora wild isolate as a model strain for the study of fungus-nematode interactions and demonstrates that trap formation is a fitness character in generalist predators of the nematode-trapping fungus family.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

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.

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.

The complete mitochondrial genome of a medicinal fungus, Tolypocladium ophioglossoides

As part of a genome sequencing project for Tolypocladium ophioglossoides (T. ophioglossoides), a complete mitochondrial (mt) genome was assembled as a single circular dsDNA of 35,159 bp. Conserved genes including the large and small rRNA subunits, 25 tRNA and 15 protein-coding genes, were identified. In addition, 4 non-conserved open reading frames (ncORFs) in the intergenic and intronic regions were also identified. Transcription analyses using RNA-seq validated the expression of most conserved genes and ncORFs. Seven introns (groups I and II) were found within conserved genes, accounting for 21% of the mitogenome. All structural genes are located on one strand and are apparently transcribed in one direction. The complete mt genomes of T. ophioglossoides would be useful for future investigation of genetics, evolution and clinical identification of Tolypocladium species.

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.