Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Genome sequence of the ectophytic fungus Ramichloridium luteum reveals unique evolutionary adaptations to plant surface niche

BMC Genomics volume 18, Article number: 729 (2017) Cite this article

Abstract

Background

Ectophytic fungi occupy the waxy plant surface, an extreme environment characterized by prolonged desiccation, nutrient limitation, and exposure to solar radiation. The nature of mechanisms that facilitate adaptation to this environment remains unclear. In this study, we sequenced the complete genome of an ectophytic fungus, Ramichloridium luteum, which colonizes the surface of apple fruit, and carried out comparative genomic and transcriptome analysis.

Results

The R. luteum genome was 28.18 Mb and encoded 9466 genes containing 1.85% repetitive elements. Compared with cell-penetrating pathogens, genes encoding plant cell wall degrading enzymes (PCWDEs), PTH11-like G protein-coupled receptors (GPCRs) and effectors were drastically reduced. In contrast, genes encoding cutinases and secretory lipases were strikingly expanded, and four of nine secretory lipases were probably acquired by horizontal gene transfer from Basidiomycota. Transcriptomic analysis revealed elevated expression of genes involved in cuticle degradation (cutinase, secretory lipase) and stress responses (melanin biosynthesis, aquaporins, lysozymes and HOG pathway).

Conclusions

Taken together, our results highlight genomic features associated with evolution of surface niche adaptation by the ectophytic fungus R. luteum, namely the contraction of PCWDEs, PTH11-like GPCRs and effectors, and the expansion of cuticle degradation and stress tolerance.

Background

Plant-colonizing microorganisms occupy two broad types of ecological niche: interior and surface [1]. Invasive phytopathogens, symbionts and endophytes establish populations in the intracellular or apoplastic spaces, which constitute the interior niche [2, 3]. Epiphytic and ectophytic organisms, on the other hand, colonize primarily the surface niche. Ectophytes differ from epiphytes in that they interact with and absorb nutrients from extracellular components of plant surfaces [4,5,6].

Spanu [2] concluded that ancestral saprophytic fungi evolved into cell-penetrating plant pathogens under suitable conditions. Adaptations required for invasive plant pathogenesis include expansion and diversification of plant cell wall depolymerases, effectors and metabolic phytotoxins [2]. In contrast, sooty blotch and flyspeck (SBFS) fungi, which descended from invasive plant pathogens, colonize epicuticular layers of plants as ectophytes [6]. They can extract required nutrients from living hosts without killing cells and tissues. During the evolutionary transition from invasive plant pathogenesis to an ectophytic niche [6, 7], most pathogenicity-related genes except some associated with cuticle penetration were lost in SBFS fungus Peltaster fructicola [6]. Transition from invasive pathogenesis to the ectophytic niche represents an important but little-studied pathway in plant-fungus co-evolution.

SBFS can cause substantial economic damage worldwide on fruit of apple, pear, persimmon and many other cultivated crops by forming superficial dark-colored blemishes which downgrade fruit for fresh-market sale and accelerate fruit desiccation during cold storage [8]. SBFS fungi include more than 80 species, most of which belong to the order Capnodiales within the class Dothideomycetes [8]. SBFS species are polyphyletic, exhibiting convergent evolution of their specialized ectophytic colonization lifestyle [9,10,11,12,13,14].

The plant surface micro-environment is characterized by direct sunlight exposure, a hydrophobic waxy surface, and scarcity of available water and nutrients. To survive, ectophytic fungi must evolve strategies to deal with these abiotic stresses. In addition, ectophytes must also compete with other surface-dwelling microorganisms for nutrient absorption and evade or counteract the toxic effects of plant-derived antimicrobial metabolites [15]. Although living as an ectophyte could reduce elicitation of host countermeasures, ectophytes must cope with the harsh environment at the plant surface. During surface colonization, SBFS fungi produce darkly pigmented (melanized) hyphae, which can enable them to resist extreme temperatures, desiccation, and ionizing radiation [16]. Xu et al. [6] sequenced the genome of the SBFS fungus P. fructicola and found that genes controlling 1,8-dihydroxynaphthalene (DHN)-melanin biosynthesis are highly expressed on the apple surface. However, the full range of adaptations by SBFS fungi to plant surfaces remains unclear. Moreover, it is unclear whether such genomic features are general characteristics of the SBFS fungal complex.

In this study, we carried out genomic and transcriptomic sequencing of a SBFS fungus in the genus Ramichloridium, R. luteum G.Y. Sun, H.Y. Li & Crous [13]. The genus Ramichloridium Stahel ex de Hoog accommodates a wide range of species with diverse lifestyles [13, 17]. Batzer et al. [9] provided evidence that Ramichloridium is in Chaetothyriales whereas Peltaster is in Dothideales, indicating that Ramichloridium spp. differ from P. fructicola at the level of taxonomic order. Morphologically, species in genus Ramichloridium also differ notably from those in Peltaster [17, 18]. The main goals of this work were to understand genomic features underlying the evolutionary adaptation of R. luteum to the plant surface niche and to compare adaptions for ectophytic colonization to those of P. fructicola.

Results

General genome features

The final genome assembly of R. luteum contained 155 scaffolds with a total length of 28.18 Mb and a coverage of 161-fold. R. luteum genome size was considerably smaller than the average genome size of Ascomycota (36.91 Mb) [19], but over 55% larger than the SBFS fungus P. fructicola (genome size 18.14 Mb) [6]. The genome scaffold N50 was 437 kb, the largest scaffold was 1.7 Mb, and there were 108 scaffolds longer than 3 kb. Seven assembly scaffolds ended with the typical fungal telomere repeat (TTAGGG)n. The completeness of the genome assembly was assessed by the CEGMA, which gave 93.95% full coverage and 95.56% partial coverage. The GC content of the genome was 54.72% and the genome contained 22 rRNAs and 78 tRNAs (Additional file 1: Table S1).

The repeat content of R. luteum was 1.85% as determined by RepeatModeler and RepeatMasker with very low repetitive content of transposable elements (TEs, 0.8%). The repeat sequence profile indicated a significant repeat-induced point mutation (RIP) in the R. luteum genome, where the TpA/ApT ratio was 1.20 and the (CpA + TpG)/(ApC + GpT) ratio was 0.83. The R. luteum genome contained the DNA methyltransferase RID gene required for RIP.

The R. luteum genome encoded 9466 predicted proteins, among which 5816 (61.44%), 5857 (61.87%), 3301 (34.87%), and 8775 (92.70%) could be annotated based on Gene ontology (GO) (Additional file 2: Figure S1), Clusters of Orthologous Groups (COG) (Additional file 3: Figure S2), Kyoto Encyclopaedia of Genes and Genomes (KEGG), and NCBI nr databases, respectively.

Phylogenomic analysis and protein family evolution

Genomes of 13 species of fungi with various lifestyles including symbionts, free-living saprophytes and invasive phytopathogens were selected for comparison. A total of 1656 single-copy orthologs conserved across all compared species (Additional file 1: Table S2) were selected to construct the phylogenomic tree. All tree branches received 100% bootstrap support, indicating a highly reliable topological structure (Fig. 1). The results showed that the Ramichloridium lineage was evolutionarily closer to the penetrating plant pathogen Zymoseptoria tritici than to the SBFS species P. fructicola. The r8s analysis showed that the divergence time between R. luteum and Z. tritici was around 167 million years ago (mya), whereas that between R. luteum and P. fructicola was around 225 mya. This phylogenetic pattern indicates that the ectophytic lifestyles of R. luteum and P. fructicola have evolved independently. CAFE analysis showed that protein families including cutinase, secretory lipase, fatty acid hydroxylase, and hydrophobic surface binding protein A (HsbA) were expanded and protein families including PWCDEs, proteins containing cysteine-rich domain (CFEM domain), LRR_8 domain and zf-MYND domain involved in protein-protein interactions were contracted in R. luteum (Fig. 2).

