Sequencing genomes from mixed DNA samples - evaluating the metagenome skimming approach in lichenized fungi

Symbionts out of sync: Decoupled physiological responses are widespread and ecologically important in lichen associations

A core vulnerability in symbioses is the need for coordination between the symbiotic partners, which are often assumed to be closely physiologically integrated. We critically re-examine this assumed integration between symbionts in lichen symbioses, recovering a long overlooked yet fundamental physiological asymmetry in carbon balance. We examine the physiological, ecological, and transcriptional basis of this asymmetry in the lichen Evernia mesomorpha. This carbon balance asymmetry depends on hydration source and aligns with climatic range limits. Differences in gene expression across the E. mesomorpha symbiosis suggest that the physiologies of the primary lichen symbionts are decoupled. Furthermore, we use gas exchange data to show that asymmetries in carbon balance are widespread and common across evolutionarily disparate lichen associations. Using carbon balance asymmetry as an example, we provide evidence for the wide-ranging importance of physiological asymmetries in symbioses.

The mitogenomes characterization of two Peltigera species (Peltigera elisabethae and Peltigera polydactylon) and comparative mitogenomic analyses of six Peltigera

In the present study, the mitochondrial genomes of Peltigera elisabethae and P. polydactylon were sequenced and assembled. The two mitogenomes were composed of circular DNA molecules, with sizes of 64,034 bp and 59,208 bp, respectively. Comparative analysis showed that the genome size, GC content, GC skew, and AT skew varied between the two mitochondrial genomes. In codon analysis, phenylalanine (Phe), isoleucine (Ile), and leucine (Leu) were most frequently used in six Peltigera genomes. Evolutionary analysis showed that all 14 protein-coding genes (PCGs) were subject to purifying selection in the six Peltigera species. Regarding gene rearrangement, the PCGs of Peltigera had the same gene sequence and gene content, and a few intron sequences and spacer sequences were rearranged in Peltigera. In the phylogenetic analysis, we used Bayesian and ML methods to construct a phylogenetic tree. Two phylogenetic trees with consistent topology with high support indicate that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships. The present study enriches the mitochondrial genome data of Peltigera and promotes further understanding of the genetics and evolution of the Peltigera genus.

DL-TODA: A Deep Learning Tool for Omics Data Analysis

Metagenomics is a technique for genome-wide profiling of microbiomes; this technique generates billions of DNA sequences called reads. Given the multiplication of metagenomic projects, computational tools are necessary to enable the efficient and accurate classification of metagenomic reads without needing to construct a reference database. The program DL-TODA presented here aims to classify metagenomic reads using a deep learning model trained on over 3000 bacterial species. A convolutional neural network architecture originally designed for computer vision was applied for the modeling of species-specific features. Using synthetic testing data simulated with 2454 genomes from 639 species, DL-TODA was shown to classify nearly 75% of the reads with high confidence. The classification accuracy of DL-TODA was over 0.98 at taxonomic ranks above the genus level, making it comparable with Kraken2 and Centrifuge, two state-of-the-art taxonomic classification tools. DL-TODA also achieved an accuracy of 0.97 at the species level, which is higher than 0.93 by Kraken2 and 0.85 by Centrifuge on the same test set. Application of DL-TODA to the human oral and cropland soil metagenomes further demonstrated its use in analyzing microbiomes from diverse environments. Compared to Centrifuge and Kraken2, DL-TODA predicted distinct relative abundance rankings and is less biased toward a single taxon.

Linking Lichen Metabolites to Genes: Emerging Concepts and Lessons from Molecular Biology and Metagenomics

