Deletion of the sex-determining gene SXI1α enhances the spread of mitochondrial introns in Cryptococcus neoformans

Distribution and Polymorphisms of Group I Introns in Mitochondrial Genes from Cryptococcus neoformans and Cryptococcus gattii

The species complexes Cryptococcus neoformans and Cryptococcus gattii are the causative agents of cryptococcosis. Virulence and susceptibility to antifungals may vary within each species according to the fungal genotype. Therefore, specific and easily accessible molecular markers are required to distinguish cryptic species and/or genotypes. Group I introns are potential markers for this purpose because they are polymorphic concerning their presence and sequence. Therefore, in this study, we evaluated the presence of group I introns in the mitochondrial genes cob and cox1 in different Cryptococcus isolates. Additionally, the origin, distribution, and evolution of these introns were investigated by phylogenetic analyses, including previously sequenced introns for the mtLSU gene. Approximately 80.5% of the 36 sequenced introns presented homing endonucleases, and phylogenetic analyses revealed that introns occupying the same insertion site form monophyletic clades. This suggests that they likely share a common ancestor that invaded the site prior to species divergence. There was only one case of heterologous invasion, probably through horizontal transfer to C. decagattii (VGIV genotype) from another fungal species. Our results showed that the C. neoformans complex has fewer introns compared to C. gattii. Additionally, there is significant polymorphism in the presence and size of these elements, both among and within genotypes. As a result, it is impossible to differentiate the cryptic species using a single intron. However, it was possible to differentiate among genotypes within each species complex, by combining PCRs of mtLSU and cox1 introns, for C. neoformans species, and mtLSU and cob introns for C. gattii species.

Insight into the evolutionary and domesticated history of the most widely cultivated mushroom Agaricus bisporus via mitogenome sequences of 361 global strains

Agaricus bisporus is the most widely cultivated edible mushroom in the world with a only around three hundred years known history of cultivation. Therefore, it represents an ideal organism not only to investigate the natural evolutionary history but also the understanding on the evolution going back to the early era of domestication. In this study, we generated the mitochondrial genome sequences of 352 A. bisporus strains and 9 strains from 4 closely related species around the world. The population mitogenomic study revealed all A. bisporus strains can be divided into seven clades, and all domesticated cultivars present only in two of those clades. The molecular dating analysis showed this species origin in Europe on 4.6 Ma and we proposed the main dispersal routes. The detailed mitogenome structure studies showed that the insertion of the plasmid-derived dpo gene caused a long fragment (MIR) inversion, and the distributions of the fragments of dpo gene were strictly in correspondence with these seven clades. Our studies also showed A. bisporus population contains 30 intron distribution patterns (IDPs), while all cultivars contain only two IDPs, which clearly exhibit intron loss compared to the others. Either the loss occurred before or after domestication, that could suggest that the change facilitates their adaptation to the cultivated environment.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

Regulatory basis for reproductive flexibility in a meningitis-causing fungal pathogen

Pathogenic fungi of the genus Cryptococcus can undergo two sexual cycles, involving either bisexual diploidization (after fusion of haploid cells of different mating type) or unisexual diploidization (by autodiploidization of a single cell). Here, we construct a gene-deletion library for 111 transcription factor genes in Cryptococcus deneoformans, and explore the roles of these regulatory networks in the two reproductive modes. We show that transcription factors crucial for bisexual syngamy induce the expression of known mating determinants as well as other conserved genes of unknown function. Deletion of one of these genes, which we term FMP1, leads to defects in bisexual reproduction in C. deneoformans, its sister species Cryptococcus neoformans, and the ascomycete Neurospora crassa. Furthermore, we show that a recently evolved regulatory cascade mediates pre-meiotic unisexual autodiploidization, supporting that this reproductive process is a recent evolutionary innovation. Our findings indicate that genetic circuits with different evolutionary ages govern hallmark events distinguishing unisexual and bisexual reproduction in Cryptococcus.

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.

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.

The mitochondrial genome of Ophiostoma himal-ulmi and comparison with other fungi causing Dutch elm disease

