Comparative genomics of the extremophile Cryomyces antarcticus and other psychrophilic Dothideomycetes

Over a billion years of fungal evolution has enabled representatives of this kingdom to populate almost all parts of planet Earth and to adapt to some of its most uninhabitable environments including extremes of temperature, salinity, pH, water, light, or other sources of radiation. Cryomyces antarcticus is an endolithic fungus that inhabits rock outcrops in Antarctica. It survives extremes of cold, humidity and solar radiation in one of the least habitable environments on Earth. This fungus is unusual because it produces heavily melanized, meristematic growth and is thought to be haploid and asexual. Due to its growth in the most extreme environment, it has been suggested as an organism that could survive on Mars. However, the mechanisms it uses to achieve its extremophilic nature are not known. Comparative genomics can provide clues to the processes underlying biological diversity, evolution, and adaptation. This effort has been greatly facilitated by the 1000 Fungal Genomes project and the JGI MycoCosm portal where sequenced genomes have been assembled into phylogenetic and ecological groups representing different projects, lifestyles, ecologies, and evolutionary histories. Comparative genomics within and between these groups provides insights into fungal adaptations, for example to extreme environmental conditions. Here, we analyze two Cryomyces genomes in the context of additional psychrophilic fungi, as well as non-psychrophilic fungi with diverse lifestyles selected from the MycoCosm database. This analysis identifies families of genes that are expanded and contracted in Cryomyces and other psychrophiles and may explain their extremophilic lifestyle. Higher GC contents of genes and of bases in the third positions of codons may help to stabilize DNA under extreme conditions. Numerous smaller contigs in C. antarcticus suggest the presence of an alternative haplotype that could indicate the sequenced isolate is diploid or dikaryotic. These analyses provide a first step to unraveling the secrets of the extreme lifestyle of C. antarcticus.

1 Introduction

Life in extreme cold and high-altitude environments like Antarctica and the Alps may be severely restricted, but it is not entirely absent. Under these conditions, rocks are the principal habitat for microorganisms, which are known as epilithic when growing on rock surfaces or endolithic when growing inside the substrate. Endolithic organisms inhabiting rock are categorized into chasmoendoliths, which inhabit fissures and cracks, and cryptoendoliths, which inhabit porous rock substrates (Friedmann, 1982). Cryptoendoliths comprise filamentous fungi, unicellular green algae, and bacteria. Fungi growing under these conditions are known as microcolonial fungi (MCF) or rock-inhabiting fungi (RIF). Most of these organisms belong to the order Dothideomycetes but some are placed in the Eurotiomycetes or remain taxonomically undetermined (Ruibal et al., 2009; Selbmann et al., 2015).

Cryptoendolithic species in genera such as Cryomyces and Friedmanniomyces have thick, melanized cell walls, produce exopolysaccharides (extracellular polymeric substances, or EPS) outside the hyphae, and have other adaptations to harsh environments (Selbmann et al., 2005). Physiological adaptations such as the production of trehalose and cryoprotectant sugars, polyols, lipids/fatty acids, antifreeze proteins (Robinson, 2001), and hydrolases and oxidoreductases (Duarte et al., 2018; Coleine et al., 2020) play an important role in the adaptation of these organisms to cold. In general, psychrophilic black fungi show increased protein expression when grown at 1 °C; Friedmanniomyces endolithicus, one of the species included in this work, showed increased numbers of high-molecular-weight proteins in comparison to a mesophilic reference species (Tesei et al., 2012), suggesting substantial reprogramming of protein expression at low versus high temperatures. These and other physiological adaptations to extreme environments should be reflected in genome modifications.

The production of proteins and enzymes underlies the physiological mechanism of cold tolerance in fungi. Observed changes occur in both gene regulation and as mutational changes to individual proteins. One of the simplest adaptations is the expansion or contraction of paralogous groups. Changes in gene number provide a rapid means to modify proteome composition, enzymatic activity, and the constituents of the membrane, all of which are important in cold/stress adaptation. However, these mechanisms are complex and not fully understood (Maggi et al., 2013).

Adaptation to cold environments may require the alteration of proteins such as those related to membrane transport and lipid biosynthesis to maintain membrane flexibility at low temperatures. A proteome comparison can identify proteins that are significantly altered in cold-adapted fungi. However, studies of proteome composition mostly have been limited to the comparison of only a few species and conditions at a time (Tesei et al., 2012), limiting our ability to see general proteomic patterns related to adaptation and survival in extreme conditions. Identification of such protein and enzyme patterns in the proteomes of cold-adapted fungi can identify proteins that are responsible for physiological mechanisms of cold tolerance.

Expansion or contraction of orthologous protein families is one of the simplest ways for organisms to adapt to an extreme environment. Expansion of specific groups has been identified in many species (Su et al., 2016; Stajich, 2017; Paiva et al., 2024). Previous analysis of a draft genome of C. antarcticus, sequenced with Ion PI, identified little that stood out, causing the authors to conclude there was “nothing special in the specialist” (Sterflinger et al., 2014). However, that analysis included comparisons with only five other species, only three of which were in the Dothideomycetes, and the assembly was very fragmented. In this work, we compare the proteome composition of 52 species in the class Dothideomycetes, 19 of which are extremophiles (cold, acidophilic, ethanol vapor, xerophilic, or salt tolerant). The selected species include 11 psychrophiles which are the primary targets of our analysis. The non-psychrophilic species include lichens, yeasts, parasitic species, opportunistic animal pathogens, and plant pathogens from selected Dothideomycetes groups. We examine and compare the expansion and contraction of orthologous groups to test the hypothesis that families of genes associated with adaptation may be altered in these stress-tolerant Dothideomycetes.

Based on previously published work (reviewed by Gostinčar et al., 2023), we expect to see changes in genes related to pigment production, cell wall biosynthesis, membrane transporters, production of compatible solutes (such as glycerol), fatty acid synthesis, and response to reactive oxygen species (ROS). In addition, specific groups necessary for survival in other environments, for instance genes involved in lignin degradation and carbohydrate‐active enzymes (CAZymes) required for pathogenicity to plants, may be dispensable in environments where the corresponding substrates are unavailable. By synthesizing results across a broad range of species, we hope to identify general principles of adaptation to extremely cold environments and determine whether these mechanisms are also found in fungal species adapted to other types of stresses.

2 Materials and methods

2.1 Genome sequencing, assembly, and annotation

Cryomyces antarcticus isolate 116301 was obtained from the CBS culture collection (Utrecht, the Netherlands). The fungus was cultured on 2% potato dextrose agar plates at 4 C for more than 6 months. Mycelia were scraped off the surface of the plates and lyophilized. DNA was extracted using Qiagen’s DNeasy kit and RNA was extracted using the RNeasy kit. Its genome was sequenced using the PacBio platform. For the PacBio library, 5 μg of genomic DNA was sheared to >10 kb using Covaris g-Tubes. The sheared DNA was treated with exonuclease to remove single-stranded ends followed by end repair and ligation of blunt adapters using the SMRTbell Template Prep Kit 1.0 (Pacific Biosciences). The library was purified with AMPure PB beads. PacBio Sequencing primer was then annealed to the SMRTbell template library and sequencing polymerase was bound to them using the Sequel II Binding kit 1.0. The prepared SMRTbell template libraries were then sequenced on a Pacific Biosystems Sequel II sequencer using 8M v1 SMRT cells and Version 1.0 sequencing chemistry with 1×900 min sequencing movie run times. Filtered subread data were processed with the JGI QC pipeline to remove artifacts and organellar reads. Filtered Circular Consensus Sequence (CCS) reads were assembled with Flye version 2.7.1-b1590 (Kolmogorov et al., 2019) to generate an assembly and polished with the PacBio SMRTLINK (v8.0.0.80529) using the gcpp –algorithm arrow command.

