Open Access

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

  • Yan-Yan Wang1, 2,
  • Bin Liu3,
  • Xin-Yu Zhang1,
  • Qi-Ming Zhou1,
  • Tao Zhang1, 2,
  • Hui Li1, 2,
  • Yu-Fei Yu1,
  • Xiao-Ling Zhang1,
  • Xi-Yan Hao3,
  • Meng Wang3,
  • Lei Wang3Email author and
  • Jiang-Chun Wei1Email author
Contributed equally
BMC Genomics201415:34

DOI: 10.1186/1471-2164-15-34

Received: 7 October 2013

Accepted: 14 January 2014

Published: 17 January 2014

Abstract

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.

Keywords

Mycobiont Phycobiont Lichenization Symbiosis Symbiosis-related gene Photosynthetic products

Background

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 [1]. 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 [2].

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 [3]. 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 [6]. 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 [9].

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 [1113]. 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 [1417] 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 [22]. 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 [12], 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 [25]. 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 [26]. 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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig1_HTML.jpg
Figure 1

Endocarpon pusillum . A. The lichen E. pusillum. B. The isolated mycobiont and phycobiont [26]. C. Cross section of a perithecium with hymenial algal cells inside. D. Cross section of a thallus under scanning electron microscopy (SEM). E. The algal layer (SEM). F. An algal cell is clasped and surrounded by some hyphae (SEM).

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) [29]. 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).
Table 1

Main features of the Endocarpon pusillum genome

Features

Endocarpon pusillum

Assembly size/Mb

37.5

Scaffold N50/kb

178

Coverage/fold

78

G + C content

46.01%

GC Exonic

51.73%

GC Intronic

47.05%

Repeat rate

1.68%

Protein-coding genes

9,285

Gene density

250.8 per Mbp

Exons per genes

2.53

tRNAs

72

rRNAs

19

SM (Secondary Metabolism) genes

28

TE

15%

https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig2_HTML.jpg
Figure 2

Genome of Endocarpon pusillum . A-E represent the circles from the outside to inside. A: Scaffolds of the genome, filtered by size (≥ 280 Kb). B: ORFs/genes. C: Gene density represented as number of genes per 100 kb (non-overlapping, window size 100 kb). D: Percentage coverage of repetitive sequences (non-overlapping windows, window size 100 kb). E: GC content was estimated by the percent GC in 100 kb non-overlapping windows.

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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig3_HTML.jpg
Figure 3

Functional classification for the Endocarpon pusillum genome. A. FunCat database classification in the second level. B. KOG database classification. C. The top 10 Gene Ontology biological process. D. The top 10 Gene Ontology molecular functions.

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 [30] 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 [9].
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig4_HTML.jpg
Figure 4

Phylogenetic analysis of 15 fungi based on 1,893 orthologous genes. The times of divergence of clades are indicated in millions of years by arrows. The symbols correspond to certain model fungi (open circles), mycorrhizal fungi (filled circles), plant pathogen (filled triangles), animal pathogen (open triangles) and fungi symbiosis with fungi (open squares). C. cinerea and L. bicolor belonging to the Basidiomycota were used as an outgroup.

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 [31]. According to the method proposed by Margolin et al.[32], sequences with a high TA/AT ratio (>0.89) and a low (CA + TG)/(AC + GT) ratio (< 1.03) are thought to indicate RIP [33].

A quantitative alignment-based method, RIPCAL [34], 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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig5_HTML.jpg
Figure 5

Proportion of genes in multi-gene families. A. The graph displays the proportion of genes in multigene families from E. pusillum and other fungi. B. Histogram of paralogous gene numbers with different levels of amino acid similarity in E. pusillum and other fungi.

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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig6_HTML.jpg
Figure 6

Configuration of the MAT locus in Endocarpon pusillum .

Secreted proteins

Secreted proteins are important effectors that modulate the interaction between pathogenic microbes and hosts [3537], and many small secreted proteins (SSPs) are virulence factors [38]. 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).

Secondary metabolism

To date, several polyketide synthase (PKS) genes have been identified in different lichen-forming fungi [3942], 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 [43]. 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.
Table 2

