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

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

The Plot Thickens: Haploid and Triploid-Like Thalli, Hybridization, and Biased Mating Type Ratios in Letharia

The study of the reproductive biology of lichen fungal symbionts has been traditionally challenging due to their complex lifestyles. Against the common belief of haploidy, a recent genomic study found a triploid-like signal in Letharia. Here, we infer the genome organization and reproduction in Letharia by analyzing genomic data from a pure culture and from thalli, and performing a PCR survey of the MAT locus in natural populations. We found that the read count variation in the four Letharia specimens, including the pure culture derived from a single sexual spore of L. lupina, is consistent with haploidy. By contrast, the L. lupina read counts from a thallus' metagenome are triploid-like. Characterization of the mating-type locus revealed a conserved heterothallic configuration across the genus, along with auxiliary genes that we identified. We found that the mating-type distributions are balanced in North America for L. vulpina and L. lupina, suggesting widespread sexual reproduction, but highly skewed in Europe for L. vulpina, consistent with predominant asexuality. Taken together, we propose that Letharia fungi are heterothallic and typically haploid, and provide evidence that triploid-like individuals are hybrids between L. lupina and an unknown Letharia lineage, reconciling classic systematic and genetic studies with recent genomic observations.

Towards a Systems Biology Approach to Understanding the Lichen Symbiosis: Opportunities and Challenges of Implementing Network Modelling

Lichen associations, a classic model for successful and sustainable interactions between micro-organisms, have been studied for many years. However, there are significant gaps in our understanding about how the lichen symbiosis operates at the molecular level. This review addresses opportunities for expanding current knowledge on signalling and metabolic interplays in the lichen symbiosis using the tools and approaches of systems biology, particularly network modelling. The largely unexplored nature of symbiont recognition and metabolic interdependency in lichens could benefit from applying a holistic approach to understand underlying molecular mechanisms and processes. Together with ‘omics’ approaches, the application of signalling and metabolic network modelling could provide predictive means to gain insights into lichen signalling and metabolic pathways. First, we review the major signalling and recognition modalities in the lichen symbioses studied to date, and then describe how modelling signalling networks could enhance our understanding of symbiont recognition, particularly leveraging omics techniques. Next, we highlight the current state of knowledge on lichen metabolism. We also discuss metabolic network modelling as a tool to simulate flux distribution in lichen metabolic pathways and to analyse the co-dependence between symbionts. This is especially important given the growing number of lichen genomes now available and improved computational tools for reconstructing such models. We highlight the benefits and possible bottlenecks for implementing different types of network models as applied to the study of lichens.

Quantitative Assessment of Nucleocytoplasmic Large DNA Virus and Host Interactions Predicted by Co-occurrence Analyses

Nucleocytoplasmic large DNA viruses (NCLDVs) are highly diverse and abundant in marine environments. However, the knowledge of their hosts is limited because only a few NCLDVs have been isolated so far. Taking advantage of the recent large-scale marine metagenomics census, in silico host prediction approaches are expected to fill the gap and further expand our knowledge of virus-host relationships for unknown NCLDVs. In this study, we built co-occurrence networks of NCLDVs and eukaryotic taxa to predict virus-host interactions using Tara Oceans sequencing data. Using the positive likelihood ratio to assess the performance of host prediction for NCLDVs, we benchmarked several co-occurrence approaches and demonstrated an increase in the odds ratio of predicting true positive relationships 4-fold compared to random host predictions. To further refine host predictions from high-dimensional co-occurrence networks, we developed a phylogeny-informed filtering method, Taxon Interaction Mapper, and showed it further improved the prediction performance by 12-fold. Finally, we inferred virophage-NCLDV networks to corroborate that co-occurrence approaches are effective for predicting interacting partners of NCLDVs in marine environments.IMPORTANCE NCLDVs can infect a wide range of eukaryotes, although their life cycle is less dependent on hosts compared to other viruses. However, our understanding of NCLDV-host systems is highly limited because few of these viruses have been isolated so far. Co-occurrence information has been assumed to be useful to predict virus-host interactions. In this study, we quantitatively show the effectiveness of co-occurrence inference for NCLDV host prediction. We also improve the prediction performance with a phylogeny-guided method, which leads to a concise list of candidate host lineages for three NCLDV families. Our results underpin the usage of co-occurrence approaches for the metagenomic exploration of the ecology of this diverse group of viruses.