Lichen secondary metabolites have tremendous pharmaceutical and industrial potential. Although more than 1000 metabolites have been reported from lichens, less than 10 have been linked to the genes coding them. The current biosynthetic research focuses strongly on linking molecules to genes as this is fundamental to adapting the molecule for industrial application. Metagenomic-based gene discovery, which bypasses the challenges associated with culturing an organism, is a promising way forward to link secondary metabolites to genes in non-model, difficult-to-culture organisms. This approach is based on the amalgamation of the knowledge of the evolutionary relationships of the biosynthetic genes, the structure of the target molecule, and the biosynthetic machinery required for its synthesis. So far, metagenomic-based gene discovery is the predominant approach by which lichen metabolites have been linked to their genes. Although the structures of most of the lichen secondary metabolites are well-documented, a comprehensive review of the metabolites linked to their genes, strategies implemented to establish this link, and crucial takeaways from these studies is not available. In this review, I address the following knowledge gaps and, additionally, provide critical insights into the results of these studies, elaborating on the direct and serendipitous lessons that we have learned from them.

Metagenomics Shines Light on the Evolution of "Sunscreen" Pigment Metabolism in the Teloschistales (Lichen-Forming Ascomycota)

Fungi produce a vast number of secondary metabolites that shape their interactions with other organisms and the environment. Characterizing the genes underpinning metabolite synthesis is therefore key to understanding fungal evolution and adaptation. Lichenized fungi represent almost one-third of Ascomycota diversity and boast impressive secondary metabolites repertoires. However, most lichen biosynthetic genes have not been linked to their metabolite products. Here we used metagenomic sequencing to survey gene families associated with production of anthraquinones, UV-protectant secondary metabolites present in various fungi, but especially abundant in a diverse order of lichens, the Teloschistales (class Lecanoromycetes, phylum Ascomycota). We successfully assembled 24 new, high-quality lichenized-fungal genomes de novo and combined them with publicly available Lecanoromycetes genomes from taxa with diverse secondary chemistry to produce a whole-genome tree. Secondary metabolite biosynthetic gene cluster (BGC) analysis showed that whilst lichen BGCs are numerous and highly dissimilar, core enzyme genes are generally conserved across taxa. This suggests metabolite diversification occurs via re-shuffling existing enzyme genes with novel accessory genes rather than BGC gains/losses or de novo gene evolution. We identified putative anthraquinone BGCs in our lichen dataset that appear homologous to anthraquinone clusters from non-lichenized fungi, suggesting these genes were present in the common ancestor of the subphylum Pezizomycotina. Finally, we identified unique transporter genes in Teloschistales anthraquinone BGCs that may explain why these metabolites are so abundant and ubiquitous in these lichens. Our results support the importance of metagenomics for understanding the secondary metabolism of non-model fungi such as lichens.

Recent advances in research for potential utilization of unexplored lichen metabolites

Several research studies have shown that lichens are productive organisms for the synthesis of a broad range of secondary metabolites. Lichens are a self-sustainable stable microbial ecosystem comprising an exhabitant fungal partner (mycobiont) and at least one or more photosynthetic partners (photobiont). The successful symbiosis is responsible for their persistence throughout time and allows all the partners (holobionts) to thrive in many extreme habitats, where without the synergistic relationship they would be rare or non-existent. The ability to survive in harsh conditions can be directly correlated with the production of some unique metabolites. Despite the potential applications, these unique metabolites have been underutilised by pharmaceutical and agrochemical industries due to their slow growth, low biomass availability and technical challenges involved in their artificial cultivation. However, recent development of biotechnological tools such as molecular phylogenetics, modern tissue culture techniques, metabolomics and molecular engineering are opening up a new opportunity to exploit these compounds within the lichen holobiome for industrial applications. This review also highlights the recent advances in culturing the symbionts and the computational and molecular genetics approaches of lichen gene regulation recognized for the enhanced production of target metabolites. The recent development of multi-omics novel biodiscovery strategies aided by synthetic biology in order to study the heterologous expressed lichen-derived biosynthetic gene clusters in a cultivatable host offers a promising means for a sustainable supply of specialized metabolites.

Tissue Cultivation, Preparation, and Extraction of High Molecular Weight DNA for Single-Molecule Genome Sequencing of Plant-Associated Fungi

