Multiomics in the central Arctic Ocean for benchmarking biodiversity change.

Multiomics approaches need to be applied in the central Arctic Ocean to benchmark biodiversity change and to identify novel species and their genes. As part of MOSAiC, EcoOmics will therefore be essential for conservation and sustainable bioprospecting in one of the least explored ecosystems on Earth.

The genome of Naegleria gruberi illuminates early eukaryotic versatility.

Genome sequences of diverse free-living protists are essential for understanding eukaryotic evolution and molecular and cell biology. The free-living amoeboflagellate Naegleria gruberi belongs to a varied and ubiquitous protist clade (Heterolobosea) that diverged from other eukaryotic lineages over a billion years ago. Analysis of the 15,727 protein-coding genes encoded by Naegleria's 41 Mb nuclear genome indicates a capacity for both aerobic respiration and anaerobic metabolism with concomitant hydrogen production, with fundamental implications for the evolution of organelle metabolism. The Naegleria genome facilitates substantially broader phylogenomic comparisons of free-living eukaryotes than previously possible, allowing us to identify thousands of genes likely present in the pan-eukaryotic ancestor, with 40% likely eukaryotic inventions. Moreover, we construct a comprehensive catalog of amoeboid-motility genes. The Naegleria genome, analyzed in the context of other protists, reveals a remarkably complex ancestral eukaryote with a rich repertoire of cytoskeletal, sexual, signaling, and metabolic modules.

Diploid-dominant life cycles characterize the early evolution of Fungi.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Speciation Underpinned by Unexpected Molecular Diversity in the Mycorrhizal Fungal Genus Pisolithus.

The mutualistic ectomycorrhizal (ECM) fungal genus Pisolithus comprises 19 species defined to date which colonize the roots of >50 hosts worldwide suggesting that substantial genomic and functional evolution occurred during speciation. To better understand this intra-genus variation, we undertook a comparative multi-omic study of nine Pisolithus species sampled from North America, South America, Asia, and Australasia. We found that there was a small core set of genes common to all species (13%), and that these genes were more likely to be significantly regulated during symbiosis with a host than accessory or species-specific genes. Thus, the genetic "toolbox" foundational to the symbiotic lifestyle in this genus is small. Transposable elements were located significantly closer to gene classes including effector-like small secreted proteins (SSPs). Poorly conserved SSPs were more likely to be induced by symbiosis, suggesting that they may be a class of protein that tune host specificity. The Pisolithus gene repertoire is characterized by divergent CAZyme profiles when compared with other fungi, both symbiotic and saprotrophic. This was driven by differences in enzymes associated with symbiotic sugar processing, although metabolomic analysis suggest that neither copy number nor expression of these genes is sufficient to predict sugar capture from a host plant or its metabolism in fungal hyphae. Our results demonstrate that intra-genus genomic and functional diversity within ECM fungi is greater than previously thought, underlining the importance of continued comparative studies within the fungal tree of life to refine our focus on pathways and evolutionary processes foundational to this symbiotic lifestyle.

A scaffolded and annotated reference genome of giant kelp (Macrocystis pyrifera).

Macrocystis pyrifera (giant kelp), is a brown macroalga of great ecological importance as a primary producer and structure-forming foundational species that provides habitat for hundreds of species. It has many commercial uses (e.g. source of alginate, fertilizer, cosmetics, feedstock). One of the limitations to exploiting giant kelp's economic potential and assisting in giant kelp conservation efforts is a lack of genomic tools like a high quality, contiguous reference genome with accurate gene annotations. Reference genomes attempt to capture the complete genomic sequence of an individual or species, and importantly provide a universal structure for comparison across a multitude of genetic experiments, both within and between species. We assembled the giant kelp genome of a haploid female gametophyte de novo using PacBio reads, then ordered contigs into chromosome level scaffolds using Hi-C. We found the giant kelp genome to be 537 MB, with a total of 35 scaffolds and 188 contigs. The assembly N50 is 13,669,674 with GC content of 50.37%. We assessed the genome completeness using BUSCO, and found giant kelp contained 94% of the BUSCO genes from the stramenopile clade. Annotation of the giant kelp genome revealed 25,919 genes. Additionally, we present genetic variation data based on 48 diploid giant kelp sporophytes from three different Southern California populations that confirms the population structure found in other studies of these populations. This work resulted in a high-quality giant kelp genome that greatly increases the genetic knowledge of this ecologically and economically vital species.

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales.

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

PhycoCosm, a comparative algal genomics resource.