Fig. 1
figure 1

A maximum likelihood phylogenetic tree constructed from concatenated alignment of 1656 single-copy orthologs conserved across all 13 species using RAxML with the best-fit model of LG + I + G + F. Bootstrap values based on 100 replicates are shown for indicated nodes. Red color represents the divergence time in million years. Boxes represent the life style of selected organisms. E, ectophytes; H, hemibiotrophs; N, necrotrophs; S, saprotrophs; Y, symbionts

Full size image
Fig. 2
figure 2

Significant expansions and contractions of protein family in R. luteum. Evolution analysis was processed by CAFE. Numbers of protein family members in each genome are shown. Fungal names: BC = Botrytis cinerea; CH = Colletotrichum higginsianum; LB = Laccaria bicolor; MO = Magnaporthe oryzae; ZT = Zymoseptoria tritici; NP = Neofusicoccum parvum; NC = Neurospora crassa; PF = Peltaster fructicola; PT = Pyrenophora tritici-repentis; RL = Ramichloridium luteum; SS = Sclerotinia sclerotiorum; TR = Trichoderma reesei; CG = Cenococcum geophilum. Overrepresented (+3 to 0) and underrepresented (0 to −3) numbers are depicted as Z-scores for each line in heatmap. The boxes on the nethermost represent the life style of the selected organisms. E, ectophytes; H, hemibiotrophs; N, necrotrophs; S, saprotrophs; Y, symbionts

Full size image

Contraction of PTH11-like GPCRs and candidate effectors

The PTH11 G protein-coupled receptors (GPCRs) are a family of cell surface integral membrane proteins required for pathogenicity [20, 21]. Proteins in this family contain seven transmembrane regions and an amino-terminal extracellular CFEM domain. The number of predicted CFEM proteins varied from 3 in R. luteum to 20 in Magnaporthe oryzae or Colletotrichum higginsianum (Fig. 2). R. luteum had fewer CFEM proteins than P. fructicola (5 members) and Z. tritici (8 members). CAFE analysis confirmed the CFEM contraction in R. luteum (P = 0.003). The three R. luteum proteins contained 0–1 transmembrane helices, ruling out the possibility that they are PTH11-like GPCRs.

Effector proteins of plant pathogenic fungi are low-molecular-weight, typically cysteine-rich, secreted compounds that are key virulence factors modulating host defense reactions and cellular activities to facilitate infection [22]. Effector is also necessary for the mutualistic symbionts facilitating hyphal entry into the root and required for microbial-microbial growth competitions [23,24,25]. Candidate effectors (putatively secreted, cysteine-rich, <200 aa) ranged from 34 in the ectophytic fungus P. fructicola to 234 in the cell-penetrating species C. higginsianum (Fig. 3). R. luteum had 59 such proteins, more than the penetrating symbiont Dothidiomycete Cenococcum geophilum (45), similar to the saprotrophic fungi Neurospora crassa (57) and Trichoderma reesei (62), but significantly fewer than penetrating plant pathogens (ranging from 89 to 234) and the penetrating symbiont Basidiomycota Laccaria bicolor (119). The hemibiotrophic plant pathogen Z. tritici had 95 candidate effectors, more than both R. luteum and P. fructicola but fewer than all host-penetrating species except for the necrotrophic fungus Botrytis cinerea (89). Overall, R. luteum has a significantly reduced effector gene content compared with penetrating phytopathogens.

Fig. 3
figure 3

Numbers of genes encoding putative candidate effectors, CAZY modules (AAs, auxiliary activities; GHs, glycoside hydrolases; CBMs, carbohydrate-binding modules), PCWDEs, putative cutinases and secretory lipases identified in the genomes of R. luteum and 12 additional fungal species included in this study. Boxes on the left represent the life style of the selected organisms. E, ectophytes; H, hemibiotrophs; N, necrotrophs; S, saprotrophs; Y, symbionts

Full size image

Fortynine of the 59 candidate effectors in R. luteum contained unknown domains, whereas only 10 contained a PFAM domain, indicating an enrichment of genes with putatively unknown function (Additional file 1: Table S3). Eight genes among 59 candidate effectors were up-regulated during ectophytic growth. One of 8 up-regulated genes, Ramle7664 is taxon-specific since it is found only in R. luteum. The up-regulation of expression of putative effectors is conceivably required for adaptation to a surface niche, such as warding off microbial antagonists.

Loss of plant cell wall-degrading enzymes

The R. luteum genome contained 474 carbohydrate active enzyme (CAZY) modules (Additional file 1: Table S4), which was more than P. fructicola (363), L. bicolor (427) and C. geophilum (436), but fewer than other plant pathogenic and saprophytic fungal species (493 to 839). Furthermore, CAZY annotation demonstrated that R. luteum contains a much smaller set of glycoside hydrolases (GHs) (163 vs. average 218, **, P < 0.005, t-test), auxiliary activities (AAs) (47 vs. average 75, **, P = 0.008, t-test), and carbohydrate binding modules (CBMs) (36 vs. average 58, *, P = 0.011, t-test) compared to the average number for all 13 species (Fig. 3). The strong reduction of GHs, AAs and CBMs in R. luteum suggested a reduced capacity to degrade plant and fungal cell walls. Plant cell walls are composed mostly of cellulose, hemicellulose and pectin, and pathogens attack them by secreting PCWDEs [26]. We further found that R. luteum and P. fructicola had similar numbers of PCWDEs, 56 and 53 respectively, from CAZY genes, which were much fewer than Z. tritici (102 members) and other all non-SBFS fungal species with various lifestyles (**, P < 0.01, t-test) (Fig. 3). Numbers of PCWDE genes in SBFS fungi were fewer than in the penetrating symbiont fungi L. bicolor and C. geophilum (70 and 82, respectively). Hierarchical clustering analysis showed that the PCWDE profile of R. luteum is most closely related to that of P. fructicola and C. geophilum, followed by Z. tritici [27] (Fig. 4). CAFE analysis showed that GH28 (0 vs. average 7.2) and GH43 (0 vs. 8.4) were most strikingly reduced in R. luteum (P < 0.05), and there was one GH28 and one GH43 in P. fructicola (Fig. 2), indicating that putative enzymes involved in the degradation of pectin were contracted in SBFS fungi. Only 5 genes involved in pectin degradation in R. luteum and 4 genes in P. fructicola were found (Fig. 4). These numbers are less than for penetrating symbiont fungi (17 to 21) and penetrating phytopathogens (16 to 82). Fungal enzymes involved in pectin backbone degradation mainly included pectate lyases belonging to PL1, PL3, PL4, and PL9 families and hydrolases belonging to GH family 28 [28]. R. luteum lacked all GH28 genes for enzymes that degrade pectin backbones, such as endopolygalacturonase (PGA), exopolygalacturonase (PGX), endorhamnogalacturonase (RHG), exorhamnogalacturonase (RGX). Moreover, R. luteum had only one pectin lyase, Ramle7389 (PEL: PL1), which was expressed similarly in planta and in vitro. R. luteum grew more slowly on a medium containing apple pectin as the sole nutrient source than that containing glucose; this result correlated with a reduced repertoire for pectin degradation (Additional file 5: Figure S4). The ectophytic fungi P. fructicola had one enzyme (GH28) involved in pectin lyase backbones and Z. tritici had five genes (two GH28, two PL1 and one PL3), which was fewer than for symbiotic fungi (av. = 7) and penetrating plant pathogens (av. = 25) (Fig. 4).

