Characterization of a transcriptome from a non-model organism, Cladonia rangiferina, the grey reindeer lichen, using high-throughput next generation sequencing and EST sequence data

Comparative Transcriptomic and Proteomic Analyses Provide New Insights into the Tolerance to Cyclic Dehydration in a Lichen Phycobiont

Desiccation tolerance (DT) is relatively frequent in non-vascular plants and green algae. However, it is poorly understood how successive dehydration/rehydration (D/R) cycles shape their transcriptomes and proteomes. Here, we report a comprehensive analysis of adjustments on both transcript and protein profiles in response to successive D/R cycles in Coccomyxa simplex (Csol), isolated from the lichen Solorina saccata. A total of 1833 transcripts and 2332 proteins were differentially abundant as a consequence of D/R; however, only 315 of these transcripts/proteins showed similar trends. Variations in both transcriptomes and proteomes along D/R cycles together with functional analyses revealed an extensive decrease in transcript and protein levels during dehydration, most of them involved in gene expression, metabolism, substance transport, signalling and folding catalysis, among other cellular functions. At the same time, a series of protective transcripts/proteins, such as those related to antioxidant defence, polyol metabolism and autophagy, was upregulated during dehydration. Overall, our results show a transient decrease in most cellular functions as a result of drying and a gradual reactivation of specific cell processes to accommodate the hydration status along successive D/R cycles. This study provides new insights into key mechanisms involved in the DT of Csol and probably other dehydration-tolerant microalgae. In addition, functionally characterising the high number of genes/proteins of unknown functions found in this study may lead to the discovery of new DT mechanisms.

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.

CPPred: coding potential prediction based on the global description of RNA sequence

The rapid and accurate approach to distinguish between coding RNAs and ncRNAs has been playing a critical role in analyzing thousands of novel transcripts, which have been generated in recent years by next-generation sequencing technology. Previously developed methods CPAT, CPC2 and PLEK can distinguish coding RNAs and ncRNAs very well, but poorly distinguish between small coding RNAs and small ncRNAs. Herein, we report an approach, CPPred (coding potential prediction), which is based on SVM classifier and multiple sequence features including novel RNA features encoded by the global description. The CPPred can better distinguish not only between coding RNAs and ncRNAs, but also between small coding RNAs and small ncRNAs than the state-of-the-art methods due to the addition of the novel RNA features. A recent study proposes 1335 novel human coding RNAs from a large number of RNA-seq datasets. However, only 119 transcripts are predicted as coding RNAs by the CPPred. In fact, almost all proposed novel coding RNAs are ncRNAs (91.1%), which is consistent with previous reports. Remarkably, we also reveal that the global description of encoding features (T2, C0 and GC) plays an important role in the prediction of coding potential.

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.

CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features

With advances in next-generation sequencing technologies, numerous novel transcripts in a large number of organisms have been identified. With the goal of fast, accurate assessment of the coding ability of RNA transcripts, we upgraded the coding potential calculator CPC1 to CPC2. CPC2 runs ∼1000 times faster than CPC1 and exhibits superior accuracy compared with CPC1, especially for long non-coding transcripts. Moreover, the model of CPC2 is species-neutral, making it feasible for ever-growing non-model organism transcriptomes. A mobile-friendly web server, as well as a downloadable standalone package, is freely available at http://cpc2.cbi.pku.edu.cn.

The CRISPR/Cas9 System for Targeted Genome Engineering in Free-Living Fungi: Advances and Opportunities for Lichenized Fungi