Advances in Understanding of Desiccation Tolerance of Lichens and Lichen-Forming Algae

Lichens are symbiotic associations (holobionts) established between fungi (mycobionts) and certain groups of cyanobacteria or unicellular green algae (photobionts). This symbiotic association has been essential in the colonization of terrestrial dry habitats. Lichens possess key mechanisms involved in desiccation tolerance (DT) that are constitutively present such as high amounts of polyols, LEA proteins, HSPs, a powerful antioxidant system, thylakoidal oligogalactolipids, etc. This strategy allows them to be always ready to survive drastic changes in their water content. However, several studies indicate that at least some protective mechanisms require a minimal time to be induced, such as the induction of the antioxidant system, the activation of non-photochemical quenching including the de-epoxidation of violaxanthin to zeaxanthin, lipid membrane remodeling, changes in the proportions of polyols, ultrastructural changes, marked polysaccharide remodeling of the cell wall, etc. Although DT in lichens is achieved mainly through constitutive mechanisms, the induction of protection mechanisms might allow them to face desiccation stress in a better condition. The proportion and relevance of constitutive and inducible DT mechanisms seem to be related to the ecology at which lichens are adapted to.

Establishment of Agrobacterium tumefaciens-Mediated Transformation of Cladonia macilenta, a Model Lichen-Forming Fungus

Despite the fascinating biology of lichens, such as the symbiotic association of lichen-forming fungi (mycobiont) with their photosynthetic partners and their ability to grow in harsh habitats, lack of genetic tools manipulating mycobiont has hindered studies on genetic mechanisms underpinning lichen biology. Thus, we established an Agrobacterium tumefaciens-mediated transformation (ATMT) system for genetic transformation of a mycobiont isolated from Cladonia macilenta. A set of combinations of ATMT conditions, such as input biomass of mycobiont, co-cultivation period with Agrobacterium cells, and incubation temperature, were tested to identify an optimized ATMT condition for the C. macilenta mycobiont. As a result, more than 10 days of co-cultivation period and at least 2 mg of input biomass of the mycobiont were recommended for an efficient ATMT, owing to extremely slow growth rate of mycobionts in general. Moreover, we examined T-DNA copy number variation in a total of 180 transformants and found that 88% of the transformants had a single copy T-DNA insertion. To identify precise T-DNA insertion sites that interrupt gene function in C. macilenta, we performed TAIL-PCR analyses for selected transformants. A hypothetical gene encoding ankyrin repeats at its C-terminus was interrupted by T-DNA insertion in a transformant producing dark-brown colored pigment. Although the identification of the pigment awaits further investigation, this proof-of-concept study demonstrated the feasibility of use of ATMT in construction of a random T-DNA insertion mutant library in mycobionts for studying genetic mechanisms behind the lichen symbiosis, stress tolerance, and secondary metabolite biosynthesis.

Predicted Input of Uncultured Fungal Symbionts to a Lichen Symbiosis from Metagenome-Assembled Genomes

