Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus

The first two mitochondrial genomes from Apiotrichum reveal mitochondrial evolution and different taxonomic assignment of Trichosporonales

Apiotrichum is a diverse anamorphic basidiomycetous yeast genus, and its mitogenome characterization has not been revealed. In this study, we assembled two Apiotrichum mitogenomes and compared them with mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. The mitogenomes of Apiotrichum gracile and A. gamsii comprised circular DNA molecules, with sizes of 34,648 bp and 38,096 bp, respectively. Intronic regions were found contributed the most to the size expansion of A. gamsii mitogenome. Comparative mitogenomic analysis revealed that 6.85-38.89% of nucleotides varied between tRNAs shared by the two Apiotrichum mitogenomes. The GC content of all core PCGs in A. gamsii was lower than that of A. gracile, with an average low value of 4.97%. The rps3 gene differentiated the most among Agaricomycotina, Pucciniomycotina and Ustilaginomycotina species, while nad4L gene was the most conserved in evolution. The Ka/Ks values for cob and rps3 genes were > 1, indicating the two genes may be subjected to positive selection in Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. Frequent intron loss/gain events and potential intron transfer events have been detected in evolution of Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. We further detected large-scale gene rearrangements between the 19 mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina, and fifteen of the 17 mitochondrial genes shared by Apiotrichum varied in gene arrangements. Phylogenetic analyses based on maximum likelihood and Bayesian inference methods using a combined mitochondrial gene dataset revealed different taxonomic assignment of two Apiotrichum species, wherein A. gamsii had a more closely relationship with Trichosporon asahii. This study served as the first report on mitogenomes from the genus Apiotrichum, which promotes the understanding of evolution, genomics, and phylogeny of Apiotrichum.

The first two mitochondrial genomes for the genus Ramaria reveal mitochondrial genome evolution of Ramaria and phylogeny of Basidiomycota

In the present study, we assembled and analyzed the mitogenomes of two Ramaria species. The assembled mitogenomes of Ramaria cfr. rubripermanens and R. rubella were circularized, with sizes of 126,497 bp and 143,271 bp, respectively. Comparative mitogenome analysis showed that intron region contributed the most (contribution rate, 43.74%) to the size variations of Ramaria mitogenomes. The genetic contents, gene length, tRNAs, and codon usages of the two Ramaria mitogenomes varied greatly. In addition, the evolutionary rates of different core protein coding genes (PCGs) in Phallomycetidae mitogenomes varied. We detected large-scale gene rearrangements between Phallomycetidae mitogenomes, including gene displacement and tRNA doubling. A total of 4499 bp and 7746 bp aligned fragments were detected between the mitochondrial and nuclear genomes of R. cfr. rubripermanens and R. rubella, respectively, indicating possible gene transferring events. We further found frequent intron loss/gain and potential intron transfer events in Phallomycetidae mitogenomes during the evolution, and the mitogenomes of R. rubella contained a novel intron P44. Phylogenetic analyses using both Bayesian inference (BI) and Maximum Likelihood (ML) methods based on a combined mitochondrial gene dataset obtained an identical and well-supported phylogenetic tree for Basidiomycota, wherein R. cfr. rubripermanens and Turbinellus floccosus are sister species. This study served as the first report on mitogenomes from the genus Ramaria, which provides a basis for understanding the evolution, genetics, and taxonomy of this important fungal group.

Comparative Mitogenomic Analysis Reveals Intraspecific, Interspecific Variations and Genetic Diversity of Medical Fungus Ganoderma