Fig. 4
figure 4

Heatmap showing the distribution of putative plant cell wall degrading enzymes among different fungi. Numbers of family members in each genome are shown. Fungal names: BC = Botrytis cinerea; CH = Colletotrichum higginsianum; LB = Laccaria bicolor; MO = Magnaporthe oryzae; ZT = Zymoseptoria tritici; NP = Neofusicoccum parvum; NC = Neurospora crassa; PF = Peltaster fructicola; PT = Pyrenophora tritici-repentis; RL = Ramichloridium luteum; SS = Sclerotinia sclerotiorum; TR = Trichoderma reesei; CG = Cenococcum geophilum. Overrepresented (+3 to 0) and underrepresented (0 to −3) numbers are depicted as Z-scores for each line in heatmap. Boxes on the nethermost represent the life style of the selected organisms. E, ectophytes; H, hemibiotrophs; N, necrotrophs; S, saprotrophs; Y, symbionts. Boxes on the right represent the CAZY substrates. P, pectin; H, hemicellulose; C, cellulose

Full size image

R. luteum and P. fructicola were also reduced in regard to β-1,4-xylosidases (BXL: GH3, GH43) that release D-xyloses from non-reducing ends of xylooligosaccharides (Fig. 4). The numbers of GH3 and GH43 in R. luteum and P. fructicola were much lower than in Z. tritici or the average of cell-penetrating plant pathogens (10 and 9 vs. 26 and 31, respectively). Cellulose is a polymer of β-1,4-linked D-glucose residues, the degradation of which is mediated by three classes of enzymes in synergy: β-1,4-endoglucanases (EGL: GH5, 7 and 12), exoglucanases/cellobiohydrolases (CBH: GH6 and GH7), and β-glucosidase (BGL: GH1 and GH3) [28]. Cellulase genes of R. luteum, P. fructicola, Z. tritici and the average of penetrating plant pathogens numbered 22, 16, 30 and 62 respectively. Interestingly, R. luteum and P. fructicola processed only EGL and BGL, whereas Z. tritici processed all three classes.

Collectively, reduction in PCWDEs content appeared to be an ancient trait for the phylogenetic clade containing R. luteum, P. fructicola, and Z. tritici. However, the fact that R. luteum and P. fructicola contained many fewer PCWDEs than Z. tritici suggested a secondary round of PCWDE gene reduction associated with adaptation to an ectophytic niche. In addition, RNA-seq data showed that only two secreted PCWDE-coding genes, Ramle6973 and Ramle1916 (GH3, for releasing of terminal glucoses from shorter oligosaccharides), were moderately expressed and up-regulated in planta (Additional file 1: Table S5).

Expansion of cutinase and secretory lipase

Cutinase and secretory lipase are key enzymes degrading plant cutin and wax derivatives, and are presumably important for surface niche adaptations [29,30,31]. CAFE analysis showed that both categories of enzymes were significantly expanded in the R. luteum genome (P < 0.05) (Fig. 2). Of a total of 13 fungal species surveyed in this work, R. luteum, together with the rice blast fungus M. oryzae, encoded the highest number of putative cutinases (PF01083, 18 members), which was far more than the average number (9). The difference was even more striking in the Dothideomycetes, in which the previously highest documented number was 14 in Cochliobolus sativus and C. heterostrophus (http://pfam.xfam.org). Interestingly, this increase in cutinase copy number was not observed with the SBFS fungus P. fructicola (6) or in the non-SBFS Z. tritici (6). Phylogenetic analysis identified 11 clusters out of 112 putative cutinases (Fig. 5a). R. luteum cutinases were distributed among clusters 1, 6, 7, 8 and 11. Gene expansion was noted in clusters 1, 6 and 7 (Fig. 5b). Clusters 8 and 11 contained two species-specific cutinase genes, Ramle3391 and Ramle3581. Interestingly, cutinases Ramle3581 (cluster 11), Ramle3560 (cluster 7) and Ramle3575 (cluster 6) form a tandem repeat cluster (Fig. 5c), indicating local duplication events over the course of cutinase evolution in this species. Meanwhile, this cluster was surrounded by genes encoding α/β-hydrolase and fatty acid hydroxylase, and there was also a Gypsy retroelement (~2 kb in length) in the 5′-flanking region of Ramle3581 (<4 kb).

Fig. 5
figure 5

Striking expansion of putative cutinase genes in R. luteum. a phylogenetic tree of predicted cutinases among the 13 compared fungal species. b gene copy number variation among the phylogenetic clusters defined in A. Numbers of putative cutinase genes in each genome are shown. Overrepresented (+3 to 0) and underrepresented (0 to −3) numbers are depicted as Z-scores in heatmap. c local tandem duplication of cutinase genes in R. luteum. Botci1 = Botrytis cinerea; Maggr1 = Magnaporthe oryzae; Mycgr3 = Zymoseptoria tritici; Neopa1 = Neofusicoccum parvum; Pelfr = Peltaster fructicola; Ramle = Ramichloridium luteum; Sclsc = Sclerotinia sclerotiorum; Trire1 = Trichoderma reesei; Lacbi2 = Laccaria bicolor; Pyrtr1 = Pyrenophora tritici-repentis; Colhi1 = Colletotrichum higginsianum; Neucr2 = Neurospora crassa. E, ectophytes; H, hemibiotrophs; N, necrotrophs; S, saprotrophs; Y, symbionts

Full size image

HsbA is a secreted protein which facilitates the binding between cutinase and hydrophobic surfaces, thus facilitating cutin degradation [32]. The genome of R. luteum contained significantly more HsbA than penetrating plant pathogen species (14 vs. 1–7) (Fig. 2), and 13 of these proteins were predicted to be secreted.

The R. luteum genome encoded nine secreted lipases (domain PF03583), which is exceeds the number in P. fructicola and Z. tritici (4 and 2, respectively) and is significantly higher than the average number found in all fungal species in this study (2, ***, P < 0.001, t-test,). Phylogenetically, lipases from the 13 compared species were divided into six clusters (Fig. 6a), with clusters 3 and 5 being specific to the SBFS fungi, R. luteum and P. fructicola. Interestingly, genes in cluster 5 were most similar to Basidiomycota sequences based on a BLASTp search (Additional file 1: Table S6). Ramle0713 was selected for illustration. When the protein was searched against the NCBI nr database with BLASTp, 46 hits (more than 200 alignment scores), all belonging to Basidiomycota, were applicable for building phylogenies. The results of phylogenetic analysis, shown in Fig. 6b, strongly suggest that cluster 5 genes have been transferred from Basidiomycota through horizontal gene transfer.

Fig. 6
figure 6

Expansion of secretory lipases in R. luteum and their horizontal transfer from Basidiomycota. a phylogenetic tree of secretory lipases among the 13 compared fungal species. b phylogenetic tree constructed with R. luteum proteins in cluster 5 together with their best BLASTp hits in NCBI nr database. Botci1 = Botrytis cinerea; Maggr1 = Magnaporthe oryzae; Mycgr3 = Zymoseptoria tritici; Neopa1 = Neofusicoccum parvum; Pelfr = Peltaster fructicola; Ramle = Ramichloridium luteum; Sclsc = Sclerotinia sclerotiorum; Trire1 = Trichoderma reesei

Full size image

Based on RNA-seq data, two secretory lipases (Ramle1714 and Ramle9316) and one cutinase (Ramle6800) were highly expressed and up-regulated at least two-fold in planta. In addition, one HsbA homolog (Ramle4244) was significantly up-regulated by ~3-fold at the colonization stage. During ectophytic colonization of the apple fruit surface, we observed an apparent disappearance of waxy crystals around the hyphae network (Fig. 7a). Later, clearing zones devoid of waxy crystals started to form within the hyphae network (Fig. 7b). The microscopic observations suggested that R. luteum could degrade apple cuticle efficiently.

Fig. 7
figure 7

R. luteum could efficiently degrade apple cuticles by scanning electron microscopy (SEM) observation. (a) SEM microscopy showing the disappearance of waxy crystals around hyphal network. (b) SEM microscopy showing the formation of a clearing zone surrounding ectophytic hyphae. Scale bars: a = 10 μm; b = 20 μm

Full size image

Fungal surface adaptation and stress responses

Microorganisms on surface niches are directly exposed to atmospheric extremes and sunlight. Furthermore, the hydrophobic waxy cuticle that covers plant epidermal cells reduces evaporation of water inside the fruit and the leaching of substances from plants [33, 34], thus resulting in an oligotrophic environment [35]. Reproduction on the plant surface requires adaptation to this habitat, including the ability to withstand abiotic and biotic stresses such as scarce nutrients and water, UV irradiation, and antagonism by other microorganisms. The following transcriptomic analyses focus on expression of traits of R. luteum in its ectophytic niche.

a) Melanin Melanin can protect fungi from UV irradiation, desiccation, high temperatures, and strong oxidants [36]. On fruit surfaces, R. luteum hyphae exposed to air are heavily melanized. In most Ascomycota species, melanin is synthesized via the 1,8-DHN biosynthesis pathway [36]. The R. luteum genome encodes a complete 1,8-DHN biosynthetic pathway, which includes PKS1 (Ramle2903), 1,3,6,8-tetrahydroxynaphthalene reductase (Ramle8102), scytalone dehydratase (Ramle4854), and 1,3,8-trihydroxynaphthalene reductase (Ramle6035). RNA-Seq data indicated that these genes were expressed similarly in planta and in vitro, with PKS1 expressed at an intermediate level whereas the other three genes were expressed at a high level (Fig. 8).