Basidiomycete yeasts have recently been reported as stably associated secondary fungal symbionts of many lichens, but their role in the symbiosis remains unknown. Attempts to sequence their genomes have been hampered both by the inability to culture them and their low abundance in the lichen thallus alongside two dominant eukaryotes (an ascomycete fungus and chlorophyte alga). Using the lichen Alectoria sarmentosa, we selectively dissolved the cortex layer in which secondary fungal symbionts are embedded to enrich yeast cell abundance and sequenced DNA from the resulting slurries as well as bulk lichen thallus. In addition to yielding a near-complete genome of the filamentous ascomycete using both methods, metagenomes from cortex slurries yielded a 36- to 84-fold increase in coverage and near-complete genomes for two basidiomycete species, members of the classes Cystobasidiomycetes and Tremellomycetes. The ascomycete possesses the largest gene repertoire of the three. It is enriched in proteases often associated with pathogenicity and harbors the majority of predicted secondary metabolite clusters. The basidiomycete genomes possess ∼35% fewer predicted genes than the ascomycete and have reduced secretomes even compared with close relatives, while exhibiting signs of nutrient limitation and scavenging. Furthermore, both basidiomycetes are enriched in genes coding for enzymes producing secreted acidic polysaccharides, representing a potential contribution to the shared extracellular matrix. All three fungi retain genes involved in dimorphic switching, despite the ascomycete not being known to possess a yeast stage. The basidiomycete genomes are an important new resource for exploration of lifestyle and function in fungal–fungal interactions in lichen symbioses.

Phytochemical Investigation of New Algerian Lichen Species: Physcia Mediterranea Nimis

The present study provides new data concerning the chemical characterisation of Physcia mediterranea Nimis, a rare Mediterranean species belonging to the family Physciaceae. The phytochemical screening was carried out using GC-MS, HPLC-ESI-MS-MS, and NMR techniques. Hot extraction of n-hexane was carried out, followed by separation of the part insoluble in methanol: wax (WA-hex), from the part soluble in methanol (ME-hex). GC-MS analysis of the ME-hex part revealed the presence of methylbenzoic acids such as sparassol and atraric acid and a diterpene with a kaurene skeleton which has never been detected before in lichen species. Out of all the compounds identified by HPLC-ESI-MS-MS, sixteen compounds are common between WA-hex and ME-hex. Most are aliphatic fatty acids, phenolic compounds and depsides. The wax part is characterised by the presence of atranorin, a depside of high biological value. Proton 1H and carbon 13C NMR have confirmed its identification. Atranol, chloroatranol (depsides compound), Ffukinanolide (sesquiterpene lactones), leprolomin (diphenyl ether), muronic acid (triterpenes), and ursolic acid (triterpenes) have also been identified in ME-hex. The results suggested that Physcia mediterranea Nimis is a valuable source of bioactive compounds that could be useful for several applications as functional foods, cosmetics, and pharmaceuticals.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

ViralRecall-A Flexible Command-Line Tool for the Detection of Giant Virus Signatures in 'Omic Data

Giant viruses are widespread in the biosphere and play important roles in biogeochemical cycling and host genome evolution. Also known as nucleo-cytoplasmic large DNA viruses (NCLDVs), these eukaryotic viruses harbor the largest and most complex viral genomes known. Studies have shown that NCLDVs are frequently abundant in metagenomic datasets, and that sequences derived from these viruses can also be found endogenized in diverse eukaryotic genomes. The accurate detection of sequences derived from NCLDVs is therefore of great importance, but this task is challenging owing to both the high level of sequence divergence between NCLDV families and the extraordinarily high diversity of genes encoded in their genomes, including some encoding for metabolic or translation-related functions that are typically found only in cellular lineages. Here, we present ViralRecall, a bioinformatic tool for the identification of NCLDV signatures in ‘omic data. This tool leverages a library of giant virus orthologous groups (GVOGs) to identify sequences that bear signatures of NCLDVs. We demonstrate that this tool can effectively identify NCLDV sequences with high sensitivity and specificity. Moreover, we show that it can be useful both for removing contaminating sequences in metagenome-assembled viral genomes as well as the identification of eukaryotic genomic loci that derived from NCLDV. ViralRecall is written in Python 3.5 and is freely available on GitHub: https://github.com/faylward/viralrecall.

Population genomics of a reindeer lichen species from North American lichen woodlands

PREMISE