Expanded and lineage-specific protein families in Endocarpon pusillum

Protein family ID

Number of proteins in E. pusillum

p-value

HMMPfam description

group10696

20

0

Hypothetical protein CIMG_07368

group10024

15

0

Ankyrin and HET domain-containing protein

group14580

14

0

Ankyrin repeat

group10771

13

0

P-loop containing nucleoside triphosphate hydrolases

group11690

13

0

MFS general substrate transporter

group10071

11

0

Putative transposase

group10011

10

0

Hypothetical protein NECHADRAFT_87861

group14903

8

0

Heterokaryon incompatibility protein

group17019

7

1E-07

DNA/RNA polymerases

group10674

10

3E-07

Magnesium transport protein CorA, transmembrane region

group10674

10

3E-07

Hypothetical protein AN9301.2

group10038

8

5E-07

Hypothetical protein SS1G_06795

group10340

8

5E-07

Conserved hypothetical protein

group15154

6

5E-07

Hypothetical protein CHGG_10108

group16626

6

5E-07

Hypothetical protein

group18347

6

5E-07

Hypothetical protein

group10002

9

2.4E-06

Ribonuclease H-like

group10015

9

2.4E-06

Ankyrin repeat

group10097

7

4.6E-06

Ankyrin repeat

group10201

7

4.6E-06

ARM repeat

group10224

7

4.6E-06

Heterokaryon incompatibility protein

group10843

7

4.6E-06

Predicted protein

group11854

7

4.6E-06

Cytochrome P450

group10000

11

0.000001

P-loop containing nucleoside triphosphate hydrolases

group10705

5

0.000005

Hypothetical protein

group11683

5

0.000005

Dimeric alpha + beta barrel

group11909

5

0.000005

Unnamed protein product

group13886

5

0.000005

Ankyrin repeat

group15448

5

0.000005

FabD/lysophospholipase-like

group16624

5

0.000005

Protein kinase-like (PK-like)

group16962

5

0.000005

Hypothetical protein NFIA_004500

group17139

5

0.000005

Cytochrome P450

group18335

5

0.000005

Predicted protein

group18344

5

0.000005

Conserved hypothetical protein

group19016

5

0.000005

Hypothetical protein

group19022

5

0.000005

Hypothetical protein

group10704

8

1.78E-05

Conserved hypothetical protein

group10239

6

3.81E-05

P-loop containing nucleoside triphosphate hydrolases

group10495

6

3.81E-05

Unnamed protein product

group10677

6

3.81E-05

Unnamed protein product

group11482

6

3.81E-05

Predicted protein

group11606

6

3.81E-05

Hypothetical protein AOR_1_1448144

group15933

6

3.81E-05

Purine and uridine phosphorylases

group16361

6

3.81E-05

Hypothetical protein SS1G_13793

group16513

6

3.81E-05

Conserved hypothetical protein

group17376

6

3.81E-05

Hypothetical protein

group10014

9

4.84E-05

POZ domain

group10175

9

4.84E-05

Hypothetical protein CIMG_05613

group10238

4

4.94E-05

WD40 repeat-like

group11337

4

4.94E-05

TPR-like

group12140

4

4.94E-05

WD40 repeat-like

group16268

4

4.94E-05

Hypothetical protein MGYG_01629

group17166

4

4.94E-05

Predicted protein

group17167

4

4.94E-05

Alpha/beta-Hydrolases

group17374

4

4.94E-05

Hypothetical protein FRAAL5044

group17553

4

4.94E-05

Protein kinase-like (PK-like)

group17955

4

4.94E-05

Conserved hypothetical protein

group18350

4

4.94E-05

Hypothetical protein CHGG_08502

group19010

4

4.94E-05

NTF2-like

Lineage-specific protein families

Protein family ID

Number of proteins in E. pusillum

p-value

Pfam description

group19669

4

4.94E-05

Ribonuclease H-like

group19844

4

4.94E-05

Transposase

group19847

4

4.94E-05

Hypothetical protein

group19848

4

4.94E-05

Hypothetical protein

group19850

4

