BACKGROUND
arbuscular mycorrhizal fungi  live in an obligate symbiosis with the roots of around two-thirds of all studied land plant species, i.e.
angiosperms, gymnosperms, pteridophytes and some bryophytes  <cit> , and in one known case with cyanobacteria  <cit> .
this interaction has a major impact on the entire soil ecosystem and plays a crucial role in agricultural systems by increasing plant tolerance to biotic and abiotic stresses  <cit> .
in the field, amf form mycelial networks providing extensive connections between roots of different plant species and mediating interactions among plants.
these networks significantly improve the rhizosphere soil stability  <cit> .
the major benefit for plants upon mycorrhization is a more efficient nutrient uptake from the soil .
in exchange, amf obtain photoassimilates from their hosts  <cit> .
all amf belong to the phylum glomeromycota which has been divided into four orders: the glomerales, the diversisporales, and the two ‘ancestral’ lineages, paraglomerales and archaeosporales  <cit> .
they may have started to diverge over  <dig> mya.
the earliest spore fossils have been found together with the first land plants  and resemble present amf morphological structures  <cit> .
it has been suggested that amf played a crucial role for the adaptation of phototrophs to the terrestrial environment  <cit> .
the am symbiosis persisted morphologically unchanged throughout the complete evolutionary development within the plant phylum from haploid gametophytes to diploid sporophytes  <cit> .
additionally, one single amf species can often be used to inoculate dicotyledons, monocotyledons and ferns, and one plant species can be mycorrhized by several amf species.
therefore, amf are considered not to be host specific and there is no evidence for evolution of host specificity  <cit> .
this finding is extremely remarkable considering the obligate biotrophic life style of amf.
thus, the mechanism of plant infection and colonization seems to be ancient and conserved within amf.
especially during infection, pathogenic and symbiotic plant-microbe interactions show striking similarities, suggesting commonalities in the underlying regulation.
for example, pathogenic and am fungi develop analogous feeding structures, haustoria and arbuscules, respectively  <cit> .
transcriptome profiles of plant cells hosting these structures indicate activated auxin signalling and increased plant metabolism in response to both, pathogenic and symbiotic fungi  <cit> .
the transcriptome sequencing project of the model am fungus, rhizophagus irregularis daom <dig>  showed that only a limited set of cell wall degrading enzymes is expressed during invasive growth, presumably aimed to avoid a major release of polysaccharide fragments, thereby, their detection by the plant immune system  <cit> .
this is analogous to the gene expression patterns in obligate biotrophic pathogens such as the fungus blumeria graminis and the oomycete albugo laibachii  <cit> .
recent studies of plant-mutualistic ectomycorrhizal fungi interactions point in the same direction: genes encoding effector proteins were shown to play a key role in host colonization by controlling the plant immune system  <cit> .
an important factor in plant-microbe interactions are microbial effector proteins released to alter plant cell structure or function, allowing successful infection by suppressing the host defence response  <cit> .
in plant pathogenic fungi and oomycetes two classes of effectors are known.
the first class comprises apoplastic effectors that are secreted into the extracellular space.
they exert their function by interfering with the plant defence in the apoplast and contain small cysteine-rich proteins  and inhibitors of plant extracellular hydrolases including proteases   <cit> , glucanases and plant cell wall degrading enzymes  <cit> .
the second class encompasses intracellular effectors that are translocated into the host cell and target different subcellular compartments.
these effectors may be recognized inside the plant cell by intracellular nb-lrr  proteins encoded by r resistance genes, which results in the induction of effector-triggered immunity   <cit> .
most of the identified eukaryotic pathogenic effectors do not contain domains or homologies to proteins with known function; therefore their roles remain unclear.
although in oomycete pathogens the rxlr motif has been implicated in translocation of effectors into the host cell, to date no universal host-transportation motif has been identified for fungi  <cit> .
fungal and oomycete effector genes seem to be at the forefront of the arms race between host and pathogen: their high rate of nonsynonymous sequence substitutions is a strong signature of positive selection.
this leads to a high amino acid polymorphism, so that, from the  <dig> predicted rxlr genes in the oomycete phytophthora infestans, only  <dig> belong to the “core orthologs” present in other phytophthora species  <cit> .
in contrast to pathogenic systems, there is a vast gap of knowledge on effector proteins in symbiotic fungal-plant interactions  <cit> .
there are indications that also here they may play a role in modulating the plant immune response and consequently the efficiency of fungal infection and colonization of plant cells  <cit> .
kloppholz and colleagues  <cit>  discovered in the am fungus r. irregularis the secreted effector protein, sp <dig>  influencing fungal development within plant roots.
sp <dig> interacts with the pathogenesis-related transcription factor erf <dig> in the medicago truncatula nucleus, and its constitutive expression in roots led to increased mycorrhization and reduced expression of plant defence genes.
since infection and colonization processes appear to be conserved within amf, we would also expect effectors playing a role during am symbiosis to be conserved among amf species.
the level of conservation could then be indicative of whether they play a major role or have instead a rather supplemental function during the symbiosis establishment.
therefore, this scenario would be the exact opposite situation as in pathogenic plant-microbe interactions.
in this study, we predicted and compared the effector candidate repertoire of the two amf, r. irregularis and rhizophagus clarus.
to this end, we have used the published r. irregularis genome  <cit>  and a r. clarus genome draft generated by illumina high-throughput sequencing.
for most of the predicted effector candidates, we were able to identify homologs in both rhizophagus species.
in addition, for several of them we could postulate potential roles in symbiosis by means of gene ontology.
methods
r. clarus monoxenic culture
r. clarus in vitro culture was established as described by declerck et al.
<cit> .
chicory ri-t dna transformed root organ culture  was cultivated with a single-spore isolate of r. clarus mucl <dig> on one side of two-compartment petri dishes containing mrs medium with  <dig>  % phytagel .
fungal hyphae, but not chicory roots, were allowed to grow over the plastic barrier to the second compartment.
plates were incubated in the dark at 26 °c for up to six months until the development of secondary spores.
extraradical mycelia and secondary spores were recovered from the root-free plate compartments by dissolving the medium in 10 mm citrate buffer  <cit> .
the fungal material was collected under sterile conditions  <cit>  and stored at −20 °c or −80 °c until used.
r. clarus genomic dna sequencing
genomic dna was extracted from 100 mg spores and hyphae using the dneasy plant mini kit  following manufacturer’s instructions.
the genomic sequencing library was constructed with the nextera dna sample preparation kit , which allows for low gdna input  and does not need prior dna shearing.
the library was quality controlled by analysis on an agilent  <dig> bioanalyzer with agilent high sensitivity dna kit  for fragment sizes of around 200–400 bp.
sequencing on a miseq sequencer   was performed in the genomics service unit , yielding about  <dig> mio reads and  <dig> gb of primary sequence.
an initial draft with approx.
53-fold coverage of the r. clarus genome was assembled using clc genomics workbench v <dig>   with the following parameters: word size, 24; bubble size, 50; mismatch cost, 2; insertion cost, 3; deletion cost, 3; length fraction,  <dig> ; similarity fraction,  <dig> ; mapping mode, map reads back to contigs; minimum contig length,  <dig> ; mismatch cost,  <dig>  contaminating sequences were identified by analysis with the metagenomics server, mg-rast  <cit> .
80 % of the contigs passed quality control  and from these, around 25 % were recognized as bacterial contamination based on sequences available in databases.
contaminating sequences were removed by mapping against genomes of bacillus cereus , bradyrhizobium elkanii , bradyrhizobium japonicum , legionella pneumophila , and sorangium cellulosum , which were closely related to the main bacterial contaminants.
in addition, all contigs with a gc content > 50 % were eliminated.
84 % of the raw sequence information could be remapped onto the final assembly.
the draft assembly statistics are summarized in additional file 1: table s <dig>  completeness of the draft genome assembly was assessed using the busco software tool  <cit>  for genomes with the fungal gene set and with default parameters .
bioinformatic prediction of the r. irregularis effectome
the r. irregularis daom <dig> secretome was predicted from the initial putative proteome  available at the doe joint genome institute  database by use of a series of domain and protein structure prediction programs.
the  <dig>  proteins carrying n-terminal signal peptides  were predicted using the well validated signalp  <dig>  and predisi programs and sub-cellular targeting using targetp  <cit>  with default parameters.
to ensure stringent standards, only proteins predicted being secreted by at least two methods were considered further.
this prediction was followed by exclusion of proteins likely to be retained into the plasma membrane in which a transmembrane domain  was forecasted by tmhmm  <dig> c and phobius  <cit>  with default parameters.
proteins in which predicted tmds overlapped with signal peptides were not excluded.
a ps-scan analysis with the prosite motif ‘ps00014’  filtered proteins that are likely to permanently reside in the lumen of endoplasmic reticulum  <cit> .
the final number of proteins belonging to the r. irregularis secretome was  <dig> .fig.
<dig>
in silico analysis pipeline for identification of the putative effector repertoire from the genomes of two amf, r. irregularis and r. clarus
effector candidates were identified from the r. irregularis secreted proteins that fulfil at least one of the previously mentioned criteria: proteins with nuclear localization signal .
to find candidate effectors that target the plant nucleus, the r. irregularis secretome was searched for proteins containing a nls using wolf psort, nucleo and eslpred <dig>  software with default parameters .
we identified  <dig> predicted proteins containing at least one nls .
small and cysteine rich  proteins.
several known effectors from fungal and oomycetal plant pathogen are characterized by a length of less than  <dig> amino acids and a cysteine content of more than 3 %  <cit> .
in order to identify candidates achieving this criterion the cysteine content and length of r. irregularis secreted proteins were calculated using clc genomics workbench.
within the secretome, we found a total of  <dig> proteins fulfilling this condition .
repeats containing proteins .
internal repeats are often found in oomycetal and fungal effectors  <cit> .
using the t-reks algorithm and xstream software with default parameters  <cit>  tandem repeat proteins were selected within the r. irregularis secretome.
sixty secreted proteins showed internal repeats .
in silico search for effector candidate homologs in r. clarus
proteins of r. irregularis matching at least one effector criterium were subsequently used to predict candidate effectors in r. clarus .
approx.
95 %  of r. irregularis effector candidates showed blastp  detectable hits with minimum similarity of 40 % within the r. clarus initial genome draft.
even though these proteins share sequence similarity they do not necessarily have a similar topology.
therefore, we applied the same selection criteria that were used on the r. irregularis secretome, filtering out proteins with tmds and endoplasmic reticulum retention signals.
also, an additional control for the presence of sp, nls and internal repeats was performed.
effector protein characterization by blast2go
to characterized predicted effector proteins the blast2go software was used  <cit> .
at first, protein sequences encoded by effector candidate genes were compared with the non-redundant sequence database.
analysis were performed using the blastp algorithm with expect-value lower than  <dig> e− <dig> recording max.
<dig> hits.
proteins with significant hits were classified into gene ontology categories by ‘mapping’.
data from both analyses were merged into the annotation.
additionally, presence of conserved protein domains or motifs was determined using pfam, superfamily  <dig>  software and ncbi conserved domain search with default parameters  <cit> .
RESULTS
defining the rhizophagus effector repertoire
to identify and classify effector candidates of the rhizophagus genus we constructed a bioinformatics pipeline using current knowledge of the properties of validated fungal and oomycete plant pathogen effectors .
it was observed that known pathogen effector proteins from fungi and oomycetes are secreted and fulfil at least one of the following criteria:  they contain a nuclear localization signal ,  they are small and cysteine rich ,  they contain internal repeats , or  they show similarity to haustorial  expressed proteins  <cit> .
for r. irregularis so far no inventory of arbuscular expressed proteins exist.
therefore, we restricted the effector candidate identification to the first three criteria.
effector candidates need to carry an n-terminal signal peptide  for secretion and should not be retained in the membrane or endoplasmic reticulum of the secreting organism.
the r. irregularis daom <dig> secretome was predicted from the initial putative proteome  available at the doe joint genome institute  database by use of a series of domain and protein structure prediction programs.
the predicted r. irregularis secretome contains  <dig> proteins .
the final number of r. irregularis predicted secreted proteins that fulfil at least one of the above-mentioned effector criteria was  <dig>  representing 23 % of all r. irregularis putatively secreted proteins .
in addition, we found  <dig> proteins fulfilling two but none that matched all three of the effector criteria.
brief gene ontology category analyses by protfun  <dig>  server  <cit>  of both, secretome and effectome, showed enrichment in hormone, stress response, growth factor, and immune response related proteins .
when comparing secretome and effectome, only three differences are apparent.
the 40 % increased number of hormone-related proteins observed in the effectome correlates with the finding that in am the root architecture is modified by essential, hormone-mediated pathways that might be controlled by amf  <cit> .
a reduction by 45 % of proteins involved in immune response between effectome and secretome might be indicative of the symbiotic interaction between amf and plant.
the same could be true for the relative reduction of proposed receptor proteins by 60 % in the effectome.
receptors play an important role in all kind of interactions by detecting and recognising signals from pathogenic or symbiotic organisms present in the environment  <cit> .
the effectome is conserved
given the central importance of secreted proteins in determining the outcome of plant-microbe encounters in other systems, it is likely that amf effectors also have a crucial function.
however, knowledge about amf effectors and their function during symbiosis is limited and to date only one effector has been functionally characterised  <cit> .
we wanted to address two questions:  how conserved is the effectome between two closely related amf species?
and , does the identity of the conserved putative effectors shed light onto their potential function during the amf colonization process?
for this purpose we generated a draft genome of r. clarus mucl <dig> by high-throughput sequencing resulting in  <dig>  scaffolds with a total of 94 mb, in size similar to the 101 mb gloin <dig> assembly of r. irregularis daom <dig>  <cit> .
an analysis of genome completeness using the busco software tool  <cit>  showed comparable levels for both genomes .
therefore, we estimate the complete r. clarus mucl <dig> genome to be slightly smaller than the 150 mb total genome size predicted for r. irregularis daom <dig>  <cit> .
to determine the level of conservation within the amf symbiotic effectome, we compared protein sequences encoded by selected housekeeping genes and proteins involved in symbiosis from r. irregularis and r. clarus and calculated the level of identity and similarity by blastp   .
the identity of selected proteins varied between 84 % for ammonium transporter  <dig> to 99 % for the conserved elongation factor 1-alpha, whereas the similarity varied between  <dig> and 99 %, respectively.
the high levels of protein identity and similarity confirm the close relationship between these two species of the same genus, which is in agreement with current knowledge  <cit> .
subsequently, we searched within the r. clarus genome for homologs of the r. irregularis proteins matching at least one effector criterion.
approx.
95 %  of r. irregularis effector candidates showed blastp  detectable hits within the r. clarus genome draft .
the surprisingly high number of putative effector protein homologs present in both rhizophagus species clearly speaks in favour of a conserved role during the infection process and demonstrates the potential of a comparative approach.
although the r. clarus draft genome is not complete and therefore incomplete orf information in the homolog sets were to be expected, analysis of the overlapping effectome could still reveal potential roles in symbiosis.
for 30 % of the effector candidates predicted in r. irregularis we found only truncated homolog sequences in r. clarus.
this can only be partly explained by genome incompleteness, because the two genome assemblies were of comparable quality .
a further 25 % of candidates in r. clarus had to be omitted because the inaccurate annotation of the r. irregularis proteome resulted in selection of false positive candidates in the effectome, as properly annotated r. clarus homologs of these candidates did not fulfil effectome criteria.
additionally, even though r. clarus effector protein homologs share sequence similarity with that of r. irregularis, they do not necessarily have a similar topology.
therefore, after applying the same selection criteria used on the r. irregularis secretome, where proteins predicted to hold sp, nls, rcp but lacking tmd and er retention signal were kept, the r. clarus candidate number was reduced by 13 %.
the final number of full-length sequence candidates fulfilling effector selection criteria present in r. clarus as well as in r. irregularis genome was  <dig>  which is approx.
32 % of all candidate homologs present in both species .
the protein similarity between the homologs was 40–98 %.
within the  <dig> r. clarus effector candidates we identified only five proteins that fulfilled two of the effector criteria and none that matched all three of them .
furthermore, the similarity between all effector candidates from both rhizophagus species was ranged from 51 % for nls/scr_ <dig>  to 97 % for nls_ <dig> as determined by the blastp  algorithm .
from these,  <dig> conserved effector candidates show more than 90 % similarity between r. clarus and r. irregularis, which falls within the same range as the similarity between housekeeping genes .
more than 44 %  of the conserved candidates belong to the nls group .
in addition, the nls group shows the highest conservation level  in comparison to the scr and rcp groups .
these findings suggest that the effectors translocated to the plant nucleus with possible effect on the transcriptional regulation in plant cells might be essential for am symbiosis.
the highly conserved effector candidates identified here between different amf species are probably playing a fundamental role during fungus-plant interaction also explaining the degree of unspecific interaction with all potential hosts.
in contrast, effector candidates showing high diversification may have minor host specific functions.table  <dig> the r. clarus effector candidates with predicted function
potential functional features of r. clarus effector candidates
to deduct putative functions of r. clarus and r. irregularis predicted effector candidates the blast2go software with automated annotation was used  <cit> .
the obtained information was used to predict putative cellular functions of rhizophagus effector candidates .
only  <dig> of r. clarus effector candidates  showed sequence similarity  to documented proteins with defined or predicted function, and  <dig>  contain known domains or motifs .
predicted domain and homology hits with an e-value higher than 10− <dig> were ignored  <cit> .
current research showed that a number of fungal and oomycete effectors despite being under high diversifying or positive selective pressure still carry recognizable pfam domains, which might be helpful in prediction of their biological function  <cit> .
thus, for  <dig>  of the r. clarus effector candidates the putative function could be determined, hence their potential role in am symbiosis was defined and subdivided into three main functional groups: signalling, cell wall modification, and transcriptional regulation.
however, the majority of r. clarus effector candidates are novel sequences showing no significant homology to known proteins in other organisms or sequence motifs , which is in agreement with previous studies highlighting the evolutionary diverse nature of fungal effector groups  <cit> .
based on computational approaches, we established effector candidate subcategories and speculated on putative functions after careful confirmation of presence of active sites and completeness of the predicted domains.fig.
<dig> illustration of the various functional modules of r. clarus effector candidates
secreted proteins likely to target plant signalling
secreted effectors can alter plant defence by interfering with signal transduction  <cit> .
within the r. clarus effectome we detected eight effector candidates that might affect plant signalling.
within the r. clarus effectome, candidate scr_ <dig> shows homology to the thermoanaerobacter sp.
peptidoglycan-binding lysm protein and contains a lysm domain .
plant receptor kinases containing a lysm  domain are able to detect chitin-like components or chitin fragments  <cit> .
chitin, an n-acetyl-glucosamine  polymer, is a main component of the fungal cell wall  <cit> .
chitin recognition activates signalling pathways, which lead to plant immune defence  <cit> .
due to this, it is suggested that many plant pathogens secrete effectors that contain lysm domains in order to inhibit the release of chitin fragments from fungal cell walls or to bind these fragments, thus, preventing their recognition by plant receptors.
for instance, the rice blast fungus magnaporthe oryzae secretes the lysm domain containing effector, slp <dig>  which binds chitin oligosaccharides.
this suppresses plant defence response induced through the recognition of chitin fragments by the rice receptor cebip   <cit> .
candidate scr_ <dig> might play a similar role in am symbiosis.
the second effector candidate nls_ <dig> contains a f-box-like domain .
proteins containing this domain are components of the e <dig> ubiquitin ligases called scf complexes  which control many intracellular regulatory processes by catalysing the ubiquitination of proteins that targets proteins for degradation  <cit> .
intense analysis showed that plants defence response pathways might be under control of an scf complex  <cit> .
many pathogenic effectors use a strategy of ubiquitination in order to regulate plant immunity.
the avrpiz-t effector of the m. oryzae binds to and destabilizes the rice ring e <dig> ubiquitin ligase apip <dig> to suppress chitin-induced plant immunity  <cit> .
the third effector candidate nls_ <dig> shows homology to the ochotona princeps nitrogen permease regulator 3-like protein  isoform x <dig> and contains a npr <dig> domain  typical of this protein family.
npr <dig> forms a heterodimer complex with npr <dig>  protein, which mediates an amino acid starvation signal to torc <dig>   <cit> .
torc <dig>  an essential serine/threonine protein kinase, incorporates signals from diverse pathways in the cell.
the homology of candidate nls_ <dig> to npr <dig> protein suggests its interference with signal transduction most likely during amino acid deficiency conditions.
a similar role in interfering with plant signalling pathways might have the fourth candidate nls_ <dig> in which a protein kinases catalytic domain, pkc  was detected.
pkc domains can transfer a phosphate group from atp to serine/threonine or tyrosine residues on protein substrates.
pks often function as components of signal transduction pathways in which one kinase activates another kinase.
this serial action transfers a signal within a cell, which results in cellular responses  <cit> .
two additional candidates might also influence the plant signalling cascades by being involved in hydrolytic processes.
the candidate nls_ <dig> shows homology to the dichomitus squalens plc-like  phosphodiesterase and has a glycerophosphodiester phosphodiesterase domain, gdpd .
plc, by hydrolysing phosphodiester bonds, integrates signalling cascades critical for a variety of cellular and physiological functions  <cit> .
the candidate rcp_ <dig> shows homology to the grosmannia clavigera phospholipase d <dig>  which catalyse the cleavage of phospholipids into fatty acids and other lipophilic substances, like phosphatidic acid, which is an established intracellular signalling lipid  <cit> .
candidate number seven, nls_ <dig>  has homology to the glycosyltransferase family  <dig> protein from colletotrichum higginsianum and encompasses a complex bifunctional domain, which includes a trehalose-6-phosphate synthase  domain  and a had-superfamily hydrolase domain, subfamily iib .
tps catalyses the synthesis of α,α- <dig> -trehalose-6-phosphate from glucose-6-phosphate using a udp-glucose donor and is a key enzyme in the trehalose synthesis pathway  <cit> .
furthermore, trehalose was reported to have a role as a stress protectant in many organisms across kingdoms  <cit> .
remarkably, in pathogens trehalose is associated with several mechanisms involved in cell wall integrity, regulation of metabolism, and successful infection  <cit> .
absence of a trehalose metabolic apparatus is correlated with lower pathogen virulence.
however, yet the underlying mechanism remains unravelled.
for instance, in m. oryzae tps <dig> regulates transcriptional effectors linked to virulence factors via nadph-binding during appressorium-mediated rice infection  <cit> .
consequently, the putative effector nls_ <dig> might take part in the colonization step of amf symbiosis by being involved in trehalose synthesis.
the last candidate, scr_ <dig>  contains a ml  domain, which is involved in specific lipid recognition.
it is predicted that ml domain containing proteins through interacting with specific lipids can mediate diverse biological functions, like lipid metabolism, host response to pathogen components like lipopolysaccharides, and other cellular functions  <cit> .
this suggests that the scr_ <dig> candidate, by interacting with lipids, can interfere with the plant lipid signalling; however the direct mechanism of this effect needs to be clarified.
genes encoding putative secreted enzymes engaged in cell wall modification
we identified four r. clarus effector candidates encoding putative secreted enzymes involved in cell wall modification.
they were divided in two subgroups depending whether their putative biological functions concern plant or fungal cell wall.
enzymes involved in degradation of the plant cell wall
most of the necrotrophic fungi degrade plant tissues by producing a broad set of carbohydrate-active enzymes specifically focused on plant polysaccharide degradation  <cit> .
transcriptomic profiling of r. irregularis daom <dig> predicted  <dig> genes encoding carbohydrate-active enzymes  <cit> .
this amount is relatively low in comparison to obligate biotrophic pathogens and ectomycorrhizal symbionts  <cit> .
also, only one putative enzyme involved in plant cell wall modification was found in the r. clarus effector repertoire.
this suggests that during am symbiosis the plant cell wall is not completely degraded, but instead a localized corruption of its structure might occur to make it more flexible, possibly for fungus penetration into the plant cell.
the single effector of this subgroup, candidate nls/rcp_ <dig> shows homology to the rhizoctonia solani galactose oxidase and contains a glyoxal oxidase domain .
glyoxal oxidase has a copper active site remarkably similar to that found in another fungal metalloenzyme, galactose oxidase  <cit> .
glyoxal oxidase is a copper metalloenzyme that catalyses the oxidation of aldehydes to carboxylic acids and is an essential component of the extracellular lignin degradation pathways  <cit> .
although lignin is present only in small amounts in the plant cell wall  <cit> , it has a strong impact on the rigidity of its structure.
lignin degradation might result in a more flexible cell wall necessary for fungal penetration.
enzymes involved in modification of the fungal cell wall
the effectome also contains three putative enzymes involved in chitin degradation.
this suggests constant remodelling of the fungal cell wall, likely an adaptation to substantial differences in environmental conditions in the soil and inside of plant tissues.
similarly, in cryptococcus neoformans it was shown that chitosan, the product of chitin deacetylation, helps to sustain cell integrity by maintaining normal capsule width during vegetative growth  <cit> .
the first effector candidate in r. clarus, rcp_ <dig>  contains a glycosyl hydrolase family  <dig> domain  and shows homology to the neofusicoccu parvum class i alpha-mannosidase enzyme.
this enzyme is involved in the cleavage of the terminal mannose from asn-linked oligosaccharides, which takes part at the early step of hydrolysis of n-glycan such as chitin  <cit> .
two other r. clarus effector candidates might function in fungal cell wall remodelling by chitin deacetylation  <cit> .
candidate rcp_ <dig> shows homology to xylase/chitin deacetylase enzyme from cryptococcus gatti and contains a chitin deacetylase domain .
the other candidate, rcp_ <dig>  shows a significant homology to rhizoctonia solani bifunctional xylanase/deacetylase.
the presence of hevein/chitin binding domain  and a catalytic nodb homology domain of the carbohydrate esterase  <dig> superfamily , which typically occur in chitin deacetylases  <cit> , suggest a chitin degradation function of both candidates.
however, candidates with putative chitin deacetylase function might also have another purpose.
it was shown that chitin deacetylase plays a significant role in protecting pathogenic fungal hyphae from being lysed by plant secreted chitinases, specifically in the cases of the wheat stem rust fungus puccinia graminis f. sp.
tritici and the broad bean rust fungus uromyces fabae  <cit> .
furthermore, it was suggested that exposed fungal chitin polymers are partially de-n-acetylated by chitin deacetylases during the infection of the host, this with purpose to avoid plant antimicrobial hydrolases by affinity modulation  <cit> .
genes encoding proteins with putative transcriptional regulatory function
during the pathogenic or symbiotic infection thousands of genes are up- or down-regulated in the plant as revealed by several transcriptome analyses.
expression of genes encoding enzymes involved in protein degradation, cell wall modification, and secondary metabolite biosynthesis is induced in order to activate plant immune response  <cit> .
however, pathogens are able to secrete the so-called transcription activation-like effectors, which target dna and are involved in transcriptional regulation of a set of plant host genes, for instance to prevent defence response  <cit> .
we found only one putative candidate within the r. clarus effectome that might regulate plant gene expression during establishment of the am symbiosis.
this is the nuclear localized  candidate, nls_ <dig>  which consists of a pwwp domain , named after a conserved pro-trp-trp-pro motif.
pwwp domain-containing proteins bind to chromatin and control its packing resulting in transcription and replication regulation.
moreover, pwwp proteins are also involved in dna damage response.
while dna damage occurs, they regulate chromatin condensation and ability of dna repair factors, like 53bpi and ctip to bind dna  <cit> .
similarly, candidate nls_ <dig> might be involved in transcription regulation by chromatin remodelling or condensation.
effector candidates with unknown function
six of the r. clarus effector candidates show homologies to known proteins or contain identified domains or motifs.
however, the similarity does not allow for prediction of a specific function.
one of the candidates, nls_ <dig>  has an aaa  domain and shows significant homology to a protozoa protein containing a denn  domain .
however, the homology is mostly based on the aaa domain region as a denn domain was not detected.
the aaa superfamily is a large and functionally diverse superfamily composed of members that are essential for a variety of cellular functions, like protein unfolding and degradation, dna recombination, replication and repair, and signal transduction  <cit> .
candidates scr_ <dig> and rcp_ <dig> show homologies to putative conserved fungal proteins with unknown function.
both of them lack of any functional domain; therefore, their putative function cannot be predicted.
candidate nls_ <dig> shows homology to the trametes versicolor p-loop protein containing a nucleoside triphosphate hydrolase domain .
nevertheless, this homology is partial and does not include the characteristics for this family of p-loop ntpase domain proteins.
in addition, two other candidates, scr_ <dig> and rcp_ <dig>  show homology to the laccari bicolor carbohydrate binding module family  <dig> protein .
proteins belonging to this family contain a module of 60– <dig> residues  with a chitin binding function  <cit> .
however, once more, this module is absent in both candidates; hence, despite the homology the putative function cannot be predicted.
CONCLUSIONS
symbiotic interaction between amf and plant host requires a well-adjusted dialog between the two partners in order to function properly.
we hypothesise that one of the key components of this communication is the fungal effectome.
in this study, we used a bioinformatic pipeline to predict effector candidates from two amf belonging to the rhizophagus genus, r. irregularis and r. clarus.
our in silico pipeline revealed a list of several candidate effector proteins that create a valuable source of information to elucidate the mechanism of plant infection and colonization by am fungi.
our in silico characterization of the amf effectome focused on r. clarus candidates.
the majority of the r. clarus effectors are novel proteins with no significant homology to known sequences or motifs in other organisms.
this finding is in agreement to previous reports enhancing the evolution diverse nature of fungal and oomycete effectors  <cit> .
however, for a small number of candidates putative functions during plant-fungus symbiosis might be predicted.
the set of putative effector candidates encoded by the r. clarus effectome appears to be suited for a fungus that wants to establish a symbiotic rather than a pathogenic relation with the host.
the presence of only one candidate encoding a putative plant cell wall degrading enzyme  in the r. clarus effectome is in agreement with earlier reports, in which it was shown that several obligate biotrophic pathogens and ectomycorrhizal symbionts have a decreased repertoire of pcwdes compared to saprophytes  <cit> .
in comparison to pathogenic fungi, amf seem not to degrade the plant cells but rather grow between them or inside of them while trying to avoid recognition by the plant defence.
however, localized plant cell wall modelling probably does occur to weaken and make it more flexible during fungal penetration and passage from cell to cell.
the lack of the typical pathogenic effectors, such as extracellular lipases, proteases, nucleases and phytases, additionally confirms the non-pathogenic status of amf  <cit> .
the relationship between symbiotic fungi and plants requires a major rearrangement of the plant signalling pathways for which fungi have evolved a range of secreted proteins.
for instance, fungal effectors control auxin signalling within plant cells to alter plant root architecture and favour colonization  <cit> .
thus, the presence of candidates in the r. clarus effectome that might target plant signalling in order to allow fungal entry was expected.
reprogramming of the plant defence could be used to avoid recognition.
alternatively, effectors might affect transport processes across the plant cell membrane.
this was previously shown for symbiotic fungi which encode effectors to control the transfer of glucose into the apoplastic space of the root where it is easily accessible for fungus  <cit> .
a remarkable aspect of our work is presence of a large portion of the effectors that are not species-specific, but conserved between different rhizophagus species.
this group of highly conserved effector proteins may play a fundamental role during fungus-plant interaction and has first priority for future functional characterization.
availability of supporting data
the sequences of the predicted effector candidates have been deposited in genbank under accession numbers: ku <dig> - ku <dig> .
other supporting data are included in an additional file  <dig>
additional file
additional file 1: figure s <dig>  gene ontology category analyses by protfun  <dig>  server  <cit>  of r. irregularis daom <dig> secretome and effectome.
table s <dig>
r. clarus draft genome assembly statistics in comparison to published r. irregularis genome assemblies.
table s <dig>  the r. irregularis putative effectome.
table s <dig>  identity and similarity levels between r. irregularis and r. clarus based on sequences of selected proteins .
table s <dig>  the r. clarus effector candidates.
abbreviations
aaaatpases associated with diverse cellular activities
amarbuscular mycorrhizal
amfarbuscular mycorrhizal fungi
cebipchitin elicitor binding protein
denndifferentially expressed in neoplastic vs. normal cells
erendoplasmic reticulum
etieffector-triggered immunity
gdpdglycerophosphodiester phosphodiesterase domain
glcnacn-acetyl-glucosamine
lysmlysine motif
m. oryzaemagnaporthe oryzae
miomillion
mlmd-2-related lipid-recognition
myamillion years ago
nb-lrrnucleotide-binding site leucine-rich repeat
nlsnuclear localization signal
npr2nitrogen permease regulator 2
npr3nitrogen permease regulator 3-like protein
pcwdeplant cell wall degrading enzyme
pkcprotein kinases catalytic domain
plc-likephospholipase c-like
p-loop ntpasep-loop protein with a nucleoside triphosphate hydrolase domain
pwwppro-trp-trp-pro
qcquality control
r. clarusrhizophagus clarus
r. irregularisrhizophagus irregularis
rcprepeats containing proteins
scrsmall cysteine-rich
spsignal peptide
sp7secreted effector protein 7
tmdtransmembrane domain
torc1target of rapamycin kinase
tpstrehalose-6-phosphate synthase
competing interests
the authors declare that they have no competing interests.
authors’ contributions
conceived and designed the experiments: ab, kst.
genome sequencing and assembly: ab.
bioinformatics and data analysis: kst.
wrote manuscript: kst, ab.
both authors read and approved the final manuscript.

 arbuscular mycorrhizal fungi live in an obligate symbiosis with the roots of around two-thirds of all studied land plant species, i.e. angiosperms, gymnosperms, pteridophytes and some bryophytes, and in one known case with cyanobacteria. this interaction has a major impact on the entire soil ecosystem and plays a crucial role in agricultural systems by increasing plant tolerance to biotic and abiotic stresses. in the field, amf form mycelial networks providing extensive connections between roots of different plant species and mediating interactions among plants. these networks significantly improve the rhizosphere soil stability. the major benefit for plants upon mycorrhization is a more efficient nutrient uptake from the soil. in exchange, amf obtain photoassimilates from their hosts. all amf belong to the phylum glomeromycota which has been divided into four orders: the glomerales, the diversisporales, and the two ‘ancestral’ lineages, paraglomerales and archaeosporales. they may have started to diverge over mya. the earliest spore fossils have been found together with the first land plants and resemble present amf morphological structures. it has been suggested that amf played a crucial role for the adaptation of phototrophs to the terrestrial environment. the am symbiosis persisted morphologically unchanged throughout the complete evolutionary development within the plant phylum from haploid gametophytes to diploid sporophytes. additionally, one single amf species can often be used to inoculate dicotyledons, monocotyledons and ferns, and one plant species can be mycorrhized by several amf species. therefore, amf are considered not to be host specific and there is no evidence for evolution of host specificity. this finding is extremely remarkable considering the obligate biotrophic life style of amf. thus, the mechanism of plant infection and colonization seems to be ancient and conserved within amf. especially during infection, pathogenic and symbiotic plant-microbe interactions show striking similarities, suggesting commonalities in the underlying regulation. for example, pathogenic and am fungi develop analogous feeding structures, haustoria and arbuscules, respectively. transcriptome profiles of plant cells hosting these structures indicate activated auxin signalling and increased plant metabolism in response to both, pathogenic and symbiotic fungi. the transcriptome sequencing project of the model am fungus, rhizophagus irregularis daom showed that only a limited set of cell wall degrading enzymes is expressed during invasive growth, presumably aimed to avoid a major release of polysaccharide fragments, thereby, their detection by the plant immune system. this is analogous to the gene expression patterns in obligate biotrophic pathogens such as the fungus blumeria graminis and the oomycete albugo laibachii. recent studies of plant-mutualistic ectomycorrhizal fungi interactions point in the same direction: genes encoding effector proteins were shown to play a key role in host colonization by controlling the plant immune system. an important factor in plant-microbe interactions are microbial effector proteins released to alter plant cell structure or function, allowing successful infection by suppressing the host defence response. in plant pathogenic fungi and oomycetes two classes of effectors are known. the first class comprises apoplastic effectors that are secreted into the extracellular space. they exert their function by interfering with the plant defence in the apoplast and contain small cysteine-rich proteins and inhibitors of plant extracellular hydrolases including proteases, glucanases and plant cell wall degrading enzymes. the second class encompasses intracellular effectors that are translocated into the host cell and target different subcellular compartments. these effectors may be recognized inside the plant cell by intracellular nb-lrr proteins encoded by r resistance genes, which results in the induction of effector-triggered immunity. most of the identified eukaryotic pathogenic effectors do not contain domains or homologies to proteins with known function; therefore their roles remain unclear. although in oomycete pathogens the rxlr motif has been implicated in translocation of effectors into the host cell, to date no universal host-transportation motif has been identified for fungi. fungal and oomycete effector genes seem to be at the forefront of the arms race between host and pathogen: their high rate of nonsynonymous sequence substitutions is a strong signature of positive selection. this leads to a high amino acid polymorphism, so that, from the predicted rxlr genes in the oomycete phytophthora infestans, only belong to the “core orthologs” present in other phytophthora species. in contrast to pathogenic systems, there is a vast gap of knowledge on effector proteins in symbiotic fungal-plant interactions. there are indications that also here they may play a role in modulating the plant immune response and consequently the efficiency of fungal infection and colonization of plant cells. kloppholz and colleagues discovered in the am fungus r. irregularis the secreted effector protein, sp influencing fungal development within plant roots. sp interacts with the pathogenesis-related transcription factor erf in the medicago truncatula nucleus, and its constitutive expression in roots led to increased mycorrhization and reduced expression of plant defence genes. since infection and colonization processes appear to be conserved within amf, we would also expect effectors playing a role during am symbiosis to be conserved among amf species. the level of conservation could then be indicative of whether they play a major role or have instead a rather supplemental function during the symbiosis establishment. therefore, this scenario would be the exact opposite situation as in pathogenic plant-microbe interactions. in this study, we predicted and compared the effector candidate repertoire of the two amf, r. irregularis and rhizophagus clarus. to this end, we have used the published r. irregularis genome and a r. clarus genome draft generated by illumina high-throughput sequencing. for most of the predicted effector candidates, we were able to identify homologs in both rhizophagus species. in addition, for several of them we could postulate potential roles in symbiosis by means of gene ontology. methods r. clarus monoxenic culture r. clarus in vitro culture was established as described by declerck et al.. chicory ri-t dna transformed root organ culture was cultivated with a single-spore isolate of r. clarus mucl on one side of two-compartment petri dishes containing mrs medium with  % phytagel. fungal hyphae, but not chicory roots, were allowed to grow over the plastic barrier to the second compartment. plates were incubated in the dark at 26 °c for up to six months until the development of secondary spores. extraradical mycelia and secondary spores were recovered from the root-free plate compartments by dissolving the medium in 10 mm citrate buffer. the fungal material was collected under sterile conditions and stored at −20 °c or −80 °c until used. r. clarus genomic dna sequencing genomic dna was extracted from 100 mg spores and hyphae using the dneasy plant mini kit following manufacturer’s instructions. the genomic sequencing library was constructed with the nextera dna sample preparation kit, which allows for low gdna input and does not need prior dna shearing. the library was quality controlled by analysis on an agilent bioanalyzer with agilent high sensitivity dna kit for fragment sizes of around 200–400 bp. sequencing on a miseq sequencer was performed in the genomics service unit, yielding about mio reads and gb of primary sequence. an initial draft with approx. 53-fold coverage of the r. clarus genome was assembled using clc genomics workbench v with the following parameters: word size, 24; bubble size, 50; mismatch cost, 2; insertion cost, 3; deletion cost, 3; length fraction,; similarity fraction,; mapping mode, map reads back to contigs; minimum contig length,; mismatch cost, contaminating sequences were identified by analysis with the metagenomics server, mg-rast. 80 % of the contigs passed quality control and from these, around 25 % were recognized as bacterial contamination based on sequences available in databases. contaminating sequences were removed by mapping against genomes of bacillus cereus, bradyrhizobium elkanii, bradyrhizobium japonicum, legionella pneumophila, and sorangium cellulosum, which were closely related to the main bacterial contaminants. in addition, all contigs with a gc content > 50 % were eliminated. 84 % of the raw sequence information could be remapped onto the final assembly. the draft assembly statistics are summarized in additional file 1: table s completeness of the draft genome assembly was assessed using the busco software tool for genomes with the fungal gene set and with default parameters. bioinformatic prediction of the r. irregularis effectome the r. irregularis daom secretome was predicted from the initial putative proteome available at the doe joint genome institute database by use of a series of domain and protein structure prediction programs. the proteins carrying n-terminal signal peptides were predicted using the well validated signalp and predisi programs and sub-cellular targeting using targetp with default parameters. to ensure stringent standards, only proteins predicted being secreted by at least two methods were considered further. this prediction was followed by exclusion of proteins likely to be retained into the plasma membrane in which a transmembrane domain was forecasted by tmhmm c and phobius with default parameters. proteins in which predicted tmds overlapped with signal peptides were not excluded. a ps-scan analysis with the prosite motif ‘ps00014’ filtered proteins that are likely to permanently reside in the lumen of endoplasmic reticulum. the final number of proteins belonging to the r. irregularis secretome was.fig. in silico analysis pipeline for identification of the putative effector repertoire from the genomes of two amf, r. irregularis and r. clarus effector candidates were identified from the r. irregularis secreted proteins that fulfil at least one of the previously mentioned criteria: proteins with nuclear localization signal. to find candidate effectors that target the plant nucleus, the r. irregularis secretome was searched for proteins containing a nls using wolf psort, nucleo and eslpred software with default parameters. we identified predicted proteins containing at least one nls. small and cysteine rich proteins. several known effectors from fungal and oomycetal plant pathogen are characterized by a length of less than amino acids and a cysteine content of more than 3 %. in order to identify candidates achieving this criterion the cysteine content and length of r. irregularis secreted proteins were calculated using clc genomics workbench. within the secretome, we found a total of proteins fulfilling this condition. repeats containing proteins. internal repeats are often found in oomycetal and fungal effectors. using the t-reks algorithm and xstream software with default parameters tandem repeat proteins were selected within the r. irregularis secretome. sixty secreted proteins showed internal repeats. in silico search for effector candidate homologs in r. clarus proteins of r. irregularis matching at least one effector criterium were subsequently used to predict candidate effectors in r. clarus. approx. 95 % of r. irregularis effector candidates showed blastp detectable hits with minimum similarity of 40 % within the r. clarus initial genome draft. even though these proteins share sequence similarity they do not necessarily have a similar topology. therefore, we applied the same selection criteria that were used on the r. irregularis secretome, filtering out proteins with tmds and endoplasmic reticulum retention signals. also, an additional control for the presence of sp, nls and internal repeats was performed. effector protein characterization by blast2go to characterized predicted effector proteins the blast2go software was used. at first, protein sequences encoded by effector candidate genes were compared with the non-redundant sequence database. analysis were performed using the blastp algorithm with expect-value lower than e− recording max. hits. proteins with significant hits were classified into gene ontology categories by ‘mapping’. data from both analyses were merged into the annotation. additionally, presence of conserved protein domains or motifs was determined using pfam, superfamily software and ncbi conserved domain search with default parameters.  defining the rhizophagus effector repertoire to identify and classify effector candidates of the rhizophagus genus we constructed a bioinformatics pipeline using current knowledge of the properties of validated fungal and oomycete plant pathogen effectors. it was observed that known pathogen effector proteins from fungi and oomycetes are secreted and fulfil at least one of the following criteria: they contain a nuclear localization signal, they are small and cysteine rich, they contain internal repeats, or they show similarity to haustorial expressed proteins. for r. irregularis so far no inventory of arbuscular expressed proteins exist. therefore, we restricted the effector candidate identification to the first three criteria. effector candidates need to carry an n-terminal signal peptide for secretion and should not be retained in the membrane or endoplasmic reticulum of the secreting organism. the r. irregularis daom secretome was predicted from the initial putative proteome available at the doe joint genome institute database by use of a series of domain and protein structure prediction programs. the predicted r. irregularis secretome contains proteins. the final number of r. irregularis predicted secreted proteins that fulfil at least one of the above-mentioned effector criteria was representing 23 % of all r. irregularis putatively secreted proteins. in addition, we found proteins fulfilling two but none that matched all three of the effector criteria. brief gene ontology category analyses by protfun server of both, secretome and effectome, showed enrichment in hormone, stress response, growth factor, and immune response related proteins. when comparing secretome and effectome, only three differences are apparent. the 40 % increased number of hormone-related proteins observed in the effectome correlates with the finding that in am the root architecture is modified by essential, hormone-mediated pathways that might be controlled by amf. a reduction by 45 % of proteins involved in immune response between effectome and secretome might be indicative of the symbiotic interaction between amf and plant. the same could be true for the relative reduction of proposed receptor proteins by 60 % in the effectome. receptors play an important role in all kind of interactions by detecting and recognising signals from pathogenic or symbiotic organisms present in the environment. the effectome is conserved given the central importance of secreted proteins in determining the outcome of plant-microbe encounters in other systems, it is likely that amf effectors also have a crucial function. however, knowledge about amf effectors and their function during symbiosis is limited and to date only one effector has been functionally characterised. we wanted to address two questions: how conserved is the effectome between two closely related amf species? and, does the identity of the conserved putative effectors shed light onto their potential function during the amf colonization process? for this purpose we generated a draft genome of r. clarus mucl by high-throughput sequencing resulting in scaffolds with a total of 94 mb, in size similar to the 101 mb gloin assembly of r. irregularis daom. an analysis of genome completeness using the busco software tool showed comparable levels for both genomes. therefore, we estimate the complete r. clarus mucl genome to be slightly smaller than the 150 mb total genome size predicted for r. irregularis daom. to determine the level of conservation within the amf symbiotic effectome, we compared protein sequences encoded by selected housekeeping genes and proteins involved in symbiosis from r. irregularis and r. clarus and calculated the level of identity and similarity by blastp. the identity of selected proteins varied between 84 % for ammonium transporter to 99 % for the conserved elongation factor 1-alpha, whereas the similarity varied between and 99 %, respectively. the high levels of protein identity and similarity confirm the close relationship between these two species of the same genus, which is in agreement with current knowledge. subsequently, we searched within the r. clarus genome for homologs of the r. irregularis proteins matching at least one effector criterion. approx. 95 % of r. irregularis effector candidates showed blastp detectable hits within the r. clarus genome draft. the surprisingly high number of putative effector protein homologs present in both rhizophagus species clearly speaks in favour of a conserved role during the infection process and demonstrates the potential of a comparative approach. although the r. clarus draft genome is not complete and therefore incomplete orf information in the homolog sets were to be expected, analysis of the overlapping effectome could still reveal potential roles in symbiosis. for 30 % of the effector candidates predicted in r. irregularis we found only truncated homolog sequences in r. clarus. this can only be partly explained by genome incompleteness, because the two genome assemblies were of comparable quality. a further 25 % of candidates in r. clarus had to be omitted because the inaccurate annotation of the r. irregularis proteome resulted in selection of false positive candidates in the effectome, as properly annotated r. clarus homologs of these candidates did not fulfil effectome criteria. additionally, even though r. clarus effector protein homologs share sequence similarity with that of r. irregularis, they do not necessarily have a similar topology. therefore, after applying the same selection criteria used on the r. irregularis secretome, where proteins predicted to hold sp, nls, rcp but lacking tmd and er retention signal were kept, the r. clarus candidate number was reduced by 13 %. the final number of full-length sequence candidates fulfilling effector selection criteria present in r. clarus as well as in r. irregularis genome was which is approx. 32 % of all candidate homologs present in both species. the protein similarity between the homologs was 40–98 %. within the r. clarus effector candidates we identified only five proteins that fulfilled two of the effector criteria and none that matched all three of them. furthermore, the similarity between all effector candidates from both rhizophagus species was ranged from 51 % for nls/scr_ to 97 % for nls_ as determined by the blastp algorithm. from these, conserved effector candidates show more than 90 % similarity between r. clarus and r. irregularis, which falls within the same range as the similarity between housekeeping genes. more than 44 % of the conserved candidates belong to the nls group. in addition, the nls group shows the highest conservation level in comparison to the scr and rcp groups. these findings suggest that the effectors translocated to the plant nucleus with possible effect on the transcriptional regulation in plant cells might be essential for am symbiosis. the highly conserved effector candidates identified here between different amf species are probably playing a fundamental role during fungus-plant interaction also explaining the degree of unspecific interaction with all potential hosts. in contrast, effector candidates showing high diversification may have minor host specific functions.table the r. clarus effector candidates with predicted function potential functional features of r. clarus effector candidates to deduct putative functions of r. clarus and r. irregularis predicted effector candidates the blast2go software with automated annotation was used. the obtained information was used to predict putative cellular functions of rhizophagus effector candidates. only of r. clarus effector candidates showed sequence similarity to documented proteins with defined or predicted function, and contain known domains or motifs. predicted domain and homology hits with an e-value higher than 10− were ignored. current research showed that a number of fungal and oomycete effectors despite being under high diversifying or positive selective pressure still carry recognizable pfam domains, which might be helpful in prediction of their biological function. thus, for of the r. clarus effector candidates the putative function could be determined, hence their potential role in am symbiosis was defined and subdivided into three main functional groups: signalling, cell wall modification, and transcriptional regulation. however, the majority of r. clarus effector candidates are novel sequences showing no significant homology to known proteins in other organisms or sequence motifs, which is in agreement with previous studies highlighting the evolutionary diverse nature of fungal effector groups. based on computational approaches, we established effector candidate subcategories and speculated on putative functions after careful confirmation of presence of active sites and completeness of the predicted domains.fig. illustration of the various functional modules of r. clarus effector candidates secreted proteins likely to target plant signalling secreted effectors can alter plant defence by interfering with signal transduction. within the r. clarus effectome we detected eight effector candidates that might affect plant signalling. within the r. clarus effectome, candidate scr_ shows homology to the thermoanaerobacter sp. peptidoglycan-binding lysm protein and contains a lysm domain. plant receptor kinases containing a lysm domain are able to detect chitin-like components or chitin fragments. chitin, an n-acetyl-glucosamine polymer, is a main component of the fungal cell wall. chitin recognition activates signalling pathways, which lead to plant immune defence. due to this, it is suggested that many plant pathogens secrete effectors that contain lysm domains in order to inhibit the release of chitin fragments from fungal cell walls or to bind these fragments, thus, preventing their recognition by plant receptors. for instance, the rice blast fungus magnaporthe oryzae secretes the lysm domain containing effector, slp which binds chitin oligosaccharides. this suppresses plant defence response induced through the recognition of chitin fragments by the rice receptor cebip. candidate scr_ might play a similar role in am symbiosis. the second effector candidate nls_ contains a f-box-like domain. proteins containing this domain are components of the e ubiquitin ligases called scf complexes which control many intracellular regulatory processes by catalysing the ubiquitination of proteins that targets proteins for degradation. intense analysis showed that plants defence response pathways might be under control of an scf complex. many pathogenic effectors use a strategy of ubiquitination in order to regulate plant immunity. the avrpiz-t effector of the m. oryzae binds to and destabilizes the rice ring e ubiquitin ligase apip to suppress chitin-induced plant immunity. the third effector candidate nls_ shows homology to the ochotona princeps nitrogen permease regulator 3-like protein isoform x and contains a npr domain typical of this protein family. npr forms a heterodimer complex with npr protein, which mediates an amino acid starvation signal to torc. torc an essential serine/threonine protein kinase, incorporates signals from diverse pathways in the cell. the homology of candidate nls_ to npr protein suggests its interference with signal transduction most likely during amino acid deficiency conditions. a similar role in interfering with plant signalling pathways might have the fourth candidate nls_ in which a protein kinases catalytic domain, pkc was detected. pkc domains can transfer a phosphate group from atp to serine/threonine or tyrosine residues on protein substrates. pks often function as components of signal transduction pathways in which one kinase activates another kinase. this serial action transfers a signal within a cell, which results in cellular responses. two additional candidates might also influence the plant signalling cascades by being involved in hydrolytic processes. the candidate nls_ shows homology to the dichomitus squalens plc-like phosphodiesterase and has a glycerophosphodiester phosphodiesterase domain, gdpd. plc, by hydrolysing phosphodiester bonds, integrates signalling cascades critical for a variety of cellular and physiological functions. the candidate rcp_ shows homology to the grosmannia clavigera phospholipase d which catalyse the cleavage of phospholipids into fatty acids and other lipophilic substances, like phosphatidic acid, which is an established intracellular signalling lipid. candidate number seven, nls_ has homology to the glycosyltransferase family protein from colletotrichum higginsianum and encompasses a complex bifunctional domain, which includes a trehalose-6-phosphate synthase domain and a had-superfamily hydrolase domain, subfamily iib. tps catalyses the synthesis of α,α- -trehalose-6-phosphate from glucose-6-phosphate using a udp-glucose donor and is a key enzyme in the trehalose synthesis pathway. furthermore, trehalose was reported to have a role as a stress protectant in many organisms across kingdoms. remarkably, in pathogens trehalose is associated with several mechanisms involved in cell wall integrity, regulation of metabolism, and successful infection. absence of a trehalose metabolic apparatus is correlated with lower pathogen virulence. however, yet the underlying mechanism remains unravelled. for instance, in m. oryzae tps regulates transcriptional effectors linked to virulence factors via nadph-binding during appressorium-mediated rice infection. consequently, the putative effector nls_ might take part in the colonization step of amf symbiosis by being involved in trehalose synthesis. the last candidate, scr_ contains a ml domain, which is involved in specific lipid recognition. it is predicted that ml domain containing proteins through interacting with specific lipids can mediate diverse biological functions, like lipid metabolism, host response to pathogen components like lipopolysaccharides, and other cellular functions. this suggests that the scr_ candidate, by interacting with lipids, can interfere with the plant lipid signalling; however the direct mechanism of this effect needs to be clarified. genes encoding putative secreted enzymes engaged in cell wall modification we identified four r. clarus effector candidates encoding putative secreted enzymes involved in cell wall modification. they were divided in two subgroups depending whether their putative biological functions concern plant or fungal cell wall. enzymes involved in degradation of the plant cell wall most of the necrotrophic fungi degrade plant tissues by producing a broad set of carbohydrate-active enzymes specifically focused on plant polysaccharide degradation. transcriptomic profiling of r. irregularis daom predicted genes encoding carbohydrate-active enzymes. this amount is relatively low in comparison to obligate biotrophic pathogens and ectomycorrhizal symbionts. also, only one putative enzyme involved in plant cell wall modification was found in the r. clarus effector repertoire. this suggests that during am symbiosis the plant cell wall is not completely degraded, but instead a localized corruption of its structure might occur to make it more flexible, possibly for fungus penetration into the plant cell. the single effector of this subgroup, candidate nls/rcp_ shows homology to the rhizoctonia solani galactose oxidase and contains a glyoxal oxidase domain. glyoxal oxidase has a copper active site remarkably similar to that found in another fungal metalloenzyme, galactose oxidase. glyoxal oxidase is a copper metalloenzyme that catalyses the oxidation of aldehydes to carboxylic acids and is an essential component of the extracellular lignin degradation pathways. although lignin is present only in small amounts in the plant cell wall, it has a strong impact on the rigidity of its structure. lignin degradation might result in a more flexible cell wall necessary for fungal penetration. enzymes involved in modification of the fungal cell wall the effectome also contains three putative enzymes involved in chitin degradation. this suggests constant remodelling of the fungal cell wall, likely an adaptation to substantial differences in environmental conditions in the soil and inside of plant tissues. similarly, in cryptococcus neoformans it was shown that chitosan, the product of chitin deacetylation, helps to sustain cell integrity by maintaining normal capsule width during vegetative growth. the first effector candidate in r. clarus, rcp_ contains a glycosyl hydrolase family domain and shows homology to the neofusicoccu parvum class i alpha-mannosidase enzyme. this enzyme is involved in the cleavage of the terminal mannose from asn-linked oligosaccharides, which takes part at the early step of hydrolysis of n-glycan such as chitin. two other r. clarus effector candidates might function in fungal cell wall remodelling by chitin deacetylation. candidate rcp_ shows homology to xylase/chitin deacetylase enzyme from cryptococcus gatti and contains a chitin deacetylase domain. the other candidate, rcp_ shows a significant homology to rhizoctonia solani bifunctional xylanase/deacetylase. the presence of hevein/chitin binding domain and a catalytic nodb homology domain of the carbohydrate esterase superfamily, which typically occur in chitin deacetylases, suggest a chitin degradation function of both candidates. however, candidates with putative chitin deacetylase function might also have another purpose. it was shown that chitin deacetylase plays a significant role in protecting pathogenic fungal hyphae from being lysed by plant secreted chitinases, specifically in the cases of the wheat stem rust fungus puccinia graminis f. sp. tritici and the broad bean rust fungus uromyces fabae. furthermore, it was suggested that exposed fungal chitin polymers are partially de-n-acetylated by chitin deacetylases during the infection of the host, this with purpose to avoid plant antimicrobial hydrolases by affinity modulation. genes encoding proteins with putative transcriptional regulatory function during the pathogenic or symbiotic infection thousands of genes are up- or down-regulated in the plant as revealed by several transcriptome analyses. expression of genes encoding enzymes involved in protein degradation, cell wall modification, and secondary metabolite biosynthesis is induced in order to activate plant immune response. however, pathogens are able to secrete the so-called transcription activation-like effectors, which target dna and are involved in transcriptional regulation of a set of plant host genes, for instance to prevent defence response. we found only one putative candidate within the r. clarus effectome that might regulate plant gene expression during establishment of the am symbiosis. this is the nuclear localized candidate, nls_ which consists of a pwwp domain, named after a conserved pro-trp-trp-pro motif. pwwp domain-containing proteins bind to chromatin and control its packing resulting in transcription and replication regulation. moreover, pwwp proteins are also involved in dna damage response. while dna damage occurs, they regulate chromatin condensation and ability of dna repair factors, like 53bpi and ctip to bind dna. similarly, candidate nls_ might be involved in transcription regulation by chromatin remodelling or condensation. effector candidates with unknown function six of the r. clarus effector candidates show homologies to known proteins or contain identified domains or motifs. however, the similarity does not allow for prediction of a specific function. one of the candidates, nls_ has an aaa domain and shows significant homology to a protozoa protein containing a denn domain. however, the homology is mostly based on the aaa domain region as a denn domain was not detected. the aaa superfamily is a large and functionally diverse superfamily composed of members that are essential for a variety of cellular functions, like protein unfolding and degradation, dna recombination, replication and repair, and signal transduction. candidates scr_ and rcp_ show homologies to putative conserved fungal proteins with unknown function. both of them lack of any functional domain; therefore, their putative function cannot be predicted. candidate nls_ shows homology to the trametes versicolor p-loop protein containing a nucleoside triphosphate hydrolase domain. nevertheless, this homology is partial and does not include the characteristics for this family of p-loop ntpase domain proteins. in addition, two other candidates, scr_ and rcp_ show homology to the laccari bicolor carbohydrate binding module family protein. proteins belonging to this family contain a module of 60– residues with a chitin binding function. however, once more, this module is absent in both candidates; hence, despite the homology the putative function cannot be predicted.  symbiotic interaction between amf and plant host requires a well-adjusted dialog between the two partners in order to function properly. we hypothesise that one of the key components of this communication is the fungal effectome. in this study, we used a bioinformatic pipeline to predict effector candidates from two amf belonging to the rhizophagus genus, r. irregularis and r. clarus. our in silico pipeline revealed a list of several candidate effector proteins that create a valuable source of information to elucidate the mechanism of plant infection and colonization by am fungi. our in silico characterization of the amf effectome focused on r. clarus candidates. the majority of the r. clarus effectors are novel proteins with no significant homology to known sequences or motifs in other organisms. this finding is in agreement to previous reports enhancing the evolution diverse nature of fungal and oomycete effectors. however, for a small number of candidates putative functions during plant-fungus symbiosis might be predicted. the set of putative effector candidates encoded by the r. clarus effectome appears to be suited for a fungus that wants to establish a symbiotic rather than a pathogenic relation with the host. the presence of only one candidate encoding a putative plant cell wall degrading enzyme in the r. clarus effectome is in agreement with earlier reports, in which it was shown that several obligate biotrophic pathogens and ectomycorrhizal symbionts have a decreased repertoire of pcwdes compared to saprophytes. in comparison to pathogenic fungi, amf seem not to degrade the plant cells but rather grow between them or inside of them while trying to avoid recognition by the plant defence. however, localized plant cell wall modelling probably does occur to weaken and make it more flexible during fungal penetration and passage from cell to cell. the lack of the typical pathogenic effectors, such as extracellular lipases, proteases, nucleases and phytases, additionally confirms the non-pathogenic status of amf. the relationship between symbiotic fungi and plants requires a major rearrangement of the plant signalling pathways for which fungi have evolved a range of secreted proteins. for instance, fungal effectors control auxin signalling within plant cells to alter plant root architecture and favour colonization. thus, the presence of candidates in the r. clarus effectome that might target plant signalling in order to allow fungal entry was expected. reprogramming of the plant defence could be used to avoid recognition. alternatively, effectors might affect transport processes across the plant cell membrane. this was previously shown for symbiotic fungi which encode effectors to control the transfer of glucose into the apoplastic space of the root where it is easily accessible for fungus. a remarkable aspect of our work is presence of a large portion of the effectors that are not species-specific, but conserved between different rhizophagus species. this group of highly conserved effector proteins may play a fundamental role during fungus-plant interaction and has first priority for future functional characterization. availability of supporting data the sequences of the predicted effector candidates have been deposited in genbank under accession numbers: ku - ku. other supporting data are included in an additional file  additional file additional file 1: figure s gene ontology category analyses by protfun server of r. irregularis daom secretome and effectome. table s r. clarus draft genome assembly statistics in comparison to published r. irregularis genome assemblies. table s the r. irregularis putative effectome. table s identity and similarity levels between r. irregularis and r. clarus based on sequences of selected proteins. table s the r. clarus effector candidates. abbreviations aaaatpases associated with diverse cellular activities amarbuscular mycorrhizal amfarbuscular mycorrhizal fungi cebipchitin elicitor binding protein denndifferentially expressed in neoplastic vs. normal cells erendoplasmic reticulum etieffector-triggered immunity gdpdglycerophosphodiester phosphodiesterase domain glcnacn-acetyl-glucosamine lysmlysine motif m. oryzaemagnaporthe oryzae miomillion mlmd-2-related lipid-recognition myamillion years ago nb-lrrnucleotide-binding site leucine-rich repeat nlsnuclear localization signal npr2nitrogen permease regulator 2 npr3nitrogen permease regulator 3-like protein pcwdeplant cell wall degrading enzyme pkcprotein kinases catalytic domain plc-likephospholipase c-like p-loop ntpasep-loop protein with a nucleoside triphosphate hydrolase domain pwwppro-trp-trp-pro qcquality control r. clarusrhizophagus clarus r. irregularisrhizophagus irregularis rcprepeats containing proteins scrsmall cysteine-rich spsignal peptide sp7secreted effector protein 7 tmdtransmembrane domain torc1target of rapamycin kinase tpstrehalose-6-phosphate synthase competing interests the authors declare that they have no competing interests. authors’ contributions conceived and designed the experiments: ab, kst. genome sequencing and assembly: ab. bioinformatics and data analysis: kst. wrote manuscript: kst, ab. both authors read and approved the final manuscript. 