Lichens are one of the main structural components of plant communities in the North American boreal biome. They play a pivotal role in lichen woodlands, a large ecosystem situated north of the closed-crown forest zone, and south of the forest-tundra zone. In Eastern Canada (Quebec), there is a remnant LW found 500 km south of its usual distribution range, in the Parc National des Grands-Jardins, originated mainly because of wildfires. We inferred the origin of the lichen Cladonia stellaris from this LW and assessed its genetic diversity in a postfire succession.

METHODS

We genotyped 122 individuals collected across a latitudinal gradient in Quebec. Using the software Stacks, we compared four different approaches of locus selection and single-nucleotide polymorphism calling. We identified the best fitting approach to investigate population structure and estimate genetic diversity of C. stellaris.

RESULTS

Populations in southern Quebec are not genetically different from those of northern LWs. The species consists of at least four phylogenetic lineages with elevated levels of genetic diversity and low co-ancestry. In Parc National des Grands-Jardins, we reported high values of genetic diversity not related with time since fire disturbance and low genetic differentiation among populations with different fire histories.

CONCLUSIONS

This first population genomic study of C. stellaris is an important step forward to understand the origin and biogeographic patterns of lichen woodlands in North America. Our findings also contribute to the understanding of the effect of postfire succession on the genetic structure of the species.

An Overview of Genomics, Phylogenomics and Proteomics Approaches in Ascomycota

Fungi are among the most successful eukaryotes on Earth: they have evolved strategies to survive in the most diverse environments and stressful conditions and have been selected and exploited for multiple aims by humans. The characteristic features intrinsic of Fungi have required evolutionary changes and adaptations at deep molecular levels. Omics approaches, nowadays including genomics, metagenomics, phylogenomics, transcriptomics, metabolomics, and proteomics have enormously advanced the way to understand fungal diversity at diverse taxonomic levels, under changeable conditions and in still under-investigated environments. These approaches can be applied both on environmental communities and on individual organisms, either in nature or in axenic culture and have led the traditional morphology-based fungal systematic to increasingly implement molecular-based approaches. The advent of next-generation sequencing technologies was key to boost advances in fungal genomics and proteomics research. Much effort has also been directed towards the development of methodologies for optimal genomic DNA and protein extraction and separation. To date, the amount of proteomics investigations in Ascomycetes exceeds those carried out in any other fungal group. This is primarily due to the preponderance of their involvement in plant and animal diseases and multiple industrial applications, and therefore the need to understand the biological basis of the infectious process to develop mechanisms for biologic control, as well as to detect key proteins with roles in stress survival. Here we chose to present an overview as much comprehensive as possible of the major advances, mainly of the past decade, in the fields of genomics (including phylogenomics) and proteomics of Ascomycota, focusing particularly on those reporting on opportunistic pathogenic, extremophilic, polyextremotolerant and lichenized fungi. We also present a review of the mostly used genome sequencing technologies and methods for DNA sequence and protein analyses applied so far for fungi.

The Mosaic Architecture of NRPS-PKS in the Arbuscular Mycorrhizal Fungus Gigaspora margarita Shows a Domain With Bacterial Signature

As obligate biotrophic symbionts, arbuscular mycorrhizal fungi (AMF) live in association with most land plants. Among them, Gigaspora margarita has been deeply investigated because of its peculiar features, i.e., the presence of an intracellular microbiota with endobacteria and viruses. The genome sequencing of this fungus revealed the presence of some hybrid non-ribosomal peptide synthases-polyketide synthases (NRPS-PKS) that have been rarely identified in AMF. The aim of this study is to describe the architecture of these NRPS-PKS sequences and to understand whether they are present in other fungal taxa related to G. margarita. A phylogenetic analysis shows that the ketoacyl synthase (KS) domain of one G. margarita NRPS-PKS clusters with prokaryotic sequences. Since horizontal gene transfer (HGT) has often been advocated as a relevant evolutionary mechanism for the spread of secondary metabolite genes, we hypothesized that a similar event could have interested the KS domain of the PKS module. The bacterial endosymbiont of G. margarita, Candidatus Glomeribacter gigasporarum (CaGg), was the first candidate as a donor, since it possesses a large biosynthetic cluster involving an NRPS-PKS. However, bioinformatics analyses do not confirm the hypothesis of a direct HGT from the endobacterium to the fungal host: indeed, endobacterial and fungal sequences show a different evolution and potentially different donors. Lastly, by amplifying a NRPS-PKS conserved fragment and mining the sequenced AMF genomes, we demonstrate that, irrespective of the presence of CaGg, G. margarita, and some other related Gigasporaceae possess such a sequence.

