Core Protein-coding Genes Research Articles

Assembly, Annotation, and Comparative Analysis of Mitochondrial Genomes in Trichoderma

Trichoderma is a widely studied ascomycete fungal genus, including more than 400 species. However, genetic information on Trichoderma is limited, with most species reporting only DNA barcodes. Mitochondria possess their own distinct DNA that plays a pivotal role in molecular function and evolution. Here, we report 42 novel mitochondrial genomes (mitogenomes) combined with 18 published mitogenomes of Trichoderma. These circular mitogenomes exhibit sizes of 26,276–94,608 bp, typically comprising 15 core protein-coding genes (PCGs), 2 rRNAs, and 16–30 tRNAs; however, the number of endonucleases and hypothetical proteins encoded in the introns of PCGs increases with genome size enlargement. According to the result of phylogenetic analysis of the whole mitogenome, these strains diverged into six distinct evolutionary branches, supported by the phylogeny based on 2830 single-copy nuclear genes. Comparative analysis revealed that dynamic Trichoderma mitogenomes exhibited variations in genome size, gene number, GC content, tRNA copy, and intron across different branches. We identified three mutation hotspots near the regions encoding nad3, cox2, and nad5 that caused major changes in the mitogenomes. Evolutionary analysis revealed that atp9, cob, nad4L, nad5, and rps3 have been influenced by positive selection during evolution. This study provides a valuable resource for exploring the important roles of the genetic and evolutionary dynamics of Trichoderma mitogenome in the adaptive evolution of biocontrol fungi.

Comparative Mitogenomics Provides Valuable Insights for the Phylogeny and New DNA Barcodes of Ganoderma

Ganoderma is the most important genus in the family Ganodermataceae; many species have attracted much attention and widely cultivated because of their medicinal values, but so far, not a sequenced mitogenome derived from dikaryon strains has been explicitly recorded. Herein, four novel mitogenomes of commonly cultivated Ganoderma (G. leucocontextum H4, G. lucidum G6, G. sinense MZ96 and G. tsugae SS) were de novo assembled and given detail functional annotations. Collinearity analysis revealed that the four mitogenomes shared 82.93–92.02% similarity with their corresponding reference mitogenomes at the nucleotide level. A total of 15 core protein-coding genes (PCGs), along with rrnL and rrnS (mtLSU and mtSSU) were chosen as potential candidates for constructing their individual phylogenetic trees. These trees were compared with those derived from the concatenated sequences of 15 core PCGs. And finally, we found that the atp9 and nad4L were the most reliable markers for the phylogenetic analysis of Ganoderma and chosen as standard sequences to generate new DNA barcodes. This finding was further verified by comparing it against almost all available Ganoderma mitogenomes in the NCBI, with Trametes versicolor (Polyporaceae) and Rigidoporus microporus (Meripilaceae) as two outgroups. A total of 52 mitogenomes from three families were highly conserved, with identical gene lengths for atp9 (222 bp) and nad4L (267 bp). These genes were capable of distinguish distinctly different various species, which are grouped into separate clades within the phylogenetic trees. The closest related clades (I and II), including at least 30 samples of the three classical taxonomic species (G. lingzhi, G. sichuanense and G. lucidum), differed in only one SNP. The single base mutation rate increased with the evolutionary divergence of the phylogenetic clades, from two to three SNPs in earlier clades (e.g., clade IV containing G. leucocontextum) to five to six SNPs in later clades (e.g., clade X containing G. sinense). Despite these variations between species, the atp9 and nad4L genes of Ganoderma mitogenomes consistently encoded the same ATP synthase F0 subunit c (73 aa) and NADH dehydrogenase subunit 4L (88 aa). These two genes have been identified as reliable markers of new DNA barcodes, offering valuable insights and contributing significantly to understanding the evolutionary relationships and phylogeny of the Ganoderma genus and even the Ganodermataceae family.

A pan-genomic analysis based multi-epitope vaccine development by targeting Stenotrophomonas maltophilia using reverse vaccinology method: an in-silico approach.