Extraction of high-quality, high molecular weight DNA is a critical step for sequencing an organism's genome. For fungi, DNA extraction is often complicated by co-precipitation of secondary metabolites, the most destructive being polysaccharides, polyphenols, and melanin. Different DNA extraction protocols and clean-up methods have been developed to address challenging materials and contaminants; however, the method of fungal cultivation and tissue preparation also plays a critical role to limit the production of inhibitory compounds prior to extraction. Here, we provide protocols and guidelines for (i) fungal tissue cultivation and processing with solid media containing a cellophane overlay or in liquid media, (ii) DNA extraction with customized recommendations for taxonomically and ecologically diverse plant-associated fungi, and (iii) assessing DNA quantity and quality for downstream genome sequencing with single-molecule technology such as PacBio.

Genome mining as a biotechnological tool for the discovery of novel biosynthetic genes in lichens

Natural products (NPs) and their derivatives are a major contributor to modern medicine. Historically, microorganisms such as bacteria and fungi have been instrumental in generating drugs and lead compounds because of the ease of culturing and genetically manipulating them. However, the ever-increasing demand for novel drugs highlights the need to bioprospect previously unexplored taxa for their biosynthetic potential. Next-generation sequencing technologies have expanded the range of organisms that can be explored for their biosynthetic content, as these technologies can provide a glimpse of an organism's entire biosynthetic landscape, without the need for cultivation. The entirety of biosynthetic genes can be compared to the genes of known function to identify the gene clusters potentially coding for novel products. In this study, we mine the genomes of nine lichen-forming fungal species of the genus Umbilicaria for biosynthetic genes, and categorize the biosynthetic gene clusters (BGCs) as "associated product structurally known" or "associated product putatively novel". Although lichen-forming fungi have been suggested to be a rich source of NPs, it is not known how their biosynthetic diversity compares to that of bacteria and non-lichenized fungi. We found that 25%-30% of biosynthetic genes are divergent as compared to the global database of BGCs, which comprises 1,200,000 characterized biosynthetic genes from plants, bacteria, and fungi. Out of 217 BGCs, 43 were highly divergant suggesting that they potentially encode structurally and functionally novel NPs. Clusters encoding the putatively novel metabolic diversity comprise polyketide synthases (30), non-ribosomal peptide synthetases (12), and terpenes (1). Our study emphasizes the utility of genomic data in bioprospecting microorganisms for their biosynthetic potential and in advancing the industrial application of unexplored taxa. We highlight the untapped structural metabolic diversity encoded in the lichenized fungal genomes. To the best of our knowledge, this is the first investigation identifying genes coding for NPs with potentially novel properties in lichenized fungi.

Physodic acid sensitizes LNCaP prostate cancer cells to TRAIL-induced apoptosis

In spite of the extensive research for developing new therapies, prostate cancer is still one of the major human diseases with poor prognosis and high mortality. Therefore, with the aim of identifying novel agents with antigrowth and pro-apoptotic activity on prostate cancer cells, in the present study, we evaluated the effect of lichen secondary metabolite physodic acid on cell growth in human prostate cancer cells. In addition, we tested the apoptotic activity of physodic acid on TRAIL-resistant LNCaP cells in combination with TRAIL. The cell viability was measured using MTT assay. LDH release, a marker of membrane breakdown, was also measured. For the detection of apoptosis, the evaluation of DNA fragmentation and caspase-3 activity assay were employed. The expression of proteins was detected by Western blot analysis. It was observed that physodic acid showed a dose-response relationship in the range of 12.5-50 μM concentrations in LNCaP and DU-145 cells, activating an apoptotic process. In addition, physodic acid sensitizes LNCaP cells to TRAIL-induced apoptosis. The combination of physodic acid with other anti-prostate cancer therapies could be considered a promising strategy that warrants further investigations.

Prediction of Genes That Function in Methanogenesis and CO2 Pathways in Extremophiles