Ganoderma species are widely distributed in the world with high diversity. Some species are considered to be pathogenic fungi while others are used as traditional medicine in Asia. In this study, we sequenced and assembled four Ganoderma complete mitogenomes, including G. subamboinense s118, G. lucidum s37, G. lingzhi s62, and G. lingzhi s74. The sizes of the four mitogenomes ranged from 50,603 to 73,416 bp. All Ganoderma specimens had a full set of core protein-coding genes (PCGs), and the rps3 gene of Ganoderma species was detected to be under positive or relaxed selection. We found that the non-conserved PCGs, which encode RNA polymerases, DNA polymerases, homing endonucleases, and unknown functional proteins, are dynamic within and between Ganoderma species. Introns were thought to be the main contributing factor in Ganoderma mitogenome size variation (p < 0.01). Frequent intron loss/gain events were detected within and between Ganoderma species. The mitogenome of G. lucidum s26 gained intron P637 in the cox3 gene compared with the other two G. lucidum mitogenomes. In addition, some rare introns in Ganoderma were detected in distinct Basidiomycetes, indicating potential gene transfer events. Comparative mitogenomic analysis revealed that gene arrangements also varied within and between Ganoderma mitogenomes. Using maximum likelihood and Bayesian inference methods with a combined mitochondrial gene dataset, phylogenetic analyses generated identical, well-supported tree topologies for 71 Agaricomycetes species. This study reveals intraspecific and interspecific variations of the Ganoderma mitogenomes, which promotes the understanding of the origin, evolution, and genetic diversity of Ganoderma species.

Genetic Differentiation and Widespread Mitochondrial Heteroplasmy among Geographic Populations of the Gourmet Mushroom Thelephora ganbajun from Yunnan, China

The mitochondrial genomes are generally considered non-recombining and homoplasmic in nature. However, our previous study provided the first evidence of extensive and stable mitochondrial heteroplasmy in natural populations of the basidiomycete fungus Thelephora ganbajun from Yunnan province, China. The heteroplasmy was characterized by the presence of two types of introns residing at adjacent but different sites in the cytochrome oxidase subunits I (cox1) gene within an individual strain. However, the frequencies of these two introns among isolates from different geographical populations and the implications for the genetic structure in natural populations have not been investigated. In this study, we analyzed DNA sequence variation at the internal transcribed spacer (ITS) regions of the nuclear ribosomal RNA gene cluster among 489 specimens from 30 geographic locations from Yunnan and compared that variation with distribution patterns of the two signature introns in the cox1 gene that are indicative of heteroplasmy in this species. In our samples, evidence for gene flow, abundant genetic diversity, and genotypic uniqueness among geographic samples in Yunnan were revealed by ITS sequence variation. While there was insignificant positive correlation between geographic distance and genetic differentiation among the geographic samples based on ITS sequences, a moderate significant correlation was found between ITS sequence variation, geographical distance of sampling sites, and distribution patterns of the two heteroplasmic introns in the cox1 gene. Interestingly, there was a significantly negative correlation between the copy numbers of the two co-existing introns. We discussed the implications of our results for a better understanding of the spread of stable mitochondrial heteroplasmy, mito-nuclear interactions, and conservation of this important gourmet mushroom.

First two mitochondrial genomes for the order Filobasidiales reveal novel gene rearrangements and intron dynamics of Tremellomycetes

In the present study, two mitogenomes from the Filobasidium genus were assembled and compared with other Tremellomycetes mitogenomes. The mitogenomes of F. wieringae and F. globisporum both comprised circular DNA molecules, with sizes of 27,861 bp and 71,783 bp, respectively. Comparative mitogenomic analysis revealed that the genetic contents, tRNAs, and codon usages of the two Filobasidium species differed greatly. The sizes of the two Filobasidium mitogenomes varied greatly with the introns being the main factor contributing to mitogenome expansion in F. globisporum. Positive selection was observed in several protein-coding genes (PCGs) in the Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina species, including cob, cox2, nad2, and rps3 genes. Frequent intron loss/gain events were detected to have occurred during the evolution of the Tremellomycetes mitogenomes, and the mitogenomes of 17 species from Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina have undergone large-scale gene rearrangements. Phylogenetic analyses based on Bayesian inference and the maximum likelihood methods using a combined mitochondrial gene set generated identical and well-supported phylogenetic trees, wherein Filobasidium species had close relationships with Trichosporonales species. This study, which is the first report on mitogenomes from the order Filobasidiales, provides a basis for understanding the genomics, evolution, and taxonomy of this important fungal group.