Antibiotic resistance in bacteria leads to high mortality rates and healthcare costs, a significant concern for public health. A colonizer of the human respiratory system, Stenotrophomonas maltophilia is frequently associated with hospital-acquired infections in individuals with cystic fibrosis, cancer, and other chronic illnesses. The importance of this study is underscored by its capacity to meet the critical demand for effective preventive strategies against this pathogen, particularly among susceptible groups of cystic fibrosis and those undergoing cancer treatment. In this study, we engineered a multi-epitope vaccine targeting S. maltophilia through genomic analysis, reverse vaccination strategies, and immunoinformatic techniques by examining a total of 81 complete genomes of S. maltophilia strains. Our investigation revealed 1945 core protein-coding genes alongside their corresponding proteomic sequences, with 191 of these genes predicted to exhibit virulence characteristics. Out of the filtered proteins, three best antigenic proteins were selected for epitope prediction while seven epitopes each from CTL, HTL, and B cell were chosen for vaccine development. The vaccine was refined and validated, showing highly antigenic and desirable physicochemical features. Molecular docking assessments revealed stable binding with TLR-4. Molecular dynamic simulation demonstrated stable dynamics with minor alterations. The originality of this investigation is rooted in the thorough techniques aimed at designing a vaccine that directly targets S. maltophilia, a microorganism of considerable clinical relevance that currently lacks an available vaccine. This study not only responds to a pressing public health crisis but also lays the groundwork for subsequent research endeavors focused on the prevention of S. maltophilia outbreaks. Further evidence from studies in mice models is needed to confirm immune protection against S. maltophilia.

Characterization of the complete mitochondrial genome of the medical fungus Ganoderma resinaceum Boud., 1889 (Polyporales: Ganodermataceae)

The medical mushroom Ganoderma resinaceum Boud., 1889, is of great interest in pharmacy due to its diverse functional active ingredients. However, the mitochondrial genome of G. resinaceum remains unexplored. Here, we present the complete mitochondrial genome of G. resinaceum, which spans 67,458 bp and has a GC content of 25.65%. This genome encompasses 15 core protein-coding genes, 8 independent ORFs, 15 intronic ORFs, 27 tRNAs, and 2 rRNA genes. Through phylogenetic analysis using Bayesian inference (BI), we elucidated the evolutionary relationships among 34 Basidiomycota fungi, revealing distinct clades and indicating a close relationship between G. resinaceum and G. subamboinense.

Characterization and phylogenetic analysis of the Talaromyces liani (kamyschko) Yilmaz, Frisvad & Samson, 2014 (Eurotiales: trichocomaceae) mitochondrial genome

The filamentous fungus Talaromyces liani (Kamyschko) Yilmaz, Frisvad & Samson, 2014, has attracted considerable interest in biotechnology due to its diverse industrial applications and physiological characteristics. However, the mitochondrial genome of T. liani remains uncharacterized. Here, we present the complete mitochondrial genome of T. liani, comprising 38,000 bp with a GC content of 24.61%. This genome includes 15 core protein-coding genes, 4 independent ORFs, 6 intronic ORFs, 26 tRNAs, and 2 rRNA genes. Phylogenetic analysis using Bayesian inference (BI) revealed the evolutionary relationships among 15 fungi from Eurotiales, strongly supporting distinct clades and indicating that T. liani most closely related to T. pinophilus.

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.

The complete mitochondrial genome of the basidiomycetous fungus, Tinctoporellus epimiltinus strain RS1

Tinctoporellus epimiltinus is widely known as a wood-decaying fungus. In the present study, we identified the complete mitochondrial genome of this species using next-generation sequencing technology. Our findings revealed that the genomic structure is a circular molecule with a size of 51,878 bp. Consistent with most Basidiomycota species, it consists of 14 core protein-coding genes, one ribosomal protein gene (rps3), 26 transfer RNA genes, and small and large ribosomal RNA (rns and rnl) genes. Seven additional open reading frames were identified. These included two sequences similar to DNA polymerases, an endonuclease-like sequence, and four hypothetical proteins. The mitochondrial genome exhibited a nucleotide composition of A (36.24%), C (12.04%), G (13.18%), and T (38.55%), resulting in a 25.21% GC content. A phylogenetic tree constructed using the combined mitochondrial gene dataset provided insight into the phylogenetic relationships of this species within the context of Basidiomycota and its members.

The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny

In the present study, three mitogenomes from the Bipolaris genus (Bipolaris maydis, B. zeicola, and B. oryzae) were assembled and compared with the other two reported Bipolaris mitogenomes (B. oryzae and B. sorokiniana). The five mitogenomes were all circular DNA molecules, with lengths ranging from 106,403 bp to 135,790 bp. The mitogenomes of the five Bipolaris species mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). The PCG length, AT skew and GC skew showed large variability among the 13 PCGs in the five mitogenomes. Across the 13 core PCGs tested, nad6 had the least genetic distance among the 16 Pleosporales species we investigated, indicating that this gene was highly conserved. In addition, the Ka/Ks values for all 12 core PCGs (excluding rps3) were < 1, suggesting that these genes were subject to purifying selection. Comparative mitogenomic analyses indicate that introns were the main factor contributing to the size variation of Bipolaris mitogenomes. The introns of the cox1 gene experienced frequent gain/loss events in Pleosporales species. The gene arrangement and collinearity in the mitogenomes of the five Bipolaris species were almost highly conserved within the genus. Phylogenetic analysis based on combined mitochondrial gene datasets showed that the five Bipolaris species formed well-supported topologies. This study is the first report on the mitogenomes of B. maydis and B. zeicola, as well as the first comparison of mitogenomes among Bipolaris species. The findings of this study will further advance investigations into the population genetics, evolution, and genomics of Bipolaris species.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Saccharomycopsis fibuligera (lindner) Klocker 1907 (saccharomycetales: saccharomycopsidaceae)

Saccharomycopsis fibuligera (Lindner) Klocker 1907 is frequently employed in the fermentation of metabolites such as citric acid, ethanol, mannitol, and pyruvate. Its heat tolerance and alcohol-producing capabilities during fermentation make it a desirable option for bread and wine production. To date, the mitochondrial genome of S. fibuligera has not been sequenced. In the present study, we obtained the full mitochondrial genome of S. fibuligera, which is 57,302 bp long and has a GC content of 24.40%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 21 intronic ORFs, 25 tRNAs, and 2 rRNA genes. By utilizing the Bayesian inference phylogenetic method, we constructed phylogenetic trees for 24 Saccharomycotina fungi, which indicated that S. fibuligera is closely related to S. capsularis.

Complete mitochondrial genome sequence and phylogenetic analysis of Tylopilus brunneirubens (Boletales, Basidiomycota)

Tylopilus brunneirubens is a common species in southern China. It is known for brown to dark brown pileus, white context turning reddish brown or rust brown when touched and distinct reticulation on the upper stem. However, little is known about its mitochondrial genome and its relationship with other boletes. Our analysis revealed that the mitochondrial genome of this species is a circular DNA molecule that spans 32,389 bp. It contains 15 core protein-coding genes, 24 transfer RNA genes, and two ribosomal RNA genes. The base composition of the mitochondrial genome is as follows: A (37.20%), C (11.32%), G (12.48%), and T (39.00%), with a GC content of 23.80%. Furthermore, a phylogenetic tree based on 24 mitochondrial genomes provided valuable insights into the phylogenetic relationships of Tylopilus brunneirubens with other boletes for the first time.

Mitochondrial genomic characteristics and phylogenetic analysis of a brewing fungus, Rhizopus microsporus Tiegh. 1875 (Mucorales: Rhizopodaceae)

Rhizopus microsporus Tiegh. 1875 is widely used in a variety of industries, such as brewing, wine making, baking, and medicine production, as it has the capability to break down proteins and generate surface-active agents. To date, the mitochondrial genome features of early evolved fungi from the Rhizopus genus have not been extensively studied. Our research obtained a full mitochondrial genome of R. microsporus species, which was 43,837 bp in size and had a GC content of 24.93%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 7 intronic ORFs, 24 tRNAs, and 2 rRNA genes. Through the use of the BI phylogenetic inference method, we were able to create phylogenetic trees for 25 early differentiation fungi which strongly supported the major clades; this indicated that R. microsporus is most closely related to Rhizopus oryzae.

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.

The second complete mitochondrial genome of Capillidium rhysosporum within the family Capillidiaceae, Entomophthorales

The complete mitochondrial genome of the entomophthoroid fungus Capillidium rhysosporum (strain no.: ATCC 12588) was sequenced using next-generation sequencing technology. The assembled circular genome has a length of 46,756 base pairs with a GC content of 27.06%. Gene prediction identified 15 core protein-coding genes (PCGs), two rRNA genes, and 27 tRNA genes. Phylogenetic analysis confirmed that C. rhysosporum belongs to the Zoopagomycota clade and is closely related to C. heterosporum. This study presents the second complete mitochondrial genome within the family Capillidiaceae, contributing to the mitochondrial DNA database of entomophthoroid fungi.

