Analysis of fungal genomes reveals commonalities of intron gain or loss and functions in intron-poor species

Retracing the evolution of Pneumocystis species, with a focus on the human pathogen Pneumocystis jirovecii

SUMMARYEvery human being is presumed to be infected by the fungus Pneumocystis jirovecii at least once in his or her lifetime. This fungus belongs to a large group of species that appear to exclusively infect mammals, with P. jirovecii being the only one known to cause disease in humans. The mystery of P. jirovecii origin and speciation is just beginning to unravel. Here, we provide a review of the major steps of P. jirovecii evolution. The Pneumocystis genus likely originated from soil or plant-associated organisms during the period of Cretaceous ~165 million years ago and successfully shifted to mammals. The transition coincided with a substantial loss of genes, many of which are related to the synthesis of nutrients that can be scavenged from hosts or cell wall components that could be targeted by the mammalian immune system. Following the transition, the Pneumocystis genus cospeciated with mammals. Each species specialized at infecting its own host. Host specialization is presumably built at least partially upon surface glycoproteins, whose protogene was acquired prior to the genus formation. P. jirovecii appeared at ~65 million years ago, overlapping with the emergence of the first primates. P. jirovecii and its sister species P. macacae, which infects macaques nowadays, may have had overlapping host ranges in the distant past. Clues from molecular clocks suggest that P. jirovecii did not cospeciate with humans. Molecular evidence suggests that Pneumocystis speciation involved chromosomal rearrangements and the mounting of genetic barriers that inhibit gene flow among species.

How to overcome constraints imposed by microsporidian genome features to ensure gene prediction?

Since the advent of sequencing techniques and due to their continuous evolution, it has become easier and less expensive to obtain the complete genome sequence of any organism. Nevertheless, to elucidate all biological processes governing organism development, quality annotation is essential. In genome annotation, predicting gene structure is one of the most important and captivating challenges for computational biology. This aspect of annotation requires continual optimization, particularly for genomes as unusual as those of microsporidia. Indeed, this group of fungal-related parasites exhibits specific features (highly reduced gene sizes, sequences with high rate of evolution) linked to their evolution as intracellular parasites, requiring the implementation of specific annotation approaches to consider all these features. This review aimed to outline these characteristics and to assess the increasingly efficient approaches and tools that have enhanced the accuracy of gene prediction for microsporidia, both in terms of sensitivity and specificity. Subsequently, a final part will be dedicated to postgenomic approaches aimed at reinforcing the annotation data generated by prediction software. These approaches include the characterization of other understudied genes, such as those encoding regulatory noncoding RNAs or very small proteins, which also play crucial roles in the life cycle of these microorganisms.

Improved recovery and annotation of genes in metagenomes through the prediction of fungal introns

Metagenomics provides a tool to assess the functional potential of environmental and host-associated microbiomes based on the analysis of environmental DNA: assembly, gene prediction and annotation. While gene prediction is straightforward for most bacterial and archaeal taxa, it has limited applicability in the majority of eukaryotic organisms, including fungi that contain introns in gene coding sequences. As a consequence, eukaryotic genes are underrepresented in metagenomics datasets and our understanding of the contribution of fungi and other eukaryotes to microbiome functioning is limited. Here, we developed a machine intelligence-based algorithm that predicts fungal introns in environmental DNA with reasonable precision and used it to improve the annotation of environmental metagenomes. Intron removal increased the number of predicted genes by up to 9.1% and improved the annotation of several others. The proportion of newly predicted genes increased with the share of eukaryotic genes in the metagenome and-within fungal taxa-increased with the number of introns per gene. Our approach provides a tool named SVMmycointron for improved metagenome annotation, especially of microbiomes with a high proportion of eukaryotes. The scripts described in the paper are made publicly available and can be readily utilized by microbiome researchers analysing metagenomics data.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians

BACKGROUND

Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution.

RESULTS

At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths.

CONCLUSIONS

Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

Functional genome annotation and transcriptome analysis of Pseudozyma hubeiensis BOT-O, an oleaginous yeast that utilizes glucose and xylose at equal rates