Studies using whole genome sequencing, computational and gene expression, targeted genome engineering techniques for generating site-specific sequence alterations through non-homologous end joining (NHEJ) by genomic double-strand break (DSB) repair pathway with high precision, resulting in gene inactivation have elucidated the complexity of gene expression, and metabolic pathways in fungi. These tools and the data generated are crucial for precise generation of fungal products such as enzymes, secondary metabolites, antibiotics etc. Artificially engineered molecular scissors, zinc finger nucleases (ZFNs), Transcriptional activator-like effector nucleases (TALENs; that use protein motifs for DNA sequence recognition in the genome) and CRISPR associated protein 9 (Cas9;CRISPR/Cas9) system (RNA-DNA recognition) are being used in achieving targeted genome modifications for modifying traits in free-living fungal systems. Here, we discuss the recent research breakthroughs and developments which utilize CRISPR/Cas9 in the metabolic engineering of free-living fungi for the biosynthesis of secondary metabolites, enzyme production, antibiotics and to develop resistance against post-harvest browning of edible mushrooms and fungal pathogenesis. We also discuss the potential and advantages of using targeted genome engineering in lichenized fungal (mycobiont) cultures to enhance their growth and secondary metabolite production in vitro can be complemented by other molecular approaches.

Lichen Biosynthetic Gene Clusters. Part I. Genome Sequencing Reveals a Rich Biosynthetic Potential

Lichens are symbionts of fungi and algae that produce diverse secondary metabolites with useful properties. Little is known of lichen natural product biosynthesis because of the challenges of working with lichenizing fungi. We describe the first attempt to comprehensively profile the genetic secondary metabolome of a lichenizing fungus. An Illumina platform combined with the Antibiotics and Secondary Metabolites Analysis Shell (FungiSMASH, version 4.0) was used to sequence and annotate assembled contigs of the fungal partner of Cladonia uncialis. Up to 48 putative gene clusters are described comprising type I and type III polyketide synthases (PKS), nonribosomal peptide synthetases (NRPS), hybrid PKS-NRPS, and terpene synthases. The number of gene clusters revealed by this work dwarfs the number of known secondary metabolites from C. uncialis, suggesting that lichenizing fungi have an unexplored biosynthetic potential.

The cost of surviving nitrogen excess: energy and protein demand in the lichen Cladonia portentosa as revealed by proteomic analysis

Different nitrogen forms affect different metabolic pathways in lichens. In particular, the most relevant changes in protein expression were observed in the fungal partner, with NO ₃^- mostly affecting the energetic metabolism and NH ₄⁺ affecting transport and regulation of proteins and the energetic metabolism much more than NO ₃^- did. Excess deposition of reactive nitrogen is a well-known agent of stress for lichens, but which symbiont is most affected and how, remains a mystery. Using proteomics can expand our understanding of stress effects on lichens. We investigated the effects of different doses and forms of reactive nitrogen, with and without supplementary phosphorus and potassium, on the proteome of the lichen Cladonia portentosa growing in a 'real-world' simulation of nitrogen deposition. Protein expression changed with the nitrogen treatments but mostly in the fungal partner, with NO₃^- mainly affecting the energetic metabolism and NH₄⁺ also affecting the protein synthesis machinery. The photobiont mainly responded overexpressing proteins involved in energy production. This suggests that in response to nitrogen stress, the photobiont mainly supports the defensive mechanisms initiated by the mycobiont with an increased energy production. Such surplus energy is then used by the cell to maintain functionality in the presence of NO₃^-, while a futile cycle of protein production can be hypothesized to be induced by NH₄⁺ excess. External supply of potassium and phosphorus influenced differently the responses of particular enzymes, likely reflecting the many processes in which potassium exerts a regulatory function.

Lichenized Fungi and the Evolution of Symbiotic Organization

Lichen symbioses comprise a fascinating relationship between algae and fungi. The lichen symbiotic lifestyle evolved early in the evolution of ascomycetes and is also known from a few basidiomycetes. The ascomycete lineages have diversified in the lichenized stage to give rise to a tremendous variety of morphologies. Their thalli are often internally complex and stratified for optimized integration of algal and fungal metabolisms. Thalli are frequently colonized by specific nonlichenized fungi and occasionally also by other lichens. Microscopy has revealed various ways these fungi interact with their hosts. Besides the morphologically recognizable diversity of the lichen mycobionts and lichenicolous (lichen-inhabiting) fungi, many other microorganisms including other fungi and bacterial communities are now detected in lichens by culture-dependent and culture-independent approaches. The application of multi-omics approaches, refined microscopic techniques, and physiological studies has added to our knowledge of lichens, not only about the taxa involved in the lichen interactions, but also about their functions.