The complete mitochondrial genome of Penicillium expansum: Insights into the fungal evolution and phylogeny

Penicillium expansum is an important phytopathogenic fungus that causes blue mold disease. In this study, the novel mitochondrial genome of P. expansum was sequenced, assembled, annotated, and compared with the previously published Penicillium mitogenomes. P. expansum mitogenome is composed of circular DNA molecules with a genome size of 25,496 bp. It encodes 16 protein-encoding genes (PCGs), two rRNA genes, and 25 tRNA genes. Comparative analysis with six other Penicillium species revealed that gene length, GC content, AT skew, and GC skew were variable among the core protein-coding genes. The Penicillium species' gene synteny analysis identified several gene rearrangements. Among the core 15 PCGs, atp8 had the lowest K2P genetic distance, which shows that this gene is highly conserved. The Ka/Ks value of most PCGs was less than 1, which shows that these genes have undergone purifying selection. Phylogenetic analysis based on 14 concatenated core mitochondrial genes revealed that P. expansum shares a close relationship with P. solitum. This study served as a first report on the complete mitochondrial genome of P. expansum and its comparative analysis that will contribute to population genetics and rapid evolutionary studies among Penicillium species.

The complete mitochondrial genome of Bauhinia variegata (Leguminosae)

The mitogenome of Bauhinia variegate was assembled and characterized in this study. The mitogenome size was 437,271 bp, and its GC content was 45.5%. 36 protein-coding genes, 17 tRNAs and 3 rRNAs were annotated in the mitogenome. A total of 12 MTPTs, ranging from 71 bp to 3562 bp, were identified in the mitogenome and covered 1.46% (6373 bp) of the mitogenome. Phylogenetic analysis of 15 species of Leguminosae based on 23 core protein-coding genes showed that B. variegata was sister to Tylosema esculentum, another member from the subfamily Cercidoideae. The mitogenome of B. variegata provides a valuable genetic resource for further phylogenetic studies of this family.

Complete mitochondrial genome sequence of Aureoboletus raphanaceus (Boletales, Basidiomycota)

Aureoboletus raphanaceus is a member of boletoid mushroom, which is named after its distinctive radish smell. The mitochondrial genome and phylogenetic relationships with other boletes need to be investigated to gain a comprehensive understanding of it. In this study, we sequenced the mitochondrial genome of A. raphanaceus using next-generation sequencing technology and found that its mitochondrial genome is a circular DNA molecule measuring 42,157 bp. It consists of 15 core protein-coding genes, 27 transfer RNA genes, and two ribosomal RNA genes. The mitochondrial genome had a base composition of A (39.89%), C (11.06%), G (11.67%), and T (37.38%), with a GC content of 22.73%. A phylogenetic tree based on 22 mitochondrial genomes was constructed, which provided the first insights into the phylogenetic relationships of this species with related boletes.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes.

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown. In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains. The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species. This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

Analysis of synonymous codon usage patterns in mitochondrial genomes of nine Amanita species.

Codon basis is a common and complex natural phenomenon observed in many kinds of organisms. In the present study, we analyzed the base bias of 12 mitochondrial core protein-coding genes (PCGs) shared by nine Amanita species. The results showed that the codons of all Amanita species tended to end in A/T, demonstrating the preference of mitochondrial codons of Amanita species for a preference for this codon. In addition, we detected the correlation between codon base composition and the codon adaptation index (CAI), codon bias index (CBI), and frequency of optimal codons (FOP) indices, indicating the influence of base composition on codon bias. The average effective number of codons (ENC) of mitochondrial core PCGs of Amanita is 30.81, which is <35, demonstrating the strong codon preference of mitochondrial core PCGs of Amanita. The neutrality plot analysis and PR2-Bias plot analysis further demonstrated that natural selection plays an important role in Amanita codon bias. In addition, we obtained 5-10 optimal codons (ΔRSCU > 0.08 and RSCU > 1) in nine Amanita species, and GCA and AUU were the most widely used optimal codons. Based on the combined mitochondrial sequence and RSCU value, we deduced the genetic relationship between different Amanita species and found large variations between them. This study promoted the understanding of synonymous codon usage characteristics and evolution 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.

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.