The mitochondrial genome of Ophiostoma himal-ulmi, a species endemic to the Western Himalayas and one of the fungi that cause Dutch elm disease, has been sequenced and characterized. The mitochondrial genome was compared with other available genomes for members of the Ophiostomatales, including other agents of Dutch elm disease (Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi, and Ophiostoma novo-ulmi subspecies americana), and it was observed that gene synteny is highly conserved, and variability among members of the fungi that cause Dutch-elm disease is primarily due to the number of intron insertions. Among the fungi that cause Dutch elm disease that we examined, O. himal-ulmi has the largest mitochondrial genomes (ranging from 94 934 to 111 712 bp), owing to the expansion of the number of introns.

Mitochondrial Inheritance in Phytopathogenic Fungi-Everything Is Known, or Is It?

Mitochondria are important organelles in eukaryotes that provide energy for cellular processes. Their function is highly conserved and depends on the expression of nuclear encoded genes and genes encoded in the organellar genome. Mitochondrial DNA replication is independent of the replication control of nuclear DNA and as such, mitochondria may behave as selfish elements, so they need to be controlled, maintained and reliably inherited to progeny. Phytopathogenic fungi meet with special environmental challenges within the plant host that might depend on and influence mitochondrial functions and services. We find that this topic is basically unexplored in the literature, so this review largely depends on work published in other systems. In trying to answer elemental questions on mitochondrial functioning, we aim to introduce the aspect of mitochondrial functions and services to the study of plant-microbe-interactions and stimulate phytopathologists to consider research on this important organelle in their future projects.

Mitochondrial Genome Polymorphisms in the Human Pathogenic Fungus Cryptococcus neoformans

The Cryptococcus complex consists of at least seven evolutionary divergent lineages and causes ∼200,000 fatal human infections each year worldwide. The dominant lineage is Cryptococcus neoformans which consists of three haploid clades VNI, VNII, and VNB, their haploid hybrids, and various diploids derived from intra- and inter-clade mating events. In this study, we analyzed the mitogenomes of 184 strains of C. neoformans. Our analyses revealed that all 184 mitogenomes contained the same 15 protein-coding genes in the same gene order. However, their mitogenome sizes varied between 24,740 and 31,327 bp, primarily due to differences in the number and size of mitochondrial introns. Twelve nucleotide sites within five mitochondrial genes were found to contain introns in at least one of the 184 strains, ranging from 2 to 7 introns within each mitogenome. The concatenated mitochondrial exon sequences of the 15 protein-coding genes and two rRNA genes showed that VNI, VNII, and VNB strains were separated into distinct clades or sub-clades, largely consistent with results based on nuclear genome SNPs. However, several novel findings were observed. First, one strain of the VNB clade contained mitogenome exon sequences identical to the main VNI mitogenome type but was distant to other VNB mitogenomes. Second, hybrids among clades VNI, VNII, and VNB identified based on their nuclear genome SNPs contained mitogenomes from different clades, with evidence of their mitogenomes inherited from either the MAT a or the MAT α parents. Third, the eight diploid VNB (C. neoformans) × VNIV (C. deneoformans) hybrids contained recombinant mitogenomes. Fourth, analyses of intron distribution and the paired exon-intron phylogenies for each of the 12 exon-intron pairs suggested frequent gains and losses of mitochondrial introns during the evolution of C. neoformans. The combined mitogenome exon-based phylogeny and intron distributions suggested that clades VNI, VNII and VNB could be further divided into sub-clades. Together, our results revealed a dynamic evolution of mitochondrial genomes in this important human fungal pathogen.

Comparative Analyses of Mitochondrial Genomes Provide Evolutionary Insights Into Nematode-Trapping Fungi

Predatory fungi in Orbiliaceae (Ascomycota) have evolved a diversity of trapping devices that enable them to trap and kill nematodes, other small animals, and protozoans. These trapping devices include adhesive hyphae, adhesive knobs, adhesive networks, constricting rings, and non-constricting rings. Their diversity and practical importance have attracted significant attention from biologists, making them excellent model organisms for studying adaptative evolution and as biological control agents against parasitic nematodes. The putative origins and evolutionary relationships among these carnivorous fungi have been investigated using nuclear protein-encoding genes, but their patterns of mitogenome relationships and divergences remain unknown. Here we analyze and compare the mitogenomes of 12 fungal strains belonging to eight species, including six species representing all four types of nematode trapping devices and two from related but non-predatory fungi. All 12 analyzed mitogenomes were of circular DNA molecules, with lengths ranging from 146,101 bp to 280,699 bp. Gene synteny analysis revealed several gene rearrangements and intron transfers among the mitogenomes. In addition, the number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among these mitogenomes. The increased number and total size of introns were the main contributors to the length differences among the mitogenomes. Phylogenetic analyses of the protein-coding genes indicated that mitochondrial and nuclear genomes evolved at different rates, and signals of positive selection were found in several genes involved in energy metabolism. Our study provides novel insights into the evolution of nematode-trapping fungi and shall facilitate further investigations of this ecologically and agriculturally important group of fungi.