4.94E-05

Protein kinase-like (PK-like)

group19853

4

4.94E-05

Ankyrin repeat

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[23].

Similarly, G protein-coupled receptors (GPCRs), which sense molecules outside the cell and activate signal transduction pathways inside the cell [44], 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).
Table 3

Evolutionary changes of GPCR protein families in Endocarpon pusillum

GPCRs family

Number of protein in E. pusillum

Evolutionary changes in E. pusillum

Annotation

Class A

6

Gain

Rhodopsin-like family

Class B

2

Gain

Secretin receptor

Class D

5

Gain

Fungal mating pheromone receptor

Class E

1

Gain

cAMP receptor

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 [4547]. 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 [48], 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 [49], 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 [49]. There are multiple group I introns present in the nuclear rDNA of many lichen-forming fungi [5053], and they are considered to decrease the growth rate of lichen-forming fungi by interfering with rRNA maturation [54]. 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 [57]. Magnesium is the most important divalent cation in cells, and is particularly important for photosynthesis [58]. 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 [61]. 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.
Table 4

Evolutionary changes of transporter protein families in Endocarpon pusillum (partial results)

Transporter family

Number of protein in E. pusillum

Evolutionary changes in E. pusillum

Annotation

1.A.11.3.3

3

Gain

Ammonium transporter-Hebeloma cylindrosporum

2.A.3.10.20

1

Gain

Lysine/arginine permease-Candida albicans

2.A.3.4.2

16

Gain

Gaba permease-Emericella nidulans-Aspergillus nidulans

2.A.3.8.4

6

Gain

High affinity methionine permease-Saccharomyces cerevisiae

2.A.39.3.3

1

Gain

Uridine permease-Saccharomyces cerevisiae

8.A.9.1.1

3

Gain

Amino acid transport related protein (RBAT)-Oryctolagus cuniculus

2.A.1.1.10

0

Loss

Maltose permease MAL6T-Saccharomyces cerevisiae

2.A.1.1.11

1

Loss

General alpha-glucoside permease-Saccharomyces cerevisiae

2.A.1.1.33

0

Loss

Hexose transporter-Kluyveromyces lactis

2.A.1.1.38

2

Loss

Sugar transporter STL1-Saccharomyces cerevisiae

2.A.1.1.39

3

Loss

High-affinity glucose transporter-Kluyveromyces lactis

2.A.1.1.57

1

Loss

Monosaccharide transporter-Aspergillus niger

2.A.1.1.58

0

Loss

Monosaccharide transporter-Aspergillus niger

2.A.1.1.68

0

Loss

Glucose transporter/Sensor OS-Pichia stipitis

2.A.1.1.9

0

Loss

Lactose permease-Kluyveromyces lactis

2.A.1.2.23

1

Loss

Fructose facilitator-Zygosaccharomyces bailii

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 [6264]. 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 [65]. 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 [66], 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 [67]. 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 [6875]. 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 [76] 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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig7_HTML.jpg
Figure 7

Relative transcription levels of certain symbiosis-related genes in Endocarpon pusillum when cultured with Diplosphaera chodatii . The names of the genes are listed on the top of the picture. Circles represent genes encoding lectins, triangles represent N transporter genes, squares represent gene encoding enzymes of nitrogen metabolism, a heptagon represents the tetracycline resistance protein gene, hexagons represent sugar transporter genes, and pentagons represent genes encoding enzymes of sugar metabolism. The culture conditions are displayed at the right of the picture, in which the ratios of mycobiont (M) and phycobiont (A) and incubation time are given. The log2 relative expression with respect to the control group is illustrated in the heat map.

Table 5

Sample treatments to identify differentially expressed genes

Experiment group

Sample

Weight ratio (mycobiont:Phycobiont)

Incubation time (h)

Control A

Mycobiont

-

24

Control B

Mycobiont

-

72

Control C

Phycobiont

-

72

I

Mycobiont + phycobiont

10:1

24

II

Mycobiont + phycobiont

10:3

24

III

Mycobiont + phycobiont

10:1

72

IV

Mycobiont + phycobiont

10:3