Algae are a diverse, polyphyletic group of photosynthetic eukaryotes spanning nearly all eukaryotic lineages of life and collectively responsible for ∼50% of photosynthesis on Earth. Sequenced algal genomes, critical to understanding their complex biology, are growing in number and require efficient tools for analysis. PhycoCosm (https://phycocosm.jgi.doe.gov) is an algal multi-omics portal, developed by the US Department of Energy Joint Genome Institute to support analysis and distribution of algal genome sequences and other 'omics' data. PhycoCosm provides integration of genome sequence and annotation for >100 algal genomes with available multi-omics data and interactive web-based tools to enable algal research in bioenergy and the environment, encouraging community engagement and data exchange, and fostering new sequencing projects that will further these research goals.

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity.

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates-namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose-methanol-choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases-we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms.

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales.

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

Ecological generalism drives hyperdiversity of secondary metabolite gene clusters in xylarialean endophytes.

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Comparative genomics reveals dynamic genome evolution in host specialist ectomycorrhizal fungi.

While there has been significant progress characterizing the 'symbiotic toolkit' of ectomycorrhizal (ECM) fungi, how host specificity may be encoded into ECM fungal genomes remains poorly understood. We conducted a comparative genomic analysis of ECM fungal host specialists and generalists, focusing on the specialist genus Suillus. Global analyses of genome dynamics across 46 species were assessed, along with targeted analyses of three classes of molecules previously identified as important determinants of host specificity: small secreted proteins (SSPs), secondary metabolites (SMs) and G-protein coupled receptors (GPCRs). Relative to other ECM fungi, including other host specialists, Suillus had highly dynamic genomes including numerous rapidly evolving gene families and many domain expansions and contractions. Targeted analyses supported a role for SMs but not SSPs or GPCRs in Suillus host specificity. Phylogenomic-based ancestral state reconstruction identified Larix as the ancestral host of Suillus, with multiple independent switches between white and red pine hosts. These results suggest that like other defining characteristics of the ECM lifestyle, host specificity is a dynamic process at the genome level. In the case of Suillus, both SMs and pathways involved in the deactivation of reactive oxygen species appear to be strongly associated with enhanced host specificity.

Desert truffle genomes reveal their reproductive modes and new insights into plant-fungal interaction and ectendomycorrhizal lifestyle.

Desert truffles are edible hypogeous fungi forming ectendomycorrhizal symbiosis with plants of Cistaceae family. Knowledge about the reproductive modes of these fungi and the molecular mechanisms driving the ectendomycorrhizal interaction is lacking. Genomes of the highly appreciated edible desert truffles Terfezia claveryi Chatin and Tirmania nivea Trappe have been sequenced and compared with other Pezizomycetes. Transcriptomes of T. claveryi × Helianthemum almeriense mycorrhiza from well-watered and drought-stressed plants, when intracellular colonizations is promoted, were investigated. We have identified the fungal genes related to sexual reproduction in desert truffles and desert-truffles-specific genomic and secretomic features with respect to other Pezizomycetes, such as the expansion of a large set of gene families with unknown Pfam domains and a number of species or desert-truffle-specific small secreted proteins differentially regulated in symbiosis. A core set of plant genes, including carbohydrate, lipid-metabolism, and defence-related genes, differentially expressed in mycorrhiza under both conditions was found. Our results highlight the singularities of desert truffles with respect to other mycorrhizal fungi while providing a first glimpse on plant and fungal determinants involved in ecto to endo symbiotic switch that occurs in desert truffle under dry conditions.

Massive lateral transfer of genes encoding plant cell wall-degrading enzymes to the mycoparasitic fungus Trichoderma from its plant-associated hosts.

Unlike most other fungi, molds of the genus Trichoderma (Hypocreales, Ascomycota) are aggressive parasites of other fungi and efficient decomposers of plant biomass. Although nutritional shifts are common among hypocrealean fungi, there are no examples of such broad substrate versatility as that observed in Trichoderma. A phylogenomic analysis of 23 hypocrealean fungi (including nine Trichoderma spp. and the related Escovopsis weberi) revealed that the genus Trichoderma has evolved from an ancestor with limited cellulolytic capability that fed on either fungi or arthropods. The evolutionary analysis of Trichoderma genes encoding plant cell wall-degrading carbohydrate-active enzymes and auxiliary proteins (pcwdCAZome, 122 gene families) based on a gene tree / species tree reconciliation demonstrated that the formation of the genus was accompanied by an unprecedented extent of lateral gene transfer (LGT). Nearly one-half of the genes in Trichoderma pcwdCAZome (41%) were obtained via LGT from plant-associated filamentous fungi belonging to different classes of Ascomycota, while no LGT was observed from other potential donors. In addition to the ability to feed on unrelated fungi (such as Basidiomycota), we also showed that Trichoderma is capable of endoparasitism on a broad range of Ascomycota, including extant LGT donors. This phenomenon was not observed in E. weberi and rarely in other mycoparasitic hypocrealean fungi. Thus, our study suggests that LGT is linked to the ability of Trichoderma to parasitize taxonomically related fungi (up to adelphoparasitism in strict sense). This may have allowed primarily mycotrophic Trichoderma fungi to evolve into decomposers of plant biomass.

Linking secondary metabolites to gene clusters through genome sequencing of six diverse Aspergillus species.

The fungal genus of Aspergillus is highly interesting, containing everything from industrial cell factories, model organisms, and human pathogens. In particular, this group has a prolific production of bioactive secondary metabolites (SMs). In this work, four diverse Aspergillus species (A. campestris, A. novofumigatus, A. ochraceoroseus, and A. steynii) have been whole-genome PacBio sequenced to provide genetic references in three Aspergillus sections. A. taichungensis and A. candidus also were sequenced for SM elucidation. Thirteen Aspergillus genomes were analyzed with comparative genomics to determine phylogeny and genetic diversity, showing that each presented genome contains 15-27% genes not found in other sequenced Aspergilli. In particular, A. novofumigatus was compared with the pathogenic species A. fumigatus This suggests that A. novofumigatus can produce most of the same allergens, virulence, and pathogenicity factors as A. fumigatus, suggesting that A. novofumigatus could be as pathogenic as A. fumigatus Furthermore, SMs were linked to gene clusters based on biological and chemical knowledge and analysis, genome sequences, and predictive algorithms. We thus identify putative SM clusters for aflatoxin, chlorflavonin, and ochrindol in A. ochraceoroseus, A. campestris, and A. steynii, respectively, and novofumigatonin, ent-cycloechinulin, and epi-aszonalenins in A. novofumigatus Our study delivers six fungal genomes, showing the large diversity found in the Aspergillus genus; highlights the potential for discovery of beneficial or harmful SMs; and supports reports of A. novofumigatus pathogenicity. It also shows how biological, biochemical, and genomic information can be combined to identify genes involved in the biosynthesis of specific SMs.

Fungal ecological strategies reflected in gene transcription - a case study of two litter decomposers.

Microbial communities interplay with their environment through their functional traits that can be a response or an effect on the environment. Here, we explore how a functional trait-the decomposition of organic matter, can be addressed based on genetic markers and how the expression of these markers reflect ecological strategies of two fungal litter decomposer Gymnopus androsaceus and Chalara longipes. We sequenced the genomes of these two fungi, as well as their transcriptomes at different steps of Pinus sylvestris needles decomposition in microcosms. Our results highlighted that if the gene content of the two species could indicate similar potential decomposition abilities, the expression levels of specific gene families belonging to the glycoside hydrolase category reflected contrasting ecological strategies. Actually, C. longipes, the weaker decomposer in this experiment, turned out to have a high content of genes involved in cell wall polysaccharides decomposition but low expression levels, reflecting a versatile ecology compare to the more competitive G. androsaceus with high expression levels of keystone functional genes. Thus, we established that sequential expression of genes coding for different components of the decomposer machinery indicated adaptation to chemical changes in the substrate as decomposition progressed.

Evolution and comparative genomics of the most common Trichoderma species.

BACKGROUND: The growing importance of the ubiquitous fungal genus Trichoderma (Hypocreales, Ascomycota) requires understanding of its biology and evolution. Many Trichoderma species are used as biofertilizers and biofungicides and T. reesei is the model organism for industrial production of cellulolytic enzymes. In addition, some highly opportunistic species devastate mushroom farms and can become pathogens of humans. A comparative analysis of the first three whole genomes revealed mycoparasitism as the innate feature of Trichoderma. However, the evolution of these traits is not yet understood. RESULTS: We selected 12 most commonly occurring Trichoderma species and studied the evolution of their genome sequences. Trichoderma evolved in the time of the Cretaceous-Palaeogene extinction event 66 (±15) mya, but the formation of extant sections (Longibrachiatum, Trichoderma) or clades (Harzianum/Virens) happened in Oligocene. The evolution of the Harzianum clade and section Trichoderma was accompanied by significant gene gain, but the ancestor of section Longibrachiatum experienced rapid gene loss. The highest number of genes gained encoded ankyrins, HET domain proteins and transcription factors. We also identified the Trichoderma core genome, completely curated its annotation, investigated several gene families in detail and compared the results to those of other fungi. Eighty percent of those genes for which a function could be predicted were also found in other fungi, but only 67% of those without a predictable function. CONCLUSIONS: Our study presents a time scaled pattern of genome evolution in 12 Trichoderma species from three phylogenetically distant clades/sections and a comprehensive analysis of their genes. The data offer insights in the evolution of a mycoparasite towards a generalist.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata.

BACKGROUND: Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. RESULTS: A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. CONCLUSIONS: The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.