- Research article
- Open Access
Genome characteristics reveal the impact of lichenization on lichen-forming fungus Endocarpon pusillum Hedwig (Verrucariales, Ascomycota)
BMC Genomics volume 15, Article number: 34 (2014)
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.
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.
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.
A lichen is a symbiotic association of a fungus (mycobiont) and a photosynthetic partner (photobiont), which may be an alga (phycobiont) or a cyanobacterium (cyanobiont). In the association the fungus produces a thallus, or body, within which the photobionts are housed . Around 20% of all Fungi and 40% of all Ascomycota are lichen-forming. Recent estimates of global diversity suggest that there are between 17,500 and 20,000 species .
Most lichens and isolated lichen-forming fungi grow extremely slow, but the lichen-symbiosis is a very successful association as lichens can survive in almost all adverse terrestrial conditions . They are also famous for their particular secondary products, which are frequently used as antibacterial and antiviral compounds [4, 5]. The lichen-forming fungi differ from non-lichenized fungi by their adaptations to symbiosis with photobiont . This mutualistic association, as called lichenization, is one of the most important fungal lifestyles and the lichenization, considered by some researchers, has evolved many times in the phylogeny of fungi [7, 8], and also some major fungal lineages may have derived from lichen symbiotic ancestors .
The principal problem about lichenization is the necessity of fungal propagules meeting a suitable photosynthetic partner for the resynthesis of the symbiosis [3, 10]. The recognition step is complicated, involving many morphological and molecular changes. Scanning electron microscopy (SEM) has been used to investigate the changes in morphology during the early resynthesis of the lichen thallus [11–13]. However, few studies have explored the resynthesis events in lichen using molecular tools. Previous studies suggest that the mycobiont-derived lectins (sugar-specific, cell agglutinating proteins) may play a key role in recognition [14–17] between both symbionts.
The interdependent relationship between mycobiont and photobiont is the foundation of lichenization, which required for both symbionts to maintain each other. In the lichen thallus, the photobiont provides its mycobiont with photosynthetic products [18, 19], previous reports showed that lichen-forming fungi absorb polyol (ribitol, sorbitol, and erythritol) or glucose from algae or cyanobacteria, respectively [20, 21], and in most green algae lichens, the hyphae of the mycobionts wrap tightly around photobiont cells, thereby protecting the photobiont cells from a range of biotic and abiotic stress, including drought, high light and mechanical damage. The protection from light-injury is associated with secondary metabolic substances, such as melanins, produced by lichen-forming fungi . However, little is known about the signals and mechanisms that lead to symbiotic recognition and maintenance in lichen, although it can be predicted that some metabolic products and macromolecules are essential.
The whole life cycle of lichen is rarely observed in the laboratory or in nature [11, 13]. The detail of lichenization, whereby a lichen-forming fungus contacts with a compatible photosynthetic partner, recognizes it and captures it, has become a hot topic in recent studies [23, 24]. Endocarpon pusillum is the most successful cases on artificial resynthesis of the fertile (perithecia-bearing) thalli from isolated mycobiont and phycobiont , and it will be very useful to reveal the origin of symbiosis between fungi and photosynthetic organisms.
Endocarpon is a special genus, which have hymenial algal cells in their perithecia and the ascospores are discharged together with the hymenial algal cells . The systematics and physiology of E. pusillum (Figure 1) have been studied well. The mycobiont E. pusillum exhibits much stronger desiccation-tolerant than other non-lichenized fungi as it can survive for 7 months under desiccation stress in combination with starvation stress . Although a large number of fungal genomes have been published, no lichen-forming is included, and only two mitochondrial genomes from lichen-forming fungi (Peltigera membranacea and Peltigera malacea) and a transcriptome from the lichen Cladonia rangiferina have been reported until now [27, 28]. Therefore, the genome of E. pusillum was sequenced and analyzed to ascertain the biological features of this lichen-forming fungus and the traits of lichenization.
Results and discussion
General features of the genome
The genome of the lichen-forming fungus E. pusillum was sequenced to about 78-fold coverage using both 454 and Illumina technologies (Additional file 1: Table S1). All sequences were assembled into 908 scaffolds (> 2 kb; N50, 178 kb) containing 1,731 contigs, with a genome size of 37.5 Mb (Table 1), which was almost identical to the result calculated by real-time polymerase chain reaction (PCR) . The average GC content of the genome is 46.1%, and exonic region has a 4% higher GC content than the intronic region. Repetitive sequences represent 15% of the genome. A circular map was generated for the 30 largest scaffolds to illustrate the genome features more clearly (Figure 2).
The whole project has been deposited at DDBJ/EMBL/GenBank under accession number APWS00000000. A total of 9,285 protein-coding genes were predicted, and 1,479 (15.7%) of these genes have no significant matches to known proteins from public databases. A total of 2,754, 3,787 and 7,589 proteins were assigned to Gene Ontology (GO) terms, the eukaryotic orthologous groups (KOG) and functional catalogue (FunCat) databases, respectively. The distributions of the top 10 GO, KOG, and FunCat terms of the sequences are presented in Figure 3.
Phylogenetic analysis of E. pusillum
The orthologous genes from the lichen-forming fungus E. pusillum and 14 other non lichen-forming fungi whose genomes were available were identified using Inparanoid  with default parameters (Figure 4). Phylogenetic analysis was performed using 1,893 single-copy orthologous genes identified among the genomes of above mentioned 15 fungi from the subkingdom Dikarya, and a linearized phylogenetic tree was constructed with estimates of the divergence times among these taxa (Figure 4). The phylogenomic analysis shows that the E. pusillum lineage is more closely related to the human pathogen Exophiala dermatitidis, which is the anamorph species of Capronia belonging to the Chaetothyriales, and the E. pusillum belongs to the Verrucariales. Both the orders belong to the same subclass Chaetothyriomycetidae, and the divergence between the lichenized Verrucariales and nonlichenized Chaetothyriales occurred approximately 131 million years (Myr) ago (Figure 4). This result is consistence with the phylogenetic analysis for Ascomycota using ribosomal RNA .
Repeat-induced point mutation
Repeat-induced point mutation (RIP) is a gene silencing mechanism that can cause C-G to T-A mutations on repetitive DNA sequences, and the mutations from C to T mostly occur at CpA dinucleotides . According to the method proposed by Margolin et al., sequences with a high TA/AT ratio (>0.89) and a low (CA + TG)/(AC + GT) ratio (< 1.03) are thought to indicate RIP .
A quantitative alignment-based method, RIPCAL , was used to search for evidence for RIP in E. pusillum genome. The RIP indices of the repetitive sequence are 1.46927 and 0.697168, and those of non-repetitive sequence are 0.565627 and 1.42277, which indicates that the regions of repetitive sequence in E. pusillum have undergone RIP. In addition, the percentage of genes in multigene families shows that E. pusillum has a low proportion of genes in multigene families (Figure 5A). The analysis of amino acid similarities in multigene families among eight fungi indicates that there are fewer best-matching (≥ 80%) genes in the E. pusillum genome (Figure 5B), which implies that some genes in the multigene families in E. pusillum are RIP mutated. However, the function of RIP in lichen-forming fungi still remains unclear.
The mating systems and reproduction
The genes involved in the mating process were identified in E. pusillum genome. The result shows that there is a single mating type (MAT) locus containing both a MATα and HMG domain in the genome, which provides the molecular evidence for homothallic lifestyle of E. pusillum (Figure 6). Genes for conidiophore development and asexual reproduction were also identified, together with those involved in fruit body development and sexual reproduction (Additional file 1: Table S2), and genes with the function of pheromone precursor (alpha-factor like), pheromone receptor (for alpha-factor like pheromone), ascus development (rhamnogalacturonase B), and receptors preventing improper sexual development have been lost in E. pusillum. However, it need more experimental evidence to determine whether these missing genes contribute to the special reproduction mode of E. pusillum.
Secreted proteins are important effectors that modulate the interaction between pathogenic microbes and hosts [35–37], and many small secreted proteins (SSPs) are virulence factors . However, E. pusillum has fewer secreted proteins than other phytopathogenic fungi (e.g. Fusarium graminearum, Magnaporthe oryzae and Metarhizium anisopliae) according to the blast results (Additional file 1: Table S3), and the function of 66% of SSPs is unknown. Lichen-forming fungus and pathogenic fungi are different, even though they both live with another organism and absorb nutrients from their partner. The analysis of the secreted proteins suggests that the SSP-mediated interaction between symbionts in E. pusillum is weaker than the virulence in parasitism (pathogenic fungi).
To date, several polyketide synthase (PKS) genes have been identified in different lichen-forming fungi [39–42], but no non-ribosomal peptide synthetase (NRPS) genes have been found. Analysis of the lichen-forming fungus E. pusillum genome revealed 2 NRPS genes and 14 PKS genes (Additional file 1: Table S4). The information of genes for secondary metabolism is helpful to the research of the secondary metabolites of E. pusillum and the synthesis conditions.
Candidate genes for drought tolerance
At least 93 genes involving in these drought resistant mechanisms were identified in the genome of E. pusillum (Additional file 1: Table S5), which will provide us many important clues to resolve the question about why E. pusillum can survive in extreme drought environment. Recently, most studies on mechanisms about desiccation-tolerance in lichens have focused on the scavenging reactive oxygen species . However, genes involved in regulating osmotic pressure, correcting protein misfolding, and scavenging the reactive oxygen species were identified in E. pusillum genome, which were worthy of further investigation.
Multigene families expansion
Up to 1,889 protein families (containing at least two genes in all selected species) were identified from E. pusillum and 14 other fungi using Markov cluster (MCL) method. Based on BadiRate analysis, we found that 1,129 protein families were expanded and 760 protein families had undergone contraction in E. pusillum. The outlier results generated by BadiRate analysis are presented in Table 2. The expansion of these protein families might reflect specific characteristics of lichen-forming fungus E. pusillum.
Protein families involved in signal transduction
Compared with other fungi, the lichen-forming fungi are special because they can form mutualistic symbiosis with photosynthetic organisms and this specialty is reflected in its genome. Little is known about the functional genes in lichens, and most expanded protein families are function-unknown in E. pusillum (Table 2); however, it is likely that many of these expansion families are related to symbiosis. For example, the expressions of WD-repeat domain-containing proteins are up-regulated in the early developmental stages of lichen-symbiosis in Cladonia grayi.
Similarly, G protein-coupled receptors (GPCRs), which sense molecules outside the cell and activate signal transduction pathways inside the cell , are also likely to provide lichen-forming fungi with a strong ability to response to such signals, and all GPCR families have undergone obvious expansion (Table 3).
Heterokaryon incompatibility protein families
There are 261 genes in the E. pusillum genome that were annotated to have heterokaryon incompatibility (HET)-related functions, and 182 of them are homologous to those in N. crassa and Podospora anserina (Additional file 1: Table S6). HI (Heterokaryon incompatibility) is a common characteristic among filamentous fungi; it can prevent the formation of heterokaryotic cells in which two different genomes coexist. Once the fusion between two individuals with incompatible het loci occurs, the HET genes would trigger the HI reaction, which is characterized as growth inhibition, repression of asexual sporulation, hyphal compartmentation and death of the heterokaryotic cell [45–47]. However, the biological significance of HI is still unknown. The expansion of HET protein families suggests that E. pusillum is likely to have a strict regulation of its vegetation , and may represent a strategy for stable genotype by defending against the transfer of exogenous genes.
Lineage-specific protein families involved in self-splicing of insertions
The genes encoding ribonuclease H and transposase, which are predicted to be able to target some insertions, such as group I introns that can be stably integrated into the genome following the reverse splicing reaction , belong to lineage-specific multigene families in E. pusillum (Table 2).
The genomic sequencing data reveal that there are three group I introns in the large-subunit ribosomal DNA (LSU rDNA) of E. pusillum Z07020 at position 856, 2169, and 2721 (corresponding to the sequence of Pichia methylivora with NCBI accession number EU011611), and one group I intron in the small-subunit ribosomal DNA (SSU rDNA) at position 1769 (corresponding to the sequence of Saccharomyces cerevisiae with NCBI accession number Z75578).
Group I introns can self-splice themselves from an RNA transcript, requiring protein factors to facilitate the correct folding of the ribozyme core . There are multiple group I introns present in the nuclear rDNA of many lichen-forming fungi [50–53], and they are considered to decrease the growth rate of lichen-forming fungi by interfering with rRNA maturation . It is considered an adaption for lichenization: most lichens grow very slowly because the symbiosis would be disrupted if the mycobionts too many nutrients from the photobionts. Therefore, it is reasonable that the protein families involved in RNA reverse splicing have undergone expansion in a lichen-forming fungal genome that harbors abundant self-splicing introns.
Protein families involved in transport
The exchange of ions and metabolites between the mycobiont and photobiont is the foundation of the symbiotic association. For this reason, genes encoding transporter proteins are expected to exhibit some traits in the lichen-forming fungus genome. There are 10 genes belonging to the magnesium transport protein CorA family, which showed significant expansion in the E. pusillum genome (Table 2). The CorA proteins were the first family to be identified that could transport Mg2+ in bacteria [55, 56]. The mechanism of action of the CorA family has been well characterized in yeast and members of this family can transfer Mg2+ both into and out of the cell . Magnesium is the most important divalent cation in cells, and is particularly important for photosynthesis . Magnesium deficiency has been reported to affect plant photosynthesis and growth [59, 60]. In the lichen thallus, the special structure of the symbionts (Figure 1D–F) makes it difficult for algal cells to absorb magnesium from the environment. However, genes for magnesium transport protein are expanded in E. pusillum, which led us to suspect that during symbiosis the mycobiont could provide magnesium to the phycobiont to meet the needs of their life cycles.
A striking finding is that most nitrogen transporter families are expanded but most sugar transporter families have been lost from the E. pusillum genome (Table 4). Though the mechanisms of substance transfer in lichen are not completely understood, it has been reported that the mycobiont absorbs certain carbohydrates generated by photosynthesis of the photobiont . The reduction of sugar transporter genes, especially those for common sugars, such as glucose and fructose, which can be utilized by many organisms, suggests that the mycobiont simply absorbs certain uncommon carbohydrates from their photosynthetic partner. Therefore, the genes for many unnecessary sugar transporters have been lost during evolution. The lichen-forming fungi do not rely on common sugars because these carbon sources can be used by numerous microorganisms and lichen cannot compete with them for its slow growth and metabolism.
By contrast, the expansion of nitrogen transporters suggests that the lichen-forming fungi need to transfer (export or import) various nitrogen sources. Nitrogen is an indispensable substance for the growth of organisms, and lichens whose photosynthetic partners are green algae must obtain inorganic and organic nitrogen [62–64]. There is no evidence for nitrogen transfer between symbionts in lichen. However, owing to the particular structure of lichen thallus in which the photobiont cells are tightly entwined by mycelia (Figure 1D–F) and the photobiont cells cannot absorb substances from their surroundings directly, it can be deduced that the necessary components for their growth, such as ions and nitrogen sources, must be transferred from the mycobiont to the photobiont. The phycobiont may be have a preference for certain forms of nitrogen; for example, in the alga-bacteria association, the algae prefer NH4+ and NO3-, which are release by bacteria when cultured in organic nitrogen conditions . Thus, the mycobiont may absorb and transform the nitrogen forms that are preferred by its photosynthetic partner. Therefore, the expansion of nitrogen transporters should be related to this physiological process. Recently, it was reported that the ammonium transporter, MEP-α, in lichen-forming fungi was obtained though horizontal gene transfer , which can partly explain the mechanism of the expansion of the nitrogen transporter superfamily. This gene appears to have been lost in some groups of lichens, such as in those that are symbiotic with nitrogen-fixing cyanobacteria or those that inhabit high-nitrogen niches . However, for the lichens whose photobiont is a green alga that cannot fix nitrogen, the expansion of nitrogen transporter families is of great significance in promoting the growth of the lichen thallus, as seen in E. pusillum.
Expression analysis of symbiosis-related genes by quantitative real-time PCR (qRT-PCR)
A qRT-PCR experiment was performed to determine the changes of some symbiosis-related genes in the lichen-forming fungus at transcription level when E. pusillum pre-contacted with its photosynthetic partner, Diplosphaera chodatii.
Symbiosis-related genes selection
Thirty-two genes encoding lectins were identified in the E. pusillum genome. They belong to four superfamilies: concanavalin A-like lectin, fucose-specific lectin, mannose-binding lectin and ricin B-related lectin. For the qRT-PCR assay, six lectin genes were selected as representatives of the four superfamilies, according to the difference in the number of transmembrane helices and signal peptides (Additional file 1: Table S7).
The expansions of nitrogen transporters and the contractions of sugar transporters imply that they are closely related to maintenance of both symbionts. Hence, certain genes encoding these transporters were chosen for the qRT-PCR assay. Among them, two ammonium transporters and one nitrate transporter that displays high homology to those in other fungi (Additional file 1: Table S8) were included. It has been reported that eight sugar transporters are able to transfer photosynthesis products in other fungi, and they are predicted to be functional at the interface of lichen association [68–75]. Thus, the homologous proteins were identified in the E. pusillum genome (Additional file 1: Table S9). The genes sharing homology with known transporters were obtained using BlastP and each of them appeared in several BlastP results; their expressions were determined by the qRT-PCR assay (Additional file 1: Table S10). Nitrogen metabolism and sugar metabolism in fungal cells would be active after nitrogen and carbohydrates were assimilated or transferred. Therefore, genes encoding glutamate synthase and nitrite reductase (Additional file 1: Table S8), and that encoding the Golgi GDP-mannose transporter, which is functional in the glycosylation of secreted proteins  and five enzymes involved in galactose and nucleotide sugar metabolism (Additional file 1: Table S10) were also included. The gene encoding the tetracycline resistance protein was chosen as a control, because it recognizes and exports tetracycline from the cell [77, 78] and is expected to show no expression change under the experimental conditions in the present study.
Differences in expression levels of genes involved in symbiosis under co-culture conditions
From the results presented in Figure 7, the expression levels of most genes increased only under the condition for experimental group IV (Table 5) (weight ratio of lichen-forming fungus and phycobiont is 10:3, and culture time is 72 hours), which indicated that the contact time and biomass ratio of both symbionts affects their recognition and nutrient transfer significantly.
Three lectin genes (F481_01961, F481_04092, F481_02882) showed significant increases (log2 fold change > 2) in their expression level in experiment group IV compared with the control group B (only the mycobiont was cultured in BBM for 72 h) (Table 5; Figure 7). Interestingly, the common aspect of these lectins is that there are no predicted signal peptides in their amino acid sequences (Additional file 1: Table S7) (http://elm.eu.org) . Lectins are sugar-binding proteins that play an important roles in cell recognition . The first lectin discovered in lichens was that found in Peltigera canina and Peltigera polydactyla. Thereafter, some studies showed that lectins secreted by lichen-forming fungi may be involved in recognizing their compatible algae [15, 16, 81]. Lectins, as surface proteins, can directly contact compatible phycobionts, and the contact requires receptor sites on the surface of the phycobiont to which the lectins binds [82, 83]. Recent studies on lectins in P. membranacea indicated that the expression of galectin lec-1 was influenced by the presence of the phycobiont , and the evolution of galectin lec-2 was driven by interaction with different strains of the phycobiont Nostoc. Galectins, together with cellular slime mold lectin discoidin I, do not have signal peptides in eukaryotic cells and both of them have intracellular and extracellular localizations [86, 87]. Our study showed that lectins without signal peptides could play a major role in symbiotic recognition and provided a practical guide for screening for lectins participating in the interaction between lichen-forming fungi and their photosynthetic partners.
An ammonium transporter (F481_01640) and a nitrate transporter (F481_07695) were up-regulated (log2 fold change > 2) in experimental group IV (Figure 7). BBM contains nitrate ions; therefore, the up-regulation of the nitrate transporter (F481_07695) is expected. However, the up-regulation of the ammonium transporter (F481_01640) suggests that ammonium could be produced from nitrate in E. pusillum and transferred to its photosynthetic partner, because there is no ammonium ion in BBM. The differences in the expressions of glutamate synthase (F481_06302) and nitrite reductase (F481_07075) indicate that nitrogen metabolism is altered in the lichen-forming fungus when it meets its compatible alga. It has been reported that algae from different lichens can utilize organic and inorganic nitrogen [62–64]. Thus, we can conclude that the pattern of nitrogen transfer in E. pusillum is that the mycobiont absorbs various nitrogen sources from the environment and can convert them to different forms in vivo to meet the algae’s demand for nitrogen, based on the expansion of nitrogen transporter families in E. pusillum genome and the result of co-culture experiments.
The lichen-forming fungus cannot obtain any organic nutrition from BBM. However, six of the symbiosis-related sugar transporter genes (F481_00818, F481_07878, F481_08427, F481_05999, F481_06586, and F481_02699) show significant differences in expression in experiment group IV (Figure 7). This suggests that some carbohydrates are released into the BBM by D. chodatii, and the transporters encoded by these genes may be responsible for transferring these carbohydrates into the fungal cells. The upregulation of the expression of the Golgi GDP-mannose transporter (F481_08390) under experimental group IV conditions, bearing in mind that a previous study proved that this transporter participates in glycosylation , indicated that the sugar metabolic pathway is active in E. pusillum after obtaining carbohydrates. Furthermore, the enzymes (F481_03176, F481_07291, F481_06583, and F481_00358) involved in galactose and nucleotide sugar metabolism also showed high expressions levels, which demonstrated that the carbohydrates produced by the alga and absorbed by the fungus are indeed used for fungal cell life activities.
Native and non-native sugars utilized by E. pusillum
The carbohydrates used for growth and metabolism of lichen-forming fungus originate from the photosynthetic products of its photosynthetic partner, however, the forms of the photosynthetic products absorbed by different lichen-forming fungi vary depending on their different photobionts. It has been report that lichen-forming fungi absorb many polyols or glucose from algae or cyanobacteria, respectively [20, 21], and the photosynthetic product transferred from D. chodatii to E. pusillum was sorbitol . Additionally, previous studies showed that some monosaccharides, such as glucose, and disaccharides, such as trehalose and sucrose, were transferred between non-lichenized symbionts [68–75]. Because some homology structures are frequently found in different sugar transporters from diverse fungi, the specificities of these transporters may be low. Therefore, 11 different carbohydrates were used to confirm the function of sugar transporters in E. pusillum and the potential carbohydrates transferred into this fungus.
The most significant up-regulation in expression of the examined genes was in the sample with trehalose for 24 hours. However, the high expression levels were not maintained over the next 48 hours (Figure 8). This suggests that trehalose cannot be absorbed by the lichen-forming fungus as an energy source, but is likely to be a signal molecule for the activation of the gene expression in E. pusillum. Trehalose is involved in many functions besides osmoprotection , and acts as a signal molecule that is possibly exported from bacteria to plants to regulate the carbon and nitrogen metabolism of plants in plant-bacteria interactions . Trehalose can be converted into trehalose-6-phosphate (T6P), which inhibits hexokinase activity and regulates glycolysis. Over expression of the trehalose 6-phospate synthase (TPS) gene induces the expression of several genes involved in stress tolerance, and carbon and nitrogen metabolism [90, 91]. In the alga-invertebrate association, trehalose produced by the alga was shown to be present in certain symbiotic interfaces . Although no research has shown a link between trehalose and lichen symbiosis interaction, our results suggest that trehalose may act as a signal molecule and would play an important role in the symbiotic interaction in E. pusillum.
Four carbohydrates (glucose, sucrose, sorbitol, and mannitol) increased the transcription levels of sugar transporters or enzymes involved in sugar metabolism in E. pusillum at 24 and 72 hours (Figure 8). These results suggest that these carbohydrates maintain the metabolism of the lichen-forming fungus over a long period, which implies that E. pusillum may absorb them as its carbon sources through some common sugar transporters. Although the transcription levels of some genes slightly fluctuated at 24 and 72 hours in the presence other sugars (arabinose, erythritol, fructose, galactose, ribitol and ribose) in the BBM, no genes exhibited a clear rising trend compared with the control group, especially those involved in sugar metabolism (Figure 8). This indicates that these sugars cannot be utilized as carbon sources by E. pusillum.
The transcription levels of some sugar transporters, such as F481_05999, F481_06586, and F481_2699, were up-regulated when various sugars were added to the BBM, indicating that they can transfer more than one carbohydrate into E. pusillum. Some genes showed different changes in their transcription levels, either in trend or extent, compared with those in the co-culture experiment. For example, sugar transporter F481_02469 showed no differential expression in the co-culture experiment, but was up-regulated slightly under some sugars treatments. In addition, the transcription levels of sugar transporters F481_05999, F481_06586, and F481_02699 were up-regulated more obviously in the co-culture experiment (Figures 7 and 8). These phenomena imply that there is some induction mechanism between the symbionts; i.e., the expressions of certain genes in a lichen-forming fungus are influenced by its photosynthetic partner to absorb organic carbon effectively.
Ion chromatography was used to detect glucose, sucrose, sorbitol and mannitol in the filtrates of experiment group IV, and controls B and C, to determine whether they could be released by D. chodatii and to determine the conditions under which this alga secretes organic carbon into the BBM. Sorbitol (12.6 mg/100 mL), glucose (0.8 mg/100 mL) and sucrose (0.8 mg/100 mL) have been detected in group IV (Figure 9). Thus, sorbitol showed the highest accumulation; however, there were other, unidentified, carbohydrates, which suggested that the lichen-forming fungus E. pusillum preferentially absorbs sorbitol. By contrast, no mannitol was found in this sample. There were no detectable carbohydrates in the control group B and C (the phycobiont and the lichen-forming fungus cultured respectively in BBM), which demonstrated that the pure cultured D. chodatii does not release carbon. Previous experiments for the phycobionts of Ramalina crassa and Ramalina subbreviuscula indicated that the pure cultured phycobiont released ribitol to the medium . However, Hawksworth (1984) suggested that the photobiont had lost its ability to release carbohydrate after isolation from the lichen thallus, and that the lichen-forming fungus exerted some specific control on the photobiont cells that lead them to secrete carbohydrates . Our result implies that E. pusillum can control the sugar export of D. chodatii, and that the biomass ratio and contact time are crucial for the interaction between the two symbionts.
Although mannitol is not released by D. chodatii, it can still be utilized by E. pusillum; it is likely to be because mannitol and sorbitol are isomers. This suggests that E. pusillum has the potential to use non-natural carbon sources. This phenomenon is likely to occur in other lichen-forming fungi, which would contribute to discovering the mechanism of the algal switch [93–95] and could provide clues for the evolution of lichenization.
Approximately 40% of all Ascomycota are lichen-forming; thus, lichenization is regarded as one of the most important fungal lifestyles . Hence, genomic information for lichen-forming fungi would expand the knowledge of fungi. In the present study, we report, for the first time, the characteristic of the lichen-forming fungal genome, which displays many features that are different to other fungi.
This is the first study to report that the lichen-forming fungal genome have undergone RIP. Genes for mating system, secondary metabolism, and the drought-related mechanisms were indentified in E. pusillum genome, which are worth being investigated in the future. The evolution analysis of multigene families indicated the expansion and contraction in E. pusillum genome reveal the effect of lichenization on lichen-forming fungi.
Co-culture experiments suggest that the lectins without signal peptides would be likely to play an essential role in the recognition of lichen symbiosis, and one of the most striking findings in these experiments is that an appropriate weigh ratio of lichen-forming fungus and its photosynthetic partner and sufficient contact time are vital for their recognition and mutual influence. We also confirmed that the most important natural carbon source for E. pusillum is sorbitol transferred from D. chodatii; however, this lichen-forming fungus can also use other non-natural carbohydrates under the pure culture condition.
A mycobiont-phycobiont interplay model is shown as Figure 10. The model reflects aspects of the recognition and interaction of the lichen thallus in E. pusillum and is likely to be applicable to other lichens, especially those whose photobionts are algae. This study provides a valuable genomic resource for future research in screening functional genes including drought-tolerance genes from lichens and would be useful for investigating the formation and divergence on the functional biology between lichenized and nonlichenized fungi.
The lichen-forming fungal strain of E. pusillum Z07020 (HMAS-L-300199) was isolated by a single-spore discharge from the perithecium of lichen E. pusillum collected from Shapotou Desert Research Station (SDRS) of the Chinese Academy of Sciences (CAS) in the Tengger Desert of northern China . The isolates were grown on 1.5% water agar for 1–2 weeks, and then cultured at room temperature after transfer to potato dextrose liquid medium.
Genome sequencing and assembly
The genome of the lichen-forming fungus E. pusillum was sequenced using high-throughput next-generation sequencing technology and the sequencing platforms were Roche 454 and Illumina Solexa systems. Genomic libraries containing 8-kb inserts were constructed and 1,394,086 paired-end reads (281.9 Mb) were generated using the 454 Roche GS FLX system. The Illumina adaptors were ligated onto the genomic DNA fragments, and DNA fragments with estimated sizes of 0.5 kb to 3 kb were selected using gel-electrophoresis. Libraries were PCR-amplified using Phusion polymerase. Sequencing libraries were denatured with sodium hydroxide and diluted in hybridization buffer for loading onto a single lane of an Illumina GA flow cell. A Solexa sequencer generated the mate-paired reads (7,155,072 reads, 716 Mb) and paired-end reads (18,176,986 reads, 1818 Mb). Solexa sequencing paired-end reads and mate-paired reads were assembled by SOAPdenovo , which adopts the de Bruijn graph data structure to construct contigs.
Gene prediction and annotation
We used Augustus , GeneID , and GeneMark-ES  programs to predict the gene models for E. pusillum. A final set of gene models was selected by EvidenceModeler . An ab initio prediction was carried on using the annotated information of A. fumigatus as a reference. All predicted gene models were subjected to GO , KOG , FunCat  and Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis . Protein domains were predicted using InterProScan  against various domain libraries (HMMPfam, superfamily, HMMTigr, and HMMSmart). Repetitive elements were identified by blasting against the RepeatMasker library (http://www.repeatmasker.org/). Non-coding RNAs were predicted according to the Rfam database , and tRNAs were predicted using tRNAscan-SE . Pseudogenes and rRNAs were designated using PseudoGene  and RNAmmer  respectively.
Orthologous gene and phylogeny analysis of E. pusillum
The sequences of corresponding orthologous genes from 15 fungi were aligned using ClustalW . A maximum parsimony (MP) phylogenomic tree was created using the concatenated amino acid sequences in phylogenetic analysis using parsimony (PAUP) , and a bootstrap analysis with 1000 replications was performed to evaluate the reliability of the phylogenetic tree. The divergence time between species was estimated using the r8s method . The time of divergence between Ascomycota and Basidiomycota is set as 500 Myr, and the time between Pezizomycotina and Saccharomycotina was 350 Myr .
Multigene families and evolution analysis
The families that were absent from the most recent common ancestor were chosen to analyze the evolution of protein families . Multigene families were identified using the MCL method . Whole genome blast analyses against Transporter Classification Database (http://www.tcdb.org/tcdb/) and GPCRDB (http://www.gpcr.org/7tm/) database were performed to identify genes exhibiting difference between E. pusillum and other fungi. All proteins in the genome were blast searched against the database with an e-value cutoff of ≤ e-10 and at least 40% identity over 60% coverage. A new software, BadiRate , was used in the evolutionary analysis of multigene families to estimate rates of gene gain and loss. Families with p-values less than 0.01 were considered to have experienced significant expansion or contraction.
To detect the sexual cycle of E. pusillum, genes involved in the mating process, incompatibility, ascomata and conidiophore development, and HET were identified using BlastP against related genes from A. niger, A. nidulans , N. crassa  and P. anserine [120, 121].
Analysis of genes involved in secondary metabolism
Genes encoding PKS, NRPS, and NRPS-PKS hybrid genes in the genome of E. pusillum were analyzed with the program SMURF (http://www.jcvi.org/smurf/index.php). Modulation analysis and domain extraction of different NRPS or PKS proteins were conducted by Blast searching against the SBSPKS database .
The potential secreted proteins of E. pusillum and other fungi, including F. graminearum, M. oryzae, and M. anisopliae, were predicted by SignaIP 3.0 analysis using a hidden Markov model [http://www.cbs.dtu.dk/services/SignalP/].
RNA was extracted from cultured lichen-forming fungus E. pusillum. cDNA synthesis and relative quantitative RT-PCR were carried out as described previously . For each treatment, qRT-PCR was performed in an Applied Biosystems 7500 Real-Time PCR system (Applied Biosystems, USA). The data were analyzed using the 2-ΔCt method.
Gene-stability measure of reference genes using geNORM
Vandesomepele et al. (2002) developed an algorithm named geNORM that determines the expression stability of reference genes . The calculated gene-stability measure (M) relies on the principle that the expression ratio of two ideal internal reference genes is identical in all samples, regardless of the experimental condition or cell type. Ten genes whose expressions showed no difference in comparative transcriptome data between control and drought-stress conditions (unpublished data) were chosen to determine the best reference genes (Additional file 1: Table S11). Ultimately, we chose the tetracycline resistance protein (F481_ 05245) gene as the reference gene.
Treatment of samples for co-culture experiments
Four treatments were carried out, comprising two weight ratios (10:1 and 10:3) for the lichen-forming fungus and phycobiont and two culture times (24 h and 72 h) (Table 5), to investigate whether the ratio and contact time between the lichen-forming fungus and phycobiont has an effect on gene expression in E. pusillum when both symbionts are incubated together.
The experiment was performed in Bold’s basal medium (BBM) (does not contain carbohydrate), and under 12 h illuminations per 24 h, which allow the algal cells to produce carbohydrates by photosynthesis.
Treatment of samples for sugar transfer experiments
To confirm the capability of E. pusillum to utilize carbohydrates, polyols (mannitol, sorbitol, ribitol, and erythritol), monosaccharides (glucose, fructose, arabinose, ribose, and galactose) and disaccharides (trehalose and sucrose), respectively, were added to BBM containing only the lichen-forming fungus, and the transcript levels of genes involved in sugar transport and metabolism were determined after culturing for 24 and 72 hrs.
Determination of photosynthetic products by ion chromatography
Ion chromatography was used to determine glucose, sucrose, sorbitol and mannitol in co-cultured samples. The analytes were separated on a CarboPac™ PA1 (4 mm × 250 mm) anion exchange column using 200 mmol/L NaOH as mobile phase at flow rate of 1.0 mL/min and detected with a pulsed amperometric detector.
Availability of supporting data
The data sets supporting the results of this article are included within the article (and its additional files).
Bold’s basal medium
G-protein coupled receptor
Mycological herbarium of microbiology institute,the chinese academy of sciences
Kyoto encyclopedia of genes and genomes database
Eukaryotic clusters of orthologous groups
Large-subunit ribosomal DNA
National center for biotechnology information
Major facilitator superfamily
Non-ribosomal peptide synthetase
Polymerase chain reaction
Quantitative real-time PCR
Repeat-induced point mutation
Scanning electron microscopy
Small-subunit ribosomal DNA
Trehalose 6-phospate synthase
Small secret protein.
Ahmadjian V: The lichen symbiosis. 1993, New York: John Wiley
Kirk PMCP, Minter DW, Stalpers JA: Dictionary of the fungi. 2008, Wallingford: CABI, 10
Nash TH: Lichen biology. 2nd edn. Cambridge. 2008, New York: Cambridge University Press
Ingolfsdottir K: Usnic acid. Phytochemistry. 2002, 61: 729-736. 10.1016/S0031-9422(02)00383-7.
Burkholder PR, Evans AW, McVeigh I, Thornton HK: Antibiotic activity of lichens. P Natl Acad Sci USA. 1944, 30: 250-10.1073/pnas.30.9.250.
Honegger R: The lichen symbiosis—what is so spectacular about it?. The Lichenologist. 1998, 30: 193-212.
Wedin M, Döring H, Gilenstam G: Saprotrophy and lichenization as options for the same fungal species on different substrata: environmental plasticity and fungal lifestyles in the Stictis–Conotrema complex. New Phytol. 2004, 164: 459-465. 10.1111/j.1469-8137.2004.01198.x.
Gargas A, DePriest PT, Grube M, Tehler A: Multiple origins of lichen symbioses in fungi suggested by SSU rDNA phylogeny. Science-New York then Washington. 1995, 1492-1492.
Lutzoni F, Pagel M, Reeb V: Major fungal lineages are derived from lichen symbiotic ancestors. Nature. 2001, 411: 937-940. 10.1038/35082053.
Rikkinen J: Molecular studies on cyanobacterial diversity in lichen symbioses. Lichens: from genome to ecosystems in a changing world MycoKeys. 2013, 6: 3-32.
Ahmadjian V, Jacobs JB, Russell LA: Scanning electron-microscope study of early lichen synthesis. Science. 1978, 200: 1062-1064. 10.1126/science.200.4345.1062.
Ahmadjian VHH: The culture and synthesis of Endocarpon pusillum and Staurothele clopima. Lichenologist. 1970, 4: 13-
Trembley ML, Ringli C, Honegger R: Morphological and molecular analysis of early stages in the resynthesis of the lichen Baeomyces rufus. Mycol Res. 2002, 106: 768-776. 10.1017/S0953756202006081.
Guillot J, Konska G: Lectins in higher fungi. Biochem Syst Ecol. 1997, 25: 203-230. 10.1016/S0305-1978(96)00110-X.
Kardish NSL, Flemminger N, Galun M: Lectin from the lichen Nephroma laevigatum. Localization and function. Symbiosis. 1991, 11: 47-62.
Andrews GA, Chavey PS, Smith JE: Reactivity of lichen lectins with blood typed canine erythrocytes. Res Vet Sci. 1992, 53: 315-319. 10.1016/0034-5288(92)90132-L.
Lockhart C, Rowell P, Stewart W: Phytohaemagglutinins from the nitrogen-fixing lichens Peltigera canina and P. polydactyla. FEMS Microbiology Letters. 1978, 3: 127-130. 10.1111/j.1574-6968.1978.tb01899.x.
Richardson DHS, Hill DJ, Smith DC: Lichen physiology. New Phytol. 1968, 67: 469-486. 10.1111/j.1469-8137.1968.tb05476.x.
Komiya T, Shibata S: Polyols produced by the cultured phyco-and mycobionts of some Ramalina species. Phytochemistry. 1971, 10: 695-699. 10.1016/S0031-9422(00)97135-8.
Hill DJ, Ahmadjia V: Relationship between carbohydrate movement and symbiosis in lichens with Green-Algae. Planta. 1972, 103: 267-10.1007/BF00386850.
Hawksworth DL, Hill DJ: The lichen-forming fungi. Blackie. 1984, 5: 62-69.
Gauslaa Y, Solhaug KA: Fungal melanins as a sun screen for symbiotic green algae in the lichen Lobaria pulmonaria. Oecologia. 2001, 126: 462-471. 10.1007/s004420000541.
Joneson S, Armaleo D, Lutzoni F: Fungal and algal gene expression in early developmental stages of lichen-symbiosis. Mycologia. 2011, 103: 291-306. 10.3852/10-064.
Suzanne Joneson FL: Compatibility and thigmotropism in the lichen symbiosis: a reappraisal. Symbiosis. 2009, 47: 109-115. 10.1007/BF03182294.
Awasthi DD: A hand book of lichens. 2000, Bishen Singh Mahendra Pal Singh: Dehradun
Zhang T, Wei J: Survival analyses of symbionts isolated from Endocarpon pusillum Hedwig to desiccation and starvation stress. Sci China Life Sci. 2011, 54: 480-489. 10.1007/s11427-011-4164-z.
Xavier BB, Miao VPW, Jonsson ZO, Andresson OS: Mitochondrial genomes from the lichenized fungi Peltigera membranacea and Peltigera malacea: Features and phylogeny. Fungal Biol-Uk. 2012, 116: 802-814. 10.1016/j.funbio.2012.04.013.
Junttila SM, Rudd S: 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. BMC Genomics. 2012, 13: 575-10.1186/1471-2164-13-575.
Wang YY, Zhang T, Zhou QM, Wei JC: Construction and characterization of a full-length cDNA library from mycobiont of Endocarpon pusillum (lichen-forming Ascomycota). World J Microb Biot. 2011, 27: 2873-2884. 10.1007/s11274-011-0768-5.
Ostlund G, Schmitt T, Forslund K, Kostler T, Messina DN, Roopra S, Frings O, Sonnhammer ELL: InParanoid 7: new algorithms and tools for eukaryotic orthology analysis. Nucleic Acids Res. 2010, 38: D196-D203. 10.1093/nar/gkp931.
Cambareri EB, Jensen BC, Schabtach E, Selker EU: Repeat-Induced G-C to a-T Mutations in Neurospora. Science. 1989, 244: 1571-1575. 10.1126/science.2544994.
Margolin BS, Garrett-Engele PW, Stevens JN, Fritz DY, Garrett-Engele C, Metzenberg RL, Selker EU: A methylated Neurospora 5S rRNA pseudogene contains a transposable element inactivated by repeat-induced point mutation. Genetics. 1998, 149: 1787-1797.
Selker EU, Tountas NA, Cross SH, Margolin BS, Murphy JG, Bird AP, Freitag M: The methylated component of the Neurospora crassa genome. Nature. 2003, 422: 893-897. 10.1038/nature01564.
Hane JK, Oliver RP: RIPCAL: a tool for alignment-based analysis of repeat-induced point mutations in fungal genomic sequences. Bmc Bioinformatics. 2008, 9: 478-490. 10.1186/1471-2105-9-478.
Hogenhout SA, Van der Hoorn RAL, Terauchi R, Kamoun S: Emerging concepts in effector biology of plant-associated organisms. Mol Plant Microbe In. 2009, 22: 115-122. 10.1094/MPMI-22-2-0115.
Plett JM, Kemppainen M, Kale SD, Kohler A, Legue V, Brun A, Tyler BM, Pardo AG, Martin F: A secreted effector protein of laccaria bicolor is required for symbiosis development. Curr Biol. 2011, 21: 1197-1203. 10.1016/j.cub.2011.05.033.
Kloppholz S, Kuhn H, Requena N: A secreted fungal effector of glomus intraradices promotes symbiotic biotrophy. Curr Biol. 2011, 21: 1204-1209. 10.1016/j.cub.2011.06.044.
Deller S, Hammond-Kosack KE, Rudd JJ: The complex interactions between host immunity and non-biotrophic fungal pathogens of wheat leaves. J plant physiol. 2011, 168: 63-71. 10.1016/j.jplph.2010.05.024.
Armaleo D, Sun XM, Culberson C: Insights from the first putative biosynthetic gene cluster for a lichen depside and depsidone. Mycologia. 2011, 103: 741-754. 10.3852/10-335.
Wang Y, Kim JA, Cheong YH, Joshi Y, Koh YJ, Hur JS: Isolation and characterization of a reducing polyketide synthase gene from the lichen-forming fungus usnea longissima. J Microbiol. 2011, 49: 473-480. 10.1007/s12275-011-0362-4.
Kim JA, Hong SG, Cheong YH, Koh YJ, Hur JS: A new reducing polyketide synthase gene from the lichen-forming fungus Cladonia metacorallifera. Mycologia. 2012, 104: 362-370. 10.3852/11-001.
Wang Y, Kim JA, Cheong YH, Koh YJ, Hur JS: Isolation and characterization of a non-reducing polyketide synthase gene from the lichen-forming fungus Usnea longissima. Mycol Prog. 2012, 11: 75-83. 10.1007/s11557-010-0730-1.
Kranner I, Beckett R, Hochman A, Nash TH: Desiccation-tolerance in lichens: a review. Bryologist. 2008, 111: 576-593. 10.1639/0007-2745-111.4.576.
Gilman AG: G-Proteins - Transducers of receptor-generated signals. Annu Rev Biochem. 1987, 56: 615-649. 10.1146/annurev.bi.56.070187.003151.
Glass NL, Jacobson DJ, Shiu PK: The genetics of hyphal fusion and vegetative incompatibility in filamentous ascomycete fungi. Annu Rev Genet. 2000, 34: 165-186. 10.1146/annurev.genet.34.1.165.
Saupe SJ, Clave C, Begueret J: Vegetative incompatibility in filamentous fungi: Podospora and Neurospora provide some clues. Curr Opin Microbiol. 2000, 3: 608-612. 10.1016/S1369-5274(00)00148-X.
Fedorova ND, Badger JH, Robson GD, Wortman JR, Nierman WC: Comparative analysis of programmed cell death pathways in filamentous fungi. BMC Genomics. 2005, 6: 177-191. 10.1186/1471-2164-6-177.
Saupe SJ: Molecular genetics of heterokaryon incompatibility in filamentous ascomycetes. Microbiol Mol Biol R. 2000, 64: 489-10.1128/MMBR.64.3.489-502.2000.
Haugen P, Simon DM, Bhattacharya D: The natural history of group I introns. Trends Genet. 2005, 21: 111-119. 10.1016/j.tig.2004.12.007.
Gargas A, Taylor JW: Phylogeny of discomycetes and early radiations of the apothecial ascomycotina inferred from Ssu Rdna sequence data. Exp Mycol. 1995, 19: 7-15. 10.1006/emyc.1995.1002.
Bhattacharya D, Friedl T, Helms G: Vertical evolution and intragenic spread of lichen-fungal group I introns. J Mol Evol. 2002, 55: 74-84. 10.1007/s00239-001-2305-x.
Gargas A, DePriest PT, Taylor JW: Positions of multiple insertions in SSU rDNA of lichen-forming fungi. Mol Biol Evol. 1995, 12: 208-218.
Hibbett DS: Phylogenetic evidence for horizontal transmission of group I introns in the nuclear ribosomal DNA of mushroom-forming fungi. Mol Biol Evol. 1996, 13: 903-917. 10.1093/oxfordjournals.molbev.a025658.
Depriest PT: Molecular innovations in lichen systematics - the use of ribosomal and intron nucleotide-sequences in the cladonia-chlorophaea complex. Bryologist. 1993, 96: 314-325. 10.2307/3243859.
Smith RL, Banks JL, Snavely MD, Maguire ME: Sequence and topology of the CorA magnesium transport systems of Salmonella typhimurium and Escherichia coli. Identification of a new class of transport protein. J Biol Chem. 1993, 268: 14071-14080.
Kehres DG, Maguire ME: Structure, properties and regulation of magnesium transport proteins. Biometals. 2002, 15: 261-270. 10.1023/A:1016078832697.
Liu GJ, Martin DK, Gardner RC, Ryan PR: Large Mg2 + −dependent currents are associated with the increased expression of ALR1 in Saccharomyces cerevisiae. FEMS Microbiology Letters. 2002, 213: 231-237. 10.1111/j.1574-6968.2002.tb11311.x.
Gardner RC: Genes for magnesium transport. Curr Opin Plant Biol. 2003, 6: 263-267. 10.1016/S1369-5266(03)00032-3.
Zhao H, Zhou Q, Zhou M, Li C, Gong X, Liu C, Qu C, Si W, Hong F: Magnesium deficiency results in damage of nitrogen and carbon cross-talk of maize and improvement by cerium addition. Biol Trace Elem Res. 2012, 148: 102-109. 10.1007/s12011-012-9340-x.
Mcswain BD, Tsujimoto HY, Arnon DI: Photochemical activity and components of membrane preparations from blue-green-algae.4. Effects of magnesium and chloride-Ions on light-induced electron-transport in membrane fragments from a blue-green-alga. Biochimica Et Biophysica Acta. 1976, 423: 313-322. 10.1016/0005-2728(76)90188-2.
Smith D, Muscatine L, Lewis D: Carbohydrate movement from autotrophs to heterotrophs in parasitic and mutualistic symbiosis. Biol Rev Camb Philos Soc. 1969, 44: 17-90. 10.1111/j.1469-185X.1969.tb00821.x.
Crittenden PD: The effect of oxygen deprivation on inorganic nitrogen uptake in an Antarctic macrolichen. Lichenologist. 1996, 28: 347-354.
Crittenden PD: Nutrient exchange in an Antarctic macrolichen during summer snowfall snow melt events. New Phytol. 1998, 139: 697-707. 10.1046/j.1469-8137.1998.00236.x.
Dahlman L, Persson J, Palmqvist K, Nasholm T: Organic and inorganic nitrogen uptake in lichens. Planta. 2004, 219: 459-467.
Berman T, Chava S: Algal growth on organic compounds as nitrogen sources. J Plankton Res. 1999, 21: 1423-1437. 10.1093/plankt/21.8.1423.
McDonald TR, Dietrich FS, Lutzoni F: Multiple horizontal gene transfers of ammonium transporters/ammonia permeases from prokaryotes to eukaryotes: toward a new functional and evolutionary classification. Mol Biol Evol. 2012, 29: 51-60. 10.1093/molbev/msr123.
McDonald TR, Mueller O, Dietrich FS, Lutzoni F: High-throughput genome sequencing of lichenizing fungi to assess gene loss in the ammonium transporter/ammonia permease gene family. Bmc Genomics. 2013, 14: 225-10.1186/1471-2164-14-225.
Wahl R, Wippel K, Goos S, Kämper J, Sauer N: A novel high-affinity sucrose transporter is required for virulence of the plant pathogen Ustilago maydis. Plos Biol. 2010, 8: e1000303-10.1371/journal.pbio.1000303.
Fang WG: St Leger RJ: Mrt, a gene unique to fungi, encodes an oligosaccharide transporter and facilitates rhizosphere competency in metarhizium robertsii. Plant Physiology. 2010, 154: 1549-1557. 10.1104/pp.110.163014.
Helber N, Wippel K, Sauer N, Schaarschmidt S, Hause B, Requena N: A versatile monosaccharide transporter that operates in the arbuscular mycorrhizal fungus glomus sp is crucial for the symbiotic relationship with plants. Plant Cell. 2011, 23: 3812-3823. 10.1105/tpc.111.089813.
Voegele RT, Struck C, Hahn M, Mendgen K: The role of haustoria in sugar supply during infection of broad bean by the rust fungus Uromyces fabae. P Natl Acad Sci USA. 2001, 98: 8133-8138. 10.1073/pnas.131186798.
Doehlemann G, Molitor F, Hahn M: Molecular and functional characterization of a fructose specific transporter from the gray mold fungus Botrytis cinerea. Fungal Genet Biol. 2005, 42: 601-610. 10.1016/j.fgb.2005.03.001.
Wiese J, Kleber R, Hampp R, Nehls U: Functional characterization of the Amanita muscaria monosaccharide transporter, Am Mst1. Plant Biology. 2000, 2: 278-282. 10.1055/s-2000-12984.
Polidori E, Ceccaroli P, Saltarelli R, Guescini M, Menotta M, Agostini D, Palma F, Stocchi V: Hexose uptake in the plant symbiotic ascomycete Tuber borchii Vittadini: biochemical features and expression pattern of the transporter TBHXT1. Fungal Genet Biol. 2007, 44: 187-198. 10.1016/j.fgb.2006.08.001.
Schussler A, Martin H, Cohen D, Fitz M, Wipf D: Characterization of a carbohydrate transporter from symbiotic glomeromycotan fungi. Nature. 2006, 444: 933-936. 10.1038/nature05364.
Abe M, Noda Y, Adachi H, Yoda K: Localization of GDP-mannose transporter in the Golgi requires retrieval to the endoplasmic reticulum depending on its cytoplasmic tail and coatomer. J Cell Sci. 2004, 117: 5687-5696. 10.1242/jcs.01491.
Schnappinger D, Hillen W: Tetracyclines: antibiotic action, uptake, and resistance mechanisms. Arch Microbiol. 1996, 165: 359-369. 10.1007/s002030050339.
Walsh C: Molecular mechanisms that confer antibacterial drug resistance. Nature. 2000, 406: 775-781. 10.1038/35021219.
Gould CM, Diella F, Via A, Puntervoll P, Gemund C, Chabanis-Davidson S, Michael S, Sayadi A, Bryne JC, Chica C, et al: ELM: the status of the 2010 eukaryotic linear motif resource. Nucleic Acids Res. 2010, 38: D167-D180. 10.1093/nar/gkp1016.
Sharon N, Lis H: Lectins as cell recognition molecules. Science. 1989, 246: 227-234. 10.1126/science.2552581.
Sacristan M, Millanes AM, Legaz ME, Vicente C: A lichen lectin specifically binds to the alpha-1,4-Polygalactoside moiety of urease located in the cell wall of homologous algae. Plant Signal Behav. 2006, 1: 23-27. 10.4161/psb.1.1.2276.
Galun MKN: Lectins as determinants of symbiotic specificity in lichens. Cryptog Bot. 1995, 5: 144-148.
Vivas M, Sacristan M, Legaz ME, Vicente C: The cell recognition model in chlorolichens involving a fungal lectin binding to an algal ligand can be extended to cyanolichens. Plant Biol (Stuttg). 2010, 12: 615-621.
Miao VPW, Manoharan SS, Snaebjarnarson V, Andresson OS: Expression of lec-1, a mycobiont gene encoding a galectin-like protein in the lichen Peltigera membranacea. Symbiosis. 2012, 57: 23-31. 10.1007/s13199-012-0175-1.
Manoharan SS, Miao VP, Andresson OS: LEC-2, a highly variable lectin in the lichen. Symbiosis. 2012, 58: 91-98. 10.1007/s13199-012-0206-y.
Liener IE, Sharon N, Goldstein IJ: The Lectins: properties, functions, and applications in biology and medicine. 1986, Orlando: Academic
Barondes SH, Cooper DN, Gitt MA, Leffler H: Galectins. Structure and function of a large family of animal lectins. J Biol Chem. 1994, 269: 20807-20810.
Iturriaga G, Suarez R, Nova-Franco B: Trehalose metabolism: from osmoprotection to signaling. Int J Mol Sci. 2009, 10: 3793-3810. 10.3390/ijms10093793.
Suarez R, Wong A, Ramirez M, Barraza A, Orozco Mdel C, Cevallos MA, Lara M, Hernandez G, Iturriaga G: Improvement of drought tolerance and grain yield in common bean by overexpressing trehalose-6-phosphate synthase in rhizobia. Mol Plant Microbe Interact. 2008, 21: 958-966. 10.1094/MPMI-21-7-0958.
Blazquez MA, Lagunas R, Gancedo C, Gancedo JM: Trehalose-6-Phosphate, a new regulator of yeast glycolysis that inhibits hexokinases. Febs Letters. 1993, 329: 51-54. 10.1016/0014-5793(93)80191-V.
Muller J, Boller T, Wiemken A: Trehalose and trehalase in plants: recent developments. Plant Science. 1995, 112: 1-9. 10.1016/0168-9452(95)04218-J.
Pardy RL, Spargo B, Crowe JH: Release of trehalose by symbiotic algae. Symbiosis. 1989, 7: 149-158.
Nelsen MP, Gargas A: Dissociation and horizontal transmission of codispersing lichen symbionts in the genus Lepraria (Lecanorales: Stereocaulaceae). New Phytol. 2008, 177: 264-275.
Piercey-Normore MD, DePriest PT: Algal switching among lichen symbioses. American Journal of Botany. 2001, 88: 1490-1498. 10.2307/3558457.
Piercey-Normore MD: The lichen-forming ascomycete Evernia mesomorpha associates with multiple genotypes of Trebouxia jamesii. New Phytol. 2006, 169: 331-344. 10.1111/j.1469-8137.2005.01576.x.
Mats Wedin HDGG: Saprotrophy and lichenization as options for the same fungal species on different substrata: environmental plasticity and fungal lifestyles in the Stictis–Conotrema complex. New Phytol. 2004, 164: 7-10.1111/j.1469-8137.2004.01180.x.
Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y: SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012, 1: 1-6. 10.1186/2047-217X-1-1.
Stanke M, Morgenstern B: AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints. Nucleic Acids Res. 2005, 33: W465-W467. 10.1093/nar/gki458.
Parra G, Blanco E, Guigó R: Geneid in drosophila. Genome research. 2000, 10: 511-515. 10.1101/gr.10.4.511.
Ter-Hovhannisyan V, Lomsadze A, Chernoff YO, Borodovsky M: Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome research. 2008, 18: 1979-1990. 10.1101/gr.081612.108.
Haas BJ, Salzberg SL, Zhu W, Pertea M, Allen JE, Orvis J, White O, Buell CR, Wortman JR: Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008, 9: R7-10.1186/gb-2008-9-1-r7.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Sonnhammer ELL, Koonin EV: Orthology, paralogy and proposed classification for paralog subtypes. Trends Genet. 2002, 18: 619-620. 10.1016/S0168-9525(02)02793-2.
Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, Tetko I, Guldener U, Mannhaupt G, Munsterkotter M, Mewes HW: The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004, 32: 5539-5545. 10.1093/nar/gkh894.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32: D277-D280. 10.1093/nar/gkh063.
Zdobnov EM, Apweiler R: InterProScan–an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.
Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Eddy SR, Bateman A: Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2005, 33: D121-D124.
Lowe TM, Eddy SR: tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.
Karro JE, Yan Y, Zheng D, Zhang Z, Carriero N, Cayting P, Harrrison P, Gerstein M: Pseudogene. org: a comprehensive database and comparison platform for pseudogene annotation. Nucleic Acids Res. 2007, 35: D55-D60. 10.1093/nar/gkl851.
Lagesen K, Hallin P, Rødland EA, Stærfeldt H-H, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35: 3100-3108. 10.1093/nar/gkm160.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W - Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Baxevanis AD: Current protocols in bioinformatics. In Book Current protocols in bioinformatics (Editor ed.^eds.). 2003, City: Wiley
Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-302. 10.1093/bioinformatics/19.2.301.
Padovan ACB, Sanson GFO, Brunstein A, Briones MRS: Fungi evolution revisited: application of the penalized likelihood method to a bayesian fungal phylogeny provides a new perspective on phylogenetic relationships and divergence dates of ascomycota groups. J Mol Evol. 2005, 60: 726-735. 10.1007/s00239-004-0164-y.
Li H, Coghlan A, Ruan J, Coin LJ, Heriche JK, Osmotherly L, Li RQ, Liu T, Zhang Z, Bolund L, et al: TreeFam: a curated database of phylogenetic trees of animal gene families. Nucleic Acids Res. 2006, 34: D572-D580. 10.1093/nar/gkj118.
Enright AJ, Van Dongen S, Ouzounis CA: An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002, 30: 1575-1584. 10.1093/nar/30.7.1575.
Librado P, Vieira FG, Rozas J: BadiRate: estimating family turnover rates by likelihood-based methods. Bioinformatics. 2012, 28: 279-281. 10.1093/bioinformatics/btr623.
Galagan JE, Calvo SE, Cuomo C, Ma LJ, Wortman JR, Batzoglou S, Lee SI, Basturkmen M, Spevak CC, Clutterbuck J, et al: Sequencing of Aspergillus nidulans and comparative analysis with A-fumigatus and A-oryzae. Nature. 2005, 438: 1105-1115. 10.1038/nature04341.
Glass NL, Kaneko I: Fatal attraction: Nonself recognition and heterokaryon incompatibility in filamentous fungi. Eukaryot Cell. 2003, 2: 1-8. 10.1128/EC.2.1.1-8.2003.
Dementhon K, Saupe SJ, Clave C: Characterization of IDI-4, a bZIP transcription factor inducing autophagy and cell death in the fungus Podospora anserina. Mol Microbiol. 2004, 53: 1625-1640. 10.1111/j.1365-2958.2004.04235.x.
Saupe S, Descamps C, Turcq B, Begueret J: Inactivation of the podospora-anserina vegetative incompatibility locus Het-C, whose product resembles a glycolipid transfer protein, drastically impairs ascospore production. P Natl Acad Sci USA. 1994, 91: 5927-5931. 10.1073/pnas.91.13.5927.
Anand S, Prasad MVR, Yadav G, Kumar N, Shehara J, Ansari MZ, Mohanty D: SBSPKS: structure based sequence analysis of polyketide synthases. Nucleic Acids Res. 2010, 38: W487-W496. 10.1093/nar/gkq340.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: research0034-
The whole project has been deposited at DDBJ/EMBL/GenBank under accession number APWS00000000.
This research was supported by the National Natural Science Foundation of China (No. 31070018, 81171524, 31270003, 31000010 and 31070360), the National 973 Program of China Grant (2013CB733904) and the Knowledge Innovation Program of the Chinese Academy of Sciences (KSCX2-EW-J-6).
There are no competing financial interests associated with this work.
YYW performed the biological assays, analyzed the genomic and qRT-PCR data and wrote the draft manuscript. BL was in charge of the genome sequencing and gene annotation. XYZ provided the analysis of genome data. QMZ took part in writing the draft manuscript, designed the biological experiments and analyzed the data. TZ isolated and determined the strain of lichen-forming fungus E. pusillum and its phycobiont D. chodatii, and optimized their culture for the study. HL provided the genes concerning sexual cycle of fungi and constructed two color figures (Figures 1 and 10). YFY and XLZ participated in the analysis of genome data. XYH and MW assisted with the genome sequencing. LW initiated the study and finalized the manuscript. JCW initiated and designed the study, wrote parts of the manuscript, and revised and finalized the manuscript. All authors read and approved the final version of the manuscript.
Yan-Yan Wang, Bin Liu, Xin-Yu Zhang, Qi-Ming Zhou contributed equally to this work.
Electronic supplementary material
Additional file 1:Genomic analysis of Endocarpon pusillum . The file contains additional information on genomic properties and qRT-PCR assays, comprising 10 tables provided in separate excel sheets. Table S1 summarizes the main features of the primary sequence data. Table S2 provides information on genes encoding proteins involved in sexual and asexual reproduction in E. pusillum. Table S3 is a comparison of the number of secreted proteins between E. pusillum and other phytopathogenic fungi. Table S4 lists the domain structures of predicted E. pusillum PKS and NRPS genes. Table S5 lists the genes involving in drought resistant mechanisms in E. pusillum. Table S6 provides information on genes involved in heterokaryon incompatibility in E. pusillum. Table S7 lists the lectins measured by qRT-PCR analysis. Table S8 lists genes that encode nitrogen transporters or proteins involved in nitrogen metabolism in E. pusillum, whose transcriptions were determined by qRT-PCR analysis. Table S9 lists the homologous genes involved in symbiotic fungal sucrose and monosaccharide transporters in E. pusillum. Table S10 lists genes that encode sugar transporters or proteins involved in sugar metabolism in E. pusillum, whose transcriptions were determined by qRT-PCR analysis. Table S11 lists the candidate reference genes for qRT-PCR analysis in E. pusillum. (XLSX 41 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.