Pseudozyma hubeiensis is a basidiomycete yeast that has the highly desirable traits for lignocellulose valorisation of being equally efficient at utilization of glucose and xylose, and capable of their co-utilization. The species has previously mainly been studied for its capacity to produce secreted biosurfactants in the form of mannosylerythritol lipids, but it is also an oleaginous species capable of accumulating high levels of triacylglycerol storage lipids during nutrient starvation. In this study, we aimed to further characterize the oleaginous nature of P. hubeiensis by evaluating metabolism and gene expression responses during storage lipid formation conditions with glucose or xylose as a carbon source. The genome of the recently isolated P. hubeiensis BOT-O strain was sequenced using MinION long-read sequencing and resulted in the most contiguous P. hubeiensis assembly to date with 18.95 Mb in 31 contigs. Using transcriptome data as experimental support, we generated the first mRNA-supported P. hubeiensis genome annotation and identified 6540 genes. 80% of the predicted genes were assigned functional annotations based on protein homology to other yeasts. Based on the annotation, key metabolic pathways in BOT-O were reconstructed, including pathways for storage lipids, mannosylerythritol lipids and xylose assimilation. BOT-O was confirmed to consume glucose and xylose at equal rates, but during mixed glucose-xylose cultivation glucose was found to be taken up faster. Differential expression analysis revealed that only a total of 122 genes were significantly differentially expressed at a cut-off of |log₂ fold change| ≥ 2 when comparing cultivation on xylose with glucose, during exponential growth and during nitrogen-starvation. Of these 122 genes, a core-set of 24 genes was identified that were differentially expressed at all time points. Nitrogen-starvation resulted in a larger transcriptional effect, with a total of 1179 genes with significant expression changes at the designated fold change cut-off compared with exponential growth on either glucose or xylose.

Psychogenetic, genetic and epigenetic mechanisms in Candida auris: Role in drug resistance

In recent years, we are facing the challenge of drug resistance emergence in fungi. The availability of limited antifungals and development of multi-drug resistance in fungal pathogens has become a serious concern in the past years in the health sector. Although several cellular, molecular, and genetic mechanisms have been proposed to explain the drug resistance mechanism in fungi, but a complete understanding of the molecular and genetic mechanisms is still lacking. Besides the genetic mechanism, epigenetic mechanisms are pivotal in the fungal lifecycle and disease biology. However, very little is understood about the role of epigenetic mechanisms in the emergence of multi-drug resistance in fungi, especially in Candida auris (C. auris). The current narrative review summaries the clinical characteristics, genomic organization, and molecular/genetic/epigenetic mechanisms underlying the emergence of drug resistance in C. auris. A very few studies have attempted to evaluate the role of epigenetic mechanisms in C. auris. Furthermore, advanced genetic tools such as the CRISP-Cas9 system can be utilized to elucidate the epigenetic mechanisms and their role in the emergence of multi-drug resistance in C. auris.

Gene Structure Evolution of the Short-Chain Dehydrogenase/Reductase (SDR) Family

SDR (Short-chain Dehydrogenases/Reductases) are one of the oldest and heterogeneous superfamily of proteins, whose classification is problematic because of the low percent identity, even within families. To get clearer insights into SDR molecular evolution, we explored the splicing site organization of the 75 human SDR genes across their vertebrate and invertebrate orthologs. We found anomalous gene structures in members of the human SDR7C and SDR42E families that provide clues of retrogene properties and independent evolutionary trajectories from a common invertebrate ancestor. The same analyses revealed that the identity value between human and invertebrate non-allelic variants is not necessarily associated with the homologous gene structure. Accordingly, a revision of the SDR nomenclature is proposed by including the human SDR40C1 and SDR7C gene in the same family.

The International Space Station Environment Triggers Molecular Responses in Aspergillus niger

Due to immense phenotypic plasticity and adaptability, Aspergillus niger is a cosmopolitan fungus that thrives in versatile environments, including the International Space Station (ISS). This is the first report of genomic, proteomic, and metabolomic alterations observed in A. niger strain JSC-093350089 grown in a controlled experiment aboard the ISS. Whole-genome sequencing (WGS) revealed that ISS conditions, including microgravity and enhanced irradiation, triggered non-synonymous point mutations in specific regions, chromosomes VIII and XII of the JSC-093350089 genome when compared to the ground-grown control. Proteome analysis showed altered abundance of proteins involved in carbohydrate metabolism, stress response, and cellular amino acid and protein catabolic processes following growth aboard the ISS. Metabolome analysis further confirmed that space conditions altered molecular suite of ISS-grown A. niger JSC-093350089. After regrowing both strains on Earth, production of antioxidant—Pyranonigrin A was significantly induced in the ISS-flown, but not the ground control strain. In summary, the microgravity and enhanced irradiation triggered unique molecular responses in the A. niger JSC-093350089 suggesting adaptive responses.