Fig. 8
figure 8

Proposed surface niche adaptation model for the ectophytic SBFS fungus R. luteum. Plant cell wall degrading enzymes (PCWDEs), PTH11-like G protein-coupled receptors (GPCRs) and effectors were drastically reduced. In contrast, genes encoding cutinases, secretory lipases and GH25-lysozymes were strikingly expanded. Genes involved in cuticle degradation (cutinase, secretory lipase) and stress responses (melanin biosynthesis, aquaporins, lysozymes and HOG pathway) were activated by R. luteum colonizing in its ectophytic niche. For heatmaps, two columns represent different treatments, i.e., inoculation on apple fruit (in planta) and growth on artificial media (in vitro), and each row is marked with the name of one gene. The colored scale bar of expression levels is divided into three grades: low (0 < FPKM <10, including 0), medium (10 < FPKM <80), and high (80 < FPKM)

Full size image

b) Lysozyme The CAZY GH25 family exhibits lysozyme-mediated lytic activity, which attacks peptidoglycans in the cell walls of several bacterial species [37]. Interestingly, the R. luteum genome encoded three GH25 proteins whereas other fungal genomes encode one at most. All three lysozymes (Ramle1917, Ramle4441 and Ramle5289) were up-regulated in planta, increasing 8-, 15- and 6-fold, respectively (Fig. 8). These putative lysozymes may give the R. luteum a competitive advantage over bacteria on the fruit surface.

c) Aquaporin Aquaporins are integral membrane proteins which form water-transporting pores on membrane surface and play critical roles in osmotic stress responses [38]. The R. luteum genome encoded three aquaporin-like proteins, two of which (Ramle5938 and Ramle6492) were up-regulated in planta, with the induction fold being 3 and 25 respectively (Fig. 8).

d) High osmolarity glycerol (HOG) pathway. The HOG pathway is an important mitogen-activated protein kinase (MAPK) signaling pathway that is able to detect, transduce signals and respond to osmotic and oxidative stresses [39, 40]. The HOG MAPK pathway was complete in the R. luteum genome (Additional file 1: Table S7). The MAPK kinase Pbs2p (Ramle3536) and the downstream bZip transcription factor AtfA (Ramle1430) were up-regulated by approximately two- and three-fold, respectively, in planta compared to in vitro. Several putative AtfA downstream targets were strongly up-regulated in planta; examples include three heat shock proteins (Ramle2917, Ramle1564 and Ramle5908, exhibiting fold increases of 718, 74, and 37, respectively), antioxidative genes like glutathione peroxidase (Ramle0653, 2-fold increase), and two catalases (Ramle0422 and Ramle1388, fold increases of 143 and 37, respectively) (Fig. 8).

The above results suggested that melanin, the HOG pathway and aquaporin play key roles in protecting ectophytic fungi from UV irradiation, desiccation, high temperatures and oxidants, and that putative lysozymes may help them to compete with other microorganisms. For adaptation to an oligotrophic environment, greatly expanded cutinase and lipase production could increase the efficiency of degradation of the hydrophobic waxy cuticle to enhance access to sugar-rich fruit leachates.

Discussion

The plant cuticle is an extracellular lipophilic biopolymer made up mainly of cutin and epicuticular wax [41, 42]. Cutinase hydrolyzes primarily cutin whereas lipase acts mainly on epicuticular waxes [29]. Our comparative genomic analysis revealed that cutinase and secretory lipase are strikingly expanded in the R. luteum genome, which is in accordance with our microscopic observation that this fungus can actively degrade the cuticle during apple surface colonization. The R. luteum genome contains three cutinases forming tandem duplications, which occur near (<4 kb) the repeat retrotransposon Gypsy, indicating that local gene duplication contributes to gene family expansion. Furthermore, based on our phylogenetic analysis R. luteum may have acquired secretory lipases from Basidiomycota fungi via horizontal gene transfer.

In contrast with R. luteum, histological studies of the SBFS fungi P. fructicola, Leptodontidium elatius, Zygophiala jamaicensis and Z. wisconsinensis suggested that these fungi do not efficiently degrade the wax and cuticle layers during colonization [6, 43,44,45]. It is therefore possible that different SBFS fungi have evolved different capacities for degrading plant cuticles despite the fact that they are all well-adapted to the surface niche. Moreover, the R. luteum and P. fructicola genomes differ strikingly in the number of putative cutinases (18 vs. 6) and secretory lipases (9 vs. 4).

On the other hand, despite their long evolutionary separation (over 200 mya), the R. luteum and P. fructicola genomes share many signatures. For example, both genomes were drastically reduced in size, repeat elements, genes involved in plant cell wall degradation, PTH11-like GPCRs and effectors in comparison to host cell-penetrating species. The substantial divergence in phylogeny, together with similarity in genome reduction, indicate that their genome reductions are the result of convergent evolution. The most closely related penetrating phytopathogenic species, Z. tritici, also showed reduction in these features when compared with other fungi. Thus, the reductions are phylogenetically related and likely are results of contractions from the genomes their most recent common ancestor. On the other hand, even when compared with Z. tritici, the reduction of PCWDEs, effectors, and PTH11-like GPCRs was still substantial in both R. luteum and P. fructicola. This suggests that there was a subsequent round of virulence-related gene content reduction in the two SBFS fungi that was associated with niche adaptation.