Cell Communications among Microorganisms, Plants, and Animals: Origin, Evolution, and Interplays

Cellular communications play pivotal roles in multi-cellular species, but they do so also in uni-cellular species. Moreover, cells communicate with each other not only within the same individual, but also with cells in other individuals belonging to the same or other species. These communications occur between two unicellular species, two multicellular species, or between unicellular and multicellular species. The molecular mechanisms involved exhibit diversity and specificity, but they share common basic features, which allow common pathways of communication between different species, often phylogenetically very distant. These interactions are possible by the high degree of conservation of the basic molecular mechanisms of interaction of many ligand-receptor pairs in evolutionary remote species. These inter-species cellular communications played crucial roles during Evolution and must have been positively selected, particularly when collectively beneficial in hostile environments. It is likely that communications between cells did not arise after their emergence, but were part of the very nature of the first cells. Synchronization of populations of non-living protocells through chemical communications may have been a mandatory step towards their emergence as populations of living cells and explain the large commonality of cell communication mechanisms among microorganisms, plants, and animals.

In vitro resynthesis of lichenization reveals the genetic background of symbiosis-specific fungal-algal interaction in Usnea hakonensis

BACKGROUND

Symbiosis is central to ecosystems and has been an important driving force of the diversity of life. Close and long-term interactions are known to develop cooperative molecular mechanisms between the symbiotic partners and have often given them new functions as symbiotic entities. In lichen symbiosis, mutualistic relationships between lichen-forming fungi and algae and/or cyanobacteria produce unique features that make lichens adaptive to a wide range of environments. Although the morphological, physiological, and ecological uniqueness of lichens has been described for more than a century, the genetic mechanisms underlying this symbiosis are still poorly known.

RESULTS

This study investigated the fungal-algal interaction specific to the lichen symbiosis using Usnea hakonensis as a model system. The whole genome of U. hakonensis, the fungal partner, was sequenced by using a culture isolated from a natural lichen thallus. Isolated cultures of the fungal and the algal partners were co-cultured in vitro for 3 months, and thalli were successfully resynthesized as visible protrusions. Transcriptomes of resynthesized and natural thalli (symbiotic states) were compared to that of isolated cultures (non-symbiotic state). Sets of fungal and algal genes up-regulated in both symbiotic states were identified as symbiosis-related genes.

CONCLUSION

From predicted functions of these genes, we identified genetic association with two key features fundamental to the symbiotic lifestyle in lichens. The first is establishment of a fungal symbiotic interface: (a) modification of cell walls at fungal-algal contact sites; and (b) production of a hydrophobic layer that ensheaths fungal and algal cells;. The second is symbiosis-specific nutrient flow: (a) the algal supply of photosynthetic product to the fungus; and (b) the fungal supply of phosphorous and nitrogen compounds to the alga. Since both features are widespread among lichens, our result may indicate important facets of the genetic basis of the lichen symbiosis.

Amyloid Signaling in Filamentous Fungi and Bacteria

Amyloids are implicated in many protein misfolding diseases. Amyloid folds, however, also display a range of functional roles particularly in the microbial world. The templating ability of these folds endows them with specific properties allowing their self-propagation and protein-to-protein transmission in vivo. This property, the prion principle, is exploited by specific signaling pathways that use transmission of the amyloid fold as a way to convey information from a receptor to an effector protein. I describe here amyloid signaling pathways involving fungal nucleotide binding and oligomerization domain (NOD)-like receptors that were found to control nonself recognition and programmed cell death processes. Studies on these fungal amyloid signaling motifs stem from the characterization of the fungal [Het-s] prion protein and have led to the identification in fungi but also in multicellular bacteria of several distinct families of signaling motifs, one of which is related to RHIM [receptor-interacting protein (RIP) homotypic interaction motif], an amyloid motif regulating mammalian necroptosis.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus. Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Complete, high-quality genomes from long-read metagenomic sequencing of two wolf lichen thalli reveals enigmatic genome architecture