An Artificial Pathway for N-Hydroxy-Pipecolic Acid Production From L-Lysine in Escherichia coli

N-hydroxy-pipecolic acid (NHP) is a hydroxylated product of pipecolic acid and an important systemic acquired resistance signal molecule. However, the biosynthesis of NHP does not have a natural metabolic pathway in microorganisms. Here, we designed and constructed a promising artificial pathway in Escherichia coli for the first time to produce NHP from biomass-derived lysine. This biosynthesis route expands the lysine catabolism pathway and employs six enzymes to sequentially convert lysine into NHP. This artificial route involves six functional enzyme coexpression: lysine α-oxidase from Scomber japonicus (RaiP), glucose dehydrogenase from Bacillus subtilis (GDH), Δ1-piperideine-2-carboxylase reductase from Pseudomonas putida (DpkA), lysine permease from E. coli (LysP), flavin-dependent monooxygenase (FMO1), and catalase from E. coli (KatE). Moreover, different FMO1s are used to evaluate the performance of the produce NHP. A titer of 111.06 mg/L of NHP was yielded in shake flasks with minimal medium containing 4 g/L of lysine. By this approach, NHP has so far been produced at final titers reaching 326.42 mg/L by 48 h in a 5-L bioreactor. To the best of our knowledge, this is the first NHP process using E. coli and the first process to directly synthesize NHP by microorganisms. This study lays the foundation for the development and utilization of renewable resources to produce NHP in microorganisms.

Characterization of the complete mitochondrial genome of Aspergillus terricola (Aspergillaceae, Eurotiales), isolated from soy sauce fermentation system

In the present study, the complete mitochondrial genome of Aspergillus terricola É.J. Marchal 1893 was sequenced and assembled. The mitochondrial genome of A. terricola was composed of circular DNA molecules, with a total size of 28,689 bp. The GC content of the A. terricola mitochondrial genome was 26.34%. A total of 18 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, and 26 transfer RNA (tRNA) genes were detected in the A. terricola mitochondrial genome. Phylogenetic analysis based on the combined mitochondrial gene dataset indicated that the A. terricola exhibited a close relationship with A. parasiticus.

Characterization of the first complete chloroplast genome of Amaranthus hybridus (Caryophyllales: Amaranthaceae) with phylogenetic implications

In the present study, the complete chloroplast genome of Amaranthus hybridus was sequenced and assembled. The complete chloroplast genome of Amaranthus hybridus is 150,709 in size, with the GC content of 36.56%. The chloroplast genome of Amaranthus hybridus contained 86 protein-coding genes (PCGs), eight ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. Phylogenetic analysis based on combined chloroplast gene dataset indicated that the Amaranthus hybridus exhibited a close relationship with A. hypochondriacus and A. caudatus.

Coproduction of 5-Aminovalerate and δ-Valerolactam for the Synthesis of Nylon 5 From L-Lysine in Escherichia coli

The compounds 5-aminovalerate and δ-valerolactam are important building blocks that can be used to synthesize bioplastics. The production of 5-aminovalerate and δ-valerolactam in microorganisms provides an ideal source that reduces the cost. To achieve efficient biobased coproduction of 5-aminovalerate and δ-valerolactam in Escherichia coli, a single biotransformation step from L-lysine was constructed. First, an equilibrium mixture was formed by L-lysine α-oxidase RaiP from Scomber japonicus. In addition, by adjusting the pH and H₂O₂ concentration, the titers of 5-aminovalerate and δ-valerolactam reached 10.24 and 1.82 g/L from 40 g/L L-lysine HCl at pH 5.0 and 10 mM H₂O₂, respectively. With the optimized pH value, the δ-valerolactam titer was improved to 6.88 g/L at pH 9.0 with a molar yield of 0.35 mol/mol lysine. The ratio of 5AVA and δ-valerolactam was obviously affected by pH value. The ratio of 5AVA and δ-valerolactam could be obtained in the range of 5.63:1-0.58:1 at pH 5.0-9.0 from the equilibrium mixture. As a result, the simultaneous synthesis of 5-aminovalerate and δ-valerolactam from L-lysine in Escherichia coli is highly promising. To our knowledge, this result constitutes the highest δ-valerolactam titer reported by biological methods. In summary, a commercially implied bioprocess developed for the coproduction of 5-aminovalerate and δ-valerolactam using engineered Escherichia coli.