Symbiont modulates expression of specific gene categories in Angomonas deanei

Trypanosomatids are parasites that cause disease in humans, animals, and plants. Most are non-pathogenic and some harbor a symbiotic bacterium. Endosymbiosis is part of the evolutionary process of vital cell functions such as respiration and photosynthesis. Angomonas deanei is an example of a symbiont-containing trypanosomatid. In this paper, we sought to investigate how symbionts influence host cells by characterising and comparing the transcriptomes of the symbiont-containing A. deanei (wild type) and the symbiont-free aposymbiotic strains. The comparison revealed that the presence of the symbiont modulates several differentially expressed genes. Empirical analysis of differential gene expression showed that 216 of the 7625 modulated genes were significantly changed. Finally, gene set enrichment analysis revealed that the largest categories of genes that downregulated in the absence of the symbiont were those involved in oxidation-reduction process, ATP hydrolysis coupled proton transport and glycolysis. In contrast, among the upregulated gene categories were those involved in proteolysis, microtubule-based movement, and cellular metabolic process. Our results provide valuable information for dissecting the mechanism of endosymbiosis in A. deanei.

New features of desiccation tolerance in the lichen photobiont Trebouxia gelatinosa are revealed by a transcriptomic approach

Trebouxia is the most common lichen-forming genus of aero-terrestrial green algae and all its species are desiccation tolerant (DT). The molecular bases of this remarkable adaptation are, however, still largely unknown. We applied a transcriptomic approach to a common member of the genus, T. gelatinosa, to investigate the alteration of gene expression occurring after dehydration and subsequent rehydration in comparison to cells kept constantly hydrated. We sequenced, de novo assembled and annotated the transcriptome of axenically cultured T. gelatinosa by using Illumina sequencing technology. We tracked the expression profiles of over 13,000 protein-coding transcripts. During the dehydration/rehydration cycle c. 92 % of the total protein-coding transcripts displayed a stable expression, suggesting that the desiccation tolerance of T. gelatinosa mostly relies on constitutive mechanisms. Dehydration and rehydration affected mainly the gene expression for components of the photosynthetic apparatus, the ROS-scavenging system, Heat Shock Proteins, aquaporins, expansins, and desiccation related proteins (DRPs), which are highly diversified in T. gelatinosa, whereas Late Embryogenesis Abundant Proteins were not affected. Only some of these phenomena were previously observed in other DT green algae, bryophytes and resurrection plants, other traits being distinctive of T. gelatinosa, and perhaps related to its symbiotic lifestyle. Finally, the phylogenetic inference extended to DRPs of other chlorophytes, embryophytes and bacteria clearly pointed out that DRPs of chlorophytes are not orthologous to those of embryophytes: some of them were likely acquired through horizontal gene transfer from extremophile bacteria which live in symbiosis within the lichen thallus.

Exploring functional contexts of symbiotic sustain within lichen-associated bacteria by comparative omics