Gaet’ale (GAL) and Mud’ara (MUP) are two hypersaline ponds located in the Danakil Depression recharged by underground water from the surrounding highlands. These two ponds have different pH, salinity, and show variation in the concentration of many ionic components. Metagenomic analysis concludes that GAL is dominated by bacteria as in the case of the other hypersaline and acidic ponds in the Danakil Depression. However, Archaea dominated the ponds of MUP. In the current study, the application of SEED and KEGG helped to map the ordered steps of specific enzyme catalyzed reaction in converting CO₂ into cell products. We predict that highly efficient and light-independent carbon fixation involving phosphoenolpyruvate carboxylase takes place in MUP. On the contrary, genes encoding enzymes involved in hydrogenotrophic and acetoclastic methanogenesis appeared solely in ponds of GAL, implying the biological source of the hazardous methane gas in that environment. Based on the investigation of the sources of the genes of interest, it is clear that cooperative interactions between members of the two communities and syntrophic metabolism is the main strategy adapted to utilize inorganic carbon as a carbon source in both MUP and GAL. This insight can be used to design biotechnological applications of microbial communities in production of methane biogas or to minimize CO₂ emissions.

Depside and Depsidone Synthesis in Lichenized Fungi Comes into Focus through a Genome-Wide Comparison of the Olivetoric Acid and Physodic Acid Chemotypes of Pseudevernia furfuracea

Primary biosynthetic enzymes involved in the synthesis of lichen polyphenolic compounds depsides and depsidones are non-reducing polyketide synthases (NR-PKSs), and cytochrome P450s. However, for most depsides and depsidones the corresponding PKSs are unknown. Additionally, in non-lichenized fungi specific fatty acid synthases (FASs) provide starters to the PKSs. Yet, the presence of such FASs in lichenized fungi remains to be investigated. Here we implement comparative genomics and metatranscriptomics to identify the most likely PKS and FASs for olivetoric acid and physodic acid biosynthesis, the primary depside and depsidone defining the two chemotypes of the lichen Pseudevernia furfuracea. We propose that the gene cluster PF33-1_006185, found in both chemotypes, is the most likely candidate for the olivetoric acid and physodic acid biosynthesis. This is the first study to identify the gene cluster and the FAS likely responsible for olivetoric acid and physodic acid biosynthesis in a lichenized fungus. Our findings suggest that gene regulation and other epigenetic factors determine whether the mycobiont produces the depside or the depsidone, providing the first direct indication that chemotype diversity in lichens can arise through regulatory and not only through genetic diversity. Combining these results and existing literature, we propose a detailed scheme for depside/depsidone synthesis.

Microbial Communities of Cladonia Lichens and Their Biosynthetic Gene Clusters Potentially Encoding Natural Products

Lichens have been widely used in traditional medicine, especially by indigenous communities worldwide. However, their slow growth and difficulties in the isolation of lichen symbionts and associated microbes have hindered the pharmaceutical utilisation of lichen-produced compounds. Advances in high-throughput sequencing techniques now permit detailed investigations of the complex microbial communities formed by fungi, green algae, cyanobacteria, and other bacteria within the lichen thalli. Here, we used amplicon sequencing, shotgun metagenomics, and in silico metabolomics together with compound extractions to study reindeer lichens collected from Southern Finland. Our aim was to evaluate the potential of Cladonia species as sources of novel natural products. We compared the predicted biosynthetic pathways of lichen compounds from isolated genome-sequenced lichen fungi and our environmental samples. Potential biosynthetic genes could then be further used to produce secondary metabolites in more tractable hosts. Furthermore, we detected multiple compounds by metabolite analyses, which revealed connections between the identified biosynthetic gene clusters and their products. Taken together, our results contribute to metagenomic data studies from complex lichen-symbiotic communities and provide valuable new information for use in further biochemical and pharmacological studies.

Lichens as a repository of bioactive compounds: an open window for green therapy against diverse cancers

Lichens, algae and fungi-based symbiotic associations, are sources of many important secondary metabolites, such as antibiotics, anti-inflammatory, antioxidants, and anticancer agents. Wide range of experiments based on in vivo and in vitro studies revealed that lichens are a rich treasure of anti-cancer compounds. Lichen extracts and isolated lichen compounds can interact with all biological entities currently identified to be responsible for tumor development. The critical ways to control the cancer development include induction of cell cycle arrests, blocking communication of growth factors, activation of anti-tumor immunity, inhibition of tumor-friendly inflammation, inhibition of tumor metastasis, and suppressing chromosome dysfunction. Also, lichen-based compounds induce the killing of cells by the process of apoptosis, autophagy, and necrosis, that inturn positively modulates metabolic networks of cells against uncontrolled cell division. Many lichen-based compounds have proven to possess potential anti-cancer activity against a wide range of cancer cells, either alone or in conjunction with other anti-cancer compounds. This review primarily emphasizes on an updated account of the repository of secondary metabolites reported in lichens. Besides, we discuss the anti-cancer potential and possible mechanism of the most frequently reported secondary metabolites derived from lichens.

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.