72

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) [79]. Lectins are sugar-binding proteins that play an important roles in cell recognition [80]. The first lectin discovered in lichens was that found in Peltigera canina and Peltigera polydactyla[17]. 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 [84], and the evolution of galectin lec-2 was driven by interaction with different strains of the phycobiont Nostoc[85]. 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 [6264]. 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 [76], 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 [20]. Additionally, previous studies showed that some monosaccharides, such as glucose, and disaccharides, such as trehalose and sucrose, were transferred between non-lichenized symbionts [6875]. 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 [88], 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 [89]. 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 [92]. 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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig8_HTML.jpg
Figure 8

Relative transcription levels of sugar transfer-related genes in Endocarpon pusillum after supplementation with different carbohydrates. The names of the genes are listed on the top of the picture. Heptagons represent the gene encoding the tetracycline transporter, hexagons represent sugar transporter genes, and pentagons represent genes encoding enzymes of sugar metabolism. The culture time is displayed on the left and sugars added into BBM are listed on the right. The log2 relative expression with respect to the control group is illustrated in the heat map.

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 [19]. 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 [21]. 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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig9_HTML.jpg
Figure 9

Ion chromatogram of co-cultured sample from experimental group IV. The sample from experiment group IV (mycobiont: phycobiont = 10:3; culture for 72 hrs) was diluted 25 times before four carbohydrates (glucose, mannitol, sorbitol and sucrose) were measured using ion chromatography.

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 [9395] and could provide clues for the evolution of lichenization.

Conclusions

Approximately 40% of all Ascomycota are lichen-forming; thus, lichenization is regarded as one of the most important fungal lifestyles [96]. 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.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-34/MediaObjects/12864_2013_Article_5660_Fig10_HTML.jpg
Figure 10

Summary of the interactions between symbionts in lichen Endocarpon pusillum . The lectins lacking signal peptides, which are shown as ellipse marked with L, act as recognition factors in the direct contact between E. pusillum and D. chodatii. After it is captured by the compatible fungal partner, the photosynthetic products sorbitol (S), glucose (G) and sucrose (G-F) are released from the phycobiont D. chodatii. These sugars are then absorbed through sugar transporters of E. pusillum and converted into glucose or fructose to support fungal metabolism. In addition, there are a small amount of other uncertainty carbon sources (X) released by the phycobiont, and the mycobiont can also utilize mannitol (M) from environment. At the same time, organic and inorganic nitrogen, together with various ions in the environment, are transferred into E. pusillum. In fungal cells, these substances may be converted into the forms preferred by the phycobiont and are then delivered to D. chodatii.

Methods

Fungal strains

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 [26]. 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 [97], which adopts the de Bruijn graph data structure to construct contigs.

Gene prediction and annotation

We used Augustus [98], GeneID [99], and GeneMark-ES [100] programs to predict the gene models for E. pusillum. A final set of gene models was selected by EvidenceModeler [101]. 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 [102], KOG [103], FunCat [104] and Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis [105]. Protein domains were predicted using InterProScan [106] 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 [107], and tRNAs were predicted using tRNAscan-SE [108]. Pseudogenes and rRNAs were designated using PseudoGene [109] and RNAmmer [110] respectively.

Orthologous gene and phylogeny analysis of E. pusillum

The sequences of corresponding orthologous genes from 15 fungi were aligned using ClustalW [111]. A maximum parsimony (MP) phylogenomic tree was created using the concatenated amino acid sequences in phylogenetic analysis using parsimony (PAUP) [112], 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 [113]. The time of divergence between Ascomycota and Basidiomycota is set as 500 Myr, and the time between Pezizomycotina and Saccharomycotina was 350 Myr [114].

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 [115]. Multigene families were identified using the MCL method [116]. 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 [117], 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.

Mating type

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 [118], N. crassa [119] 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 [122].

Secreted proteins

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/].

Quantitative RT-PCR

RNA was extracted from cultured lichen-forming fungus E. pusillum. cDNA synthesis and relative quantitative RT-PCR were carried out as described previously [28]. 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 [123]. 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).

Notes

Abbreviations

BBM: 