Intron losses and gains in the nematodes

BACKGROUND

The evolution of spliceosomal introns has been widely studied among various eukaryotic groups. Researchers nearly reached the consensuses on the pattern and the mechanisms of intron losses and gains across eukaryotes. However, according to previous studies that analyzed a few genes or genomes, Nematoda seems to be an eccentric group.

RESULTS

Taking advantage of the recent accumulation of sequenced genomes, we extensively analyzed the intron losses and gains using 104 nematode genomes across all the five Clades of the phylum. Nematodes have a wide range of intron density, from less than one to more than nine per kbp coding sequence. The rates of intron losses and gains exhibit significant heterogeneity both across different nematode lineages and across different evolutionary stages of the same lineage. The frequency of intron losses far exceeds that of intron gains. Five pieces of evidence supporting the model of cDNA-mediated intron loss have been observed in ten Caenorhabditis species, the dominance of the precise intron losses, frequent loss of adjacent introns, high-level expression of the intron-lost genes, preferential losses of short introns, and the preferential losses of introns close to 3'-ends of genes. Like studies in most eukaryotic groups, we cannot find the source sequences for the limited number of intron gains detected in the Caenorhabditis genomes.

CONCLUSIONS

These results indicate that nematodes are a typical eukaryotic group rather than an outlier in intron evolution.

Evolutionary Morphogenesis of Sexual Fruiting Bodies in Basidiomycota: Toward a New Evo-Devo Synthesis

The development of sexual fruiting bodies is one of the most complex morphogenetic processes in fungi. Mycologists have long been fascinated by the morphological and developmental diversity of fruiting bodies; however, evolutionary developmental biology of fungi still lags significantly behind that of animals or plants. Here, we summarize the current state of knowledge on fruiting bodies of mushroom-forming Basidiomycota, focusing on phylogenetic and developmental biology. Phylogenetic approaches have revealed a complex history of morphological transformations and convergence in fruiting body morphologies. Frequent transformations and convergence is characteristic of fruiting bodies in contrast to animals or plants, where main body plans are highly conserved. At the same time, insights into the genetic bases of fruiting body development have been achieved using forward and reverse genetic approaches in selected model systems. Phylogenetic and developmental studies of fruiting bodies have each yielded major advances, but they have produced largely disjunct bodies of knowledge. An integrative approach, combining phylogenetic, developmental, and functional biology, is needed to achieve a true fungal evolutionary developmental biology (evo-devo) synthesis for fungal fruiting bodies.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Modeling in yeast how rDNA introns slow growth and increase desiccation tolerance in lichens

We connect ribosome biogenesis to desiccation tolerance in lichens, widespread symbioses between specialized fungi (mycobionts) and unicellular phototrophs. We test whether the introns present in the nuclear ribosomal DNA of lichen mycobionts contribute to their anhydrobiosis. Self-splicing introns are found in the rDNA of several eukaryotic microorganisms, but most introns populating lichen rDNA are unable to self-splice, being either catalytically impaired group I introns, or spliceosomal introns ectopically present in rDNA. Although the mycobiont’s splicing machinery removes all introns from rRNA, Northern analysis indicates delayed post-transcriptional removal during rRNA processing, suggesting interference with ribosome assembly. To study the effects of lichen introns in a model system, we used CRISPR to introduce a spliceosomal rDNA intron from the lichen fungus Cladonia grayi into all nuclear rDNA copies of Saccharomyces cerevisiae, which lacks rDNA introns. Three intron-bearing yeast mutants were constructed with the intron inserted either in the 18S rRNA genes, the 25S rRNA genes, or in both. The mutants removed the introns correctly but had half the rDNA genes of the wildtype, grew 4.4–6 times slower, and were 40–1700 times more desiccation tolerant depending on intron position and number. Intracellular trehalose, a disaccharide implicated in desiccation tolerance, was detected at low concentration. Our data suggest that the interference of the splicing machinery with ribosome assembly leads to fewer ribosomes and proteins and to slow growth and increased desiccation tolerance in the yeast mutants. The relevance of these findings for slow growth and desiccation tolerance in lichens is discussed.