Using target enrichment sequencing to study the higher-level phylogeny of the largest lichen-forming fungi family: Parmeliaceae (Ascomycota)

Parmeliaceae is the largest family of lichen-forming fungi with a worldwide distribution. We used a target enrichment data set and a qualitative selection method for 250 out of 350 genes to infer the phylogeny of the major clades in this family including 81 taxa, with both subfamilies and all seven major clades previously recognized in the subfamily Parmelioideae. The reduced genome-scale data set was analyzed using concatenated-based Bayesian inference and two different Maximum Likelihood analyses, and a coalescent-based species tree method. The resulting topology was strongly supported with the majority of nodes being fully supported in all three concatenated-based analyses. The two subfamilies and each of the seven major clades in Parmelioideae were strongly supported as monophyletic. In addition, most backbone relationships in the topology were recovered with high nodal support. The genus Parmotrema was found to be polyphyletic and consequently, it is suggested to accept the genus Crespoa to accommodate the species previously placed in Parmotrema subgen. Crespoa. This study demonstrates the power of reduced genome-scale data sets to resolve phylogenetic relationships with high support. Due to lower costs, target enrichment methods provide a promising avenue for phylogenetic studies including larger taxonomic/specimen sampling than whole genome data would allow.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

IMA Genome - F13: Draft genome sequences of Ambrosiella cleistominuta, Cercospora brassicicola, C. citrullina, Physcia stellaris, and Teratosphaeria pseudoeucalypti

Draft genomes of the fungal species Ambrosiella cleistominuta, Cercospora brassicicola, C. citrullina, Physcia stellaris, and Teratosphaeria pseudoeucalypti are presented. Physcia stellaris is an important lichen forming fungus and Ambrosiella cleistominuta is an ambrosia beetle symbiont. Cercospora brassicicola and C. citrullina are agriculturally relevant plant pathogens that cause leaf-spots in brassicaceous vegetables and cucurbits respectively. Teratosphaeria pseudoeucalypti causes severe leaf blight and defoliation of Eucalyptus trees. These genomes provide a valuable resource for understanding the molecular processes in these economically important fungi.

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.

Bacterial line of defense in Dirinaria lichen from two different ecosystems: First genomic insights of its mycobiont Dirinaria sp. GBRC AP01

Lichens have been widely studied for their symbiotic properties and for the secondary metabolites production by its fungal symbiont. Recent molecular studies have confirmed coexistence of bacteria along with the fungal and algal symbionts. Direct nucleic acid study by -omics approaches is providing better insights into their structural and functional dynamics. However, genomic analysis of individual members of lichen is difficult by the conventional approach. Hence, genome assembly from metagenome data needs standardization in the eukaryotic system like lichens. The present study aimed at metagenomic characterization of rock associated lichen Dirinaria collected from Kutch and Dang regions of Gujarat, followed by genome reconstruction and annotation of the mycobiont Dirinaria. The regions considered in the study are eco-geographically highly variant. The results revealed higher alpha diversity in the dry region Kutch as compared to the tropical forest associated lichen from Dang. Ascomycota was the most abundant eukaryote while Proteobacteria dominated the bacterial population. There were 23 genera observed only in the Kutch lichen (KL) and one genus viz., Candidatus Vecturithrix unique to the Dang lichen (DL). The exclusive bacterial genera in the Kutch mostly belonged to groups reported for stress tolerance and earlier isolated from lithobionts of extreme niches. The assembled data of KL & DL were further used for genome reconstruction of Dirinaria sp. using GC and tetra-pentamer parameters and reassembly that resulted into a final draft genome of 31.7 Mb and 9556 predicted genes. Twenty-eight biosynthesis gene clusters were predicted that included genes for polyketide, indole and terpene synthesis. Association analysis of bacteria and mycobiont revealed 8 pathways specific to bacteria with implications in lichen symbiosis and environment interaction. The study provides the first draft genome of the entire fungal Dirinaria genus and provides insights into the Dirinaria lichen metagenome from Gujarat region.