The first complete mitochondrial genome of mushroom Leucoagaricus naucinus (Agaricaceae, Agaricales) and insights into its phylogeny

Leucoagaricus naucinus (Fr.) Singer is a mycorrhizal fungus widely distributed in the northern Hemisphere. In the present study, the complete mitochondrial genome of Leucoagaricus naucinus was sequenced, assembled, and annotated. The L. naucinus mitochondrial genome was composed of circular DNA molecules, with the total size of 61,434 bp. The GC content of the L. naucinus mitochondrial genome was 26.07%. A total of 30 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 26 transfer RNA (tRNA) genes were detected in the L. naucinus mitochondrial genome. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the L. naucinus exhibited a close relationship with Agaricus bisporus.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

Induction of Apoptosis in HeLa Cells by a Novel Peptide from Fruiting Bodies of Morchella importuna via the Mitochondrial Apoptotic Pathway

Morels (Morchella spp.) are a genus of edible fungi with important economic and medicinal value. In this study, a novel peptide (MIPP) was extracted from the fruiting bodies of Morchella importuna using gel filtration chromatography. Structural analysis showed that the molecular mass of MIPP is 831 Da, and it has a simple amino acid sequence: Ser-Leu-Ser-Leu-Ser-Val-Ala-Arg. To explore the antitumor activity of MIPP, the effect of MIPP on HeLa cell apoptosis and the underlying preventative mechanisms were investigated. Results showed that MIPP reduced the viability of HeLa cells in a concentration-dependent manner. TUNEL analysis and flow cytometric examination showed that MIPP decreased cell proliferation via a mitochondrial-dependent pathway, as manifested by downregulation of Bcl-2/Bax, promotion of the movement of cytochrome C from the mitochondria to the cytoplasm, and triggering of caspase-9 and caspase-3. Therefore, MIPP may be a promising tumor-preventive agent, especially in human cervical cancer.

Characterization of the complete mitochondrial genome of Peniophora lycii (Russulales: Peniophoraceae) with its phylogenetic analysis

Peniophora lycii is a resupinate lichen-like species distributed all over the world. In the present study, we sequenced and assembled the complete mitochondrial genome of Peniophora lycii. The size of the mitochondrial genome of P. lycii was 38,296 bp, with a GC content of 25.89%. Twenty protein-coding genes, 2 ribosomal RNA genes, and 24 transfer RNA genes were identified in the mitochondrial genome of P. lycii. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the mitochondrial genome of P. lycii exhibited a close relationship with that of Heterobasidion irregulare.

The first complete mitochondrial genome of edible and medicinal fungus Chroogomphus rutilus (Gomphidiaceae, Boletales) and insights into its phylogeny

In the present study, we assembled and annotated the complete mitochondrial genome of Chroogomphus rutilus. The complete mitochondrial genome of C. rutilus was composed of circular DNA molecules, with a size of 37,508 bp. The GC content of the C. rutilus mitogenome was 22.82%. A total of 18 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, and 24 transfer RNA (tRNA) genes were detected in the C. rutilus mitogenome. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the C. rutilus exhibited a close relationship with species from the genus Rhizopogon. This study served as the first report on the complete mitochondrial genome from the family Gomphidiaceae, which will promote the understanding of phylogeny, evolution, and taxonomy of this important fungal species.