Symbioses represent a frequent and successful lifestyle on earth and lichens are one of their classic examples. Recently, bacterial communities were identified as stable, specific and structurally integrated partners of the lichen symbiosis, but their role has remained largely elusive in comparison to the well-known functions of the fungal and algal partners. We have explored the metabolic potentials of the microbiome using the lung lichen Lobaria pulmonaria as the model. Metagenomic and proteomic data were comparatively assessed and visualized by Voronoi treemaps. The study was complemented with molecular, microscopic and physiological assays. We have found that more than 800 bacterial species have the ability to contribute multiple aspects to the symbiotic system, including essential functions such as (i) nutrient supply, especially nitrogen, phosphorous and sulfur, (ii) resistance against biotic stress factors (that is, pathogen defense), (iii) resistance against abiotic factors, (iv) support of photosynthesis by provision of vitamin B12, (v) fungal and algal growth support by provision of hormones, (vi) detoxification of metabolites, and (vii) degradation of older parts of the lichen thallus. Our findings showed the potential of lichen-associated bacteria to interact with the fungal as well as algal partner to support health, growth and fitness of their hosts. We developed a model of the symbiosis depicting the functional multi-player network of the participants, and argue that the strategy of functional diversification in lichens supports the longevity and persistence of lichens under extreme and changing ecological conditions.

Genome characteristics reveal the impact of lichenization on lichen-forming fungus Endocarpon pusillum Hedwig (Verrucariales, Ascomycota)

BACKGROUND

Lichen is a classic mutualistic organism and the lichenization is one of the fungal symbioses. The lichen-forming fungus Endocarpon pusillum is living in symbiosis with the green alga Diplosphaera chodatii Bialsuknia as a lichen in the arid regions.

RESULTS

454 and Illumina technologies were used to sequence the genome of E. pusillum. A total of 9,285 genes were annotated in the 37.5 Mb genome of E. pusillum. Analyses of the genes provided direct molecular evidence for certain natural characteristics, such as homothallic reproduction and drought-tolerance. Comparative genomics analysis indicated that the expansion and contraction of some protein families in the E. pusillum genome reflect the specific relationship with its photosynthetic partner (D. chodatii). Co-culture experiments using the lichen-forming fungus E. pusillum and its algal partner allowed the functional identification of genes involved in the nitrogen and carbon transfer between both symbionts, and three lectins without signal peptide domains were found to be essential for the symbiotic recognition in the lichen; interestingly, the ratio of the biomass of both lichen-forming fungus and its photosynthetic partner and their contact time were found to be important for the interaction between these two symbionts.

CONCLUSIONS

The present study lays a genomic analysis of the lichen-forming fungus E. pusillum for demonstrating its general biological features and the traits of the interaction between this fungus and its photosynthetic partner D. chodatii, and will provide research basis for investigating the nature of its drought resistance and symbiosis.

Whole transcriptome characterization of the effects of dehydration and rehydration on Cladonia rangiferina, the grey reindeer lichen

BACKGROUND

Lichens are symbiotic organisms with a fungal and an algal or a cyanobacterial partner. Lichens inhabit some of the harshest climates on earth and most lichen species are desiccation-tolerant. Lichen desiccation-tolerance has been studied at the biochemical level and through proteomics, but the underlying molecular genetic mechanisms remain largely unexplored. The objective of our study was to examine the effects of dehydration and rehydration on the gene expression of Cladonia rangiferina.

RESULTS

Samples of C. rangiferina were collected at several time points during both the dehydration and rehydration process and the gene expression intensities were measured using a custom DNA microarray. Several genes, which were differentially expressed in one or more time points, were identified. The microarray results were validated using qRT-PCR analysis. Enrichment analysis of differentially expressed transcripts was also performed to identify the Gene Ontology terms most associated with the rehydration and dehydration process.

CONCLUSIONS

Our data identify differential expression patterns for hundreds of genes that are modulated during dehydration and rehydration in Cladonia rangiferina. These dehydration and rehydration events clearly differ from each other at the molecular level and the largest changes to gene expression are observed within minutes following rehydration. Distinct changes are observed during the earliest stage of rehydration and the mechanisms not appear to be shared with the later stages of wetting or with drying. Several of the most differentially expressed genes are similar to genes identified in previous studies that have investigated the molecular mechanisms of other desiccation-tolerant organisms. We present here the first microarray experiment for any lichen species and have for the first time studied the genetic mechanisms behind lichen desiccation-tolerance at the whole transcriptome level.