Bold’s basal medium

bp: 

Base pair

FunCat: 

Functionally annotation

GO: 

Gene ontology

GPCR: 

G-protein coupled receptor

HET: 

Heterokaryon incompatibility

HI: 

Heterokaryon incompatibility

HMAS: 

Mycological herbarium of microbiology institute,the chinese academy of sciences

KEGG: 

Kyoto encyclopedia of genes and genomes database

KOG: 

Eukaryotic clusters of orthologous groups

LSUrDNA: 

Large-subunit ribosomal DNA

NCBI: 

National center for biotechnology information

MAT: 

Mating-type

MCL: 

Markov cluster

MFS: 

Major facilitator superfamily

Myr: 

million years

NRPS: 

Non-ribosomal peptide synthetase

PCR: 

Polymerase chain reaction

PKS: 

Polyketide synthase

qRT-PCR: 

Quantitative real-time PCR

RIP: 

Repeat-induced point mutation

SEM: 

Scanning electron microscopy

SSUrDNA: 

Small-subunit ribosomal DNA

T6P: 

Trehalose-6-phosphate

TPS: 

Trehalose 6-phospate synthase

SSP: 

Small secret protein.

Declarations

Accession number

The whole project has been deposited at DDBJ/EMBL/GenBank under accession number APWS00000000.

Acknowledgments

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).

Authors’ Affiliations

(1)
State Key Laboratory of Mycology, Institute of Microbiology, Chinese Academy of Sciences
(2)
University of Chinese Academy of Sciences
(3)
TEDA School of Biological Sciences and Biotechnology, Nankai University