The first complete mitochondrial genome of Macalpinomyces bursus (Ustilaginales: Ustilaginaceae) and insights into its phylogeny

In the present study, the complete mitochondrial genome of Macalpinomyces bursus (Berk.) Vanky 2002 was sequenced and assembled. The complete mitochondrial genome of M. bursus was 49,024 bp in length, with the GC content of 30.4%. The M. bursus mitochondrial genome contained 27 protein-coding genes, 2 ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the M. bursus exhibited a close relationship with species from the genera Ustilago, Sporisorium, and Anthracocystis.

The first complete chloroplast genome of Fagopyrum leptopodum (Diels) Hedberg (Caryophyllales: Polygonaceae) with phylogenetic implications

In the present study, we sequenced and assembled the complete chloroplast genome of Fagopyrum leptopodum (Diels) Hedberg. The chloroplast genome of F. leptopodum was composed of 85 protein-coding genes, 8 ribosomal RNA genes, and 37 transfer RNA genes. The F. leptopodum chloroplast genome is 159,375 bp in length, with a GC content of 37.81%. Phylogenetic analysis based on the combined chloroplast gene dataset indicated that the F. leptopodum exhibited a close relationship with Fagopyrum luojishanense.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

Characterization of the complete chloroplast genome of Hordeum vulgare L. var. trifurcatum with phylogenetic analysis

In the present study, the complete chloroplast genome of Hordeum vulgare L. var. trifurcatum was sequenced, assembled and compared with closely related species. The chloroplast genome of Hordeum vulgare L. var. trifurcatum was composed of 84 protein-coding genes (PCG), 8 ribosomal RNA (rRNA) genes, and 38 transfer RNA (tRNA) genes. The Hordeum vulgare L. var. trifurcatum chloroplast genome is 136,485 bp in size, with the GC content of 38.32%. Phylogenetic analysis based on the combined chloroplast gene dataset indicated that the Hordeum vulgare L. var. trifurcatum exhibited a close relationship with Hordeum vulgare subsp. spontaneum and Hordeum vulgare subsp. vulgare.

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.

Genome sequence and transcriptome profiles of pathogenic fungus Paecilomyces penicillatus reveal its interactions with edible fungus Morchella importuna

Paecilomyces penicillatus is one of the pathogens of morels, which greatly affects the yield and quality of Morchella spp.. In the present study, we de novo assembled the genome sequence of the fungus P. penicillatus SAAS_ppe1. We analyzed the transcriptional profile of P. penicillatus SAAS_ppe1 infection of Morchella importuna at different stages (3 days and 6 days after infection) and the response of M. importuna using the transcriptome. The assembled genome sequence of P. penicillatus SAAS_ppe1 was 39.78 Mb in length (11 scaffolds; scaffold N50, 6.50 Mb), in which 99.7% of the expected genes were detected. A total of 7.48% and 19.83% clean transcriptional reads from the infected sites were mapped to the P. penicillatus genome at the early and late stages of infection, respectively. There were 3,943 genes differently expressed in P. penicillatus at different stages of infection, of which 24 genes had increased expression with the infection and infection stage, including diphthamide biosynthesis, aldehyde reductase, and NAD (P)H-hydrate epimerase (P < 0.05). Several genes had variable expression trends at different stages of infection, indicating P. penicillatus had diverse regulation patterns to infect M. importuna. GO function, involving cellular components, and KEGG pathways, involving glycerolipid metabolism, and plant-pathogen interaction were significantly enriched during infection by P. penicillatus. The expression of ten genes in M. importuna increased during the infection and infection stage, and these may regulate the response of M. importuna to P. penicillatus infection. This is the first comprehensive study on P. penicillatus infection mechanism and M. importuna response mechanism, which will lay a foundation for understanding the fungus-fungus interactions, gene functions, and variety breeding of pathogenic and edible fungi.