we hypothesise that also for symbiotic-plant interactions the secreted effectome of the fungus is a major component of communication and the conservation level of effector proteins between amf species may be indicative whether they play a fundamental role. in this study, we used a bioinformatics pipeline to predict and compare the effector candidate repertoire of the two amf species, rhizophagus irregularis and rhizophagus clarus. a remarkable aspect of our work is presence of a large portion of the effector proteins involved in symbiosis, which are not unique to each fungi or plant species, but shared along the glomeromycota phylum. our in silico pipeline revealed a list of r. irregularis candidate effector genes that create a valuable information source to elucidate the mechanism of plant infection and colonization by fungi during amf symbiotic interaction. for 95 % of r. irregularis candidates we found homologs in a r. clarus genome draft generated by illumina high-throughput sequencing. amfeffectorin silico pipelineplantrhizophagussymbiosishttp://dx.doi.org/ /501100001659deutsche forschungsgemeinschaftbr 3527/1-1brachmann andreas http://dx.doi.org/http://dx.doi.org/ /501100001659deutsche forschungsgemeinschaftpa 493/11-1issue-copyright-statement© the author 2016 therefore, we state that this group of highly conserved effector proteins between amf species may play a fundamental role during fungus-plant interaction. identification and functional characterization of the primary amf effectors that regulate symbiotic development will help in understanding the mechanisms of fungus-plant interaction. while most of the candidate effectors show no homologies to known domains or proteins, the candidates with homologies point to potential roles in signal transduction, cell wall modification or transcription regulation. recent studies of plant-pathogen interactions showed that effector proteins play a key role in host colonization by controlling the plant immune system. interestingly, 9 % of the predicted effectors are at least as conserved between the two rhizophagus species as proteins with housekeeping functions. arbuscular mycorrhizal fungi form an ecologically important symbiosis with more than two thirds of studied land plants.
 arbuscular mycorrhizal fungi form an ecologically important symbiosis with more than two thirds of studied land plants. recent studies of plant-pathogen interactions showed that effector proteins play a key role in host colonization by controlling the plant immune system. we hypothesise that also for symbiotic-plant interactions the secreted effectome of the fungus is a major component of communication and the conservation level of effector proteins between amf species may be indicative whether they play a fundamental role.  in this study, we used a bioinformatics pipeline to predict and compare the effector candidate repertoire of the two amf species, rhizophagus irregularis and rhizophagus clarus. our in silico pipeline revealed a list of r. irregularis candidate effector genes that create a valuable information source to elucidate the mechanism of plant infection and colonization by fungi during amf symbiotic interaction. while most of the candidate effectors show no homologies to known domains or proteins, the candidates with homologies point to potential roles in signal transduction, cell wall modification or transcription regulation. a remarkable aspect of our work is presence of a large portion of the effector proteins involved in symbiosis, which are not unique to each fungi or plant species, but shared along the glomeromycota phylum. for 95 % of r. irregularis candidates we found homologs in a r. clarus genome draft generated by illumina high-throughput sequencing. interestingly, 9 % of the predicted effectors are at least as conserved between the two rhizophagus species as proteins with housekeeping functions. therefore, we state that this group of highly conserved effector proteins between amf species may play a fundamental role during fungus-plant interaction.  we hypothesise that in symbiotic interactions the secreted effectome of the fungus might be an important component of communication. identification and functional characterization of the primary amf effectors that regulate symbiotic development will help in understanding the mechanisms of fungus-plant interaction. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords amfeffectorin silico pipelineplantrhizophagussymbiosishttp://dx.doi.org/ /501100001659deutsche forschungsgemeinschaftbr 3527/1-1brachmann andreas http://dx.doi.org/http://dx.doi.org/ /501100001659deutsche forschungsgemeinschaftpa 493/11-1issue-copyright-statement© the author 2016 