Draft genome sequences of five Calonectria species from Eucalyptus plantations in China, Celoporthe dispersa, Sporothrix phasma and Alectoria sarmentosa

Draft genome sequences of five Calonectria species [including Calonectria aciculata, C. crousiana, C. fujianensis, C. honghensis and C. pseudoturangicola], Celoporthe dispersa, Sporothrix phasma and Alectoria sarmentosa are presented. Species of Calonectria are the causal agents of Eucalyptus leaf blight disease, threatening the growth and sustainability of Eucalyptus plantations in China. Celoporthe dispersa is the causal agent of stem canker in native Syzygium cordatum and exotic Tibouchina granulosa in South Africa. Sporothrix phasma was first discovered in the infructescences of Protea laurifolia and Protea neriifolia in South Africa. Alectoria sarmentosa is fruticose lichen belongs to the alectorioid clade of the family Parmeliaceae. The availability of these genome sequences will facilitate future studies on the systematics, population genetics, and genomics of these fungi.

Improved metagenomic analysis with Kraken 2

Although Kraken's k-mer-based approach provides a fast taxonomic classification of metagenomic sequence data, its large memory requirements can be limiting for some applications. Kraken 2 improves upon Kraken 1 by reducing memory usage by 85%, allowing greater amounts of reference genomic data to be used, while maintaining high accuracy and increasing speed fivefold. Kraken 2 also introduces a translated search mode, providing increased sensitivity in viral metagenomics analysis.

Improved metagenomic analysis with Kraken 2

Although Kraken’s k-mer-based approach provides a fast taxonomic classification of metagenomic sequence data, its large memory requirements can be limiting for some applications. Kraken 2 improves upon Kraken 1 by reducing memory usage by 85%, allowing greater amounts of reference genomic data to be used, while maintaining high accuracy and increasing speed fivefold. Kraken 2 also introduces a translated search mode, providing increased sensitivity in viral metagenomics analysis.

Fungal evolution: major ecological adaptations and evolutionary transitions

Fungi are a highly diverse group of heterotrophic eukaryotes characterized by the absence of phagotrophy and the presence of a chitinous cell wall. While unicellular fungi are far from rare, part of the evolutionary success of the group resides in their ability to grow indefinitely as a cylindrical multinucleated cell (hypha). Armed with these morphological traits and with an extremely high metabolical diversity, fungi have conquered numerous ecological niches and have shaped a whole world of interactions with other living organisms. Herein we survey the main evolutionary and ecological processes that have guided fungal diversity. We will first review the ecology and evolution of the zoosporic lineages and the process of terrestrialization, as one of the major evolutionary transitions in this kingdom. Several plausible scenarios have been proposed for fungal terrestralization and we here propose a new scenario, which considers icy environments as a transitory niche between water and emerged land. We then focus on exploring the main ecological relationships of Fungi with other organisms (other fungi, protozoans, animals and plants), as well as the origin of adaptations to certain specialized ecological niches within the group (lichens, black fungi and yeasts). Throughout this review we use an evolutionary and comparative-genomics perspective to understand fungal ecological diversity. Finally, we highlight the importance of genome-enabled inferences to envision plausible narratives and scenarios for important transitions.

Multiple historical processes obscure phylogenetic relationships in a taxonomically difficult group (Lobariaceae, Ascomycota)