A reduced complement of PCWDEs is meaningful for pathogenic and mutualistic interactions of biotrophs and hemibiotrophs to avoid triggering plant defense mechanisms and establish a biotrophic lifestyle [2, 46]. Ectomycorrhizal fungi colonizing the extracellular spaces also benefit from PCWDEs to penetrate the roots, although they exhibited a general reduction of such genes [46, 47]. Interestingly, numbers of PCWDEs in R. luteum and P. fructicola were fewer than for the Dothidiomycete species C. geophilum and even the symbiont Basidiomycete L. bicolor. Even when compared to other ectomycorrhizal fungi in Agaricomycetes [48], the number of CAZY associated with plant cell wall degradation are fewer in the SBFS species. This result may indicate evolutionary adaptations to the ectophytic lifestyle; i.e., colonizing epicuticular layers of plants and extracting nutrients from the surface of living hosts without penetrating and killing cells.

Our data also highlight the importance of stress response genes for ectophytic colonization by R. luteum. On the plant surface, microorganisms must deal with ultraviolet radiation and low or fluctuating water availability. Moreover, microorganisms must compete with other surface-inhabiting organisms. In bacteria, protection mechanisms that were shown experimentally to be important in epiphytic fitness involve the production of pigments and the activation of DNA repair mechanisms, catalases and superoxide dismutases [35, 49,50,51]. In the present study, we found a high level of expression of genes involved in melanin biosynthesis, aquaporins, heat shock proteins, catalases and peroxidase processing by R. luteum. Interestingly, R. luteum expanded three CH25-lysozymes during its evolution, which is more than any of the other 12 fungal species we used for comparison, suggesting a unique mechanism to deal with biotic stress at the plant surface.

Conclusion

In this study, we reported the genome sequence of Ramichloridium luteum, one of the representative SBFS species on apple in China [13]. Comparing the genome and transcriptome of this ectophytic species with P. fructicola, a first model of surface colonization, as well as those of fungi that occupy other ecological niches, has yielded new insight into the evolutionary adaptability and life strategy of ectophytic fungi. R. luteum and P. fructicola share many genomic similarities related to surface niche adaptation. For instance, compared with penetrating plant-pathogenic species and penetrating symbiont fungi, both genomes showed a dramatic reduction in pathogenicity-related genes. Our data supports the ectophytic phytopathogen evolutionary model proposed by Xu et al. [6], a pathway from interior to exterior colonizers. In fact, colonization by stealth is strategic for an ectophyte in order to reduce the risk of triggering defensive countermeasures by the host. Moreover, both SBFS fungi share similarity in genes controlling 1,8-DHN-melanin biosynthesis, which were expressed at a high level on the plant surface. On the other hand, R. luteum has strikingly more genes encoding cutinase and secretory lipase than P. fructicola, suggesting a scenario whereby different SBFS fungi have evolved distinct strategies to ensure successful surface niche adaptation. Generating and analyzing additional SBFS genomes will be necessary for discerning common genome features associated with their surface niche adaptations.

Methods

Fungal strain and culture conditions

Strain CPC 18961 of R. luteum was originally isolated from the surface of an apple fruit showing SBFS symptoms near Weifang City, Shandong Province, China. The culture was single-spore purified, maintained on potato dextrose agar (PDA) medium at 25°C, and stored in glycerol (15%) at -80°C in the Fungal Laboratory of Northwest A&F University, Yangling, Shaanxi Province, China.

DNA isolation, genome sequencing and assembly

Highly purified total genomic DNA was isolated from the fungal mycelia collected from a 2-week old PDA culture following the modified cetyltrimethyl ammonium bromide protocol [52]. The genome of R. luteum was sequenced with the Illumina HiSeq™ 2000 platform. The insertion size of the sequencing library was 360 bp and the sequencing strategy was 100 bp pair-ends. Filtered clean reads were assembled into scaffolds using the ABySS assembler v. 1.3.5 [53]. The completeness of assembly was assessed using CEGMA v. 2.4 [54]. GapFiller software was used to further fill the gaps and generate scaffolds [55].

Gene prediction and genome annotation