References

  1. Ahmadjian V: The lichen symbiosis. 1993, New York: John WileyGoogle Scholar
  2. Kirk PMCP, Minter DW, Stalpers JA: Dictionary of the fungi. 2008, Wallingford: CABI, 10Google Scholar
  3. Nash TH: Lichen biology. 2nd edn. Cambridge. 2008, New York: Cambridge University PressGoogle Scholar
  4. Ingolfsdottir K: Usnic acid. Phytochemistry. 2002, 61: 729-736. 10.1016/S0031-9422(02)00383-7.Google Scholar
  5. 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.Google Scholar
  6. Honegger R: The lichen symbiosis—what is so spectacular about it?. The Lichenologist. 1998, 30: 193-212.Google Scholar
  7. 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.Google Scholar
  8. 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.Google Scholar
  9. Lutzoni F, Pagel M, Reeb V: Major fungal lineages are derived from lichen symbiotic ancestors. Nature. 2001, 411: 937-940. 10.1038/35082053.PubMedGoogle Scholar
  10. Rikkinen J: Molecular studies on cyanobacterial diversity in lichen symbioses. Lichens: from genome to ecosystems in a changing world MycoKeys. 2013, 6: 3-32.Google Scholar
  11. 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.PubMedGoogle Scholar
  12. Ahmadjian VHH: The culture and synthesis of Endocarpon pusillum and Staurothele clopima. Lichenologist. 1970, 4: 13-Google Scholar
  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.Google Scholar
  14. Guillot J, Konska G: Lectins in higher fungi. Biochem Syst Ecol. 1997, 25: 203-230. 10.1016/S0305-1978(96)00110-X.Google Scholar
  15. Kardish NSL, Flemminger N, Galun M: Lectin from the lichen Nephroma laevigatum. Localization and function. Symbiosis. 1991, 11: 47-62.Google Scholar
  16. 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.PubMedGoogle Scholar
  17. 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.Google Scholar
  18. Richardson DHS, Hill DJ, Smith DC: Lichen physiology. New Phytol. 1968, 67: 469-486. 10.1111/j.1469-8137.1968.tb05476.x.Google Scholar
  19. 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.Google Scholar
  20. Hill DJ, Ahmadjia V: Relationship between carbohydrate movement and symbiosis in lichens with Green-Algae. Planta. 1972, 103: 267-10.1007/BF00386850.PubMedGoogle Scholar
  21. Hawksworth DL, Hill DJ: The lichen-forming fungi. Blackie. 1984, 5: 62-69.Google Scholar
  22. 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.Google Scholar
  23. 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.PubMedGoogle Scholar
  24. Suzanne Joneson FL: Compatibility and thigmotropism in the lichen symbiosis: a reappraisal. Symbiosis. 2009, 47: 109-115. 10.1007/BF03182294.Google Scholar
  25. Awasthi DD: A hand book of lichens. 2000, Bishen Singh Mahendra Pal Singh: DehradunGoogle Scholar
  26. 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.PubMedGoogle Scholar
  27. 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.Google Scholar
  28. 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.PubMed CentralPubMedGoogle Scholar
  29. 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.Google Scholar
  30. 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.PubMed CentralPubMedGoogle Scholar
  31. 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.PubMedGoogle Scholar
  32. 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.PubMed CentralPubMedGoogle Scholar
  33. 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.PubMedGoogle Scholar
  34. 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.PubMed CentralPubMedGoogle Scholar
  35. 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.Google Scholar
  36. 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.PubMedGoogle Scholar
  37. 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.PubMedGoogle Scholar
  38. 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.PubMedGoogle Scholar
  39. 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.PubMedGoogle Scholar
  40. 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.PubMedGoogle Scholar
  41. 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.PubMedGoogle Scholar
  42. 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.Google Scholar
  43. 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.Google Scholar
  44. Gilman AG: G-Proteins - Transducers of receptor-generated signals. Annu Rev Biochem. 1987, 56: 615-649. 10.1146/annurev.bi.56.070187.003151.PubMedGoogle Scholar
  45. 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.PubMedGoogle Scholar
  46. 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.PubMedGoogle Scholar
  47. 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.PubMed CentralPubMedGoogle Scholar
  48. 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.Google Scholar
  49. 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.PubMedGoogle Scholar
  50. 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.PubMedGoogle Scholar
  51. 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.PubMedGoogle Scholar
  52. Gargas A, DePriest PT, Taylor JW: Positions of multiple insertions in SSU rDNA of lichen-forming fungi. Mol Biol Evol. 1995, 12: 208-218.PubMedGoogle Scholar
  53. 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.PubMedGoogle Scholar
  54. 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.Google Scholar
  55. 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.PubMedGoogle Scholar
  56. Kehres DG, Maguire ME: Structure, properties and regulation of magnesium transport proteins. Biometals. 2002, 15: 261-270. 10.1023/A:1016078832697.PubMedGoogle Scholar
  57. 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.PubMedGoogle Scholar
  58. Gardner RC: Genes for magnesium transport. Curr Opin Plant Biol. 2003, 6: 263-267. 10.1016/S1369-5266(03)00032-3.PubMedGoogle Scholar
  59. 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.PubMedGoogle Scholar
  60. 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.PubMedGoogle Scholar
  61. 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.PubMedGoogle Scholar
  62. Crittenden PD: The effect of oxygen deprivation on inorganic nitrogen uptake in an Antarctic macrolichen. Lichenologist. 1996, 28: 347-354.Google Scholar
  63. 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.Google Scholar
  64. Dahlman L, Persson J, Palmqvist K, Nasholm T: Organic and inorganic nitrogen uptake in lichens. Planta. 2004, 219: 459-467.PubMedGoogle Scholar
  65. 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.Google Scholar
  66. 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.PubMed CentralPubMedGoogle Scholar
  67. 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.PubMed CentralPubMedGoogle Scholar
  68. 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.PubMed CentralPubMedGoogle Scholar
  69. 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.PubMed CentralPubMedGoogle Scholar
  70. 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.PubMed CentralPubMedGoogle Scholar
  71. 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.Google Scholar
  72. 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.PubMedGoogle Scholar
  73. 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.Google Scholar
  74. 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.PubMedGoogle Scholar
  75. 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.PubMedGoogle Scholar
  76. 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.PubMedGoogle Scholar
  77. Schnappinger D, Hillen W: Tetracyclines: antibiotic action, uptake, and resistance mechanisms. Arch Microbiol. 1996, 165: 359-369. 10.1007/s002030050339.PubMedGoogle Scholar
  78. Walsh C: Molecular mechanisms that confer antibacterial drug resistance. Nature. 2000, 406: 775-781. 10.1038/35021219.PubMedGoogle Scholar
  79. 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.PubMed CentralPubMedGoogle Scholar
  80. Sharon N, Lis H: Lectins as cell recognition molecules. Science. 1989, 246: 227-234. 10.1126/science.2552581.PubMedGoogle Scholar
  81. 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.PubMed CentralPubMedGoogle Scholar
  82. Galun MKN: Lectins as determinants of symbiotic specificity in lichens. Cryptog Bot. 1995, 5: 144-148.Google Scholar
  83. 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.Google Scholar
  84. 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.Google Scholar
  85. 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.PubMed CentralPubMedGoogle Scholar
  86. Liener IE, Sharon N, Goldstein IJ: The Lectins: properties, functions, and applications in biology and medicine. 1986, Orlando: AcademicGoogle Scholar
  87. 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.PubMedGoogle Scholar
  88. Iturriaga G, Suarez R, Nova-Franco B: Trehalose metabolism: from osmoprotection to signaling. Int J Mol Sci. 2009, 10: 3793-3810. 10.3390/ijms10093793.PubMed CentralPubMedGoogle Scholar
  89. 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.PubMedGoogle Scholar
  90. 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.PubMedGoogle Scholar
  91. 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.Google Scholar
  92. Pardy RL, Spargo B, Crowe JH: Release of trehalose by symbiotic algae. Symbiosis. 1989, 7: 149-158.Google Scholar
  93. Nelsen MP, Gargas A: Dissociation and horizontal transmission of codispersing lichen symbionts in the genus Lepraria (Lecanorales: Stereocaulaceae). New Phytol. 2008, 177: 264-275.PubMedGoogle Scholar
  94. Piercey-Normore MD, DePriest PT: Algal switching among lichen symbioses. American Journal of Botany. 2001, 88: 1490-1498. 10.2307/3558457.PubMedGoogle Scholar
  95. 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.PubMedGoogle Scholar
  96. 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.Google Scholar
  97. 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.Google Scholar
  98. 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.PubMed CentralPubMedGoogle Scholar
  99. Parra G, Blanco E, Guigó R: Geneid in drosophila. Genome research. 2000, 10: 511-515. 10.1101/gr.10.4.511.PubMed CentralPubMedGoogle Scholar
  100. 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.PubMed CentralPubMedGoogle Scholar
  101. 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.PubMed CentralPubMedGoogle Scholar
  102. 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.PubMed CentralPubMedGoogle Scholar
  103. 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.PubMedGoogle Scholar
  104. 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.PubMed CentralPubMedGoogle Scholar
  105. 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.PubMed CentralPubMedGoogle Scholar
  106. 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.PubMedGoogle Scholar
  107. 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.PubMed CentralPubMedGoogle Scholar
  108. 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.PubMed CentralPubMedGoogle Scholar
  109. 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.PubMed CentralPubMedGoogle Scholar
  110. 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.PubMed CentralPubMedGoogle Scholar
  111. 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.PubMed CentralPubMedGoogle Scholar
  112. Baxevanis AD: Current protocols in bioinformatics. In Book Current protocols in bioinformatics (Editor ed.^eds.). 2003, City: WileyGoogle Scholar
  113. 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.PubMedGoogle Scholar
  114. 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.PubMedGoogle Scholar
  115. 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.PubMed CentralPubMedGoogle Scholar
  116. 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.PubMed CentralPubMedGoogle Scholar
  117. Librado P, Vieira FG, Rozas J: BadiRate: estimating family turnover rates by likelihood-based methods. Bioinformatics. 2012, 28: 279-281. 10.1093/bioinformatics/btr623.PubMedGoogle Scholar
  118. 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.PubMedGoogle Scholar
  119. 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.PubMed CentralPubMedGoogle Scholar
  120. 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.PubMedGoogle Scholar
  121. 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.Google Scholar
  122. 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.PubMed CentralPubMedGoogle Scholar
  123. 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-PubMed CentralPubMedGoogle Scholar

Copyright

© Wang et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.