Fungal genomes display incredible levels of complexity and diversity, and are exceptional study systems for genome evolution. Here we used the Oxford Nanopore MinION sequencing platform to generate high-quality fungal genomes from complex metagenomic samples of lichen thalli. We sequenced two wolf lichens using one flow cell per sample, generating 17.1 Gbps for Letharia lupina and 14.3 Gbps for Letharia columbiana. The resulting L. lupina genome is one of the most contiguous lichen genomes available to date, with 49.2 Mbp contained on 31 contigs. The L. columbiana genome, while less contiguous, is still relatively high quality, with 52.3 Mbp on a total of 161 contigs. Each thallus for both species contained multiple distinct haplotypes, a phenomenon that has rarely been empirically demonstrated. The Oxford Nanopore sequencing technologies are robust and effective when applied to complex symbioses, and have the potential to fundamentally transform our understanding of fungal genetics.

What Is in Umbilicaria pustulata? A Metagenomic Approach to Reconstruct the Holo-Genome of a Lichen

Lichens are valuable models in symbiosis research and promising sources of biosynthetic genes for biotechnological applications. Most lichenized fungi grow slowly, resist aposymbiotic cultivation, and are poor candidates for experimentation. Obtaining contiguous, high-quality genomes for such symbiotic communities is technically challenging. Here, we present the first assembly of a lichen holo-genome from metagenomic whole-genome shotgun data comprising both PacBio long reads and Illumina short reads. The nuclear genomes of the two primary components of the lichen symbiosis-the fungus Umbilicaria pustulata (33 Mb) and the green alga Trebouxia sp. (53 Mb)-were assembled at contiguities comparable to single-species assemblies. The analysis of the read coverage pattern revealed a relative abundance of fungal to algal nuclei of ∼20:1. Gap-free, circular sequences for all organellar genomes were obtained. The bacterial community is dominated by Acidobacteriaceae and encompasses strains closely related to bacteria isolated from other lichens. Gene set analyses showed no evidence of horizontal gene transfer from algae or bacteria into the fungal genome. Our data suggest a lineage-specific loss of a putative gibberellin-20-oxidase in the fungus, a gene fusion in the fungal mitochondrion, and a relocation of an algal chloroplast gene to the algal nucleus. Major technical obstacles during reconstruction of the holo-genome were coverage differences among individual genomes surpassing three orders of magnitude. Moreover, we show that GC-rich inverted repeats paired with nonrandom sequencing error in PacBio data can result in missing gene predictions. This likely poses a general problem for genome assemblies based on long reads.

Lichens and associated fungi from Glacier Bay National Park, Alaska