In the age of next-generation sequencing, the number of loci available for phylogenetic analyses has increased by orders of magnitude. But despite this dramatic increase in the amount of data, some phylogenomic studies have revealed rampant gene-tree discordance that can be caused by many historical processes, such as rapid diversification, gene duplication, or reticulate evolution. We used a target enrichment approach to sample 400 single-copy nuclear genes and estimate the phylogenetic relationships of 13 genera in the lichen-forming family Lobariaceae to address the effect of data type (nucleotides and amino acids) and phylogenetic reconstruction method (concatenation and species tree approaches). Furthermore, we examined datasets for evidence of historical processes, such as rapid diversification and reticulate evolution. We found incongruence associated with sequence data types (nucleotide vs. amino acid sequences) and with different methods of phylogenetic reconstruction (species tree vs. concatenation). The resulting phylogenetic trees provided evidence for rapid and reticulate evolution based on extremely short branches in the backbone of the phylogenies. The observed rapid and reticulate diversifications may explain conflicts among gene trees and the challenges to resolving evolutionary relationships. Based on divergence times, the diversification at the backbone occurred near the Cretaceous-Paleogene (K-Pg) boundary (65 Mya) which is consistent with other rapid diversifications in the tree of life. Although some phylogenetic relationships within the Lobariaceae family remain with low support, even with our powerful phylogenomic dataset of up to 376 genes, our use of target-capturing data allowed for the novel exploration of the mechanisms underlying phylogenetic and systematic incongruence.

Whole-Genome Sequence Data Uncover Widespread Heterothallism in the Largest Group of Lichen-Forming Fungi

Fungal reproduction is regulated by the mating-type (MAT1) locus, which typically comprises two idiomorphic genes. The presence of one or both allelic variants at the locus determines the reproductive strategy in fungi—homothallism versus heterothallism. It has been hypothesized that self-fertility via homothallism is widespread in lichen-forming fungi. To test this hypothesis, we characterized the MAT1 locus of 41 genomes of lichen-forming fungi representing a wide range of growth forms and reproductive strategies in the class Lecanoromycetes, the largest group of lichen-forming fungi. Our results show the complete lack of genetic homothallism suggesting that lichens evolved from a heterothallic ancestor. We argue that this may be related to the symbiotic lifestyle of these fungi, and may be a key innovation that has contributed to the accelerated diversification rates in this fungal group.

Differential gene expression associated with fungal trophic shifts along the senescence gradient of the moss Dicranum scoparium

Bryophytes harbour microbiomes, including diverse communities of fungi. The molecular mechanisms by which perennial mosses interact with these fungal partners along their senescence gradients are unknown, yet this is an ideal system to study variation in gene expression associated with trophic state transitions. We investigated differentially expressed genes of fungal communities and their host Dicranum scoparium across its naturally occurring senescence gradient using a metatranscriptomic approach. Higher activity of fungal nutrient-related (carbon, nitrogen, phosphorus and sulfur) transporters and Carbohydrate-Active enZyme (CAZy) genes was detected toward the bottom, partially decomposed, layer of the moss. The most prominent variation in the expression levels of fungal nutrient transporters was from inorganic nitrogen-related transporters, whereas the breakdown of organonitrogens was detected as the most enriched gene ontology term for the host D. scoparium, for those transcripts having higher expression in the partially decomposed layer. The abundance of bacterial rRNA transcripts suggested that more living members of Cyanobacteria are associated with the photosynthetic layer of D. scoparium, while members of Rhizobiales are detected throughout the gametophytes. Plant genes for specific fungal-plant communication, including defense responses, were differentially expressed, suggesting that different genetic pathways are involved in plant-microbe crosstalk in photosynthetic tissues compared to partially decomposed tissues.

Biosynthetic Gene Content of the 'Perfume Lichens' Evernia prunastri and Pseudevernia furfuracea

Lichen-forming fungi produce a vast number of unique natural products with a wide variety of biological activities and human uses. Although lichens have remarkable potential in natural product research and industry, the molecular mechanisms underlying the biosynthesis of lichen metabolites are poorly understood. Here we use genome mining and comparative genomics to assess biosynthetic gene clusters and their putative regulators in the genomes of two lichen-forming fungi, which have substantial commercial value in the perfume industry, Evernia prunastri and Pseudevernia furfuracea. We report a total of 80 biosynthetic gene clusters (polyketide synthases (PKS), non-ribosomal peptide synthetases and terpene synthases) in E. prunastri and 51 in P. furfuracea. We present an in-depth comparison of 11 clusters, which show high homology between the two species. A ketosynthase (KS) phylogeny shows that biosynthetic gene clusters from E. prunastri and P. furfuracea are widespread across the Fungi. The phylogeny includes 15 genomes of lichenized fungi and all fungal PKSs with known functions from the MIBiG database. Phylogenetically closely related KS domains predict not only similar PKS architecture but also similar cluster architecture. Our study highlights the untapped biosynthetic richness of lichen-forming fungi, provides new insights into lichen biosynthetic pathways and facilitates heterologous expression of lichen biosynthetic gene clusters.