For the Cryomyces antarcticus transcriptome, plate-based RNA sample prep was performed using the PerkinElmer Sciclone NGS robotic liquid handling system and Illumina’s TruSeq Stranded mRNA HT sample prep kit, using poly-A selection of mRNA following the protocol outlined by Illumina in their user guide (https://support.illumina.com/sequencing/sequencing_kits/truseq-stranded-mrna.html) and with the following conditions: total RNA starting material was 1 μg per sample and 8 cycles of PCR were used for library amplification. The prepared library was then quantified using KAPA Biosystem’s next-generation sequencing library qPCR kit (Roche) and run on a Roche LightCycler 480 real-time PCR instrument. The quantified library was then multiplexed with other libraries, and the pool of libraries was prepared for sequencing on the Illumina NovaSeq 6000 sequencing platform using NovaSeq XP v1 reagent kits, S4 flow cell, following a 2×150 indexed run recipe. Raw RNAseq reads were filtered and trimmed using the BBDuk kmer matching (kmer=25) and phred trimming method set at Q6. Following trimming, reads under the length threshold were removed (minimum length of 25 bases or 1/3 of the original read length - whichever is longer). Filtered reads were used as input for de novo assembly of RNA contigs using Trinity version 2.8.5 with parameters –normalize_reads and –jaccard_clip.

The genome was annotated using the JGI Annotation pipeline and made available at the JGI MycoCosm genome portal (http://mycocosm.jgi.doe.gov/Cryan3) along with tools for interactive comparative analysis (Grigoriev et al., 2014). Both the assembly and annotation have been deposited to GenBank under accession JBDRTV000000000.

2.2 Acquisition of additional genome data

The proteomes of 52 fungal species were obtained from the JGI MycoCosm portal (https://mycocosm.jgi.doe.gov/mycocosm/; Grigoriev et al., 2014) (Table 1). The majority of the species belong to the subclass Dothideomycetidae within the Dothideomycetes. Members of this subclass were chosen for their close phylogenetic relationship to Cryomyces antarcticus and C. minteri, as indicated by the phylogenetic tree available in MycoCosm (https://mycocosm.jgi.doe.gov/mycocosm/species-tree/tree;j7jS2o?organism=dothideomycetes). The group of psychrophilic fungi (https://mycocosm.jgi.doe.gov/Psychrophilic_fungi) consists of 11 species, including our target species C. antarcticus and C. minteri, as well as the species Aureobasidium subglaciare, Extremus antarcticus, Friedmanniomyces simplex, F. endolithicus, Hortaea thailandica, Rachicladosporium sp., Rachicladosporium antarcticum, and two isolates of the genus Alternaria from the order Pleosporales, Alternaria sp. UNIPAMPA012 and Alternaria sp. UNIPAMPA017. Both these species were isolated from healthy leaves of Antarctic hair grass (https://mycocosm.jgi.doe.gov/Altsp012_1/Altsp012_1.home.html).

Table 1

Table 1. Summary information about the Dothideomycetes fungal species included in the comparative genomic analyses.

Among the selected species, we include representatives of different lifestyles: endophytic, endolithic, molds, yeasts, opportunistic human pathogens, parasites, and plant pathogens (Table 1). The non-psychrophilic group contains fungal pathogens of a diverse range of hosts including the sugar beet pathogen Cercospora berteroae, the corn pathogen Cercospora zeae-maydis, the tomato and grape pathogens Cladosporium fulvum and Elsinoë ampelina, the banana pathogens Pseudocercospora (formerly Mycosphaerella) eumusae, P. (Mycosphaerella) fijiensis, and P. musae, the wheat pathogen Zymoseptoria tritici (formerly Mycosphaerella graminicola) and its close relatives Z. ardabiliae, Z. brevis, and Z. pseudotritici (isolated from wild grass species). We also included the forest pathogens Dothistroma septosporum, Pseudocercospora ulei, Sphaerulina musiva (formerly Septoria musiva; teleomorph: Mycosphaerella populorum), Sphaerulina populicola (formerly Septoria populicola), and Teratosphaeria nubilosa. In addition, other extremophiles are included in the analysis to compare the survival strategies under different extreme environmental conditions: the acidophilic fungi Acidomyces richmondensis and Hortaea acidophila; the heat-tolerant Aureobasidium namibiae; the ethanol vapor and heat-tolerant species Baudoinia compniacensis and Zasmidium cellare; and the salt-tolerant species Hortaea werneckii and Cladosporium sphaerospermum. Furthermore, six species of lichens (Acarospora strigata, Cladonia grayi, Dibaeis baeomyces, Graphis scripta, Umbilicaria pustulata and Lobaria pulmonaria) from the class Lecanoromycetes, and one endophytic species (Xylona heveae) from the Xylonomycetes, were included as outgroups.

The genome sizes of fungal species in the study were calculated, and an ANOVA test was conducted using the “car” package in R (Fox and Weisberg, 2019). To determine statistically significant differences between the genome sizes of fungal groups, Tukey’s test was subsequently performed using the “stats” base package in R (R Core Team, 2024).

2.3 GC contents of the genomes, coding sequences and third-codon position

GC content was analyzed using the stats.sh function of bbtools (Bushnell, 2014) for determining the overall GC contents of the whole-genome sequences, and the EMBOSS application “cusp” (Rice et al., 2000) for calculating the GC content of protein-coding genes and the third-codon position (GC3). GC contents of genomic contigs and scaffolds of C. antarcticus and C. minteri were calculated using countgc.sh of the bbtools package (Bushnell, 2014). Plots representing the GC distribution were created using ggplot2 (Wickham, 2016) in R. Statistically significant differences in genome-wide GC content, GC3 and GC content of protein-coding genes between the groups of fungi were studied with a MANOVA test using the “car” (Fox and Weisberg, 2019) package in R. Tukey’s test was conducted for multiple group comparison using the base R “stats” package.

2.4 Genome-wide analysis of repeat-induced point mutation

Genome assemblies of species included in this study were used to conduct RIP analyses using the RIPper software (https://github.com/TheRIPper-Fungi/TheRIPper/) (van Wyk et al., 2019). This program calculates the RIP indices and the average GC contents of whole-genome sequence information using a sliding-window approach. Large RIP-affected genomic regions (LRARs) and the index values associated with them are also determined. We used the default parameters for genome-wide RIP index calculations consisting of 1000-bp windows with a 500-bp step size. Stringent parameters were used (product index value > 1.15; substrate index value < 0.75; RIP composite index > 0) to identify genomic regions affected by RIP. RIPper considers a window as RIP positive if all three indices indicate RIP activity. Only regions of more than 4000 consecutive bp that are affected by RIP are considered as LRARs (van Wyk et al., 2019). The total estimated genome-wide RIP percentage and the composite index for LRARs are shown along with a rooted species tree that was constructed using OrthoFinder v. 2.5.5 (Emms and Kelly, 2019) and the proteome files of each species. The species phylogenetic tree along with the RIP analysis results were visualized using the Interactive Tree of Life (iTOL v. 5.0; Letunic and Bork, 2021).

2.5 Identification of orthogroups

OrthoFinder v. 2.5.5 (Emms and Kelly, 2019) was used to identify orthologous gene families using DIAMOND (Buchfink et al., 2021) and MCL graph clustering (Enright et al., 2002) for sequence comparison and STAG (Emms and Kelly, 2018) as the default species tree method. For each of the 25,761 orthogroups, gene counts per species were normalized by converting to standard normal deviates using the 50% trimmed mean and standard deviation of each orthogroup. Expanded and contracted orthogroups were identified by comparing the normalized scores for the 11 psychrophiles and the 41 non-psychrophiles using Welch’s t-test statistic,

$$\mathit{\text{t=}}\frac{c_p-{\mathit{\text{c}}}_{\mathit{\text{n}}}}{\sqrt{{\mathit{\text{s}}}_{\mathit{\text{p}}}^2+{\mathit{\text{s}}}_{\mathit{\text{n}}}^2}}$$

where c is the mean normalized count, s is the standard deviation, and n is the number of sequences. Subscripts p and n indicate psychrophiles and non-psychrophiles, respectively. The 50 most expanded and 50 most contracted groups were selected for further analysis.

2.6 Correlational matrix and heatmap

The R built-in stats functions “cor” (R Core Team, 2023) and “reshape2” (Wickham, 2007) packages were used to build a correlation matrix using the normalized counts. We used the package “factoextra” (Kassambara and Mundt, 2020) in R to conduct a principal component analysis (PCA) using the correlational matrix to determine clusters of species explained by the orthogroups. The packages “ggplot2” and “ggcorrplot” (R Core Team, 2023) were used to visualize the associations in the matrix.

The R package pheatmap (v1.0.12) (Kolde, 2019) was used to construct a heatmap of the selected expanded and contracted orthogroups. Correlation was used as the distance measure for both orthogroups and species, and groups were clustered using Ward’s method (Ward Jr., 1963).

2.7 Functional annotation of orthogroups

Orthogroup annotations were obtained by Blast comparison to the UniRef50 database (The UniProt Consortium, 2023) with a maximum E-value of 10⁻¹⁰ and minimum of 40% identity. To identify domain annotations and GO terms associated with orthogroups, we used Interproscan (v 5.66-98.0) and InterPro (v 98.0) (Paysan-Lafosse et al., 2023). CAZymes were annotated using dbCAN2 (Yin et al., 2012).

3 Results

3.1 Genome features of psychrophilic fungi

Genome sizes of the psychrophilic fungi range from 23.88 to 51.66 Mbp, with C. antarcticus being the largest in the list (Figure 1). The number of encoded proteins varies from 8,778 to 18,892, with 9,404 in C. antarcticus (Figure 2). Repeat content in the two Cryomyces species varied from 8.4% in C. minteri to 12.5% in C. antarcticus, or about 2 Mbp, which also helps to explain the disparity in genome sizes between these species. A VISTA dot plot of the C. antarcticus contigs revealed numerous diagonal lines (Figure 3; Supplementary Figure 1) indicating that the sequenced isolate may be a polymorphic dikaryon or possibly a diploid. The VISTA dot plot displays blue and red lines, which represent off-diagonal matching contigs or their parts, indicating a partially duplicated genome (Supplementary Figure 1). Scaffold 6 comprises entire scaffolds 62 and 67, while scaffold 5 is composed of scaffold 61 (label not shown but it precedes 62), inverted scaffold 66, and scaffold 72 (Figure 3).

Figure 1

Figure 1. Genome size distribution of fungi with different lifestyles. Horizontal lines indicate the means, boxes show the group sizes, the vertical lines represent the range of the data, and the dots are outliers. Lifestyle groups are color coded according to the legend on the right. Seven fungal species that belong to both endolithic and psychrophilic lifestyles are represented in both groups in the figure. The lifestyle Acidophilic + Heat tolerant encompasses the two ethanol vapor-tolerant fungi included in the study.

Figure 2

Figure 2. Summary information about the genomes of extremophilic fungi in MycoCosm, (A) Genome sizes compared to the numbers of encoded proteins for the extremophilic fungi. Each group is color coded according to the legend on the right. Species framed indicate the psychrophiles that are also endolithic. (B, C) MycoCosm group pages for for (B) psychrophilic (https://mycocosm.jgi.doe.gov/Psychrophilic_fungi) and (C) endolithic (https://mycocosm.jgi.doe.gov/Endolithic_fungi) groups of fungi with sequenced and annotated genomes showing genome assembly sizes, gene counts, and genome publication status.

Figure 3

Figure 3. VISTA dot plot of selected Cryomyces antarcticus contigs from the whole-genome DNA-based self alignment (Dubchak, 2007) – see Supplementary Figure 1 for complete picture. The red and blue diagonal lines are nucleotide matches between different contigs, suggesting DNA duplication.

To gain additional insights into the likely diploidy or possible dikaryotic ploidy of C. antarcticus, we generated a k-mer plot using a k-mer length of 31 (Supplementary Figure 2). The k-mer frequency plot reveals two peaks, providing further evidence of diploidy or dikaryosis, or even a recent whole-genome duplication. The taller peak indicated by a long arrow represents k-mers found at a higher frequency (Supplementary Figure 2), common to both chromosome sets. The second, smaller peak indicated by a small arrow contains k-mers that are unique to one of the homologous chromosomes or have a lower frequency because they are not duplicated across the homologous chromosomes (Ranallo-Benavidez et al., 2020). The small peaks observed in the k-mer frequency plot normally indicate repeats in the genome (Lo and Chain, 2014).

We also compared the genome sizes of the psychrophilic species with those of other lifestyle groups included in our analysis (Figure 1). The comparison revealed significant variations in genome size within individual lifestyles. Plant pathogens generally exhibit the largest genome sizes among the Dothideomycetes species included in the study. Notably, psychrophilic (including endolithic) and halotolerant fungi exhibit similar ranges of genome sizes, while other extremophiles classified as acidophilic and heat tolerant have reduced genome sizes compared to the other fungi. Although the ANOVA test revealed significant differences between the lifestyle groups, Tukey’s test did not find any statistically significant differences between specific pairs of lifestyle groups. This suggests that, while there is an overall significant difference between groups, the specific pairwise differences are not strong enough to be considered significant after adjusting for multiple comparisons. Therefore, genome size appears not to be a reliable predictor of a given lifestyle and cannot be used to differentiate psychrophilic fungi from other groups.

There is a positive correlation between genome size and the number of encoded proteins among extremophilic fungi (Figure 2), where larger genome sizes corresponded to higher protein counts. However, there was a departure from this trend for six psychrophilic fungi, which deviated from the linear regression line. Notably, F. endolithicus, Rachicladosporium sp., and R. antarcticum showed higher-than-expected protein counts relative to their genome sizes, while C. antarcticus showed lower-than-expected protein counts. Despite this difference, we did not observe a consistent pattern for genome/proteome ratio among all psychrophilic or extremophilic fungi when compared to those with other lifestyles. This suggests that the discrepancies in genome sizes compared to proteome size for some psychrophilic fungi may be unique to the individual species rather than indicative of general mechanisms of adaptation to extreme environments.

3.2 GC contents in Cryomyces genomes

The C. antarcticus genome exhibited a bimodal distribution of GC content across contigs (Figure 4). The main peak of GC content contains a higher proportion of contigs with GC content values between 42% and 60%, while a small peak of low-GC content between 0% and 32% contains a smaller number of contigs. This may be associated with some residual contamination or the presence of RIP mutations in repetitive regions, giving rise to low GC content. In C. minteri, we did not observe a clear bi-modal distribution of GC content as we did for C. antarcticus, potentially reflecting its lower proportion of repetitive DNA. We then analyzed the genome-wide RIP mutations in the genomes and found that both C. minteri and C. antarcticus showed a similar percentage of genome-wide RIP mutations of 4.36% and 4.04%, respectively, as calculated with the RIPper tool (Figure 5, Supplementary Table 1). These values are the highest for total estimated genome-wide RIP among the psychrophilic species in this study. Taken together, these results suggest that the genomes of these two psychrophilic species might have been affected by RIP to some extent and at a higher level than other psychrophilic fungi. However, once again we did not observe a consistent pattern of genome-wide RIP or associated metrics across all psychrophilic fungal genomes. This suggests that this mechanism may not represent a general strategy for adaptation to cold stress.

Figure 4

Figure 4. GC contents in scaffolds of (A) Cryomyces minteri and (B) C. antarcticus. The y axis “Fraction of scaffolds” represents the fraction of the total complete scaffolds in each species that have a given GC content.

Figure 5

Figure 5. Phylogenetic tree produced by OrthoFinder showing the percentage of each genome affected by RIP and the composite index for large RIP-affected genomic regions (LRARs) in each species. Tree leaves are labeled by JGI portal ids (see Table 1 for full species names). Extremophilic classification encompasses heat tolerant + halotolerant fungi.

3.3 GC content in whole-genome, protein-coding regions and GC3 is higher in psychrophilic and extremophilic fungi

We grouped extremophilic fungi, excluding psychrophiles, into a single category encompassing acidophilic, heat-tolerant, and halotolerant fungi. Next, we examined whether significant variations exist in GC contents across whole genomes, protein-coding regions, and GC content in the third-codon position (GC3) between psychrophilic fungi and other fungal groups. The overall GC content in C. antarcticus and C. minteri is 53.31%. and 53.89%, respectively. While these values are consistent with those of fungi with other lifestyles, a broader analysis reveals some trends. Specifically, when considering the entire group of psychrophilic and extremophilic fungi, we observed that GC content in the psychrophiles (average = 54.38%) and extremophiles (average=53.47%) is higher than that in plant pathogens (average=50.38%), lichens (average = 47.74%) and fungi with other lifestyles (average = 51.02%). The MANOVA test found significant differences in genome-wide GC content, GC3 and GC content of coding sequences. Tukey’s test revealed significant differences in the GC content of the whole genome, with psychrophiles exhibiting higher GC content compared to plant pathogens and lichens, and the group of other extremophiles also showing significant differences compared to lichens (Figure 6). Further examination using Tukey’s test revealed that the GC contents of coding sequences and the GC3 in psychrophilic fungi are significantly higher than those in lichens, opportunistic fungi, and endophytic fungi, with the latter two categorized as “other” (Figure 6). For extremophiles, significant differences were found in GC contents of coding sequences and GC3 only compared to lichens (Figure 6). Overall, these results suggest that the genome-wide GC content may be a distinctive feature in psychrophilic fungi, potentially contributing to their mechanisms for tolerance to cold environments, as indicated by significant differences observed compared to those for plant pathogens and lichens.

Figure 6

Figure 6. Plot of GC contents for coding sequences (CDS), third-position bases of genes (GC3) and the total genomes (Genome) for fungi in the different lifestyle groups, which are colored according to the legend on the right. Horizontal lines indicate the means, boxes show the group sizes, the vertical lines represent the range of the data, and the dots are outliers. Red star symbols denote groups with statistically significant differences as determined by Tukey’s HSD test. Lifestyle groups are color coded according to the legend on the right. The category ‘Other Extremophiles’ includes acidophilic, heat tolerant and halotolerant fungi.

3.4 OrthoFinder results

OrthoFinder assigned 616,448 genes from 52 species to 25,761 orthogroups. The number of orthogroups with all species present was 1,968, revealing the core genome among the tested fungal species. We identified 40 species-specific orthogroups for C. antarcticus, all of which have at least 2 paralogs in the species, and one orthogroup that has 4 members in C. antarcticus. In the case of C. minteri, we identified 148 species-specific genes, all of which are present in at least 2 copies, and in some cases even 4 to 8 copies. Overall, we found 233 genus-specific orthogroups in Cryomyces.

3.5 Expansion and contraction events in gene families distinguish psychrophilic and extremophilic fungi from those with other lifestyles

The total number of expanded and contracted orthogroups in the fungal species revealed more expansions and fewer contractions in gene families for psychrophilic fungi, while for fungi with other lifestyles, for instance some plant pathogens and lichens, the trend shows more contracted than expanded orthogroups in general (Figure 7). The heatmap clustering analysis reveals three distinct groups within the fungal species, which are due to the correlations between expanded and contracted orthogroups across all 52 species (Figure 8).

Figure 7

Figure 7. Numbers of orthogroups that are expanded (purple bars) or contracted (green bars) relative to a basal point of zero, no expanded/contracted. Groups are considered to be expanded or contracted in t > abs(0.75). Species are indicated by JGI portal ids (see Table 1 for full species names).

Figure 8

Figure 8. Heatmap of the top expanded and contracted orthogroups in the target psychrophilic species compared to the background fungal species.

The first group encompasses the psychrophilic species, which notably are clustered with other extremophilic fungi, including the acidophilic, halotolerant, heat tolerant, and the ethanol vapor-tolerant B. compniacensis. This implies the existence of a common survival strategy for withstanding environmental conditions associated with cold temperatures among psychrophilic species. These strategies parallel those employed by extremophilic fungi to overcome a diverse array of stresses associated with life under extreme conditions. The heatmap shows the 50 most contracted and expanded orthogroups based on t-score values. There are 30 orthogroups that are expanded in the psychrophilic species but generally contracted in the fungal plant pathogens and lichen species. On the other hand, 49 orthogroups are consistently contracted in the psychrophilic species and expanded in the plant pathogens; while 16 out of these 49 orthogroups are expanded in the third group which encompasses mostly the lichens, along with endophytic, opportunistic pathogens, yeasts and molds (Figure 8).

3.6 Functional annotation of expanded gene families in psychrophilic species reveals diversification of mechanisms for cold tolerance

We divided the expanded orthogroups into five general categories (Table 2) based on their predicted protein functions. We found CAZymes, including glycosyl hydrolases and transferases. Furthermore, enzymes with diverse functions underwent expansion, including carboxylases, peptidases, oxidoreductases, alcohol dehydrogenases, phytanoyl-CoA dioxygenase, and hydrolases. Additionally, we found gene families that function as transporters including at least three orthogroups annotated as major facilitator superfamily (MFS) transporters. Some expanded orthogroups in the psychrophilic species have functional annotations related to DNA and RNA regulation, such as helicases, Zn(2)-C6 transcription factor, G-patch proteins, and Gfd2/YDR514C-like proteins. In addition, a fifth group of annotations consists of stress-related proteins, including cell wall integrity and stress response (WSC) family and aquaporins. Finally, a few orthogroups were uncharacterized or contain domains of unknown function.

Table 2

Table 2. Functional annotations of the top 30 orthogroups that were expanded in psychrophilic species based on t-score analysis.

3.7 Gene families contracted in psychrophilic fungi are mostly required for interaction with plants and microorganisms

Annotation analysis of the gene families that were contracted in psychrophilic fungi revealed that the same orthogroups were often expanded in both psychrophiles (section 3.6) and plant pathogens. Our analysis identified a diverse array of annotations, including proteins involved in pathogenicity, stress response, effector production, cell wall integrity, and production of secondary metabolites (Table 3). Interestingly, the Ecp2 effector gene family, which is widely distributed among Dothideomycetes and involved in host-pathogen interactions, is expanded in the plant pathogens but contracted in psychrophiles. Several orthogroups are associated with cell wall integrity, featuring annotations such as cell wall proteins, hydrophobic surface-binding proteins, and acetyltransferases. In terms of stress response, we found gene families annotated as methyltransferases, and NACHT domain-containing proteins. Furthermore, multiple gene families associated with pathogenicity were identified, including F-box proteins, extracellular membrane proteins, oxidases, and enzymes for cell wall degradation and nutrient uptake. A remarkable finding was the expansion of the Zn(2)-C6 fungal-type domain, which can act as an aflatoxin regulatory protein according to its Interpro annotation.

Table 3

Table 3. Functional annotations of the top 49 orthogroups that were contracted in psychrophilic species based on t-score analysis.

4 Discussion

Previous analyses did not reveal much that was unique in the draft genome of a different isolate of C. antarcticus, and concluded with the question, “Nothing special in the specialist?” (Sterflinger et al., 2014). The previously identified haploid genome size and total GC contents were similar to our results. However, that genome was generated with short-read technology that yielded a very fragmented assembly containing over 12,400 contigs, compared to our PacBio assembly of isolate CBS 116301 of only 262 contigs. By greatly increasing the number of genomes from fungi with different lifestyles for comparative analyses it was possible to identify systematically expanded or contracted gene families in C. antarcticus and other psychrophiles that provide clues about their adaptations to extreme cold.

An unexpected result from our analysis is that many of the smaller contigs appear to represent an alternative haplotype to the major contigs, indicating that the sequenced isolate was likely diploid or dikaryotic. This is surprising because these fungi are believed to be haploid and asexual (Nieuwenhuis and James, 2016; Selbmann et al., 2020). The presence of an alternate “haplotigs” set has been reported in phased diploid genome assemblies of plant pathogenic strains of Candida albicans (Hamlin et al., 2019). In this species, considered a strict diploid, the genome assembly is composed of two sets of contigs; a primary contigs set, and an alternate haplotigs set. The latter one represents alternate haplotypes in the regions for which phasing was achieved for the genome assembly. The N50 metric for the primary contigs in C. albicans approached 2 million bp, indicating substantial contiguity. In contrast, the N50 values for the haplotigs ranged from 170,000 bp to 254,000 bp, indicating that contigs representing an alternative haplotype are shorter in length (Hamlin et al., 2019), as was also seen in Cryomyces antarcticus. Additional evidence of diploidy in C. antarcticus was obtained through a k-mer frequency plot, which displayed two major peaks of different heights (Supplementary Figure 2). This bimodal distribution is characteristic of a diploid genome (Ranallo-Benavidez et al., 2020). Similar findings have been observed in the entomopathogenic fungal genus Metarhizium, typically a haploid organism, with an instance of diploidy described in the species M. acridum (Nielsen et al., 2021).

Whether the expanded genome of C. antarcticus represents a stage in the process of sexual reproduction, a balanced heterokaryon, or some other phenomenon such as whole-genome duplication is not known and should be the subject of additional research. Future experiments might include fluorescence flow cytometry analysis on multiple isolates of C. antarcticus. This analysis would determine whether some isolates contain twice the cellular DNA compared to a haploid control strain or a closely related haploid fungal organism (Lin et al., 2009). Additionally, 4,6-diamidino-2-phenylindole (DAPI) nuclear staining using fluorescence microscopy could be applied to test for cells with a single nucleus, confirming diploidy. Alternatively, identifying two independent haploid nuclei would indicate dikaryosis (Lin et al., 2009). The DAPI technique has traditionally been used as a rapid method to determine ploidy in fungal cells and was first assessed in the model yeast Schizosaccharomyces pombe (Talbot et al., 1988). A further way to test for ploidy would be to assess the genome composition of single-ascospore isolates. If single-ascospore isolates have half the DNA content and show no evidence of duplicated sequences then the original isolate was likely diploid. If they have the same genomic composition that would support the hypothesis of a whole-genome, autopolyploid duplication or possibly an allopolyploidization event from hybridization with a closely related species. However, these experiments would depend on generating the sexual stage, which so far remains unknown.

Due to differences in sequencing technologies caution should be applied when comparing contig-level genome assemblies, particularly when assessing features of genome evolution such as RIP, GC contents, and genome size. We observed a significant genome size difference between C. minteri and C. antarcticus, which is larger than expected given their close phylogenetic relationship, but may be explained by differences in sequencing and assembly methods, ploidy, and repeat content. Pepijn and Pellicer (2020) showed that genome assembly size depends on the software used. In that analysis, the genome of two isolates of Leucocoprinus sp. and one of Leucoagaricus barssii were assembled using four software packages (ABySS 2.0.0, SGA 0.10.15, SPAdes 3.7.1 and SOAPdenovo 2.04) showing considerable differences in the total lengths of the assemblies, ranging from 41 to 172 Mb for the same isolate. A similar difference in size was observed between the genomes of the closely related pathogens Epichloe typhina and E. clarkia, where size discrepancies have been attributed to the impact of RIP on genome expansion (Treindl et al., 2021). Repeat-Induced Point mutation (RIP) functions as a defense mechanism in fungal genomes that counteracts the activity of transposable elements to minimize their proliferation (Gladyshev, 2017). This defense strategy induces C-to-T transition mutations within duplicated DNA sequences throughout the genome, resulting in the formation of GC-depleted regions and AT-rich sequences (Hood et al., 2005; Testa et al., 2016). RIP occurs during the sexual stage in haploid nuclei after fertilization but prior to meiotic DNA replication (de Wit et al., 2012). Our analysis revealed traces of RIP in the C. minteri and C. antarcticus genomes, covering roughly 4% of each genome. While the genome-wide RIP percentages observed in both Cryomyces species do not rank among the highest reported for fungi in the class Dothideomycetes (Nagel et al., 2021), they do represent the highest within the psychrophilic fungi in our analysis and indicate that these species must undergo sexual reproduction at least occasionally, since RIP acts only during the sexual cycle (Galagan and Selker, 2004).

In our analysis, we did not investigate the presence of sexual mating type genes. However, a previous study of the C. antarcticus genome by Sterflinger et al. (2014) identified 32 mating genes and 34 mating proteins using BlastP searches of the annotated proteins against all available mating proteins in the GenBank database of fungal species. The homologs of these mating genes were identified in the classes Eurotiomycetes, Sordariomycetes, and Dothideomycetes (Sterflinger et al., 2014). This finding, along with the likely diploidy of C. antarcticus, suggests its potential capability for sexual reproduction. However, to confirm a sexual cycle in this fungus, the functionality of the mating type genes needs further investigation (Usher, 2019). Additionally, we observed a maximum heterozygosity level of 1.15% in C. antarcticus based on the kmer frequency analysis, indicating some degree of genetic recombination and, consequently, possible sexual reproduction assuming the fungus is diploid. If it is not a diploid, then the low level of heterozygosity supports whole-genome duplication by an autopolyploidization event rather than allopolyploidy resulting from hybridization. Despite the evidence of a possible sexual stage, the teleomorph of C. antarcticus remains undiscovered under laboratory conditions. This absence of experimental evidence might be due to the silencing of mating genes or the presence of cryptic sexual cycles triggered by specific environmental conditions, or infrequent sexual reproduction events (Sterflinger et al., 2014; Usher, 2019).

Research on RIP and TE content in psychrophilic and extremophilic fungi is scarce. One study, focusing on the psychrophilic fungus M. psychrophile, revealed minimal repeat sequence content (Su et al., 2016). Currently, most of the research on RIP and TEs predominantly centers on plant pathogens. In fact, among the species included in our analysis, plant pathogens exhibited the highest RIP percentages; this was especially evident in the pathogens P. ulei, P. fijiensis, P. musae, C. fulvum and P. eumusae, which all show more than 35% RIP. A high percentage of RIP has been reported previously in other fungal plant pathogens in the Dothideomycetes (de Wit et al., 2012; Ohm et al., 2012; Isaza et al., 2016).

While research on extremophiles in this area is limited, the impact of RIP has been extensively explored in Dothideomycetes plant pathogens including Z. tritici and its sister species, where a significant proportion of TEs exhibit RIP-like signatures (Lorrain et al., 2021). Analyses of Z. tritici isolates have revealed a higher incidence of RIP and TEs in accessory chromosomes compared to core chromosomes (Badet et al., 2020; Lorrain et al., 2021). Notably, Z. tritici maintains most of its effector genes on accessory chromosomes, often adjacent to TEs, indicating the crucial role of TEs in shaping the evolution of pathogenicity genes in fungi (Lorrain et al., 2021). Drawing parallels with plant pathogens, where key virulence genes are subject to RIP-induced changes, it is conceivable that RIP serves as a driver of gene evolution in extremophilic fungi, facilitating adaptation to extreme conditions. Given the necessity for environmental adaptation in extremophiles, similar to the importance of virulence in plant pathogens, RIP may contribute significantly to the evolutionary dynamics of genes involved in survival mechanisms. However, this hypothesis needs to be evaluated with a deeper analysis of large population-genomics datasets of Cryomyces species to test for RIP and where it occurs relative to chromosomal position.

GC content provides a simple overall indicator of nucleotide composition and can reveal variations in genomic stability and resistance to environmental stressors (Testa et al., 2016; Hu et al., 2022). At low temperatures, base-paired DNA and RNA structures are stabilized. To allow essential processes like replication and transcription to continue normally, one might expect to see lower GC content in psychrophiles. Surprisingly, we observed higher GC contents in the genomes, CDS and GC3s of the psychrophilic fungi compared to mesophiles. A Tukey’s test revealed statistically significant differences in genome-wide GC content between psychrophiles and two major comparison groups: plant pathogens and lichens. This suggests that the GC content of the genome might be a distinctive feature of psychrophilic fungi. Previous analyses reported similar results. For instance, the GC contents in psychrophilic yeasts are higher than they are in those that are mesophilic (Liu et al., 2023). In the fungus Mrakia psychrophila, overall GC content surpassed that of mesophilic fungi such as C. neoformans. This result extended to the GC contents of coding sequences (CDS), which were 56.5% and 56.8% in the two psychrophilic species M. psychrophila and D. cryoxerica, respectively. A similar pattern, with even more pronounced differences, was observed for the third bases of codons (GC3) in both M. psychrophila and D. cryoxerica, compared to C. neoformans (Su et al., 2016).

GC base pairs in DNA are more stable than AT base pairs (See Ghosh et al., 2020a and Mo, 2006). This is not due to hydrogen bonding, which contributes minimally to DNA stability, but rather to Van der Waals and electronic interactions between the stacked base-pairs in the DNA helix (Ghosh et al., 2020a). Entropic effects play a huge role in DNA helix stability (which is why DNA denatures at high temperatures), with the result that DNA helices are more stable at low temperatures. What then could cause the observed increase in GC content in psychrophiles? Recent results show that solvent viscosity decreases helix stability, possibly due to molecular crowding effects (Ghosh et al., 2020b; Stellwagen and Stellwagen, 2019). High viscosity also destabilizes initiation of RNA stem formation and increases the destabilizing effect of terminal AU base-pairs (Ghosh et al., 2020b). Fungal psychrophiles produce cryoprotectants such as glycine betaine, trehalose, glycerol, sucrose, sarcosin, mannitol, sorbitol, 2-polyalkanoates, and exopolysaccharides (Chen et al., 2022; Pacheco et al., 2019; Ghobakhlou et al., 2015; Goordial et al., 2016; Su et al., 2016; Deegenaars and Watson, 1997). Such cryoprotectants depress the freezing point and decrease water loss, both important to survival in the cold, dry Antarctic environment. These compounds probably also increase the viscosity of the intracellular environment, leading to DNA helix destabilization, and destabilization of RNA helix initiation. This may provide an explanation for the otherwise non-intuitive increase in genome-wide GC contents in psychrophiles.

Expansions and contractions of gene families may provide additional clues about adaptations to an extreme environment. For example, in the psychrophilic species we found expansions of helicases and G-patch domain-containing proteins, which are involved in multiple processes for genome regulation. In other organisms, DEAD-box RNA helicases are involved in regulation of nucleic acid binding, modulation of RNA/DNA: DNA interactions, and RNA metabolism, including translation initiation and ribosome biogenesis due to their capability to unwind the duplex RNA in eukaryotes (Lim et al., 2000; Taschuk and Cherry, 2020). In psychrophilic microorganisms RNA helicases are differentially expressed at low temperatures (Feller, 2013) and they can act as nucleic acid chaperones that help the microorganisms adapt to cold conditions (Xu et al., 2024). In particular, DEAD-box RNA helicases can presumably enable some bacteria to survive cold shock and grow at low temperatures, due to their role in removing cold-stabilized secondary structures in mRNA (Lim et al., 2000). In the Antarctic archaeon M. burtonii, an RNA helicase gene was highly transcribed during cell growth at 4°C, revealing the importance of these proteins for low-temperature adaptation (Lim et al., 2000). G-patch proteins are involved in regulation and activation of helicases. The first characterized protein with a G-patch domain that regulates an RNA helicase was Spp2 (Roy et al., 1995). The in vivo interaction of Spp2 was first documented in budding yeast, where a DEAH/RHA helicase named Prp2 interacts physically with Spp2, and this interaction is required for the activation of Prp2 function. Also, a helicase (Prp43) that is required in both splicing and ribosome biogenesis, is activated by a G-patch protein (Ntr1) which forms a stable complex with another G-patch protein (Ntr2). This complex recruits the helicase Prp43 targeting its helicase activity for spliceosome dissociation (Robert-Paganin et al., 2015).

Another gene family that was expanded in the psychrophiles was Zn(2)-C6 transcription factors, which are exclusively found in fungi. They can act as repressors, activators, or both, depending on the target genes (MacPherson et al., 2006). The regulatory scope of these transcription factors is broad, including activation of sterol biosynthesis (Vik and Rine, 2001), modulation of multidrug sensitivity (Coste et al., 2004), and regulation of genes involved in galactose metabolism (Breunig, 2000). Previous analyses have highlighted the versatility of Zn(2)-C6 transcription factors. It is particularly interesting that, in Colletotrichum lagenarium, Cmr1p positively regulates the expression of genes involved in melanin biosynthesis, which is crucial for successful invasion of host plants during anthracnose of cucumbers (Tsuji et al., 2000). In rock-inhabiting psychrophilic fungi, melanin biosynthesis is crucial for survival because melanin pigments protect the cells from extreme temperatures, solar irradiation, extreme osmotic and pH conditions, and provide tolerance to toxic metals (Selbmann et al., 2015). Moreover, Stb5, a Zn(2)-C6 transcription factor in Saccharomyces cerevisiae, can switch between activator and repressor roles under oxidative stress conditions (Larochelle et al., 2006). It binds to and regulates the expression of genes in the pentose phosphate pathway and those involved in NADPH production, which is essential for oxidative stress resistance. In fact, the activation of mechanisms for protection against oxidative stress plays a crucial role in psychrophilic fungi. It has been shown that a cold environment induces the production of reactive oxygen species (ROS), which can cause severe damage to DNA, lipids and proteins, affecting survival of the cell (Kostadinova et al., 2017). The expansion of genes for Zn(2)-C6 transcription factors in C. antarcticus and other psychrophilic extremophiles may help to counter this damage.

Another expanded group includes major facilitator superfamily (MFS) transporters which likely are involved in drug resistance and amino acid uptake. The MFS are secondary transporters that have primary roles in antifungal resistance and nutrition transport. They can indirectly control membrane potential by changing membrane lipid homeostasis, and regulate internal pH and the stress-response machinery in fungi (Su et al., 2016; Chen et al., 2017). MFS transporters are expanded in extremophilic fungal species such as Mrakia psychrophila and Penicillium funiculosum where they are required for accumulation of nutrients from the environment and resistance to acidic stress, respectively (Xu et al., 2014; Su et al., 2016). In Alternaria alternata, an MFS transporter named AaMFS19 is required for cellular resistance to oxidative stress, including H₂O₂, KO₂, and other oxygen-generating compounds (Chen et al., 2017).

Some enzymes play a role in alternative pathways for energy production in fungi. We observed expansion of gene families encoding alcohol dehydrogenases. The expression of this family of enzymes has been reported, along with increased expression of glycolytic genes that feed into ethanol production via fermentation, in the Antarctic yeast Rodothorula frigidialcoholis growing at 0°C (Touchette et al., 2022). This suggests that R. frigidialcoholis adapts to cold temperatures by redirecting pentose phosphate pathway molecules to ethanol fermentation. This switch from respiration to fermentative metabolism to produce energy contributes to its survival in extreme cryoenvironments (Touchette et al., 2022). In Mucor lusitanicus, a knock-out strain of the alcohol dehydrogenase gene (adh1), showed reduction in lipid and fatty acid content, which was restored by supplementing the medium with external ethanol, suggesting that alcohol dehydrogenases are also involved in lipid biosynthesis (Gutiérrez-Corona et al., 2023). The expansion of genes for alcohol dehydrogenases in psychrophiles might be involved either in the activation of alternative pathways for energy production, or in specific stages of lipid biosynthesis, both of which could be alternative mechanisms of adaptation to life in extremely cold environments.

The expansion of CAZymes in the psychrophilic species may facilitate cell wall modifications and energy storage. Fungi living in cold environments maintain cell integrity by remodeling the cell wall. CAZymes control the repair and modification of cell wall components necessary to overcome extreme low temperatures and desiccation. For example, in Candida albicans the architecture of the cell wall undergoes alterations in response to environmental conditions, resulting in changes to its composition and structure in newly synthesized cells (Ene et al., 2015). In previous analyses, extremotolerant species have been shown to encode more members of the CAZyme family glycosyl hydrolase 32 (GH32) compared to opportunistic, non-opportunistic, and pathogenic species (Gostinčar et al., 2018). This family of CAZymes contains invertases and other enzymes involved in energy storage and recovery. In addition, studies have concluded that fungal cell wall elasticity is critical for survival under conditions of osmotic shock (Ene et al., 2015). Structural changes that improve cell wall elasticity involve cross-linking between cell wall macromolecules, which is catalyzed by the action of carbohydrate-active cell wall-remodeling enzymes, including hydrolases, transglycosidases, and transferases. These include the Gas-like family of β-1,3-glucanosyltransferases, and a family of chitin-glucanosyltransferases, which encompasses family GH16 (Ene et al., 2015; Ibe and Munro, 2021). In the fungal pathogen B. cinerea, the GT2 family contributes to the biogenesis and remodeling of the cell wall. This family contains the synthases for production of chitin, which plays a structural role in the fungal cell wall (Blandenet et al., 2022). In our analysis, we observed expansion of gene families of AB hydrolases (OG0006140), glycoside hydrolase family 1 (OG0006802), and the alpha-mannosyltransferase protein (OG0007561), which may be involved in cell wall modifications in response to cold temperatures.

Other expanded families, such as the phytanoyl-CoA dioxygenase proteins, are peroxisomal enzymes that catalyze the first step of phytanic acid alpha-oxidation. In Streptomyces coelicolor, an ectoine dioxygenase is involved in the synthesis of 5-hydroxyectoine, a compatible solute that is synthesized upon exposure to high salinity (Bursy et al., 2008). This compound helps organisms to survive extreme osmotic stress by acting as a highly soluble organic osmolyte. Bursy et al. (2008) also found that synthesis of 5-hydroxyectoine by S. coelicolor is a stress response that protects against conditions of high salt and heat.

In contrast to gene families that were expanded in the psychrophiles, others were reduced, especially in comparison to the plant pathogens, in which gene families involved in microbiome interactions and defense or immune responses often were expanded. These families were contracted in psychrophilic fungi, most likely because plant pathogens require more proteins associated with microorganism and plant interactions, or interactions with other living things in general, while psychrophilic fungi are more isolated and face bigger challenges from abiotic and environmental factors rather than from other living organisms.

Most of the orthogroups that are contracted in the psychrophiles are proteins related to signaling and the membrane. They show diverse functions from degradation, response to stress, and synthesis of new compounds. Acetyltransferase (GNAT), heme haloperoxidase, and peroxidase are common domains in the contracted orthogroups. These proteins are involved in lignin and chitin degradation, which mostly are not used by cold-tolerant fungi (Leung et al., 2011). However, the insect pathogen Lecanicillium muscarium, which was isolated from Antarctica, is a high producer of cold-tolerant enzymes, including those for hydrolyzing chitin, as expected for a fungus that feeds on arthropods (Fenice, 2016). Several fungi isolated from King George Island, Antarctica, including the Dothideomycetes Penidiella kurandae, showed amylase and cellulase activities, even when grown at 4°C (Krishnan et al., 2016).

Other reduced orthogroups in the psychrophilic fungi are associated with virulence/pathogenicity and toxin production in plant pathogens such as polyketide synthase, glycoside hydrolase, ribonuclease ribotoxin, aflatoxin regulatory protein, and pyrroline-5-carboxylate reductase. A previous comparative genomic analysis of 20 dothideomycetous and eurotiomycetous black fungi showed a lack of specialized virulence traits in extremotolerant and opportunistic fungi (Gostinčar et al., 2018). This is consistent with our observed contraction of these orthogroups.

We observed expansion of the NACHT-NTPase and P-loop NTPase N-terminal domain gene families in the group of plant pathogens, which suggests involvement of these proteins in controlling a variety of biotic interactions, not strictly limited to immune responses. The nucleotide-binding oligomerization domain (NOD)-like receptors (NLRs) are signal-transducing proteins that control innate immunity in plants, animals and fungi (Dyrka et al., 2014). In fungi, NACHT domain proteins, a subtype of NLRs, control immunity and other nonpathogenic biotic interactions, such as symbiotic relationships with microbiomes (Chu and Mazmanian, 2013). The HET-E protein of Podospora anserina, which first defined the NACHT domain, is involved in fungal non-self recognition and programmed cell death due to heterokaryon incompatibility (Saupe et al., 1995; Koonin and Aravind, 2000). Additionally, studies have shown that NACHT domain proteins are specifically expressed during mycorrhizal symbiosis in certain fungi (Martin et al., 2008). These interactions could include fungal pathogenicity and symbiotic interactions (such as ECM formation, endophytic growth, lichen formation, or microbiome interactions).

In support of this idea, some studies suggest that variation in the number of NLR homologs among different fungal species might be associated with the diversity of ecological niches and lifestyles (Dyrka et al., 2014). For example, highly versatile pathogens like Fusarium species or mycoparasitic Trichoderma species tend to have large NLR repertoires, possibly for protection against microbial competitors. However, the ethanol-tolerant B. compniacensis only has a small repertoire of these NLRs, which might result from inhabiting a restrictive niche (Dyrka et al., 2014). Similarly, our results suggest that extremophilic fungi which inhabit environments with a restricted number of microbial competitors and potential hosts, might not require specialized mechanisms for biotic interactions, which explains the contraction of the NACHT domain-containing protein family.

The contraction of hydrophobic surface-binding proteins among the psychrophiles may be due to the involvement of hydrophobins in water-sensing mechanisms during spore germination. Hydrophobins contribute to the surface hydrophobicity needed for conidial dispersal, and protection against a host defense system (Tanaka et al., 2022), all of which are essential processes for pathogenic or symbiotic interactions but are not required in non-pathogenic, psychrophilic fungi. In plant pathogens like M. grisea and F. graminearum, hydrophobins are required for production of infective structures, penetration of the water-air interface of the host and attachment to hydrophobic surfaces (Talbot et al., 1993; Quarantin et al., 2019). These properties are not required by saprobic, extremophilic psychrophiles and this is reflected in their proteome complements.

In addition to the gene families discussed above, we found that genes involved in stress responses against UV irradiation, oxidative stress, and antifungal agents, are also contracted. This finding is remarkable given that the Antarctic continent receives the highest levels of UV radiation on Earth. In 2020, the continent experienced the highest UV irradiances in more than two decades, a consequence of the recurring Antarctic ozone hole that appears every spring (Cordero et al., 2022). The endolithic lifestyle of most psychrophilic fungi in this study might explain the contraction of gene families involved in stress responses to UV radiation. The lithic environment provides a refuge against harsh weather conditions such as high UV radiation and water scarcity. For example, expansion of hyphal growth into deeper rock layers helps to protect the cells of rock-inhabiting fungi from UV radiation (Gorbushina, 2007). Consequently, psychrophiles living as endolithics can withstand severe radiation levels and other stress forms. Indeed, endolithic communities consist of some of the most extremophilic and resilient microorganisms known to date (Coleine et al., 2021). These microorganisms migrate into fractures or in pore spaces where the necessary nutrient, moisture, and light levels are sufficient for survival (Stivaletta and Barbieri, 2008). Thus, the endolithic lifestyle not only provides protection against radiation but also serves as a strategy to avoid water stress and access a more humid environment in some of the driest and coldest places on Earth. A future analysis that includes a broader range of extremophilic fungal species with both epilithic and endolithic lifestyles is of interest to determine the differences in gene families that are expanded and contracted in endolithics compared to epilithics. We likely would not observe the contraction of gene families associated with UV stress response in a group predominantly composed of epilithic or epiphytic psychrophiles. This finding suggests that evolution and adaptation may have more to do with a dynamic and diversified set of strategies, as opposed to using only a specific set of proteins and mechanisms to cope with extreme environments.

While at first glance there was nothing special in the specialist (Sterflinger et al., 2014), upon a closer look with a larger cohort of fungi plus a much-improved genome sequence, expanded and contracted gene families could be identified and provide clues about likely adaptations to extreme conditions. This analysis highlights the importance of a large number of genomes available for comparison and validates the goals of programs like the JGI 1000 Fungal Genomes project (Grigoriev et al., 2014) which seeks to generate the requisite sequences. In a sense, Sterflinger et al. (2014) were correct in that there were not radical changes in gene composition that distinguish psychrophiles from other fungi. Instead, the changes most likely reflect a slow adaptation of C. antarcticus and its relatives over long periods of evolutionary time as gene families expanded or contracted and homologs changed function as they were selected for increasing tolerance to cold stress. Clearly much additional work is needed to understand the expression of the identified gene families under extreme conditions to fully elucidate the ability of Cryomyces to grow under conditions that are more similar to those found on Mars (Onofri et al., 2004; 2015) than they are to most areas of Earth.

Data availability statement

The datasets presented in this study can be found in online repositories. This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession JBDRTV000000000. The version described in this paper is version JBDRTV010000000.

Author contributions

SG-G: Writing – original draft, Formal analysis, Data curation, Visualization, Investigation, Methodology. WS-H: Writing – original draft, Methodology, Investigation, Formal analysis, Data curation. SH: Writing – review & editing, Resources, Methodology, Investigation, Formal analysis, Data curation. KL: Writing – review & editing, Resources, Formal analysis. JE: Writing – review & editing, Resources. NK: Writing – review & editing, Resources. AL: Formal analysis, Writing – review & editing, Resources. KB: Project administration, Writing – review & editing. SG: Methodology, Investigation, Conceptualization, Supervision, Writing – review & editing, Writing – original draft, Project administration. MG: Writing – original draft, Visualization, Supervision, Formal analysis, Data curation, Methodology, Investigation. IG: Writing – review & editing, Resources, Project administration, Conceptualization, Visualization, Supervision, Methodology, Investigation.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded in part by USDA-ARS research projects 5020-21220-019-00D and 5020-21220-014-00D and by DOE under contract DE-AC02-05CH11231.

Acknowledgments

We thank Dr. Filipe Victoria, Dr. Fabiola Lucini, and Dr. Francis Martin for providing access to the unpublished genome data of Alternaria sp. UNIPAMPA012 and Alternaria sp. UNIPAMPA017 produced within the 1000 Fungal Genomes project, Dr. Laura Selbmann, the Italian National Program for Antarctic Researches (PNRA) and the Italian National Antarctic Museum (MNA) for providing access to the unpublished genome data of Extremus antarcticus MNA-CCFEE and Capnodiales sp, Dr. Elizabeth Arnold for providing access to the unpublished genome data of Myriangiaceae sp. NC1570 and Teratosphaeriaceae sp. NC1134, Dr. Ólafur Andrésson for providing access to the unpublished genome data of Lobaria pulmonaria, and Dr. Eduardo Mizubuti for providing access to the unpublished genome data of Pseudocercospora ulei. The genome sequence data were generated by the US Department of Energy Joint Genome Institute in collaboration with the user community. The work (proposal: 10.46936/10.25585/60001019) was conducted by the U.S. Department of Energy Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science User Facility, supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231. We thank Jessica R. Cavaletto for the time-consuming process of extracting nucleic acids from the extremely slow growing and uncooperative C. antarcticus. We also thank the reviewers for their very helpful suggestions for improving the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ffunb.2024.1418145/full#supplementary-material

Supplementary Figure 1 | VISTA dot plot of Cryomyces antarcticus contigs against themselves based on the whole-genome DNA alignment (Dubchak, 2007). The red and blue diagonal lines are nucleotide matches between sequences on different contigs, indicating DNA duplication. The main diagonal of self alignment has been removed for clarity. A higher quality version of the supplementary figure can be accessed at https://mycocosm.jgi.doe.gov/vista_embed/?viewMode=dotPlot&organism=Cryan3&?&run=8474-crZbaJ90&xdset=3156&ydset=3714&cutoff=1000

Supplementary Figure 2 | K-mer frequency plot of the circular consensus sequencing data from Cryomyces antarcticus using a k-mer length of 31 nucleotides. Two peaks indicating diploidy or dikaryosis are shown: a taller peak with high-frequency k-mers and a smaller peak with low-frequency k-mers that are unique to one heterologous chromosome. The other small peaks possibly indicate repeats in the genome.

Supplementary Table 1 | Repeat-induced point mutation (RIP) in the genomes used for comparative analysis.

References

Arango Isaza R. E., Diaz-Trujillo C., Dhillon B., Aerts A., Carlier J., Crane C. F., et al. (2016). Combating a Global Threat to a Clonal Crop: Banana Black Sigatoka Pathogen Pseudocercospora fijiensis (Synonym Mycosphaerella fijiensis) Genomes Reveal Clues for Disease Control. PLoS Genet. 12 (8), e1005876. doi: 10.1371/journal.pgen.1005876. Erratum in: PLoS Genet. 2016 Oct 4;12(10):e1006365. doi: 10.1371/journal.pgen.1006365

PubMed Abstract | Crossref Full Text | Google Scholar

Armaleo D., Müller O., Lutzoni F., Andrésson Ó. S., Blanc G., Bode H. B., et al. (2019). The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata. BMC Genomics 20 (1), 605. doi: 10.1186/s12864-019-5629-x

PubMed Abstract | Crossref Full Text | Google Scholar

Badet T., Oggenfuss U., Abraham L., McDonald B.A., Croll D. (2020). A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici. BMC Biol. 18, 1–18. doi: 10.1186/s12915-020-0744-3

PubMed Abstract | Crossref Full Text | Google Scholar

Blandenet M., Gonçalves I. R., Rascle C., Dupuy J.-W., Gillet F.-X., Poussereau N., et al. (2022). Evidencing new roles for the glycosyl-transferase cps1 in the phytopathogenic fungus botrytis cinerea. J. Fungi (Basel) 8, 899. doi: 10.3390/jof8090899

PubMed Abstract | Crossref Full Text | Google Scholar

Breunig K. D. (2000). ‘Regulation of Transcription Activation by Gal4p’. Food technol. biotechnol. 38 (4), 287–293.

Google Scholar

Buchfink B., Reuter K., Drost H.-G. (2021). Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 18, 366–368. doi: 10.1038/s41592-021-01101-x

PubMed Abstract | Crossref Full Text | Google Scholar

Bursy J., Kuhlmann A. U., Pittelkow M., Hartmann H., Jebbar M., Pierik A. J., et al. (2008). Synthesis and uptake of the compatible solutes ectoine and 5-hydroxyectoine by streptomyces coelicolor A3(2) in response to salt and heat stresses. Appl. Environ. Microbiol. 74, 7286–7296. doi: 10.1128/AEM.00768-08

PubMed Abstract | Crossref Full Text | Google Scholar

Bushnell B. (2014). “BBTools software package.”. Available at: https://sourceforge.net/projects/bbmap/.

Google Scholar

Chang T. C., Salvucci A., Crous P. W., Stergiopoulos I. (2016). Comparative Genomics of the Sigatoka Disease Complex on Banana Suggests a Link between Parallel Evolutionary Changes in Pseudocercospora fijiensis and Pseudocercospora eumusae and Increased Virulence on the Banana Host. PLoS Genet. 12 (8), e1005904. doi: 10.1371/journal.pgen.1005904

PubMed Abstract | Crossref Full Text | Google Scholar

Chen L.-H., Tsai H.-C., Yu P.-L., Chung K.-R. (2017). A major facilitator superfamily transporter-mediated resistance to oxidative stress and fungicides requires yap1, skn7, and MAP kinases in the citrus fungal pathogen alternaria alternata. PloS One 12, e0169103. doi: 10.1371/journal.pone.0169103

PubMed Abstract | Crossref Full Text | Google Scholar

Chen W., Zhang X., Li S., Cui J., Yang X., Zhang Q. (2022). The MAP-kinase HOG1 controls cold adaptation in rhodosporidium kratochvilovae by promoting biosynthesis of polyunsaturated fatty acids and glycerol. Curr. Microbiol. 79. doi: 10.1007/s00284-022-02957-8

Crossref Full Text | Google Scholar

Chu H., Mazmanian S. K. (2013). Innate immune recognition of the microbiota promotes host-microbial symbiosis. Nat. Immunol. 14, 668–675. doi: 10.1038/ni.2635

PubMed Abstract | Crossref Full Text | Google Scholar

Coleine C., Masonjones S., Selbmann L., Zucconi L., Onofri S., Pacelli C., et al. (2017). Draft Genome Sequences of the Antarctic Endolithic Fungi Rachicladosporium antarcticum CCFEE 5527 and Rachicladosporium sp. CCFEE 5018. Genome Announc. 5 (27), e00397-17. doi: 10.1128/genomeA.00397-17

PubMed Abstract | Crossref Full Text | Google Scholar

Coleine C., Masonjones S., Sterflinger K., Onofri S., Selbmann L., Stajich J. E. (2020). Peculiar genomic traits in the stress-adapted cryptoendolithic Antarctic fungus Friedmanniomyces endolithicus. Fungal Biology Fungal Adaption to Stress 124, 458–467. doi: 10.1016/j.funbio.2020.01.005

Crossref Full Text | Google Scholar

Coleine C., Stajich J. E., de los Ríos A., Selbmann L. (2021). Beyond the extremes: Rocks as ultimate refuge for fungi in drylands. Mycologia 113, 108–133. doi: 10.1080/00275514.2020.1816761

PubMed Abstract | Crossref Full Text | Google Scholar

Cordero R. R., Feron S., Damiani A., Redondas A., Carrasco J., Sepúlveda E., et al. (2022). Persistent extreme ultraviolet irradiance in Antarctica despite the ozone recovery onset. Sci. Rep. 12, 1266. doi: 10.1038/s41598-022-05449-8

PubMed Abstract | Crossref Full Text | Google Scholar

Coste A. T., Karababa M., Ischer F., Bille J., Sanglard D. (2004). TAC1, transcriptional activator of CDR genes, is a new transcription factor involved in the regulation of candida albicans ABC transporters CDR1 and CDR2. Eukaryot Cell 3, 1639–1652. doi: 10.1128/EC.3.6.1639-1652.2004

PubMed Abstract | Crossref Full Text | Google Scholar

Dal Grande F., Sharma R., Meiser A., Rolshausen G., Büdel B., Mishra B., et al. (2017). Adaptive differentiation coincides with local bioclimatic conditions along an elevational cline in populations of a lichen-forming fungus. BMC Evol. Biol. 17, 93. doi: 10.1186/s12862-017-0929-8

PubMed Abstract | Crossref Full Text | Google Scholar

Deegenaars M. L., Watson K. (1997). Stress proteins and stress tolerance in an Antarctic, psychrophilic yeast, Candida psychrophila. FEMS Microbiol. Lett. 151, 191–196. doi: 10.1111/j.1574-6968.1997.tb12569.x

PubMed Abstract | Crossref Full Text | Google Scholar

de Jonge R., Ebert M. K., Huitt-Roehl C. R., Pal P., Suttle J. C., Spanner R. E., et al. (2018). Gene cluster conservation provides insight into cercosporin biosynthesis and extends production to the genus Colletotrichum. Proc. Natl. Acad. Sci. U.S.A. 115, E5459–E5466. doi: 10.1073/pnas.1712798115

PubMed Abstract | Crossref Full Text | Google Scholar

de Wit P. J. G. M., van der Burgt A., Ökmen B., Stergiopoulos I., Abd-Elsalam K. A., Aerts A. L., et al. (2012). The Genomes of the Fungal Plant Pathogens Cladosporium fulvum and Dothistroma septosporum Reveal Adaptation to Different Hosts and Lifestyles But Also Signatures of Common Ancestry. PloS Genet. 8, e1003088. doi: 10.1371/journal.pgen.1003088

PubMed Abstract | Crossref Full Text | Google Scholar

Duarte A. W. F., dos Santos J. A., Vianna M. V., Vieira J. M. F., Mallagutti V. H., Inforsato F. J., et al. (2018). Cold-adapted enzymes produced by fungi from terrestrial and marine antarctic environments. Crit. Rev. Biotechnol. 38, 600–619. doi: 10.1080/07388551.2017.1379468

PubMed Abstract | Crossref Full Text | Google Scholar

Dubchak I. (2007). Comparative analysis and visualization of genomic sequences using VISTA browser and associated computational tools. Methods Mol. Biol. 395, 3–16. doi: 10.1007/978-1-59745-514-5_1

PubMed Abstract | Crossref Full Text | Google Scholar

Dyrka W., Lamacchia M., Durrens P., Kobe B., Daskalov A., Paoletti M., et al. (2014). Diversity and variability of NOD-like receptors in fungi. Genome Biol. Evol. 6, 3137–3158. doi: 10.1093/gbe/evu251

PubMed Abstract | Crossref Full Text | Google Scholar

Emms D. M., Kelly S. (2018). STAG: species tree inference from all genes. bioRxiv p, 267914. doi: 10.1101/267914

Crossref Full Text | Google Scholar

Emms D. M., Kelly S. (2019). OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20, 238. doi: 10.1186/s13059-019-1832-y

PubMed Abstract | Crossref Full Text | Google Scholar

Ene I. V., Walker L. A., Schiavone M., Lee K. K., Martin-Yken H., Dague E., et al. (2015). Cell wall remodeling enzymes modulate fungal cell wall elasticity and osmotic stress resistance. mBio 6, e00986–e00915. doi: 10.1128/mBio.00986-15

PubMed Abstract | Crossref Full Text | Google Scholar

Enright A. J., Van Dongen S., Ouzounis C. A. (2002). An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 30, 1575–1584. doi: 10.1093/nar/30.7.1575

PubMed Abstract | Crossref Full Text | Google Scholar

Feller G. (2013). Psychrophilic enzymes: from folding to function and biotechnology. Scientifica, 512840. doi: 10.1155/2013/512840

PubMed Abstract | Crossref Full Text | Google Scholar

Fenice M. (2016). The psychrotolerant antarctic fungus lecanicillium muscarium CCFEE 5003: A powerful producer of cold-tolerant chitinolytic enzymes. Molecules 21, 447. doi: 10.3390/molecules21040447

PubMed Abstract | Crossref Full Text | Google Scholar

Fox J., Weisberg S. (2019). An R Companion to Applied Regression, Third edition (Thousand Oaks CA: Sage). Available at: https://socialsciences.mcmaster.ca/jfox/Books/Companion/.

Google Scholar

Friedmann E. I. (1982). Endolithic microorganisms in the antarctic cold desert. Science 215, 1045–1053. doi: 10.1126/science.215.4536.1045

PubMed Abstract | Crossref Full Text | Google Scholar

Galagan J. E., Selker E. U. (2004). RIP: the evolutionary cost of genome defense. Trends Genet. 20, pp.417–pp.423. doi: 10.1016/j.tig.2004.07.007

Crossref Full Text | Google Scholar

Gazis R., Kuo A., Riley R., LaButti K., Lipzen A., Lin J., et al. (2016). The genome of xylona heveae provides a window into fungal endophytism. Fungal Biol. 120, 26–42. doi: 10.1016/j.funbio.2015.10.002

PubMed Abstract | Crossref Full Text | Google Scholar

Ghobakhlou A. F., Johnston A., Harris L., Antoun H., Laberge S. (2015). Microarray transcriptional profiling of Arctic Mesorhizobium strain N33 at low temperature provides insights into cold adaption strategies. BMC Genomics 16. doi: 10.1186/s12864-015-1611-4

Crossref Full Text | Google Scholar

Ghosh S., Takahashi S., Banerjee D., Ohyama T., Endoh T., Tateishi-Karimata H., et al. (2020b). Nearest-neighbor parameters for the prediction of RNA duplex stability in diverse in vitro and cellular-like crowding conditions. Nucleic Acids Res. 51, 4101–4111. doi: 10.1093/nar/gkad020

Crossref Full Text | Google Scholar

Ghosh S., Takahashi S., Ohyama T., Endoh T., Tateishi-Karimata H., Sugimoto N. (2020a). Nearest-neighbor parameters for predicting DNA duplex stability in diverse molecular crowding conditions. Proc. Natl. Acad. Sci. U S A. 117, 14194–14201. doi: 10.1073/pnas.1920886117

PubMed Abstract | Crossref Full Text | Google Scholar

Gladyshev E. (2017). Repeat-induced point mutation (RIP) and other genome defense mechanisms in fungi. Microbiol. Spectr. 5. doi: 10.1128/microbiolspec.FUNK-0042–2017

PubMed Abstract | Crossref Full Text | Google Scholar

Grandaubert J., Bhattacharyya A., Stukenbrock E. H. (2015). RNA-seq-Based Gene Annotation and Comparative Genomics of Four Fungal Grass Pathogens in the Genus Zymoseptoria Identify Novel Orphan Genes and Species-Specific Invasions of Transposable Elements. G3 (Bethesda) 5 (7), 1323–33. doi: 10.1534/g3.115.017731

PubMed Abstract | Crossref Full Text | Google Scholar

Goodwin S. B., M’Barek S. B., Dhillon B., Wittenberg A. H. J., Crane C. F., Hane J. K., et al. (2011). Finished genome of the fungal wheat pathogen mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PloS Genet. 7, 1002070. doi: 10.1371/journal.pgen.1002070

Crossref Full Text | Google Scholar

Goordial J., Raymond-Bouchard I., Zolotarov Y., de Bethencourt L., Ronholm J., Shapiro N., et al. (2016). Cold adaptive traits revealed by comparative genomic analysis of the eurypsychrophile Rhodococcus sp. JG3 isolated from high elevation McMurdo Dry Valley permafrost, Antarctica. FEMS Microbiol. Ecol. 92, fiv154. doi: 10.1093/femsec/fiv154

PubMed Abstract | Crossref Full Text | Google Scholar

Gorbushina A. A. (2007). Life on the rocks. Environ. Microbiol. 9, 1613–1631. doi: 10.1111/j.1462-2920.2007.01301.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gostinčar C., Ohm R. A., Kogej T., Sonjak S., Turk M., Zajc J., et al. (2014). Genome sequencing of four Aureobasidium pullulans varieties: biotechnological potential, stress tolerance, and description of new species. BMC Genomics 15, 549. doi: 10.1186/1471-2164-15-549

PubMed Abstract | Crossref Full Text | Google Scholar

Gostinčar C., Stajich J. E., Gunde-Cimerman N. (2023). Extremophilic and extremotolerant fungi. Curr. Biol. 33, R752–R756. doi: 10.1016/j.cub.2023.06.011

PubMed Abstract | Crossref Full Text | Google Scholar

Gostinčar C., Zajc J., Lenassi M., Plemenitaš A., de Hoog S., Al-Hatmi A. M. S., et al. (2018). Fungi between extremotolerance and opportunistic pathogenicity on humans. Fungal Diversity 93, 195–213. doi: 10.1007/s13225-018-0414-8

Crossref Full Text | Google Scholar

Grigoriev I. V., Nikitin R., Haridas S., Kuo A., Ohm R., Otillar R., et al. (2014). MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res. 42, D699–D704. doi: 10.1093/nar/gkt1183

PubMed Abstract | Crossref Full Text | Google Scholar

Gutiérrez-Corona J. F., González-Hernández G. A., Padilla-Guerrero I. E., Olmedo-Monfil V., Martínez-Rocha A. L., Patiño-Medina J. A., et al. (2023). Fungal alcohol dehydrogenases: physiological function, molecular properties, regulation of their production, and biotechnological potential. Cells 12, 2239. doi: 10.3390/cells12182239

PubMed Abstract | Crossref Full Text | Google Scholar

Hamlin J. A. P., Dias G. B., Bergman C. M., Bensasson D. (2019). Phased diploid genome assemblies for three strains of candida albicans from oak trees. G3: Genes|Genomes|Genetics , 9, 3547–3554. doi: 10.1534/g3.119.400486

Crossref Full Text | Google Scholar

Haridas S., Albert R., Binder M., Bloem J., LaButti K., Salamov A., et al. (2020). 101 dothideomycetes genomes: A test case for predicting lifestyles and emergence of pathogens. Stud. Mycol 96, 141–153. doi: 10.1016/j.simyco.2020.01.003

PubMed Abstract | Crossref Full Text | Google Scholar

Hood M. E., Katawczik M., Giraud T. (2005). Repeat-induced point mutation and the population structure of transposable elements in microbotryum violaceum. Genetics , 170, 1081–1089. doi: 10.1534/genetics.105.042564

PubMed Abstract | Crossref Full Text | Google Scholar

Hu E.-Z., Lan X.-R., Liu Z.-L., Gao J., Niu D.-K. (2022). A positive correlation between GC content and growth temperature in prokaryotes. BMC Genomics 23, 110. doi: 10.1186/s12864-022-08353-7

PubMed Abstract | Crossref Full Text | Google Scholar

Ibe C., Munro C. A. (2021). Fungal cell wall proteins and signaling pathways form a cytoprotective network to combat stresses. J. Fungi 7, 739. doi: 10.3390/jof7090739

Crossref Full Text | Google Scholar

Isaza R. E. A., Diaz-Trujillo C., Dhillon B., Aerts A., Carlier J., Crane C. F., et al. (2016). Combating a Global Threat to a Clonal Crop: Banana Black Sigatoka Pathogen Pseudocercospora Fijiensis (Synonym Mycosphaerella Fijiensis) Genomes Reveal Clues for Disease Control. PloS Genet. 12, e1005876. doi: 10.1371/journal.pgen.1005876

PubMed Abstract | Crossref Full Text | Google Scholar

Kassambara A., Mundt F. (2020) Factoextra: Extract and visualize the results of multivariate data analyses. Available at: https://CRAN.R-project.org/package=factoextra.

Google Scholar

Kolde R. (2019) Pheatmap: Pretty heatmaps. Available at: https://CRAN.R-project.org/package=pheatmap.

Google Scholar

Kolmogorov M., Yuan J., Lin Y., Pevzner P. A. (2019). Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 37, 540–546. doi: 10.1038/s41587-019-0072-8

PubMed Abstract | Crossref Full Text | Google Scholar

Koonin E. V., Aravind L. (2000). The NACHT family – a new group of predicted NTPases implicated in apoptosis and MHC transcription activation. Trends Biochem. Sci. 25, 223–224. doi: 10.1016/S0968-0004(00)01577-2

PubMed Abstract | Crossref Full Text | Google Scholar

Kostadinova N., Tosi S., Spassova B., Angelova M. (2017). Comparison of the oxidative stress response of two Antarctic fungi to different growth temperatures. Polish Polar Res. 38, 393–408. doi: 10.1515/popore-2017-0015

Crossref Full Text | Google Scholar

Krishnan A., Convey P., Gonzalez-Rocha G., Alias S. A. (2016). Production of extracellular hydrolase enzymes by fungi from King George Island. Polar Biol. 39, 65–76. doi: 10.1007/s00300-014-1606-7

Crossref Full Text | Google Scholar

Larochelle M., Drouin S., Robert F., Turcotte B. (2006). Oxidative stress-activated zinc cluster protein stb5 has dual activator/repressor functions required for pentose phosphate pathway regulation and NADPH production. Mol. Cell. Biol. 26, 6690–6701. doi: 10.1128/MCB.02450-05

PubMed Abstract | Crossref Full Text | Google Scholar

Letunic I., Bork P. (2021). Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296. doi: 10.1093/nar/gkab301

PubMed Abstract | Crossref Full Text | Google Scholar

Lenassi M., Gostinčar C., Jackman S., Turk M., Sadowski I., Nislow C., et al. (2013). Whole genome duplication and enrichment of metal cation transporters revealed by de novo genome sequencing of extremely halotolerant black yeast Hortaea werneckii. PLoS One 8 (8), e71328. doi: 10.1371/journal.pone.0071328

PubMed Abstract | Crossref Full Text | Google Scholar

Leung G., Robson G. D., Robinson C. H. (2011). ) Characterization of cold-tolerant fungi from a decomposing High Arctic moss. Soil Biol. Biochem. 43, 1975–1979. doi: 10.1016/j.soilbio.2011.05.003

Crossref Full Text | Google Scholar

Lim J., Thomas T., Cavicchioli R. (2000). Low temperature regulated DEAD-box RNA helicase from the antarctic archaeon, Methanococcoides burtonii. J. Mol. Biol. 297, 553–567. doi: 10.1006/jmbi.2000.3585

PubMed Abstract | Crossref Full Text | Google Scholar

Lin X., Patel S., Litvintseva A. P., Floyd A., Mitchell T. G., Heitman J. (2009). Diploids in the Cryptococcus neoformans Serotype A Population Homozygous for the α Mating Type Originate via Unisexual Mating. PloS Pathog. 5, e1000283. doi: 10.1371/journal.ppat.1000283

PubMed Abstract | Crossref Full Text | Google Scholar

Liu H., Zheng G., Chen Z., Ding X., Wu J., Zhang H., et al. (2023). Psychrophilic yeasts: insights into their adaptability to extremely cold environments. Genes 14, 158. doi: 10.3390/genes14010158

PubMed Abstract | Crossref Full Text | Google Scholar

Lo C.-C., Chain P. S. G. (2014). Rapid evaluation and quality control of next generation sequencing data with FaQCs. BMC Bioinf. 15, 366. doi: 10.1128/MMBR.00015-06

Crossref Full Text | Google Scholar

Lorrain C., Feurtey A., Möller M., Haueisen J., Stukenbrock E. (2021). Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses. G3 Genes|Genomes|Genetics 11, jkab068. doi: 10.1093/g3journal/jkab068

Crossref Full Text | Google Scholar

MacPherson S., Larochelle M., Turcotte B. (2006). A fungal family of transcriptional regulators: the zinc cluster proteins. Microbiol. Mol. Biol. Rev. 70, 583–604. doi: 10.1128/MMBR.00015-06

PubMed Abstract | Crossref Full Text | Google Scholar

Maggi O., Tosi S., Angelova M., Lagostina E., Fabbri A. A., Pecoraro L., et al. (2013). Adaptation of fungi, including yeasts, to cold environments. Plant Biosyst. - Int. J. Dealing all Aspects Plant Biol. 147, 247–258. doi: 10.1080/11263504.2012.753135

Crossref Full Text | Google Scholar

Martin F., Aerts A., Ahrén D., Brun A., Danchin E. G. J., Duchaussoy F., et al. (2008). The genome of Laccaria bicolor provides insights into mycorrhizal symbiosis. Nature 452, 88–92. doi: 10.1038/nature06556

PubMed Abstract | Crossref Full Text | Google Scholar

McDonald T. R., Mueller O., Dietrich F. S., Lutzoni F. (2013). High-throughput genome sequencing of lichenizing fungi to assess gene loss in the ammonium transporter/ammonia permease gene family. BMC Genomics 14, 225. doi: 10.1186/1471-2164-14-225

PubMed Abstract | Crossref Full Text | Google Scholar

Mo Y. (2006). Probing the nature of hydrogen bonds in DNA base pairs. J. Mol. Modeling 12, 665–672. doi: 10.1007/s00894-005-0021-y

Crossref Full Text | Google Scholar

Mosier A. C., Miller C. S., Frischkorn K. R., Ohm R. A., Li Z., LaButti K., et al. (2016). Fungi contribute critical but spatially varying roles in nitrogen and carbon cycling in acid mine drainage. Front. Microbiol. 7. doi: 10.3389/fmicb.2016.00238

Crossref Full Text | Google Scholar

Nagel J. H., Wingfield M. J., Slippers B. (2021). Increased abundance of secreted hydrolytic enzymes and secondary metabolite gene clusters define the genomes of latent plant pathogens in the Botryosphaeriaceae. BMC Genomics 22, 589. doi: 10.1186/s12864-021-07902-w

PubMed Abstract | Crossref Full Text | Google Scholar

Ng K. P., Yew S. M., Chan C. L., Soo-Hoo T. S., Na S. L., Hassan H., et al. (2012). Sequencing of Cladosporium sphaerospermum, a dematiaceous fungus isolated from blood culture. Eukaryot Cell 11, 705–706. doi: 10.1128/EC.00081-12

PubMed Abstract | Crossref Full Text | Google Scholar

Nielsen K. N., Salgado J. F. M., Natsopoulou M. E., Kristensen T., Stajich J. E., De Fine Licht H. H. (2021). Diploidy within a haploid genus of entomopathogenic fungi. Genome Biol. Evol. 13, evab158. doi: 10.1093/gbe/evab158

PubMed Abstract | Crossref Full Text | Google Scholar

Nieuwenhuis B. P. S., James T. Y. (2016). The frequency of sex in fungi. Philos. Trans. R. Soc. B: Biol. Sci. 371, 20150540. doi: 10.1098/rstb.2015.0540

Crossref Full Text | Google Scholar

Ohm R. A., Feau N., Henrissat B., Schoch C. L., Horwitz B. A., Barry K. W., et al. (2012). Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog 8, e1003037. doi: 10.1371/journal.ppat.1003037

PubMed Abstract | Crossref Full Text | Google Scholar

Onofri S., de Vera J.-P., Zucconi L., Selbmann L., Scalzi G., Venkateswaran K. J., et al. (2015). Survival of Antarctic cryptoendolithic fungi in simulated Martian conditions on board the International Space Station. Astrobiology 15, 2015. doi: 10.1089/ast.2015.1324

Crossref Full Text | Google Scholar

Onofri S., Selbmann L., Zucconi L., Pagano S. (2004). Antarctic microfungi as models for exobiology. Planetary Space Sci. 52, 229–237. doi: 10.1016/j.pss.2003.08.019

Crossref Full Text | Google Scholar

Pacheco N., Orellana-Saez M., Pepczynska M., Enrione J., Bassas-Galia M., Borrero-de Acuña J. M., et al. (2019). Exploiting the natural poly(3-hydroxyalkanoates) production capacity of Antarctic Pseudomonas strains: from unique phenotypes to novel biopolymers. J. Ind. Microbiol. Biotechnol. 46, 1139–1153. doi: 10.1007/s10295-019-02186-2

PubMed Abstract | Crossref Full Text | Google Scholar

Paiva D. S., Fernandes L., Portugal A., Trovão J. (2024). First Genome Sequence of the Microcolonial Black Fungus Saxispiralis lemnorum MUM 23.14: Insights into the Unique Genomic Traits of the Aeminiaceae Family. Microorganisms 12, 104. doi: 10.3390/microorganisms12010104

PubMed Abstract | Crossref Full Text | Google Scholar

Paysan-Lafosse T., Blum M., Stein J., Devraj R., Gueye K., Gupta D., et al. (2023). InterPro in 2023: new tools and resources for the functional analysis of protein sequences. Nucleic Acids Res. 51 (D1), D418–D427. doi: 10.1093/nar/gkac1051

PubMed Abstract | Crossref Full Text | Google Scholar

Pepijn W. K., Pellicer J. (2020). Genome size versus genome assemblies: are the genomes truly expanded in polyploid fungal symbionts? Genome Biol. Evol. 12, 2384–2390. doi: 10.1093/gbe/evaa217

PubMed Abstract | Crossref Full Text | Google Scholar

Peter M., Kohler A., Ohm R. A., Kuo A., Krützmann J., Morin E., et al. (2016). Ectomycorrhizal ecology is imprinted in the genome of the dominant symbiotic fungus Cenococcum geophilum. Nat. Commun. 7, 12662. doi: 10.1038/ncomms12662

PubMed Abstract | Crossref Full Text | Google Scholar

Quarantin A., Hadeler B., Kröger C., Schäfer W., Favaron F., Sella L., et al. (2019). Different hydrophobins of fusarium graminearum are involved in hyphal growth, attachment, water-air interface penetration and plant infection. Front. Microbiol. 10, 751. doi: 10.3389/fmicb.2019.00751

PubMed Abstract | Crossref Full Text | Google Scholar

Ranallo-Benavidez T. R., Jaron K. S., Schatz M. C. (2020). GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat. Commun. 11, 1432. doi: 10.1038/s41467-020-14998-3

PubMed Abstract | Crossref Full Text | Google Scholar

R Core Team. (2023). R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria). Available at: https://www.R-project.org/.

Google Scholar

R Core Team. (2024). _R: A language and environment for statistical computing_ (Vienna, Austria: R Foundation for Statistical Computing). Available at: https://www.R-project.org/.

Google Scholar

Rice P., Longden I., Bleasby A. (2000). EMBOSS: the european molecular biology open software suite. Trends Genet. 16, 276–277. doi: 10.1016/S0168-9525(00)02024-2

PubMed Abstract | Crossref Full Text | Google Scholar

Robert-Paganin J., Réty S., Leulliot N. (2015). Regulation of DEAH/RHA helicases by G-patch proteins. BioMed. Res. Int., 931857. doi: 10.1155/2015/931857

PubMed Abstract | Crossref Full Text | Google Scholar

Robinson C. H. (2001). Cold adaptation in Arctic and Antarctic fungi. New Phytol. 151, 341–353. doi: 10.1046/j.1469-8137.2001.00177.x

Crossref Full Text | Google Scholar

Roy J., Kim K., Maddock J. R., Anthony J. G., Woolford J. L. (1995). The final stages of spliceosome maturation require Spp2p that can interact with the DEAH box protein Prp2p and promote step 1 of splicing. RNA 1, 375–390.

PubMed Abstract | Google Scholar

Ruibal C., Gueidan C., Selbmann L., Gorbushina A. A., Crous P. W., Groenewald J. Z., et al. (2009). Phylogeny of rock-inhabiting fungi related to Dothideomycetes. Stud. Mycology 64, 123–133. doi: 10.3114/sim.2009.64.06

Crossref Full Text | Google Scholar

Saupe S., Turcq B., Bégueret J. (1995). A gene responsible for vegetative incompatibility in the fungus Podospora anserina encodes a protein with a GTP-binding motif and Gβ homologous domain. Gene 162, 135–139. doi: 10.1016/0378-1119(95)00272-8

PubMed Abstract | Crossref Full Text | Google Scholar

Selbmann L., Hoogde G. S., Mazzaglia A., Friedmann E., Onofri S. (2005). Fungi at the edge of life: cryptoendolithic black fungi from Antarctic desert. Stud Mycol. 51, 1–32.

Google Scholar

Selbmann L., Zucconi L., Isola D., Onofri S. (2015). Rock black fungi: excellence in the extremes, from the Antarctic to space. Curr. Genet. 61, 335–345. doi: 10.1007/s00294-014-0457-7

PubMed Abstract | Crossref Full Text | Google Scholar

Selbmann L., Benkő Z., Coleine C., de Hoog S., Donati C., Druzhinina I., et al. (2020). Shed light in the daRk lineagES of the fungal tree of life—STRES. Life 10, 362. doi: 10.3390/life10120362

PubMed Abstract | Crossref Full Text | Google Scholar

Stajich J. E. (2017). Fungal genomes and insights into the evolution of the kingdom. Microbiol. Spectr. 5. doi: 10.1128/microbiolspec.FUNK-0055–2016

Crossref Full Text | Google Scholar

Stellwagen N. C., Stellwagen E. (2019). DNA Thermal Stability Depends Solvent Viscosity. J. Phys. Chem. B. 123, 3649–3657. doi: 10.1021/acs.jpcb.9b01217

PubMed Abstract | Crossref Full Text | Google Scholar

Sterflinger K., Lopandic K., Pandey R. V., Blasi B., Kriegner A. (2014). Nothing special in the specialist? draft genome sequence of cryomyces antarcticus, the most extremophilic fungus from antarctica. PloS One 9 (10), e109908. doi: 10.1371/journal.pone.0109908

PubMed Abstract | Crossref Full Text | Google Scholar

Stivaletta N., Barbieri R. (2008). “Endoliths in terrestrial arid environments: implications for astrobiology,” in From fossils to astrobiology: records of life on earth and search for extraterrestrial biosignatures. Eds. Seckbach J., Walsh M. (Springer Netherlands, Dordrecht), 319–333. doi: 10.1007/978-1-4020-8837-7_15

Crossref Full Text | Google Scholar

Stukenbrock E. H., Christiansen F. B., Hansen T. T., Dutheil J. Y., Schierup M. H. (2012). Fusion of two divergent fungal individuals led to the recent emergence of a unique widespread pathogen species. Proc. Natl. Acad. Sci. U.S.A. 109, 10954–10959. doi: 10.1073/pnas.1201403109

PubMed Abstract | Crossref Full Text | Google Scholar

Su Y., Jiang X., Wu W., Wang M., Hamid M. I., Xiang M., et al. (2016). Genomic, transcriptomic, and proteomic analysis provide insights into the cold adaptation mechanism of the obligate psychrophilic fungus mrakia psychrophila. G3: Genes|Genomes|Genetics 6, 3603. doi: 10.1534/g3.116.033308

Crossref Full Text | Google Scholar

Talbot N. J., Ebbole D. J., Hamer J. E. (1993). Identification and characterization of MPG1, a gene involved in pathogenicity from the rice blast fungus Magnaporthe grisea. Plant Cell 5, 1575–1590. 10.1105/tpc.5.11.1575

PubMed Abstract | Google Scholar

Talbot N. J., Rawlins D., Coddington A. (1988). A rapid method for ploidy determination in fungal cells. Curr. Genet. 14, 51–52. doi: 10.1007/BF00405853

Crossref Full Text | Google Scholar

Tanaka T., Terauchi Y., Yoshimi A., Abe K. (2022). Aspergillus hydrophobins: physicochemical properties, biochemical properties, and functions in solid polymer degradation. Microorganisms 10, 1498. doi: 10.3390/microorganisms10081498

PubMed Abstract | Crossref Full Text | Google Scholar

Taschuk F., Cherry S. (2020). DEAD-box helicases: sensors, regulators, and effectors for antiviral defense. Viruses 12, 181. doi: 10.3390/v12020181

PubMed Abstract | Crossref Full Text | Google Scholar

Tesei D., Marzban G., Zakharova K., Isola D., Selbmann L., Sterflinger K. (2012). Alteration of protein patterns in black rock inhabiting fungi as a response to different temperatures. Fungal Biol. 116, 932–940. doi: 10.1016/j.funbio.2012.06.004

PubMed Abstract | Crossref Full Text | Google Scholar

Testa A. C., Oliver R. P., Hane J. K. (2016). OcculterCut: A comprehensive survey of AT-rich regions in fungal genomes. Genome Biol. Evol. 8, 2044–2064. doi: 10.1093/gbe/evw121

PubMed Abstract | Crossref Full Text | Google Scholar

Teixeira M. M., Moreno L. F., Stielow B. J., Muszewska A., Hainaut M., Gonzaga L., et al. (2017). Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, ascomycota). Stud. Mycol 86, 1–28. doi: 10.1016/j.simyco.2017.01.001

PubMed Abstract | Crossref Full Text | Google Scholar

The UniProt Consortium (2023). UniProt: the universal protein knowledgebase in 2023. Nucleic Acids Res. 51, D523–D531. doi: 10.1093/nar/gkac1052

PubMed Abstract | Crossref Full Text | Google Scholar

Touchette D., Altshuler I., Gostinčar C., Zalar P., Raymond-Bouchard I., Zajc J., et al. (2022). Novel Antarctic yeast adapts to cold by switching energy metabolism and increasing small RNA synthesis. ISME J. 16, 221–232. doi: 10.1038/s41396-021-01030-9

PubMed Abstract | Crossref Full Text | Google Scholar

Treindl A. D., Stapley J., Winter D. J., Cox M. P., Leuchtmann A. (2021). Chromosome-level genomes provide insights into genome evolution, organization and size in Epichloe fungi. Genomics 113, 4267–4275. doi: 10.1016/j.ygeno.2021.11.009

PubMed Abstract | Crossref Full Text | Google Scholar

Tsuji G., Kenmochi Y., Takano Y., Sweigard J., Farrall L., Furusawa I., et al. (2000). Novel fungal transcriptional activators, Cmr1p of Colletotrichum lagenarium and Pig1p of Magnaporthe grisea, contain Cys2His2 zinc finger and Zn(II)2Cys6 binuclear cluster DNA-binding motifs and regulate transcription of melanin biosynthesis genes in a developmentally specific manner. Mol. Microbiol. 38, 940–954. doi: 10.1046/j.1365-2958.2000.02181.x

PubMed Abstract | Crossref Full Text | Google Scholar

Usher J. (2019). The mechanisms of mating in pathogenic fungi-A plastic trait. Genes (Basel). 10, 831. doi: 10.3390/genes10100831

PubMed Abstract | Crossref Full Text | Google Scholar

van Wyk S., Harrison C. H., Wingfield B. D., Vos L. D., Merwe N. A., Steenkamp E. T. (2019). The RIPper, a web-based tool for genome-wide quantification of Repeat-Induced Point (RIP) mutations. PeerJ 7, e7447. doi: 10.7717/peerj.7447

PubMed Abstract | Crossref Full Text | Google Scholar

Vik Å., Rine J. (2001). Upc2p and ecm22p, dual regulators of sterol biosynthesis in saccharomyces cerevisiae. Mol. Cell. Biol. 21, 6395–6405. doi: 10.1128/MCB.21.19.6395-6405.2001

PubMed Abstract | Crossref Full Text | Google Scholar

Ward J. J. H. Jr. (1963). Hierarchical grouping to optimize an objective function. J. Am. Stat. Assoc. 58, pp. doi: 10.1080/01621459.1963.10500845

Crossref Full Text | Google Scholar

Wickham H. (2007). Reshaping data with the reshape package. J. Stat. Software 21 (12), 1–20. doi: 10.18637/jss.v021.i12

Crossref Full Text | Google Scholar

Wickham H. (2016). ggplot2: Elegant graphics for data analysis (New York: Springer-Verlag). Available at: https://ggplot2.tidyverse.org.

Google Scholar

Xu X., Chen J., Xu H., Li D. (2014). Role of a major facilitator superfamily transporter in adaptation capacity of Penicillium funiculosum under extreme acidic stress. Fungal Genet. Biol. 69, 75–83. doi: 10.1016/j.fgb.2014.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

Xu H., Xu D., Liu Y. (2024). Molecular biology applications of psychrophilic enzymes: adaptations, advantages, expression, and prospective. Appl. Biochem. Biotechnol. doi: 10.1007/s12010-023-04810-5

PubMed Abstract | Crossref Full Text | Google Scholar

Yin Y., Mao X., Yang J., Chen X., Mao F., Xu Y. (2012). dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 40, W445–W451. doi: 10.1093/nar/gks479

PubMed Abstract | Crossref Full Text | Google Scholar