Lichens are widely acknowledged to be a key component of high latitude ecosystems. However, the time investment needed for full inventories and the lack of taxonomic identification resources for crustose lichen and lichenicolous fungal diversity have hampered efforts to fully gauge the depth of species richness in these ecosystems. Using a combination of classical field inventory and extensive deployment of chemical and molecular analysis, we assessed the diversity of lichens and associated fungi in Glacier Bay National Park, Alaska (USA), a mixed landscape of coastal boreal rainforest and early successional low elevation habitats deglaciated after the Little Ice Age. We collected nearly 5000 specimens and found a total of 947 taxa, including 831 taxa of lichen-forming and 96 taxa of lichenicolous fungi together with 20 taxa of saprotrophic fungi typically included in lichen studies. A total of 98 species (10.3% of those detected) could not be assigned to known species and of those, two genera and 27 species are described here as new to science: Atrophysma cyanomelanos gen. et sp. nov., Bacidina circumpulla, Biatora marmorea, Carneothele sphagnicola gen. et sp. nov., Cirrenalia lichenicola, Corticifraga nephromatis, Fuscidea muskeg, Fuscopannaria dillmaniae, Halecania athallina, Hydropunctaria alaskana, Lambiella aliphatica, Lecania hydrophobica, Lecanora viridipruinosa, Lecidea griseomarginata, L. streveleri, Miriquidica gyrizans, Niesslia peltigerae, Ochrolechia cooperi, Placynthium glaciale, Porpidia seakensis, Rhizocarpon haidense, Sagiolechia phaeospora, Sclerococcum fissurinae, Spilonema maritimum, Thelocarpon immersum, Toensbergia blastidiata and Xenonectriella nephromatis. An additional 71 'known unknown' species are cursorily described. Four new combinations are made: Lepra subvelata (G. K. Merr.) T. Sprib., Ochrolechia minuta (Degel.) T. Sprib., Steineropsis laceratula (Hue) T. Sprib. & Ekman and Toensbergia geminipara (Th. Fr.) T. Sprib. & Resl. Thirty-eight taxa are new to North America and 93 additional taxa new to Alaska. We use four to eight DNA loci to validate the placement of ten of the new species in the orders Baeomycetales, Ostropales, Lecanorales, Peltigerales, Pertusariales and the broader class Lecanoromycetes with maximum likelihood analyses. We present a total of 280 new fungal DNA sequences. The lichen inventory from Glacier Bay National Park represents the second largest number of lichens and associated fungi documented from an area of comparable size and the largest to date in North America. Coming from almost 60°N, these results again underline the potential for high lichen diversity in high latitude ecosystems.

Understanding transport processes in lichen, Azolla-cyanobacteria, ectomycorrhiza, endomycorrhiza, and rhizobia-legume symbiotic interactions

Intimate interactions between photosynthetic and non-photosynthetic organisms require the orchestrated transfer of ions and metabolites between species. We review recent progress in identifying and characterizing the transport proteins involved in five mutualistic symbiotic interactions: lichens, Azolla–cyanobacteria, ectomycorrhiza, endomycorrhiza, and rhizobia–legumes. This review focuses on transporters for nitrogen and carbon and other solutes exchanged in the interactions. Their predicted functions are evaluated on the basis of their transport mechanism and prevailing transmembrane gradients of H ⁺ and transported substrates. The symbiotic interactions are presented in the assumed order from oldest to most recently evolved.

Morphological and Molecular Characterization of Three Endolichenic Isolates of Xylaria (Xylariaceae), from Cladonia curta Ahti & Marcelli (Cladoniaceae)

Endophyte biology is a branch of science that contributes to the understanding of the diversity and ecology of microorganisms that live inside plants, fungi, and lichen. Considering that the diversity of endolichenic fungi is little explored, and its phylogenetic relationship with other lifestyles (endophytism and saprotrophism) is still to be explored in detail, this paper presents data on axenic cultures and phylogenetic relationships of three endolichenic fungi, isolated in laboratory. Cladonia curta Ahti & Marcelli, a species of lichen described in Brazil, is distributed at three sites in the Southeast of the country, in mesophilous forests and the Cerrado. Initial hyphal growth of Xylaria spp. on C. curta podetia started four days after inoculation and continued for the next 13 days until the hyphae completely covered the podetia. Stromata formation and differentiation was observed, occurring approximately after one year of isolation and consecutive subculture of lineages. Phylogenetic analyses indicate lineages of endolichenic fungi in the genus Xylaria, even as the morphological characteristics of the colonies and anamorphous stromata confirm this classification. Our preliminary results provide evidence that these endolichenic fungi are closely related to endophytic fungi, suggesting that the associations are not purely incidental. Further studies, especially phylogenetic analyses using robust multi-locus datasets, are needed to accept or reject the hypothesis that endolichenic fungi isolated from Xylaria spp. and X. berteri are conspecific.