Evidence of Extensive Intraspecific Noncoding Reshuffling in a 169-kb Mitochondrial Genome of a Basidiomycetous Fungus

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes’ intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

The Pheromone and Pheromone Receptor Mating-Type Locus Is Involved in Controlling Uniparental Mitochondrial Inheritance in Cryptococcus

Mitochondria are inherited uniparentally during sexual reproduction in the majority of eukaryotic species studied, including humans, mice, and nematodes, as well as many fungal species. Mitochondrial uniparental inheritance (mito-UPI) could be beneficial in that it avoids possible genetic conflicts between organelles with different genetic backgrounds, as recently shown in mice, and it could prevent the spread of selfish genetic elements in the mitochondrial genome. Despite the prevalence of observed mito-UPI, the underlying mechanisms and the genes involved in controlling this non-Mendelian inheritance are poorly understood in many species. In Cryptococcus neoformans, a human pathogenic basidiomyceteous fungus, mating types (MATα and MATa) are defined by alternate alleles at the single MAT locus that evolved from fusion of the two MAT loci (P/R encoding pheromones and pheromone receptors, and HD encoding homeodomain transcription factors) that are the ancestral state in the basidiomycota. Mitochondria are inherited uniparentally from the MATa parent in C. neoformans, and this requires the SXI1α and SXI2a HD factors encoded by MAT However, there is evidence that additional genes contribute to the control of mito-UPI in Cryptococcus Here, we show that in C. amylolentus, a sibling species of C. neoformans with unlinked P/R and HD MAT loci, mito-UPI is controlled by the P/R locus and is independent of the HD locus. Consistently, by replacing the MATα alleles of the pheromones (MF) and pheromone receptor (STE3) with the MATa alleles, we show that these P/R locus-defining genes indeed affect mito-UPI in C. neoformans during sexual reproduction. Additionally, we show that during early stages of C. neoformans sexual reproduction, conjugation tubes are always produced by the MATα cells, resulting in unidirectional migration of the MATα nucleus into the MATa cell during zygote formation. This process is controlled by the P/R locus and could serve to physically restrict movement of MATα mitochondria in the zygotes, and thereby contribute to mito-UPI. We propose a model in which both physical and genetic mechanisms function in concert to prevent the coexistence of mitochondria from the two parents in the zygote, and subsequently in the meiotic progeny, thus ensuring mito-UPI in pathogenic Cryptococcus, as well as in closely related nonpathogenic species. The implications of these findings are discussed in the context of the evolution of mito-UPI in fungi and other more diverse eukaryotes.

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha-proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein-coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome-wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad-scale genome streamlining via loss of protein-coding genes as well as noncoding, parasitic DNA elements.

Fungal mitochondrial genomes and genetic polymorphisms

Mitochondria are the powerhouses of eukaryotic cells, responsible for ATP generation and playing a role in a diversity of cellular and organismal functions. Different from the majority of other intracellular membrane structures, mitochondria contain their own genetic materials that are capable of independent replication and inheritance. In this mini-review, we provide brief summaries of fungal mitochondrial genome structure, size, gene content, inheritance, and genetic variation. We pay special attention to the relative genetic polymorphisms of the mitochondrial vs nuclear genomes at the population level within individual fungal species. Among the 20 species/groups of species reviewed here, there is a range of variation among genes and species in the relative nuclear and mitochondrial genetic polymorphisms. Interestingly, most (15/20) showed a greater genetic diversity for nuclear genes and genomes than for mitochondrial genes and genomes, with the remaining five showing similar or slower nuclear genome genetic variations. This fungal pattern is different from the dominant pattern in animals, but more similar to that in plants. At present, the mechanisms for the variations among fungal species and the overall low level of mitochondrial sequence polymorphisms are not known. The increasing availability of population genomic data should help us reveal the potential genetic and ecological factors responsible for the observed variations.