GeneMark-ET was first used to predict gene structures based on RNA-seq data [56]. Obtained gene models were used to train Augustus v. 3.1 [57]. Predicted gene models from GeneMark-ET and Augustus, as well as the homology proteins of Capnodiales, were combined in MAKER2 [58]. InterProScan and Pfam were used to identify protein domains [59, 60]. Protein encoding genes were annotated through BLASTp search against the GO [61], COG [62], KEGG [63], and non-redundant NCBI databases, at the threshold of E-value ≤1 × 10−5. Repeat sequences were identified by RepeatMasker v. 4.0.5 (http://www.repeatmasker.org) and RepeatModeler v. 1.0.7 [64], and non-coding RNA was predicted by rRNAmmer v. 1.2 (http://www.cbs.dtu.dk/services/RNAmmer), tRNAscan-SE v. 1.3.1 [65] and Rfam (http://rfam.xfam.org). For calculation of RIP indices, dinucleotide frequencies were determined using the RIPCAL program [66].

Functional annotation of predicted genes

Genes encoding putative CAZY were identified using the hmmscan program by searching the fungal proteomes with the family-specific HMM profiles of CAZY downloaded from the dbCAN database [67]. CAZY families corresponding to PCWDEs were in accordance with previous research [68]. Melanin synthase genes were identified by searching with Cochliobolus heterostrophus homologs [69]. Stress-responsive proteins were identified by BLASTp search with queries from other fungal species such as Saccharomyces cerevisiae, S. pombe and Aspergillus nidulans [40].

Candidate secreted proteins have a secretion signal as determined by SignalP v. 4.1 [70] and have no transmembrane domain as determined by TMHMM 2.0 (http://www.cbs.dtu.dk/services/TMHMM), and finally having extracellular score > 15 as predicted by WoLF-PSort v. 0.2 software [71]. Small secreted proteins are defined here as proteins that are smaller than 200 aa and were labeled as ‘cysteine rich’ when the percentage of cysteine residue was at least twice as high as the average percentage in entire proteomes [72].

Phylogenomic analysis

OrthoMCL v. 2.0.9 was used to identify ortholog pairs among compared genomes. The cutoff E-value was set as 1 × 10−5 [73]. To construct a genome-based phylogenic tree, single-copy ortholog pairs were aligned with MAFFT v. 7 (http://mafft.cbrc.jp/alignment/server), conserved sites in the alignments were further extracted with Gblocks v. 0.91b using the default parameters [74], and the dataset was used for maximum likelihood tree construction in RAxML [75] with the LG + I + G + F amino acid substitution model selected by ProtTest v. 3.4 [76]. Divergence times between species were estimated using the PL method with r8s [77]. The estimated Sordariomycota divergence (290–380 million years ago) was used for calibration [78]. MEGA v. 7.0 was used to generate maximum likelihood phylogeny for cutinase and secretory lipase with the JTT (Jones, Taylor, and Thorton) amino acid substitution model [79, 80]. Statistical support for phylogenetic grouping was assessed by 1000 bootstrap re-samplings. The program CAFE (Computational Analysis of gene Family Evolution) v. 3 [81] was used to analyze gene family size evolution, and the Viterbi algorithm in the CAFE program was used to identify families going through significant expansions/contractions at each branch using a cutoff of P < 0.05.

Inoculation experiment

The R. luteum isolate CPC 18961 was used to inoculate apple fruit in the field following the described protocol [9]. Prior to inoculation, single-conidium-derived isolates were cultured on PDA for 1 month. Five mycelial colonies were harvested, excessive agar was cut away, colonies were transferred to five 1.5-ml plastic centrifuge tubes (one colony per tube) to which 600 μl sterile deionized water (SDW) was added, and the mixture was shaken on a vortex oscillator (Model QL-901, Kylin-Bell Lab Instruments Company Limited, Haimen City, Jiangsu Province, China) for 60 s. After that, they were refrigerated until use, which was within 2 h of preparation [45].

Immature apple fruit (cv. Fuji; 5–6 cm in diameter) were selected for in situ inoculation. Fruit surfaces were sterilized by swabbing with 70% ethanol, air-dried, and then swabbed with inoculum suspension using sterilized brushes. In total, 10 fruit were inoculated and 10 fruit were treated only with SDW as a control. All 20 fruits were placed in sealed transparent polyethylene bags (20 × 15 cm) (The CLOROX Company, Oakland, California); two corners (1.5 × 1.5 cm) of each bag were cut off each bag to facilitate aeration. After approximately 3 months of incubation (temperature range from 11o C to 37°C), apples with visible colonies of R. luteum were harvested for total RNA extraction and microscopic examination.

Transcriptome analysis

Mycelia from artificial inoculations onto apple fruit in the field were used as ‘in planta’ samples for RNA sequencing. For ‘in vitro’ mycelium, R. luteum was cultured in liquid Czapek-Dox medium [82] for 7 d at 25°C with a shake speed of 150 rpm. For ‘in planta’ and ‘in vitro’ treatments, three independent biological samples were collected for RNA extraction and sequencing. The average cDNA library was 350 bp. The libraries were paired-end sequenced on a flow cell using an Illumina HiSeq™ 2500 sequencing platform with a read length of 125 bp. After filtration, reads were mapped to the genome sequence using TopHat v. 2.0.9 with the parameter “-g 1” [83]. Reads mapped to multiple locations were discarded. In total, 29,060,627 RNA-seq reads were mapped to the R. luteum genome (Additional file 1: Table S8). A gene expression profile was created using Cufflinks v. 2.2.1 [84]. FPKM (fragments per kilobase of transcript per million fragments mapped) values were used to determine Pearson’s correlation coefficient among biological replicates (0.84 to 0.99) (Additional file 4: Figure S3a). We performed principal component analysis (PCA) to assess biological variability among all samples (Additional file 4: Figure S3b); FPKM of transcripts with a significant P-value (<0.05) and a greater than twofold change (log2) in transcript level were considered to be differently expressed. All P-values were corrected for false discoveries resulting from multiple hypothesis testing using the Benjamini-Hochberg procedure.

Scanning electron microscopy

Epicarps with sooty blotch colonies were peeled off the fruit and cut into 5 mm × 5 mm pieces, then air dried for 48 h. Dehydrated samples were mounted on stubs, sputter coated with gold-palladium [85] and viewed using a HITACHI S-3400N scanning electron microscopy operating at 5 KV.

Apple pectin utilization

For determining physiological correlates of gene reduction associated with pectin degradation, we grew the SBFS fungi on media containing pectin as a sole source of carbon. We used water-agar (Bacto™ Agar, USA) medium as the basal medium, amended with 0.2% apple pectin (Sigma-Aldrich, USA) as a sole carbon source and 0.2% glucose as a sole carbon source (Guangdong Guanghua Sci-Tech Co. Ltd., China), respectively. For strain growth, 1 μl of conidial suspension (concentration: 1 × 10−6 CFU/ml) was inoculated onto each plate; the culture was incubated at 25°C for 10 d and then photographed. Water-agar medium was set as the internal control medium. There were six plates per treatment.

Abbreviations

1,8-DHN:

1,8-dihydroxynaphthalene

BLAST:

Basic Local Alignment Search Tool

CAFE:

Computational Analysis of gene Family Evolution

COG:

Clusters of Orthologous Groups

FPKM:

Fragments Per Kilobase of transcript per Million fragments mapped

GPCRs:

G Protein-Coupled Receptors

HOG:

High Osmolarity Glycerol

HsbA:

Hydrophobic surface binding protein A

KEGG:

Kyoto Encyclopaedia of Genes and Genomes

MAPK:

Mitogen-Activated Protein Kinase

MEGA:

Molecular Evolutionary Genetic Analysis

PCWDEs:

Plant Cell Wall Degrading Enzymes

RIP:

Repeat-Induced Point mutation

SBFS:

Sooty Blotch and FlySpeck

SEM:

Scanning Electron Microscopy

References

  1. Yu X, Lund SP, Scott RA, Greenwald JW, Records AH, Nettleton D, et al. Transcriptional responses of Pseudomonas syringae to growth in epiphytic versus apoplastic leaf sites. Proc Natl Acad Sci U S A. 2013;110:E425–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Spanu PD. The genomics of obligate (and nonobligate) biotrophs. Annu Rev Phytopathol. 2012;50:91–109.

    Article  CAS  PubMed  Google Scholar 

  3. Fesel PH, Zuccaro A. Dissecting endophytic lifestyle along the parasitism/mutualism continuum in Arabidopsis. Curr Opin Microbiol. 2016;32:103–12.

    Article  PubMed  Google Scholar 

  4. Thomas J, Mathew K. Biochemical changes in the Lawsonia inermis L. infected with Asterina lawsoniae Henn. & Nyn. Curr Res Environ Appl Mycol. 2014;4:217–20.

    Google Scholar 

  5. Flores-Palacios A. Does structural parasitism by epiphytes exist? A case study between Tillandsia recurvata and Parkinsonia praecox. Plant Biol. 2016;18:463–70.

    Article  CAS  PubMed  Google Scholar 

  6. Xu C, Chen H, Gleason ML, Xu JR, Liu HQ, Zhang R, et al. Peltaster fructicola genome reveals evolution from an invasive phytopathogen to an ectophytic parasite. Sci Rep. 2016;6:22926.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ismail SI, Batzer JC, Harrington TC, Crous PW, Lavrov DV, Li HY, et al. Ancestral state reconstruction infers phytopathogenic origins of sooty blotch and flyspeck fungi on apple. Mycologia. 2016;108:292–302.

    Article  PubMed  Google Scholar 

  8. Gleason ML, Batzer JC, Sun G, Zhang R, Arias MMD, Sutton TB, et al. A new view of sooty blotch and flyspeck. Plant Dis. 2011;95:368–83.

    Article  Google Scholar 

  9. Batzer JC, Gleason ML, Harrington TC, Tiffany LH. Expansion of the sooty blotch and flyspeck complex on apples based on analysis of ribosomal DNA gene sequences and morphology. Mycologia. 2005;97:1268–86.

    Article  CAS  PubMed  Google Scholar 

  10. Arias MMD, Batzer JC, Harrington TC, Wong AW, Bost SC, Cooley DR, et al. Diversity and biogeography of sooty blotch and flyspeck fungi on apple in the eastern and midwestern United States. Phytopathology. 2010;100:345–55.

    Article  Google Scholar 

  11. Yang HL, Sun GY, Batzer JC, Crous PW, Groenewald JZ, Gleason ML. Novel fungal genera and species associated with the sooty blotch and flyspeck complex on apple in China and the USA. Persoonia. 2010;24:29–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Li HY, Sun G, Batzer JC, Crous PW, Groenewald JZ, Karakaya A, et al. Scleroramularia gen. nov. associated with sooty blotch and flyspeck of apple and pawpaw from the Northern Hemisphere. Fungal Divers. 2011;46:53–66.

    Article  Google Scholar 

  13. Li HY, Sun GY, Zhai XR, Batzer JC, Mayfield DA, Crous PW, et al. Dissoconiaceae associated with sooty blotch and flyspeck on fruits in China and the United States. Persoonia. 2012;28:113–25.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Mayfield DA, Karakaya A, Batzer JC, Blaser JM, Gleason ML. Diversity of sooty blotch and flyspeck fungi from apples in northeastern Turkey. Eur J Plant Pathol. 2013;135:805–15.

    Article  Google Scholar 

  15. Huang M, Sanchez-Moreiras AM, Abel C, Sohrabi R, Lee S, Gershenzon J, et al. The major volatile organic compound emitted from Arabidopsis thaliana flowers, the sesquiterpene (E)-β-caryophyllene, is a defense against a bacterial pathogen. New Phytol. 2012;193:997–1008.

    Article  CAS  PubMed  Google Scholar 

  16. Kershaw MJ, Talbot NJ. Genome-wide functional analysis reveals that infection-associated fungal autophagy is necessary for rice blast disease. Proc Natl Acad Sci U S A. 2009;106:15967–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Arzanlou M, Groenewald JZ, Gams W, Braun U, Shin HD, Crous PW. Phylogenetic and morphotaxonomic revision of Ramichloridium and allied genera. Stud Mycol. 2007;58:57–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Sydow H, Sydow P. Beitrag zur Kenntniss der Pilzflora der Philippinen-Inseln. Ann Mycol. 1917;15:165–268.

    Google Scholar 

  19. Mohanta TK, Bae H. The diversity of fungal genome. Biol Proced Online. 2015;17:1.

    Article  CAS  Google Scholar 

  20. DeZwaan TM, Carroll AM, Valent B, Sweigard JA. Magnaporthe grisea Pth11p is a novel plasma membrane protein that mediates appressorium differentiation in response to inductive substrate cues. Plant Cell. 1999;11:2013–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kulkarni RD, Thon MR, Pan H, Dean RA. Novel G-protein-coupled receptor-like proteins in the plant pathogenic fungus Magnaporthe grisea. Genome Biol. 2005;6:R24.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Stergiopoulos I, de Wit PJ. Fungal effector proteins. Annu Rev Phytopathol. 2009;47:233–63.

    Article  CAS  PubMed  Google Scholar 

  23. Plett JM, Martin F. Blurred boundaries: lifestyle lessons from ectomycorrhizal fungal genomes. Trends Genet. 2011;27:14–22.

    Article  CAS  PubMed  Google Scholar 

  24. Rovenich H, Boshoven JC, Thomma BP. Filamentous pathogen effector functions: of pathogens, hosts and microbiomes. Curr Opin Plant Biol. 2014;20:96–103.

    Article  CAS  PubMed  Google Scholar 

  25. Nunes-Alves C. Bacterial secretion: A trans-kingdom effector. Nat Rev Microbiol. 2014;12:460–1.

    Article  CAS  PubMed  Google Scholar 

  26. Underwood W. The plant cell wall: a dynamic barrier against pathogen invasion. Front Plant Sci. 2012;3:85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Goodwin SB, M'barek SB, Dhillon B, Wittenberg AH, Crane CF, Hane JK, et al. Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet. 2011;7:e1002070.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. van den Brink J, de Vries RP. Fungal enzyme sets for plant polysaccharide degradation. Appl Microbiol Biotechnol. 2011;91:1477–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Feng J, Liu G, Selvaraj G, Hughes GR, Wei Y. A secreted lipase encoded by LIP1 is necessary for efficient use of saturated triglyceride lipids in Fusarium graminearum. Microbiology. 2005;151:3911–21.

    Article  CAS  PubMed  Google Scholar 

  30. Feng J, Wang F, Liu G, Greenshields D, Shen W, Kaminskyj S, et al. Analysis of a Blumeria graminis-secreted lipase reveals the importance of host epicuticular wax components for fungal adhesion and development. Mol Plant Microbe In. 2009;22:1601–10.

    Article  CAS  Google Scholar 

  31. Nawrath C. Unraveling the complex network of cuticular structure and function. Curr Opin Plant Biol. 2006;9:281–7.

    Article  CAS  PubMed  Google Scholar 

  32. Ohtaki S, Maeda H, Takahashi T, Yamagata Y, Hasegawa F, Gomi K, et al. Novel hydrophobic surface binding protein, HsbA, produced by Aspergillus oryzae. Appl Environ Microb. 2006;72:2407–13.

    Article  CAS  Google Scholar 

  33. Tukey HB. The leaching of substances from plants. Annu Rev Plant Physiol. 1970;21:305–24.

    Article  CAS  Google Scholar 

  34. Schönherr J. Resistance of plant surfaces to water loss: transport properties of cutin, suberin and associated lipids. In: Lange OL, Nobel PS, Osmond CB, Ziegler H, editors. Encyclopedia of Plant Physiology. Berlin: Springer; 1982. p. 153–79.

    Google Scholar 

  35. Vorholt JA. Microbial life in the phyllosphere. Nat Rev Microbiol. 2012;10:828–40.

    Article  CAS  PubMed  Google Scholar 

  36. Schumacher J. DHN melanin biosynthesis in the plant pathogenic fungus Botrytis cinerea is based on two developmentally regulated key enzyme (PKS)-encoding genes. Mol Microbiol. 2016;99:729–48.

    Article  CAS  PubMed  Google Scholar 

  37. Moraïs S, Cockburn DW, Ben-David Y, Koropatkin NM, Martens EC, Duncan SH, et al. Lysozyme activity of the Ruminococcus champanellensis cellulosome. Environ Microbiol. 2016;18:5112–22.

    Article  PubMed  Google Scholar 

  38. Alleva K, Chara O, Amodeo G. Aquaporins: another piece in the osmotic puzzle. FEBS Lett. 2012;586:2991–9.

    Article  CAS  PubMed  Google Scholar 

  39. Zhao X, Mehrabi R, Xu JR. Mitogen-activated protein kinase pathways and fungal pathogenesis. Eukaryot Cell. 2007;6:1701–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Miskei M, Karányi Z, Pócsi I. Annotation of stress-response proteins in the aspergilli. Fungal Genet Biol. 2009;46:S105–20.

    Article  CAS  PubMed  Google Scholar 

  41. Kolattukudy PE. Polyesters in higher plants. Adv Biochem Eng Biotechnol. 2001;71:1–49.

    CAS  PubMed  Google Scholar 

  42. Samuels L, Kunst L, Jetter R. Sealing plant surfaces: cuticular wax formation by epidermal cells. Annu Rev Plant Biol. 2008;59:683–707.

    Article  CAS  PubMed  Google Scholar 

  43. Nasu H, Kunoh H. Scanning electron microscopy of flyspeck of apple, pear, Japanese persimmon, plum, Chinese quince, and pawpaw. Plant Dis. 1987;71:361–4.

    Article  Google Scholar 

  44. Belding RD, Sutton TB, Blankenship SM, Young E. Relationship between apple fruit epicuticular wax and growth of Peltaster fructicola and Leptodontidium elatius, two fungi that cause sooty blotch disease. Plant Dis. 2000;84:767–72.

    Article  CAS  Google Scholar 

  45. Gao L, Sun G, Zhang R, Gleason ML. Secondary spread of Zygophiala wisconsinensis on the surface of apple fruit. Eur J Plant Pathol. 2014;139:117–24.

    Article  Google Scholar 

  46. Kohler A, Kuo A, Nagy LG, Morin E, Barry KW, Buscot F, et al. Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat Genet. 2015;47:410–5.

    Article  CAS  PubMed  Google Scholar 

  47. Peter M, Kohler A, Ohm RA, Kuo A, Krützmann J, Morin E, et al. Ectomycorrhizal ecology is imprinted in the genome of the dominant symbiotic fungus Cenococcum geophilum. Nat Commun. 2016;7:12662.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Floudas D, Binder M, Riley R, Barry K, Blanchette RA, Henrissat B, et al. The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. Science. 2012;336:1715–9.

    Article  CAS  PubMed  Google Scholar 

  49. Jacobs JL, Carroll TL, Sundin GW. The role of pigmentation, ultraviolet radiation tolerance, and leaf colonization strategies in the epiphytic survival of phyllosphere bacteria. Microb Ecol. 2005;49:104–13.

    Article  CAS  PubMed  Google Scholar 

  50. Gunasekera TS, Sundin GW. Role of nucleotide excision repair and photoreactivation in the solar UVB radiation survival of Pseudomonas syringae pv. syringae B728a. J Appl Microbiol. 2006;100:1073–83.

    Article  CAS  PubMed  Google Scholar 

  51. Gourion B, Rossignol M, Vorholt JA. A proteomic study of Methylobacterium extorquens reveals a response regulator essential for epiphytic growth. Proc Natl Acad Sci U S A. 2006;103:13186–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8:4321–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23:1061–7.

    Article  CAS  PubMed  Google Scholar 

  55. Nadalin F, Vezzi F, Policriti A. GapFiller: a de novo assembly approach to fill the gap within paired reads. BMC Bioinform. 2012;13:S8.

    Article  Google Scholar 

  56. Lomsadze A, Burns PD, Borodovsky M. Integration of mapped RNA-Seq reads into automatic training of eukaryotic gene finding algorithm. Nucleic Acids Res. 2014;42:e119.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Stanke M, Morgenstern B. AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints. Nucleic Acids Res. 2005;33(Suppl 2):W465–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Cantarel BL, Korf I, Robb SM, Parra G, Ross E, Moore B, et al. MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2008;18:188–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Zdobnov EM, Apweiler R. InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17:847–8.

    Article  CAS  PubMed  Google Scholar 

  60. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42:D222–30.

    Article  CAS  PubMed  Google Scholar 

  61. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinform. 2003;4:41.

    Article  Google Scholar 

  63. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nuclei Acids Res. 1999;27:29–34.

    Article  CAS  Google Scholar 

  64. Saha S, Bridges S, Magbanua ZV, Peterson DG. Empirical comparison of ab initio repeat finding programs. Nucleic Acids Res. 2008;36:2284–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Hane J, Oliver R. RIPCAL: a tool for alignment-based analysis of repeat-induced point mutations in fungal genomic sequences. BMC Bioinform. 2008;9:478.

    Article  Google Scholar 

  67. Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40:W445–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Gan P, Ikeda K, Irieda H, Narusaka M, O'Connell RJ, Narusaka Y. Comparative genomic and transcriptomic analyses reveal the hemibiotrophic stage shift of Colletotrichum fungi. New Phytol. 2013;197:1236–49.

    Article  CAS  PubMed  Google Scholar 

  69. Eliahu N, Igbaria A, Rose MS, Horwitz BA, Lev S. Melanin biosynthesis in the maize pathogen Cochliobolus heterostrophus depends on two mitogen-activated protein kinases, Chk1 and Mps1, and the transcription factor Cmr1. Eukaryot Cell. 2007;6:421–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.

    Article  CAS  PubMed  Google Scholar 

  71. Paul H, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, et al. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35:585–7.

    Article  Google Scholar 

  72. Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, Barry KW, et al. Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 2012;8:e1003037.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Castresama J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540–52.

    Article  Google Scholar 

  75. Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.

    Article  CAS  PubMed  Google Scholar 

  76. Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27:1164–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Taylor JW, Berbee ML. Dating divergences in the Fungal Tree of Life: review and new analyses. Mycologia. 2006;98:838–49.

    Article  PubMed  Google Scholar 

  78. Lücking R, Huhndorf S, Pfister DH, Plata ER, Lumbsch HT. Fungi evolved right on track. Mycologia. 2009;101:810–22.

    Article  PubMed  Google Scholar 

  79. Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Bioinformatics. 1992;8:275–82.

    Article  CAS  Google Scholar 

  80. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    Article  CAS  PubMed  Google Scholar 

  81. Han MV, Thomas GWC, Lugo-Martinez J, Hahn MW. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013;30:1987–97.

    Article  CAS  PubMed  Google Scholar 

  82. Raper KB, Thom C. A manual of the penicillia. Maryland: Williams & Wilkins Press; 1949. p. 875.

    Google Scholar 

  83. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Kang ZS, Zingen-Sell I, Buchenauer H. Infection of wheat spikes by Fusarium avenaceum and alterations of cell wall components in the infected tissue. Eur J Plant Pathol. 2005;111:19–28.

    Article  Google Scholar 

Download references

Acknowledgments

We would like to thank the anonymous reviewers for their kind and helpful comments on the original manuscript.

Funding

This work was supported by National Natural Science Foundation of China (31371887), the 111 Project from Education Ministry of China (B07049), and China Agriculture Research System (CARS-28).

Availability of data and materials

This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession MTSC00000000 (version MTSC01000000). The complete transcriptome data sets are available at the Gene Expression Omnibus at the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/geo) and derived data sets as series accession number GSE100245.

Author information

Authors and Affiliations

  1. State Key Laboratory of Crop Stress Biology in Arid Areas and College of Plant Protection, Northwest A&F University, Yangling, Shaanxi Province, 712100, China

    Bo Wang, Xiaofei Liang, Rong Zhang & Guangyu Sun

  2. Department of Plant Pathology and Microbiology, Iowa State University, Ames, IA, 50011, USA

    Mark L. Gleason

Authors
  1. Bo Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiaofei Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Mark L. Gleason
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Rong Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Guangyu Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

GS initiated and coordinated the project. BW, GS and RZ conceived and planned the research. BW performed bioinformatics analysis and microscopic observation. BW wrote the paper. MLG, XL and GS revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Guangyu Sun.

Ethics declarations

Ethics approval and consent to participate

Strain CPC 18961 of R. luteum in this study was originally isolated from the surface of an apple fruit showing SBFS symptoms near Weifang City, Shandong Province, China. Inoculate apple fruits were collected from an orchard in Yangling, Shaanxi Province, China. Collecting such samples are in accordance with local regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Genome characteristics of R. luteum. Table S2. Species used in this comparative study. Table S3. The expression of 59 candidate effectors in R. luteum. Green fills represent the up-regulated genes on ectophytic growth. Table S4. Number of carbohydrate-active enzyme modules of R. luteum and 12 other fungi according to the CAZY database. GHs = Glycosyl hydrolases; GTs = Glycosyl transferases; CBMs = Carbohydrate-binding modules; CEs = Carbohydrate esterase; AAs = Auxiliary activities; PLs = Polysaccharide lyases. Table S5. Up-regulated genes involved in PCWDEs of R. luteum on the apple surface. Table S6. Nine secretory lipases show best-hit relationships with the fungi in NCBI nr database. Table S7. Orthologues and putative orthologues in response regulator-dependent osmotic and oxidative stress-response regulatory in R. luteum. Table S8. General features of data got by the high throughput transcriptome analysis. (XLS 72 kb)

Additional file 2: Figure S1.

GO functional classification of the R. luteum genome. In total, there are 5816 genes (61.44%) that have functional assignments. (TIFF 1334 kb)

Additional file 3: Figure S2.

Distribution of COG function annotation of the R. luteum genome. In total, there are 5857 genes (61.87%) that have functional assignments. (TIFF 700 kb)

Additional file 4: Figure S3.

Relationships of biological samples in RNA-seq experiment. (a) Pearson’s correlation outcomes; (b) PCA plot. (TIFF 543 kb)

Additional file 5: Figure S4.

Comparative growth profiling of R. luteum on two carbon substrates. (a) On water-agar medium; (b) On water-agar medium containing apple pectin as a sole carbon source; (c) On water-agar medium containing glucose as a sole carbon source. (TIFF 5055 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, B., Liang, X., Gleason, M.L. et al. Genome sequence of the ectophytic fungus Ramichloridium luteum reveals unique evolutionary adaptations to plant surface niche. BMC Genomics 18, 729 (2017). https://doi.org/10.1186/s12864-017-4118-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-4118-3

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us