Limited Correlation of Shotgun Metagenomics Following Host Depletion and Routine Diagnostics for Viruses and Bacteria in Low Concentrated Surrogate and Clinical Samples

The etiologic cause of encephalitis, meningitis or meningo-encephalitis is unknown in up to 70% of cases. Clinical shotgun metagenomics combined with host depletion is a promising technique to identify infectious etiologies of central nervous system (CNS) infections. We developed a straightforward eukaryotic host nucleic acid depletion method that preserves intact viruses and bacteria for subsequent shotgun metagenomics screening of clinical samples, focusing on cerebrospinal fluid (CSF). A surrogate CSF sample for a CNS infection paradigm was used to evaluate the proposed depletion method consisting of selective host cell lysis, followed by enzymatic degradation of the liberated genomic DNA for final depletion with paramagnetic beads. Extractives were subjected to reverse transcription, followed by whole genome amplification and next generation sequencing. The effectiveness of the host depletion method was demonstrated in surrogate CSF samples spiked with three 1:100 dilutions of Influenza A H3N2 virus (qPCR Ct-values 20.7, 28.8, >42/negative). Compared to the native samples, host depletion increased the amount of the virus subtype reads by factor 7127 and 132, respectively, while in the qPCR negative sample zero vs. 31 (1.4E-4 %) virus subtype reads were detected (native vs. depleted). The workflow was applied to thirteen CSF samples of patients with meningo-/encephalitis (two bacterial, eleven viral etiologies), a serum of an Andes virus infection and a nose swab of a common cold patient. Unlike surrogate samples, host depletion of the thirteen human CSF samples and the nose swab did not result in more reads indicating presence of damaged pathogens due to, e.g., host immune response. Nevertheless, previously diagnosed pathogens in the human CSF samples (six viruses, two bacteria), the serum, and the nose swab (Human rhinovirus A31) were detected in the depleted and/or the native samples. Unbiased evaluation of the taxonomic profiles supported the diagnosed pathogen in two native CSF samples and the native and depleted serum and nose swab, while detecting various contaminations that interfered with pathogen identification at low concentration levels. In summary, damaged pathogens and contaminations complicated analysis and interpretation of clinical shotgun metagenomics data. Still, proper consideration of these issues may enable future application of metagenomics for clinical diagnostics.

Comparing Sanger sequencing and high-throughput metabarcoding for inferring photobiont diversity in lichens

The implementation of HTS (high-throughput sequencing) approaches is rapidly changing our understanding of the lichen symbiosis, by uncovering high bacterial and fungal diversity, which is often host-specific. Recently, HTS methods revealed the presence of multiple photobionts inside a single thallus in several lichen species. This differs from Sanger technology, which typically yields a single, unambiguous algal sequence per individual. Here we compared HTS and Sanger methods for estimating the diversity of green algal symbionts within lichen thalli using 240 lichen individuals belonging to two species of lichen-forming fungi. According to HTS data, Sanger technology consistently yielded the most abundant photobiont sequence in the sample. However, if the second most abundant photobiont exceeded 30% of the total HTS reads in a sample, Sanger sequencing generally failed. Our results suggest that most lichen individuals in the two analyzed species, Lasallia hispanica and L. pustulata, indeed contain a single, predominant green algal photobiont. We conclude that Sanger sequencing is a valid approach to detect the dominant photobionts in lichen individuals and populations. We discuss which research areas in lichen ecology and evolution will continue to benefit from Sanger sequencing, and which areas will profit from HTS approaches to assessing symbiont diversity.